From 16da62546b15062dadd9086f5c8ef565e6f30d7a Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Sat, 14 Dec 2013 13:43:17 +0100
Subject: [PATCH] [Cleanup]
---
src/one-body-sample.cc | 6 ++----
1 file changed, 2 insertions(+), 4 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 7d356f56..54b50bc3 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -295,9 +295,7 @@ int main(int argc, char *argv[]) {
Vector a_initial(fineVertexCount);
a_initial = 0.0;
{
- /* We solve Au + Cv + Ma + Psi(v) = ell, thus
- Ma = - (Au + Cv + Psi(v) - ell)
- */
+ // We solve Ma = ell - [Au + Cv + Psi(v)]
Vector accelerationRHS(fineVertexCount);
{
accelerationRHS = 0.0;
@@ -305,8 +303,8 @@ int main(int argc, char *argv[]) {
Arithmetic::addProduct(accelerationRHS, C, v_initial);
// NOTE: We assume differentiability of Psi at v0 here!
myGlobalNonlinearity->addGradient(v_initial, accelerationRHS);
- accelerationRHS -= ell;
accelerationRHS *= -1.0;
+ accelerationRHS += ell;
}
solveLinearProblem(noNodes, M, accelerationRHS, a_initial, MNorm,
parset.sub("a0.solver"));
--
GitLab