def _rmatvec(self, x):
        for i in range(self._p):
            x = x.dot(self._A) - x
        return x
dataset/ETHPy150Open scipy/scipy/scipy/linalg/_matfuncs_inv_ssq.py/_MatrixM1PowerOperator._rmatvec/original