1 SUBROUTINE zwrsk(ZRR, ZRI, FNU, KODE, N, YR, YI, NZ, CWR, CWI,
12 DOUBLE PRECISION ACT, ACW, ALIM, ASCLE, CINUI, CINUR, CSCLR, CTI,
13 * ctr, cwi, cwr, c1i, c1r, c2i, c2r, elim, fnu, pti, ptr, ract,
14 * sti,
str, tol, yi, yr, zri, zrr, xzabs,
d1mach
15 INTEGER I, KODE, N, NW, NZ
23 CALL
zbknu(zrr, zri, fnu, kode, 2, cwr, cwi, nw, tol, elim, alim)
25 CALL
zrati(zrr, zri, fnu, n, yr, yi, tol)
32 IF (kode.EQ.1) go
to 10
42 acw = xzabs(cwr(2),cwi(2))
45 IF (acw.GT.ascle) go
to 20
50 IF (acw.LT.ascle) go
to 30
63 ptr =
str*c1r - sti*c1i
64 pti =
str*c1i + sti*c1r
67 ctr = zrr*ptr - zri*pti
68 cti = zrr*pti + zri*ptr
75 cinur = ptr*ctr - pti*cti
76 cinui = ptr*cti + pti*ctr
81 ptr =
str*cinur - sti*cinui
82 cinui =
str*cinui + sti*cinur
subroutine zrati(ZR, ZI, FNU, N, CYR, CYI, TOL)
subroutine zwrsk(ZRR, ZRI, FNU, KODE, N, YR, YI, NZ, CWR, CWI, TOL, ELIM, ALIM)
double precision function d1mach(i)
F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T const F77_REAL const F77_REAL F77_REAL &F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T const F77_DBLE F77_DBLE &F77_RET_T const F77_REAL F77_REAL &F77_RET_T F77_REAL F77_REAL &F77_RET_T F77_DBLE F77_DBLE &F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE * d
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 zbknu(ZR, ZI, FNU, KODE, N, YR, YI, NZ, TOL, ELIM, ALIM)