5 * Revision 1.1.1.1 1996/04/01 15:01:48 mclareni
10 #if defined(CERNLIB_DOUBLE)
11 FUNCTION DWIG9J(A,B,C,P,Q,R,X,Y,Z)
12 IMPLICIT DOUBLE PRECISION (A-H,O-Z)
14 PARAMETER (R1 = 1, HF = R1/2)
27 IF(IA .LT. 0 .OR. IB .LT. 0 .OR. IC .LT. 0 .OR.
28 1 IP .LT. 0 .OR. IQ .LT. 0 .OR. IR .LT. 0 .OR.
29 2 IX .LT. 0 .OR. IY .LT. 0 .OR. IZ .LT. 0) GO TO 99
30 J0=MAX(ABS(IA-IZ),ABS(IY-IP),ABS(IB-IR))
31 J1=MIN(IA+IZ,IY+IP,IB+IR)
36 H=DWIG6J(A,P,X,Y,Z,AJ)
37 IF(H .NE. 0) H=H*DWIG6J(B,Q,Y,P,AJ,R)
38 IF(H .NE. 0) H=H*DWIG6J(C,R,Z,AJ,A,B)