character*240 record
character*240 string
external fgvalue,optsave
+ common /lbfgstat/ status,niter,ncalls
c
c
c initialize some values to be used below
c
do while (.not. done)
niter = niter + 1
+c write (jout,*) "LBFGS niter",niter
muse = min(niter-1,msav)
m = m + 1
if (m .gt. msav) m = 1
c perform line search along the new conjugate direction
c
status = blank
+c write (jout,*) "Before search"
call search (nvar,f,g,x0,p,f_move,angle,ncalls,fgvalue,status)
+c write (jout,*) "After search"
c
c update variables based on results of this iteration
c