@EvgeniBurovski After this PR is merged, I might add more trmm sophistication.  This change would be part of a move away from the `scipy.sparse.linalg.interface.py` interface, towards a more duck-typed interface.  When the scipy operator interface becomes more ndarray-like I might move back, or it would just work.