From 8c5f0a566e6536f0a2b8e14cd678446430a9721d Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 20 Nov 2007 18:24:59 +0000
Subject: [PATCH] bugfix: damped update of orientation part of lambda was wrong

[[Imported from SVN: r1793]]
---
 dirneucoupling.cc | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index 0dc6222b..af825e19 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -430,7 +430,9 @@ int main (int argc, char *argv[]) try
             lambda.r[j] = (1-dirichletDamping) * lambda.r[j] 
                 + dirichletDamping * (referenceInterface.r[j] + averageInterface.r[j]);
 
-        lambda.q = Quaternion<double>::interpolate(lambda.q, averageInterface.q, dirichletDamping);
+        lambda.q = Quaternion<double>::interpolate(lambda.q, 
+                                                   referenceInterface.q.mult(averageInterface.q), 
+                                                   dirichletDamping);
 
         std::cout << "Lambda: " << lambda << std::endl;
 
-- 
GitLab