Skip to content

FLUPS v1.0 - scalar and vector solver for 2-D and 3-D unbounded Poisson

Compare
Choose a tag to compare
@DGCaprace DGCaprace released this 12 Dec 22:45
· 713 commits to release since this release

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.

⚠️ This comes with an important change while computing the FFT'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