fix of beuler res

exanauts · Apr 19, 2018 · 28637d9 · 28637d9
1 parent 4e3055d
commit 28637d9
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/include/ad.hpp b/include/ad.hpp
@@ -489,7 +489,7 @@ template <class T> void integrate(pVector<T> &x) {
     adlinsolve<T>(J, y);
 #endif
     pVector<T> res(dim);
-    res = Jold * y + yold;
+    res = Jold * y - yold;
 #ifdef DBUG
     cout << "Norm(y): " << y.norm() << endl << " y: " << y << endl;
     cout << "|Ax - b| " << res.norm() << endl;