1 SUBROUTINE dsubsp(NMAX, N, A, B, Z, FTEST, EPS, NDIM, FAIL, IND)
2 INTEGER NMAX, N, FTEST, NDIM, IND(n)
4 DOUBLE PRECISION A(nmax,n), B(nmax,n), Z(nmax,n), EPS
37 INTEGER L, LS, LS1, LS2, L1, LL, NUM, IS, L2I, L2K, I, K, II,
39 DOUBLE PRECISION S, P, D, ALPHA, BETA
55 IF (a(l1,l).EQ.0.) go
to 10
59 s = (a(l,l)*b(l1,l1)+a(l1,l1)*b(l,l)-a(l1,l)*b(l,l1))/d
60 p = (a(l,l)*a(l1,l1)-a(l,l1)*a(l1,l))/d
61 is = ftest(ls,alpha,beta,s,p)
65 is = ftest(ls,a(l,l),b(l,l),s,p)
67 IF (is.EQ.1) ndim = ndim + ls
74 IF (ind(i).GT.0) go
to 90
79 IF (ind(k).LT.0) go
to 50
95 CALL
exchqz(nmax, n, a, b, z, l, ls1, ls2, eps, fail)
97 ind(ifirst+1) = ind(ifirst)
100 90 l2i = l2i + ind(i)
subroutine exchqz(NMAX, N, A, B, Z, L, LS1, LS2, EPS, FAIL)
subroutine dsubsp(NMAX, N, A, B, Z, FTEST, EPS, NDIM, FAIL, IND)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to