same for lmdiff

Author: jacobwilliams

test/test_lmder.f90

@@ -97,9 +97,9 @@ program test_lmder
                 factor = ten*factor
 
                 ! compare with previously generated solutions:
-                if (    info_original(ic)<5 .and. &   ! ignore any where the original minpack failed
-                        any(abs( solution(ic) - x)>tol .and. &
-                        abs((solution(ic) - x)/(solution(ic))) > solution_reltol)) then
+                if (info_original(ic)<5 .and. &   ! ignore any where the original minpack failed
+                    any(abs( solution(ic) - x)>tol .and. &
+                    abs((solution(ic) - x)/(solution(ic))) > solution_reltol)) then
                     write(nwrite,'(A)') 'Failed case'
                     write(nwrite, '(//5x, a//(5x, 5d15.7))') 'Expected x: ', solution(ic)
                     write(nwrite, '(/5x, a//(5x, 5d15.7))')  'Computed x: ', x

test/test_lmdif.f90

@@ -24,6 +24,10 @@ program test_lmdif
     integer,dimension(ncases),parameter :: ms      = [10,50,10,50,10,50,2,3,4,2,15,11,16,31,31,31,10,10,20,8,8,9,10,10,30,40,33,65]
     integer,dimension(ncases),parameter :: ntriess = [1,1,1,1,1,1,3,3,3,3,3,3,2,3,3,3,1,1,3,3,1,1,1,3,1,1,1,1]
 
+    integer,dimension(*),parameter :: info_original = [1,1,1,1,1,1,2,2,2,2,2,2,5,5,5,1,1,1,1,1,1,1,1,5,1,4,&
+                                                       1,1,1,1,2,2,2,2,2,2,1,5,1,5,1,1,1,1,2,1,1,2,1,2,2,1,1]
+                                                       !! original `info` from the original minpack
+
     integer :: i, ic, info, k, m, n, NFEv, NJEv, NPRob, ntries, icase, lwa
     real(wp) :: factor, fnorm1, fnorm2
     integer :: ma(60), na(60), nf(60), nj(60), np(60), nx(60)
@@ -86,8 +90,9 @@ program test_lmdif
                 factor = ten*factor
 
                 ! compare with previously generated solutions:
-                if (any(abs( solution(ic) - x)>tol .and. &
-                        abs((solution(ic) - x)/(solution(ic))) > solution_reltol)) then
+                if (info_original(ic)<5 .and. &   ! ignore any where the original minpack failed
+                    any(abs( solution(ic) - x)>tol .and. &
+                    abs((solution(ic) - x)/(solution(ic))) > solution_reltol)) then
                     write(nwrite,'(A)') 'Failed case'
                     write(nwrite, '(//5x, a//(5x, 5d15.7))') 'Expected x: ', solution(ic)
                     write(nwrite, '(/5x, a//(5x, 5d15.7))')  'Computed x: ', x