From 9edaad6c2d13f75ae5fc261c7177ac451d8eac54 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 14 Feb 2013 16:07:48 +0000
Subject: [PATCH] Simplify formula for a third derivative
[[Imported from SVN: r9191]]
---
dune/gfe/hyperbolichalfspacepoint.hh | 5 +----
1 file changed, 1 insertion(+), 4 deletions(-)
diff --git a/dune/gfe/hyperbolichalfspacepoint.hh b/dune/gfe/hyperbolichalfspacepoint.hh
index 96ef67f4..aa77c3ce 100644
--- a/dune/gfe/hyperbolichalfspacepoint.hh
+++ b/dune/gfe/hyperbolichalfspacepoint.hh
@@ -387,10 +387,7 @@ public:
} else if (i==N-1 and j==N-1 and k==N-1) {
- dFdqdqdq[i][j][k] = -2.0/Dune::Power<3>::eval(q[N-1])
- - (2*p[N-1]*p[N-1]*q[N-1] - p[N-1]*q[N-1]*q[N-1]) / (p[N-1]*p[N-1]*Dune::Power<4>::eval(q[N-1]))
- + 2 * (p[N-1]-q[N-1]) / (p[N-1]*Dune::Power<3>::eval(q[N-1]))
- - 3 * diffNormSquared / (p[N-1]*Dune::Power<4>::eval(q[N-1]));
+ dFdqdqdq[i][j][k] = (2.0*p[N-1] -1 -4*p[N-1] -3*diffNormSquared)/(p[N-1]*Dune::Power<3>::eval(q[N-1]));
}
--
GitLab