~maddevelopers/mg5amcnlo/2.9.4

« back to all changes in this revision

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

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

Show diffs side-by-side

added

removed

vendor/CutTools/src/qcdloop/qlbox8.f

 

1

      subroutine qlbox8(Y,musq,Ires)

 

2

      implicit none

 

3

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

 

4

 

 

5

c                 [                            s12     msq - p4sq ]

 

6

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

 

7

c                 [                             2          2      ]

 

8

c                 [                                               ]

 

9

c                 [                                    msq - s23  ]

 

10

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

 

11

c                 [                                        2      ]

 

12

c            y8 = [                                               ]

 

13

c                 [     s12                            msq - p3sq ]

 

14

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

 

15

c                 [      2                                 2      ]

 

16

c                 [                                               ]

 

17

c                 [ msq - p4sq  msq - s23  msq - p3sq             ]

 

18

c                 [ ----------  ---------  ----------     msq     ]

 

19

c                 [     2           2          2                  ]

 

20

      include 'qlconstants.f'

 

21

      integer iep

 

22

      double precision msq,si,p3sqbar,p4sqbar,tabar,musq,Y(4,4),

 

23

     .     r1,r2,ieps1,ieps2

 

24

      double complex Ires(-2:0),qllnrat,wlogs,qlLi2omrat,qlLi2omx2,

 

25

     . wlogp3,wlogp4,dilog34,dilog4,dilog3

 

26

 

 

27

      msq=Y(4,4)

 

28

      tabar=2d0*Y(2,4)

 

29

      si=2d0*Y(1,3)

 

30

      p3sqbar=2d0*Y(3,4)

 

31

      p4sqbar=2d0*Y(1,4)

 

32

      wlogs=qllnrat(si,musq)

 

33

      wlogp3=qllnrat(p3sqbar,tabar)

 

34

      wlogp4=qllnrat(p4sqbar,tabar)

 

35

 

 

36

      dilog3=qlLi2omrat(p3sqbar,tabar)

 

37

      dilog4=qlLi2omrat(p4sqbar,tabar)

 

38

      dilog34=qlLi2omx2(p3sqbar,p4sqbar,si,msq)

 

39

 

 

40

      call qlratreal(p3sqbar,si,r1,ieps1)

 

41

      call qlratreal(p4sqbar,msq,r2,ieps2)

 

42

       

 

43

      Ires(-2)=cone

 

44

      Ires(-1)=wlogp3+wlogp4-wlogs

 

45

 

 

46

      Ires( 0)=-2d0*dilog3-2d0*dilog4-dilog34

 

47

     . -dcmplx(pisq/6d0)+0.5d0*(qllnrat(si,musq)**2-qllnrat(si,msq)**2)

 

48

     . +2d0*qllnrat(si,musq)*qllnrat(tabar,msq)

 

49

     . -qllnrat(p3sqbar,musq)*qllnrat(p3sqbar,msq)

 

50

     . -qllnrat(p4sqbar,musq)*qllnrat(p4sqbar,msq)

 

51

      do iep=-2,0

 

52

      Ires(iep)=Ires(iep)/dcmplx(si*tabar)

 

53

      enddo

 

54

      return

 

55

      end

 

56

 

• Older »

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