Mentions légales du service

Skip to content
Snippets Groups Projects
Commit 312391de authored by Mathieu Faverge's avatar Mathieu Faverge
Browse files

Silent bug issue reletaed to a = a in hetra1 kernels

parent 6d852b8d
No related branches found
No related tags found
1 merge request!87Sonarqube /hotfix
Hide whitespace changes
Inline Side-by-side
coreblas/compute/core_ztsmlq_hetra1.c
+ 4
− 2
View file @ 312391de
...@@ -143,8 +143,9 @@ int CORE_ztsmlq_hetra1( MORSE_enum side, MORSE_enum trans, ...@@ -143,8 +143,9 @@ int CORE_ztsmlq_hetra1( MORSE_enum side, MORSE_enum trans,
/* in-place transposition of A1 */ /* in-place transposition of A1 */
for (j = 0; j < n1; j++){ for (j = 0; j < n1; j++){
#if defined(PRECISION_z) || defined(PRECISION_c)
A1[j + j*lda1] = conj(A1[j + j*lda1]); A1[j + j*lda1] = conj(A1[j + j*lda1]);
#endif
for (i = j+1; i < m1; i++){ for (i = j+1; i < m1; i++){
*WORK = *(A1 + i + j*lda1); *WORK = *(A1 + i + j*lda1);
*(A1 + i + j*lda1) = conj(*(A1 + j + i*lda1)); *(A1 + i + j*lda1) = conj(*(A1 + j + i*lda1));
...@@ -159,8 +160,9 @@ int CORE_ztsmlq_hetra1( MORSE_enum side, MORSE_enum trans, ...@@ -159,8 +160,9 @@ int CORE_ztsmlq_hetra1( MORSE_enum side, MORSE_enum trans,
/* in-place transposition of A1 */ /* in-place transposition of A1 */
for (j = 0; j < n1; j++){ for (j = 0; j < n1; j++){
#if defined(PRECISION_z) || defined(PRECISION_c)
A1[j + j*lda1] = conj(A1[j + j*lda1]); A1[j + j*lda1] = conj(A1[j + j*lda1]);
#endif
for (i = j+1; i < m1; i++){ for (i = j+1; i < m1; i++){
*WORK = *(A1 + i + j*lda1); *WORK = *(A1 + i + j*lda1);
*(A1 + i + j*lda1) = conj(*(A1 + j + i*lda1)); *(A1 + i + j*lda1) = conj(*(A1 + j + i*lda1));
......
coreblas/compute/core_ztsmqr_hetra1.c
+ 4
− 2
View file @ 312391de
...@@ -145,8 +145,9 @@ int CORE_ztsmqr_hetra1( MORSE_enum side, MORSE_enum trans, ...@@ -145,8 +145,9 @@ int CORE_ztsmqr_hetra1( MORSE_enum side, MORSE_enum trans,
/* in-place transposition of A1 */ /* in-place transposition of A1 */
for (j = 0; j < n1; j++){ for (j = 0; j < n1; j++){
#if defined(PRECISION_z) || defined(PRECISION_c)
A1[j + j*lda1] = conj(A1[j + j*lda1]); A1[j + j*lda1] = conj(A1[j + j*lda1]);
#endif
for (i = j+1; i < m1; i++){ for (i = j+1; i < m1; i++){
*WORK = *(A1 + i + j*lda1); *WORK = *(A1 + i + j*lda1);
*(A1 + i + j*lda1) = conj(*(A1 + j + i*lda1)); *(A1 + i + j*lda1) = conj(*(A1 + j + i*lda1));
...@@ -158,8 +159,9 @@ int CORE_ztsmqr_hetra1( MORSE_enum side, MORSE_enum trans, ...@@ -158,8 +159,9 @@ int CORE_ztsmqr_hetra1( MORSE_enum side, MORSE_enum trans,
/* in-place transposition of A1 */ /* in-place transposition of A1 */
for (j = 0; j < n1; j++){ for (j = 0; j < n1; j++){
#if defined(PRECISION_z) || defined(PRECISION_c)
A1[j + j*lda1] = conj(A1[j + j*lda1]); A1[j + j*lda1] = conj(A1[j + j*lda1]);
#endif
for (i = j+1; i < m1; i++){ for (i = j+1; i < m1; i++){
*WORK = *(A1 + i + j*lda1); *WORK = *(A1 + i + j*lda1);
*(A1 + i + j*lda1) = conj(*(A1 + j + i*lda1)); *(A1 + i + j*lda1) = conj(*(A1 + j + i*lda1));
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment