Hi, I’m a student studying computational physics. I’m now solving very large sparse linear system on GPU, and those problmes have very weak diagonal dominance, so old solvers for double precision arithmetics, such as BiCGSTAB do not work for it; solution diverges. So I’ve tried to find more attractive solvers and I found Algebraic Multigrid method (AMG) and the amgX library CUDA supports. But I cannot find how to bind the library into the Fortran language. Can anybody help me wrap the library?