5 SUBROUTINE dheqr (A, LDA, N, Q, INFO, IJOB)
67 INTEGER LDA, N, INFO, IJOB
68 DOUBLE PRECISION A(lda,*), Q(*)
69 INTEGER I, IQ, J, K, KM1, KP1, NM1
70 DOUBLE PRECISION C, S, T, T1, T2
72 IF (ijob .GT. 1) go
to 70
88 IF (km1 .LT. 1) go
to 20
96 a(j+1,k) = s*t1 + c*t2
105 IF (t2 .NE. 0.0d0) go
to 30
112 s = -1.0d0/
sqrt(1.0d0+t*t)
117 c = 1.0d0/
sqrt(1.0d0+t*t)
123 IF (a(k,k) .EQ. 0.0d0) info = k
143 a(k+1,n) = s*t1 + c*t2
152 IF (t2 .NE. 0.0d0) go
to 110
159 s = -1.0d0/
sqrt(1.0d0+t*t)
164 c = 1.0d0/
sqrt(1.0d0+t*t)
171 IF (a(n,n) .EQ. 0.0d0) info = n
subroutine dheqr(A, LDA, N, Q, INFO, IJOB)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
OCTAVE_EXPORT octave_value_list return the value of the option it must match the dimension of the state and the relative tolerance must also be a vector of the same length tem it must match the dimension of the state and the absolute tolerance must also be a vector of the same length The local error test applied at each integration step is xample roup abs(local error in x(i))<
octave_value sqrt(void) const