mesh_3_volume = Mesh_3.GetVolume()
volume_error_3 = abs(mesh_2_volume-mesh_3_volume)/mesh_2_volume
-assert abs(mesh_3_volume-mesh_2_volume)/mesh_2_volume < 1e-10
+assert math.isclose(mesh_3_volume,mesh_2_volume,rel_tol=1e-7)
# Check viscous layers with imprinting
for name in faces_imprinted: