Fixed: When the elasticity solver is mixed, use an additional MADOF array

defined on the first basis for explicit phase field. Don't use shared grid
if the phase field solver is used and the elasticity solver is mixed.

Author: akva2

SIMDynElasticity.h

@@ -336,6 +336,12 @@ class SIMDynElasticity : public Sim
     return result;
   }
 
+  //! \brief Adds an additional MADOF array if required.
+  bool preprocessB() override
+  {
+    return this->mixedProblem() ? this->addMADOF(1,1,false) : true;
+  }
+
 private:
   std::string energFile; //!< File name for global energy output
 

SIMExplPhaseField.C

@@ -79,8 +79,7 @@ bool SIMExplPhaseField::solveStep (TimeStep& tp, bool)
 bool SIMExplPhaseField::saveStep (const TimeStep& tp, int& nBlock)
 {
   if (tp.step%opt.saveInc == 0 && opt.format >= 0)
-    if (myOwner->writeGlvS1(phaseField,++myStep,nBlock,tp.time.t,
-                            "phase",6,1,true) < 0) return false;
+    return myOwner->writeGlvS(phaseField,"phase",++myStep,nBlock);
 
   return true;
 }

SIMFractureDynamics.h

@@ -61,7 +61,14 @@ class SIMFracture : public Coupling<SolidSolver,PhaseSolver>
   //! \brief Initializes and sets up field dependencies.
   virtual void setupDependencies()
   {
-    this->S1.registerDependency(&this->S2,"phasefield",1);
+    if (this->S1.mixedProblem() && this->S2.getNoParamDim() == 0)
+      // The solid solver is mixed and we use a mesh-less explicit phase field.
+      // Then set up the dependency to the phase field using an additional
+      // MADOF array defined for a scalar field on the first basis.
+      this->S1.registerDependency("phasefield",&this->S2);
+    else
+      this->S1.registerDependency(&this->S2,"phasefield");
+
     // The tensile energy is defined on integration points and not nodal points.
     // It is a global buffer array across all patches in the model.
     // Use an explicit call instead of normal couplings for this.

SIMPhaseField.h

@@ -54,7 +54,7 @@ template<class Dim> class SIMPhaseField : public Dim
   explicit SIMPhaseField(Dim* gridOwner = nullptr, size_t n = 2) : Dim(1)
   {
     Dim::myHeading = "Cahn-Hilliard solver";
-    if (gridOwner && gridOwner->createFEMmodel())
+    if (gridOwner && !gridOwner->mixedProblem() && gridOwner->createFEMmodel())
       this->clonePatches(gridOwner->getFEModel(),gridOwner->getGlob2LocMap());
 
     eps_d0 = refTol = 0.0;