← Back to branch summary

~maddevelopers/mg5amcnlo/2.9.4

« back to all changes in this revision

Viewing changes to vendor/CutTools/src/qcdloop/npointes.f

pass to v2.0.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
vendor/CutTools/src/qcdloop/npointes.f
 
1
        program testff
 
2
        include 'ff.h'
 
3
*       call ffinit
 
4
*       lwrite = .TRUE.
 
5
        call form
 
6
        end
 
7
      subroutine FORM
 
8
      implicit DOUBLE PRECISION (A-H,O-Z)
 
9
      DOUBLE COMPLEX TADMU
 
10
      DOUBLE COMPLEX AA0,B0,B0PM,B1,B1PM,B2,C0,C1,C2,C3,D0,D1,D2,D3,D4
 
11
      COMMON PX(6),RM(4),DEL
 
12
      COMMON B0,B0PM,B1,B1PM,B2(2),C0,C1(2),C2(4),C3(6)
 
13
      COMMON D0,D1(3),D2(7),D3(13),D4(22)
 
14
      DEL=0.d0
 
15
      XM=105.6d0**2
 
16
      TADMU=AA0(XM,DEL)
 
17
      PRINT 1,TADMU
 
18
    1 FORMAT(1H ,4E20.13)
 
19
      PX(1)=-90000.d0**2
 
20
      RM(1)=105.6d0**2
 
21
      RM(2)=0.d0
 
22
      CALL NPOIN(2)
 
23
      PRINT 1,B0,B0PM,B1PM,B2(1)
 
24
      PX(1)=-105.6d0**2
 
25
      PX(2)=PX(1)
 
26
      PX(5)=-90000.d0**2
 
27
      RM(1)=80000.d0**2
 
28
      RM(2)=0.d0
 
29
      RM(3)=80000.d0**2
 
30
      CALL NPOIN(3)
 
31
      PRINT 1,C0,C1(1)
 
32
      PX(1)=-0.511d0**2
 
33
      PX(2)=PX(1)
 
34
      PX(3)=-105.6d0**2
 
35
      PX(4)=PX(3)
 
36
      PX(5)=-80000.d0**2
 
37
      PX(6)=20000.d0**2
 
38
      RM(1)=90000.d0**2
 
39
      RM(2)=0.511d0**2
 
40
      RM(3)=90000.d0**2
 
41
      RM(4)=105.6d0**2
 
42
      CALL NPOIN(4)
 
43
      PRINT 1,D0,D1(1)
 
44
      P12=-0.511d0**2
 
45
      P22=-105.6d0**2
 
46
      P1P2=0.5d0*(20000.d0**2-P12-P22)
 
47
      P10=40000.d0
 
48
      P20=-40000.d0
 
49
      DELE=2000.d0
 
50
      PM2=0.003d0**2
 
51
      BREM=ALIJ(P22,P12,P1P2,P20,P10,DELE,PM2)
 
52
      PRINT 1,BREM
 
53
      call ffexi
 
54
      STOP
 
55
C
 
56
C  THE FOLLOWING OUTPUT IS EXPECTED    *********************************
 
57
C 0.                   .9156199386460E+06
 
58
C-.3100623399361E+02 -.2054369006935E+03 -.5269961187649E-14 -.1218492318111E-08
 
59
C .7255214122160E-20  .6092365170853E-09 -.1033542556010E+02 -.6793071440282E+02
 
60
C-.2272160404013E-21  .1739461728623E-08  .3037109802926E-22 -.1286491875423E-08
 
61
C-.2462904245037E-31  .4952559280190E-18  .1165014267077E-31 -.4571732002934E-18
 
62
C-.8203023744960E-05
 
63
C
 
64
C  NUMBERS OF ORDER  E-20  ARE TO BE CONSIDERED ZERO.
 
65
C
 
66
C ******************************************************************************
 
67
C
 
68
      END