156 SUBROUTINE sorbdb5( M1, M2, N, X1, INCX1, X2, INCX2, Q1, LDQ1, Q2,
157 $ ldq2, work, lwork, info )
165 INTEGER incx1, incx2, info, ldq1, ldq2, lwork, m1, m2,
169 REAL q1(ldq1,*), q2(ldq2,*), work(*), x1(*), x2(*)
176 parameter( one = 1.0e0, zero = 0.0e0 )
179 INTEGER childinfo, i,
j
198 ELSE IF( m2 .LT. 0 )
THEN
200 ELSE IF( n .LT. 0 )
THEN
202 ELSE IF( incx1 .LT. 1 )
THEN
204 ELSE IF( incx2 .LT. 1 )
THEN
206 ELSE IF( ldq1 .LT. max( 1, m1 ) )
THEN
208 ELSE IF( ldq2 .LT. max( 1, m2 ) )
THEN
210 ELSE IF( lwork .LT. n )
THEN
214 IF( info .NE. 0 )
THEN
215 CALL
xerbla(
'SORBDB5', -info )
221 CALL
sorbdb6( m1, m2, n, x1, incx1, x2, incx2, q1, ldq1, q2, ldq2,
222 $ work, lwork, childinfo )
226 IF(
snrm2(m1,x1,incx1) .NE. zero
227 $ .OR.
snrm2(m2,x2,incx2) .NE. zero )
THEN
242 CALL
sorbdb6( m1, m2, n, x1, incx1, x2, incx2, q1, ldq1, q2,
243 $ ldq2, work, lwork, childinfo )
244 IF(
snrm2(m1,x1,incx1) .NE. zero
245 $ .OR.
snrm2(m2,x2,incx2) .NE. zero )
THEN
261 CALL
sorbdb6( m1, m2, n, x1, incx1, x2, incx2, q1, ldq1, q2,
262 $ ldq2, work, lwork, childinfo )
263 IF(
snrm2(m1,x1,incx1) .NE. zero
264 $ .OR.
snrm2(m2,x2,incx2) .NE. zero )
THEN
subroutine sorbdb5(M1, M2, N, X1, INCX1, X2, INCX2, Q1, LDQ1, Q2, LDQ2, WORK, LWORK, INFO)
SORBDB5
subroutine xerbla(SRNAME, INFO)
XERBLA
real function snrm2(N, X, INCX)
SNRM2
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine sorbdb6(M1, M2, N, X1, INCX1, X2, INCX2, Q1, LDQ1, Q2, LDQ2, WORK, LWORK, INFO)
SORBDB6