From c86acdf9928fc0ad59d63aa20a6300510f0c4864 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 16 Nov 2011 10:00:37 +0100
Subject: [PATCH] The linear update is negative!
---
src/one-body-sample.cc | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index db4fe8ec..fd1f4bf0 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -233,10 +233,10 @@ int main(int argc, char *argv[]) {
b2 = b1;
b3 = b1;
- // b += linear update
- stiffnessMatrix.umv(u1_diff_old, b1);
- stiffnessMatrix.umv(u2_diff_old, b2);
- stiffnessMatrix.umv(u3_diff_old, b3);
+ // b -= linear update
+ stiffnessMatrix.mmv(u1_diff_old, b1);
+ stiffnessMatrix.mmv(u2_diff_old, b2);
+ stiffnessMatrix.mmv(u3_diff_old, b3);
// {{{ Assemble terms for the nonlinearity
--
GitLab