diff --git a/benchmarks/3d/CubeAnisoBL.geo b/benchmarks/3d/CubeAnisoBL.geo
new file mode 100644
index 0000000000000000000000000000000000000000..40265e9da40720c9026ac2a1f0983cb70e63b4f1
--- /dev/null
+++ b/benchmarks/3d/CubeAnisoBL.geo
@@ -0,0 +1,55 @@
+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[1] = BoundaryLayer;
+Field[1].FacesList = {23};
+Field[1].hfar = 0.1;
+Field[1].hwall_n = 0.005;
+Field[1].hwall_t = 0.1;
+Field[1].thickness = 0.04;
+Background Field = 1;