From 008842b624939349828b22435498dabf94db1ca7 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sun, 3 Apr 2011 08:45:23 +0000
Subject: [PATCH] bugfix in
 thirdDerivativeOfDistanceSquaredWRTFirst1AndSecond2Argument.  Now it produces
 the same result as the fd approximation at least in for S^1

[[Imported from SVN: r7066]]
---
 dune/gfe/unitvector.hh | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/dune/gfe/unitvector.hh b/dune/gfe/unitvector.hh
index 202f4137..205a07de 100644
--- a/dune/gfe/unitvector.hh
+++ b/dune/gfe/unitvector.hh
@@ -257,7 +257,7 @@ public:
                  + secondDerivativeOfArcCosSquared(sp) * Tensor3<double,N,N,N>::product(PpPq,pProjected)
                  + secondDerivativeOfArcCosSquared(sp) * Tensor3<double,N,N,N>::product(PpPqp,Pq)
                  - secondDerivativeOfArcCosSquared(sp) * Tensor3<double,N,N,N>::product(qProjected,Pq_sp)
-                 - derivativeOfArcCosSquared(sp)       * Tensor3<double,N,N,N>::product(PpPq,qProjected);
+                 - derivativeOfArcCosSquared(sp)       * Tensor3<double,N,N,N>::product(qProjected,Pq);
                
         return result;
     }
-- 
GitLab