← Back to branch summary

~maddevelopers/mg5amcnlo/2.0.2_onshellsubtract_newresh

« back to all changes in this revision

Viewing changes to vendor/IREGI/src/oneloop/example/example.f

  • Committer: Marco Zaro
  • Date: 2015-04-09 16:12:12 UTC
  • mfrom: (250.22.93 2.3)
  • Revision ID: marco.zaro@gmail.com-20150409161212-t7umsvk7p4z983rt
merged with 2.3 r345

Show diffs side-by-side

added added

removed removed

Lines of Context:
vendor/IREGI/src/oneloop/example/example.f
 
1
      call test
 
2
      end
 
3
 
 
4
 
 
5
      subroutine test
 
6
!  ********************************************************************
 
7
!  ********************************************************************
 
8
      use avh_olo
 
9
      implicit none
 
10
      integer :: nn
 
11
      character(4) :: next
 
12
      real(kind(1d0)) :: mu
 
13
!
 
14
      call ltini !|LTVSNge26=yesyes
 
15
!#      call ffini !|LTVSNge26=yesno
 
16
!      call olo_onshell( 1d-6 )
 
17
!      call olo_unit( 6 ,'printall' )
 
18
!
 
19
      next = 'next'
 
20
      do while (next.eq.'next')
 
21
        read(5,*) nn ,mu
 
22
        call olo_scale( mu )
 
23
        call setmudim( mu*mu ) !|LT=yes
 
24
        if     (nn.eq.2) then
 
25
          call test_2
 
26
        elseif (nn.eq.3) then
 
27
          call test_3
 
28
        else!if(nn.eq.4) then
 
29
          call test_4
 
30
        endif
 
31
        read(5,*) next
 
32
      enddo
 
33
      end
 
34
 
 
35
 
 
36
      subroutine test_4
 
37
!  ********************************************************************
 
38
!  ********************************************************************
 
39
      use avh_olo
 
40
      implicit none
 
41
      complex(kind(1d0)) :: p1,p2,p3,p4,p12,p23,m1,m2,m3,m4
 
42
      complex(kind(1d0)) :: zolo(0:2),zLT,D0C
 
43
!
 
44
      read(5,*) p1
 
45
      read(5,*) p2
 
46
      read(5,*) p3
 
47
      read(5,*) p4
 
48
      read(5,*) p12
 
49
      read(5,*) p23
 
50
      read(5,*) m1
 
51
      read(5,*) m2
 
52
      read(5,*) m3
 
53
      read(5,*) m4
 
54
      write(6,101) 'p1 : (',real(p1),',',aimag(p1),')'
 
55
      write(6,101) 'p2 : (',real(p2),',',aimag(p2),')'
 
56
      write(6,101) 'p3 : (',real(p3),',',aimag(p3),')'
 
57
      write(6,101) 'p4 : (',real(p4),',',aimag(p4),')'
 
58
      write(6,101) 'p12: (',real(p12),',',aimag(p12),')'
 
59
      write(6,101) 'p23: (',real(p23),',',aimag(p23),')'
 
60
      write(6,101) 'm1 : (',real(m1),',',aimag(m1),')'
 
61
      write(6,101) 'm2 : (',real(m2),',',aimag(m2),')'
 
62
      write(6,101) 'm3 : (',real(m3),',',aimag(m3),')'
 
63
      write(6,101) 'm4 : (',real(m4),',',aimag(m4),')'
 
64
  101 format(a6,d31.24,a1,d31.24,a1)
 
65
!
 
66
      call olo( zolo ,p1,p2,p3,p4,p12,p23 ,m1,m2,m3,m4 )
 
67
!
 
68
      zLT = 0
 
69
      if (zolo(1).eq.0)                               !|LT=yes
 
70
     &  zLT = D0C( p1,p2,p3,p4,p12,p23 ,m1,m2,m3,m4 ) !|LT=yes
 
71
!
 
72
      write(6,'(a10)') ' eps^( 0):'
 
73
      if (zLT.ne.0)
 
74
     &write(6,'(a8,2d24.16)') '     LT:',zLT
 
75
      write(6,'(a8,2d24.16)') '    olo:',zolo(0)
 
76
      write(6,'(a10)') ' eps^(-1):'
 
77
      write(6,'(a8,2d24.16)') '    olo:',zolo(1)
 
78
      write(6,'(a10)') ' eps^(-2):'
 
79
      write(6,'(a8,2d24.16)') '    olo:',zolo(2)
 
80
      write(6,*)
 
81
!
 
82
      end
 
83
   
 
84
      subroutine test_3
 
85
!  ********************************************************************
 
86
!  ********************************************************************
 
87
      use avh_olo
 
88
      implicit none
 
89
      complex(kind(1d0)) :: p1,p2,p3,m1,m2,m3 ,zolo(0:2),zLT,C0C
 
90
!
 
91
      read(5,*) p1
 
92
      read(5,*) p2
 
93
      read(5,*) p3
 
94
      read(5,*) m1
 
95
      read(5,*) m2
 
96
      read(5,*) m3
 
97
      write(6,101) 'p1 : (',real(p1),',',aimag(p1),')'
 
98
      write(6,101) 'p2 : (',real(p2),',',aimag(p2),')'
 
99
      write(6,101) 'p3 : (',real(p3),',',aimag(p3),')'
 
100
      write(6,101) 'm1 : (',real(m1),',',aimag(m1),')'
 
101
      write(6,101) 'm2 : (',real(m2),',',aimag(m2),')'
 
102
      write(6,101) 'm3 : (',real(m3),',',aimag(m3),')'
 
103
  101 format(a6,d31.24,a1,d31.24,a1)
 
104
!
 
105
      call olo( zolo ,p1,p2,p3 ,m1,m2,m3 )
 
106
!
 
107
      zLT = 0
 
108
      if (zolo(1).eq.0) zLT = C0C( p1,p2,p3 ,m1,m2,m3 ) !|LT=yes
 
109
!
 
110
      write(6,'(a10)') ' eps^( 0):'
 
111
      if (zLT.ne.0)
 
112
     &write(6,'(a8,2d24.16)') '     LT:',zLT
 
113
      write(6,'(a8,2d24.16)') '    olo:',zolo(0)
 
114
      write(6,'(a10)') ' eps^(-1):'
 
115
      write(6,'(a8,2d24.16)') '    olo:',zolo(1)
 
116
      write(6,'(a10)') ' eps^(-2):'
 
117
      write(6,'(a8,2d24.16)') '    olo:',zolo(2)
 
118
      write(6,*)
 
119
!
 
120
      end
 
121
 
 
122
      subroutine test_2
 
123
!  ********************************************************************
 
124
!  ********************************************************************
 
125
      use avh_olo
 
126
      implicit none
 
127
      complex(kind(1d0)) :: p1,m1,m2 ,b0(0:2),b1(0:2),b00(0:2),b11(0:2)
 
128
      complex(kind(1d0)) :: rslt(0:2)
 
129
      complex(kind(1d0)) :: b0LT,b1LT,b00LT,b11LT,b0C,b1C,b00C,b11C
 
130
!
 
131
      read(5,*) p1
 
132
      read(5,*) m1
 
133
      read(5,*) m2
 
134
      write(6,101) 'p1 : (',real(p1),',',aimag(p1),')'
 
135
      write(6,101) 'm1 : (',real(m1),',',aimag(m1),')'
 
136
      write(6,101) 'm2 : (',real(m2),',',aimag(m2),')'
 
137
  101 format(a6,d31.24,a1,d31.24,a1)
 
138
!
 
139
      call olo( rslt ,p1 ,m1,m2 )
 
140
      call olo( b11,b00,b1,b0 ,p1 ,m1,m2 )
 
141
      b0LT  = B0C(  p1 ,m1,m2 ) !|LT=yes
 
142
      b1LT  = B1C(  p1 ,m1,m2 ) !|LT=yes
 
143
      b00LT = B00C( p1 ,m1,m2 ) !|LT=yes
 
144
      b11LT = B11C( p1 ,m1,m2 ) !|LT=yes
 
145
!
 
146
      write(6,'(a10)') ' eps^( 0):'
 
147
      write(6,'(a8,2d24.16)') 'b0   LT:',b0LT !|LT=yes
 
148
      write(6,'(a8,2d24.16)') '    olo:',rslt(0)
 
149
      write(6,'(a8,2d24.16)') '    olo:',b0(0)
 
150
      write(6,'(a8,2d24.16)') 'b1   LT:',b1LT !|LT=yes
 
151
      write(6,'(a8,2d24.16)') '    olo:',b1(0)
 
152
      write(6,'(a8,2d24.16)') 'b00  LT:',b00LT !|LT=yes
 
153
      write(6,'(a8,2d24.16)') '    olo:',b00(0)
 
154
      write(6,'(a8,2d24.16)') 'b11  LT:',b11LT !|LT=yes
 
155
      write(6,'(a8,2d24.16)') '    olo:',b11(0)
 
156
      write(6,'(a10)') ' eps^(-1):'
 
157
      write(6,'(a8,2d24.16)') '    olo:',b0(1)
 
158
      write(6,'(a10)') ' eps^(-2):'
 
159
      write(6,'(a8,2d24.16)') '    olo:',b0(2)
 
160
      write(6,*)
 
161
!
 
162
      end