Possible bug fix for points in intersections.

wyldckat · Oct 6, 2013 · 42c3b96 · 42c3b96
1 parent 4ff24f5
commit 42c3b96
Showing 1 changed file with 19 additions and 1 deletion.
diff --git a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoSnapDriverFeature.C b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoSnapDriverFeature.C
@@ -1660,8 +1660,26 @@ void Foam::autoSnapDriver::determineFeatures
                 // Mark point on the nearest feature edge. Note that we
                 // only search within the surrounding since the plane
                 // reconstruction might find a feature where there isn't one.
-                const point estimatedPt(pt + patchAttraction[pointI]);
+                point estimatedPt(pt + patchAttraction[pointI]);
 
+                if (multiRegionFeatureSnap)
+                {
+                    pointIndexHit multiPatchPt
+                    (
+                        findMultiPatchPoint
+                        (
+                            pt,
+                            pointFacePatchID[pointI],
+                            pointFaceCentres[pointI]
+                        )
+                    );
+                    if (multiPatchPt.hit())
+                    {
+                      estimatedPt = (estimatedPt+multiPatchPt.hitPoint())/scalar(2.0);
+                    }
+                }
+
+
                 // Determine nearest point on feature edge. Store constraint
                 // (calculated from feature edge, alternative would be to
                 //  use constraint calculated from both surfaceNormals)