Skip to content

Commit c1b96f8

Browse files
committed
Fix double complex cases
1 parent bd02149 commit c1b96f8

File tree

2 files changed

+7
-4
lines changed

2 files changed

+7
-4
lines changed

SRC/zlaic1.f

Lines changed: 5 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -252,7 +252,9 @@ SUBROUTINE ZLAIC1( JOB, J, X, SEST, W, GAMMA, SESTPR, S, C )
252252
*
253253
SINE = -( ALPHA / ABSEST ) / T
254254
COSINE = -( GAMMA / ABSEST ) / ( ONE+T )
255-
TMP = DBLE( SQRT( SINE*DCONJG( SINE )+COSINE*DCONJG( COSINE ) ) )
255+
TMP = DBLE( SQRT( SINE * DCONJG( SINE )
256+
$ + COSINE * DCONJG( COSINE ) ) )
257+
256258
S = SINE / TMP
257259
C = COSINE / TMP
258260
SESTPR = SQRT( T+ONE )*ABSEST
@@ -354,7 +356,8 @@ SUBROUTINE ZLAIC1( JOB, J, X, SEST, W, GAMMA, SESTPR, S, C )
354356
COSINE = -( GAMMA / ABSEST ) / ( ONE+T )
355357
SESTPR = SQRT( ONE+T+FOUR*EPS*EPS*NORMA )*ABSEST
356358
END IF
357-
TMP = DBLE( SQRT( SINE*DCONJG( SINE )+COSINE*DCONJG( COSINE ) ) )
359+
TMP = DBLE( SQRT( SINE * DCONJG( SINE )
360+
$ + COSINE * DCONJG( COSINE ) ) )
358361
S = SINE / TMP
359362
C = COSINE / TMP
360363
RETURN

SRC/zpptrf.f

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -189,8 +189,8 @@ SUBROUTINE ZPPTRF( UPLO, N, AP, INFO )
189189
*
190190
* Compute U(J,J) and test for non-positive-definiteness.
191191
*
192-
AJJ = DBLE( AP( JJ ) ) - DBLE( ZDOTC( J-1, AP( JC ), 1, AP( JC ),
193-
$ 1 ) )
192+
AJJ = DBLE( AP( JJ ) ) - DBLE( ZDOTC( J-1,
193+
$ AP( JC ), 1, AP( JC ), 1 ) )
194194
IF( AJJ.LE.ZERO ) THEN
195195
AP( JJ ) = AJJ
196196
GO TO 30

0 commit comments

Comments
 (0)