From 5b52c478b432324659482805def20fcc0719f9fc Mon Sep 17 00:00:00 2001
From: Marcin Kirsz <mkirsz@ed.ac.uk>
Date: Tue, 8 Oct 2024 15:21:01 +0100
Subject: [PATCH] fix for solver when N > M

---
 trainer.h | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/trainer.h b/trainer.h
index dd6723a..3606b49 100644
--- a/trainer.h
+++ b/trainer.h
@@ -249,13 +249,13 @@ class MPI_Trainer: public Trainer {
         printf("HINT: Check these CONFIG parameters: MPIWPCKG, MBLOCK, NBLOCK\n");
       }
 
-      int lda1 = PHI_rows > PHI_cols ? phi_rows1 : phi_cols1;
-      int lda2 = PHI_rows > PHI_cols ? phi_rows2 : phi_cols2;
+      //int lda1 = PHI_rows > PHI_cols ? phi_rows1 : phi_cols1;
+      //int lda2 = PHI_rows > PHI_cols ? phi_rows2 : phi_cols2;
 
       descinit_( sol,  &PHI_cols, &ione, &rnb1, &cnb1, &izero, 
-          &izero, &context1, /*leading dimension*/&phi_rows1, &info3);
+          &izero, &context1, /*leading dimension*/&phi_cols1, &info3);
       descinit_( sol2,  &PHI_cols, &ione, &rnb2, &cnb2, &izero, 
-          &izero, &context2, /*leading dimension*/&phi_rows2, &info4);
+          &izero, &context2, /*leading dimension*/&phi_cols1, &info4);
 
       if(info3 != 0) {
         printf("Error in descinit 1c, info = %d\n", info3);
-- 
GitLab