stop if no descent

Cra2yDavid · Feb 2, 2018 · 7ad8757 · 7ad8757
1 parent 306ce09
commit 7ad8757
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/PoleFigureAnalysis/@MLSSolver/calcODF.m b/PoleFigureAnalysis/@MLSSolver/calcODF.m
@@ -46,7 +46,8 @@
   function do_iterate
 
     solver.initIter;
-    for iter = 1:7 %solver.iterMax
+    lasterr = inf;
+    for iter = 1:solver.iterMax
 
       % one step
       solver.doIter;
@@ -56,6 +57,8 @@
         err(i) = norm(solver.u{i}) ./ norm(solver.pf.allI{i}(:));
       end
       fprintf(format,xnum2str(iter,[],2),err);  
+      if (lasterr-err)/err < 0.05, break; end
+      lasterr = err;
     end
   end
 % -------------------------------------------------------------------