From 300ebd02b07458adadcba8c5ceb872a453971c01 Mon Sep 17 00:00:00 2001 From: Johannes Pfeifer Date: Sat, 12 Mar 2016 19:43:08 +0100 Subject: [PATCH] solve1: Move condition that could never be reached into if-clause --- matlab/solve1.m | 15 +++++++-------- 1 file changed, 7 insertions(+), 8 deletions(-) diff --git a/matlab/solve1.m b/matlab/solve1.m index 937909d09..5d5241fce 100644 --- a/matlab/solve1.m +++ b/matlab/solve1.m @@ -118,7 +118,13 @@ for its = 1:maxit if check > 0 den = max([f;0.5*nn]) ; if max(abs(g).*max([abs(x(j2)') ones(1,nn)])')/den < tolmin - return + if max(abs(x(j2)-xold(j2))./max([abs(x(j2)') ones(1,nn)])') < tolx + disp (' ') + disp (['SOLVE: Iteration ' num2str(its)]) + disp (['Convergence on dX.']) + disp (x) + return + end else disp (' ') disp (['SOLVE: Iteration ' num2str(its)]) @@ -127,13 +133,6 @@ for its = 1:maxit return end - if max(abs(x(j2)-xold(j2))./max([abs(x(j2)') ones(1,nn)])') < tolx - disp (' ') - disp (['SOLVE: Iteration ' num2str(its)]) - disp (['Convergence on dX.']) - disp (x) - return - end elseif max(abs(fvec)) < tolf return end