From 7018c1b001bd79dff9d7af842b7f48c66ec03dcb Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Mon, 14 Oct 2024 18:56:44 +0200 Subject: [PATCH] Fix leading dimension for B (Reference-LAPACK PR 1064) --- lapack-netlib/TESTING/EIG/dget52.f | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/lapack-netlib/TESTING/EIG/dget52.f b/lapack-netlib/TESTING/EIG/dget52.f index 68196f5f3b..b662bb6945 100644 --- a/lapack-netlib/TESTING/EIG/dget52.f +++ b/lapack-netlib/TESTING/EIG/dget52.f @@ -293,7 +293,7 @@ SUBROUTINE DGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHAR, BCOEFR = SCALE*SALFR CALL DGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC ), 1, $ ZERO, WORK( N*( JVEC-1 )+1 ), 1 ) - CALL DGEMV( TRANS, N, N, -BCOEFR, B, LDA, E( 1, JVEC ), + CALL DGEMV( TRANS, N, N, -BCOEFR, B, LDB, E( 1, JVEC ), $ 1, ONE, WORK( N*( JVEC-1 )+1 ), 1 ) ELSE * @@ -323,16 +323,16 @@ SUBROUTINE DGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHAR, * CALL DGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC ), 1, $ ZERO, WORK( N*( JVEC-1 )+1 ), 1 ) - CALL DGEMV( TRANS, N, N, -BCOEFR, B, LDA, E( 1, JVEC ), + CALL DGEMV( TRANS, N, N, -BCOEFR, B, LDB, E( 1, JVEC ), $ 1, ONE, WORK( N*( JVEC-1 )+1 ), 1 ) - CALL DGEMV( TRANS, N, N, BCOEFI, B, LDA, E( 1, JVEC+1 ), + CALL DGEMV( TRANS, N, N, BCOEFI, B, LDB, E( 1, JVEC+1 ), $ 1, ONE, WORK( N*( JVEC-1 )+1 ), 1 ) * CALL DGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC+1 ), $ 1, ZERO, WORK( N*JVEC+1 ), 1 ) - CALL DGEMV( TRANS, N, N, -BCOEFI, B, LDA, E( 1, JVEC ), + CALL DGEMV( TRANS, N, N, -BCOEFI, B, LDB, E( 1, JVEC ), $ 1, ONE, WORK( N*JVEC+1 ), 1 ) - CALL DGEMV( TRANS, N, N, -BCOEFR, B, LDA, E( 1, JVEC+1 ), + CALL DGEMV( TRANS, N, N, -BCOEFR, B, LDB, E( 1, JVEC+1 ), $ 1, ONE, WORK( N*JVEC+1 ), 1 ) END IF END IF