24 axpy(x_dot, -1., x_prev);
25 scale(x_dot, 1. / dt);
26 matMult(M, x_dot, residuum);
27 matMultAdd(K, x, residuum, residuum);
28 axpy(residuum, -1., b);
Global vector based on Eigen vector.
GlobalVector computeResiduum(double const dt, GlobalVector const &x, GlobalVector const &x_prev, GlobalMatrix const &M, GlobalMatrix const &K, GlobalVector const &b)