# Questions tagged [fftw]

The tag has no usage guidance.

34 questions
Filter by
Sorted by
Tagged with
9 views

### How to use RODFT00 and REDFT00

I have some difficulty in implementing RODFT00 and REDFT00. I want to use them for fluid simulations. I would really appreciate ...
22 views

...
31 views

### How to use discrete cosine and discrete sine transforms in fftw

I work on fluid-related simulations. I have used FFT for fluid simulation. I want to use discrete cosine transform (DCT) and discrete sine transform (DST) to transform my velocity field to wavenumbers....
110 views

### MKL/FFTW performance of batch 1-D FFTs

MKL and FFTW offer 1-D FFTs that can operate on many inputs simultaneously - in other words, they can batch-transform the columns of some input matrix. Is the performance of these multi-transforms ...
48 views

### FFTW on subarray with MPI

With the guru interfaces of FFTW, I can apply transforms only to parts of a multidimensional array by modifying the fftw_iodim ...
70 views

### Discrete sine and cosine transform for mixed derivatives

Using sine and cosine transforms to solve Poisson's equation with Dirichlet boundary conditions seem quite standard nowadays (see, e.g., here or Table 2 in this paper). In the case of Poisson's ...
179 views

### Use of multidimensional FFTW and normalisation factor

I am using the FFTW MPI in C. I have a simple question. Quoting from fftw.org The multi-dimensional transforms of FFTW, in general, compute simply the separable product of the given 1d ...
190 views

233 views

### MPI support for discrete Fourier transform (DFT) in Python

I am looking for a discrete Fourier transform (DFT) library that can be run with MPI on Python. Usually, in other languages (C, Fortran) FFTW is used. There's a Python wrapper for FFTW called pyFFTW, ...
196 views

### Solving a 3D (almost radial) convolution with FFT

I have a 3D integral that is almost a radial convolution of the form $$\int d^{3}k'h(\mathbf{k'})g(|\mathbf{k-k'}|)$$ and I am looking for a fast and efficient algorithm (e.g. FFT) to solve it ...
10k views

### C++ libraries for Fast Fourier Transform in high precision

I am looking for a C++ library for Fast Fourier Transform (FFT) in high precision (e.g., using high precision real data types similar to mpfr_t in MPFR or ...
284 views

### How to obtain values in physical space for a given spectrum?

My question falls under purview of turbulent flows. I want to add an initial perturbation, for which I have a given energy spectrum (say$E(k)=ak^4e^{-bk^2}$). The steps involved in getting these ...
3k views

### computing turbulent energy spectrum from isotropic turbulence flow field in a box

I have my 3 dimensional velocity flow-field u, v and w at a given instant of time from DNS using pseudo-spectral method. I need to calculate the energy spectrum ( in Fourier space ) as a function of ...
129 views

### Computing Fourier representation of space dependent advection operator via FFT

Consider the following equation on the circle: $$\dfrac{\partial p(x,t)}{\partial t} = a(x)\dfrac{\partial p(x,t)}{\partial x} \equiv L(p) \enspace ,$$ where $L$ is the operator acting on $p(x,t)$. ...
59 views

### Broadening spectral data by using FFT's

I obtain numerical discrete data of the form $$S_{raw}(\omega) = \sum_{j}w_{j} \delta(\omega-\omega_{j})$$ to compare the result with experimental data the delta peaks need to be broadened ...
607 views

### radially averaged power spectrum of a binary image does not look like power law

It has been known that Fourier power spectrum somehow obeys power law therefore the slope of the spectrum can be used to calculate the fractal dimension of an image. Many people have used it for ...
255 views

### precision loss in non-trigonometric, periodic functions using FFTW and NaNs after marching forward in time (Fortran)

I have developed a pseudospectral solver of the Navier-Stokes equations using FFTW. I tested my formulation of right hand sides (RHS) of the NS equations against standard trigonometric functions (...
459 views

### FFT parallel processing in MPI

I am working now in Beowulf Cluster and parallel processing, I want code for Fast Fourier transfer functions written in any language, e.g., C/C++. Without using FFTW library based on message passing ...
653 views

### Time array from frequency array in FFT using Python

I have done a Fourier transformation of two signals (in time) $S_1(t),S_2(t)$ using numpy's fft which will give me $S_1(f),S_2(f)$. The corresponding frequency grid ...
74 views

### Which version of FFTW

At the moment, I'm installing the "GADGET-2" application. In its documentation it says "Note that the MPI-capable version 2.x of FFTW is required, the new version 3 lacks MPI capability at this point. ...
462 views

### DST using FFT routine

Please can you help me with my problem? On Wikipedia, in article Discrete sine transform, this is written (chapter Computation): "Although the direct application of these formulas would require O(N2) ...
80 views

### A function as a sum of serie of modified FFT coeff. of another function - multiplied by sum number

I solve such a problem. Lets have a function $Y=\sum_{k=-\infty}^\infty i\hat Y e^{ik\pi y}$ and then I have a function which is defined as $X=\sum_{k=-\infty}^\infty ik^2\hat Y e^{ik\pi y}$. I ...
2k views

### MPI-based Implementations of FFT

In a numerical computation, I am required to take a multi-dimensional FFT on a distributed-memory cluster. The data is currently distributed using a distributed array in PETSc (DMDA). I initial ...
151 views

### FFT - function only in sine series? Can be done with MKL / Lapack? [duplicate]

please can I ask, how one can make from function sine series (Fourier transform) with MKL? I can do "normal" exponential FFT with MKL (Lapack of course), how can I say that I want only sine series? ...
182 views

### Discrete convolution

please can I ask a bit stupid question? Let say I need to solve an equation in a form $\frac{\partial X}{\partial t}=\sum_k M_k * X_{n-k}$ How can I do the discrete convolution numerically? I will say ...
2k views

### Fast (approximate) evaluation of Chebyshev polynomial

Is there a preferred way how to implement a fast (approximate) evaluation of the Chebyshev interpolation polynomial on uniform grid (given the function values at the Chebyshev nodes)? My problem is ...
To use the Fast Fourier Transform (FFT) on uniformly sampled data, e.g. in connection with PDE solvers, it is well known that the FFT is an $\mathcal{O}(n\log(n)$) algorithm. How well do the FFT scale ...