99 REAL FUNCTION sqrt11( M, K, A, LDA, TAU, WORK, LWORK )
107 INTEGER k, lda, lwork, m
110 REAL a( lda, * ), tau( * ), work( lwork )
117 parameter( zero = 0.0e0, one = 1.0e0 )
141 IF( lwork.LT.m*m+m )
THEN
142 CALL
xerbla(
'SQRT11', 7 )
151 CALL
slaset(
'Full', m, m, zero, one, work, m )
155 CALL
sorm2r(
'Left',
'No transpose', m, m, k, a, lda, tau, work,
156 $ m, work( m*m+1 ), info )
160 CALL
sorm2r(
'Left',
'Transpose', m, m, k, a, lda, tau, work, m,
161 $ work( m*m+1 ), info )
164 work( (
j-1 )*m+
j ) = work( (
j-1 )*m+
j ) - one
168 $ (
REAL( m )*
slamch(
'Epsilon' ) )
subroutine slaset(UPLO, M, N, ALPHA, BETA, A, LDA)
SLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values...
REAL function slamch(CMACH)
SLAMCH
REAL function sqrt11(M, K, A, LDA, TAU, WORK, LWORK)
SQRT11
subroutine xerbla(SRNAME, INFO)
XERBLA
REAL function slange(NORM, M, N, A, LDA, WORK)
SLANGE returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value ...
subroutine sorm2r(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, INFO)
SORM2R multiplies a general matrix by the orthogonal matrix from a QR factorization determined by sge...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j