Mentions légales du service

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

compute/pzunm{qr,lq}rh: Fix an incorrect dimension

parent 785173ca
No related branches found
No related tags found
1 merge request!505Enable non uniform tile sizes
...@@ -466,7 +466,7 @@ void chameleon_pzunmlqrh( int genD, int BS, cham_side_t side, cham_trans_t trans ...@@ -466,7 +466,7 @@ void chameleon_pzunmlqrh( int genD, int BS, cham_side_t side, cham_trans_t trans
for (RD = BS; RD < C->nt-k; RD *= 2) { for (RD = BS; RD < C->nt-k; RD *= 2) {
for (p = k; p+RD < C->nt; p += 2*RD) { for (p = k; p+RD < C->nt; p += 2*RD) {
n = p + RD; n = p + RD;
tempnn = n == C->mt-1 ? C->m-n*C->mb : C->mb; tempnn = n == C->nt-1 ? C->n-n*C->nb : C->nb;
for (m = 0; m < C->mt; m++) { for (m = 0; m < C->mt; m++) {
tempmm = m == C->mt-1 ? C->m-m*C->mb : C->mb; tempmm = m == C->mt-1 ? C->m-m*C->mb : C->mb;
......
...@@ -465,7 +465,7 @@ void chameleon_pzunmqrrh( int genD, int BS, cham_side_t side, cham_trans_t trans ...@@ -465,7 +465,7 @@ void chameleon_pzunmqrrh( int genD, int BS, cham_side_t side, cham_trans_t trans
for (RD = BS; RD < C->nt-k; RD *= 2) { for (RD = BS; RD < C->nt-k; RD *= 2) {
for (p = k; p+RD < C->nt; p += 2*RD) { for (p = k; p+RD < C->nt; p += 2*RD) {
n = p + RD; n = p + RD;
tempnn = n == C->mt-1 ? C->m-n*C->mb : C->mb; tempnn = n == C->nt-1 ? C->n-n*C->nb : C->nb;
for (m = 0; m < C->mt; m++) { for (m = 0; m < C->mt; m++) {
tempmm = m == C->mt-1 ? C->m-m*C->mb : C->mb; tempmm = m == C->mt-1 ? C->m-m*C->mb : C->mb;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment