Allow passing of user data for callback procedure

Author: awvwgk

src/minpack.f90

@@ -87,11 +87,12 @@ program test
       end do
       end program test
 
-      subroutine fcn(n,x,Fvec,Iflag)
+      subroutine fcn(n,x,Fvec,Iflag,udata)
       implicit none
 
       integer n , Iflag
       double precision x(n) , Fvec(n)
+      class(*), optional :: udata
 !     **********
 !
 !     THE CALLING SEQUENCE OF FCN SHOULD BE IDENTICAL TO THE

test/test_hybrd.f90

@@ -93,11 +93,12 @@ program test
       end do
       end program test
 
-      subroutine fcn(n,x,Fvec,Fjac,Ldfjac,Iflag)
+      subroutine fcn(n,x,Fvec,Fjac,Ldfjac,Iflag,udata)
       implicit none
 
       integer n , Ldfjac , Iflag
       double precision x(n) , Fvec(n) , Fjac(Ldfjac,n)
+      class(*), optional :: udata
 !     **********
 !
 !     THE CALLING SEQUENCE OF FCN SHOULD BE IDENTICAL TO THE

test/test_hybrj.f90

@@ -94,7 +94,7 @@ program test
 !  FUNCTION AND JACOBIAN SUBROUTINES SSQFCN AND SSQJAC WITH
 !  THE APPROPRIATE VALUE OF PROBLEM NUMBER (NPROB).
 
-    subroutine fcn(m, n, x, Fvec, Fjac, Ldfjac, Iflag)
+    subroutine fcn(m, n, x, Fvec, Fjac, Ldfjac, Iflag, udata)
 
         implicit none
 
@@ -112,6 +112,7 @@ subroutine fcn(m, n, x, Fvec, Fjac, Ldfjac, Iflag)
         real(wp), intent(in) :: x(n) !! independant variable vector
         real(wp), intent(inout) :: fvec(m) !! value of function at `x`
         real(wp), intent(inout) :: fjac(ldfjac, n) !! jacobian matrix at `x`
+        class(*), intent(inout), optional :: udata
 
         select case (iflag)
         case (1)
@@ -994,4 +995,4 @@ end subroutine ssqfcn
 
 !*****************************************************************************************
     end program test
-!*****************************************************************************************
+!*****************************************************************************************

test/test_lmder.f90

@@ -96,11 +96,12 @@ program test
       end do
       end program test
 
-      subroutine fcn(m,n,x,Fvec,Iflag)
+      subroutine fcn(m,n,x,Fvec,Iflag,udata)
       implicit none
 
       integer m , n , Iflag
       double precision x(n) , Fvec(m)
+      class(*), optional :: udata
 !     **********
 !
 !     THE CALLING SEQUENCE OF FCN SHOULD BE IDENTICAL TO THE

test/test_lmdif.f90

@@ -103,11 +103,12 @@ program test
    end program test
 
 
-      subroutine fcn(m,n,x,Fvec,Fjrow,Iflag)
+      subroutine fcn(m,n,x,Fvec,Fjrow,Iflag,udata)
       implicit none
 
       integer m , n , Iflag
       double precision x(n) , Fvec(m) , Fjrow(n)
+      class(*), optional :: udata
 !     **********
 !
 !     THE CALLING SEQUENCE OF FCN SHOULD BE IDENTICAL TO THE