edu.emory.mathcs.jtransforms.fft

Class RealFFTUtils_3D

java.lang.Object

edu.emory.mathcs.jtransforms.fft.RealFFTUtils_3D

public class RealFFTUtils_3D
extends Object

This is a set of utility methods for R/W access to data resulting from a call to the Fourier transform of real data. Memory optimized methods, namely

• DoubleFFT_3D.realForward(double[])
• DoubleFFT_3D.realForward(double[][][])
• FloatFFT_3D.realForward(float[])
• FloatFFT_3D.realForward(float[][][])

are implemented to handle this case specifically. However, packing of the data in the data array is somewhat obscure. This class provides methods for direct access to the data, without the burden of all necessary tests.

Example for Fourier Transform of real, double precision 1d data

   DoubleFFT_3D fft = new DoubleFFT_2D(slices, rows, columns);
   double[] data = new double[2 * slices * rows * columns];
   ...
   fft.realForwardFull(data);
   data[(s1 * rows + r1) * 2 * columns + c1] = val1;
   val2 = data[(s2 * rows + r2) * 2 * columns + c2];
 

is equivalent to

   DoubleFFT_3D fft = new DoubleFFT_3D(slices, rows, columns);
   RealFFTUtils_3D unpacker = new RealFFTUtils_3D(slices, rows, columns);
   double[] data = new double[slices * rows * columns];
   ...
   fft.realForward(data);
   unpacker.pack(val1, s1, r1, c1, data);
   val2 = unpacker.unpack(s2, r2, c2, data, 0);
 

Even (resp. odd) values of c correspond to the real (resp. imaginary) part of the Fourier mode.

Example for Fourier Transform of real, double precision 3d data

   DoubleFFT_3D fft = new DoubleFFT_3D(slices, rows, columns);
   double[][][] data = new double[slices][rows][2 * columns];
   ...
   fft.realForwardFull(data);
   data[s1][r1][c1] = val1;
   val2 = data[s2][r2][c2];
 

is equivalent to

   DoubleFFT_3D fft = new DoubleFFT_3D(slices, rows, columns);
   RealFFTUtils_3D unpacker = new RealFFTUtils_3D(slices, rows, columns);
   double[][][] data = new double[slices][rows][columns];
   ...
   fft.realForward(data);
   unpacker.pack(val1, s1, r1, c1, data);
   val2 = unpacker.unpack(s2, r2, c2, data, 0);
 

Even (resp. odd) values of c correspond to the real (resp. imaginary) part of the Fourier mode.

Author:

Sébastien Brisard

Constructor Summary

Constructors

Constructor and Description
RealFFTUtils_3D(int slices, int rows, int columns) Creates a new instance of this class.

Method Summary

Modifier and Type	Method and Description
int	getIndex(int s, int r, int c) Returns the 1d index of the specified 3d Fourier mode.
void	pack(double val, int s, int r, int c, double[][][] packed) Sets the specified Fourier mode of the transformed data.
void	pack(double val, int s, int r, int c, double[] packed, int pos) Sets the specified Fourier mode of the transformed data.
void	pack(float val, int s, int r, int c, float[][][] packed) Sets the specified Fourier mode of the transformed data.
void	pack(float val, int s, int r, int c, float[] packed, int pos) Sets the specified Fourier mode of the transformed data.
double	unpack(int s, int r, int c, double[][][] packed) Returns the specified Fourier mode of the transformed data.
double	unpack(int s, int r, int c, double[] packed, int pos) Returns the specified Fourier mode of the transformed data.
float	unpack(int s, int r, int c, float[][][] packed) Returns the specified Fourier mode of the transformed data.
float	unpack(int s, int r, int c, float[] packed, int pos) Returns the specified Fourier mode of the transformed data.

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Constructor Detail

RealFFTUtils_3D

public RealFFTUtils_3D(int slices,
                       int rows,
                       int columns)

Creates a new instance of this class. The size of the underlying DoubleFFT_3D or FloatFFT_3D must be specified.

Parameters:

slices - number of slices

rows - number of rows

columns - number of columns

Method Detail

getIndex

public int getIndex(int s,
                    int r,
                    int c)

Returns the 1d index of the specified 3d Fourier mode. In other words, if packed contains the transformed data following a call to DoubleFFT_3D.realForwardFull(double[]) or FloatFFT_3D.realForward(float[]), then the returned value index gives access to the [s][r][c] Fourier mode

• if index == Integer.MIN_VALUE, then the Fourier mode is zero,
• if index >= 0, then the Fourier mode is packed[index],
• if index < 0, then the Fourier mode is -packed[-index],

Parameters:

s - the slice index

r - the row index

c - the column index

Returns:

the value of index

pack

public void pack(double val,
                 int s,
                 int r,
                 int c,
                 double[] packed,
                 int pos)

Sets the specified Fourier mode of the transformed data. The data array results from a call to DoubleFFT_3D.realForward(double[]).

Parameters:

val - the new value of the [s][r][c] Fourier mode

s - the slice index

r - the row index

c - the column index

packed - the transformed data

pos - index of the first element in array packed

pack

public void pack(double val,
                 int s,
                 int r,
                 int c,
                 double[][][] packed)

Sets the specified Fourier mode of the transformed data. The data array results from a call to DoubleFFT_3D.realForward(double[][][]).

Parameters:

val - the new value of the [s][r][c] Fourier mode

s - the slice index

r - the row index

c - the column index

packed - the transformed data

pack

public void pack(float val,
                 int s,
                 int r,
                 int c,
                 float[] packed,
                 int pos)

Sets the specified Fourier mode of the transformed data. The data array results from a call to FloatFFT_3D.realForward(float[]).

Parameters:

val - the new value of the [s][r][c] Fourier mode

s - the slice index

r - the row index

c - the column index

packed - the transformed data

pos - index of the first element in array packed

pack

public void pack(float val,
                 int s,
                 int r,
                 int c,
                 float[][][] packed)

Sets the specified Fourier mode of the transformed data. The data array results from a call to FloatFFT_3D.realForward(float[][][]).

Parameters:

val - the new value of the [s][r][c] Fourier mode

s - the slice index

r - the row index

c - the column index

packed - the transformed data

unpack

public double unpack(int s,
                     int r,
                     int c,
                     double[] packed,
                     int pos)

Returns the specified Fourier mode of the transformed data. The data array results from a call to DoubleFFT_3D.realForward(double[]).

Parameters:

s - the slice index

r - the row index

c - the column index

packed - the transformed data

pos - index of the first element in array packed

Returns:

the value of the [s][r][c] Fourier mode

unpack

public double unpack(int s,
                     int r,
                     int c,
                     double[][][] packed)

Returns the specified Fourier mode of the transformed data. The data array results from a call to DoubleFFT_3D.realForward(double[][][]) .

Parameters:

s - the slice index

r - the row index

c - the column index

packed - the transformed data

Returns:

the value of the [s][r][c] Fourier mode

unpack

public float unpack(int s,
                    int r,
                    int c,
                    float[] packed,
                    int pos)

Returns the specified Fourier mode of the transformed data. The data array results from a call to FloatFFT_3D.realForward(float[]) .

Parameters:

s - the slice index

r - the row index

c - the column index

packed - the transformed data

pos - index of the first element in array packed

Returns:

the value of the [s][r][c] Fourier mode

unpack

public float unpack(int s,
                    int r,
                    int c,
                    float[][][] packed)

Returns the specified Fourier mode of the transformed data. The data array results from a call to FloatFFT_3D.realForward(float[][][]) .

Parameters:

s - the slice index

r - the row index

c - the column index

packed - the transformed data

Returns:

the value of the [s][r][c] Fourier mode