corrections of NavierStokes_TH_Multiphase, NavierStokes_TaylorHodd and CahnHilliard

extensions/base_problems/CahnHilliard.cc

@@ -124,8 +124,12 @@ void CahnHilliard::solveInitialProblem(AdaptInfo *adaptInfo)
       HL_SignedDistTraverse reinit("reinit", getMesh()->getDim());
       reinit.calcSignedDistFct(adaptInfo, &tmp);
     
+#ifndef HAVE_PARALLEL_DOMAIN_AMDIS
       Recovery recovery(L2_NORM, 1);
       recovery.recoveryUh(&tmp, *prob->getSolution()->getDOFVector(0));
+#else
+      prob->getSolution()->getDOFVector(0)->interpol(tmp);
+#endif
     }
     
 

extensions/base_problems/NavierStokes_TH_MultiPhase.cc

@@ -46,8 +46,14 @@ NavierStokes_TH_MultiPhase::NavierStokes_TH_MultiPhase(const std::string &name_,
 
 NavierStokes_TH_MultiPhase::~NavierStokes_TH_MultiPhase()
 {
-  if (densityPhase) { delete densityPhase; densityPhase = NULL; }
-  if (viscosityPhase) { delete viscosityPhase; viscosityPhase = NULL; }
+  if (densityPhase) { 
+    delete densityPhase; 
+    densityPhase = NULL; 
+  }
+  if (viscosityPhase) { 
+    delete viscosityPhase; 
+    viscosityPhase = NULL; 
+  }
 }
 
 

extensions/base_problems/NavierStokes_TaylorHood.cc

@@ -44,8 +44,8 @@ NavierStokes_TaylorHood::NavierStokes_TaylorHood(const std::string &name_, bool
 
   // type of non-linear term: 
   //   0... u^old*grad(u_i^old), 
-  //   1... u'*grad(u_i^old), 
-  //   2... u^old*grad(u'_i)
+  //   1... u*grad(u_i^old), 
+  //   2... u^old*grad(u_i)
   //   3... (1)+(2)-(0)
   Initfile::get(self::name + "->non-linear term", nonLinTerm);