def _matmat(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._matmat/operator `-` -> `*`