Commit 5268b5e6 authored by Andrew Heather's avatar Andrew Heather
Browse files

BUG: Corrected no-op in reverse transform multi-D fft - see #813

parent ef882e54
......@@ -185,16 +185,16 @@ void Foam::fft::transform
const label N = field.size();
fftw_complex in[N], out[N];
forAll(field, i)
// If reverse transform : renumber before transform
if (dir == REVERSE_TRANSFORM)
{
in[i][0] = field[i].Re();
in[i][1] = field[i].Im();
fftRenumber(field, nn);
}
// If backward transform : renumber before transform
if (dir == FFTW_BACKWARD)
forAll(field, i)
{
fftRenumber(field, nn);
in[i][0] = field[i].Re();
in[i][1] = field[i].Im();
}
// Create the plan
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment