LAPACK dlarfg question
We use LAPACK in a Finite Element program for structural mechanics
applications. It has been and still is a task of ours to port it from
Unix to a Linux cluster. However, we've stumbled upon some real
difficulties.
The code currently gets infinitely "stuck" in the dlarfg routine. To be
even more precise, this code is similar to where it enters and never
returns:
BETA = SIGN( DLAPY2( ALPHA, XNORM ), ALPHA )
SAFMIN = DLAMCH( 'S' ) / DLAMCH( 'E' )
IF ( ABS( BETA ).LT.SAFMIN) THEN
DO WHILE ( ABS( BETA ).LT.SAFMIN )
RSAFMN = ONE / SAFMIN
KNT = 0
KNT = KNT + 1
CALL DSCAL( N1, RSAFMN, X, INCX )
BETA = BETA*RSAFMN
ALPHA = ALPHA*RSAFMN
END DO
Beta is forever equal to a small number and SAFMIN to 1. Have you
encountered this problem before and can you offer any advice? Thanks in
advance!!
Regards,
Håkan Lane
Chalmers
(Remove CUT., SPAM. and PLEASE. for email address)
