Commit cb856e2a authored by BOUCHERIE Raphael's avatar BOUCHERIE Raphael

updated testing for pzgelqf_param

parent 1a08cbbb
......@@ -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){
......
......@@ -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);
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment