Go to the source code of this file.
|
subroutine | mneig (A, NDIMA, N, MITS, WORK, PRECIS, IFAULT) |
|
◆ mneig()
subroutine mneig |
( |
dimension(ndima,*) |
A, |
|
|
|
NDIMA, |
|
|
|
N, |
|
|
|
MITS, |
|
|
dimension(*) |
WORK, |
|
|
|
PRECIS, |
|
|
|
IFAULT |
|
) |
| |
Definition at line 10 of file mneig.f.
References a, b, c, e, f, h, i, j, m, n, o, parameter(), pt, r, and z.
Referenced by mnpsdf().
23 IMPLICIT DOUBLE PRECISION (
a-
h,
o-
z)
27 dimension
a(ndima,*),work(*)
43 IF(gl .GT. tol)
GO TO 30
52 IF(
f .GE. zero) gl = -gl
62 40 gl = gl+
a(
j,k)*
a(
i,k)
68 45 gl = gl+
a(k,
j)*
a(
i,k)
78 a(
j,k) =
a(
j,k)-
f*work(
n+k)-gl*
a(
i,k)
88 IF(work(
i) .EQ. zero .OR. l .EQ. 0)
GO TO 100
93 80 gl = gl+
a(
i,k)*
a(k,
j)
100 IF(l .EQ. 0)
GO TO 110
112 130 work(i0) = work(i0+1)
118 h = precis*(abs(work(l))+abs(work(
n+l)))
125 IF(abs(work(
n+
m)) .LE.
b)
GO TO 150
129 150
IF(
m .EQ. l)
GO TO 205
131 160
IF(
j .EQ. mits)
RETURN 134 pt = (work(l+1)-work(l))/(two*work(
n+l))
138 IF(
pt .LT. zero) pr=
pt-
r 140 h = work(l)-work(
n+l)/pr
142 170 work(
i) = work(
i)-
h 155 IF(abs(
pt) .GE. abs(work(
n+
i)))
GO TO 180
159 work(
n+
j) = s*work(
n+
i)*
r 168 190
pt =
c*work(
i)-s*gl
169 work(
j) =
h+s*(
c*gl+s*work(
i))
178 IF(abs(work(
n+l)) .GT.
b)
GO TO 160
180 205 work(l) = work(l)+
f 188 IF(work(
j) .GE.
pt)
GO TO 220
194 IF(k .EQ.
i)
GO TO 240
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
dE dx *! Nuc Int sampling table e
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz r
! 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 t endmap map ! pt before pz is set real pt
real(4), dimension(:), allocatable, save h
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data m
dE dx *! Nuc Int sampling table f
dE dx *! Nuc Int sampling table c