/* end if residual raises again */
#ifdef RESIDUAL
//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) {
+ printf(" %.1f%%", (prev_residual-curr_residual)/prev_residual*100);
+ fflush(stdout);
+ if (curr_residual == 0.0) {
+ printf(" 0");
+ fflush(stdout);
+ break;
+ }
+ if ((prev_residual-curr_residual)/prev_residual < target_residual_change)
+ break;
}
- if (target_residual_change && (prev_residual-curr_residual)/prev_residual < target_residual_change)
- break;
prev_residual = curr_residual;
#endif
}