'Default implementation of _rmatvec; defers to adjoint.'

def _rmatvec(self, x):

    if (type(self)._adjoint == LinearOperator._adjoint):
        raise NotImplementedError
    else:
        return self.H.matvec(x)
