The fully-vectorial MR-FEM-CSI algorithm has been implemented in C/C++ for distributed memory systems using MPI under the PETSc framework .
A major benefit of PETSc is that parallel storage and operations can be transparent to the implementation: once vectors such as [[chi].
and Aztec approaches, object-oriented design methodologies are at the basis of the Diffpack project [Bruaset and Langtangen 1997].
7: Iteratively solve the system using either the KSP AGMRES (Algorithm 1) or the PETSc
built-in KSP GMRES [41, Algorithm 4].
Knyazev, Preconditioned eigensolver LOBPCG in hypre and PETSc
, in Proceedings of the 16th International Conference on Domain Decomposition Methods 2005, O.
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
SMITH, Making automatic differentiation truly automatic: coupling PETSc
with ADIC, Future Generation Computer Systems, 21 (2005), pp.
TT to obtain the finite element discretization and PETSc
for the numerical linear algebra (routines for matrix-vector multiplication, scalar products, etc.