diff options
| author | julie <julielangou@users.noreply.github.com> | 2010-11-07 18:30:58 +0000 |
|---|---|---|
| committer | julie <julielangou@users.noreply.github.com> | 2010-11-07 18:30:58 +0000 |
| commit | 811fbdb838936721133eb83b769027ba50a5c55d (patch) | |
| tree | 8af7b3e42e70127ea1d77f8656cb4aec60b718dc /SRC/zunbdb.f | |
| parent | e9718dbb16036b1ef3684a9d396954a72d3c25a0 (diff) | |
| download | lapack-811fbdb838936721133eb83b769027ba50a5c55d.tar.gz lapack-811fbdb838936721133eb83b769027ba50a5c55d.tar.bz2 lapack-811fbdb838936721133eb83b769027ba50a5c55d.zip | |
Fix a minor problem on csd - intrinsic function is CMPLX and not COMPLEX (same for double precision)
Diffstat (limited to 'SRC/zunbdb.f')
| -rw-r--r-- | SRC/zunbdb.f | 50 |
1 files changed, 25 insertions, 25 deletions
diff --git a/SRC/zunbdb.f b/SRC/zunbdb.f index 08298387..71aa02da 100644 --- a/SRC/zunbdb.f +++ b/SRC/zunbdb.f @@ -208,7 +208,7 @@ * .. * .. Intrinsic Functions INTRINSIC ATAN2, COS, MAX, MIN, SIN - INTRINSIC COMPLEX, DCONJG + INTRINSIC DCMPLX, DCONJG * .. * .. Executable Statements .. * @@ -280,19 +280,19 @@ DO I = 1, Q * IF( I .EQ. 1 ) THEN - CALL ZSCAL( P-I+1, COMPLEX( Z1, 0.0D0 ), X11(I,I), 1 ) + CALL ZSCAL( P-I+1, DCMPLX( Z1, 0.0D0 ), X11(I,I), 1 ) ELSE - CALL ZSCAL( P-I+1, COMPLEX( Z1*COS(PHI(I-1)), 0.0D0 ), + CALL ZSCAL( P-I+1, DCMPLX( Z1*COS(PHI(I-1)), 0.0D0 ), $ X11(I,I), 1 ) - CALL ZAXPY( P-I+1, COMPLEX( -Z1*Z3*Z4*SIN(PHI(I-1)), + CALL ZAXPY( P-I+1, DCMPLX( -Z1*Z3*Z4*SIN(PHI(I-1)), $ 0.0D0 ), X12(I,I-1), 1, X11(I,I), 1 ) END IF IF( I .EQ. 1 ) THEN - CALL ZSCAL( M-P-I+1, COMPLEX( Z2, 0.0D0 ), X21(I,I), 1 ) + CALL ZSCAL( M-P-I+1, DCMPLX( Z2, 0.0D0 ), X21(I,I), 1 ) ELSE - CALL ZSCAL( M-P-I+1, COMPLEX( Z2*COS(PHI(I-1)), 0.0D0 ), + CALL ZSCAL( M-P-I+1, DCMPLX( Z2*COS(PHI(I-1)), 0.0D0 ), $ X21(I,I), 1 ) - CALL ZAXPY( M-P-I+1, COMPLEX( -Z2*Z3*Z4*SIN(PHI(I-1)), + CALL ZAXPY( M-P-I+1, DCMPLX( -Z2*Z3*Z4*SIN(PHI(I-1)), $ 0.0D0 ), X22(I,I-1), 1, X21(I,I), 1 ) END IF * @@ -314,14 +314,14 @@ $ DCONJG(TAUP2(I)), X22(I,I), LDX22, WORK ) * IF( I .LT. Q ) THEN - CALL ZSCAL( Q-I, COMPLEX( -Z1*Z3*SIN(THETA(I)), 0.0D0 ), + CALL ZSCAL( Q-I, DCMPLX( -Z1*Z3*SIN(THETA(I)), 0.0D0 ), $ X11(I,I+1), LDX11 ) - CALL ZAXPY( Q-I, COMPLEX( Z2*Z3*COS(THETA(I)), 0.0D0 ), + CALL ZAXPY( Q-I, DCMPLX( Z2*Z3*COS(THETA(I)), 0.0D0 ), $ X21(I,I+1), LDX21, X11(I,I+1), LDX11 ) END IF - CALL ZSCAL( M-Q-I+1, COMPLEX( -Z1*Z4*SIN(THETA(I)), 0.0D0 ), + CALL ZSCAL( M-Q-I+1, DCMPLX( -Z1*Z4*SIN(THETA(I)), 0.0D0 ), $ X12(I,I), LDX12 ) - CALL ZAXPY( M-Q-I+1, COMPLEX( Z2*Z4*COS(THETA(I)), 0.0D0 ), + CALL ZAXPY( M-Q-I+1, DCMPLX( Z2*Z4*COS(THETA(I)), 0.0D0 ), $ X22(I,I), LDX22, X12(I,I), LDX12 ) * IF( I .LT. Q ) @@ -360,7 +360,7 @@ * DO I = Q + 1, P * - CALL ZSCAL( M-Q-I+1, COMPLEX( -Z1*Z4, 0.0D0 ), X12(I,I), + CALL ZSCAL( M-Q-I+1, DCMPLX( -Z1*Z4, 0.0D0 ), X12(I,I), $ LDX12 ) CALL ZLACGV( M-Q-I+1, X12(I,I), LDX12 ) CALL ZLARFGP( M-Q-I+1, X12(I,I), X12(I,I+1), LDX12, @@ -381,7 +381,7 @@ * DO I = 1, M - P - Q * - CALL ZSCAL( M-P-Q-I+1, COMPLEX( Z2*Z4, 0.0D0 ), + CALL ZSCAL( M-P-Q-I+1, DCMPLX( Z2*Z4, 0.0D0 ), $ X22(Q+I,P+I), LDX22 ) CALL ZLACGV( M-P-Q-I+1, X22(Q+I,P+I), LDX22 ) CALL ZLARFGP( M-P-Q-I+1, X22(Q+I,P+I), X22(Q+I,P+I+1), @@ -401,21 +401,21 @@ DO I = 1, Q * IF( I .EQ. 1 ) THEN - CALL ZSCAL( P-I+1, COMPLEX( Z1, 0.0D0 ), X11(I,I), + CALL ZSCAL( P-I+1, DCMPLX( Z1, 0.0D0 ), X11(I,I), $ LDX11 ) ELSE - CALL ZSCAL( P-I+1, COMPLEX( Z1*COS(PHI(I-1)), 0.0D0 ), + CALL ZSCAL( P-I+1, DCMPLX( Z1*COS(PHI(I-1)), 0.0D0 ), $ X11(I,I), LDX11 ) - CALL ZAXPY( P-I+1, COMPLEX( -Z1*Z3*Z4*SIN(PHI(I-1)), + CALL ZAXPY( P-I+1, DCMPLX( -Z1*Z3*Z4*SIN(PHI(I-1)), $ 0.0D0 ), X12(I-1,I), LDX12, X11(I,I), LDX11 ) END IF IF( I .EQ. 1 ) THEN - CALL ZSCAL( M-P-I+1, COMPLEX( Z2, 0.0D0 ), X21(I,I), + CALL ZSCAL( M-P-I+1, DCMPLX( Z2, 0.0D0 ), X21(I,I), $ LDX21 ) ELSE - CALL ZSCAL( M-P-I+1, COMPLEX( Z2*COS(PHI(I-1)), 0.0D0 ), + CALL ZSCAL( M-P-I+1, DCMPLX( Z2*COS(PHI(I-1)), 0.0D0 ), $ X21(I,I), LDX21 ) - CALL ZAXPY( M-P-I+1, COMPLEX( -Z2*Z3*Z4*SIN(PHI(I-1)), + CALL ZAXPY( M-P-I+1, DCMPLX( -Z2*Z3*Z4*SIN(PHI(I-1)), $ 0.0D0 ), X22(I-1,I), LDX22, X21(I,I), LDX21 ) END IF * @@ -444,14 +444,14 @@ CALL ZLACGV( M-P-I+1, X21(I,I), LDX21 ) * IF( I .LT. Q ) THEN - CALL ZSCAL( Q-I, COMPLEX( -Z1*Z3*SIN(THETA(I)), 0.0D0 ), + CALL ZSCAL( Q-I, DCMPLX( -Z1*Z3*SIN(THETA(I)), 0.0D0 ), $ X11(I+1,I), 1 ) - CALL ZAXPY( Q-I, COMPLEX( Z2*Z3*COS(THETA(I)), 0.0D0 ), + CALL ZAXPY( Q-I, DCMPLX( Z2*Z3*COS(THETA(I)), 0.0D0 ), $ X21(I+1,I), 1, X11(I+1,I), 1 ) END IF - CALL ZSCAL( M-Q-I+1, COMPLEX( -Z1*Z4*SIN(THETA(I)), 0.0D0 ), + CALL ZSCAL( M-Q-I+1, DCMPLX( -Z1*Z4*SIN(THETA(I)), 0.0D0 ), $ X12(I,I), 1 ) - CALL ZAXPY( M-Q-I+1, COMPLEX( Z2*Z4*COS(THETA(I)), 0.0D0 ), + CALL ZAXPY( M-Q-I+1, DCMPLX( Z2*Z4*COS(THETA(I)), 0.0D0 ), $ X22(I,I), 1, X12(I,I), 1 ) * IF( I .LT. Q ) @@ -482,7 +482,7 @@ * DO I = Q + 1, P * - CALL ZSCAL( M-Q-I+1, COMPLEX( -Z1*Z4, 0.0D0 ), X12(I,I), 1 ) + CALL ZSCAL( M-Q-I+1, DCMPLX( -Z1*Z4, 0.0D0 ), X12(I,I), 1 ) CALL ZLARFGP( M-Q-I+1, X12(I,I), X12(I+1,I), 1, TAUQ2(I) ) X12(I,I) = ONE * @@ -498,7 +498,7 @@ * DO I = 1, M - P - Q * - CALL ZSCAL( M-P-Q-I+1, COMPLEX( Z2*Z4, 0.0D0 ), + CALL ZSCAL( M-P-Q-I+1, DCMPLX( Z2*Z4, 0.0D0 ), $ X22(P+I,Q+I), 1 ) CALL ZLARFGP( M-P-Q-I+1, X22(P+I,Q+I), X22(P+I+1,Q+I), 1, $ TAUQ2(P+I) ) |