diff --git a/dirneucoupling.cc b/dirneucoupling.cc index d672df29f2acdf444312b596d4e5eb95fded9645..c0d9dcd8c4bfd3aea87a6c5306e6b7c0c1bb4105 100644 --- a/dirneucoupling.cc +++ b/dirneucoupling.cc @@ -356,8 +356,6 @@ int main (int argc, char *argv[]) try VectorType oldSolution3d = x3d; RodSolutionType oldSolutionRod = rodX; - RigidBodyMotion<3> averageInterface; - if (ddType=="FixedPointIteration") { // ////////////////////////////////////////////////// @@ -428,6 +426,8 @@ int main (int argc, char *argv[]) try // Extract new interface position and orientation // /////////////////////////////////////////////////////////// + RigidBodyMotion<3> averageInterface; + computeAverageInterface(interfaceBoundary[toplevel], x3d, averageInterface); //averageInterface.r[0] = averageInterface.r[1] = 0;