Skip to content
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

Jump to bottom

Implement an iterative Poisson solver #69

Draft
wants to merge 99 commits into
base: main
Choose a base branch
from
Draft

Implement an iterative Poisson solver #69

wants to merge 99 commits into from

Conversation

pbartholomew08
Copy link
Member

This is based off the Poisson branch #32 and adds the interface to PETSc as an alternative Poisson solver.

semi-h and others added 30 commits February 19, 2024 11:46
@semi-h
Add a Poisson solver selection mechanism.
7a342be
@semi-h
Implement a dedicated class for the FFT based Poisson solver.
bcd1056
@semi-h
Store compact scheme coefficients in tdsops class.
e6c488b
@semi-h
Store domain length and dx in dirps_t.
d4c6831
@semi-h
Change the way we initialise poisson_fft class.
d0d11d9
@semi-h
Construct the waves array for all periodic BCs.
9cbdf95
@semi-h
Create cufft plans.
9d36b6b
@semi-h
Add postprocessing in spectral space.
bcbaf18
@pbartholomew08
Merge branch 'poisson' into 44-poisson-solver-iterative-solver
2b648d2
@pbartholomew08
Ignore emacs project configuration file
11a6960
@pbartholomew08
Merge branch 'main' into 44-poisson-solver-iterative-solver
fad59f1
@semi-h
Add a Poisson solver selection mechanism.
257b4ee
@semi-h
Implement a dedicated class for the FFT based Poisson solver.
87f0763
@semi-h
Store compact scheme coefficients in tdsops class.
dc3fe5c
@semi-h
Store domain length and dx in dirps_t.
7d67a0e
@semi-h
Change the way we initialise poisson_fft class.
bbe17d6
@semi-h
Construct the waves array for all periodic BCs.
266c3a9
@semi-h
Create cufft plans.
b1a3589
@semi-h
Add postprocessing in spectral space.
b138b40
@semi-h
Store plans at poisson_fft class level.
4953d87
@semi-h
Add forward FFT transforms.
4cde72d
@semi-h
Add backward FFT transforms.
178c920
@semi-h
Fix allocation size of the waves array on GPU.
7893007
@semi-h
Change shape of the c_x/y/z arrays.
de51f15
@semi-h
Add postprocessing kernel call.
994e83e
@semi-h
Add complex reordering functions.
f609516
@semi-h
Call complex reordering functions in forward and backward FFTs.
1f5564a
@semi-h
No need to initialise backends with dirps any more.
add29cc
@semi-h
Cleanups.
b982e52
@semi-h
Cleanup and style fix.
8e1f39b
pbartholomew08 and others added 21 commits July 3, 2024 14:20
@pbartholomew08
Merge branch 'main' into 44-poisson-solver-iterative-solver
842fecb
@pbartholomew08
Rewrite Poisson/CG test to support refinement
72e43f8
@pbartholomew08
Poisson CG test passes for varying field
ccb78f7 
Note the tolerance was raised to 1.0e-6
@pbartholomew08
Improve CG/Poisson test
b1bd12f 
Check that constant field solution is zero, check that varying
(periodic) field converges with 6th order
@pbartholomew08
Ensure test status is always set in Poisson/CG Laplacian test
7c79fe0
@pbartholomew08
Simplify the Poisson_CG interface
389541d 
The mesh object is available via the backend object
@pbartholomew08
Update the Poisson object
a928c1f
@pbartholomew08
Fix error in backend
fea6d3e
@pbartholomew08
Pass backend (target) through Poisson CG
022015a
@pbartholomew08
Allocate backend in CG test
49fd0fa
@pbartholomew08
Give up passing matrix context and use a module global
96cbf32 
This is the pragmatic solution, but it SUCKS
@pbartholomew08
Remove debug prints
5b9300e
@pbartholomew08
Refine matrix through test
bbc4176
@pbartholomew08
Rewrite Poisson CG to run in DIR_Z orientation
80185d6
@pbartholomew08
Add a 7-point star preconditioner (WIP)
fd1b6e9
@pbartholomew08
Analytical solution and RHS were back to front
b2220fe
@pbartholomew08
Correct order of setting mat options
2bef93b
@pbartholomew08
Check if PETSc is initialised (for convergence testing, etc)
8dcbfad
@pbartholomew08
Finalise PETSc in test
f3824e8
@pbartholomew08
Initialise PETSc properly
fb84bc9
@pbartholomew08
Set Null space on Poisson matrices
60b2726
@pbartholomew08 pbartholomew08 force-pushed the 44-poisson-solver-iterative-solver branch from 7223542 to 60b2726 Compare July 24, 2024 08:58
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Reviewers
No reviews
Assignees

@pbartholomew08 pbartholomew08

Labels
None yet
Projects
None yet
Milestone
No milestone
Development

Successfully merging this pull request may close these issues.

Poisson solver (iterative solver)
2 participants
@pbartholomew08 @semi-h