COSMOS v7.655  COSMOSv7655
(AirShowerMC)
Functions/Subroutines
mncuve.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine mncuve (FCN, FUTIL)
 

Function/Subroutine Documentation

◆ mncuve()

subroutine mncuve ( external  FCN,
external  FUTIL 
)

Definition at line 10 of file mncuve.f.

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

Referenced by mncont(), and mnexcm().

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 Makes sure that the current point is a local
25 CC minimum and that the error matrix exists,
26 CC or at least something good enough for MINOS and MNCONT
27 CC
28 *
29 * $Id: d506cm.inc,v 1.1.1.1 1996/03/07 14:31:32 mclareni Exp $
30 *
31 * $Log: d506cm.inc,v $
32 * Revision 1.1.1.1 1996/03/07 14:31:32 mclareni
33 * Minuit
34 *
35 *
36 *
37 *
38 * d506cm.inc
39 *
40  parameter(mne=100 , mni=50)
41  parameter(mnihl=mni*(mni+1)/2)
42  CHARACTER*10 cpnam
43  COMMON
44  1/mn7nam/ cpnam(mne)
45  2/mn7ext/ u(mne) ,alim(mne) ,blim(mne)
46  3/mn7err/ erp(mni) ,ern(mni) ,werr(mni) ,globcc(mni)
47  4/mn7inx/ nvarl(mne) ,niofex(mne),nexofi(mni)
48  5/mn7int/ x(mni) ,xt(mni) ,dirin(mni)
49  6/mn7fx2/ xs(mni) ,xts(mni) ,dirins(mni)
50  7/mn7der/ grd(mni) ,g2(mni) ,gstep(mni) ,gin(mne) ,dgrd(mni)
51  8/mn7fx3/ grds(mni) ,g2s(mni) ,gsteps(mni)
52  9/mn7fx1/ ipfix(mni) ,npfix
53  a/mn7var/ vhmat(mnihl)
54  b/mn7vat/ vthmat(mnihl)
55  c/mn7sim/ p(mni,mni+1),pstar(mni),pstst(mni) ,pbar(mni),prho(mni)
56 C
57  parameter(maxdbg=10, maxstk=10, maxcwd=20, maxp=30, maxcpt=101)
58  parameter(zero=0.0, one=1.0, half=0.5)
59  COMMON
60  d/mn7npr/ maxint ,npar ,maxext ,nu
61  e/mn7iou/ isysrd ,isyswr ,isyssa ,npagwd ,npagln ,newpag
62  e/mn7io2/ istkrd(maxstk) ,nstkrd ,istkwr(maxstk) ,nstkwr
63  f/mn7tit/ cfrom ,cstatu ,ctitl ,cword ,cundef ,cvrsn ,covmes
64  g/mn7flg/ isw(7) ,idbg(0:maxdbg) ,nblock ,icomnd
65  h/mn7min/ amin ,up ,edm ,fval3 ,epsi ,apsi ,dcovar
66  i/mn7cnv/ nfcn ,nfcnmx ,nfcnlc ,nfcnfr ,itaur,istrat,nwrmes(2)
67  j/mn7arg/ word7(maxp)
68  k/mn7log/ lwarn ,lrepor ,limset ,lnolim ,lnewmn ,lphead
69  l/mn7cns/ epsmac ,epsma2 ,vlimlo ,vlimhi ,undefi ,bigedm,updflt
70  m/mn7rpt/ xpt(maxcpt) ,ypt(maxcpt)
71  n/mn7cpt/ chpt(maxcpt)
72  o/mn7xcr/ xmidcr ,ymidcr ,xdircr ,ydircr ,ke1cr ,ke2cr
73  CHARACTER ctitl*50, cword*(maxcwd), cundef*10, cfrom*8,
74  + cvrsn*6, covmes(0:3)*22, cstatu*10, chpt*1
75  LOGICAL lwarn, lrepor, limset, lnolim, lnewmn, lphead
76  EXTERNAL fcn,futil
77  IF (isw(4) .LT. 1) THEN
78  WRITE (isyswr,'(/A,A)')
79  + ' FUNCTION MUST BE MINIMIZED BEFORE CALLING ',cfrom
80  apsi = epsi
81  CALL mnmigr(fcn,futil)
82  ENDIF
83  IF (isw(2) .LT. 3) THEN
84  CALL mnhess(fcn,futil)
85  IF (isw(2) .LT. 1) THEN
86  CALL mnwarn('W',cfrom,'NO ERROR MATRIX. WILL IMPROVISE.')
87  DO 555 i=1,npar
88  ndex = i*(i-1)/2
89  DO 554 j=1,i-1
90  ndex = ndex + 1
91  554 vhmat(ndex) = 0.
92  ndex = ndex + 1
93  IF (g2(i) .LE. zero) THEN
94  wint = werr(i)
95  iext = nexofi(i)
96  IF (nvarl(iext) .GT. 1) THEN
97  CALL mndxdi(x(i),i,dxdi)
98  IF (abs(dxdi) .LT. .001) THEN
99  wint = .01
100  ELSE
101  wint = wint/abs(dxdi)
102  ENDIF
103  ENDIF
104  g2(i) = up/wint**2
105  ENDIF
106  vhmat(ndex) = 2./g2(i)
107  555 CONTINUE
108  isw(2) = 1
109  dcovar = 1.
110  ELSE
111  CALL mnwerr
112  ENDIF
113  ENDIF
114  RETURN
parameter
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
z
nodes z
Definition: cstdatmos0.fNew2.h:93
up
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/
e
dE dx *! Nuc Int sampling table e
Definition: cblkMuInt.h:130
i
nodes i
Definition: cstdatmos0.fNew2.h:93
mnwerr
subroutine mnwerr
Definition: mnwerr.f:10
mnmigr
subroutine mnmigr(FCN, FUTIL)
Definition: mnmigr.f:25
nrl_atmos::h
real(4), dimension(:), allocatable, save h
Definition: cNRLAtmos.f:28
xs
! 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
p
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
Definition: cblkHeavy.h:7
j
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
Definition: cblkHeavy.h:36
mnhess
subroutine mnhess(FCN, FUTIL)
Definition: mnhess.f:10
o
struct ob o[NpMax]
Definition: Zprivate.h:34
mndxdi
subroutine mndxdi(PINT, IPAR, DXDI)
Definition: mndxdi.f:10
d
dE dx *! Nuc Int sampling table d
Definition: cblkMuInt.h:130
m
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data m
Definition: ZlatfitBD.h:35
nrl_atmos::a
real(4), save a
Definition: cNRLAtmos.f:20
g
dE dx *! Nuc Int sampling table g
Definition: cblkMuInt.h:130
nrl_atmos::b
real(4), save b
Definition: cNRLAtmos.f:21
maxp
!onst int maxp
Definition: Zprivate.h:3
mnwarn
subroutine mnwarn(COPT, CORG, CMES)
Definition: mnwarn.f:10
n
integer n
Definition: Zcinippxc.h:1
softenpik::half
integer, parameter half
Definition: csoftenPiK.f:108
x
! 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
f
dE dx *! Nuc Int sampling table f
Definition: cblkMuInt.h:130
c
dE dx *! Nuc Int sampling table c
Definition: cblkMuInt.h:130
Here is the call graph for this function:
Here is the caller graph for this function: