Skip to content

Use FFT Poisson solver on fully refined levels #449

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
ChunYen-Chen wants to merge 8 commits into
base: main
Choose a base branch
from
Open

Use FFT Poisson solver on fully refined levels #449

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
8381411
Add `FFTW_Inited[]` and `GreenK_Inited[]`
ChunYen-Chen Mar 11, 2025
47d2131
Refactor `Init_FFTW()` and `End_FFTW()`
ChunYen-Chen Mar 14, 2025
0c44daa
Apply FFT Poisson solver to fully refined level
ChunYen-Chen Mar 14, 2025
31443e5
Rename `CPU_ExtPotSolver_BaseLevel.cpp` to
ChunYen-Chen Mar 14, 2025
daca378
Merge branch 'master' into full_fft_possion
ChunYen-Chen Jun 24, 2025
d09a23a
Bugfix
ChunYen-Chen Aug 25, 2025
485e5c2
Optimzie CellFactor
ChunYen-Chen Aug 25, 2025
1a00693
Optimize gravity part in EvolveLevel
ChunYen-Chen Aug 26, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 3 additions & 1 deletion include/Global.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -181,7 +181,9 @@ extern double AveDensity_Init; // initial average mass density (in al
extern int Pot_ParaBuf; // number of parallel buffers to exchange potential for the
// Poisson/Gravity solvers and the potential refinement
extern int Rho_ParaBuf; // number of parallel buffers to exchange density for the Poisson solver
extern real *GreenFuncK;
extern bool FFTW_Inited[NLEVEL];
extern bool GreenFuncK_Inited[NLEVEL];
extern real *GreenFuncK[NLEVEL];
extern double GFUNC_COEFF0;
extern double DT__GRAVITY;
extern double NEWTON_G;
Expand Down
17 changes: 9 additions & 8 deletions include/Prototype.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -260,14 +260,15 @@ void Init_ByRestart_HDF5( const char *FileName );
#endif
#ifdef SUPPORT_FFTW
void End_FFTW();
void Init_FFTW();
void Init_FFTW( const int lv );
void Patch2Slab( real *VarS, real *SendBuf_Var, real *RecvBuf_Var, long *SendBuf_SIdx, long *RecvBuf_SIdx,
int **List_PID, int **List_k, long *List_NSend_Var, long *List_NRecv_Var,
const int *List_z_start, const int local_nz, const int FFT_Size[], const int NRecvSlice,
const double PrepTime, const long TVar, const bool InPlacePad, const bool ForPoisson, const bool AddExtraMass );
const double PrepTime, const long TVar, const bool InPlacePad, const bool ForPoisson,
const bool AddExtraMass, const int lv );
void Slab2Patch( const real *VarS, real *SendBuf, real *RecvBuf, const int SaveSg, const long *List_SIdx,
int **List_PID, int **List_k, long *List_NSend, long *List_NRecv, const int local_nz, const int FFT_Size[],
const int NSendSlice, const long TVar, const bool InPlacePad );
const int NSendSlice, const long TVar, const bool InPlacePad, const int lv );
#endif // #ifdef SUPPORT_FFTW
void Microphysics_Init();
void Microphysics_End();
Expand Down Expand Up @@ -405,12 +406,12 @@ void CPU_PoissonGravitySolver( const real h_Rho_Array [][RHO_NXT][RHO_NXT][RH
const bool SelfGravity, const OptExtPot_t ExtPot, const OptExtAcc_t ExtAcc,
const double TimeNew, const double TimeOld, const real MinEint,
const bool UseWaveFlag );
void CPU_ExtPotSolver_BaseLevel( const ExtPot_t Func, const double AuxArray_Flt[], const int AuxArray_Int[],
const real Table[], void **GenePtr,
const double Time, const bool PotIsInit, const int SaveSg );
void CPU_ExtPotSolver_FullyRefinedLevel( const ExtPot_t Func, const double AuxArray_Flt[], const int AuxArray_Int[],
const real Table[], void **GenePtr,
const double Time, const bool PotIsInit, const int SaveSg, const int lv );
#ifdef SUPPORT_FFTW
void CPU_PoissonSolver_FFT( const real Poi_Coeff, const int SaveSg, const double PrepTime );
void Init_GreenFuncK();
void CPU_PoissonSolver_FFT( const real Poi_Coeff, const int SaveSg, const double PrepTime, const int lv );
void Init_GreenFuncK( const int lv );
#endif
void End_MemFree_PoissonGravity();
void Gra_AdvanceDt( const int lv, const double TimeNew, const double TimeOld, const double dt,
Expand Down
Loading