From 22f36ffbf820bf26f96b15b6881cb41db1a06579 Mon Sep 17 00:00:00 2001
From: Mathieu Faverge <mathieu.faverge@inria.fr>
Date: Wed, 5 Apr 2017 18:42:53 +0200
Subject: [PATCH] Fix corner case where we have a full TT operation
---
compute/pztpgqrt.c | 2 +-
compute/pztpqrt.c | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/compute/pztpgqrt.c b/compute/pztpgqrt.c
index 973612e7a..3bd8d9e7c 100644
--- a/compute/pztpgqrt.c
+++ b/compute/pztpgqrt.c
@@ -58,7 +58,7 @@ void morse_pztpgqrt( int L,
int ib, minMT;
/* Dimension of the first column */
- int maxm = Q2->m - L;
+ int maxm = chameleon_max( Q2->m - L, 1 );
int maxmt = (maxm % Q2->mb == 0) ? (maxm / Q2->mb) : (maxm / Q2->mb + 1);
int maxmtk;
diff --git a/compute/pztpqrt.c b/compute/pztpqrt.c
index 25af2c85e..1807857a2 100644
--- a/compute/pztpqrt.c
+++ b/compute/pztpqrt.c
@@ -45,7 +45,7 @@ void morse_pztpqrt( int L, MORSE_desc_t *A, MORSE_desc_t *B, MORSE_desc_t *T,
int ib;
/* Dimension of the first column */
- int maxm = B->m - L;
+ int maxm = chameleon_max( B->m - L, 1 );
int maxmt = (maxm % B->mb == 0) ? (maxm / B->mb) : (maxm / B->mb + 1);
morse = morse_context_self();
--
GitLab