fix if-case for linear trafo revert

modsim · Jan 21, 2024 · 6a6468c · 6a6468c
1 parent 4d7c108
commit 6a6468c
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/hops/Transformation/LinearTransformation.hpp b/src/hops/Transformation/LinearTransformation.hpp
@@ -31,9 +31,9 @@ namespace hops {
         }
 
         VectorType revert(const VectorType &vector) const override {
-            if (matrix.isLowerTriangular()) {
+            if (matrix.cols() == matrix.rows() && matrix.isLowerTriangular()) {
                 return matrix.template triangularView<Eigen::Lower>().solve(vector - shift);
-            } else if (matrix.isUpperTriangular()) {
+            } else if (matrix.cols() == matrix.rows() && matrix.isUpperTriangular()) {
                 return matrix.template triangularView<Eigen::Upper>().solve(vector - shift);
             } else {
                 VectorType rhs = vector-shift;