mnread.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	mnread (FCN, IFLGIN, IFLGUT, FUTIL)

Function/Subroutine Documentation

mnread()

subroutine mnread	(	external	FCN,
			IFLGIN,
			IFLGUT,
		external	FUTIL
	)

Definition at line 10 of file mnread.f.

References a, b, c, d, e, f, false, g, h, softenpik::half, i, j, m, maxp, mncomd(), mnmatu(), mnpars(), mnprin(), mnseti(), mnstin(), n, o, p, parameter(), true, up(), x, xs, and z.

Referenced by minuit(), and mnintr().

*
* $Id: d506dp.inc,v 1.1.1.1 1996/03/07 14:31:32 mclareni Exp $
*
* $Log: d506dp.inc,v $
* Revision 1.1.1.1  1996/03/07 14:31:32  mclareni
* Minuit
*
*
*
*
* d506dp.inc
*
C ************ DOUBLE PRECISION VERSION *************
      IMPLICIT DOUBLE PRECISION (a-h,o-z)
CC        Called from MINUIT.  Reads all user input to MINUIT.
CC     This routine is highly unstructured and defies normal logic.
CC
CC     IFLGIN indicates the function originally requested:
CC           = 1: read one-line title
CC             2: read parameter definitions
CC             3: read MINUIT commands
CC
CC     IFLGUT= 1: reading terminated normally
CC             2: end-of-data on input
CC             3: unrecoverable read error
CC             4: unable to process parameter requests
CC             5: more than 100 incomprehensible commands
CC internally,
CC     IFLGDO indicates the subfunction to be performed on the next
CC         input record: 1: read a one-line title
CC                       2: read a parameter definition
CC                       3: read a command
CC                       4: read in covariance matrix
CC     for example, when IFLGIN=3, but IFLGDO=1, then it should read
CC       a title, but this was requested by a command, not by MINUIT.
CC
*
* $Id: d506cm.inc,v 1.1.1.1 1996/03/07 14:31:32 mclareni Exp $
*
* $Log: d506cm.inc,v $
* Revision 1.1.1.1  1996/03/07 14:31:32  mclareni
* Minuit
*
*
*
*
* d506cm.inc
*
      parameter(mne=100 , mni=50)
      parameter(mnihl=mni*(mni+1)/2)
      CHARACTER*10 cpnam
      COMMON
     1/mn7nam/ cpnam(mne)
     2/mn7ext/ u(mne)     ,alim(mne)  ,blim(mne)
     3/mn7err/ erp(mni)   ,ern(mni)   ,werr(mni)  ,globcc(mni)
     4/mn7inx/ nvarl(mne) ,niofex(mne),nexofi(mni)
     5/mn7int/ x(mni)     ,xt(mni)    ,dirin(mni)
     6/mn7fx2/ xs(mni)    ,xts(mni)   ,dirins(mni)
     7/mn7der/ grd(mni)   ,g2(mni)    ,gstep(mni) ,gin(mne) ,dgrd(mni)
     8/mn7fx3/ grds(mni)  ,g2s(mni)   ,gsteps(mni)
     9/mn7fx1/ ipfix(mni) ,npfix
     a/mn7var/ vhmat(mnihl)
     b/mn7vat/ vthmat(mnihl)
     c/mn7sim/ p(mni,mni+1),pstar(mni),pstst(mni) ,pbar(mni),prho(mni)
C
      parameter(maxdbg=10, maxstk=10, maxcwd=20, maxp=30, maxcpt=101)
      parameter(zero=0.0,  one=1.0,   half=0.5)
      COMMON
     d/mn7npr/ maxint ,npar   ,maxext ,nu
     e/mn7iou/ isysrd ,isyswr ,isyssa ,npagwd ,npagln ,newpag
     e/mn7io2/ istkrd(maxstk) ,nstkrd ,istkwr(maxstk) ,nstkwr
     f/mn7tit/ cfrom  ,cstatu ,ctitl  ,cword  ,cundef ,cvrsn ,covmes
     g/mn7flg/ isw(7) ,idbg(0:maxdbg) ,nblock ,icomnd
     h/mn7min/ amin   ,up     ,edm    ,fval3  ,epsi   ,apsi  ,dcovar
     i/mn7cnv/ nfcn   ,nfcnmx ,nfcnlc ,nfcnfr ,itaur,istrat,nwrmes(2)
     j/mn7arg/ word7(maxp)
     k/mn7log/ lwarn  ,lrepor ,limset ,lnolim ,lnewmn ,lphead
     l/mn7cns/ epsmac ,epsma2 ,vlimlo ,vlimhi ,undefi ,bigedm,updflt
     m/mn7rpt/ xpt(maxcpt)    ,ypt(maxcpt)
     n/mn7cpt/ chpt(maxcpt)
     o/mn7xcr/ xmidcr ,ymidcr ,xdircr ,ydircr ,ke1cr  ,ke2cr
      CHARACTER ctitl*50, cword*(maxcwd), cundef*10, cfrom*8,
     +          cvrsn*6,  covmes(0:3)*22, cstatu*10, chpt*1
      LOGICAL   lwarn, lrepor, limset, lnolim, lnewmn, lphead
      EXTERNAL fcn,futil
      CHARACTER  crdbuf*80, cupbuf*10
      CHARACTER cpromt(3)*40, clower*26, cupper*26
      LOGICAL leof
      DATA cpromt/' ENTER MINUIT TITLE, or "SET INPUT n" : ',
     +            ' ENTER MINUIT PARAMETER DEFINITION:     ',
     +            ' ENTER MINUIT COMMAND:                  '/
C
      DATA clower/'abcdefghijklmnopqrstuvwxyz'/
      DATA cupper/'ABCDEFGHIJKLMNOPQRSTUVWXYZ'/
C
      iflgut = 1
      iflgdo = iflgin
      leof = .false.
      incomp = 0
C                                           . . . . read next record
   10 CONTINUE
      IF (isw(6) .EQ. 1) THEN
           WRITE (isyswr,'(A)') cpromt(iflgdo)
           IF (iflgdo .EQ. 2)  lphead = .false.
      ENDIF
      crdbuf = '   '
      READ (isysrd,'(A)',err=500,end=45)  crdbuf
C
C                 CUPBUF is the first few characters in upper case
      cupbuf(1:10) = crdbuf(1:10)
      DO 12 i= 1, 10
      IF (crdbuf(i:i) .EQ. '''') GO TO 13
         DO 11 ic= 1, 26
         IF (crdbuf(i:i) .EQ. clower(ic:ic)) cupbuf(i:i)=cupper(ic:ic)
   11    CONTINUE
   12 CONTINUE
   13 CONTINUE
C                                           . .   preemptive commands
      leof = .false.
      IF (index(cupbuf,'*EOF') .EQ. 1)    THEN
         WRITE (isyswr,'(A,I3)') ' *EOF ENCOUNTERED ON UNIT NO.',isysrd
         lphead = .true.
         GO TO 50
         ENDIF
      IF (index(cupbuf,'SET INP') .EQ. 1)    THEN
         icomnd = icomnd + 1
         WRITE (isyswr, 21) icomnd,crdbuf(1:50)
   21    FORMAT (' **********'/' **',i5,' **',a/' **********')
         lphead = .true.
         GO TO 50
         ENDIF
      GO TO 80
C                                    . . hardware EOF on current ISYSRD
   45 crdbuf = '*EOF '
      WRITE (isyswr,'(A,I3)') ' END OF DATA ON UNIT NO.',isysrd
C                                     or SET INPUT command
   50 CONTINUE
         CALL mnstin(crdbuf,ierr)
         IF (ierr .EQ. 0)  GO TO 10
         IF (ierr .EQ. 2)  THEN
            IF (.NOT. leof) THEN
               WRITE (isyswr,'(A,A/)') ' TWO CONSECUTIVE EOFs ON ',
     +              'PRIMARY INPUT FILE WILL TERMINATE EXECUTION.'
               leof = .true.
               GO TO 10
            ENDIF
         ENDIF
         iflgut = ierr
         GO TO 900
   80 IF (iflgdo .GT. 1) GO TO 100
C                            read title        . . . . .   IFLGDO = 1
C              if title is 'SET TITLE', skip and read again
      IF (index(cupbuf,'SET TIT') .EQ. 1)  GO TO 10
      CALL mnseti(crdbuf(1:50))
      WRITE (isyswr,'(1X,A50)')  ctitl
      WRITE (isyswr,'(1X,78(1H*))')
         lphead = .true.
      IF (iflgin .EQ. iflgdo)  GO TO 900
      iflgdo = iflgin
      GO TO 10
C                            data record is not a title.
  100 CONTINUE
      IF (iflgdo .GT. 2)  GO TO 300
C                          expect parameter definitions.   IFLGDO = 2
C              if parameter def is 'PARAMETER', skip and read again
      IF (index(cupbuf,'PAR') .EQ. 1)  GO TO 10
C              if line starts with SET TITLE, read a title first
      IF (index(cupbuf,'SET TIT') .EQ. 1)  THEN
         iflgdo = 1
         GO TO 10
         ENDIF
C                      we really have parameter definitions now
      CALL mnpars(crdbuf,icondp)
      IF (icondp .EQ. 0)  GO TO 10
C          format error
      IF (icondp .EQ. 1)  THEN
         IF (isw(6) .EQ. 1)  THEN
           WRITE (isyswr,'(A)') ' FORMAT ERROR.  IGNORED.  ENTER AGAIN.'
           GO TO 10
         ELSE
           WRITE (isyswr,'(A)') ' ERROR IN PARAMETER DEFINITION'
           iflgut = 4
           GO TO 900
         ENDIF
      ENDIF
C                     ICONDP = 2            . . . end parameter requests
      IF (isw(5).GE.0 .AND. isw(6).LT.1) WRITE (isyswr,'(4X,75(1H*))')
      lphead = .true.
      IF (iflgin .EQ. iflgdo)  GO TO 900
      iflgdo = iflgin
      GO TO 10
C                                              . . . . .   IFLGDO = 3
C                                           read commands
  300 CONTINUE
      CALL mncomd(fcn,crdbuf,icondn,futil)
CC     ICONDN = 0: command executed normally
CC              1: command is blank, ignored
CC              2: command line unreadable, ignored
CC              3: unknown command, ignored
CC              4: abnormal termination (e.g., MIGRAD not converged)
CC              5: command is a request to read PARAMETER definitions
CC              6: 'SET INPUT' command
CC              7: 'SET TITLE' command
CC              8: 'SET COVAR' command
CC              9: reserved
CC             10: END command
CC             11: EXIT or STOP command
CC             12: RETURN command
      IF (icondn .EQ. 2 .OR. icondn .EQ. 3) THEN
         incomp = incomp + 1
         IF (incomp .GT. 100) THEN
            iflgut = 5
            GO TO 900
            ENDIF
         ENDIF
C                         parameter
      IF (icondn .EQ. 5)  iflgdo = 2
C                         SET INPUT
      IF (icondn .EQ. 6)  GO TO 50
C                         SET TITLE
      IF (icondn .EQ. 7)  iflgdo = 1
C                                        . . . . . . . . . . set covar
      IF (icondn .EQ. 8) THEN
         icomnd = icomnd + 1
         WRITE (isyswr,405) icomnd,crdbuf(1:50)
  405    FORMAT (1h ,10(1h*)/' **',i5,' **',a)
         WRITE (isyswr, '(1H ,10(1H*))' )
         npar2 = npar*(npar+1)/2
         READ (isysrd,420,err=500,end=45)  (vhmat(i),i=1,npar2)
  420    FORMAT (bn,7e11.4,3x)
         isw(2) = 3
         dcovar = 0.0
         IF (isw(5) .GE. 0)  CALL mnmatu(1)
         IF (isw(5) .GE. 1)  CALL mnprin(2,amin)
         GO TO 10
         ENDIF
      IF (icondn .LT. 10) GO TO 10
      GO TO 900
C                                              . . . . error conditions
  500 iflgut = 3
  900 RETURN

Here is the call graph for this function:

Here is the caller graph for this function: