diff --git a/benchmarks/boolean/aneurysm.stl b/benchmarks/boolean/aneurysm.stl
new file mode 100644
index 0000000000000000000000000000000000000000..4aedb8c65d6a0f6e930b75de870293d894cd68aa
Binary files /dev/null and b/benchmarks/boolean/aneurysm.stl differ
diff --git a/benchmarks/boolean/aneurysmBL.py b/benchmarks/boolean/aneurysmBL.py
new file mode 100644
index 0000000000000000000000000000000000000000..68b0b95745f1c5efc18e8a38e6a9dedb29b3f64e
--- /dev/null
+++ b/benchmarks/boolean/aneurysmBL.py
@@ -0,0 +1,22 @@
+from dgpy import *
+
+GmshSetOption('General', 'Verbosity', 99.0) 
+GmshSetOption('Mesh', 'SaveAll', 1.0) 
+GmshSetOption('Mesh', 'Algorithm', 1.0)
+
+g =  GModel()
+g.setFactory('Gmsh')
+
+g.load("aneurysm.stl");
+g.createTopologyFromMesh();
+
+face = g.getFaceByTag(1);
+newface  = g.extrudeBoundaryLayer(face, 4, 0.5, 0, -1)
+newface2 = g.extrudeBoundaryLayer(face, 4, -0.5, 1, -1)
+
+print "*** face = %d " % face.tag()
+print "*** new face = %d newface2 = %d " % (newface.tag(), newface2.tag())
+
+g.mesh(2)
+g.save("aneurysmBL.msh")
+