lmcom.f90

Go to the documentation of this file.

c NLSL Version 1.9.0 beta 2/13/15
c----------------------------------------------------------------------
c                       =========================
c                             module LMCOM
c                       =========================
c     Work arrays for Levenberg-Marquardt least squares fitting
c     procedure and parameters controlling termination criteria.
c
c     NOTE: This module uses the nlsdim.inc module.
c
c----------------------------------------------------------------------
c
      module lmcom
      use nlsdim
      implicit none
c
      double precision, save ::
     #                 fjac(mxpt,mxjcol), fvec(mxpt), x(mxjcol),
     #                 diag(mxjcol), qtf(mxjcol), corr(mxjcol,mxjcol),
     #                 work1(mxjcol), work2(mxjcol), work3(mxjcol),
     #                 work4(mxpt), gnvec(mxjcol), gradf(mxjcol),
     #                 tcov(mxtv,mxtv)
c
      integer, save :: ipvt(mxjcol), nprint, lmflag
c
      integer, parameter :: nflmpr=4, nilmpr=4
c
      double precision, target, save  :: flmprm(nflmpr)
      double precision, pointer, save :: ftol, gtol, xtol, factor
c
      integer, target, save  :: ilmprm(nilmpr)
      integer, pointer, save :: maxev, maxitr, mode, info
c
      contains
c
      subroutine lmcom_init
      implicit none
c
      ftol   => flmprm(1)
      gtol   => flmprm(2)
      xtol   => flmprm(3)
      factor => flmprm(4)
c
      maxev  => ilmprm(1)
      maxitr => ilmprm(2)
      mode   => ilmprm(3)
      info   => ilmprm(4)
c
      end subroutine lmcom_init
c
      end module lmcom