1 SUBROUTINE ssubsp(NMAX, N, A, B, Z, FTEST, EPS, NDIM, FAIL, IND)
2 INTEGER NMAX, N, FTEST, NDIM, IND(n)
4 REAL 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 REAL 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
sexchqz(nmax, n, a, b, z, l, ls1, ls2, eps, fail)
97 ind(ifirst+1) = ind(ifirst)
100 90 l2i = l2i + ind(i)
subroutine sexchqz(NMAX, N, A, B, Z, L, LS1, LS2, EPS, FAIL)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
subroutine ssubsp(NMAX, N, A, B, Z, FTEST, EPS, NDIM, FAIL, IND)