Adding discrete cosine and sine transforms
There are possible applications of adding pseudospectral solvers with boundary conditions and to use it for postprocessing non-periodic fields. We could start by using mpi4py-fft's API which spectralDNS uses. Thoughts?