naginterfaces.library.sum.fft_​hermitian_​3d

naginterfaces.library.sum.fft_hermitian_3d(n1, n2, n3, y)

fft_hermitian_3d computes the three-dimensional inverse discrete Fourier transform of a trivariate Hermitian sequence of complex data values.

For full information please refer to the NAG Library document for c06pz

https://support.nag.com/numeric/nl/nagdoc_30/flhtml/c06/c06pzf.html

Parameters

n1int

$n_1$, the first dimension of the transform.

n2int

$n_2$, the second dimension of the transform.

n3int

$n_3$, the third dimension of the transform.

ycomplex, array-like, shape $((n1/2+1)\times n2\times n3)$

The Hermitian sequence of complex input dataset $z$, where $z_{{\text{j}}_1{\text{j}}_2{\text{j}}_3}$ is stored in $y[{\text{j}}_3\times (n_1/2+1)n_2+{\text{j}}_2\times (n_1/2+1)+{\text{j}}_1+1-1]$, for ${\text{j}}_3=0,1,\dots ,n_3-1$, for ${\text{j}}_2=0,1,\dots ,n_2-1$, for ${\text{j}}_1=0,1,\dots ,n_1/2$.

Returns

xfloat, ndarray, shape $(n1\times n2\times n3)$

The real output dataset $\hat{x}$, where ${\hat{x}}_{{\text{k}}_1{\text{k}}_2{\text{k}}_3}$ is stored in $x[{\text{k}}_3\times n_1{n}_2+{\text{k}}_2\times n_1+{\text{k}}_1+1-1]$, for ${\text{k}}_3=0,1,\dots ,n_3-1$, for ${\text{k}}_2=0,1,\dots ,n_2-1$, for ${\text{k}}_1=0,1,\dots ,n_1-1$.

Raises

NagValueError

(errno $1$)

On entry, $n1=⟨value⟩$.

Constraint: $n1\ge 1$.

(errno $2$)

On entry, $n2=⟨value⟩$.

Constraint: $n2\ge 1$.

(errno $3$)

On entry, $n3=⟨value⟩$.

Constraint: $n3\ge 1$.

(errno $4$)

An internal error has occurred in this function. Check the function call and any array sizes. If the call is correct then please contact NAG for assistance.

Notes

fft_hermitian_3d computes the three-dimensional inverse discrete Fourier transform of a trivariate Hermitian sequence of complex data values $z_{j_1{j}_2{j}_3}$, for ${\text{j}}_3=0,1,\dots ,n_3-1$, for ${\text{j}}_2=0,1,\dots ,n_2-1$, for ${\text{j}}_1=0,1,\dots ,n_1-1$.

The discrete Fourier transform is here defined by

$${\hat{x}}_{k_1{k}_2{k}_3}=\frac{1}{\sqrt{n_1{n}_2{n}_3}}\sum _0^{n_1-1}\sum _0^{n_2-1}\sum _0^{n_3-1}{z}_{j_1{j}_2{j}_3}\times exp\left(2\pi i(\frac{j_1{k}_1}{n_1}+\frac{j_2{k}_2}{n_2}+\frac{j_3{k}_3}{n_3})\right)\text{,}$$

where $k_1=0,1,\dots ,n_1-1$, $k_2=0,1,\dots ,n_2-1$ and $k_3=0,1,\dots ,n_3-1$. (Note the scale factor of $\frac{1}{\sqrt{n_1{n}_2{n}_3}}$ in this definition.)

Because the input data satisfies conjugate symmetry (i.e., $z_{j_1{j}_2{j}_3}$ is the complex conjugate of $z_{(n_1-j_1)(n_2-j_2)(n_3-j_3)}$), the transformed values ${\hat{x}}_{k_1{k}_2{k}_3}$ are real.

A call of fft_real_3d() followed by a call of fft_hermitian_3d will restore the original data.

This function calls fft_realherm_1d_multi_col() and fft_complex_1d_multi_row() to perform multiple one-dimensional discrete Fourier transforms by the fast Fourier transform (FFT) algorithm in Brigham (1974) and Temperton (1983).

References

Brigham, E O, 1974, The Fast Fourier Transform, Prentice–Hall

Temperton, C, 1983, Fast mixed-radix real Fourier transforms, J. Comput. Phys. (52), 340–350