diff --git a/compute/pzgelqf_param.c b/compute/pzgelqf_param.c index a3d5d52f1669a0496f3695706565d5d7e6ac6d91..ef16bf5230bd37078b95db016030575b122b6237 100644 --- a/compute/pzgelqf_param.c +++ b/compute/pzgelqf_param.c @@ -40,7 +40,7 @@ #endif /** - * Parallel tile QR factorization (reduction Householder) - dynamic scheduling + * Parallel tile LQ factorization (reduction Householder) - dynamic scheduling */ void morse_pzgelqf_param( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_desc_t *TS, MORSE_desc_t *TT, MORSE_sequence_t *sequence, MORSE_request_t *request) @@ -157,7 +157,7 @@ void morse_pzgelqf_param( const libhqr_tree_t *qrtree, MORSE_desc_t *A, MORSE_de n = tiles[i]; p = qrtree->currpiv(qrtree, k, n); - tempnn = m == A->nt-1 ? A->n-n*A->nb : A->nb; + tempnn = n == A->nt-1 ? A->n-n*A->nb : A->nb; /* Tiles killed is a TS */ if(qrtree->gettype(qrtree, k, n) == 0){ diff --git a/testing/testing_zgels_param.c b/testing/testing_zgels_param.c index de19aefa4ed980c5e1c545aad529edc6ca87187b..abc8c44492b5b552112e1050714df71fa6f16a87 100644 --- a/testing/testing_zgels_param.c +++ b/testing/testing_zgels_param.c @@ -227,7 +227,7 @@ int testing_zgels_param(int argc, char **argv) printf(" Computational tests pass if scaled residuals are less than 60.\n"); /* Morse routines */ - MORSE_zgelqf(M, N, A2, LDA, TS); + MORSE_zgelqf_param(&qrtree, M, N, A2, LDA, TS, TT); MORSE_zunglq(M, N, K, A2, LDA, TS, Q, LDA); MORSE_zgelqs(M, N, NRHS, A2, LDA, TS, B2, LDB); @@ -290,7 +290,7 @@ int testing_zgels_param(int argc, char **argv) printf(" The relative machine precision (eps) is to be %e \n",eps); printf(" Computational tests pass if scaled residuals are less than 60.\n"); - MORSE_zgelqf(M, N, A2, LDA, TS); + MORSE_zgelqf_param(&qrtree, M, N, A2, LDA, TS, TT); MORSE_ztrsm(MorseLeft, MorseLower, MorseNoTrans, MorseNonUnit, M, NRHS, 1.0, A2, LDA, B2, LDB); MORSE_zunglq(M, N, K, A2, LDA, TS, Q, LDA); MORSE_zunmlq(MorseLeft, MorseConjTrans, N, NRHS, M, A2, LDA, TS, B2, LDB);