COSMOS v7.655  COSMOSv7655
(AirShowerMC)
mnstat.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine mnstat (FMIN, FEDM, ERRDEF, NPARI, NPARX, ISTAT)
 

Function/Subroutine Documentation

◆ mnstat()

subroutine mnstat (   FMIN,
  FEDM,
  ERRDEF,
  NPARI,
  NPARX,
  ISTAT 
)

Definition at line 10 of file mnstat.f.

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

10 *
11 * $Id: d506dp.inc,v 1.1.1.1 1996/03/07 14:31:32 mclareni Exp $
12 *
13 * $Log: d506dp.inc,v $
14 * Revision 1.1.1.1 1996/03/07 14:31:32 mclareni
15 * Minuit
16 *
17 *
18 *
19 *
20 * d506dp.inc
21 *
22 C ************ DOUBLE PRECISION VERSION *************
23  IMPLICIT DOUBLE PRECISION (a-h,o-z)
24 CC User-called
25 CC Provides the user with information concerning the current status
26 CC of the current minimization. Namely, it returns:
27 CC FMIN: the best function value found so far
28 CC FEDM: the estimated vertical distance remaining to minimum
29 CC ERRDEF: the value of UP defining parameter uncertainties
30 CC NPARI: the number of currently variable parameters
31 CC NPARX: the highest (external) parameter number defined by user
32 CC ISTAT: a status integer indicating how good is the covariance
33 CC matrix: 0= not calculated at all
34 CC 1= approximation only, not accurate
35 CC 2= full matrix, but forced positive-definite
36 CC 3= full accurate covariance matrix
37 CC
38 *
39 * $Id: d506cm.inc,v 1.1.1.1 1996/03/07 14:31:32 mclareni Exp $
40 *
41 * $Log: d506cm.inc,v $
42 * Revision 1.1.1.1 1996/03/07 14:31:32 mclareni
43 * Minuit
44 *
45 *
46 *
47 *
48 * d506cm.inc
49 *
50  parameter(mne=100 , mni=50)
51  parameter(mnihl=mni*(mni+1)/2)
52  CHARACTER*10 cpnam
53  COMMON
54  1/mn7nam/ cpnam(mne)
55  2/mn7ext/ u(mne) ,alim(mne) ,blim(mne)
56  3/mn7err/ erp(mni) ,ern(mni) ,werr(mni) ,globcc(mni)
57  4/mn7inx/ nvarl(mne) ,niofex(mne),nexofi(mni)
58  5/mn7int/ x(mni) ,xt(mni) ,dirin(mni)
59  6/mn7fx2/ xs(mni) ,xts(mni) ,dirins(mni)
60  7/mn7der/ grd(mni) ,g2(mni) ,gstep(mni) ,gin(mne) ,dgrd(mni)
61  8/mn7fx3/ grds(mni) ,g2s(mni) ,gsteps(mni)
62  9/mn7fx1/ ipfix(mni) ,npfix
63  a/mn7var/ vhmat(mnihl)
64  b/mn7vat/ vthmat(mnihl)
65  c/mn7sim/ p(mni,mni+1),pstar(mni),pstst(mni) ,pbar(mni),prho(mni)
66 C
67  parameter(maxdbg=10, maxstk=10, maxcwd=20, maxp=30, maxcpt=101)
68  parameter(zero=0.0, one=1.0, half=0.5)
69  COMMON
70  d/mn7npr/ maxint ,npar ,maxext ,nu
71  e/mn7iou/ isysrd ,isyswr ,isyssa ,npagwd ,npagln ,newpag
72  e/mn7io2/ istkrd(maxstk) ,nstkrd ,istkwr(maxstk) ,nstkwr
73  f/mn7tit/ cfrom ,cstatu ,ctitl ,cword ,cundef ,cvrsn ,covmes
74  g/mn7flg/ isw(7) ,idbg(0:maxdbg) ,nblock ,icomnd
75  h/mn7min/ amin ,up ,edm ,fval3 ,epsi ,apsi ,dcovar
76  i/mn7cnv/ nfcn ,nfcnmx ,nfcnlc ,nfcnfr ,itaur,istrat,nwrmes(2)
77  j/mn7arg/ word7(maxp)
78  k/mn7log/ lwarn ,lrepor ,limset ,lnolim ,lnewmn ,lphead
79  l/mn7cns/ epsmac ,epsma2 ,vlimlo ,vlimhi ,undefi ,bigedm,updflt
80  m/mn7rpt/ xpt(maxcpt) ,ypt(maxcpt)
81  n/mn7cpt/ chpt(maxcpt)
82  o/mn7xcr/ xmidcr ,ymidcr ,xdircr ,ydircr ,ke1cr ,ke2cr
83  CHARACTER ctitl*50, cword*(maxcwd), cundef*10, cfrom*8,
84  + cvrsn*6, covmes(0:3)*22, cstatu*10, chpt*1
85  LOGICAL lwarn, lrepor, limset, lnolim, lnewmn, lphead
86  fmin = amin
87  fedm = edm
88  errdef = up
89  npari = npar
90  nparx = nu
91  istat = isw(2)
92  IF (edm .EQ. bigedm) THEN
93  fedm = up
94  ENDIF
95  IF (amin .EQ. undefi) THEN
96  fmin = 0.0
97  fedm = up
98  istat= 0
99  ENDIF
100  RETURN
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
nodes z
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data h g *is for param c g data up(2, 1)/7.0d0/
dE dx *! Nuc Int sampling table e
Definition: cblkMuInt.h:130
nodes i
real(4), dimension(:), allocatable, save h
Definition: cNRLAtmos.f:28
! constants thru Cosmos real ! if multiplied to deg radian Torad ! light velocity m sec ! infinty ! kg m2 *Togpcm2 g cm2 ! g cm2 *Tokgpm2 kg m2 ! cm *Tom m ! m *Tocm cm ! g cm3 *Tokgpm3 kg m3 ! kg m3 *Togpcm3 g cm3 ! sec *Tonsec nsec ! Tesla m ! Avogadro *A2deninv ! mfp *n * xs
Definition: Zglobalc.h:18
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
Definition: cblkHeavy.h:7
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
Definition: cblkHeavy.h:36
struct ob o[NpMax]
Definition: Zprivate.h:34
dE dx *! Nuc Int sampling table d
Definition: cblkMuInt.h:130
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data m
Definition: ZlatfitBD.h:35
real(4), save a
Definition: cNRLAtmos.f:20
dE dx *! Nuc Int sampling table g
Definition: cblkMuInt.h:130
real(4), save b
Definition: cNRLAtmos.f:21
!onst int maxp
Definition: Zprivate.h:3
integer n
Definition: Zcinippxc.h:1
integer, parameter half
Definition: csoftenPiK.f:108
! structure defining a particle at production ! Basic idea of what is to be contained in ! the particle structue is that dynamical ones should be included those derivable from the particle code ! is not included ******************************************************type fmom momentum sequence union map real e endmap map real * x
Definition: Zptcl.h:21
dE dx *! Nuc Int sampling table f
Definition: cblkMuInt.h:130
dE dx *! Nuc Int sampling table c
Definition: cblkMuInt.h:130
Here is the call graph for this function: