lminew.F

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file lminew.F

C **********************************************************************

      SUBROUTINE lminew

      IMPLICIT NONE

C...This is the MINUIT routine MINNEW.
C     C        THIS IS THE MAIN PROGRAM, DISGUISED AS A SUBROUTINE FOR

CC        REASONS OF COMPATIBILITY BETWEEN SYSTEMS.    IT INITIALIZES
CC        SOME CONSTANTS IN COMMON INCLUDING THE LOGICAL I/O UNIT NOS.)
CC        THEN VERIFIES THAT FCN GIVES THE SAME VALUE WHEN CALLED
CC        TWICE WITH THE SAME ARGUMENTS, AND PASSES CONTROL TO LMCMND.
CC

      COMMON /lpflag/ lst3
      INTEGER lst3
      SAVE /lpflag/

      COMMON
     1/lmmine/ erp(30)  ,ern(30)
     2/lmpari/ x(15)    ,xt(15)   ,dirin(15) ,maxint     ,npar
     3/lmpare/ u(30)              ,werr(30)  ,maxext     ,nu
     4/lmlimi/ alim(30) ,blim(30) ,lcode(30) ,lcorsp(30) ,limset
     5/lmvari/ v(15,15)
     7/lmfix / ipfix(15),xs(15)   ,xts(15)   ,dirins(15) ,npfix
     7/lmfix2/ grds(15) ,g2s(15)  ,gsteps(15),aberfs(15)
     c/lmcasc/ y(16)    ,jh       ,jl
     f/lmderi/ gin(30)  ,grd(15)  ,g2(15)    ,gstep(15)  ,aberf(15)
     g/lmsimv/ p(15,16) ,pstar(15),pstst(15) ,pbar(15)   ,prho(15)
     j/lmvart/ vt(15,15)
      COMMON
     6/lmunit/ isysrd   ,isyswr   ,isyspu
     8/lmtitl/ title(13),date(2)  ,isw(7)    ,nblock
     9/lmconv/ epsi     ,apsi     ,vtest     ,nstepq     ,nfcn ,nfcnmx
     a/lmcard/ cword    ,cword2   ,cword3    ,word7(7)
     b/lmmini/ amin     ,up       ,newmin    ,itaur      ,sigma,epsmac

      REAL erp,ern
      INTEGER maxint,npar
      REAL  x,xt,dirin
      INTEGER maxext     ,nu
      REAL u,werr
      INTEGER lcode,lcorsp,limset
      REAL alim,blim
      REAL v
      INTEGER ipfix,npfix
      REAL xs,xts,dirins
      REAL grds,g2s,gsteps,aberfs
      REAL y
      INTEGER jh,jl
      REAL gin,grd,g2,gstep,aberf
      REAL p,pstar,pstst,pbar,prho
      REAL vt
*****      COMMON
      INTEGER isysrd   ,isyswr   ,isyspu
      REAL title,date
      INTEGER isw,nblock
      REAL epsi     ,apsi     ,vtest
      INTEGER nstepq     ,nfcn ,nfcnmx
      REAL cword    ,cword2   ,cword3    ,word7
      REAL amin     ,up      ,sigma,epsmac
      INTEGER newmin    ,itaur
      REAL f
      INTEGER i

C     UNIT NUMBERS FOR CARD READER, PRINTER, PUNCH
C
      isysrd = 5
      isyswr = 6
      isyspu = 7
      maxint=15
      maxext=30
C                   DETERMINE MACHINE ACCURACY EPSMAC
      epsmac = 0.5
      DO 33 i= 1, 100
      epsmac = epsmac * 0.5
      IF ((1.0+epsmac) .EQ. 1.0)  go to 35
   33 CONTINUE
      epsmac = 1.0e-6
   35 epsmac = 2.0 * epsmac
C                             . . . . . . . . .
  110 CONTINUE
      nfcn = 1
      CALL lmidat
      CALL lminto(x)
      IF(lst3.GE.5) WRITE (isyswr,120)
  120 FORMAT (/,'0FIRST ENTRY TO FCN ')
      CALL lsigmx(npar,gin,amin,u,1)
      CALL lsigmx(npar,gin,amin,u,4)
      CALL lmprin(1,amin)
      CALL lsigmx(npar,gin,f   ,u,4)
      IF  (f .NE. amin)  go to 160
      nfcn = 3
      CALL lmcmnd
      RETURN
  160 CONTINUE
      IF(lst3.GE.1) WRITE (isyswr,880) amin, f
      IF(lst3.GE.2) stop
  880 FORMAT('0FOR THE ABOVE VALUES OF THE PARAMETERS, FCN IS TIME-',
     +'DEPENDENT',/,'0F = ',e22.14,' FOR FIRST CALL',/,' F =',e22.14,
     +' FOR SECOND')
      END