mnmatu.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	mnmatu (KODE)

Function/Subroutine Documentation

mnmatu()

subroutine mnmatu

(

KODE

)

Definition at line 10 of file mnmatu.f.

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

Referenced by mnexcm(), mnmigr(), mnmnos(), mnread(), and 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        prints the covariance matrix v when KODE=1.
CC        always prints the global correlations, and
CC        calculates and prints the individual correlation coefficients
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
      dimension vline(mni)
      isw2 = isw(2)
      IF (isw2 .LT. 1)  THEN
          WRITE (isyswr,'(1X,A)')  covmes(isw2)
          GO TO 500
      ENDIF
      IF (npar .EQ. 0)  THEN
          WRITE (isyswr,'('' MNMATU: NPAR=0'')')
          GO TO 500
          ENDIF
C                                       . . . . .external error matrix
      IF (kode .EQ. 1)  THEN
         isw5 = isw(5)
         isw(5) = 2
         CALL mnemat(p,maxint)
           IF (isw2.LT.3)  WRITE (isyswr,'(1X,A)')  covmes(isw2)
         isw(5) = isw5
      ENDIF
C                                       . . . . . correlation coeffs. .
      IF (npar .LE. 1)   GO TO 500
      CALL mnwerr
C     NCOEF is number of coeff. that fit on one line, not to exceed 20
      ncoef = (npagwd-19)/6
      ncoef = min(ncoef,20)
      nparm = min(npar,ncoef)
      WRITE (isyswr, 150) (nexofi(id),id=1,nparm)
  150 FORMAT (/36h PARAMETER  correlation coefficients  /
     +         18h       no.  global   ,20i6)
      DO 200 i= 1, npar
         ix = nexofi(i)
         ndi = i*(i+1)/2
           DO 170 j= 1, npar
           m = max(i,j)
           n = min(i,j)
           ndex = m*(m-1)/2 + n
           ndj = j*(j+1)/2
  170      vline(j) = vhmat(ndex)/sqrt(abs(vhmat(ndi)*vhmat(ndj)))
         nparm = min(npar,ncoef)
         WRITE (isyswr,171)   ix, globcc(i), (vline(it),it=1,nparm)
  171    FORMAT (6x,i3,2x,f7.5,1x,20f6.3)
         IF (i.LE.nparm) GO TO 200
            DO 190 iso= 1, 10
            nsofar = nparm
            nparm = min(npar,nsofar+ncoef)
            WRITE (isyswr,181)  (vline(it),it=nsofar+1,nparm)
  181       FORMAT (19x,20f6.3)
            IF (i .LE. nparm) GO TO 192
  190       CONTINUE
  192    CONTINUE
  200 CONTINUE
      IF (isw2.LT.3)  WRITE (isyswr,'(1X,A)')  covmes(isw2)
  500 RETURN

Here is the call graph for this function:

Here is the caller graph for this function: