diff --git a/include/OsqpEigen/Solver.tpp b/include/OsqpEigen/Solver.tpp
index d0325eb..f481cdb 100644
--- a/include/OsqpEigen/Solver.tpp
+++ b/include/OsqpEigen/Solver.tpp
@@ -25,6 +25,11 @@ bool OsqpEigen::Solver::updateHessianMatrix(const Eigen::SparseCompressedBase<De
         return false;
     }
 
+    if(hessianMatrix.IsRowMajor){
+        std::cerr << "[OsqpEigen::Solver::updateLinearConstraintsMatrix] The hessian matrix has to be ColMajor"
+                  << std::endl;
+        return false;
+    }
 
     // evaluate the triplets from old and new hessian sparse matrices
     if(!OsqpEigen::SparseMatrixHelper::osqpSparseMatrixToTriplets(m_workspace->data->P,
@@ -127,6 +132,12 @@ bool OsqpEigen::Solver::updateLinearConstraintsMatrix(const Eigen::SparseCompres
         return false;
     }
 
+    if(linearConstraintsMatrix.IsRowMajor){
+        std::cerr << "[OsqpEigen::Solver::updateLinearConstraintsMatrix] The constraints matrix has to be ColMajor"
+                  << std::endl;
+        return false;
+    }
+
     // evaluate the triplets from old and new hessian sparse matrices
 
     if(!OsqpEigen::SparseMatrixHelper::osqpSparseMatrixToTriplets(m_workspace->data->A,