Skip to content

Commit

Permalink
Merge pull request #4937 from martin-frbg/lapack1064
Browse files Browse the repository at this point in the history 
Fix leading dimension for B in LAPACK tests for GGEV (Reference-LAPACK PR 1064)
  • Loading branch information
@martin-frbg
martin-frbg authored Oct 14, 2024
2 parents c979c1d + 22628f1 commit 8483a71
Show file tree
Hide file tree
Showing 4 changed files with 12 additions and 12 deletions.
2 changes: 1 addition & 1 deletion lapack-netlib/TESTING/EIG/cget52.f
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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
*
Expand Down
10 changes: 5 additions & 5 deletions lapack-netlib/TESTING/EIG/dget52.f
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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
*
Expand Down Expand Up @@ -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
Expand Down
10 changes: 5 additions & 5 deletions lapack-netlib/TESTING/EIG/sget52.f
View file
Original file line number Diff line number Diff line change
Expand Up @@ -293,7 +293,7 @@ SUBROUTINE SGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHAR,
BCOEFR = SCALE*SALFR
CALL SGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC ), 1,
$ ZERO, WORK( N*( JVEC-1 )+1 ), 1 )
CALL SGEMV( TRANS, N, N, -BCOEFR, B, LDA, E( 1, JVEC ),
CALL SGEMV( TRANS, N, N, -BCOEFR, B, LDB, E( 1, JVEC ),
$ 1, ONE, WORK( N*( JVEC-1 )+1 ), 1 )
ELSE
*
Expand Down Expand Up @@ -323,16 +323,16 @@ SUBROUTINE SGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHAR,
*
CALL SGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC ), 1,
$ ZERO, WORK( N*( JVEC-1 )+1 ), 1 )
CALL SGEMV( TRANS, N, N, -BCOEFR, B, LDA, E( 1, JVEC ),
CALL SGEMV( TRANS, N, N, -BCOEFR, B, LDB, E( 1, JVEC ),
$ 1, ONE, WORK( N*( JVEC-1 )+1 ), 1 )
CALL SGEMV( TRANS, N, N, BCOEFI, B, LDA, E( 1, JVEC+1 ),
CALL SGEMV( TRANS, N, N, BCOEFI, B, LDB, E( 1, JVEC+1 ),
$ 1, ONE, WORK( N*( JVEC-1 )+1 ), 1 )
*
CALL SGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC+1 ),
$ 1, ZERO, WORK( N*JVEC+1 ), 1 )
CALL SGEMV( TRANS, N, N, -BCOEFI, B, LDA, E( 1, JVEC ),
CALL SGEMV( TRANS, N, N, -BCOEFI, B, LDB, E( 1, JVEC ),
$ 1, ONE, WORK( N*JVEC+1 ), 1 )
CALL SGEMV( TRANS, N, N, -BCOEFR, B, LDA, E( 1, JVEC+1 ),
CALL SGEMV( TRANS, N, N, -BCOEFR, B, LDB, E( 1, JVEC+1 ),
$ 1, ONE, WORK( N*JVEC+1 ), 1 )
END IF
END IF
Expand Down
2 changes: 1 addition & 1 deletion lapack-netlib/TESTING/EIG/zget52.f
View file
Original file line number Diff line number Diff line change
Expand Up @@ -257,7 +257,7 @@ SUBROUTINE ZGET52( LEFT, N, A, LDA, B, LDB, E, LDE, ALPHA, BETA,
END IF
CALL ZGEMV( TRANS, N, N, ACOEFF, A, LDA, E( 1, JVEC ), 1,
$ CZERO, WORK( N*( JVEC-1 )+1 ), 1 )
CALL ZGEMV( TRANS, N, N, -BCOEFF, B, LDA, E( 1, JVEC ), 1,
CALL ZGEMV( TRANS, N, N, -BCOEFF, B, LDB, E( 1, JVEC ), 1,
$ CONE, WORK( N*( JVEC-1 )+1 ), 1 )
10 CONTINUE
*
Expand Down

0 comments on commit 8483a71

Please sign in to comment.