Skip to content

Adds CRSCL #839

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 5 commits into from
Jul 4, 2023
Merged

Adds CRSCL #839

Changes from 1 commit
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
831ac96
Adds CRSCL
weslleyspereira May 24, 2023
59a1772
(Last?) version of RSCL
weslleyspereira Jun 1, 2023
bfcece3
Adding zrscl.f
weslleyspereira Jun 26, 2023
c9a49d9
Adds zrscl in the build environments
weslleyspereira Jun 26, 2023
d212879
Fix bugs in ZRSCL
weslleyspereira Jun 26, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Fix bugs in ZRSCL
  • Loading branch information
@weslleyspereira
weslleyspereira committed Jun 26, 2023
commit d212879672abb6f4a1d14711b742da2dad785e6c
8 changes: 5 additions & 3 deletions SRC/zrscl.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -105,7 +105,7 @@ SUBROUTINE ZRSCL( N, A, X, INCX )
* ..
* .. External Functions ..
DOUBLE PRECISION DLAMCH
COMPLEX ZLADIV
COMPLEX*16 ZLADIV
EXTERNAL DLAMCH, ZLADIV
* ..
* .. External Subroutines ..
Expand Down Expand Up @@ -164,7 +164,8 @@ SUBROUTINE ZRSCL( N, A, X, INCX )
*
IF( (ABS( UR ).LT.SAFMIN).OR.(ABS( UI ).LT.SAFMIN) ) THEN
* This means that both alphaR and alphaI are very small.
CALL ZSCAL( N, DCMPLX( SAFMIN / UR, -SAFMIN / UI ), X, INCX )
CALL ZSCAL( N, DCMPLX( SAFMIN / UR, -SAFMIN / UI ), X,
$ INCX )
CALL ZDSCAL( N, SAFMAX, X, INCX )
ELSE IF( (ABS( UR ).GT.SAFMAX).OR.(ABS( UI ).GT.SAFMAX) ) THEN
IF( (ABSR.GT.OV).OR.(ABSI.GT.OV) ) THEN
Expand All @@ -183,7 +184,8 @@ SUBROUTINE ZRSCL( N, A, X, INCX )
UR = (SAFMIN * AR) + AI * ( (SAFMIN * AI) / AR )
UI = (SAFMIN * AI) + SAFMIN * (AR * ( AR / AI ))
END IF
CALL ZSCAL( N, DCMPLX( ONE / UR, -ONE / UI ), X, INCX )
CALL ZSCAL( N, DCMPLX( ONE / UR, -ONE / UI ), X,
$ INCX )
ELSE
CALL ZSCAL( N, DCMPLX( SAFMAX / UR, -SAFMAX / UI ),
$ X, INCX )
Expand Down