Merge pull request #4937 from martin-frbg/lapack1064
Fix leading dimension for B in LAPACK tests for GGEV (Reference-LAPACK PR 1064)
This commit is contained in:
commit
8483a71169
|
@ -256,7 +256,7 @@
|
||||||
END IF
|
END IF
|
||||||
CALL CGEMV( TRANS, N, N, ACOEFF, A, LDA, E( 1, JVEC ), 1,
|
CALL CGEMV( TRANS, N, N, ACOEFF, A, LDA, E( 1, JVEC ), 1,
|
||||||
$ CZERO, WORK( N*( JVEC-1 )+1 ), 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 )
|
$ CONE, WORK( N*( JVEC-1 )+1 ), 1 )
|
||||||
10 CONTINUE
|
10 CONTINUE
|
||||||
*
|
*
|
||||||
|
|
|
@ -293,7 +293,7 @@
|
||||||
BCOEFR = SCALE*SALFR
|
BCOEFR = SCALE*SALFR
|
||||||
CALL DGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC ), 1,
|
CALL DGEMV( 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 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 )
|
$ 1, ONE, WORK( N*( JVEC-1 )+1 ), 1 )
|
||||||
ELSE
|
ELSE
|
||||||
*
|
*
|
||||||
|
@ -323,16 +323,16 @@
|
||||||
*
|
*
|
||||||
CALL DGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC ), 1,
|
CALL DGEMV( 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 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 )
|
$ 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 )
|
$ 1, ONE, WORK( N*( JVEC-1 )+1 ), 1 )
|
||||||
*
|
*
|
||||||
CALL DGEMV( TRANS, N, N, ACOEF, A, LDA, E( 1, JVEC+1 ),
|
CALL DGEMV( 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 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 )
|
$ 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 )
|
$ 1, ONE, WORK( N*JVEC+1 ), 1 )
|
||||||
END IF
|
END IF
|
||||||
END IF
|
END IF
|
||||||
|
|
|
@ -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
|
||||||
|
|
|
@ -257,7 +257,7 @@
|
||||||
END IF
|
END IF
|
||||||
CALL ZGEMV( TRANS, N, N, ACOEFF, A, LDA, E( 1, JVEC ), 1,
|
CALL ZGEMV( TRANS, N, N, ACOEFF, A, LDA, E( 1, JVEC ), 1,
|
||||||
$ CZERO, WORK( N*( JVEC-1 )+1 ), 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 )
|
$ CONE, WORK( N*( JVEC-1 )+1 ), 1 )
|
||||||
10 CONTINUE
|
10 CONTINUE
|
||||||
*
|
*
|
||||||
|
|
Loading…
Reference in New Issue