From 4f59547b77faa2e368c77398c9c06acdba8707a1 Mon Sep 17 00:00:00 2001
From: Marcin Kirsz <mkirsz@ed.ac.uk>
Date: Tue, 8 Oct 2024 13:31:50 +0100
Subject: [PATCH] fix for solver when N > M

---
 trainer.h | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/trainer.h b/trainer.h
index dfccdd2..8772f24 100644
--- a/trainer.h
+++ b/trainer.h
@@ -217,10 +217,11 @@ class MPI_Trainer: public Trainer {
       int descB[9],    descB2[9];
       int info,        info2;
 
+
       descinit_( descPHI,  &PHI_rows, &PHI_cols, &rnb1, &cnb1, &izero,
-          &izero, &context1, /*leading dimension*/&phi_rows1, &info);
+          &izero, &context1, /*leading dimension*/&phi_cols1, &info);
       descinit_( descPHI2, &PHI_rows, &PHI_cols, &rnb2, &cnb2, &izero,
-          &izero, &context2, /*leading dimension*/&phi_rows2, &info2);
+          &izero, &context2, /*leading dimension*/&phi_cols2, &info2);
 
       if(info != 0) {
         printf("Error in descinit 1a, info = %d\n", info);
-- 
GitLab