void MatMult(const FullMatrix& A, const dTensor1& invec, dTensor1& outvec) { for (int i=1; i<=A.get_NumRows(); i++) { int jmax = A.get_NZrow(i); double tmp = 0.0; for (int j=1; j<=jmax; j++) { int ind = A.get_Index(i,j); double val = A.get_Value(i,j); tmp = tmp + val*invec.get(ind); } outvec.set(i, tmp ); } }