Hi @jdm204, nice to see your using ParallelStencil.jl in your work. The above snippet represents in this context an explicit update. However, depending on your needs, you could use various iterative-based solution approaches to achieve an implicit solution. Relying on direct-iterative approaches (“inverting a big matrix”) may be a bottleneck in 3D. Some hints about (GPU-friendly) iterative implicit approaches in this paper and this course.
2 Likes