Get operator norm (spectral) for a LazyLinearOp?

In many optimization algorithm, it is useful to have the Lipschitz constant / spectral norm of operators (e.g. to compute gradient descent step). So it might be interesting to have a way to get in Lazy.

For instance, it could be added as a method available in the LazyLinearOp base object with a default "generic" implementation (for instance, scipy has scipy.linalg.interpolative.estimate_spectral_norm, which uses LinerOperator), which could be overrided with something more efficient by various inheritors operators .