warn when cg does not converge

Author: Ahdhn

include/rxmesh/matrix/cg_mat_free_attr_solver.h

@@ -112,6 +112,11 @@ struct CGMatFreeAttrSolver : public IterativeSolver<T, Attribute<T, HandleT>>
 
             this->m_iter_taken++;
         }
+        RXMESH_WARN(
+            "CGMatFreeAttrSolver::solve() did not converge after {} iterations. Residual "
+            "= {}",
+            this->m_iter_taken,
+            delta_new);
         this->m_final_residual = delta_new;
     }
 

include/rxmesh/matrix/cg_solver.h

@@ -107,6 +107,14 @@ struct CGSolver : public IterativeSolver<T, DenseMatrix<T, DenseMatOrder>>
 
             this->m_iter_taken++;
         }
+
+
+        RXMESH_WARN(
+            "CGSolver::solve() did not converge after {} iterations. Residual "
+            "= {}",
+            this->m_iter_taken,
+            delta_new);
+
         this->m_final_residual = delta_new;
     }
 

include/rxmesh/matrix/pcg_mat_free_attr_solver.h

@@ -119,6 +119,12 @@ struct PCGMatFreeAttrSolver : public CGMatFreeAttrSolver<T, HandleT>
 
             this->m_iter_taken++;
         }
+        RXMESH_WARN(
+            "PCGMatFreeAttrSolver::solve() did not converge after {} "
+            "iterations. Residual "
+            "= {}",
+            this->m_iter_taken,
+            delta_new);
         this->m_final_residual = this->delta_new;
     }
 

include/rxmesh/matrix/pcg_solver.h

@@ -125,6 +125,13 @@ struct PCGSolver : public CGSolver<T, DenseMatOrder>
 
             this->m_iter_taken++;
         }
+        RXMESH_WARN(
+            "PCGSolver::solve() did not converge after {} "
+            "iterations. Residual "
+            "= {}",
+            this->m_iter_taken,
+            delta_new);
+
         this->m_final_residual = this->delta_new;
     }