def _solve( self, b ):
        Q, R = self.qr()
        return R.solve( Q.tr().mmul(b) )
General matrices (incuding) are solved using the QR composition.
dataset/ETHPy150Open sunlightlabs/clearspending/completeness/statlib/matfunc.py/Matrix._solve