@@ -331,6 +331,28 @@ void CRadialBasisFunctionInterpolation::SetInternalNodes(CGeometry* geometry, CC
331331 }
332332 }
333333
334+ /* --- Adding nodes on markers considered as internal nodes ---*/
335+ for (auto iMarker = 0u ; iMarker < geometry->GetnMarker (); iMarker++){
336+
337+ /* --- Check if marker is considered as internal nodes ---*/
338+ if (config->GetMarker_All_Deform_Mesh_Internal (iMarker)){
339+
340+ /* --- Loop over marker vertices ---*/
341+ for (auto iVertex = 0ul ; iVertex < geometry->nVertex [iMarker]; iVertex++) {
342+
343+ /* --- Local node index ---*/
344+ auto iNode = geometry->vertex [iMarker][iVertex]->GetNode ();
345+
346+ /* --- if not among the boundary nodes ---*/
347+ if (find_if (BoundNodes.begin (), BoundNodes.end (), [&](CRadialBasisFunctionNode* i){return i->GetIndex () == iNode;}) == BoundNodes.end ()) {
348+ internalNodes.push_back (iNode);
349+ }
350+ }
351+ }
352+ }
353+
354+
355+
334356 /* --- In case of a parallel computation, the nodes on the send/receive markers are included as internal nodes
335357 if they are not already a boundary node with known deformation ---*/
336358
@@ -347,7 +369,7 @@ void CRadialBasisFunctionInterpolation::SetInternalNodes(CGeometry* geometry, CC
347369 /* --- Local node index ---*/
348370 auto iNode = geometry->vertex [iMarker][iVertex]->GetNode ();
349371
350- // / *--- if not among the boundary nodes ---*/
372+ /* --- if not among the boundary nodes ---*/
351373 if (find_if (BoundNodes.begin (), BoundNodes.end (), [&](CRadialBasisFunctionNode* i){return i->GetIndex () == iNode;}) == BoundNodes.end ()) {
352374 internalNodes.push_back (iNode);
353375 }
0 commit comments