From 3911aa76ec0bd4ca8589fae5888630dc2eafccc6 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 14 Feb 2013 16:07:51 +0000
Subject: [PATCH] Bugfix in dFdpdq
[[Imported from SVN: r9193]]
---
dune/gfe/hyperbolichalfspacepoint.hh | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/dune/gfe/hyperbolichalfspacepoint.hh b/dune/gfe/hyperbolichalfspacepoint.hh
index 307ad7f2..192236b0 100644
--- a/dune/gfe/hyperbolichalfspacepoint.hh
+++ b/dune/gfe/hyperbolichalfspacepoint.hh
@@ -107,7 +107,9 @@ class HyperbolicHalfspacePoint
} else if (i==N-1 and j==N-1) {
- dFdpdq[i][j] = -1/(p[N-1]*p[N-1]*q[N-1]) - (p[N-1]-q[N-1]) / (p[N-1]*q[N-1]*q[N-1]) + diffNormSquared / (p[N-1]*p[N-1]*q[N-1]*q[N-1]);
+ dFdpdq[i][j] = -1/(p[N-1]*q[N-1])
+ - (p[N-1]-q[N-1]) / (p[N-1]*q[N-1]*q[N-1])
+ + (p[N-1]-q[N-1]) / (p[N-1]*p[N-1]*q[N-1]) + diffNormSquared / (2*p[N-1]*p[N-1]*q[N-1]*q[N-1]);
}
--
GitLab