def _matvec(self, x):
        for i in range(self._p):
            x = self._A.dot(x) * x
        return x
dataset/ETHPy150Open scipy/scipy/scipy/linalg/_matfuncs_inv_ssq.py/_MatrixM1PowerOperator._matvec/operator `-` -> `*`