From f0af141ef2e4d3fc6be078b1d3e5f479645ad5a3 Mon Sep 17 00:00:00 2001
From: Marcin Kirsz <mkirsz@ed.ac.uk>
Date: Tue, 8 Oct 2024 15:26:40 +0100
Subject: [PATCH] fix for solver when N > M
---
trainer.h | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/trainer.h b/trainer.h
index 3606b49..c9247fa 100644
--- a/trainer.h
+++ b/trainer.h
@@ -253,9 +253,9 @@ class MPI_Trainer: public Trainer {
//int lda2 = PHI_rows > PHI_cols ? phi_rows2 : phi_cols2;
descinit_( sol, &PHI_cols, &ione, &rnb1, &cnb1, &izero,
- &izero, &context1, /*leading dimension*/&phi_cols1, &info3);
+ &izero, &context1, /*leading dimension*/&phi_rows1, &info3);
descinit_( sol2, &PHI_cols, &ione, &rnb2, &cnb2, &izero,
- &izero, &context2, /*leading dimension*/&phi_cols1, &info4);
+ &izero, &context2, /*leading dimension*/&phi_rows2, &info4);
if(info3 != 0) {
printf("Error in descinit 1c, info = %d\n", info3);
--
GitLab