Assorted MPI Fixes and updates of interpolation classes to account for MP.

joshkellyjak · joshkellyjak · commit fd206572e495 · 2025-12-09T16:24:46.000Z
diff --git a/Common/include/interface_interpolation/CInterpolator.hpp b/Common/include/interface_interpolation/CInterpolator.hpp
@@ -131,7 +131,7 @@ class CInterpolator {
    * \note Main method that derived classes must implement.
    * \param[in] config - Definition of the particular problem.
    */
-  virtual void SetTransferCoeff(const CConfig* const* config) = 0;
+  virtual void SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) = 0;
 
   /*!
    * \brief Print information about the interpolation.
diff --git a/Common/include/interface_interpolation/CIsoparametric.hpp b/Common/include/interface_interpolation/CIsoparametric.hpp
@@ -64,7 +64,7 @@ class CIsoparametric final : public CInterpolator {
    * \brief Set up transfer matrix defining relation between two meshes
    * \param[in] config - Definition of the particular problem.
    */
-  void SetTransferCoeff(const CConfig* const* config) override;
+  void SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) override;
 
   /*!
    * \brief Print information about the interpolation.
diff --git a/Common/include/interface_interpolation/CMirror.hpp b/Common/include/interface_interpolation/CMirror.hpp
@@ -54,5 +54,5 @@ class CMirror final : public CInterpolator {
    * \brief Set up transfer matrix defining relation between two meshes
    * \param[in] config - Definition of the particular problem.
    */
-  void SetTransferCoeff(const CConfig* const* config) override;
+  void SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) override;
 };
diff --git a/Common/include/interface_interpolation/CMixingPlane.hpp b/Common/include/interface_interpolation/CMixingPlane.hpp
@@ -39,7 +39,7 @@ class CMixingPlane final : public CInterpolator {
     CMixingPlane(CGeometry**** geometry_container, const CConfig* const* config, unsigned int iZone,
                                    unsigned int jZone);
 
-    void SetTransferCoeff(const CConfig* const* config) override;
+    void SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) override;
 
     /*!
      * \brief Write interpolation details to file.
diff --git a/Common/include/interface_interpolation/CNearestNeighbor.hpp b/Common/include/interface_interpolation/CNearestNeighbor.hpp
@@ -63,7 +63,7 @@ class CNearestNeighbor final : public CInterpolator {
    * \brief Set up transfer matrix defining relation between two meshes.
    * \param[in] config - Definition of the particular problem.
    */
-  void SetTransferCoeff(const CConfig* const* config) override;
+  void SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) override;
 
   /*!
    * \brief Print interpolation statistics.
diff --git a/Common/include/interface_interpolation/CRadialBasisFunction.hpp b/Common/include/interface_interpolation/CRadialBasisFunction.hpp
@@ -56,7 +56,7 @@ class CRadialBasisFunction final : public CInterpolator {
    * \brief Set up transfer matrix defining relation between two meshes
    * \param[in] config - Definition of the particular problem.
    */
-  void SetTransferCoeff(const CConfig* const* config) override;
+  void SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) override;
 
   /*!
    * \brief Print information about the interpolation.
diff --git a/Common/include/interface_interpolation/CSlidingMesh.hpp b/Common/include/interface_interpolation/CSlidingMesh.hpp
@@ -49,7 +49,7 @@ class CSlidingMesh final : public CInterpolator {
    * \brief Set up transfer matrix defining relation between two meshes
    * \param[in] config - Definition of the particular problem.
    */
-  void SetTransferCoeff(const CConfig* const* config) override;
+  void SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) override;
 
  private:
   /*!
diff --git a/Common/src/interface_interpolation/CIsoparametric.cpp b/Common/src/interface_interpolation/CIsoparametric.cpp
@@ -37,7 +37,7 @@ using namespace GeometryToolbox;
 CIsoparametric::CIsoparametric(CGeometry**** geometry_container, const CConfig* const* config, unsigned int iZone,
                                unsigned int jZone)
     : CInterpolator(geometry_container, config, iZone, jZone) {
-  SetTransferCoeff(config);
+  SetTransferCoeff(geometry_container, config);
 }
 
 void CIsoparametric::PrintStatistics() const {
@@ -46,7 +46,7 @@ void CIsoparametric::PrintStatistics() const {
        << "  Interpolation clipped for " << ErrorCounter << " (" << ErrorRate << "%) target vertices." << endl;
 }
 
-void CIsoparametric::SetTransferCoeff(const CConfig* const* config) {
+void CIsoparametric::SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) {
   const su2double matchingVertexTol = 1e-12;  // 1um^2
 
   const int nProcessor = size;
diff --git a/Common/src/interface_interpolation/CMirror.cpp b/Common/src/interface_interpolation/CMirror.cpp
@@ -40,10 +40,10 @@ CMirror::CMirror(CGeometry**** geometry_container, const CConfig* const* config,
                        to_string(iZone) + string(" and ") + to_string(jZone) + string("."),
                    CURRENT_FUNCTION);
   }
-  SetTransferCoeff(config);
+  SetTransferCoeff(geometry_container, config);
 }
 
-void CMirror::SetTransferCoeff(const CConfig* const* config) {
+void CMirror::SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) {
   const int nProcessor = size;
 
   vector<unsigned long> allNumVertexTarget(nProcessor);
diff --git a/Common/src/interface_interpolation/CMixingPlane.cpp b/Common/src/interface_interpolation/CMixingPlane.cpp
@@ -33,16 +33,19 @@
 CMixingPlane::CMixingPlane(CGeometry**** geometry_container, const CConfig* const* config, unsigned int iZone,
                                    unsigned int jZone)
     : CInterpolator(geometry_container, config, iZone, jZone) {
-  SetTransferCoeff(config);
+  SetTransferCoeff(geometry_container, config);
 }
 
-void CMixingPlane::SetTransferCoeff(const CConfig* const* config) {
+void CMixingPlane::SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) {
     const auto nMarkerInt = config[donorZone]->GetnMarker_MixingPlaneInterface() / 2;
     const auto nDim = donor_geometry->GetnDim();
 
     const auto donor_config = config[donorZone];
     const auto target_config = config[targetZone];
 
+    const auto donor_geometry = geometry[donorZone][INST_0][MESH_0];
+    const auto target_geometry = geometry[targetZone][INST_0][MESH_0];
+
     //TODO turbo this approach only works if all the turboamchinery marker
     //    of all zones have the same amount of span wise sections.
     //TODO turbo initialization needed for the MPI routine should be place somewhere else.
@@ -52,14 +55,14 @@ void CMixingPlane::SetTransferCoeff(const CConfig* const* config) {
     targetSpans.resize(config[donorZone]->GetnMarker_MixingPlaneInterface());
 
     /*--- On the donor side ---*/
-    for (auto iMarkerInt = 0; iMarkerInt < nMarkerInt; iMarkerInt++){
+    for (auto iMarkerInt = 1; iMarkerInt < nMarkerInt + 1; iMarkerInt++){
         int markDonor = -1, markTarget = -1;
-        unsigned short donorFlag = 0, targetFlag = 0;
+        short donorFlag = 0, targetFlag = 0;
 
-        markDonor = donor_config->FindMixingPlaneInterfaceMarker(donor_config->GetnMarker_All());
+        markDonor = donor_config->FindMixingPlaneInterfaceMarker(donor_geometry->GetnMarker(), iMarkerInt);
         donorFlag = (markDonor != -1) ? donor_config->GetMarker_All_MixingPlaneInterface(markDonor) : -1;
 
-        markTarget = target_config->FindMixingPlaneInterfaceMarker(target_config->GetnMarker_All());
+        markTarget = target_config->FindMixingPlaneInterfaceMarker(target_geometry->GetnMarker(), iMarkerInt);
         targetFlag = (markTarget != -1) ? target_config->GetMarker_All_MixingPlaneInterface(markTarget) : -1;
 
 #ifdef HAVE_MPI
@@ -85,11 +88,17 @@ void CMixingPlane::SetTransferCoeff(const CConfig* const* config) {
     targetFlag= -1;
 
     for (int iSize=0; iSize<size; iSize++) {
-        if(buffMarkerDonor[iSize] >= 0.0) {
+        if(buffMarkerDonor[iSize] != -1) {
             markDonor = buffMarkerDonor[iSize];
             donorFlag = buffDonorFlag[iSize];
-            markTarget = buffMarkerDonor[iSize];
-            targetFlag = buffDonorFlag[iSize];
+            break;
+        }
+    }
+
+    for (int iSize=0; iSize<size; iSize++) {
+        if(buffMarkerTarget[iSize] != -1) {
+            markTarget = buffMarkerTarget[iSize];
+            targetFlag = buffTargetFlag[iSize];
             break;
         }
     }
@@ -98,7 +107,6 @@ void CMixingPlane::SetTransferCoeff(const CConfig* const* config) {
     delete [] buffMarkerTarget;
     delete [] buffTargetFlag;
 #endif
-
         if (markTarget == -1 || markDonor == -1) continue;
 
         nSpanDonor = donor_config->GetnSpanWiseSections();
@@ -115,9 +123,9 @@ void CMixingPlane::SetTransferCoeff(const CConfig* const* config) {
         if (nDim > 2) {
             targetSpans[iMarkerInt][nSpanTarget-1].donorSpan = nSpanDonor-1;
             targetSpans[iMarkerInt][nSpanTarget-1].coefficient = 0.0;
-            targetSpans[iMarkerInt][nSpanTarget].donorSpan = nSpanDonor;
-            targetSpans[iMarkerInt][nSpanTarget].coefficient = 0.0;
         }
+        targetSpans[iMarkerInt][nSpanTarget].donorSpan = nSpanDonor;
+        targetSpans[iMarkerInt][nSpanTarget].coefficient = 0.0;
 
         for(auto iSpanTarget = 1; iSpanTarget < nSpanTarget - 1; iSpanTarget++){
             auto &targetSpan = targetSpans[iMarkerInt][iSpanTarget];
@@ -147,31 +155,46 @@ void CMixingPlane::SetTransferCoeff(const CConfig* const* config) {
                     break;
                     
                 case LINEAR_INTERPOLATION:
-                    // Find the donor span interval that brackets the target span
-                    for (auto iSpanDonor = 1; iSpanDonor < nSpanDonor - 1; iSpanDonor++) {
-                        const auto test = abs(spanValuesTarget[iSpanTarget] - spanValuesDonor[iSpanDonor]);
-                        if(test < minDist && spanValuesTarget[iSpanTarget] > spanValuesDonor[iSpanDonor]){
-                            kSpan = iSpanDonor;
-                            minDist = test;
+                {
+                    bool printWarning = false;
+                    // Check if target span is within donor bound
+                    if (spanValuesTarget[iSpanTarget] <= spanValuesDonor[0]) {
+                        // Below hub - use hub value
+                        targetSpan.donorSpan = 0;
+                        targetSpan.coefficient = 0.0;
+                        printWarning = true;
+                    }
+                    else if (spanValuesTarget[iSpanTarget] >= spanValuesDonor[nSpanDonor - 1]) {
+                        // Above shroud - use shroud value
+                        targetSpan.donorSpan = nSpanDonor - 1;
+                        targetSpan.coefficient = 0.0;
+                        printWarning = true;
+                    }
+                    else {
+                        // Find the donor span interval that brackets the target span
+                        for (auto iSpanDonor = 0; iSpanDonor < nSpanDonor - 1; iSpanDonor++) {
+                            const auto test = abs(spanValuesTarget[iSpanTarget] - spanValuesDonor[iSpanDonor]);
+                            if(test < minDist && spanValuesTarget[iSpanTarget] > spanValuesDonor[iSpanDonor]){
+                                kSpan = iSpanDonor;
+                                minDist = test;
+                            }
                         }
+                        // Calculate interpolation coefficient
+                        coeff = (spanValuesTarget[iSpanTarget] - spanValuesDonor[kSpan]) / 
+                                (spanValuesDonor[kSpan + 1] - spanValuesDonor[kSpan]);
+                        targetSpan.donorSpan = kSpan;
+                        targetSpan.coefficient = coeff;
                     }
-                    // Calculate interpolation coefficient
-                    coeff = (spanValuesTarget[iSpanTarget] - spanValuesDonor[kSpan]) / 
-                            (spanValuesDonor[kSpan + 1] - spanValuesDonor[kSpan]);
-                    if ((coeff < 0) || (coeff > 1)) {
+                    if (printWarning) {
                         if (rank == MASTER_NODE) {
-                            cout << "Warning! Target spans exist outside the bounds of donor spans!" <<  endl;
-                            cout << "Target span " << iSpanTarget << " <- " << kSpan << " coeff = " << coeff << endl;
-                            if (iSpanTarget < nSpanTarget/2) cout << "This is likely an issue at the hub." << endl;
-                            if (iSpanTarget > nSpanTarget/2) cout << "This is likely an issue at the shroud." << endl;
-                            cout << "Setting coeff = 0.0" << endl;
+                            cout << "Warning! Target spans exist outside the bounds of donor spans! Clamping interpolator..." <<  endl;
+                            if (spanValuesTarget[iSpanTarget] <= spanValuesDonor[0]) cout << "This is an issue at the hub." << endl;
+                            if (spanValuesTarget[iSpanTarget] >= spanValuesDonor[nSpanDonor - 1]) cout << "This is an issue at the shroud." << endl;
+                            cout << "Setting coeff = 0.0 and transferring endwall value!" << endl;
                         }
-                        coeff = 0.0;
                     }
-                    targetSpan.donorSpan = kSpan;
-                    targetSpan.coefficient = coeff;
                     break;
-                    
+                }
                 default:
                     SU2_MPI::Error("MixingPlane interface option not implemented yet", CURRENT_FUNCTION);
                     break;
diff --git a/Common/src/interface_interpolation/CNearestNeighbor.cpp b/Common/src/interface_interpolation/CNearestNeighbor.cpp
@@ -33,15 +33,15 @@
 CNearestNeighbor::CNearestNeighbor(CGeometry**** geometry_container, const CConfig* const* config, unsigned int iZone,
                                    unsigned int jZone)
     : CInterpolator(geometry_container, config, iZone, jZone) {
-  SetTransferCoeff(config);
+  SetTransferCoeff(geometry_container, config);
 }
 
 void CNearestNeighbor::PrintStatistics() const {
   if (rank != MASTER_NODE) return;
   cout << "  Avg/max distance to closest donor point: " << AvgDistance << "/" << MaxDistance << endl;
 }
 
-void CNearestNeighbor::SetTransferCoeff(const CConfig* const* config) {
+void CNearestNeighbor::SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) {
   /*--- Desired number of donor points. ---*/
   const auto nDonor = max<unsigned long>(config[donorZone]->GetNumNearestNeighbors(), 1);
 
diff --git a/Common/src/interface_interpolation/CRadialBasisFunction.cpp b/Common/src/interface_interpolation/CRadialBasisFunction.cpp
@@ -47,7 +47,7 @@ extern "C" void dgemm_(const char*, const char*, const int*, const int*, const i
 CRadialBasisFunction::CRadialBasisFunction(CGeometry**** geometry_container, const CConfig* const* config,
                                            unsigned int iZone, unsigned int jZone)
     : CInterpolator(geometry_container, config, iZone, jZone) {
-  SetTransferCoeff(config);
+  SetTransferCoeff(geometry_container, config);
 }
 
 void CRadialBasisFunction::PrintStatistics() const {
@@ -98,7 +98,7 @@ su2double CRadialBasisFunction::Get_RadialBasisValue(RADIAL_BASIS type, const su
   return rbf;
 }
 
-void CRadialBasisFunction::SetTransferCoeff(const CConfig* const* config) {
+void CRadialBasisFunction::SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) {
   /*--- RBF options. ---*/
   const auto kindRBF = config[donorZone]->GetKindRadialBasisFunction();
   const bool usePolynomial = config[donorZone]->GetRadialBasisFunctionPolynomialOption();
diff --git a/Common/src/interface_interpolation/CSlidingMesh.cpp b/Common/src/interface_interpolation/CSlidingMesh.cpp
@@ -33,10 +33,10 @@
 CSlidingMesh::CSlidingMesh(CGeometry**** geometry_container, const CConfig* const* config, unsigned int iZone,
                            unsigned int jZone)
     : CInterpolator(geometry_container, config, iZone, jZone) {
-  SetTransferCoeff(config);
+  SetTransferCoeff(geometry_container, config);
 }
 
-void CSlidingMesh::SetTransferCoeff(const CConfig* const* config) {
+void CSlidingMesh::SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) {
   /* 0 - Variable declaration */
 
   /* --- General variables --- */
diff --git a/SU2_CFD/src/drivers/CMultizoneDriver.cpp b/SU2_CFD/src/drivers/CMultizoneDriver.cpp
@@ -270,7 +270,7 @@ void CMultizoneDriver::Preprocess(unsigned long TimeIter) {
     for (iZone = 0; iZone < nZone; iZone++) {
       for (unsigned short jZone = 0; jZone < nZone; jZone++){
         if(jZone != iZone && interpolator_container[iZone][jZone] != nullptr && prefixed_motion[iZone])
-          interpolator_container[iZone][jZone]->SetTransferCoeff(config_container);
+          interpolator_container[iZone][jZone]->SetTransferCoeff(geometry_container, config_container);
       }
     }
   }
diff --git a/SU2_CFD/src/interfaces/cfd/CMixingPlaneInterface.cpp b/SU2_CFD/src/interfaces/cfd/CMixingPlaneInterface.cpp
@@ -64,11 +64,11 @@ void CMixingPlaneInterface::BroadcastData_MixingPlane(const CInterpolator& inter
 
   /*--- Loop over interface markers. ---*/
 
-  for (auto iMarkerInt = 0u; iMarkerInt < donor_config->GetnMarker_MixingPlaneInterface()/2; iMarkerInt++) {
+  for (auto iMarkerInt = 1u; iMarkerInt < donor_config->GetnMarker_MixingPlaneInterface()/2 + 1; iMarkerInt++) {
 
     /*--- Find the markers containing the interface ---*/
-    short markDonor = donor_config->FindMixingPlaneInterfaceMarker(donor_geometry->GetnMarker());
-    short markTarget= target_config->FindMixingPlaneInterfaceMarker(target_geometry->GetnMarker());
+    short markDonor = donor_config->FindMixingPlaneInterfaceMarker(donor_geometry->GetnMarker(), iMarkerInt);
+    short markTarget= target_config->FindMixingPlaneInterfaceMarker(target_geometry->GetnMarker(), iMarkerInt);
 
     /*--- Check if this interface connects the two zones, if not continue. ---*/
     if(!CInterpolator::CheckInterfaceBoundary(markDonor, markTarget)) continue;
@@ -78,37 +78,38 @@ void CMixingPlaneInterface::BroadcastData_MixingPlane(const CInterpolator& inter
 
     /*--- Fill send buffers. ---*/
 
-    vector<short> sendDonorMarker(nSpanDonor);
-    vector<su2double> sendDonorVar(static_cast<size_t>(nSpanDonor) * nMixingVars);
+    vector<short> sendDonorMarker(nSpanDonor + 1);
+    vector<su2double> sendDonorVar(static_cast<size_t>(nSpanDonor + 1) * nMixingVars);
 
     if (markDonor != -1) {
-      for (auto iSpan = 0ul; iSpan < nSpanDonor; iSpan++) {
+      for (auto iSpan = 0ul; iSpan < nSpanDonor + 1; iSpan++) {
         GetDonor_Variable(donor_solution, donor_geometry, donor_config, markDonor, iSpan, 0);
         for (auto iVar = 0u; iVar < nMixingVars; iVar++) sendDonorVar[iSpan * nMixingVars + iVar] = Donor_Variable[iVar];
         sendDonorMarker[iSpan] = markDonor;
       }
     }
 #ifdef HAVE_MPI
     /*--- Gather data. ---*/
-    const auto nTotalDonors = nSpanDonor * size; // Number of donor spans across all ranks
-    const auto nSpanDonorVars = nSpanDonor * nMixingVars; // Number of variables to be transferred on each rank
+    const auto nTotalDonors = (nSpanDonor + 1) * size; // Number of donor spans across all ranks
+    const auto nSpanDonorVars = (nSpanDonor + 1) * nMixingVars; // Number of variables to be transferred on each rank
     vector<short> buffDonorMarker(nTotalDonors);
     vector<su2double> buffDonorVar(static_cast<unsigned long>(nTotalDonors) * nMixingVars); // Total number of variables to be transferred on all ranks
 
-    SU2_MPI::Allgather(sendDonorMarker.data(), nSpanDonor, MPI_SHORT,
-                      buffDonorMarker.data(), nSpanDonor, MPI_SHORT,
+    SU2_MPI::Allgather(sendDonorMarker.data(), nSpanDonor + 1, MPI_SHORT,
+                      buffDonorMarker.data(), nSpanDonor + 1, MPI_SHORT,
                       SU2_MPI::GetComm());
 
     SU2_MPI::Allgather(sendDonorVar.data(), nSpanDonorVars, MPI_DOUBLE,
                       buffDonorVar.data(), nSpanDonorVars, MPI_DOUBLE,
                       SU2_MPI::GetComm());
 
     for (auto iSize = 0; iSize < size; iSize++){
-      if (buffDonorMarker[static_cast<size_t>(iSize) * static_cast<unsigned long>(nSpanDonor)] != -1) {
-        for (auto iSpan = 0ul; iSpan < nSpanDonor; iSpan++){
+      if (buffDonorMarker[static_cast<size_t>(iSize) * static_cast<unsigned long>(nSpanDonor + 1)] != -1) {
+        for (auto iSpan = 0ul; iSpan < nSpanDonor + 1; iSpan++){
           for (auto iVar = 0u; iVar < nMixingVars; iVar++) sendDonorVar[iSpan * nMixingVars + iVar] = buffDonorVar[static_cast<size_t>(iSize) * nSpanDonorVars + iSpan * nMixingVars + iVar];
         }
-        markDonor = buffDonorMarker[static_cast<unsigned long>(iSize) * static_cast<unsigned long>(nSpanDonor)];
+        markDonor = buffDonorMarker[static_cast<unsigned long>(iSize) * static_cast<unsigned long>(nSpanDonor + 1)];
+        break; // Avoid overwriting
       }
     }
 #endif

Original file line number	Diff line number	Diff line change
`@@ -37,7 +37,7 @@ using namespace GeometryToolbox;`
`37`	`37`	`CIsoparametric::CIsoparametric(CGeometry**** geometry_container, const CConfig* const* config, unsigned int iZone,`
`38`	`38`	`unsigned int jZone)`
`39`	`39`	`: CInterpolator(geometry_container, config, iZone, jZone) {`
`40`		`- SetTransferCoeff(config);`
	`40`	`+ SetTransferCoeff(geometry_container, config);`
`41`	`41`	`}`
`42`	`42`
`43`	`43`	`void CIsoparametric::PrintStatistics() const {`
`@@ -46,7 +46,7 @@ void CIsoparametric::PrintStatistics() const {`
`46`	`46`	`<< " Interpolation clipped for " << ErrorCounter << " (" << ErrorRate << "%) target vertices." << endl;`
`47`	`47`	`}`
`48`	`48`
`49`		`-void CIsoparametric::SetTransferCoeff(const CConfig* const* config) {`
	`49`	`+void CIsoparametric::SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) {`
`50`	`50`	`const su2double matchingVertexTol = 1e-12; // 1um^2`
`51`	`51`
`52`	`52`	`const int nProcessor = size;`
Original file line number	Diff line number	Diff line change
`@@ -40,10 +40,10 @@ CMirror::CMirror(CGeometry**** geometry_container, const CConfig* const* config,`
`40`	`40`	`to_string(iZone) + string(" and ") + to_string(jZone) + string("."),`
`41`	`41`	`CURRENT_FUNCTION);`
`42`	`42`	`}`
`43`		`- SetTransferCoeff(config);`
	`43`	`+ SetTransferCoeff(geometry_container, config);`
`44`	`44`	`}`
`45`	`45`
`46`		`-void CMirror::SetTransferCoeff(const CConfig* const* config) {`
	`46`	`+void CMirror::SetTransferCoeff(CGeometry**** geometry, const CConfig* const* config) {`
`47`	`47`	`const int nProcessor = size;`
`48`	`48`
`49`	`49`	`vector<unsigned long> allNumVertexTarget(nProcessor);`