Skip to content

Commit d212879

Browse files
Fix bugs in ZRSCL
1 parent c9a49d9 commit d212879

File tree

1 file changed

+5
-3
lines changed

1 file changed

+5
-3
lines changed

SRC/zrscl.f

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -105,7 +105,7 @@ SUBROUTINE ZRSCL( N, A, X, INCX )
105105
* ..
106106
* .. External Functions ..
107107
DOUBLE PRECISION DLAMCH
108-
COMPLEX ZLADIV
108+
COMPLEX*16 ZLADIV
109109
EXTERNAL DLAMCH, ZLADIV
110110
* ..
111111
* .. External Subroutines ..
@@ -164,7 +164,8 @@ SUBROUTINE ZRSCL( N, A, X, INCX )
164164
*
165165
IF( (ABS( UR ).LT.SAFMIN).OR.(ABS( UI ).LT.SAFMIN) ) THEN
166166
* This means that both alphaR and alphaI are very small.
167-
CALL ZSCAL( N, DCMPLX( SAFMIN / UR, -SAFMIN / UI ), X, INCX )
167+
CALL ZSCAL( N, DCMPLX( SAFMIN / UR, -SAFMIN / UI ), X,
168+
$ INCX )
168169
CALL ZDSCAL( N, SAFMAX, X, INCX )
169170
ELSE IF( (ABS( UR ).GT.SAFMAX).OR.(ABS( UI ).GT.SAFMAX) ) THEN
170171
IF( (ABSR.GT.OV).OR.(ABSI.GT.OV) ) THEN
@@ -183,7 +184,8 @@ SUBROUTINE ZRSCL( N, A, X, INCX )
183184
UR = (SAFMIN * AR) + AI * ( (SAFMIN * AI) / AR )
184185
UI = (SAFMIN * AI) + SAFMIN * (AR * ( AR / AI ))
185186
END IF
186-
CALL ZSCAL( N, DCMPLX( ONE / UR, -ONE / UI ), X, INCX )
187+
CALL ZSCAL( N, DCMPLX( ONE / UR, -ONE / UI ), X,
188+
$ INCX )
187189
ELSE
188190
CALL ZSCAL( N, DCMPLX( SAFMAX / UR, -SAFMAX / UI ),
189191
$ X, INCX )

0 commit comments

Comments
 (0)