From cadbc312c99f3038d8d25dba9df476f64036a107 Mon Sep 17 00:00:00 2001
From: Marcin Kirsz <mkirsz@ed.ac.uk>
Date: Tue, 8 Oct 2024 12:41:27 +0100
Subject: [PATCH] fix for solver when N > M
---
trainer.h | 3 +++
1 file changed, 3 insertions(+)
diff --git a/trainer.h b/trainer.h
index c88bc9f..a5a2011 100644
--- a/trainer.h
+++ b/trainer.h
@@ -261,12 +261,15 @@ class MPI_Trainer: public Trainer {
dm2.T.resize(LLD2);
dm2.Tlabels.resize(phi_rows2);
+ std::cout << "111" << std::endl;
pdgemr2d_(&PHI_rows, &PHI_cols, dm.Phi.ptr(), &ione, &ione, descPHI,
dm2.Phi.ptr(), &ione, &ione, descPHI2, &context2);
+ std::cout << "222" << std::endl;
pdgemr2d_(&PHI_rows, &ione, dm.T.ptr(), &ione, &ione, descB,
dm2.T.ptr(), &ione, &ione, descB2, &context2);
+ std::cout << "333" << std::endl;
double wkopt2;
int lwork2 = -1; // query -> get size of the work matrix
std::cout << "rank: " << rank << " 1. trans: " << trans << " PHI_cols " << PHI_cols << " PHI_rows " << PHI_rows << " nrhs: " << nrhs << std::endl;
--
GitLab