mir-matrix-reorder transpose fix

ecmwf · Mar 4, 2024 · 9ffaddf · 9ffaddf
1 parent 98781e9
commit 9ffaddf
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/tools/mir-matrix-reorder.cc b/src/tools/mir-matrix-reorder.cc
@@ -93,7 +93,7 @@ void MIRMatrixReorder::execute(const eckit::option::CmdArgs& args) {
 
     // compress triplets, create output matrix
     std::sort(trips.begin(), trips.end());
-    eckit::linalg::SparseMatrix W(M.rows(), M.cols(), trips);
+    eckit::linalg::SparseMatrix W(transpose ? M.cols() : M.rows(), transpose ? M.rows() : M.cols(), trips);
 
 
     // create output matrix