1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
|
C THE CORE SUBROUTINE CALLED BY CUTTOOLS WHICH CONTAINS THE HELAS
C CALLS BUILDING THE LOOP
SUBROUTINE ML5_0_LOOPNUM(Q,RES)
C
C CONSTANTS
C
INTEGER NLOOPGROUPS
PARAMETER (NLOOPGROUPS=77)
INTEGER LOOPMAXCOEFS
PARAMETER (LOOPMAXCOEFS=70)
C These are constants related to the split orders
INTEGER NSQUAREDSO
PARAMETER (NSQUAREDSO=0)
C
C ARGUMENTS
C
COMPLEX*16 Q(0:3)
COMPLEX*16 RES
C
C GLOBAL VARIABLES
C
INTEGER ID,SQSOINDEX,RANK
COMMON/ML5_0_LOOP/ID,SQSOINDEX,RANK
COMPLEX*16 LOOPCOEFS(0:LOOPMAXCOEFS-1,NSQUAREDSO,NLOOPGROUPS)
COMMON/ML5_0_LCOEFS/LOOPCOEFS
RES=(0.0D0,0.0D0)
CALL ML5_0_EVAL_POLY(LOOPCOEFS(0,SQSOINDEX,ID),RANK,-Q,RES)
END
SUBROUTINE ML5_0_MPLOOPNUM(Q,RES)
C
C MODULE
C
INCLUDE 'cts_mprec.h'
C
C CONSTANTS
C
INTEGER NLOOPGROUPS
PARAMETER (NLOOPGROUPS=77)
INTEGER NEXTERNAL
PARAMETER (NEXTERNAL=5)
INTEGER LOOPMAXCOEFS
PARAMETER (LOOPMAXCOEFS=70)
C These are constants related to the split orders
INTEGER NSQUAREDSO
PARAMETER (NSQUAREDSO=0)
C
C ARGUMENTS
C
INCLUDE 'cts_mpc.h'
$ , INTENT(IN), DIMENSION(0:3) :: Q
INCLUDE 'cts_mpc.h'
$ , INTENT(OUT) :: RES
C
C LOCAL VARIABLES
C
COMPLEX*32 QRES
REAL*8 DUMMY(3,0:NSQUAREDSO)
REAL*16 QPP(0:3,NEXTERNAL)
COMPLEX*32 QQ(0:3)
INTEGER I,J
C
C GLOBAL VARIABLES
C
LOGICAL MP_DONE
COMMON/ML5_0_MP_DONE/MP_DONE
INTEGER ID,SQSOINDEX,RANK
COMMON/ML5_0_LOOP/ID,SQSOINDEX,RANK
COMPLEX*32 LOOPCOEFS(0:LOOPMAXCOEFS-1,NSQUAREDSO,NLOOPGROUPS)
COMMON/ML5_0_MP_LCOEFS/LOOPCOEFS
C MP_PS IS THE FIXED (POSSIBLY IMPROVED) MP PS POINT AND MP_P IS
C THE ONE WHICH CAN BE MODIFIED (I.E. ROTATED ETC.) FOR STABILITY
C PURPOSE
REAL*16 MP_PS(0:3,NEXTERNAL),MP_P(0:3,NEXTERNAL)
COMMON/ML5_0_MP_PSPOINT/MP_PS,MP_P
C ----------
C BEGIN CODE
C ----------
DO I=0,3
QQ(I) = Q(I)
ENDDO
QRES=(0.0E0_16,0.0E0_16)
IF (.NOT.MP_DONE) THEN
CALL ML5_0_MP_COMPUTE_LOOP_COEFS(MP_P,DUMMY)
MP_DONE=.TRUE.
ENDIF
CALL MP_ML5_0_EVAL_POLY(LOOPCOEFS(0,SQSOINDEX,ID),RANK,-QQ,QRES)
RES=QRES
END
SUBROUTINE ML5_0_MPLOOPNUM_DUMMY(Q,RES)
C
C MODULE
C
INCLUDE 'cts_mprec.h'
C
C ARGUMENTS
C
INCLUDE 'cts_mpc.h'
$ , INTENT(IN), DIMENSION(0:3) :: Q
INCLUDE 'cts_mpc.h'
$ , INTENT(OUT) :: RES
C
C LOCAL VARIABLES
C
COMPLEX*16 DRES
COMPLEX*16 DQ(0:3)
INTEGER I
C ----------
C BEGIN CODE
C ----------
DO I=0,3
DQ(I) = Q(I)
ENDDO
CALL ML5_0_LOOPNUM(DQ,DRES)
RES=DRES
END
|