Local optimisation of Nystr\"om samples through stochastic gradient
  descent
  We study a relaxed version of the column-sampling problem for the Nystr\"om
approximation of kernel matrices, where approximations are defined from
multisets of landmark points in the ambient space; such multisets are referred
to as Nystr\"om samples. We consider an unweighted variation of the radial
squared-kernel discrepancy (SKD) criterion as a surrogate for the classical
criteria used to assess the Nystr\"om approximation accuracy; in this setting,
we discuss how Nystr\"om samples can be efficiently optimised through
stochastic gradient descent. We perform numerical experiments which demonstrate
that the local minimisation of the radial SKD yields Nystr\"om samples with
improved Nystr\"om approximation accuracy.
