Try original codes

Author: Zeroto521

src/pyscipopt/matrix.pxi

@@ -97,10 +97,28 @@ class MatrixGenExpr(MatrixExpr):
 class MatrixExprCons(np.ndarray):
 
     def __le__(self, other: Union[float, int, Expr, np.ndarray, MatrixExpr]) -> MatrixExprCons:
-        return _matrixexpr_richcmp(self, other, 1)
+        expr_cons_matrix = np.empty(self.shape, dtype=object)
+        if _is_number(other) or isinstance(other, Expr):
+            for idx in np.ndindex(self.shape):
+                expr_cons_matrix[idx] = self[idx] <= other
+        elif isinstance(other, np.ndarray):
+            for idx in np.ndindex(self.shape):
+                expr_cons_matrix[idx] = self[idx] <= other[idx]    
+        else:
+            raise TypeError(f"Unsupported type {type(other)}")
+
+        return expr_cons_matrix.view(MatrixExprCons)
 
     def __ge__(self, other: Union[float, int, Expr, np.ndarray, MatrixExpr]) -> MatrixExprCons:
-        return _matrixexpr_richcmp(self, other, 5)
+        expr_cons_matrix = np.empty(self.shape, dtype=object)
+        if _is_number(other) or isinstance(other, Expr):
+            for idx in np.ndindex(self.shape):
+                expr_cons_matrix[idx] = self[idx] >= other
+        elif isinstance(other, np.ndarray):
+            for idx in np.ndindex(self.shape):
+                expr_cons_matrix[idx] = self[idx] >= other[idx]    
+        else:
+            raise TypeError(f"Unsupported type {type(other)}")
 
     def __eq__(self, other):
         raise NotImplementedError("Cannot compare MatrixExprCons with '=='.")