From 528c78663e709a5c63a47d52fc698353824903a3 Mon Sep 17 00:00:00 2001
From: Jean-Francois Remacle <jean-francois.remacle@uclouvain.be>
Date: Wed, 23 May 2012 07:06:07 +0000
Subject: [PATCH] hopla
---
Mesh/CenterlineField.cpp | 28 +++++++++++++++++++++-------
1 file changed, 21 insertions(+), 7 deletions(-)
diff --git a/Mesh/CenterlineField.cpp b/Mesh/CenterlineField.cpp
index 056a7d3583..a2dc1dd1a3 100644
--- a/Mesh/CenterlineField.cpp
+++ b/Mesh/CenterlineField.cpp
@@ -1184,20 +1184,34 @@ void Centerline::operator() (double x, double y, double z, SMetric3 &metr, GEnt
double lc_a = 3.5*lc;
double lc_n, lc_t;
- if (onTubularSurface){
- lc_n = lc_t = lc;
+ if ( onTubularSurface){
+ double e = radMax/5.;
+ double hn = e/50.;
+ lc_n = hn;
+ double crv = 1./radMax;
+ double oneOverD2 = .5/(lc*lc) * (1. + sqrt (1. + ( 4.*crv*crv*lc*lc*lc*lc/ (lc_n*lc_n*CTX::instance()->mesh.smoothRatio*CTX::instance()->mesh.smoothRatio))));
+
+ lc_n = lc_t = sqrt(1./oneOverD2);
}
else{
- double e = radMax/4.;
- double hn = e/10.;
+ /// thickness of the refined layer
+ double e = radMax/5.;
+ // small size
+ double hn = e/50.;
double rm = std::max(radMax-ds, radMax-e);
- lc_t = 2*M_PI*rm/nbPoints;
+
+ lc_t = lc;
+ lc_n = std::min(lc,ds*(1.3-1) + hn);
+
+ double crv = 1./radMax;
+ double oneOverD2 = .5/(lc*lc) * (1. + sqrt (1. + ( 4.*crv*crv*lc*lc*lc*lc/ (lc_n*lc_n*CTX::instance()->mesh.smoothRatio*CTX::instance()->mesh.smoothRatio))));
+ lc_t = sqrt(1./oneOverD2);
+
+
//lc_n = lc_t = lc;
//double ratio = 1.02; //1. + (lc_t-hn)/e;
//printf("ratio =%g \n", ratio);
//lc_n = ds*(ratio-1) + hn;
- if (ds < e) lc_n = hn;
- else lc_n = lc_t;
}
double lam_a = 1./(lc_a*lc_a);
double lam_n = 1./(lc_n*lc_n);
--
GitLab