ready for MPI
I will put matrices here, too.
put matrices here
call FFTW routine directly
add ERRORs