From 0a65b142cfa7169bf6271ecfff57c86cd4cb0bf9 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Wed, 12 Jul 2017 10:58:21 +0100
Subject: [PATCH] BUG: snappyHexMesh: calculate surface normals if face
 splitting. Fixes #522.

---
 .../snappyHexMesh/snappyHexMeshDriver/snappySnapDriver.C  | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriver.C b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriver.C
index 39a637d84c..b0d7e48425 100644
--- a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriver.C
+++ b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriver.C
@@ -2750,13 +2750,13 @@ void Foam::snappySnapDriver::doSnap
             motionSmoother& meshMover = meshMoverPtr();
 
 
-            // Calculate displacement at every patch point. Insert into
-            // meshMover.
-            // Calculate displacement at every patch point
+            // Calculate displacement at every patch point if we need it:
+            // - if automatic near-surface detection
+            // - if face splitting active
             pointField nearestPoint;
             vectorField nearestNormal;
 
-            if (snapParams.detectNearSurfacesSnap())
+            if (snapParams.detectNearSurfacesSnap() || doSplit)
             {
                 nearestPoint.setSize(pp.nPoints(), vector::max);
                 nearestNormal.setSize(pp.nPoints(), Zero);