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

This commit is contained in:
Martin Kroeker 2024-10-14 18:57:50 +02:00 committed by GitHub
parent 7018c1b001
commit 27ed6da331
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
1 changed files with 5 additions and 5 deletions

View File

@ -293,7 +293,7 @@
BCOEFR = SCALE*SALFR BCOEFR = SCALE*SALFR
CALL SGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC ), 1, CALL SGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC ), 1,
$ ZERO, WORK( N*( JVEC-1 )+1 ), 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 ) $ 1, ONE, WORK( N*( JVEC-1 )+1 ), 1 )
ELSE ELSE
* *
@ -323,16 +323,16 @@
* *
CALL SGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC ), 1, CALL SGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC ), 1,
$ ZERO, WORK( N*( JVEC-1 )+1 ), 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 ) $ 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 ) $ 1, ONE, WORK( N*( JVEC-1 )+1 ), 1 )
* *
CALL SGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC+1 ), CALL SGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC+1 ),
$ 1, ZERO, WORK( N*JVEC+1 ), 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 ) $ 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 ) $ 1, ONE, WORK( N*JVEC+1 ), 1 )
END IF END IF
END IF END IF