BOUT-dev
BOUT-dev copied to clipboard
Parallisation in x/z
Currently not working with FCI.
One option would be to figure out how many guard cells we need based on the parallisation and the grid.
Another option is openmp
And another is to use PETSc and implement the interpolation as a matrix routine
FCI does currently work with OpenMP.
The main issue with MPI and FCI is working out what processors the field lines end on and doing that communication. If the field lines only ever go to the processors immediately up/down or left/right, then it should be fine.
More flexible grids and meshes are really the key enabler here.