From 21f3e17aa112a5c1b6fc1f3e275968cb1439a40e Mon Sep 17 00:00:00 2001
From: Thomas Toulorge <thomas.toulorge@mines-paristech.fr>
Date: Fri, 14 Jun 2013 09:28:14 +0000
Subject: [PATCH] Fixed NACA00 level-set for non-zero y0
---
Geo/gmshLevelset.cpp | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/Geo/gmshLevelset.cpp b/Geo/gmshLevelset.cpp
index e553ede60d..53e6c0b404 100644
--- a/Geo/gmshLevelset.cpp
+++ b/Geo/gmshLevelset.cpp
@@ -1284,7 +1284,7 @@ void gLevelsetNACA00::getClosestBndPoint(double x, double y, double z,
if (xtb < tolr) xtb = tolr;
}
xb = _x0+xtb;
- yb = (y >= _y0) ? _y0+ytb : -(_y0+ytb);
+ yb = (y >= _y0) ? _y0+ytb : _y0-ytb;
const double norm = sqrt(1.+dyb*dyb);
curvRad = norm*norm*norm/fabs(ddyb);
}
--
GitLab