~maddevelopers/mg5amcnlo/dim6_eft

« back to all changes in this revision

Viewing changes to vendor/IREGI/src/qcdloop/ql/qlbox2.f

• browse files at revision 255
• compare with another revision
• download diff
• download tarball
• view history from revision 255

Show diffs side-by-side

added

removed

vendor/IREGI/src/qcdloop/ql/qlbox2.f

 

1

      subroutine qlbox2(Y,musq,Ires)

 

2

C $I_4^{D=4-2 \epsilon}(0,0,0,\pq^2;s_{12},s_{23};0,0,0,0)$}

 

3

C     One mass integral as given in egz2, Eq.(A22).

 

4

c                          [                  s12    p4sq ]

 

5

c                          [   0       0    - ---  - ---- ]

 

6

c                          [                   2      2   ]

 

7

c                          [                              ]

 

8

c                          [                         s23  ]

 

9

c                          [   0       0      0    - ---  ]

 

10

c                          [                          2   ]

 

11

c                     Y2 = [                              ]

 

12

c                          [   s12                        ]

 

13

c                          [ - ---     0      0      0    ]

 

14

c                          [    2                         ]

 

15

c                          [                              ]

 

16

c                          [   p4sq    s23                ]

 

17

c                          [ - ----  - ---    0      0    ]

 

18

c                          [    2       2                 ]

 

19

      implicit none

 

20

      double precision musq,si,ta,mp4sq,Y(4,4)

 

21

      double complex Ires(-2:0),qlLsm1,qllnrat,ctwo,fac

 

22

      parameter(ctwo=(2d0,0d0))

 

23

 

 

24

      si=2d0*Y(1,3)

 

25

      ta=2d0*Y(2,4)

 

26

      mp4sq=2d0*Y(1,4)

 

27

      fac=dcmplx(1d0/(si*ta))      

 

28

      Ires(-2)=fac*ctwo

 

29

      Ires(-1)=fac*ctwo

 

30

     . *(qllnrat(mp4sq,musq)

 

31

     .  -qllnrat(ta,musq)

 

32

     .  -qllnrat(si,musq))

 

33

      Ires( 0)=fac*(-qllnrat(mp4sq,musq)**2

 

34

     . +qllnrat(ta,musq)**2

 

35

     . +qllnrat(si,musq)**2

 

36

     . +ctwo*(qlLsm1(ta,mp4sq,si,mp4sq)

 

37

     . -qllnrat(ta,mp4sq)*qllnrat(si,mp4sq))

 

38

     . +qllnrat(mp4sq,ta)**2+qllnrat(mp4sq,si)**2

 

39

     . -qllnrat(ta,si)**2)

 

40

 

 

41

      return

 

42

      end

• Older »

Loggerhead is a web-based interface for Breezy
Version: 2.0.1