]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | |
2 | *================================================================== | |
3 | ||
4 | FUNCTION RKZPR(L,Q1,Q2) | |
5 | IMPLICIT DOUBLE PRECISION(A-H,O-Z) | |
6 | COMPLEX RKZPR | |
7 | DIMENSION Q1(0:4),Q2(0:4) | |
8 | IF(IABS(L).NE.1) THEN | |
9 | WRITE(6,*) ' RKZPR: ERROR L=',L | |
10 | STOP | |
11 | ENDIF | |
12 | C...Introduce cutoff to check that R1 and R2 not zero. | |
13 | R1=DSQRT(MAX(1D-10,Q1(0)-Q1(1))) | |
14 | R2=DSQRT(MAX(1D-10,Q2(0)-Q2(1))) | |
15 | RKZPR=CMPLX(SNGL(Q1(2)),SNGL(Q1(3)))*R2/R1 | |
16 | . -CMPLX(SNGL(Q2(2)),SNGL(Q2(3)))*R1/R2 | |
17 | IF(L.EQ.-1) RKZPR=-CONJG(RKZPR) | |
18 | END |