I would like to be able to easily switch between programs which use PETSc and the same programs which use Trilinos. The rationale for this would be to easily compare the two (see whether the solutions are the same, speed, etc.). Probably need to accomplish this with some sort of struct which takes in a template parameter and gives back a namespace or something.
Actually need to get rid of PETSC. Will hopefully used geometric multigrid in the future so that I can also get rid of Trilinos (might be some work to do that though).
I would like to be able to easily switch between programs which use PETSc and the same programs which use Trilinos. The rationale for this would be to easily compare the two (see whether the solutions are the same, speed, etc.). Probably need to accomplish this with some sort of struct which takes in a template parameter and gives back a namespace or something.