backtrack print

wortiz · Mar 8, 2024 · 7c3f8b3 · 7c3f8b3
1 parent de2532f
commit 7c3f8b3
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/mm_sol_nonlinear.c b/src/mm_sol_nonlinear.c
@@ -2067,6 +2067,7 @@ int solve_nonlinear_problem(struct GomaLinearSolverData *ams,
       if (best_norm < Epsilon[pg->imtrx][0]) {
         skip = TRUE;
       }
+      P0PRINTF("\nNewton Line Search: %f %f\n", damp, best_norm);
 
       while (!skip && (damp > min_damp)) {
         damp *= reduction_factor;
@@ -2087,6 +2088,7 @@ int solve_nonlinear_problem(struct GomaLinearSolverData *ams,
         exchange_dof(cx, dpi, R, pg->imtrx);
 
         dbl g_check = L2_norm(R, NumUnknowns[pg->imtrx]);
+        P0PRINTF("\nNewton Line Search: %f %f\n", damp, g_check);
         if (isnan(g_check)) {
           break;
         }
@@ -2101,11 +2103,13 @@ int solve_nonlinear_problem(struct GomaLinearSolverData *ams,
         g_check = 0.5 * sqrt(g_check * g_check);
 
         if (g_check <= r_check + 0.5 * slope * damp) {
+          P0PRINTF("\nNewton Line Search: stop reached %f %f %f <= %f\n", slope, damp, g_check,
+                   r_check + 0.5 * slope * damp);
           break;
         }
       }
 
-      P0PRINTF("Newton Line Search: best damping factor: %f\n", best_damp);
+      P0PRINTF("\nNewton Line Search: best damping factor: %f\n", best_damp);
       for (i = 0; i < NumUnknowns[pg->imtrx]; i++) {
         x[i] = x_save[i] - best_damp * delta_x[i];
       }