Hello,
i want to solve a pde with periodic boundary conditions on multiple processors.
There is a good way to do it with Petsc, but i would like to move to julia.
Is there a way to do this with ParallelStencil / ImplicitGlobalGrid?
I could not figure out how to make this work. Any suggestions or alternative packages are welcome!
thanks
Hi @gerhardu and welcome!
solve a pde with periodic boundary conditions on multiple processors […] Is there a way to do this with ParallelStencil / ImplicitGlobalGrid?
Yes, using ParallelStencil and ImplicitGlobalGrid, periodic boundary condition implementation is straight forward.
From this example, just modify the init_global_grid() accordingly and you should have the 3D diffusion running periodic in x,y,z dims:
init_global_grid(nx, ny, nz; periodx=1, periody=1, periodz=1)
Note that for the CPU version of the code, modify the init_global_grid() in here.
Great thanks!