From 7fd1c722a8541da565472ab2c186b03151d22b67 Mon Sep 17 00:00:00 2001
From: Emilie Marchandise <emilie.marchandise@uclouvain.be>
Date: Thu, 27 Jan 2011 10:48:57 +0000
Subject: [PATCH]

---
 .../cylinder_double_boundary_layer.geo        | 37 +++++++++++++++++++
 .../extrude/sphere_double_boundary_layer.geo  |  3 ++
 benchmarks/extrude/u_shape_boundary_layer.geo |  4 +-
 3 files changed, 42 insertions(+), 2 deletions(-)
 create mode 100644 benchmarks/extrude/cylinder_double_boundary_layer.geo

diff --git a/benchmarks/extrude/cylinder_double_boundary_layer.geo b/benchmarks/extrude/cylinder_double_boundary_layer.geo
new file mode 100644
index 0000000000..8e67b5d038
--- /dev/null
+++ b/benchmarks/extrude/cylinder_double_boundary_layer.geo
@@ -0,0 +1,37 @@
+//-*- C++ -*-
+/*********************************************       
+Cube uniformly meshed along all three axes      
+**********************************************/       
+lc = 0.2; //.05;     
+// Radius of the cylinder:
+R = 0.5; 
+// Length of the cylinder:
+L = 5; 
+
+Point(0) = {0, 0, 0, lc}; 
+Point(1) = {R, 0, 0, lc}; 
+Point(2) = {0, R, 0, lc};
+Point(3) = {-R, 0, 0, lc};
+Point(4) = {0, -R, 0, lc};
+
+Circle(1) = {1,0,2};
+Circle(2) = {2,0,3};
+Circle(3) = {3,0,4};
+Circle(4) = {4,0,1};
+
+nb_pts_edge = 5;
+nb_layers = 30;
+
+For iLine In {1:4}
+	Transfinite Line{iLine} = nb_pts_edge;
+EndFor
+
+Line Loop(1) = {1,2,3,4};
+Plane Surface(1) = {1};
+
+out[] = Extrude{0, 0, L}{Surface{1}; Layers{nb_layers};};
+
+Extrude { Surface{13,17,21,25}; Layers{5, 0.1}; Recombine; Using Index[0]; }
+Extrude { Surface{-13,-17,-21,-25}; Layers{5, 0.1}; Recombine; Using Index[1]; }
+
+
diff --git a/benchmarks/extrude/sphere_double_boundary_layer.geo b/benchmarks/extrude/sphere_double_boundary_layer.geo
index e7f91dc106..268efa27d5 100644
--- a/benchmarks/extrude/sphere_double_boundary_layer.geo
+++ b/benchmarks/extrude/sphere_double_boundary_layer.geo
@@ -37,5 +37,8 @@ Ruled Surface(26) = {25};
 Line Loop(27) = {-4,12,-6};
 Ruled Surface(28) = {27};
 
+Surface Loop(1)={14,16,18,20,22,24,26,28};
+Volume(2) = {1};
+
 Extrude { Surface{14:28:2}; Layers{5, 0.1}; Recombine; Using Index[0]; }
 Extrude { Surface{-14,-16,-18,-20,-22,-24,-26,-28}; Layers{5, 0.1}; Recombine; Using Index[1]; }
diff --git a/benchmarks/extrude/u_shape_boundary_layer.geo b/benchmarks/extrude/u_shape_boundary_layer.geo
index 83ea224f35..1c157cf1f1 100644
--- a/benchmarks/extrude/u_shape_boundary_layer.geo
+++ b/benchmarks/extrude/u_shape_boundary_layer.geo
@@ -11,8 +11,8 @@ BSpline(1) = {5, 4, 3, 2, 1};
 BSpline(2) = {1, 9, 8, 7, 5};
 Line(3) = {5, 6};
 
-Extrude { Line{1,-3}; Layers{5,0.1}; Using Index[0]; }
-Extrude { Line{2,3}; Layers{5,0.1}; Using Index[1]; }
+Extrude { Line{1,-3}; Layers{5,0.2}; Using Index[0]; }
+Extrude { Line{2,3}; Layers{5,0.2}; Using Index[1]; }
 
 // fix leading edge by hand
 Coherence Point {25, 16};
-- 
GitLab