24 axpy(residuum, -1, b);
Global vector based on Eigen vector.
void axpy(PETScVector &y, double const a, PETScVector const &x)
void scale(PETScVector &x, double const a)
void matMult(PETScMatrix const &A, PETScVector const &x, PETScVector &y)
void matMultAdd(PETScMatrix const &A, PETScVector const &v1, PETScVector const &v2, PETScVector &v3)
GlobalVector computeResiduum(GlobalVector const &x, GlobalVector const &xdot, GlobalMatrix const &M, GlobalMatrix const &K, GlobalVector const &b)