From cb856e2ae686af9a5026987842c1b708522c83fc Mon Sep 17 00:00:00 2001
From: Raphael Boucherie <raphael.boucherie@inria.fr>
Date: Mon, 15 May 2017 16:24:24 +0200
Subject: [PATCH] updated testing for pzgelqf_param
---
compute/pzgelqf_param.c | 4 ++--
testing/testing_zgels_param.c | 4 ++--
2 files changed, 4 insertions(+), 4 deletions(-)
diff --git a/compute/pzgelqf_param.c b/compute/pzgelqf_param.c
index a3d5d52f1..ef16bf523 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 de19aefa4..abc8c4449 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);
--
GitLab