--- /dev/null
+puts "========"
+puts "OCC26896"
+puts "========"
+puts ""
+#################################
+# Wrong result of bsection algorithm
+#################################
+
+restore [locate_data_file bug26896_cutshape1.brep] b1
+
+vertex v1 554.852479447 5561.24970151 -18225.8420231
+vertex v2 19830.8118639 5561.31737822 -18184.3230952
+vertex v3 19967.9643885 5620.52434157 18138.50479
+vertex v4 692.019599679 5620.80022039 18307.753443
+
+edge e1 v1 v2
+edge e2 v2 v3
+edge e3 v3 v4
+edge e8 v4 v1
+
+wire w3 e1 e2 e3 e8
+mkplane f3 w3
+
+bsection result b1 f3
+
+regexp {Tolerance +MAX=([-0-9.+eE]+)} [tolerance result] full MaxTolerance
+
+set expected_MaxTolerance 0.0023476472179396139
+set tol_abs_MaxTolerance 0.0
+set tol_rel_MaxTolerance 1.0e-4
+checkreal "MaxTolerance" ${MaxTolerance} ${expected_MaxTolerance} ${tol_abs_MaxTolerance} ${tol_rel_MaxTolerance}
+
+set nbshapes_expected "
+Number of shapes in shape
+ VERTEX : 2
+ EDGE : 1
+ WIRE : 0
+ FACE : 0
+ SHELL : 0
+ SOLID : 0
+ COMPSOLID : 0
+ COMPOUND : 1
+ SHAPE : 4
+"
+checknbshapes result -ref ${nbshapes_expected} -t -m "BSECTION"
+
+smallview
+
+donly result
+fit
+display b1 f3
+
+set only_screen_axo 1