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!