I have problems executing a programm. I know the code is working (because it’s just a free code from the internet) and I’m able to compile the code with the clang++ compiler, but if I want to execute the Programm I get the following .if !dPS .ds PS .if !dPE .de PE .lf – 1 ..

#### Category : fftw

I have a simple "complex to real" FFT function that contains a dummy variable initialized in order to avoid any data being overwritten. This function works fine when used with my serial code, but I am trying to write a parallel version of this. When I do this and run Valgrind with openmp I get ..

Recently I’ve been trying to deepen my understanding of the MPI standard and I’ve came across topologies which probably can simplify some parts of my existing code which among other things uses FFTW3 with MPI. It’s a single grid based Schrödinger equation solver, which I’m now extending for multiple grids that should communicate using inter-communicators. ..

GOAL I have 2d data that represent a real scalar field f(x,y) and I want to calculate the 1D spectrum in the fourier (q) space. My project is written mainly in C++ but there are some blocks of code in C. I am also using the FFTW3 library and perform in-place transforms. CODE I have ..

GOAL I have 2d data that represent a real scalar field f(x,y) and I want to calculate the 1D spectrum in the fourier (q) space. My project is written mainly in C++ but there are some blocks of code in C. I am also using the FFTW3 library and perform in-place transforms. CODE I have ..

This question builds up on the previous post: FFTW3 segmentation errors while using std::vector container for FFT I cannot localize the problems with my FFT routines, done repeatedly in a loop. I compile using g++ -std=c++17 smoother.cpp -lfftw3_omp -lfftw3f -lfftw3 -lm -fopenmp -lstdc++fs -o smoother.exe The binary files I’m working on are single precision and ..

I would like to ask for your advice regarding the use of std::vector container in performing FFT. In my program, I use a certain function (get_PDF_CDF) in a loop. When I perform the FFTs using std::vector containers, during a certain iteration (usually the second one), I am getting a segmentation error / memory corruption errors. ..

i’m trying to port some code from CPU to GPU that includes some FFTs. So, on CPU code some complex array is transformed using fftw_plan_many_r2r for both real and imag parts of it separately. Function foo represents R2R transform routine and called twice for each part of complex array. void foo(vector_double &evg) { int nx ..

Consider a container for two-dimensional complex-valued arrays #include <vector> #include <array> struct Array2D { typedef std::array<double, 2> complex; int X, Y; std::vector<complex> values; }; and a function that uses FFTW to compute the Fourier transform: #include <fftw3.h> void FourierTransform(const Array2D& source, Array2D *dest) { fftw_plan p = fftw_plan_dft_2d(source.X, source.Y, (fftw_complex*)source.values.data(), (fftw_complex*)dest->values.data(), FFTW_FORWARD, FFTW_ESTIMATE); fftw_execute(p); fftw_destroy_plan(p); ..

I have a function which does a Hilbert transform on a std::vector containing discrete samples. I have learned that a band pass filter on my samples can increase the quality of the Hilbert Transform. The code of the Hilbert transform function is shown below. How can I add a band pass filter to this function? ..

## Recent Comments