I am having trouble using two or more solvers from the AMGX library simultaneously. Doing this leads to undefined behavior of one of the solvers, i.e. the reported convergence rate are considerably different (often diverging) even though matrices start guess, RHS and solver configuration are identical. It seems that some objects are shared between the solvers even though they have separate handle objects (resources, config, solve, matrix, vectors). In particular “resources_destroy” cannot be called without error on any of the resoruces objects active simultaneously. I am running on a single GPU system.
Do you know the cause of the issue? Or do you have a working example of simultaneous use of several AMGX solvers?