changed: force field projections for mixed multi-patch

the nodal averaging code in SIMbase::project uses the global
node number and this breaks since since these are for all
bases, not just basis 1.

src/ASM/LR/ASMu2D.C

@@ -2464,9 +2464,6 @@ Field* ASMu2D::getProjectedField (const Vector& coefs) const
 
 Fields* ASMu2D::getProjectedFields (const Vector& coefs, size_t) const
 {
-  if (!this->separateProjectionBasis())
-    return nullptr;
-
   size_t ncmp = coefs.size() / this->getNoProjectionNodes();
   if (ncmp*this->getNoProjectionNodes() == coefs.size())
     return new LRSplineFields2D(this->getBasis(ASM::PROJECTION_BASIS),

src/ASM/LR/ASMu3D.C

@@ -2006,9 +2006,6 @@ Field* ASMu3D::getProjectedField (const Vector& coefs) const
 
 Fields* ASMu3D::getProjectedFields (const Vector& coefs, size_t) const
 {
-  if (!this->separateProjectionBasis())
-    return nullptr;
-
   size_t ncmp = coefs.size() / this->getNoProjectionNodes();
   if (ncmp*this->getNoProjectionNodes() == coefs.size())
     return new LRSplineFields3D(this->getBasis(ASM::PROJECTION_BASIS),coefs,ncmp);

src/SIM/SIMbase.C

@@ -2362,6 +2362,11 @@ bool SIMbase::evalSecondarySolution (Matrix& field, int pindx) const
 
 bool SIMbase::fieldProjections () const
 {
+  // Forced for mixed multi-patch models as the nodal averaging
+  // in SIMbase::project() fails.
+  if (myModel.size() > 1 && myModel[0]->getNoBasis() > 1)
+    return true;
+
   for (const ASMbase* pch : myModel)
     if (pch->separateProjectionBasis())
       return true;