From 07fd5db588937ff300494e1971c2fef13e080df5 Mon Sep 17 00:00:00 2001
From: Emilie Marchandise <emilie.marchandise@uclouvain.be>
Date: Thu, 10 Nov 2011 15:38:45 +0000
Subject: [PATCH]
---
benchmarks/3d/CubeAniso.geo | 56 ++++++++++++++++++++++++++++++++++++
benchmarks/3d/CubeAniso2.geo | 26 +++++++++++++++++
2 files changed, 82 insertions(+)
create mode 100644 benchmarks/3d/CubeAniso.geo
create mode 100644 benchmarks/3d/CubeAniso2.geo
diff --git a/benchmarks/3d/CubeAniso.geo b/benchmarks/3d/CubeAniso.geo
new file mode 100644
index 0000000000..69e324e12d
--- /dev/null
+++ b/benchmarks/3d/CubeAniso.geo
@@ -0,0 +1,56 @@
+Mesh.Algorithm = 7; //Bamg
+Mesh.Algorithm3D= 7; //mmg3d
+lc = 0.1;
+
+Point(1) = {0.0,0.0,0.0,lc};
+Point(2) = {1,0.0,0.0,lc};
+Point(3) = {1,1,0.0,lc};
+Point(4) = {0,1,0.0,lc};
+Line(1) = {4,3};
+Line(2) = {3,2};
+Line(3) = {2,1};
+Line(4) = {1,4};
+Line Loop(5) = {2,3,4,1};
+Plane Surface(6) = {5};
+Extrude Surface { 6, {0,0.0,1} };
+
+/*
+Field[1] = Attractor;
+Field[1].FacesList = {23};
+Field[1].NNodesByEdge = 100;
+
+Field[2] = Threshold;
+Field[2].IField = 1;
+Field[2].DistMax = 0.3;
+Field[2].DistMin = 0.05;
+Field[2].LcMax = 0.1;
+Field[2].LcMin = 0.01;
+Background Field = 2;
+
+Field[3] = BoundaryLayer;
+Field[3].IField = 1;
+Field[3].hfar = 0.1;
+Field[3].hwall_n = 0.005;
+Field[3].hwall_t = 0.1;
+Background Field = 3;
+*/
+
+/*
+Field[1] = BoundaryLayer;
+Field[1].FaceList = {23};
+Field[1].hfar = 0.01;
+Field[1].hwall_n = 0.005;
+Field[1].hwall_t = 0.1;
+Field[1].ratio = 1.0;
+Field[1].thickness = 0.5;
+Background Field = 1;
+*/
+
+Field[2] = MathEvalAniso;
+Field[2].m11 = "1./(0.1)^2";
+Field[2].m12 = "0";
+Field[2].m13 = "0";
+Field[2].m22 = "y+1/(0.02)^2";
+Field[2].m23 = "0";
+Field[2].m33 = "1/(0.1)^2";
+Background Field = 2;
diff --git a/benchmarks/3d/CubeAniso2.geo b/benchmarks/3d/CubeAniso2.geo
new file mode 100644
index 0000000000..a46a3390df
--- /dev/null
+++ b/benchmarks/3d/CubeAniso2.geo
@@ -0,0 +1,26 @@
+lc = 0.1;
+
+Mesh.RemeshParametrization=1; //(0) harmonic (1) conformal
+Mesh.RemeshAlgorithm=1; //(0) nosplit (1) automatic (2) split only with metis ///Default: 1
+
+Merge "CubeAniso.msh";
+CreateTopology;
+
+Compound Surface(106) = {6};
+Compound Surface(115) = {15};
+Compound Surface(119) = {19};
+Compound Surface(123) = {23};
+Compound Surface(127) = {27};
+Compound Surface(128) = {28};
+
+Surface Loop(100) = {106,115,119,123,127,128};
+Volume(1) = {100};
+
+Field[2] = MathEvalAniso;
+Field[2].m11 = "1./(0.1)^2";
+Field[2].m12 = "0";
+Field[2].m13 = "0";
+Field[2].m22 = "y+1/(0.02)^2";
+Field[2].m23 = "0";
+Field[2].m33 = "1/(0.1)^2";
+Background Field = 2;
--
GitLab