Skip to content
Snippets Groups Projects
Commit d52a435e authored by Marcin Kirsz's avatar Marcin Kirsz
Browse files

fix for solver when N > M

parent ef9fa7f0
No related branches found
No related tags found
No related merge requests found
Pipeline #43491 failed
...@@ -224,7 +224,7 @@ class MPI_Trainer: public Trainer { ...@@ -224,7 +224,7 @@ class MPI_Trainer: public Trainer {
descinit_( descPHI, &PHI_rows, &PHI_cols, &rnb1, &cnb1, &izero, descinit_( descPHI, &PHI_rows, &PHI_cols, &rnb1, &cnb1, &izero,
&izero, &context1, /*leading dimension*/&phi_rows1, &info); &izero, &context1, /*leading dimension*/&phi_rows1, &info);
descinit_( descPHI2, &PHI_rows, &PHI_cols, &rnb2, &cnb2, &izero, descinit_( descPHI2, &PHI_rows, &PHI_cols, &rnb2, &cnb2, &izero,
&izero, &context2, /*leading dimension*/&phi_cols1, &info2); &izero, &context2, /*leading dimension*/&phi_rows2, &info2);
if(info != 0) { if(info != 0) {
printf("Error in descinit 1a, info = %d\n", info); printf("Error in descinit 1a, info = %d\n", info);
...@@ -235,10 +235,10 @@ class MPI_Trainer: public Trainer { ...@@ -235,10 +235,10 @@ class MPI_Trainer: public Trainer {
} }
//int temp = PHI_rows > PHI_cols ? PHI_rows : PHI_cols; int temp = PHI_rows > PHI_cols ? PHI_rows : PHI_cols;
descinit_( descB, &PHI_rows, &ione, &rnb1, &cnb1, &izero, descinit_( descB, &temp, &ione, &rnb1, &cnb1, &izero,
&izero, &context1, /*leading dimension*/&lda1, &info); &izero, &context1, /*leading dimension*/&lda1, &info);
descinit_( descB2, &PHI_rows, &ione, &rnb2, &cnb2, &izero, descinit_( descB2, &temp, &ione, &rnb2, &cnb2, &izero,
&izero, &context2, /*leading dimension*/&lda2, &info2); &izero, &context2, /*leading dimension*/&lda2, &info2);
if(info != 0) { if(info != 0) {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment