diff --git a/SofaKernel/modules/SofaBaseLinearSolver/CGLinearSolver.inl b/SofaKernel/modules/SofaBaseLinearSolver/CGLinearSolver.inl
index 2963ada0295..4bd9913c1e7 100644
--- a/SofaKernel/modules/SofaBaseLinearSolver/CGLinearSolver.inl
+++ b/SofaKernel/modules/SofaBaseLinearSolver/CGLinearSolver.inl
@@ -216,8 +216,8 @@ void CGLinearSolver<TMatrix,TVector>::solve(Matrix& M, Vector& x, Vector& b)
         graph_den.push_back(den);
 
 
-        /// Break condition = THRESHOLD criterion regarding the denominator is reached
-        if( fabs(den)<f_smallDenominatorThreshold.getValue() )
+        /// Break condition = THRESHOLD criterion regarding the denominator is reached (but do at least one iteration)
+        if( fabs(den)<f_smallDenominatorThreshold.getValue() && nb_iter>1)
         {
             endcond = "threshold";
             if( verbose )