94 SUBROUTINE sget10( M, N, A, LDA, B, LDB, WORK, RESULT )
102 INTEGER lda, ldb, m, n
106 REAL a( lda, * ), b( ldb, * ), work( * )
113 parameter( one = 1.0e+0, zero = 0.0e+0 )
117 REAL anorm, eps, unfl, wnorm
127 INTRINSIC max, min, real
133 IF( m.LE.0 .OR. n.LE.0 )
THEN
138 unfl =
slamch(
'Safe minimum' )
139 eps =
slamch(
'Precision' )
143 CALL
scopy( m, a( 1, j ), 1, work, 1 )
144 CALL
saxpy( m, -one, b( 1, j ), 1, work, 1 )
145 wnorm = max( wnorm,
sasum( n, work, 1 ) )
148 anorm = max(
slange(
'1', m, n, a, lda, work ), unfl )
150 IF( anorm.GT.wnorm )
THEN
151 result = ( wnorm / anorm ) / ( m*eps )
153 IF( anorm.LT.one )
THEN
154 result = ( min( wnorm, m*anorm ) / anorm ) / ( m*eps )
156 result = min( wnorm / anorm,
REAL( M ) ) / ( m*eps )