From: Frederic Pons <frederic.pons@opencascade.com>
Date: Tue, 28 Feb 2023 16:52:59 +0000 (+0100)
Subject: Check result of Compute in examples
X-Git-Tag: V9_11_0a1~10^2
X-Git-Url: http://git.salome-platform.org/gitweb/?a=commitdiff_plain;h=796122663b68ea1cee8c90033398a34d93518d39;p=modules%2Fsmesh.git

Check result of Compute in examples
---

diff --git a/doc/examples/a3DmeshOnModified2Dmesh.py b/doc/examples/a3DmeshOnModified2Dmesh.py
index fb49f873c..33b9a47ea 100644
--- a/doc/examples/a3DmeshOnModified2Dmesh.py
+++ b/doc/examples/a3DmeshOnModified2Dmesh.py
@@ -56,7 +56,9 @@ isDone = Mesh_1.SplitQuadObject( Mesh_1, Diag13=True )
 Mesh_1.Tetrahedron()
 
 # 5) Compute 3D mesh
-Mesh_1.Compute()
+isDone = Mesh_1.Compute()
+if not isDone:
+  raise Exception("Error when computing Mesh")
 
 if salome.sg.hasDesktop():
   salome.sg.updateObjBrowser()
diff --git a/doc/examples/cartesian_algo.py b/doc/examples/cartesian_algo.py
index 584d66715..690aaff8a 100644
--- a/doc/examples/cartesian_algo.py
+++ b/doc/examples/cartesian_algo.py
@@ -25,7 +25,7 @@ coords = list(range(-100,100,10))
 cartHyp = cartAlgo.SetGrid( coords,coords,coords, 1000000)
 
 # compute the mesh
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 print("nb hexahedra",mesh.NbHexas())
 print("nb tetrahedra",mesh.NbTetras())
 print("nb polyhedra",mesh.NbPolyhedrons())
@@ -34,7 +34,7 @@ print()
 # define the grid by setting constant spacing
 cartHyp = cartAlgo.SetGrid( "10","10","10", 1000000)
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 print("nb hexahedra",mesh.NbHexas())
 print("nb tetrahedra",mesh.NbTetras())
 print("nb polyhedra",mesh.NbPolyhedrons())
@@ -44,7 +44,7 @@ print()
 # activate creation of faces
 cartHyp.SetToCreateFaces( True )
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 print("nb hexahedra",mesh.NbHexas())
 print("nb tetrahedra",mesh.NbTetras())
 print("nb polyhedra",mesh.NbPolyhedrons())
@@ -53,7 +53,7 @@ print()
 
 # enable consideration of shared faces
 cartHyp.SetToConsiderInternalFaces( True )
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 print("nb hexahedra",mesh.NbHexas())
 print("nb tetrahedra",mesh.NbTetras())
 print("nb polyhedra",mesh.NbPolyhedrons())
@@ -64,7 +64,7 @@ print()
 spaceFuns = ["5","10+10*t"]
 cartAlgo.SetGrid( [spaceFuns, [0.5]], [spaceFuns, [0.5]], [spaceFuns, [0.25]], 10 )
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 print("nb hexahedra",mesh.NbHexas())
 print("nb tetrahedra",mesh.NbTetras())
 print("nb polyhedra",mesh.NbPolyhedrons())
@@ -85,24 +85,24 @@ spc = "0.1" # spacing
 mesh = smesh_builder.Mesh( box, "custom axes")
 algo = mesh.BodyFitted()
 algo.SetGrid( spc, spc, spc, 10000 )
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 print("Default axes")
 print("   nb hex:",mesh.NbHexas())
 
 # set axes using edges of the box
 algo.SetAxesDirs( xDir, [-0.1,1,0], zDir )
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 print("Manual axes")
 print("   nb hex:",mesh.NbHexas())
 
 # set optimal orthogonal axes
 algo.SetOptimalAxesDirs( isOrthogonal=True )
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 print("Optimal orthogonal axes")
 print("   nb hex:",mesh.NbHexas())
 
 # set optimal non-orthogonal axes
 algo.SetOptimalAxesDirs( isOrthogonal=False )
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 print("Optimal non-orthogonal axes")
 print("   nb hex:",mesh.NbHexas())
diff --git a/doc/examples/create_dual_mesh.py b/doc/examples/create_dual_mesh.py
index 08dfea822..d65e9bd16 100644
--- a/doc/examples/create_dual_mesh.py
+++ b/doc/examples/create_dual_mesh.py
@@ -45,15 +45,14 @@ Mesh_1 = smesh.Mesh(Sphere_1,'Mesh_1')
 status = Mesh_1.AddHypothesis( Sphere_1, NETGEN_3D_Parameters_1 )
 NETGEN_1D_2D_3D = Mesh_1.Tetrahedron(algo=smeshBuilder.NETGEN_1D2D3D)
 isDone = Mesh_1.Compute()
-
+if not isDone:
+  raise Exception("Error when computing Mesh")
 
 # Creating Dual mesh
 dual_Mesh_1 = smesh.CreateDualMesh( Mesh_1, 'dual_Mesh_1', True)
 
-
 assert(dual_Mesh_1.NbPolyhedrons() > 0)
 assert(dual_Mesh_1.NbTetras() == 0)
 
-
 if salome.sg.hasDesktop():
   salome.sg.updateObjBrowser()
diff --git a/doc/examples/creating_meshes_ex01.py b/doc/examples/creating_meshes_ex01.py
index e5cb0e934..0a174d8a3 100644
--- a/doc/examples/creating_meshes_ex01.py
+++ b/doc/examples/creating_meshes_ex01.py
@@ -27,8 +27,7 @@ algo3D.MaxElementVolume(900.)
 
 # compute the mesh
 ret = tetra.Compute()
-if ret == 0:
-    print("problem when computing the mesh")
-else:
-    print("mesh computed")
-    pass
+if not ret:
+    raise Exception("problem when computing the mesh")
+
+print("mesh computed")
diff --git a/doc/examples/creating_meshes_ex02.py b/doc/examples/creating_meshes_ex02.py
index 676708c68..bb0358837 100644
--- a/doc/examples/creating_meshes_ex02.py
+++ b/doc/examples/creating_meshes_ex02.py
@@ -52,11 +52,11 @@ mesh.Hexahedron()
 surfaces = geom_builder.SubShapeAll(box, geom_builder.ShapeType["FACE"])
 
 # method 1: no sub-mesh is created
-mesh.Compute( surfaces[0] )
+if not mesh.Compute( surfaces[0] ): raise Exception("Error when computing Mesh")
 
 # method 2: a sub-mesh is created
 submesh = mesh.GetSubMesh( surfaces[2], "submesh 2" )
-submesh.Compute()
+if not submesh.Compute(): raise Exception("Error when computing Mesh")
 
 # compute the whole mesh
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/creating_meshes_ex03.py b/doc/examples/creating_meshes_ex03.py
index 03a1804a1..88b3d0b1a 100644
--- a/doc/examples/creating_meshes_ex03.py
+++ b/doc/examples/creating_meshes_ex03.py
@@ -35,14 +35,17 @@ mesh.Triangle(geom=Face_3)
 # get existing sub-mesh priority order: F1 -> F2 -> F3
 [[SubMesh_F1, SubMesh_F3, SubMesh_F2]] = mesh.GetMeshOrder()
 isDone = mesh.Compute()
+if not isDone: raise Exception("Error when computing Mesh")
 print("Nb elements at initial order of sub-meshes:", mesh.NbElements())
 
 # set new sub-mesh order: F2 -> F1 -> F3
 isDone = mesh.SetMeshOrder([[SubMesh_F2, SubMesh_F1, SubMesh_F3]])
 isDone = mesh.Compute()
+if not isDone: raise Exception("Error when computing Mesh")
 print("Nb elements at new order of sub-meshes:", mesh.NbElements())
 
 # compute with other sub-mesh order: F3 -> F2 -> F1
 isDone = mesh.SetMeshOrder([[SubMesh_F3, SubMesh_F2, SubMesh_F1]])
 isDone = mesh.Compute()
+if not isDone: raise Exception("Error when computing Mesh")
 print("Nb elements at another order of sub-meshes:", mesh.NbElements())
diff --git a/doc/examples/creating_meshes_ex04.py b/doc/examples/creating_meshes_ex04.py
index 3298e971c..405fcd9ee 100644
--- a/doc/examples/creating_meshes_ex04.py
+++ b/doc/examples/creating_meshes_ex04.py
@@ -41,19 +41,19 @@ hyp3 = algo_local.Arithmetic1D(1, 6)
 hyp4 = algo_local.Propagation()
 
 # compute the mesh
-tria.Compute()
+if not tria.Compute(): raise Exception("Error when computing Mesh")
 PrintMeshInfo(tria)
 
 # remove a local hypothesis
 tria.RemoveHypothesis(hyp4, edge)
 
 # compute the mesh
-tria.Compute()
+if not tria.Compute(): raise Exception("Error when computing Mesh")
 PrintMeshInfo(tria)
 
 # change the value of the 2D hypothesis
 hyp2.SetMaxElementArea(2.)
 
 # compute the mesh
-tria.Compute()
+if not tria.Compute(): raise Exception("Error when computing Mesh")
 PrintMeshInfo(tria)
diff --git a/doc/examples/creating_meshes_ex05.py b/doc/examples/creating_meshes_ex05.py
index 5c2ab863b..4f4f6d171 100644
--- a/doc/examples/creating_meshes_ex05.py
+++ b/doc/examples/creating_meshes_ex05.py
@@ -24,7 +24,7 @@ tetra.Triangle()
 tetra.Tetrahedron()
 
 # compute the mesh
-tetra.Compute()
+if not tetra.Compute(): raise Exception("Error when computing Mesh")
 
 # export the mesh in a MED file
 medFile = tempfile.NamedTemporaryFile(suffix=".med").name
diff --git a/doc/examples/creating_meshes_ex06.py b/doc/examples/creating_meshes_ex06.py
index a460eeb8b..0d706851e 100644
--- a/doc/examples/creating_meshes_ex06.py
+++ b/doc/examples/creating_meshes_ex06.py
@@ -113,7 +113,7 @@ algo.Propagation()
 hexa.Quadrangle()
 hexa.Hexahedron()
 
-hexa.Compute()
+if not hexa.Compute(): raise Exception("Error when computing Mesh")
 
 hexa.Group(group_a)
 hexa.Group(group_b)
diff --git a/doc/examples/creating_meshes_ex07.py b/doc/examples/creating_meshes_ex07.py
index 4a181c237..d0d4c8e9f 100644
--- a/doc/examples/creating_meshes_ex07.py
+++ b/doc/examples/creating_meshes_ex07.py
@@ -46,7 +46,7 @@ algo1D_1=Mesh_inf.Segment()
 algo1D_1.NumberOfSegments(10)
 algo2D_1=Mesh_inf.Quadrangle()
 algo3D_1=Mesh_inf.Hexahedron()
-Mesh_inf.Compute()
+if not Mesh_inf.Compute(): raise Exception("Error when computing Mesh")
 
 # create a group on the top face
 Gsup1=Mesh_inf.Group(Fsup1, "Sup")
@@ -59,7 +59,7 @@ algo1D_2=Mesh_sup.Segment()
 algo1D_2.NumberOfSegments(5)
 algo2D_2=Mesh_sup.Quadrangle()
 algo3D_2=Mesh_sup.Hexahedron()
-Mesh_sup.Compute()
+if not Mesh_sup.Compute(): raise Exception("Error when computing Mesh")
 
 # create a group on the top face
 Gsup2=Mesh_sup.Group(Fsup2, "Sup")
diff --git a/doc/examples/creating_meshes_ex08.py b/doc/examples/creating_meshes_ex08.py
index 71fdff73b..aa6714339 100644
--- a/doc/examples/creating_meshes_ex08.py
+++ b/doc/examples/creating_meshes_ex08.py
@@ -20,7 +20,7 @@ localAlgo = mesh.Triangle(face)
 mesh.Segment().NumberOfSegments( 3 )
 mesh.Quadrangle()
 mesh.Prism()
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # objects to copy
 fGroup = mesh.GroupOnGeom( face, "2D on face")
@@ -56,4 +56,4 @@ smallBox = geom_builder.MakeScaleAlongAxes( box, None, 1, 0.5, 0.5 )
 cutBox = geom_builder.MakeCut( box, smallBox, theName="box - smallBox" )
 
 ok, newMesh, groups, submehses, hyps, invIDs = smesh_builder.CopyMeshWithGeom( mesh, cutBox, "cutBox" )
-newMesh.Compute()
+if not newMesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_adaptive1d.py b/doc/examples/defining_hypotheses_adaptive1d.py
index d28342a19..e0f567a2a 100644
--- a/doc/examples/defining_hypotheses_adaptive1d.py
+++ b/doc/examples/defining_hypotheses_adaptive1d.py
@@ -31,4 +31,4 @@ deflection = 0.05
 mesh = smesh_builder.Mesh( shape )
 mesh.Segment().Adaptive( minSize, maxSize, deflection )
 mesh.Triangle().MaxElementArea( 300 )
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex01.py b/doc/examples/defining_hypotheses_ex01.py
index 734c7ea02..48571de5f 100644
--- a/doc/examples/defining_hypotheses_ex01.py
+++ b/doc/examples/defining_hypotheses_ex01.py
@@ -40,4 +40,4 @@ hexa.Quadrangle()
 hexa.Hexahedron()
 
 # compute the mesh
-hexa.Compute()
+if not hexa.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex02.py b/doc/examples/defining_hypotheses_ex02.py
index 29e15745d..b98b1a185 100644
--- a/doc/examples/defining_hypotheses_ex02.py
+++ b/doc/examples/defining_hypotheses_ex02.py
@@ -43,4 +43,4 @@ algo_local = hexa.Segment(e_arc)
 algo_local.Deflection1D(1.0)
 
 # compute the mesh
-hexa.Compute()
+if not hexa.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex03.py b/doc/examples/defining_hypotheses_ex03.py
index 21be409db..56be06dee 100644
--- a/doc/examples/defining_hypotheses_ex03.py
+++ b/doc/examples/defining_hypotheses_ex03.py
@@ -40,4 +40,4 @@ algo_local.StartEndLength(1, 6)
 algo_local.Propagation()
 
 # compute the mesh
-hexa.Compute()
+if not hexa.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex04.py b/doc/examples/defining_hypotheses_ex04.py
index 043ff1e44..f27cb3b4c 100644
--- a/doc/examples/defining_hypotheses_ex04.py
+++ b/doc/examples/defining_hypotheses_ex04.py
@@ -40,4 +40,4 @@ algo_local.LocalLength(2.)
 algo_local.Propagation()
 
 # compute the mesh
-hexa.Compute()
+if not hexa.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex05.py b/doc/examples/defining_hypotheses_ex05.py
index 4b2d1a75b..e18ccebb8 100644
--- a/doc/examples/defining_hypotheses_ex05.py
+++ b/doc/examples/defining_hypotheses_ex05.py
@@ -40,4 +40,4 @@ algo = tria_mesh.Triangle()
 algo.MaxElementArea(100)
 
 # compute the mesh
-tria_mesh.Compute()
+if not tria_mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex06.py b/doc/examples/defining_hypotheses_ex06.py
index aa0cb45bf..351dc795c 100644
--- a/doc/examples/defining_hypotheses_ex06.py
+++ b/doc/examples/defining_hypotheses_ex06.py
@@ -31,6 +31,6 @@ algo3D.MaxElementVolume(200.)
 # compute the mesh
 ret = tetra.Compute()
 if ret == 0:
-    print("problem when computing the mesh")
-else:
-    print("Computation succeeded")
+    raise Exception("problem when computing the mesh")
+
+print("Computation succeeded")
diff --git a/doc/examples/defining_hypotheses_ex07.py b/doc/examples/defining_hypotheses_ex07.py
index 7295778f0..3cc6860bc 100644
--- a/doc/examples/defining_hypotheses_ex07.py
+++ b/doc/examples/defining_hypotheses_ex07.py
@@ -33,4 +33,4 @@ algo2D = tria.Triangle()
 algo2D.LengthFromEdges()
 
 # compute the mesh
-tria.Compute()
+if not tria.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex08.py b/doc/examples/defining_hypotheses_ex08.py
index 43d77f1c2..d59d4430f 100644
--- a/doc/examples/defining_hypotheses_ex08.py
+++ b/doc/examples/defining_hypotheses_ex08.py
@@ -39,7 +39,7 @@ algo_local.Arithmetic1D(1, 4)
 algo_local.Propagation()
 
 # compute the mesh which contains prisms
-hexa.Compute()
+if not hexa.Compute(): raise Exception("Error when computing Mesh")
 
 # create another mesh on the box
 mesh = smesh_builder.Mesh(box, "Propagation of distribution of nodes")
@@ -60,4 +60,4 @@ algo_local.Arithmetic1D(1, 4)
 algo_local.PropagationOfDistribution()
 
 # compute the mesh which contains hexahedra only
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex09.py b/doc/examples/defining_hypotheses_ex09.py
index 9646d1707..cb02c912e 100644
--- a/doc/examples/defining_hypotheses_ex09.py
+++ b/doc/examples/defining_hypotheses_ex09.py
@@ -29,7 +29,7 @@ algo3D = hexa.Hexahedron()
 algo1D.Arithmetic1D(1, 4)
 
 # compute the mesh
-hexa.Compute()
+if not hexa.Compute(): raise Exception("Error when computing Mesh")
 
 # 2. Create a tetrahedral mesh on the box
 tetra = smesh_builder.Mesh(box, "Box : tetrahedrical mesh")
@@ -48,4 +48,4 @@ algo1D.Arithmetic1D(1, 4)
 algo2D.LengthFromEdges()
 
 # compute the mesh
-tetra.Compute()
+if not tetra.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex10.py b/doc/examples/defining_hypotheses_ex10.py
index 0a1b3d85c..92e0f1163 100644
--- a/doc/examples/defining_hypotheses_ex10.py
+++ b/doc/examples/defining_hypotheses_ex10.py
@@ -47,7 +47,7 @@ src_mesh.Segment().NumberOfSegments(9,10)
 src_mesh.Quadrangle()
 src_mesh.Hexahedron()
 src_mesh.Triangle(f1) # triangular sub-mesh
-src_mesh.Compute()
+if not src_mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Mesh the box using projection algorithms
 
@@ -75,7 +75,7 @@ proj2D.SourceFace( f2 )
 # 3D hypotheses to project prisms from the source to the target mesh
 proj3D = tgt_mesh.Projection3D()
 proj3D.SourceShape3D( box, src_mesh, v1F1, v1F2, v2F1, v2F2 )
-tgt_mesh.Compute()
+if not tgt_mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Move the source mesh to visually compare the two meshes
 src_mesh.TranslateObject( src_mesh, smesh_builder.MakeDirStruct( 210, 0, 0 ), Copy=False)
diff --git a/doc/examples/defining_hypotheses_ex11.py b/doc/examples/defining_hypotheses_ex11.py
index 86796723f..0fb36d3c6 100644
--- a/doc/examples/defining_hypotheses_ex11.py
+++ b/doc/examples/defining_hypotheses_ex11.py
@@ -29,9 +29,9 @@ geom_builder.addToStudyInFather( box, Face_2, 'Face_2' )
 src_mesh = smesh_builder.Mesh(Face_1, "Source mesh")
 src_mesh.Segment().NumberOfSegments(15)
 src_mesh.Triangle()
-src_mesh.Compute()
+if not src_mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Mesh the target mesh using the algorithm Projection1D2D
 tgt_mesh = smesh_builder.Mesh(Face_2, "Target mesh")
 tgt_mesh.Projection1D2D().SourceFace(Face_1,src_mesh)
-tgt_mesh.Compute()
+if not tgt_mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex12.py b/doc/examples/defining_hypotheses_ex12.py
index d7f9ff5b6..d391b5f21 100644
--- a/doc/examples/defining_hypotheses_ex12.py
+++ b/doc/examples/defining_hypotheses_ex12.py
@@ -37,4 +37,4 @@ Mesh_1.AddHypothesis(Fixed_points_1D_1)
 Regular_1D = Mesh_1.Segment()
 Quadrangle_2D = Mesh_1.Quadrangle()
 # Compute mesh
-Mesh_1.Compute()
+if not Mesh_1.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex13.py b/doc/examples/defining_hypotheses_ex13.py
index f1b359e50..0ef0433d8 100644
--- a/doc/examples/defining_hypotheses_ex13.py
+++ b/doc/examples/defining_hypotheses_ex13.py
@@ -24,7 +24,7 @@ radial_Quad_algo = mesh.Quadrangle(algo=smeshBuilder.RADIAL_QUAD)
 # The Radial Quadrange algorithm can work without any hypothesis
 # In this case it uses "Default Nb of Segments" preferences parameter to discretize edges
 # So by default there will be 15 segments in both radial and circular directions
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # The Radial Quadrange uses global or local 1d hypotheses if it does
 # not have its own hypotheses.
@@ -32,9 +32,9 @@ mesh.Compute()
 # So that there will be 5 radial layers and 10 circular segments
 global_Nb_Segments = mesh.Segment().NumberOfSegments(5)
 local_Nb_Segments  = mesh.Segment(circle).NumberOfSegments(10)
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Define own parameters of Radial Quadrange algorithm
 # The number of radial layers will be 4
 radial_Quad_algo.NumberOfLayers( 4 )
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex14.py b/doc/examples/defining_hypotheses_ex14.py
index accc5879b..bc7ce0b99 100644
--- a/doc/examples/defining_hypotheses_ex14.py
+++ b/doc/examples/defining_hypotheses_ex14.py
@@ -28,4 +28,4 @@ Nb_Segments_1 = Regular_1D.NumberOfSegments(10)
 Quadrangle_2D = Mesh_1.Quadrangle().TriangleVertex( triaVertex )
 
 # Compute the mesh
-Mesh_1.Compute()
+if not Mesh_1.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex15.py b/doc/examples/defining_hypotheses_ex15.py
index d1c711844..06ec83c19 100644
--- a/doc/examples/defining_hypotheses_ex15.py
+++ b/doc/examples/defining_hypotheses_ex15.py
@@ -36,8 +36,8 @@ Regular_1D_1 = Mesh_1.Segment(geom=Edge_2)
 Nb_Segments_2 = Regular_1D_1.NumberOfSegments(10)
 
 # Compute mesh (with Quadrangle Preference type)
-isDone = Mesh_1.Compute()
+if not Mesh_1.Compute(): raise Exception("Error when computing Mesh")
 
 # Change type to Reduced and compute again
 Quadrangle_Parameters_1.SetQuadType( smeshBuilder.QUAD_REDUCED )
-isDone = Mesh_1.Compute()
+if not Mesh_1.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex16.py b/doc/examples/defining_hypotheses_ex16.py
index 850bc4865..c3fb9ec5b 100644
--- a/doc/examples/defining_hypotheses_ex16.py
+++ b/doc/examples/defining_hypotheses_ex16.py
@@ -33,7 +33,7 @@ geom_builder.addToStudyInFather( boxes[1], midFace1, "middle Face")
 srcMesh = smesh_builder.Mesh(boxes[0], "source mesh") # box coloser to CS origin
 nSeg1 = srcMesh.Segment().NumberOfSegments(4)
 srcMesh.Quadrangle()
-srcMesh.Compute()
+if not srcMesh.Compute(): raise Exception("Error when computing Mesh")
 srcFaceGroup = srcMesh.GroupOnGeom( midFace0, "src faces", SMESH.FACE )
 
 # Import faces from midFace0 to the target mesh
@@ -43,8 +43,8 @@ importAlgo = tgtMesh.UseExisting2DElements(midFace1)
 import2hyp = importAlgo.SourceFaces( [srcFaceGroup] )
 tgtMesh.Segment().NumberOfSegments(3)
 tgtMesh.Quadrangle()
-tgtMesh.Compute()
+if not tgtMesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Import the whole source mesh with groups
 import2hyp.SetCopySourceMesh(True,True)
-tgtMesh.Compute()
+if not tgtMesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_ex17.py b/doc/examples/defining_hypotheses_ex17.py
index 917ddd487..9861b2ac4 100644
--- a/doc/examples/defining_hypotheses_ex17.py
+++ b/doc/examples/defining_hypotheses_ex17.py
@@ -43,7 +43,7 @@ layersHyp = algo3D.ViscousLayers(thickness,numberOfLayers,stretchFactor,
                                  ignoreFaces,           # optional
                                  groupName = groupName) # optional
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # retrieve boundary prisms created by mesh.Compute()
 boundaryGroup = mesh.GetGroupByName( layersHyp.GetGroupName() )[0]
@@ -65,7 +65,7 @@ mesh.Segment().NumberOfSegments( 5 )
 vlHyp = mesh.Triangle().ViscousLayers2D( 2, 3, 1.5,
                                          edgeIds, isEdgesToIgnore=True, # optional
                                          groupName=groupName)           # optional
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # retrieve boundary elements created by mesh.Compute()
 quadrangles = mesh.GetGroupByName( vlHyp.GetGroupName() )[0]
@@ -74,4 +74,4 @@ print( "Nb boundary quadrangles", quadrangles.Size() )
 # viscous layers will be created on 3 edges, as we pass isEdgesToIgnore=False
 vlHyp.SetEdges( edgeIds, False )
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/defining_hypotheses_len_near_vertex.py b/doc/examples/defining_hypotheses_len_near_vertex.py
index f2f024ffd..5fca17b5b 100644
--- a/doc/examples/defining_hypotheses_len_near_vertex.py
+++ b/doc/examples/defining_hypotheses_len_near_vertex.py
@@ -28,4 +28,4 @@ mesh.Hexahedron()
 # define refinement near vertices
 algo1d.LengthNearVertex( 0.2 )
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/filters_ex09.py b/doc/examples/filters_ex09.py
index d78f685ba..30bf3d9e9 100644
--- a/doc/examples/filters_ex09.py
+++ b/doc/examples/filters_ex09.py
@@ -15,7 +15,7 @@ face = geom_builder.MakeFaceHW(100, 100, 1, theName="quadrangle")
 mesh = smesh_builder.Mesh(face)
 mesh.Segment().NumberOfSegments(10)
 mesh.Triangle().MaxElementArea(25)
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # get all free borders
 filter = smesh_builder.GetFilter(SMESH.EDGE, SMESH.FT_FreeBorders)
diff --git a/doc/examples/filters_ex10.py b/doc/examples/filters_ex10.py
index 8f27a7c18..741bfe188 100644
--- a/doc/examples/filters_ex10.py
+++ b/doc/examples/filters_ex10.py
@@ -16,7 +16,7 @@ geom_builder.addToStudy( face, "quadrangle" )
 mesh = smesh_builder.Mesh(face)
 mesh.Segment().NumberOfSegments(10)
 mesh.Triangle().MaxElementArea(25)
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # get all faces with free edges
 filter = smesh_builder.GetFilter(SMESH.FACE, SMESH.FT_FreeEdges)
diff --git a/doc/examples/filters_ex16.py b/doc/examples/filters_ex16.py
index 9b9f5523e..01b2fa468 100644
--- a/doc/examples/filters_ex16.py
+++ b/doc/examples/filters_ex16.py
@@ -16,7 +16,7 @@ mesh = smesh_builder.Mesh( box, "Box" )
 mesh.Segment().NumberOfSegments(10)
 mesh.Quadrangle()
 mesh.Hexahedron()
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 # copy all elements with translation and Merge nodes
 mesh.TranslateObject( mesh, smesh_builder.MakeDirStruct( 10,0,0), Copy=True )
 mesh.MergeNodes( mesh.FindCoincidentNodes(1e-7) )
diff --git a/doc/examples/filters_ex17.py b/doc/examples/filters_ex17.py
index 43a5af866..9d6c27c75 100644
--- a/doc/examples/filters_ex17.py
+++ b/doc/examples/filters_ex17.py
@@ -16,7 +16,7 @@ mesh = smesh_builder.Mesh( box, "Box" )
 mesh.Segment().NumberOfSegments(10)
 mesh.Quadrangle()
 mesh.Hexahedron()
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 # copy all elements with translation
 mesh.TranslateObject( mesh, [10,0,0], Copy=True )
 # create  a filter to find nodes equal within tolerance of 1e-5
diff --git a/doc/examples/filters_ex18.py b/doc/examples/filters_ex18.py
index f5c0e7d3b..22369a418 100644
--- a/doc/examples/filters_ex18.py
+++ b/doc/examples/filters_ex18.py
@@ -16,7 +16,7 @@ mesh = smesh_builder.Mesh( box, "Box" )
 mesh.Segment().NumberOfSegments(10)
 mesh.Quadrangle()
 mesh.Hexahedron()
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 # copy all elements with translation and merge nodes
 mesh.TranslateObject( mesh, [10,0,0], Copy=True )
 mesh.MergeNodes( mesh.FindCoincidentNodes( 1e-5 ))
diff --git a/doc/examples/filters_ex39.py b/doc/examples/filters_ex39.py
index b66a8d0b7..587138ea0 100644
--- a/doc/examples/filters_ex39.py
+++ b/doc/examples/filters_ex39.py
@@ -34,7 +34,7 @@ v000 = geom_builder.MakeVertex( 100,100,0, theName='v000' ) # can use box sub-ve
 v001 = geom_builder.GetVertexNearPoint( box1, geom_builder.MakeVertex(100,100,100), theName='v001')
 ijkAlgo.Renumber([ smeshBuilder.BlockCS( box1, v000, v001 ) ])
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Create filters with FT_ConnectedElements criterion by pointing a domain in different ways:
 
diff --git a/doc/examples/generate_flat_elements.py b/doc/examples/generate_flat_elements.py
index b632fce21..f816da018 100644
--- a/doc/examples/generate_flat_elements.py
+++ b/doc/examples/generate_flat_elements.py
@@ -44,7 +44,7 @@ Mesh_1 = smesh_builder.Mesh(Partition_1)
 Mesh_1.Segment().NumberOfSegments(15)
 Mesh_1.Triangle().LengthFromEdges()
 Mesh_1.Tetrahedron()
-Mesh_1.Compute()
+if not Mesh_1.Compute(): raise Exception("Error when computing Mesh")
 
 # relevant groups of volumes and faces
 
diff --git a/doc/examples/grouping_elements_ex02.py b/doc/examples/grouping_elements_ex02.py
index e3ddd4e22..fff4b4041 100644
--- a/doc/examples/grouping_elements_ex02.py
+++ b/doc/examples/grouping_elements_ex02.py
@@ -33,7 +33,7 @@ quadra.Quadrangle()
 algo1D.NumberOfSegments(7) 
 
 # compute the mesh
-quadra.Compute()
+if not quadra.Compute(): raise Exception("Error when computing Mesh")
 
 # create SMESH group on the face with name "SMESHGroup1"
 aSmeshGroup1 = quadra.GroupOnGeom(face, "SMESHGroup1")
diff --git a/doc/examples/grouping_elements_ex03.py b/doc/examples/grouping_elements_ex03.py
index 2de59d734..2c592aebd 100644
--- a/doc/examples/grouping_elements_ex03.py
+++ b/doc/examples/grouping_elements_ex03.py
@@ -16,7 +16,7 @@ box = geom_builder.MakeBoxDXDYDZ(10,10,10)
 mesh = smesh_builder.Mesh(box,"Quad mesh")
 hyp1D = mesh.Segment().StartEndLength( 1, 4 )
 mesh.Quadrangle()
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # create a group on filter selecting faces of medium size
 critaria = [ \
@@ -30,7 +30,7 @@ print("Group on filter contains %s elements" % filtGroup.Size())
 # group on filter is updated if the mesh is modified
 hyp1D.SetStartLength( 2.5 )
 hyp1D.SetEndLength( 2.5 )
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 print("After mesh change, group on filter contains %s elements" % filtGroup.Size())
 
 # set a new filter defining the group
diff --git a/doc/examples/measurements_ex01.py b/doc/examples/measurements_ex01.py
index a0c9991a8..9e1053db2 100644
--- a/doc/examples/measurements_ex01.py
+++ b/doc/examples/measurements_ex01.py
@@ -27,8 +27,8 @@ mesh2 = smesh_builder.Mesh(box2, 'box2')
 mesh2.Segment().NumberOfSegments(5)
 mesh2.Triangle().MaxElementArea(20)
 
-mesh1.Compute()
-mesh2.Compute()
+if not mesh1.Compute(): raise Exception("Error when computing Mesh")
+if not mesh2.Compute(): raise Exception("Error when computing Mesh")
 
 # compute min distance from grp1 to the origin (not available yet)
 smesh_builder.MinDistance(mesh1)
diff --git a/doc/examples/measurements_ex02.py b/doc/examples/measurements_ex02.py
index 33bffbaea..e4d8338c6 100644
--- a/doc/examples/measurements_ex02.py
+++ b/doc/examples/measurements_ex02.py
@@ -27,8 +27,8 @@ mesh2 = smesh_builder.Mesh(box2, 'box2')
 mesh2.Segment().NumberOfSegments(5)
 mesh2.Triangle().MaxElementArea(20)
 
-mesh1.Compute()
-mesh2.Compute()
+if not mesh1.Compute(): raise Exception("Error when computing Mesh")
+if not mesh2.Compute(): raise Exception("Error when computing Mesh")
 
 # compute bounding box for mesh1
 mesh1.BoundingBox()
diff --git a/doc/examples/measurements_ex03.py b/doc/examples/measurements_ex03.py
index 44d2f6e71..77a59ee60 100644
--- a/doc/examples/measurements_ex03.py
+++ b/doc/examples/measurements_ex03.py
@@ -19,7 +19,7 @@ submesh_1d = mesh.Segment().NumberOfSegments(5)
 submesh_2d = mesh.Quadrangle()
 submesh_3d = mesh.Hexahedron()
 submesh_2d_face = mesh.Triangle(face)
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # create a group
 
diff --git a/doc/examples/mechanic.py b/doc/examples/mechanic.py
index 902ce7184..37a372e6f 100644
--- a/doc/examples/mechanic.py
+++ b/doc/examples/mechanic.py
@@ -83,4 +83,4 @@ mesh.Quadrangle(sub_face4)
 
 # finally compute whole mesh
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/mesh_3d.py b/doc/examples/mesh_3d.py
index 69bd3f0e1..0041b612c 100644
--- a/doc/examples/mesh_3d.py
+++ b/doc/examples/mesh_3d.py
@@ -73,7 +73,7 @@ algo3d = tetra.Tetrahedron()
 algo3d.MaxElementVolume(100) 
 
 # Compute the mesh
-tetra.Compute() 
+if not tetra.Compute(): raise Exception("Error when computing Mesh")
 
 # Create a mesh group of all triangles generated on geom faces present in faces_group
 group = tetra.Group(faces_group)
diff --git a/doc/examples/modifying_meshes_ex11.py b/doc/examples/modifying_meshes_ex11.py
index 268ec1235..4555823c8 100644
--- a/doc/examples/modifying_meshes_ex11.py
+++ b/doc/examples/modifying_meshes_ex11.py
@@ -16,7 +16,7 @@ mesh = smesh_builder.Mesh(box, 'box')
 mesh.Segment().NumberOfSegments(10)
 mesh.Triangle().MaxElementArea(5)
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 print("After Compute(): %s nodes, %s faces" % ( mesh.NbNodes(), mesh.NbFaces()))
 
 # remove nodes #246 and #255
diff --git a/doc/examples/modifying_meshes_ex15.py b/doc/examples/modifying_meshes_ex15.py
index cd058178a..5b570e396 100644
--- a/doc/examples/modifying_meshes_ex15.py
+++ b/doc/examples/modifying_meshes_ex15.py
@@ -14,7 +14,7 @@ box = geom_builder.MakeBoxDXDYDZ(200, 200, 200)
 mesh = smesh_builder.Mesh( box )
 mesh.Segment().AutomaticLength(0.1)
 mesh.Quadrangle()
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # find node at (0,0,0) which is located on a geom vertex
 node000 = None
diff --git a/doc/examples/modifying_meshes_ex23.py b/doc/examples/modifying_meshes_ex23.py
index 673569516..acfb7570d 100644
--- a/doc/examples/modifying_meshes_ex23.py
+++ b/doc/examples/modifying_meshes_ex23.py
@@ -46,7 +46,7 @@ def Mesh1D(shape1d, nbSeg, name):
   algo = mesh1d_tool.Segment()
   hyp  = algo.NumberOfSegments(nbSeg)
   isDone = mesh1d_tool.Compute()
-  if not isDone: print('Mesh ', name, ': computation failed')
+  if not isDone: raise Exception('Mesh %s: computation failed' % name)
   return mesh1d_tool
 
 # Create a mesh with six nodes, seven edges and two quadrangle faces
diff --git a/doc/examples/modifying_meshes_ex25.py b/doc/examples/modifying_meshes_ex25.py
index 39250d3a3..b01e6be15 100644
--- a/doc/examples/modifying_meshes_ex25.py
+++ b/doc/examples/modifying_meshes_ex25.py
@@ -28,7 +28,7 @@ algo1D.NumberOfSegments(3)
 Mesh_1.Quadrangle()
 
 isDone = Mesh_1.Compute()
-if not isDone: print('Mesh Mesh_1 : computation failed')
+if not isDone: raise Exception('Mesh Mesh_1 : computation failed')
 
 # build a triangle mesh on Face_2
 Mesh_2 = smesh_builder.Mesh(Face_2)
@@ -39,7 +39,7 @@ algo2D = Mesh_2.Triangle()
 algo2D.MaxElementArea(240)
 
 isDone = Mesh_2.Compute()
-if not isDone: print('Mesh Mesh_2 : computation failed')
+if not isDone: raise Exception('Mesh Mesh_2 : computation failed')
 
 # create a 2d pattern
 pattern = smesh_builder.GetPattern()
@@ -60,7 +60,7 @@ Mesh_3.Segment().NumberOfSegments(1)
 Mesh_3.Quadrangle()
 Mesh_3.Hexahedron()
 isDone = Mesh_3.Compute()
-if not isDone: print('Mesh Mesh_3 : computation failed')
+if not isDone: raise Exception('Mesh Mesh_3 : computation failed')
 
 # create a 3d pattern (hexahedrons)
 pattern_hexa = smesh_builder.GetPattern()
@@ -103,7 +103,7 @@ Mesh_4.Segment().NumberOfSegments(1)
 Mesh_4.Quadrangle()
 Mesh_4.Hexahedron()
 isDone = Mesh_4.Compute()
-if not isDone: print('Mesh Mesh_4 : computation failed')
+if not isDone: raise Exception('Mesh Mesh_4 : computation failed')
 
 # create another 3d pattern (pyramids)
 pattern_pyra = smesh_builder.GetPattern()
diff --git a/doc/examples/modifying_meshes_ex26.py b/doc/examples/modifying_meshes_ex26.py
index b837d6a9c..1a8ea2a23 100644
--- a/doc/examples/modifying_meshes_ex26.py
+++ b/doc/examples/modifying_meshes_ex26.py
@@ -23,7 +23,7 @@ Mesh.Tetrahedron()
 
 # compute mesh
 
-Mesh.Compute()
+if not Mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # convert to quadratic
 # theForce3d = 1; this results in the medium node lying at the
diff --git a/doc/examples/notebook_smesh.py b/doc/examples/notebook_smesh.py
index 8d90ad382..9e65604ee 100644
--- a/doc/examples/notebook_smesh.py
+++ b/doc/examples/notebook_smesh.py
@@ -37,7 +37,7 @@ algo3D = tetra.Tetrahedron()
 algo3D.MaxElementVolume("MaxElementVolume")
 
 # compute the mesh
-ret = tetra.Compute()
+if not tetra.Compute(): raise Exception("Error when computing Mesh")
 
 # translate the mesh
 point = SMESH.PointStruct("Offset", 0., 0.)
diff --git a/doc/examples/prism_3d_algo.py b/doc/examples/prism_3d_algo.py
index 20ab781e8..481fdad01 100644
--- a/doc/examples/prism_3d_algo.py
+++ b/doc/examples/prism_3d_algo.py
@@ -72,4 +72,4 @@ mesh.Segment(bigQuad).LocalLength( 3 )
 mesh.Triangle(bigQuad)
 
 # compute the mesh
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/quad_medial_axis_algo.py b/doc/examples/quad_medial_axis_algo.py
index 450827bd4..669cfd8ad 100644
--- a/doc/examples/quad_medial_axis_algo.py
+++ b/doc/examples/quad_medial_axis_algo.py
@@ -29,4 +29,4 @@ algo1d.NumberOfSegments( circNbSeg ) # division of circle edges
 algo2d = mesh.Quadrangle( smeshBuilder.QUAD_MA_PROJ )
 algo2d.StartEndLength( circleLen2 / circNbSeg, circleLen1 / circNbSeg ) # radial division
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/quality_controls_defl.py b/doc/examples/quality_controls_defl.py
index 42c04a315..e9c4074b0 100644
--- a/doc/examples/quality_controls_defl.py
+++ b/doc/examples/quality_controls_defl.py
@@ -20,7 +20,7 @@ mesh = smesh_builder.Mesh( Fuse, "Deflection_2D")
 algo = mesh.Segment()
 algo.LocalLength(35)
 algo = mesh.Triangle()
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # get min and max deflection
 minMax = mesh.GetMinMax( SMESH.FT_Deflection2D )
diff --git a/doc/examples/quality_controls_ex01.py b/doc/examples/quality_controls_ex01.py
index dda1aca12..76b8be68c 100644
--- a/doc/examples/quality_controls_ex01.py
+++ b/doc/examples/quality_controls_ex01.py
@@ -23,7 +23,7 @@ algo = mesh.Segment()
 algo.NumberOfSegments(5)
 algo = mesh.Triangle()
 algo.MaxElementArea(20.)
-mesh.Compute() 
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # criterion : free borders
 aFilter = smesh_builder.GetFilter(SMESH.EDGE, SMESH.FT_FreeBorders)
diff --git a/doc/examples/quality_controls_ex02.py b/doc/examples/quality_controls_ex02.py
index 3b48bf305..e4c6c5213 100644
--- a/doc/examples/quality_controls_ex02.py
+++ b/doc/examples/quality_controls_ex02.py
@@ -23,7 +23,7 @@ algo = mesh.Segment()
 algo.NumberOfSegments(5)
 algo = mesh.Triangle()
 algo.MaxElementArea(20.)
-mesh.Compute() 
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Criterion : Borders at multi-connection
 nb_conn = 2
diff --git a/doc/examples/quality_controls_ex03.py b/doc/examples/quality_controls_ex03.py
index 9958d7115..cfe71c97b 100644
--- a/doc/examples/quality_controls_ex03.py
+++ b/doc/examples/quality_controls_ex03.py
@@ -23,7 +23,7 @@ algo = mesh.Segment()
 algo.NumberOfSegments(5)
 algo = mesh.Triangle()
 algo.MaxElementArea(20.)
-mesh.Compute() 
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Criterion : Length > 3.
 length_margin = 3.
diff --git a/doc/examples/quality_controls_ex05.py b/doc/examples/quality_controls_ex05.py
index 4559cb897..008ee2fbc 100644
--- a/doc/examples/quality_controls_ex05.py
+++ b/doc/examples/quality_controls_ex05.py
@@ -18,7 +18,7 @@ idbox = geom_builder.addToStudy(box, "box")
 mesh = smesh_builder.Mesh(box, "Mesh_free_nodes")
 mesh.Segment().NumberOfSegments(10)
 mesh.Triangle().MaxElementArea(150.)
-mesh.Compute() 
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Remove some elements to obtain free nodes
 # Criterion : AREA < 80.
diff --git a/doc/examples/quality_controls_ex06.py b/doc/examples/quality_controls_ex06.py
index 7ecbc5ba3..22d6a9b10 100644
--- a/doc/examples/quality_controls_ex06.py
+++ b/doc/examples/quality_controls_ex06.py
@@ -42,7 +42,7 @@ Mesh_1 = smesh_builder.Mesh(Partition_1)
 Mesh_1.Segment().MaxSize(34.641)
 Mesh_1.Triangle()
 Mesh_1.Tetrahedron()
-Mesh_1.Compute()
+if not Mesh_1.Compute(): raise Exception("Error when computing Mesh")
 
 # create a group of free faces
 aFilter = smesh_builder.GetFilter(SMESH.FACE, SMESH.FT_FreeFaces )
diff --git a/doc/examples/quality_controls_ex07.py b/doc/examples/quality_controls_ex07.py
index ce88f3773..5b1b695fb 100644
--- a/doc/examples/quality_controls_ex07.py
+++ b/doc/examples/quality_controls_ex07.py
@@ -16,7 +16,7 @@ geom_builder.addToStudy( box, "box" )
 mesh = smesh_builder.Mesh(box)
 mesh.Segment().NumberOfSegments(3)
 mesh.Quadrangle()
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # remove 2 faces
 allFaces = mesh.GetElementsByType(SMESH.FACE)
diff --git a/doc/examples/quality_controls_ex11.py b/doc/examples/quality_controls_ex11.py
index 1b08179f1..731b3b01f 100644
--- a/doc/examples/quality_controls_ex11.py
+++ b/doc/examples/quality_controls_ex11.py
@@ -23,7 +23,7 @@ algo = mesh.Segment()
 algo.NumberOfSegments(5)
 algo = mesh.Triangle()
 algo.MaxElementArea(20.)
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Criterion : Length 2D > 5.7
 length_margin = 5.7
diff --git a/doc/examples/quality_controls_ex12.py b/doc/examples/quality_controls_ex12.py
index a76fa3543..1e875a8ca 100644
--- a/doc/examples/quality_controls_ex12.py
+++ b/doc/examples/quality_controls_ex12.py
@@ -23,7 +23,7 @@ algo = mesh.Segment()
 algo.NumberOfSegments(5)
 algo = mesh.Triangle()
 algo.MaxElementArea(20.)
-mesh.Compute() 
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Criterion : MULTI-CONNECTION 2D = 3
 nb_conn = 3
diff --git a/doc/examples/radial_prism_3d_algo.py b/doc/examples/radial_prism_3d_algo.py
index a2623d846..a84ecc829 100644
--- a/doc/examples/radial_prism_3d_algo.py
+++ b/doc/examples/radial_prism_3d_algo.py
@@ -33,4 +33,4 @@ mesh.Projection1D2D( faces[1] ).SourceFace( faces[0] ) # projection faces[0] ->
 prism_algo.NumberOfSegments( 4, 5. )
 
 # compute the mesh
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/split_biquad.py b/doc/examples/split_biquad.py
index 90d544e36..a8e51ed24 100644
--- a/doc/examples/split_biquad.py
+++ b/doc/examples/split_biquad.py
@@ -24,7 +24,7 @@ mesh = smesh_builder.Mesh( shape )
 mesh.Segment().NumberOfSegments(1)
 mesh.Quadrangle()
 mesh.Triangle( ff[1] )
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # make group of quadrangles and extrude them into a hexahedron
 quadGroup = mesh.Group( ff[0], "quads")
diff --git a/doc/examples/test_homard_adapt.py b/doc/examples/test_homard_adapt.py
index 0e250eee0..9f3672b63 100644
--- a/doc/examples/test_homard_adapt.py
+++ b/doc/examples/test_homard_adapt.py
@@ -40,7 +40,7 @@ smeshhomard.SetLogInFile(True)
 smeshhomard.SetLogFile(log_file_1)
 smeshhomard.SetRemoveLogOnSuccess(False)
 smeshhomard.SetVerboseLevel(3)
-smeshhomard.Compute()
+if not smeshhomard.Compute(): raise Exception("Error when computing Mesh")
 
 if osp.isfile(output_med_1):
   os.remove(output_med_1)
@@ -79,7 +79,7 @@ smeshhomard.SetLogInFile(True)
 smeshhomard.SetLogFile(log_file_2)
 smeshhomard.SetRemoveLogOnSuccess(True)
 smeshhomard.SetVerboseLevel(0)
-smeshhomard.Compute()
+if not smeshhomard.Compute(): raise Exception("Error when computing Mesh")
 
 if osp.isfile(output_med_2):
   print("Test Uniform refinement Case 2: Error: output med file has not been removed")
diff --git a/doc/examples/test_uniform_refinement.py b/doc/examples/test_uniform_refinement.py
index 96139e01b..b0e12e370 100644
--- a/doc/examples/test_uniform_refinement.py
+++ b/doc/examples/test_uniform_refinement.py
@@ -42,7 +42,7 @@ smeshhomard.SetLogInFile(True)
 smeshhomard.SetLogFile(log_file_1)
 smeshhomard.SetRemoveLogOnSuccess(False)
 smeshhomard.SetVerboseLevel(3)
-smeshhomard.Compute()
+if not smeshhomard.Compute(): raise Exception("Error when computing Mesh")
 
 if osp.isfile(output_med_1):
   os.remove(output_med_1)
@@ -85,7 +85,7 @@ smeshhomard.SetLogInFile(True)
 smeshhomard.SetLogFile(log_file_2)
 smeshhomard.SetRemoveLogOnSuccess(True)
 smeshhomard.SetVerboseLevel(0)
-smeshhomard.Compute()
+if not smeshhomard.Compute(): raise Exception("Error when computing Mesh")
 
 if osp.isfile(output_med_2):
   print("Test Uniform refinement Case 2: Error: output med file has not been removed")
diff --git a/doc/examples/transforming_meshes_ex03.py b/doc/examples/transforming_meshes_ex03.py
index 07c8bb86a..750040b2b 100644
--- a/doc/examples/transforming_meshes_ex03.py
+++ b/doc/examples/transforming_meshes_ex03.py
@@ -18,7 +18,7 @@ Regular_1D = Mesh1.Segment()
 Nb_Segments_1 = Regular_1D.NumberOfSegments(3)
 Nb_Segments_1.SetDistrType( 0 )
 Quadrangle_2D = Mesh1.Quadrangle()
-isDone = Mesh1.Compute()
+if not Mesh1.Compute(): raise Exception("Error when computing Mesh")
 
 #Perform scale operation for the whole mesh and creation of a new mesh
 newMesh = Mesh1.ScaleMakeMesh(Mesh1,SMESH.PointStruct(100,100,200),[0.5,0.3,0.7],True,"ScaledMesh")
diff --git a/doc/examples/transforming_meshes_ex06.py b/doc/examples/transforming_meshes_ex06.py
index 2600cb734..72fd27d40 100644
--- a/doc/examples/transforming_meshes_ex06.py
+++ b/doc/examples/transforming_meshes_ex06.py
@@ -40,7 +40,7 @@ algo1D.NumberOfSegments(6)
 algo2D = trias.Triangle()
 algo2D.LengthFromEdges()
 
-trias.Compute()
+if not trias.Compute(): raise Exception("Error when computing Mesh")
 
 # create a group of all triangles currently present in the mesh
 faceTriGroup = trias.Group( face1, "face triangles" )
@@ -49,7 +49,7 @@ faceTriGroup = trias.Group( face1, "face triangles" )
 circlemesh = smesh_builder.Mesh(circle, "Path mesh")
 algo = circlemesh.Segment()
 algo.NumberOfSegments(10)
-circlemesh.Compute()
+if not circlemesh.Compute(): raise Exception("Error when computing Mesh")
 
 # extrusion of the mesh
 trias.ExtrusionAlongPath([], circlemesh, circle, 1, MakeGroups=True )
diff --git a/doc/examples/transforming_meshes_ex07.py b/doc/examples/transforming_meshes_ex07.py
index db7a6cb7c..e54dec1a3 100644
--- a/doc/examples/transforming_meshes_ex07.py
+++ b/doc/examples/transforming_meshes_ex07.py
@@ -35,7 +35,7 @@ algo_local = mesh.Segment(edge1)
 algo_local.Arithmetic1D(1, 4)
 algo_local.Propagation()
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # sew border to side
 # FirstNodeIDOnFreeBorder, SecondNodeIDOnFreeBorder, LastNodeIDOnFreeBorder,
diff --git a/doc/examples/transforming_meshes_ex08.py b/doc/examples/transforming_meshes_ex08.py
index 3de67b9df..50249a7e1 100644
--- a/doc/examples/transforming_meshes_ex08.py
+++ b/doc/examples/transforming_meshes_ex08.py
@@ -35,7 +35,7 @@ algo_local = mesh.Segment(edge1)
 algo_local.Arithmetic1D(1, 4)
 algo_local.Propagation()
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # sew conform free borders
 # FirstNodeID1, SecondNodeID1, LastNodeID1, FirstNodeID2, SecondNodeID2
diff --git a/doc/examples/transforming_meshes_ex09.py b/doc/examples/transforming_meshes_ex09.py
index 19ab5ed95..7941b29f9 100644
--- a/doc/examples/transforming_meshes_ex09.py
+++ b/doc/examples/transforming_meshes_ex09.py
@@ -24,7 +24,7 @@ mesh = smesh_builder.Mesh(shape, "Two faces : quadrangle mesh")
 mesh.Segment().Arithmetic1D( 0.1, 0.4 )
 mesh.Segment(q1).NumberOfSegments( 5 )
 mesh.Quadrangle()
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # sew free borders
 
diff --git a/doc/examples/transforming_meshes_ex10.py b/doc/examples/transforming_meshes_ex10.py
index 36c2f5c22..317a61bb1 100644
--- a/doc/examples/transforming_meshes_ex10.py
+++ b/doc/examples/transforming_meshes_ex10.py
@@ -30,7 +30,7 @@ algo_local = mesh.Segment(EdgesList[8])
 algo_local.NumberOfSegments(4)
 algo_local.Propagation()
 
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # sew side elements
 
diff --git a/doc/examples/transforming_meshes_ex11.py b/doc/examples/transforming_meshes_ex11.py
index 5fddd3015..a08d21f52 100644
--- a/doc/examples/transforming_meshes_ex11.py
+++ b/doc/examples/transforming_meshes_ex11.py
@@ -21,7 +21,7 @@ mesh.Quadrangle()
 mesh.Hexahedron()
 
 # Compute mesh
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Duplicate nodes only
 
diff --git a/doc/examples/transforming_meshes_ex13.py b/doc/examples/transforming_meshes_ex13.py
index d8fade895..8ee0260b7 100644
--- a/doc/examples/transforming_meshes_ex13.py
+++ b/doc/examples/transforming_meshes_ex13.py
@@ -25,7 +25,7 @@ mesh = smesh_builder.Mesh( shape, "test_Reorient2D")
 mesh.AutomaticHexahedralization(0.5)
 localAlgo = mesh.Segment(faces[0])
 localAlgo.NumberOfSegments( 11 )
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 group = mesh.Group( faces[1] )
 
 vec = geom_builder.MakeVectorDXDYDZ( 1, 1, 1 )
diff --git a/doc/examples/use_existing_faces.py b/doc/examples/use_existing_faces.py
index 4cdb9f26a..ef13199b6 100644
--- a/doc/examples/use_existing_faces.py
+++ b/doc/examples/use_existing_faces.py
@@ -101,7 +101,7 @@ geom_builder.addToStudy( f2, "f2" )
 # compute 1D mesh
 mesh = smesh_builder.Mesh( box )
 mesh.Segment().NumberOfSegments( 5 )
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # compute 2D mesh
 mesh.Quadrangle()
@@ -109,8 +109,8 @@ mesh.UseExistingFaces(f1) # UseExistingFaces() allows using my2DMeshing();
 mesh.UseExistingFaces(f2) # assign UseExistingFaces() BEFORE calling my2DMeshing()!
 my2DMeshing(f1, mesh)
 my2DMeshing(f2, mesh)
-assert mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # compute 3D mesh
 mesh.Prism()
-assert mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
diff --git a/doc/examples/viewing_meshes_ex01.py b/doc/examples/viewing_meshes_ex01.py
index e75f036a3..ded69ba88 100644
--- a/doc/examples/viewing_meshes_ex01.py
+++ b/doc/examples/viewing_meshes_ex01.py
@@ -32,11 +32,11 @@ Regular_1D_1_1 = tetra.Segment(geom=Face_1)
 Nb_Segments_1 = Regular_1D_1_1.NumberOfSegments(5)
 Nb_Segments_1.SetDistrType( 0 )
 Quadrangle_2D = tetra.Quadrangle(geom=Face_1)
-isDone = tetra.Compute()
+if not tetra.Compute(): raise Exception("Error when computing Mesh")
 submesh = Regular_1D_1_1.GetSubMesh()
 
 # compute the mesh
-tetra.Compute()
+if not tetra.Compute(): raise Exception("Error when computing Mesh")
 
 # Creation of group
 group = tetra.CreateEmptyGroup( SMESH.FACE, 'Group' )
diff --git a/doc/examples/viewing_meshes_ex02.py b/doc/examples/viewing_meshes_ex02.py
index a2f2ed17a..ddaba327b 100644
--- a/doc/examples/viewing_meshes_ex02.py
+++ b/doc/examples/viewing_meshes_ex02.py
@@ -16,7 +16,7 @@ box = geom_builder.MakeBoxDXDYDZ(100,100,100)
 # Create a mesh
 mesh = smesh_builder.Mesh(box,"Mesh")
 mesh.AutomaticHexahedralization()
-mesh.Compute()
+if not mesh.Compute(): raise Exception("Error when computing Mesh")
 
 # Create a point
 x,y,z = 0, 0, 1