diff --git a/benchmarks/levelset/cube.geo b/benchmarks/levelset/cube.geo
index b69a57bbde2b4f9802c37e3123e087fffaf1e3d7..45c1076681f5f53c79436abd4f5695d24791aff7 100644
--- a/benchmarks/levelset/cube.geo
+++ b/benchmarks/levelset/cube.geo
@@ -29,4 +29,8 @@ Levelset Point(13)={{0.5,0,0},{0.5,0.5,0.},{0.5,1,0},{0.5,0,0.5},{0.5,0,0.8},{0.
 
 Levelset CutMeshTri {5};
 
+Print "cube_cut.msh";
+
+
+
 
diff --git a/benchmarks/levelset/square.geo b/benchmarks/levelset/square.geo
index e51bf058f167ce1adbde9ae91fbe232dbab8bc55..413dc38e2d855defb1a20ebe1cb56b9aa89bacb8 100644
--- a/benchmarks/levelset/square.geo
+++ b/benchmarks/levelset/square.geo
@@ -29,6 +29,7 @@ Levelset Plane (1) = {0,1,0,-7};
 
 Levelset CutMeshTri {1};
 
+Print "square_cut.msh";
 
 
 
diff --git a/benchmarks/levelset/squareP.geo b/benchmarks/levelset/squareP.geo
new file mode 100644
index 0000000000000000000000000000000000000000..aab74d4b9c3c42365f403e36213483b1c25b561a
--- /dev/null
+++ b/benchmarks/levelset/squareP.geo
@@ -0,0 +1,33 @@
+Mesh.Algorithm=1;
+Mesh.CharacteristicLengthFactor=1.5;
+
+lc=0.1;
+Point(1) = {0.1, 0.1, 0.1}; //,lc};
+Point(2) = {0.1, 10, 0.1}; //,lc};
+Point(3) = {10, 10, 0.1}; //,lc};
+Point(4) = {10, 0.1, 0.1}; //,lc};
+Line(1) = {2, 3};
+Line(2) = {3, 4};
+Line(3) = {4, 1};
+Line(4) = {1, 2};
+
+Line Loop(5) = {1, 2, 3, 4};
+Plane Surface(6) = {5};
+
+Physical Line("top") = {1};
+Physical Line("bottom") = {3};
+Physical Line("left") = {4};
+Physical Line("right") = {2};
+Physical Surface("inside") = {6};
+
+Mesh 2;
+Levelset Plane (1) = {0,1,0,-7};
+//Levelset Point (1) = {{0.1, 2, 0},{3,2,0},{9,2,0},{5,2,0}, {0.1, 2.2, -1},{3,2.5,-1},{9,2,-1},{5,2,-1}, {0.2, 2, 1},{3,2,1},{9,2,1},{5,2,1}}; 
+
+Levelset CutMeshTri {1};
+
+Print "squareP_cut.msh";
+
+
+
+