* * $Id: lminew.F,v 1.1.1.1 1996/03/08 17:40:16 mclareni Exp $ * * $Log: lminew.F,v $ * Revision 1.1.1.1 1996/03/08 17:40:16 mclareni * Lepto63 * * C ********************************************************************** SUBROUTINE LMINEW C...This is the MINUIT routine MINNEW. CC 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 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 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