Add the missing CONJG ...
authorjulie <julielangou@users.noreply.github.com>
Sun, 15 Nov 2015 02:50:39 +0000 (02:50 +0000)
committerjulie <julielangou@users.noreply.github.com>
Sun, 15 Nov 2015 02:50:39 +0000 (02:50 +0000)
TESTINGS are all good now!

TESTING/EIG/cdrvbd.f
TESTING/EIG/zdrvbd.f

index 612672df9c788610b70c45a21e0f4d9aff563d4a..16e558dc022698abfc3d8a73f7799d1b25a9f4f0 100644 (file)
 *
                   DO J=1,N
                      DO I=1,N
-                        VTSAV(J,I) = A(I,J)
+                        VTSAV(J,I) = CONJG (A(I,J))
                      END DO
                   END DO
 *
 *
 *                 Do tests 15--18
 *
-* TEST 15 NOT PASSING THE THREASHOLD
-*                  CALL CBDT01( M, N, 0, ASAV, LDA, USAV, LDU, SSAV, E,
-*     $                         VTSAV, LDVT, WORK, RWORK, RESULT( 15 ) )
+                  CALL CBDT01( M, N, 0, ASAV, LDA, USAV, LDU, SSAV, E,
+     $                         VTSAV, LDVT, WORK, RWORK, RESULT( 15 ) )
                   IF( M.NE.0 .AND. N.NE.0 ) THEN
                      CALL CUNT01( 'Columns', M, M, USAV, LDU, WORK,
      $                            LWORK, RWORK, RESULT( 16 ) )
      $            LSWORK = LWORK
                LRWORK = MAX( 7, N + 2*M)
 *
-                  CALL CLACPY( 'F', M, N, ASAV, LDA, VTSAV, LDA )
+                 CALL CLACPY( 'F', M, N, ASAV, LDA, VTSAV, LDA )
                   SRNAMT = 'CGEJSV'
                   CALL CGEJSV( 'G', 'U', 'V', 'R', 'N', 'N',
      &                   M, N, VTSAV, LDA, SSAV, USAV, LDU, A, LDVT,
 *
                   DO 133 J=1,N
                      DO 132 I=1,N
-                        VTSAV(J,I) = A(I,J)
+                        VTSAV(J,I) = CONJG (A(I,J))
   132                END DO
   133             END DO
 *
 *
 *                 Do tests 19--22
 *
-* TEST 19 NOT PASSING THE THREASHOLD
-*                  CALL CBDT01( M, N, 0, ASAV, LDA, USAV, LDU, SSAV, E,
-*     $                         VTSAV, LDVT, WORK, RWORK, RESULT( 19 ) )
+                  CALL CBDT01( M, N, 0, ASAV, LDA, USAV, LDU, SSAV, E,
+     $                         VTSAV, LDVT, WORK, RWORK, RESULT( 19 ) )
                   IF( M.NE.0 .AND. N.NE.0 ) THEN
                      CALL CUNT01( 'Columns', M, M, USAV, LDU, WORK,
      $                            LWORK, RWORK, RESULT( 20 ) )
index f390bc066fe89e0f74f11d5867962a021df868e5..06d90fee31d6165be1d3e6c51b4554b8b1f2b674 100644 (file)
 *
                   DO J=1,N
                      DO I=1,N
-                        VTSAV(J,I) = A(I,J)
+                        VTSAV(J,I) = CONJG (A(I,J))
                      END DO
                   END DO
 *
 *
 *                 Do tests 15--18
 *
-* TEST 15 NOT PASSING THE THREASHOLD
-*                  CALL ZBDT01( M, N, 0, ASAV, LDA, USAV, LDU, SSAV, E,
-*     $                         VTSAV, LDVT, WORK, RWORK, RESULT( 15 ) )
+                  CALL ZBDT01( M, N, 0, ASAV, LDA, USAV, LDU, SSAV, E,
+     $                         VTSAV, LDVT, WORK, RWORK, RESULT( 15 ) )
                   IF( M.NE.0 .AND. N.NE.0 ) THEN
                      CALL ZUNT01( 'Columns', M, M, USAV, LDU, WORK,
      $                            LWORK, RWORK, RESULT( 16 ) )
 *
                   DO 133 J=1,N
                      DO 132 I=1,N
-                        VTSAV(J,I) = A(I,J)
+                        VTSAV(J,I) = CONJG (A(I,J))
   132                END DO
   133             END DO
 *
 *
 *                 Do tests 19--22
 *
-* TEST 19 NOT PASSING THE THREASHOLD
-*                  CALL ZBDT01( M, N, 0, ASAV, LDA, USAV, LDU, SSAV, E,
-*     $                         VTSAV, LDVT, WORK, RWORK, RESULT( 19 ) )
+                  CALL ZBDT01( M, N, 0, ASAV, LDA, USAV, LDU, SSAV, E,
+     $                         VTSAV, LDVT, WORK, RWORK, RESULT( 19 ) )
                   IF( M.NE.0 .AND. N.NE.0 ) THEN
                      CALL ZUNT01( 'Columns', M, M, USAV, LDU, WORK,
      $                            LWORK, RWORK, RESULT( 20 ) )