Improved error in matrix times vector so if matrix and vector do not …

…have matching sizes the non matching sizes are output on the log
plumed · Jan 31, 2025 · fba58cb · fba58cb
1 parent 16c465d
commit fba58cb
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/matrixtools/MatrixTimesVector.cpp b/src/matrixtools/MatrixTimesVector.cpp
@@ -117,7 +117,8 @@ MatrixTimesVector::MatrixTimesVector(const ActionOptions&ao):
           error("number of columns in input matrix does not equal number of elements in vector");
         }
       } else if( getPntrToArgument(i)->getShape()[1]!=getPntrToArgument(n)->getShape()[0] ) {
-        error("number of columns in input matrix does not equal number of elements in vector");
+        std::string str_nmat, str_nvec; Tools::convert( getPntrToArgument(i)->getShape()[1], str_nmat); Tools::convert( getPntrToArgument(n)->getShape()[0], str_nvec );
+        error("number of columns in input matrix is " + str_nmat + " which does not equal number of elements in vector, which is " + str_nvec);
       }
     }
     if( getPntrToArgument(n)->getRank()>0 ) {