]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | |
2 | C*********************************************************************** | |
3 | ||
4 | SUBROUTINE PYWAUX(IAUX,EPS,WRE,WIM) | |
5 | ||
6 | C...Calculates real and imaginary parts of the auxiliary functions W1 | |
7 | C...and W2; see R. K. Ellis, I. Hinchliffe, M. Soldate and J. J. van | |
8 | C...der Bij, Nucl. Phys. B297 (1988) 221. | |
9 | COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200) | |
10 | SAVE /LUDAT1/ | |
11 | ||
12 | ASINH(X)=LOG(X+SQRT(X**2+1.)) | |
13 | ACOSH(X)=LOG(X+SQRT(X**2-1.)) | |
14 | ||
15 | IF(EPS.LT.0.) THEN | |
16 | IF(IAUX.EQ.1) WRE=2.*SQRT(1.-EPS)*ASINH(SQRT(-1./EPS)) | |
17 | IF(IAUX.EQ.2) WRE=4.*(ASINH(SQRT(-1./EPS)))**2 | |
18 | WIM=0. | |
19 | ELSEIF(EPS.LT.1.) THEN | |
20 | IF(IAUX.EQ.1) WRE=2.*SQRT(1.-EPS)*ACOSH(SQRT(1./EPS)) | |
21 | IF(IAUX.EQ.2) WRE=4.*(ACOSH(SQRT(1./EPS)))**2-PARU(1)**2 | |
22 | IF(IAUX.EQ.1) WIM=-PARU(1)*SQRT(1.-EPS) | |
23 | IF(IAUX.EQ.2) WIM=-4.*PARU(1)*ACOSH(SQRT(1./EPS)) | |
24 | ELSE | |
25 | IF(IAUX.EQ.1) WRE=2.*SQRT(EPS-1.)*ASIN(SQRT(1./EPS)) | |
26 | IF(IAUX.EQ.2) WRE=-4.*(ASIN(SQRT(1./EPS)))**2 | |
27 | WIM=0. | |
28 | ENDIF | |
29 | ||
30 | RETURN | |
31 | END |