Fix leading dimension for B (Reference-LAPACK PR 1064)

martin-frbg · web-flow · commit a659f40fe116 · 2024-10-14T18:53:30.000+02:00
diff --git a/lapack-netlib/TESTING/EIG/cget52.f b/lapack-netlib/TESTING/EIG/cget52.f
@@ -256,7 +256,7 @@ SUBROUTINE CGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHA, BETA,
          END IF
          CALL CGEMV( TRANS, N, N, ACOEFF, A, LDA, E( 1, JVEC ), 1,
      $               CZERO, WORK( N*( JVEC-1 )+1 ), 1 )
-         CALL CGEMV( TRANS, N, N, -BCOEFF, B, LDA, E( 1, JVEC ), 1,
+         CALL CGEMV( TRANS, N, N, -BCOEFF, B, LDB, E( 1, JVEC ), 1,
      $               CONE, WORK( N*( JVEC-1 )+1 ), 1 )
    10 CONTINUE
 *

Original file line number	Diff line number	Diff line change
`@@ -256,7 +256,7 @@ SUBROUTINE CGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHA, BETA,`
`256`	`256`	`END IF`
`257`	`257`	`CALL CGEMV( TRANS, N, N, ACOEFF, A, LDA, E( 1, JVEC ), 1,`
`258`	`258`	`$ CZERO, WORK( N*( JVEC-1 )+1 ), 1 )`
`259`		`- CALL CGEMV( TRANS, N, N, -BCOEFF, B, LDA, E( 1, JVEC ), 1,`
	`259`	`+ CALL CGEMV( TRANS, N, N, -BCOEFF, B, LDB, E( 1, JVEC ), 1,`
`260`	`260`	`$ CONE, WORK( N*( JVEC-1 )+1 ), 1 )`
`261`	`261`	`10 CONTINUE`
`262`	`262`	`*`