diff --git a/gendocs.sh b/gendocs.sh
index 8549cd4..2cc20b4 100755
--- a/gendocs.sh
+++ b/gendocs.sh
@@ -46,41 +46,41 @@ func_primitive="meshabox meshasphere meshanellip meshunitsphere meshacylinder
                 meshgrid5 meshgrid6 latticegrid extrudecurve meshcylinders"
 func_inquery="finddisconnsurf surfedge volface extractloops meshconn  
                 meshcentroid nodevolume elemvolume neighborelem layersurf 
-		faceneighbors edgeneighbors maxsurf flatsegment orderloopedge  
-		mesheuler bbxflatsegment surfplane surfinterior surfpart
+				faceneighbors edgeneighbors maxsurf flatsegment orderloopedge  
+				mesheuler bbxflatsegment surfplane surfinterior surfpart
                 surfseeds meshquality meshedge meshface surfacenorm nodesurfnorm
                 uniqedges uniqfaces advancefront innersurf outersurf surfvolume
-		insurface"
+				insurface"
 func_meshfix="meshcheckrepair meshreorient removedupelem 
                 removedupnodes removeisolatednode removeisolatedsurf
                 surfaceclean getintersecttri delendelem surfreorient"
 func_metch="proj2mesh dist2surf regpt2surf affinemap metchgui metchgui_one"
 func_remesh="meshresample remeshsurf smoothsurf sortmesh mergemesh 
                 meshrefine mergesurf surfboolean fillsurf highordertet
-		elemfacecenter barydualmesh meshinterp meshremap extrudesurf"
+				elemfacecenter barydualmesh meshinterp meshremap extrudesurf"
 func_polyline="slicesurf slicesurf3 polylinelen polylinesimplify polylineinterp closestnode"
 func_fileio="saveasc savedxf savestl savebinstl saveinr saveoff 
                 savesmf savesurfpoly savegts readgts savemsh
                 savevrml readasc readinr readmedit readoff readsmf
-	        readtetgen deletemeshfile mcpath mwpath savemedit
-		savejson loadjson saveubjson loadubjson loadmsgpack savemsgpack
+	            readtetgen deletemeshfile mcpath mwpath savemedit
+				savejson loadjson saveubjson loadubjson loadmsgpack savemsgpack
                 savebj loadbj savemphtxt savetetgenele savetetgennode saveabaqus
-		savenirfast readnirfast readnifti readmptiff"
+				savenirfast readnirfast readnifti readmptiff"
 func_jdata="savejmesh loadjnifti savejnifti loadnifti savenifti jdataencode 
                 jdatadecode jload jsave decodevarname encodevarname jnifticreate
                 nifticreate nii2jnii jnii2nii niicodemap niiformat savebnii savejnii"
 func_compression="zlibencode zlibdecode gzipencode gzipdecode lzmaencode 
                 lzmadecode lzipencode lzipdecode lz4encode lz4decode lz4hcencode 
-		lz4hcdecode base64decode base64encode"
+				lz4hcdecode base64decode base64encode"
 func_binimage="bwislands fillholes3d deislands2d deislands3d ndgaussian ndimfilter
-                imedge3d internalpoint smoothbinvol 
-		thickenbinvol thinbinvol maskdist"
+                imedge3d internalpoint smoothbinvol volgrow volshrink volclose volopen
+				thickenbinvol thinbinvol maskdist laplacefill"
 func_plotting="plotmesh plotsurf plottetra plotedges qmeshcut"
 func_misc="surfdiffuse volmap2mesh isoctavemesh getvarfrom raytrace linextriangle
-		getplanefrom3pt getexeext fallbackexeext iso2meshver
+				getplanefrom3pt getexeext fallbackexeext iso2meshver
                 raysurf getoptkey rotatevec3d rotmat2vec varargin2struct
                 jsonopt mergestruct orthdisk nestbracket2dim fast_match_bracket 
-		match_bracket memmapstream"
+				match_bracket memmapstream"
 
 print_title Streamlined mesh generation - shortcuts 
 print_group $func_shortcut
diff --git a/surfboolean.m b/surfboolean.m
index cfca2d9..d205a6b 100644
--- a/surfboolean.m
+++ b/surfboolean.m
@@ -118,7 +118,7 @@
             newnode(:, 4) = 1;
             newelem(:, 4) = 1;
         end
-        opstr = ' --decouple-inin 1 --shells 2'; % -q
+        opstr = ' --decouple-inin 1e-3 --shells 2'; % -q
         saveoff(node1(:, 1:3), elem1(:, 1:3), mwpath('pre_decouple1.off'));
         if (size(no, 2) ~= 3)
             opstr = ['-q --shells ' num2str(no)];
@@ -136,7 +136,7 @@
             newnode(:, 4) = 1;
             newelem(:, 4) = 1;
         end
-        opstr = [' --decouple-outout 1 --shells 2']; % -q
+        opstr = [' --decouple-outout 1e-3 --shells 2']; % -q
         saveoff(node1(:, 1:3), elem1(:, 1:3), mwpath('pre_decouple1.off'));
         if (size(no, 2) ~= 3)
             opstr = ['-q --shells ' num2str(no)];
diff --git a/test/run_iso2mesh_test.m b/test/run_iso2mesh_test.m
index 66d3e01..dad30f9 100644
--- a/test/run_iso2mesh_test.m
+++ b/test/run_iso2mesh_test.m
@@ -134,6 +134,9 @@ function run_iso2mesh_test(tests)
     test_iso2mesh('elemvolume', @savejson, unique(round_to_digits(elemvolume(no, el), 6)), '[0.083333]');
     test_iso2mesh('surfvolume', @savejson, surfvolume(no, fc), '[1]');
     test_iso2mesh('insurface', @savejson, insurface(no, fc, [1.5, -0.9, 2.1; 1, 0, 2; -1, 0 2; 1.2, -0, 2.5])', '[1,1,0,1]');
+
+    [no1, fc1] = highordertet(no, el);
+    test_iso2mesh('highordertet', @savejson, fc1, '[[1,7,3,10,8,14],[15,21,17,24,22,28],[2,3,7,11,13,14],[16,17,21,25,27,28],[1,5,7,9,10,22],[15,19,21,23,24,32],[4,7,5,17,15,22],[18,21,19,31,29,32],[2,7,6,13,12,25],[16,21,20,27,26,33],[4,6,7,16,17,25],[18,20,21,30,31,33]]');
 end
 
 %%
@@ -204,4 +207,72 @@ function run_iso2mesh_test(tests)
     test_iso2mesh('meshrefine elem sizefield node', @savejson, size(no1) - size(no), '[34,0]');
     test_iso2mesh('meshrefine elem sizefield elem', @savejson, size(el1) - size(el), '[105,0]');
     test_iso2mesh('meshrefine elem sizefield face', @savejson, size(fc1) - size(fc), '[64,1]');
+
+    [no, el] = meshgrid5(1:2, -1:0, 2:0.5:3);
+    [no, el] = removeisolatednode(no, volface(el));
+    no1 = sms(no, el, 10);
+    [no2, el2] = s2m(no1, el, 1, 100);
+    test_iso2mesh('sms laplacianhc', @savejson, sum(elemvolume(no2, el2(:, 1:4))) > 0.8, '[true]');
+
+    no1 = sms(no, el, 10, 0.5, 'laplacian');
+    [no2, el2] = s2m(no1, el, 1, 100);
+    test_iso2mesh('sms laplacian', @savejson, sum(elemvolume(no2, el2(:, 1:4))) < 0.1, '[true]');
+
+    no1 = sms(no, el, 10, 0.5, 'lowpass');
+    [no2, el2] = s2m(no1, el, 1, 100);
+    test_iso2mesh('sms lowpass', @savejson, sum(elemvolume(no2, el2(:, 1:4))) > 0.55, '[true]');
+end
+
+%%
+if (ismember('boolean', tests))
+    [no1, el1] = meshgrid5(1:2, 1:2, 1:2);
+    el1 = volface(el1);
+    [no1, el1] = removeisolatednode(no1, el1);
+    [no2, el2] = meshgrid6(1.7:4, 1.7:4, 1.7:4);
+    el2 = volface(el2);
+    [no2, el2] = removeisolatednode(no2, el2);
+
+    [no3, el3] = surfboolean(no1, el1, 'and', no2, el2);
+    [no3, el3] = meshcheckrepair(no3, el3, 'dup', 'tolerance', 1e-4);
+    [node, elem] = s2m(no3, el3, 1, 100);
+
+    test_iso2mesh('surfboolean and', @savejson, round_to_digits(sum(elemvolume(node, elem(:, 1:4))), 5), '[0.027]');
+
+    [no3, el3] = surfboolean(no1, el1, 'or', no2, el2);
+    [no3, el3] = meshcheckrepair(no3, el3, 'dup', 'tolerance', 1e-4);
+    [node, elem] = s2m(no3, el3, 1, 100);
+    test_iso2mesh('surfboolean or', @savejson, round(sum(elemvolume(node, elem(:, 1:4))) * 1000), '[8973]');
+
+    [no3, el3] = surfboolean(no1, el1, 'diff', no2, el2);
+    [no3, el3] = meshcheckrepair(no3, el3, 'dup', 'tolerance', 1e-4);
+    [node, elem] = s2m(no3, el3, 1, 100);
+    test_iso2mesh('surfboolean diff', @savejson, round_to_digits(sum(elemvolume(node, elem(:, 1:4))), 5), '[0.973]');
+
+    [no3, el3] = surfboolean(no1, el1, 'first', no2, el2);
+    [no3, el3] = meshcheckrepair(no3, el3, 'dup', 'tolerance', 1e-4);
+    [node, elem] = s2m(no3, el3, 1, 100, 'tetgen', [1.5, 1.5, 1.5]);
+    test_iso2mesh('surfboolean first region 1', @savejson, round_to_digits(sum(elemvolume(node, elem(elem(:, 5) == 1, 1:4))), 5), '[0.973]');
+    test_iso2mesh('surfboolean first region 0', @savejson, round_to_digits(sum(elemvolume(node, elem(elem(:, 5) == 0, 1:4))), 5), '[0.027]');
+
+    [no3, el3] = surfboolean(no1, el1, 'second', no2, el2);
+    [no3, el3] = meshcheckrepair(no3, el3, 'dup', 'tolerance', 1e-4);
+    [node, elem] = s2m(no3, el3, 1, 100, 'tetgen', [2.6, 2.6, 2.6]);
+    test_iso2mesh('surfboolean second region 1', @savejson, round_to_digits(sum(elemvolume(node, elem(elem(:, 5) == 1, 1:4))), 5), '[7.973]');
+    test_iso2mesh('surfboolean second region 0', @savejson, round_to_digits(sum(elemvolume(node, elem(elem(:, 5) == 0, 1:4))), 5), '[0.027]');
+
+    [no3, el3] = surfboolean(no1, el1, 'resolve', no2, el2);
+    [no3, el3] = meshcheckrepair(no3, el3, 'dup', 'tolerance', 1e-4);
+    [node, elem] = s2m(no3, el3, 1, 100, 'tetgen', [1.5, 1.5, 1.5; 2.6, 2.6, 2.6]);
+    test_iso2mesh('surfboolean resolve region 0', @savejson, round_to_digits(sum(elemvolume(node, elem(elem(:, 5) == 0, 1:4))), 5), '[0.027]');
+    test_iso2mesh('surfboolean resolve region 1', @savejson, round_to_digits(sum(elemvolume(node, elem(elem(:, 5) == 1, 1:4))), 5), '[0.973]');
+    test_iso2mesh('surfboolean resolve region 2', @savejson, round_to_digits(sum(elemvolume(node, elem(elem(:, 5) == 2, 1:4))), 5), '[7.973]');
+    test_iso2mesh('surfboolean self intersecting test', @savejson, surfboolean(no1, el1, 'self', no2, el2), '[1]');
+
+    [no3, el3] = meshgrid5(1:0.4:1.4, 1:0.4:1.4, 1:0.4:1.4);
+    el3 = volface(el3);
+    [no3, el3] = removeisolatednode(no3, el3);
+
+    [no4, el4] = surfboolean(no1, el1, 'separate', no3, el3);
+    [node, elem] = s2m(no4, el4, 1, 100, 'tetgen', [1.5, 1.5, 1.5]);
+    test_iso2mesh('surfboolean resolve region 2', @savejson, unique(elem(:, 5))', '[0,1]');
 end