From 5f269f9f643f87ec2240d459fa650b0a08289c40 Mon Sep 17 00:00:00 2001 From: julie Date: Sun, 7 Nov 2010 20:15:40 +0000 Subject: [PATCH] fix again some COMPLEX instead of CMPLX --- TESTING/EIG/ccsdts.f | 8 ++++---- TESTING/EIG/zcsdts.f | 8 ++++---- 2 files changed, 8 insertions(+), 8 deletions(-) diff --git a/TESTING/EIG/ccsdts.f b/TESTING/EIG/ccsdts.f index 1a4a798..897f4ff 100644 --- a/TESTING/EIG/ccsdts.f +++ b/TESTING/EIG/ccsdts.f @@ -180,7 +180,7 @@ END DO DO I = 1, R X(MIN(P,Q)-R+I,MIN(P,Q)-R+I) = - $ X(MIN(P,Q)-R+I,MIN(P,Q)-R+I) - COMPLEX( COS(THETA(I)), + $ X(MIN(P,Q)-R+I,MIN(P,Q)-R+I) - CMPLX( COS(THETA(I)), $ 0.0E0 ) END DO * @@ -196,7 +196,7 @@ DO I = 1, R X(P-(MIN(P,M-Q)-R)+1-I,M-(MIN(P,M-Q)-R)+1-I) = $ X(P-(MIN(P,M-Q)-R)+1-I,M-(MIN(P,M-Q)-R)+1-I) + - $ COMPLEX( SIN(THETA(R-I+1)), 0.0E0 ) + $ CMPLX( SIN(THETA(R-I+1)), 0.0E0 ) END DO * CALL CGEMM( 'No transpose', 'Conjugate transpose', M-P, Q, Q, ONE, @@ -211,7 +211,7 @@ DO I = 1, R X(M-(MIN(M-P,Q)-R)+1-I,Q-(MIN(M-P,Q)-R)+1-I) = $ X(M-(MIN(M-P,Q)-R)+1-I,Q-(MIN(M-P,Q)-R)+1-I) - - $ COMPLEX( SIN(THETA(R-I+1)), 0.0E0 ) + $ CMPLX( SIN(THETA(R-I+1)), 0.0E0 ) END DO * CALL CGEMM( 'No transpose', 'Conjugate transpose', M-P, M-Q, M-Q, @@ -226,7 +226,7 @@ DO I = 1, R X(P+(MIN(M-P,M-Q)-R)+I,Q+(MIN(M-P,M-Q)-R)+I) = $ X(P+(MIN(M-P,M-Q)-R)+I,Q+(MIN(M-P,M-Q)-R)+I) - - $ COMPLEX( COS(THETA(I)), 0.0E0 ) + $ CMPLX( COS(THETA(I)), 0.0E0 ) END DO * * Compute norm( U1'*X11*V1 - D11 ) / ( MAX(1,P,Q)*EPS2 ) . diff --git a/TESTING/EIG/zcsdts.f b/TESTING/EIG/zcsdts.f index b185c17..408e230 100644 --- a/TESTING/EIG/zcsdts.f +++ b/TESTING/EIG/zcsdts.f @@ -180,7 +180,7 @@ END DO DO I = 1, R X(MIN(P,Q)-R+I,MIN(P,Q)-R+I) = - $ X(MIN(P,Q)-R+I,MIN(P,Q)-R+I) - COMPLEX( COS(THETA(I)), + $ X(MIN(P,Q)-R+I,MIN(P,Q)-R+I) - DCMPLX( COS(THETA(I)), $ 0.0D0 ) END DO * @@ -196,7 +196,7 @@ DO I = 1, R X(P-(MIN(P,M-Q)-R)+1-I,M-(MIN(P,M-Q)-R)+1-I) = $ X(P-(MIN(P,M-Q)-R)+1-I,M-(MIN(P,M-Q)-R)+1-I) + - $ COMPLEX( SIN(THETA(R-I+1)), 0.0D0 ) + $ DCMPLX( SIN(THETA(R-I+1)), 0.0D0 ) END DO * CALL ZGEMM( 'No transpose', 'Conjugate transpose', M-P, Q, Q, ONE, @@ -211,7 +211,7 @@ DO I = 1, R X(M-(MIN(M-P,Q)-R)+1-I,Q-(MIN(M-P,Q)-R)+1-I) = $ X(M-(MIN(M-P,Q)-R)+1-I,Q-(MIN(M-P,Q)-R)+1-I) - - $ COMPLEX( SIN(THETA(R-I+1)), 0.0D0 ) + $ DCMPLX( SIN(THETA(R-I+1)), 0.0D0 ) END DO * CALL ZGEMM( 'No transpose', 'Conjugate transpose', M-P, M-Q, M-Q, @@ -226,7 +226,7 @@ DO I = 1, R X(P+(MIN(M-P,M-Q)-R)+I,Q+(MIN(M-P,M-Q)-R)+I) = $ X(P+(MIN(M-P,M-Q)-R)+I,Q+(MIN(M-P,M-Q)-R)+I) - - $ COMPLEX( COS(THETA(I)), 0.0D0 ) + $ DCMPLX( COS(THETA(I)), 0.0D0 ) END DO * * Compute norm( U1'*X11*V1 - D11 ) / ( MAX(1,P,Q)*EPS2 ) . -- 2.7.4