diff --git a/examples/bouncingBall/softSphereFalling.py b/examples/bouncingBall/softSphereFalling.py
index 58c84ef..d323f6f 100644
--- a/examples/bouncingBall/softSphereFalling.py
+++ b/examples/bouncingBall/softSphereFalling.py
@@ -34,7 +34,7 @@ def createScene(rootNode):
                 rootNode.gravity=[0,-9810, 0]
                 rootNode.dt = 0.0001
                 rootNode.addObject('FreeMotionAnimationLoop')
-                rootNode.addObject('ProjectedGaussSeidelConstraintSolver', name='GSSolver', maxIterations='10000', tolerance='1e-15')
+                rootNode.addObject('BlockGaussSeidelConstraintSolver', name='GSSolver', maxIterations='10000', tolerance='1e-15')
                 rootNode.addObject('CollisionPipeline', verbose='0')
                 rootNode.addObject('BruteForceBroadPhase', name='N2')
                 rootNode.addObject('BVHNarrowPhase')
diff --git a/examples/bouncingBall/softSphereFalling_reduced.py b/examples/bouncingBall/softSphereFalling_reduced.py
index 27aa98e..5a36f30 100644
--- a/examples/bouncingBall/softSphereFalling_reduced.py
+++ b/examples/bouncingBall/softSphereFalling_reduced.py
@@ -38,7 +38,7 @@ def createScene(rootNode):
                 rootNode.gravity=[0,-9810, 0]
                 rootNode.dt = 0.0001
                 rootNode.addObject('FreeMotionAnimationLoop')
-                rootNode.addObject('ProjectedGaussSeidelConstraintSolver', name='GSSolver', maxIterations='10000', tolerance='1e-15')
+                rootNode.addObject('BlockGaussSeidelConstraintSolver', name='GSSolver', maxIterations='10000', tolerance='1e-15')
                 rootNode.addObject('CollisionPipeline', verbose="0")
                 rootNode.addObject('BruteForceBroadPhase', name="N2")
                 rootNode.addObject('BVHNarrowPhase')
diff --git a/examples/others/SphereOnAPlane/softSphereFalling.pyscn b/examples/others/SphereOnAPlane/softSphereFalling.pyscn
index 255d951..fd107f7 100644
--- a/examples/others/SphereOnAPlane/softSphereFalling.pyscn
+++ b/examples/others/SphereOnAPlane/softSphereFalling.pyscn
@@ -32,12 +32,12 @@ def createScene(rootNode):
     lambdaCoeffs = "lambdaCoeffsSphereFine6.txt"
 
     rootNode.addObject('FreeMotionAnimationLoop')
-    #rootNode.addObject('ProjectedGaussSeidelConstraintSolver', printLog=True, tolerance="1e-6", maxIterations="500", storeLambdas=True)
+    #rootNode.addObject('BlockGaussSeidelConstraintSolver', printLog=True, tolerance="1e-6", maxIterations="500", storeLambdas=True)
     if reductionOn:
-        rootNode.addObject('ProjectedGaussSeidelConstraintSolver', name='GSSolver', maxIterations='10000', tolerance='1e-15')
+        rootNode.addObject('BlockGaussSeidelConstraintSolver', name='GSSolver', maxIterations='10000', tolerance='1e-15')
     else:
-        #rootNode.addObject('ProjectedGaussSeidelConstraintSolver', name='GSSolver', maxIterations='10000', tolerance='1e-15',storeLambdas=True, lambdaPath="lambdaStored.txt",computeConstraintForces=True)
-        rootNode.addObject('ProjectedGaussSeidelConstraintSolver', name='GSSolver', maxIterations='10000', tolerance='1e-15',computeConstraintForces=True)
+        #rootNode.addObject('BlockGaussSeidelConstraintSolver', name='GSSolver', maxIterations='10000', tolerance='1e-15',storeLambdas=True, lambdaPath="lambdaStored.txt",computeConstraintForces=True)
+        rootNode.addObject('BlockGaussSeidelConstraintSolver', name='GSSolver', maxIterations='10000', tolerance='1e-15',computeConstraintForces=True)
         rootNode.addObject('PythonScriptController', filename="lambdaDumper.py", classname="storeLambda", variables="Zetest.txt")
 
 
diff --git a/examples/others/caduceus/caduceusNG.pyscn b/examples/others/caduceus/caduceusNG.pyscn
index 3f39692..236003b 100644
--- a/examples/others/caduceus/caduceusNG.pyscn
+++ b/examples/others/caduceus/caduceusNG.pyscn
@@ -34,7 +34,7 @@ def createScene(rootNode):
 
 
     rootNode.addObject('FreeMotionAnimationLoop')
-    rootNode.addObject('ProjectedGaussSeidelConstraintSolver', printLog='0', tolerance="1e-6", maxIterations="500")
+    rootNode.addObject('BlockGaussSeidelConstraintSolver', printLog='0', tolerance="1e-6", maxIterations="500")
     rootNode.addObject('CollisionPipeline', verbose="0")
     rootNode.addObject('BruteForceBroadPhase', name="N2")
     rootNode.addObject('BVHNarrowPhase')
diff --git a/examples/others/caduceus/reduced/reduced_caduceusNG.py b/examples/others/caduceus/reduced/reduced_caduceusNG.py
index 41405e4..5e6b6b4 100644
--- a/examples/others/caduceus/reduced/reduced_caduceusNG.py
+++ b/examples/others/caduceus/reduced/reduced_caduceusNG.py
@@ -150,7 +150,7 @@ def createScene(rootNode):
     rootNode.VisualStyle.displayFlags="showForceFields"
     
     rootNode.addObject('FreeMotionAnimationLoop')
-    rootNode.addObject('ProjectedGaussSeidelConstraintSolver', printLog='0', tolerance="1e-6", maxIterations="500")
+    rootNode.addObject('BlockGaussSeidelConstraintSolver', printLog='0', tolerance="1e-6", maxIterations="500")
     rootNode.addObject('CollisionPipeline', verbose="0")
     rootNode.addObject('BruteForceBroadPhase', name="N2")
     rootNode.addObject('BVHNarrowPhase')
diff --git a/examples/others/hexaBeam/hexaBeam.pyscn b/examples/others/hexaBeam/hexaBeam.pyscn
index c257033..f73d2e4 100644
--- a/examples/others/hexaBeam/hexaBeam.pyscn
+++ b/examples/others/hexaBeam/hexaBeam.pyscn
@@ -25,7 +25,7 @@ def createScene(rootNode):
     rootNode.addObject('RequiredPlugin', name='Sofa.Component.Visual') # Needed to use components [VisualStyle]  
 
     rootNode.addObject('FreeMotionAnimationLoop')
-    rootNode.addObject('ProjectedGaussSeidelConstraintSolver', tolerance="1e-12", maxIterations="10000")
+    rootNode.addObject('BlockGaussSeidelConstraintSolver', tolerance="1e-12", maxIterations="10000")
 
     
     M1 = rootNode.addChild("M1")
diff --git a/examples/softRobots/diamond/diamondRobot.py b/examples/softRobots/diamond/diamondRobot.py
index 66897e5..165f7cf 100644
--- a/examples/softRobots/diamond/diamondRobot.py
+++ b/examples/softRobots/diamond/diamondRobot.py
@@ -55,7 +55,7 @@ def createScene(rootNode):
     rootNode.addObject('RequiredPlugin', pluginName=plugins, printLog=False)
     
     rootNode.addObject('FreeMotionAnimationLoop')
-    rootNode.addObject('ProjectedGaussSeidelConstraintSolver', tolerance="1e-6", maxIterations="1000")
+    rootNode.addObject('BlockGaussSeidelConstraintSolver', tolerance="1e-6", maxIterations="1000")
     rootNode.addObject('OglSceneFrame', style="Arrows", alignment="TopRight")
     rootNode.addObject('VisualStyle', displayFlags='showVisualModels showForceFields')
 
diff --git a/examples/softRobots/diamond/diamondRobotHyperElastic.py b/examples/softRobots/diamond/diamondRobotHyperElastic.py
index 178e230..6c7e2a5 100644
--- a/examples/softRobots/diamond/diamondRobotHyperElastic.py
+++ b/examples/softRobots/diamond/diamondRobotHyperElastic.py
@@ -64,7 +64,7 @@ def createScene(rootNode):
     rootNode.addObject('OglSceneFrame', style="Arrows", alignment="TopRight")
 
     rootNode.addObject('FreeMotionAnimationLoop')
-    rootNode.addObject('ProjectedGaussSeidelConstraintSolver', tolerance="1e-6", maxIterations="1000")
+    rootNode.addObject('BlockGaussSeidelConstraintSolver', tolerance="1e-6", maxIterations="1000")
 
     modelNode = rootNode.addChild('modelNode')
     modelNode.addObject('EulerImplicitSolver', rayleighStiffness='0.1', rayleighMass='0.1')
diff --git a/examples/softRobots/finger/finger.pyscn b/examples/softRobots/finger/finger.pyscn
index 2486624..7b03c5a 100644
--- a/examples/softRobots/finger/finger.pyscn
+++ b/examples/softRobots/finger/finger.pyscn
@@ -35,7 +35,7 @@ def createScene(rootNode):
     scene.VisualStyle.displayFlags='showBehavior'
 
     rootNode.addObject("FreeMotionAnimationLoop")
-    rootNode.addObject("ProjectedGaussSeidelConstraintSolver", maxIterations=1000, tolerance=0.001)
+    rootNode.addObject("BlockGaussSeidelConstraintSolver", maxIterations=1000, tolerance=0.001)
 
     #finger = ElasticMaterialObject(name="finger",
                                    #volumeMeshFileName="mesh/finger.vtk",
diff --git a/examples/softRobots/multiGait/multiGait.py b/examples/softRobots/multiGait/multiGait.py
index 4697e04..3b394ba 100644
--- a/examples/softRobots/multiGait/multiGait.py
+++ b/examples/softRobots/multiGait/multiGait.py
@@ -41,7 +41,7 @@ def createScene(rootNode):
     rootNode.addObject('FreeMotionAnimationLoop')
     rootNode.addObject('EulerImplicitSolver', name='odesolver', firstOrder="false", rayleighStiffness='0.1',
                        rayleighMass='0.1')
-    rootNode.addObject('ProjectedGaussSeidelConstraintSolver', printLog='0', tolerance="1e-15", maxIterations="5000")
+    rootNode.addObject('BlockGaussSeidelConstraintSolver', printLog='0', tolerance="1e-15", maxIterations="5000")
     rootNode.addObject('CollisionPipeline', verbose="0")
     rootNode.addObject('BruteForceBroadPhase', name="N2")
     rootNode.addObject('BVHNarrowPhase')
diff --git a/examples/softRobots/multiGait/reducedMultiGait-softRobot.pyscn b/examples/softRobots/multiGait/reducedMultiGait-softRobot.pyscn
index ff822af..c79fb19 100644
--- a/examples/softRobots/multiGait/reducedMultiGait-softRobot.pyscn
+++ b/examples/softRobots/multiGait/reducedMultiGait-softRobot.pyscn
@@ -219,7 +219,7 @@ def createScene(rootNode):
                 rootNode.addObject('VisualStyle', displayFlags='showVisualModels hideBehaviorModels hideCollisionModels hideBoundingCollisionModels showForceFields showInteractionForceFields hideWireframe')
 
                 rootNode.addObject('FreeMotionAnimationLoop')
-                rootNode.addObject('ProjectedGaussSeidelConstraintSolver', printLog=False, tolerance="1e-4", maxIterations="1000")
+                rootNode.addObject('BlockGaussSeidelConstraintSolver', printLog=False, tolerance="1e-4", maxIterations="1000")
                 rootNode.addObject('CollisionPipeline')
                 rootNode.addObject('BruteForceBroadPhase', name="N2")
                 rootNode.addObject('BVHNarrowPhase')