From ba9d868f51d72b8c390f777e10296f29f8a9d8c0 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 20 Aug 2007 09:38:49 +0000
Subject: [PATCH] damp rotation as well
[[Imported from SVN: r1539]]
---
dirneucoupling.cc | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index 9e6062cd..1146d64e 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -356,7 +356,8 @@ int main (int argc, char *argv[]) try
for (int j=0; j<dim; j++)
lambda.r[j] = (1-damping) * lambda.r[j] + damping * (referenceInterface.r[j] + averageInterface.r[j]);
- lambda.q = averageInterface.q.mult(referenceInterface.q);
+
+ lambda.q = Quaternion<double>::interpolate(lambda.q, averageInterface.q, damping);
// ////////////////////////////////////////////////////////////////////////
// Write the two iterates to disk for later convergence rate measurement
--
GitLab