GPU acceleration in GetFEM

Hi everyone,

I would like to know if anyone has already successfully run GetFEM on a GPU, or integrated GetFEM with a GPU-based computing framework.

My goal is to accelerate computationally intensive FEM simulations by taking advantage of GPU parallelization. In particular, I am interested in knowing whether there is currently a way to offload parts of the GetFEM computation (assembly, linear algebra, nonlinear solver, etc.) to a GPU, for example using CUDA, HIP, PETSc, or another framework.

Has anyone already experimented with this, or developed a working implementation?

I would be very interested in hearing about your experience, including any practical approach, limitations, or recommendations on where to start.

Thanks in advance!