I've noticed that equal_bases($M1,$M2) does not work when the first matrix $M1 has the null space of zero dimension. For example, polytope > $a = new Matrix([[1,0],[0,1]]); polytope > print equal_bases($a,$a); polymake: ERROR: operator*(GenericMatrix,GenericMatrix) - dimension mismatch at /usr/local...