mnlims.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	mnlims

Function/Subroutine Documentation

mnlims()

subroutine mnlims

(

)

Definition at line 25 of file mnlims.f.

References a, b, c, d, e, f, from, g, h, softenpik::half, i, j, m, maxp, mndxdi(), mnexin(), mnrset(), n, o, p, parameter(), to, up(), x, xs, and z.

Referenced by mnset().

*
* $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 MNSET
CC       Interprets the SET LIM command, to reset the parameter limits
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
C
      cfrom = 'SET LIM '
      nfcnfr = nfcn
      cstatu= 'NO CHANGE '
      i2 = word7(1)
      IF (i2 .GT. maxext .OR. i2 .LT. 0)  GO TO 900
      IF (i2 .GT. 0)  GO TO 30
C                                     set limits on all parameters
      newcod = 4
      IF (word7(2) .EQ. word7(3))  newcod = 1
      DO 20 inu= 1, nu
      IF (nvarl(inu) .LE. 0)  GO TO 20
      IF (nvarl(inu).EQ.1 .AND. newcod.EQ.1)  GO TO 20
      kint = niofex(inu)
C             see if parameter has been fixed
      IF (kint .LE. 0)  THEN
         IF (isw(5) .GE. 0)  WRITE (isyswr,'(11X,A,I3)')
     +      ' LIMITS NOT CHANGED FOR FIXED PARAMETER:',inu
         GO TO 20
      ENDIF
      IF (newcod .EQ. 1)  THEN
C            remove limits from parameter
         IF (isw(5) .GT. 0)     WRITE (isyswr,134)  inu
         cstatu = 'NEW LIMITS'
         CALL mndxdi(x(kint),kint,dxdi)
         snew = gstep(kint)*dxdi
         gstep(kint) = abs(snew)
         nvarl(inu) = 1
      ELSE
C             put limits on parameter
         alim(inu) = min(word7(2),word7(3))
         blim(inu) = max(word7(2),word7(3))
         IF (isw(5) .GT. 0) WRITE (isyswr,237)  inu,alim(inu),blim(inu)
         nvarl(inu) = 4
         cstatu = 'NEW LIMITS'
         gstep(kint) = -0.1
      ENDIF
   20 CONTINUE
      GO TO 900
C                                       set limits on one parameter
   30 IF (nvarl(i2) .LE. 0)  THEN
        WRITE (isyswr,'(A,I3,A)') ' PARAMETER ',i2,' IS NOT VARIABLE.'
        GO TO 900
      ENDIF
      kint = niofex(i2)
C                                       see if parameter was fixed
      IF (kint .EQ. 0)  THEN
         WRITE (isyswr,'(A,I3)')
     +     ' REQUEST TO CHANGE LIMITS ON FIXED PARAMETER:',i2
         DO 82 ifx= 1, npfix
         IF (i2 .EQ. ipfix(ifx)) GO TO 92
   82    CONTINUE
         WRITE (isyswr,'(A)') ' MINUIT BUG IN MNLIMS. SEE F. JAMES'
   92    CONTINUE
      ENDIF
      IF (word7(2) .NE. word7(3))  GO TO 235
C                                       remove limits
      IF (nvarl(i2) .NE. 1)  THEN
         IF (isw(5) .GT. 0)  WRITE (isyswr,134)  i2
  134    FORMAT (30h limits removed from PARAMETER  ,i4)
         cstatu = 'NEW LIMITS'
         IF (kint .LE. 0)  THEN
            gsteps(ifx) = abs(gsteps(ifx))
         ELSE
            CALL mndxdi(x(kint),kint,dxdi)
            IF (abs(dxdi) .LT. 0.01)  dxdi=0.01
            gstep(kint) = abs(gstep(kint)*dxdi)
            grd(kint) = grd(kint)*dxdi
         ENDIF
         nvarl(i2) = 1
      ELSE
         WRITE (isyswr,'(A,I3)') ' NO LIMITS SPECIFIED.  PARAMETER ',
     +        i2,' IS ALREADY UNLIMITED.  NO CHANGE.'
      ENDIF
      GO TO 900
C                                        put on limits
  235 alim(i2) = min(word7(2),word7(3))
      blim(i2) = max(word7(2),word7(3))
      nvarl(i2) = 4
      IF (isw(5) .GT. 0)   WRITE (isyswr,237)  i2,alim(i2),blim(i2)
  237 FORMAT (10h PARAMETER ,i3, 14h limits set to  ,2g15.5)
      cstatu = 'NEW LIMITS'
      IF (kint .LE. 0)  THEN
         gsteps(ifx) = -0.1
      ELSE
         gstep(kint) = -0.1
      ENDIF
C
  900 CONTINUE
      IF (cstatu .NE. 'NO CHANGE ')  THEN
        CALL mnexin(x)
        CALL mnrset(1)
      ENDIF
      RETURN

Here is the call graph for this function:

Here is the caller graph for this function: