i think we had this discussion before and decided to support non-square matrices for generality. actually petsc supports non-square matrices.