[petsc-users] 1d FFT of 2D vector (using fftw for preference)

Jed Brown jedbrown at mcs.anl.gov
Fri Jun 28 09:10:23 CDT 2013


Matthew Knepley <knepley at gmail.com> writes:

> I do not know of one other than repeated calls to the 1D transform. If
> you can find it in the FFTW API when could put it in.

I used this:

  c->p_forward = fftw_plan_guru_r2r(1, &(c->tdim), c->rank-1, c->dim, x, c->work, &redft00, flag | FFTW_PRESERVE_INPUT);

https://github.com/jedbrown/spectral-petsc/blob/master/chebyshev.c#L127
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130628/216785e8/attachment.pgp>


More information about the petsc-users mailing list