7: Iteratively solve the system using either the KSP AGMRES (Algorithm 1) or the
PETSc built-in KSP GMRES [41, Algorithm 4].
It is worth noting here that, despite all the different optimization strategies, sparse matrix products are still not standard in sparse libraries such as
PETSc [2].
SMITH, Making automatic differentiation truly automatic: coupling
PETSc with ADIC, Future Generation Computer Systems, 21 (2005), pp.