FLUPS v1.0 - scalar and vector solver for 2-D and 3-D unbounded Poisson
Changelog
In this release, we introduce the resolution of a Poisson equation in 2/3D for a vector field, each component having a different set of BC's.
In order to store only once the Green's function in the fourier domain, we have to impose a common Fourier space to every component of the vector, starting from the 0-mode to the flip-flop mode.
To achieve this, regardless of the DCT/DST/DFT used, we correct the output of the DCT by adding the flip-flop mode and the output of the DST with a shift and the addition of the 0-mode.
This operation is local, hence scalable, and is done after/before every DCT/DST computation while in the correct topology.
To be achieved in the near future
- sample code for FORTRAN
- support of single-precision domains (float vs double).
- support of FORTRAN and C way of indexing
- support of vertex-centered meshes