1 SUBROUTINE zs1s2(ZRR, ZRI, S1R, S1I, S2R, S2I, NZ, ASCLE, ALIM,
17 DOUBLE PRECISION AA, ALIM, ALN, ASCLE, AS1, AS2, C1I, C1R, S1DI,
18 * s1dr, s1i, s1r, s2i, s2r, zeroi, zeror, zri, zrr, xzabs
20 DATA zeror,zeroi / 0.0d0 , 0.0d0 /
24 IF (s1r.EQ.0.0d0 .AND. s1i.EQ.0.0d0) go
to 10
25 IF (as1.EQ.0.0d0) go
to 10
26 aln = -zrr - zrr + dlog(as1)
32 IF (aln.LT.(-alim)) go
to 10
33 CALL
xzlog(s1dr, s1di, c1r, c1i, idum)
36 CALL
xzexp(c1r, c1i, s1r, s1i)
41 IF (aa.GT.ascle)
RETURN
subroutine xzlog(AR, AI, BR, BI, IERR)
subroutine xzexp(AR, AI, BR, BI)
subroutine zs1s2(ZRR, ZRI, S1R, S1I, S2R, S2I, NZ, ASCLE, ALIM, IUF)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to