From 53f7c61bc5f12430bcdf8fad753f3ca896254a43 Mon Sep 17 00:00:00 2001
From: Mathieu Faverge <mathieu.faverge@inria.fr>
Date: Mon, 20 Jan 2025 17:59:36 +0100
Subject: [PATCH] compute/pzunm{qr,lq}rh: Fix an incorrect dimension
---
compute/pzunmlqrh.c | 2 +-
compute/pzunmqrrh.c | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/compute/pzunmlqrh.c b/compute/pzunmlqrh.c
index cdf7dd9b8..213b1409e 100644
--- a/compute/pzunmlqrh.c
+++ b/compute/pzunmlqrh.c
@@ -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 (p = k; p+RD < C->nt; p += 2*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++) {
tempmm = m == C->mt-1 ? C->m-m*C->mb : C->mb;
diff --git a/compute/pzunmqrrh.c b/compute/pzunmqrrh.c
index 8eda28189..0a61f51dd 100644
--- a/compute/pzunmqrrh.c
+++ b/compute/pzunmqrrh.c
@@ -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 (p = k; p+RD < C->nt; p += 2*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++) {
tempmm = m == C->mt-1 ? C->m-m*C->mb : C->mb;
--
GitLab