Implement an iterative Poisson solver - #69

Draft

pbartholomew08 wants to merge 132 commits into

xcompact3d:mainfrom

pbartholomew08:44-poisson-solver-iterative-solver

Commits on Feb 19, 2024

Add a Poisson solver selection mechanism.
semi-h
committed
7a342be
Implement a dedicated class for the FFT based Poisson solver.
semi-h
committed
bcd1056
Store compact scheme coefficients in tdsops class.
semi-h
committed
e6c488b
Store domain length and dx in dirps_t.
semi-h
committed
d4c6831
Change the way we initialise poisson_fft class.
semi-h
committed
d0d11d9
Construct the waves array for all periodic BCs.
semi-h
committed
9cbdf95
Create cufft plans.
semi-h
committed
9d36b6b
Add postprocessing in spectral space.
semi-h
committed
bcbaf18

Commits on Feb 20, 2024

Commits on Feb 29, 2024

Add a Poisson solver selection mechanism.
semi-h
committed
257b4ee
Implement a dedicated class for the FFT based Poisson solver.
semi-h
committed
87f0763
Store compact scheme coefficients in tdsops class.
semi-h
committed
dc3fe5c
Store domain length and dx in dirps_t.
semi-h
committed
7d67a0e
Change the way we initialise poisson_fft class.
semi-h
committed
bbe17d6
Construct the waves array for all periodic BCs.
semi-h
committed
266c3a9
Create cufft plans.
semi-h
committed
b1a3589
Add postprocessing in spectral space.
semi-h
committed
b138b40
Store plans at poisson_fft class level.
semi-h
committed
4953d87
Add forward FFT transforms.
semi-h
committed
4cde72d
Add backward FFT transforms.
semi-h
committed
178c920
Fix allocation size of the waves array on GPU.
semi-h
committed
7893007
Change shape of the c_x/y/z arrays.
semi-h
committed
de51f15
Add postprocessing kernel call.
semi-h
committed
994e83e
Add complex reordering functions.
semi-h
committed
f609516
Call complex reordering functions in forward and backward FFTs.
semi-h
committed
1f5564a
No need to initialise backends with dirps any more.
semi-h
committed
add29cc
Cleanups.
semi-h
committed
b982e52

Commits on Mar 1, 2024

Commits on Mar 4, 2024

Merge branch 'poisson' into 44-poisson-solver-iterative-solver
pbartholomew08
committed
4dcfc0f
Add configure option to select FFT or Iterative Poisson solver
pbartholomew08
committed
6e422d4
Find and link PETSc
pbartholomew08
committed
9a2756e

Commits on Mar 5, 2024

Commits on Mar 7, 2024

Merge branch 'main' into 44-poisson-solver-iterative-solver
pbartholomew08
committed
bc55cab

Commits on Mar 8, 2024

Commits on Mar 9, 2024

Move the PETSc definitions to a PETSc-specific submodule
pbartholomew08
committed
2b723cc

Commits on Mar 18, 2024

Commits on Mar 22, 2024

Merge branch 'main' into 44-poisson-solver-iterative-solver
pbartholomew08
committed
e556fdb

Commits on Apr 2, 2024

Merge branch 'main' into 44-poisson-solver-iterative-solver
pbartholomew08
committed
f8e80be
Indent consistently with fprettify
pbartholomew08
committed
332cf8e

Commits on Apr 3, 2024

Add PETSc->x3d2->PETSc operations
pbartholomew08
committed
161f8d3
Remove out of date comment
pbartholomew08
committed
966782d

Commits on Apr 30, 2024

Merge branch 'main' into 44-poisson-solver-iterative-solver
pbartholomew08
committed
f91b436

Commits on May 3, 2024

Implement test for Laplacian operator
pbartholomew08
committed
5a6a384

Commits on May 13, 2024

Make the CG/Laplacian test serial for now
pbartholomew08
committed
cdcfb82

Commits on Jun 21, 2024

Commits on Jun 24, 2024

Update PETSc/CG solver to use mesh object
pbartholomew08
committed
72b5400
Cleanup PETSc/Poisson solver
pbartholomew08
committed
40c7f15
Correct typo in CG test
pbartholomew08
committed
1a5e1dd

Commits on Jun 26, 2024

Implement OMP scalar product
pbartholomew08
committed
b03e7da
Implement vecAdd
pbartholomew08
committed
c1ea615
Add error messages to mesh dims getters
pbartholomew08
committed
73b4372
Implement OMP sum_yintox
pbartholomew08
committed
5ab5c11
Implement OMP sum_zintox
pbartholomew08
committed
154f63d
Implement vecadd and scalar product in vectorised form
pbartholomew08
committed
cf02979
Make sure inner loop indices are private
pbartholomew08
committed
c007d31
Run fprettify
pbartholomew08
committed
8076ed7

Commits on Jun 27, 2024

Merge branch '93-implement-omp-backend-for-basic-math-operators' into 44-poisson-solver-iterative-solver
pbartholomew08
committed
cd8658e

Commits on Jul 3, 2024

Commits on Jul 8, 2024

Poisson CG test passes for varying field
pbartholomew08
committed
ccb78f7
Improve CG/Poisson test
pbartholomew08
committed
b1bd12f
Ensure test status is always set in Poisson/CG Laplacian test
pbartholomew08
committed
7c79fe0

Commits on Jul 21, 2024

Simplify the Poisson_CG interface
pbartholomew08
committed
389541d
Update the Poisson object
pbartholomew08
committed
a928c1f
Fix error in backend
pbartholomew08
committed
fea6d3e
Pass backend (target) through Poisson CG
pbartholomew08
committed
022015a
Allocate backend in CG test
pbartholomew08
committed
49fd0fa
Give up passing matrix context and use a module global
pbartholomew08
committed
96cbf32
Remove debug prints
pbartholomew08
committed
5b9300e
Refine matrix through test
pbartholomew08
committed
bbc4176
Rewrite Poisson CG to run in DIR_Z orientation
pbartholomew08
committed
80185d6
Add a 7-point star preconditioner (WIP)
pbartholomew08
committed
fd1b6e9
Analytical solution and RHS were back to front
pbartholomew08
committed
b2220fe

Commits on Jul 22, 2024

Commits on Jul 24, 2024

Set Null space on Poisson matrices
pbartholomew08
committed
60b2726

Commits on Jul 25, 2024

Merge branch 'main' into 44-poisson-solver-iterative-solver
pbartholomew08
committed
9852c2d

Commits on Jul 30, 2024

Update iterative Poisson solver to latest x3d2 and begin parallelisation
pbartholomew08
committed
24fe79d
Correct solver setup - use O6 operator, O2 preconditioner
pbartholomew08
committed
2c4340c
Revert to serial preconditioner for time being
pbartholomew08
committed
fa6fc6a
Clean up test
pbartholomew08
committed
141c12e
Integrate CG Poisson with solver -- serial only
pbartholomew08
committed
da72568
Nonzero initial guess goes faster
pbartholomew08
committed
d10704c

Commits on Aug 26, 2024

Implement halo->global map for PETSc/Poisson preconditioner
pbartholomew08
committed
706058c

Commits on Sep 2, 2024

Commits on Oct 21, 2024

Resolve parallel indexing issues for PETSc
pbartholomew08
committed
b671eef

Commits on Oct 28, 2024

Test the Poisson solver in parallel
pbartholomew08
committed
10f8833

Commits on Nov 5, 2024

Generalise the Poisson evaluation test to work with generic linear operators
pbartholomew08
committed
824c631

Commits on Nov 12, 2024

Correct preconditioner indices
pbartholomew08
committed
48580f7

Commits on Dec 2, 2024

Integrate preconditioner class and CG class
pbartholomew08
committed
a724a5d
Documenting the Poisson/CG code
pbartholomew08
committed
c54191e

Commits on Feb 10, 2025

Enable building the project without PETSc
pbartholomew08
committed
079d9a6
Silence unused compiler warnings
pbartholomew08
committed
4fadb10
Enable running the high-order Poisson operator tests without PETSc
pbartholomew08
committed
c0649b9
Merge branch 'main' into 44-poisson-solver-iterative-solver
pbartholomew08
committed
816bc70
Correct Laplacian implementation by setting BCs...
pbartholomew08
committed
47ac899
Apply fprettify
pbartholomew08
committed
b3116c9
Re-apply fprettify
pbartholomew08
committed
12f0b2a
Remove submodules (CG/dummy) build
pbartholomew08
committed
e24619e

Commits on Feb 11, 2025

Remove submodules (CG/petsc) build
pbartholomew08
committed
4295cdc

Commits on Feb 13, 2025

Commits on Feb 17, 2025

Commits on Feb 18, 2025

Commits on Feb 23, 2025

Implement 27-point preconditioner based on PETSc DMDA
pbartholomew08
committed
f135b17

Commits on Feb 27, 2025

27-point stencil application working in serial
pbartholomew08
committed
a850241
Fix DMDA vector interactions
pbartholomew08
committed
4178b2e

Commits on Mar 3, 2025

Correct pointer use in PETSc solver
pbartholomew08
committed
df157cf

Commits on May 9, 2025

Implement pressure Laplacian based on face-staggered operators
pbartholomew08
committed
6e43b19

Commits on Jun 9, 2025

Merge branch 'main' into 44-poisson-solver-iterative-solver
pbartholomew08
committed
3ad375a