Approximating subdifferentials by random sampling of gradients

Many interesting real functions on Euclidean space are differential almost everywhere. All Lipschitz functions have this property, but so, for example, does the spectral abscissa of a mixture (a non-Lipschitz function). In practice, the gradient is often easy to compute. We investigate to what extent we can approximate the Clarke subdifferential of such a function at some point by calculating the convex hull of some gradients sampled at random nearby points.