From: Andreas Eversberg <jolly@eversberg.eu>
Date: Wed, 21 Oct 2015 04:27:20 +0000 (+0200)
Subject: Stop if residual is 0
X-Git-Url: http://git.eversberg.eu/gitweb.cgi?p=colorize.git;a=commitdiff_plain;h=071dfc3639d11a5710e5f0e089aa34f07667c2b2

Stop if residual is 0
---

diff --git a/lib/multigrid.c b/lib/multigrid.c
index f3097ee..a2409c8 100644
--- a/lib/multigrid.c
+++ b/lib/multigrid.c
@@ -1434,6 +1434,10 @@ int solve_mg(int nlevel, int iters, int ncycle, darray_t **values, darray_t **ma
 //printf("residual %.8f (change=%%%.4f)\n", curr_residual, (prev_residual-curr_residual)/prev_residual*100);
 		if (target_residual_change)
 			printf(" %.1f%%", (prev_residual-curr_residual)/prev_residual*100); fflush(stdout);
+		if (curr_residual == 0.0) {
+			printf(" 0"); fflush(stdout);
+			break;
+		}
 		if (target_residual_change && (prev_residual-curr_residual)/prev_residual < target_residual_change)
 			break;
 		prev_residual = curr_residual;