1 SUBROUTINE zkscl(ZRR,ZRI,FNU,N,YR,YI,NZ,RZR,RZI,ASCLE,TOL,ELIM)
12 DOUBLE PRECISION ACS, AS, ASCLE, CKI, CKR, CSI, CSR, CYI,
13 * cyr, elim, fn, fnu, rzi, rzr,
str, s1i, s1r, s2i,
14 * s2r, tol, yi, yr, zeroi, zeror, zri, zrr, xzabs,
15 * zdr, zdi, celmr, elm, helim, alas
16 INTEGER I, IC, IDUM, KK, N, NN, NW, NZ
18 DATA zeror,zeroi / 0.0d0 , 0.0d0 /
33 IF (acs.LT.(-elim)) go
to 10
34 CALL
xzlog(s1r, s1i, csr, csi, idum)
40 CALL
zuchk(csr, csi, nw, ascle, tol)
75 s2r = ckr*csr - cki*csi + s1r
76 s2i = cki*csr + ckr*csi + s1i
87 IF (acs.LT.(-elim)) go
to 25
88 CALL
xzlog(s2r, s2i, csr, csi, idum)
94 CALL
zuchk(csr, csi, nw, ascle, tol)
99 IF (ic.EQ.kk-1) go
to 40
103 IF(alas.LT.helim) go
to 30
subroutine xzlog(AR, AI, BR, BI, IERR)
subroutine zkscl(ZRR, ZRI, FNU, N, YR, YI, NZ, RZR, RZI, ASCLE, TOL, ELIM)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)
subroutine zuchk(YR, YI, NZ, ASCLE, TOL)