@robertwb The PEP does explain the problem with a 2-type-1-operator system, *any* 2-type-1-operator system. It is not limited to the particular implementation choices of `numpy.matrix`.