Bifurcation diagram 2D reaction diffusion type system on GPU

Hi all,

This is my first time posting, so forgive me if I am not following the proper etiquette.

I am currently working on computing bifurcation diagrams of 2D reaction diffusion type models (e.g. Brusselator model) with zero flux boundary conditions. I can successfully compute diagrams, but when increasing the resolution of the domain the computations become quite slow. I already tried to optimize the code by using the analytical expression for the Jacobian matrix and using a shift invariant eigensolver, but this only gets me so far.

Because of this I am interested in performing these computations on the GPU. I tried to adapt my CPU code based on this tutorial 🟠 2d Swift-Hohenberg equation (non-local) on the GPU · Bifurcation Analysis in Julia. However, this is (in my opinion ;)) not a very simple example to start with.

Hence my question: are there any more resources available for computing bifurcation diagrams on the GPU? Any other tips, tricks or general advice is also greatly appreciated.

Thank you in advance!