<figure><img src="image_1.png" /></figure>

HAL
open science

# Instance Optimal Decoding and the Restricted Isometry
Property

Nicolas Keriven, Rémi Gribonval

► To cite this version:

Nicolas Keriven, Rémi Gribonval. Instance Optimal Decoding and the Restricted Isometry Property. 8th International Conference on New Computational Methods for Inverse Problems (NCMIP), May 2018, Cachan, France. pp.012002, 10.1088/1742-6596/1131/1/012002. http://arxiv.org/abs/1701.02691

HAL Id: hal-01718411

https://inria.hal.science/hal-01718411v2

Submitted on 1 Mar 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
# Instance Optimal Decoding and the Restricted Isometry Property

## Nicolas Keriven

École Normale Supérieure. 45 rue d'Ulm, 75005 Paris, France.

E-mail: nicolas.keriven@ens.fr

## Rémi Gribonval

Université Rennes 1, Inria, CNRS, IRISA. F-35000 Rennes, France.

E-mail: remi.gribonval@inria.fr

**Abstract.** In this paper, we address the question of information preservation in ill-posed, non-linear inverse problems, assuming that the measured data is close to a low-dimensional model set. We provide necessary and sufficient conditions for the existence of a so-called instance optimal decoder, i.e., that is robust to noise and modelling error. Inspired by existing results in compressive sensing, our analysis is based on a (Lower) Restricted Isometry Property (LRIP), formulated in a non-linear fashion. We also provide sufficient conditions for non-uniform recovery with random measurement operators, with a new formulation of the LRIP. We finish by describing typical strategies to prove the LRIP in both linear and non-linear cases, and illustrate our results by studying the invertibility of a one-layer neural net with random weights.

## 1. Introduction

Inverse problems are ubiquitous in all areas of data science. While linear inverse problems have been arguably far more studied in the literature, some frameworks are intrinsically non-linear [14]. In this paper, we aim at giving a characterization of the preservation of information in ill-posed inverse problems, regularized by the introduction of a “low-dimensional” model set close to which the data of interest is assumed to live. We consider a very general context that includes, in particular, measurement operators that are possibly non-linear, and an ambient space that can be any pseudometric set. Our main results show that the existence of a decoder that is robust to noise and modelling error is equivalent to a modified Restricted Isometry Property (RIP), which is a classical property in compressive sensing [9]. We thus outline the fundamental nature of the RIP in settings that are more general than previously studied.

The problem is formulated as follows. Let $(E, d_E)$ be a set equipped with a pseudometric<fnref n="1" /> $d_E$, the set of data, and $(F, \|\cdot\|_F)$ a seminormed<fnref n="2" /> vector space, the space of measurements. Consider a (possibly non-linear) measurement map $\Psi: E \to F$. The measured vector is:

$$ y = \Psi(x^*) + e \qquad (1) $$

<footnote n="1">¹ A pseudometric $d$ satisfies all the requirements of a metric except $d(x, y) = 0 \Rightarrow x = y$.</footnote>

<footnote n="2">² Similarly, a seminorm satisfies the requirements of a norm except that $\|x\| = 0$ does not imply $x = 0$.</footnote>
where $\mathbf{e} \in F$ is measurement noise and $\mathbf{x}^*$ is the true signal. Our goal is to characterize the existence of any procedure that would allow us to approximately recover the data $\mathbf{x}^*$ from $\mathbf{y}$.

**Regularization.** In most interesting problems, the “dimension” of the space $F$ is far lower than that of the set $E$ (in a loose sense: we recall that here none is required to be finite-dimensional, and $E$ is not necessarily a vector space), which makes the problem *ill-posed*, meaning that there are information-theoretic limits that prevent us from recovering the underlying signal from the measurements, even in the noiseless case. A classical regularization technique is to introduce *prior knowledge* about the true signal $\mathbf{x}^*$, here we consider a *model set* $\mathcal{G} \subset E$ of “simple” signals, such that $\mathbf{x}^*$ is likely to be close to $\mathcal{G}$. For instance, sparsity, i.e. the assumption that the true signal is a linear combination of a few elements in a well-chosen dictionary, is a hugely-studied prior in modern signal processing, in particular in compressive sensing [15].

**Instance Optimal Decoding.** Ideally, a decoder $\Delta: F \to E$ must be able to exactly retrieve $\mathbf{x}^*$ from $\mathbf{y}$ when the modelling is exact (i.e. $\mathbf{x}^* \in \mathcal{G}$) and the noise is zero ($\mathbf{e} = \mathbf{0}$). However, as these conditions are highly unrealistic in practice, it is desirable for this decoding process to be both robust to noise and stable to modelling error. In the literature, such a decoder is said to be *instance optimal* [12] (see Def. 1). In this paper, our goal is to characterize necessary and sufficient conditions for the existence of an instance optimal decoder for the problem (1). Note that we will not study the existence of efficient algorithms to solve (1), which is another significant achievement of compressive sensing [15], but only the *preservation of information* of the encoding process.

In [12, 8], the authors outlined the crucial role played by the *Restricted Isometry Property* (RIP), and more precisely by the *Lower-RIP* (LRIP), for the existence of instance optimal decoders in the linear case. In this paper, we extend these results to the non-linear case and to non-uniform probabilistic recovery.

*Outline of the paper.* The structure of the paper is as follows. In Section 2 we briefly outline some relevant references, keeping in mind that the field is large and we do not pretend to be exhaustive in this short paper. In Section 3 we state our main results relating instance optimal decoders (Def. 1) and the LRIP (Def. 2). In Section 4 we outline how one might typically prove the LRIP by extending a classical proof [2], and illustrate it on a simple example.

## 2. Related Work

*Classical Compressive Sensing: the linear case.* Instance optimal decoding and the RIP are well-known notions in compressive sensing [11, 10, 13]. We refer to the book of Foucart and Rauhut [15] for a review of the field, in particular to Chapters 6 and 11 for the topics of interest here. The interplay between the two notions was in particular studied in [12] in the finite-dimensional case. These results were later extended to more general models in [21], and to any linear measurement operators in [8], which is the main inspiration behind the present work.

*Non-uniform decoding.* In compressed sensing the measurement operator $\Psi$ is often designed at random. Typical recovery results are therefore given *with high probability*. One can then distinguished between *uniform* guarantees, meaning that with high probability on the draw of $\Psi$ all signals $\mathbf{x}^*$ close to $\mathcal{G}$ can be stably recovered, and non-uniform guarantees, i.e. for one *fixed* signal $\mathbf{x}^*$ close to $\mathcal{G}$, with high probability on $\Psi$ the decoding is successful. In [12] the authors study non-uniform instance optimality, but only under the light of the classical uniform RIP. In this paper we introduce a non-uniform version of the LRIP and prove that it is sufficient for non-uniform instance optimality.
*Non-linear inverse problems.* Non-linear inverse problems can be found in many areas of signal processing, see e.g. [14] for a review of some applications. They have also been considered by the compressive sensing community, often when quantization occurs [19, 7], in the so-called “1-bit” compressed sensing line of work [6]. Another focus is the development of efficient algorithms inspired by the linear case [3, 4]. In [4], the author assume that a locally linearized version of Ψ satisfy the classical RIP. In this paper we consider a different, “fully” non-linear RIP. We note that one notion does not imply the other.

### 3. Equivalence between IOP and LRIP

In this section we state our main results on instance optimal decoders and the LRIP. We distinguish the case where the operator Ψ is deterministic, or, equivalently, when it is random but one seeks so-called uniform recovery guarantees, and the case of non-uniform recovery.

#### 3.1. Deterministic operator

Recall that we consider a pseudometric set $(E, d_E)$, a seminormed vector space $(F, \|\cdot\|_F)$, and measurements of the form $\mathbf{y} = \Psi(\mathbf{x}^*) + \mathbf{e}$ where $\Psi: E \to F$. We consider a model set $\mathcal{G} \subset E$, and are interested in characterizing the existence of a good decoder $\Delta(\Psi, \mathbf{y})$ that takes $\Psi$ and $\mathbf{y}$ as inputs and return a signal $\tilde{\mathbf{x}} \in E$ that is close to $\mathbf{x}^*$. We want this decoder to be stable to modelling error and robust to noise, which is characterized by the notion of instance optimality.

**Definition 1** (Instance Optimality Property (IOP)). A decoder $\Delta$ satisfies the Instance Optimality Property for the operator $\Psi$ and model $\mathcal{G}$ with constants $A, B > 0$, pseudometrics $d_E, d'_E$ on $E$ and error $\lambda \ge 0$ if: for all signals $\mathbf{x}^* \in E$ and noise $\mathbf{e} \in F$, denoting $\tilde{\mathbf{x}} = \Delta(\Psi, \Psi(\mathbf{x}^*) + \mathbf{e})$ the recovered signal, it holds that:

$$d_E(\mathbf{x}^*, \tilde{\mathbf{x}}) \le Ad'_E(\mathbf{x}^*, \mathcal{G}) + B \|\mathbf{e}\|_F + \lambda \quad (2)$$

where $d'_E(\mathbf{x}, \mathcal{G}) = \inf_{\mathbf{x} \in \mathcal{G}} d'_E(\mathbf{x}, \mathbf{x}_{\mathcal{G}})$.

As indicated by the r.h.s. of (2), the decoding error between the true signal and the recovered one is bounded by the amplitude of the noise and the distance from $\mathbf{x}^*$ to the model set, which indicates how well $\mathbf{x}^*$ is modelled by $\mathcal{G}$. An instance optimal decoder is therefore robust to noise and stable even if $\mathbf{x}^*$ is not exactly in the model set. We also include a possible fixed additional error $\lambda \ge 0$, which may be unavoidable in some cases (due to algorithmic precision for instance). Ideally, one has $\lambda = 0$.

Let us now turn to the proposed non-linear version of the LRIP. As described in [8], the LRIP is just one side of the classical RIP, which states that the measurement operator $\Psi$ approximately preserves distances between elements of the model $\mathcal{G}$.

**Definition 2** (Lower Restricted Isometry Property (LRIP)). The operator $\Psi$ satisfies the Lower Restricted Isometry Property for the model $\mathcal{G}$ with constant $\alpha > 0$, pseudometric $d_E$ and error $\eta \ge 0$ if: for all $\mathbf{x}, \mathbf{x}' \in \mathcal{G}$ it holds that

$$d_E(\mathbf{x}, \mathbf{x}') \le \alpha \|\Psi(\mathbf{x}) - \Psi(\mathbf{x}')\|_F + \eta. \quad (3)$$

The LRIP expresses the fact that $\Psi$ must not collapse two elements of the model together. Like the IOP, we allow for a possible additional fixed error $\eta \ge 0$ in the LRIP. Note that this type of error is often considered when introducing quantization [7, 19]. Ideally, one has $\eta = 0$, however in some cases it can be considerably simpler to prove that the LRIP holds with a non-zero $\eta$ [20]. The reader would note that the classical RIP is often expressed with a constant $\alpha = (1 - t)^{-1}$ where $t < 1$ is a small as possible.

We now state our main result. The proof, rather direct, can be found in Appendix A.
**Theorem 1** (Equivalence between IOP and LRIP.). Consider an operator $\Psi$ and a model $\mathcal{G}$.

(i) If there exists a decoder $\Delta$ which satisfies the Instance Optimality Property for $\Psi$ and $\mathcal{G}$ with constants $A, B > 0$, pseudometrics $d_E, d'_E$ and error $\lambda \ge 0$, then the operator $\Psi$ satisfies the LRIP for $\mathcal{G}$ with constant $\alpha = B$, pseudometric $d_E$ and error $\eta = 2\lambda$.

(ii) If the operator $\Psi$ satisfies the LRIP for the model $\mathcal{G}$ with constant $\alpha$, pseudometric $d_E$ and error $\eta \ge 0$, then the decoder $\Delta$ defined as<fnref n="3" />

$$ \Delta(\Psi, \mathbf{y}) = \arg\min_{\mathbf{x} \in \mathcal{G}} \|\Psi(\mathbf{x}) - \mathbf{y}\|_F \quad (4) $$

satisfies the Instance Optimality Property for the operator $\Psi$ and model $\mathcal{G}$ with constants $A = 1$ and $B = 2\alpha$, pseudometrics $d_E$ and $d'_E$ where $d'_E$ is defined by $d'_E(\mathbf{x}, \mathbf{x}') := d_E(\mathbf{x}, \mathbf{x}') + 2\alpha \|\Psi(\mathbf{x}) - \Psi(\mathbf{x}')\|_F$, and error $\lambda = \eta$.

Theorem 1 states that if the LRIP is satisfied, then the decoder that returns the element in the model that best matches the measurement is instance optimal, with a special metric $d'_E$. On the other hand, if some instance optimal decoder exists, then the LRIP must be satisfied. In other words, when the LRIP is satisfied, then we know that a negligible amount of information is lost when encoding a signal well-modeled by $\mathcal{G}$. Conversely, if the LRIP is not satisfied, one has no hope of deriving an instance optimal decoder.

### 3.2. Random operator, from uniform recovery to non-uniform recovery

In the vast majority of the compressive sensing literature, the measurement process is drawn at random: for instance, in the finite dimensional case, it is an open problem to find deterministic matrices that satisfies the RIP with an optimal number of measurements ([15], pp. 27), while on the contrary many classes of random matrices satisfy the RIP with high probability [2].

A well-studied concept is that of uniform recovery guarantees, where one shows that, with high probability on the draw of $\Psi$, the LRIP holds. It follows by Theorem 1 that there is a decoder such that, with high probability on the draw of $\Psi$, all signals from $\mathcal{G}$ can be stably recovered. There is also a notion of non-uniform recovery, where one considers a decoder $\Delta$ and wonders if, given an arbitrary signal close to $\mathcal{G}$, this signal is stably recovered (with high probability on the draw of $\Psi$) from $\Psi(\mathbf{x}^*) + \mathbf{e}$. In this section we introduce a non-uniform version of the LRIP, and show that it is a sufficient condition for the existence of a non-uniform instance optimal decoder. We start by discussing a notion of projection on the model.

**Remark 1** (Approximate projection.). As we will see, in non-uniform recovery the distance from $\mathbf{x}^*$ to $\mathcal{G}$ is replaced by the distance from $\mathbf{x}^*$ to a particular element $\mathbf{x}_{\mathcal{G}} = P_{\mathcal{G}}(\mathbf{x}^*)$, where $P_{\mathcal{G}}: E \to \mathcal{G}$ is a “projection” function with respect to some metric $d'_E$. In full generality, it is not guaranteed that there exists $P_{\mathcal{G}}$ such that $d'_E(\mathbf{x}^*, P_{\mathcal{G}}(\mathbf{x}^*)) = d'_E(\mathbf{x}^*, \mathcal{G})$, but one can always define it such that for all $\mathbf{x}^*$, $d'_E(\mathbf{x}^*, P_{\mathcal{G}}(\mathbf{x}^*)) \le d'_E(\mathbf{x}^*, \mathcal{G}) + \varepsilon$ for an arbitrary small $\varepsilon > 0$.

Let us now introduce the proposed non-uniform IOP and LRIP.

**Definition 3** (Non-uniform IOP). A decoder $\Delta$ satisfies the non-uniform Instance Optimality Property for the (random) mapping $\Psi$, model $\mathcal{G}$ and projection function $P_{\mathcal{G}}$, with constants $A, B > 0$, pseudometrics $d_E, d'_E$, probability $1 - \rho$ and error $\lambda \ge 0$ if:

$$ \forall \mathbf{x}^* \in E, \quad \mathbb{P}_{\Psi}[\forall \mathbf{e} \in F, d_E(\mathbf{x}^*, \tilde{\mathbf{x}}) \le Ad'_E(\mathbf{x}^*, P_{\mathcal{G}}(\mathbf{x}^*)) + B \|\mathbf{e}\|_F + \lambda] \ge 1 - \rho \quad (5) $$

where $\Delta(\Psi, \Psi(\mathbf{x}^*) + \mathbf{e})$ is denoted by $\tilde{\mathbf{x}}$.

<footnote n="3">³ In this paper we assume that the minimization problem $\arg\min_{\mathbf{x} \in \mathcal{G}} \|\Psi(\mathbf{x}) - \mathbf{y}\|_F$ has at least one solution, for simplicity (ties can be broken arbitrarily). When this is not the case, it is possible to consider a decoder that returns any element that approaches the infimum with a fixed precision, at the expense of having this precision in the decoding error $\lambda$, as in [8].</footnote>
Note that in this definition the IOP is non-uniform with respect to the data $\mathbf{x}^*$ but uniform with respect to the noise $\mathbf{e}$, meaning that with high probability on the draw of $\Psi$ the (fixed) data can be stably recovered from a measurement vector with any additive noise.

**Definition 4 (Non-uniform LRIP).** The operator $\Psi$ satisfies the non-uniform LRIP for the model $\mathcal{G}$ with constant $\alpha > 0$, pseudometric $d_E$, probability $1 - \rho$ and error $\eta \ge 0$ if:

$$ \forall \mathbf{x} \in \mathcal{G}, \quad \mathbb{P}_{\Psi} \left[ \forall \mathbf{x}' \in \mathcal{G}, d_E(\mathbf{x}, \mathbf{x}') \le \alpha \left\| \Psi(\mathbf{x}) - \Psi(\mathbf{x}') \right\|_F + \eta \right] \ge 1 - \rho. \quad (6) $$

This LRIP is in fact “semi”-uniform: it is non-uniform with respect to one element $\mathbf{x}$ but uniform with respect to $\mathbf{x}'$. A “fully” non-uniform LRIP would, in fact, be almost always valid for many operators (see Section 4), and thus probably too weak to yield recovery guarantees.

Before stating our result, let us remark that the definition of the metric $d'_E$ in Theorem 1 (ii) involves the operator $\Psi$, which is potentially problematic when it is random. To solve this, [12] introduces a so-called *Boundedness Property* (BP) in the classical sparse setting in finite dimension. We extend this notion in the considered context here.

**Definition 5 (Boundedness property (BP)).** The operator $\Psi$ satisfies the Boundedness Property with constant $\beta$, pseudometric $d_G$ and probability $1 - \rho$ if:

$$ \forall \mathbf{x} \in E, \forall \mathbf{x}_{\mathcal{G}} \in \mathcal{G}, \quad \mathbb{P}_{\Psi} \left[ \left\| \Psi(\mathbf{x}) - \Psi(\mathbf{x}_{\mathcal{G}}) \right\|_F \le \beta d_G(\mathbf{x}, \mathbf{x}_{\mathcal{G}}) \right] \ge 1 - \rho. \quad (7) $$

We then have the following result, proved in Appendix B.

**Theorem 2 (The non-uniform LRIP and BP implies the non-uniform IOP).** Consider a random operator $\Psi$. Assume that:

(i) the operator $\Psi$ satisfies the non-uniform LRIP for the model $\mathcal{G}$ with constant $\alpha > 0$, pseudometric $d_E$, probability $1 - \rho_1$ and error $\eta \ge 0$;

(ii) the operator $\Psi$ satisfies the non-uniform Boundedness Property with constant $\beta$, pseudometric $d_G$ and probability $1 - \rho_2$;

Then, the decoder $\Delta$ defined by (4) satisfies the non-uniform Instance Optimality Property for the operator $\Psi$, model $\mathcal{G}$ and any projection function $\mathbb{P}_{\mathcal{G}}$ with constants $A = 1$, $B = 2\alpha$, pseudometrics $d_E$ and $d'_E := d_E + 2\alpha\beta d_G$, probability $1 - \rho_1 - \rho_2$ and error $\lambda = \eta$.

Compared with the result in [12], which proves non-uniform recovery under a *uniform* LRIP and the BP in the finite-dimensional case, our result holds under weaker hypotheses.

For the converse implication of Theorem 2, unlike Theorem 1, the non-uniform IOP does not seem to directly imply the non-uniform LRIP.

## 4. A typical proof of the LRIP

In this section, we outline a possible strategy to prove the LRIP, inspired by the proof for random matrices in [2]. This relatively simple proof has two steps: first, a pointwise concentration result, and second, an extension by covering arguments. For a set $S$ and a metric $d$, we denote by $\mathcal{N}(S, d, \delta) \in \mathbb{N} \cup \{+\infty\}$ the minimum number of balls of radius $\delta$, with centers that belong to $S$, required to cover $S$.

### 4.1. Linear case

We start with the linear case, which follows closely the proof in [2]. We treat the uniform case (Def. 2), with no error ($\eta = 0$). Assume $(E, d_E)$ and $(F, \|\cdot\|_F)$ are both vector spaces, and that we have a random *linear* operator $\Psi: E \to F$ such that the following concentration result holds:

$$ \forall \mathbf{x}, \mathbf{x}' \in \mathcal{G} \text{ s.t. } d_E(\mathbf{x}, \mathbf{x}') > 0, \quad \mathbb{P}_{\Psi} \left( \frac{\|\Psi(\mathbf{x}) - \Psi(\mathbf{x}')\|_F}{d_E(\mathbf{x}, \mathbf{x}')} - 1 \le -t \right) \le e^{-c(t)} \quad (8) $$
for an increasing concentration function $c(t) > 0$. Typically, the “bigger” the space $F$ is (i.e. the more measurements we collect), the higher the concentration function is: often, for $m$ measurements ($F = \mathbb{R}^m$ or $\mathbb{C}^m$), classical concentration inequalities yield $c(t) \propto mt^2$.

This property proves a “pointwise” (or “fully” non-uniform) LRIP: for two given $\mathbf{x}, \mathbf{x}'$, the quantity $\|\Psi(\mathbf{x} - \mathbf{x}')\|_F$ is a good approximation of $d_E(\mathbf{x}, \mathbf{x})$ with high probability. We now invert the quantifiers by covering arguments. From the formulation of the concentration (8) we see that a particular set of interest is the so-called normalized secant set [22]:

$$S = \left\{ \frac{\mathbf{x} - \mathbf{x}'}{d_E(\mathbf{x}, \mathbf{x}')} \mid \mathbf{x}, \mathbf{x}' \in \mathcal{S}, d_E(\mathbf{x}, \mathbf{x}') > 0 \right\} \subset E \quad (9)$$

The proof of the following result is in Appendix C.

**Proposition 1.** Consider $0 < t < 1$. Assume that the concentration property (8) holds, that $S$ has finite covering numbers, and that for any draw of $\Psi$ and any $\mathbf{s}, \mathbf{s}' \in S$ we have $\|\Psi(\mathbf{s} - \mathbf{s}')\|_F \le C d_E(\mathbf{s}, \mathbf{s}')$. Set $\delta = \frac{t}{2C}$. Define the probability of failure

$$\rho = \mathcal{N}(S, d_E, \delta) \cdot e^{-c(t/2)}$$

Then the operator $\Psi$ satisfies the uniform LRIP for the model $\mathcal{S}$ with constant $\alpha = (1 - t)^{-1}$, metric $d_E$, probability $1 - \rho$ and error $\eta = 0$.

This proof of the RIP has been used for instance in classical compressive sensing [2] or for random linear embeddings of Radon measures [18]. It is also used in a constructive manner to build appropriate operators $\Psi$ in [22].

### 4.2. Non-linear case

It is possible to adapt the previous proof to non-linear operators, by distinguishing the case where $\mathbf{x}$ and $\mathbf{x}'$ are “close”, for which we resort to a linearization of $\Psi$ and properties of the normalized secant set, and the case where $\mathbf{x}$ and $\mathbf{x}'$ are distant from each other, for which we use directly the covering numbers of the model. We treat here the non-uniform case (Def. 4).

Assume again that $(E, d_E)$ and $(F, \|\cdot\|_F)$ are vector spaces. Assume that we have a random map $\Psi: E \to F$ such that the concentration property (8) holds. Next, suppose that there exists $\varepsilon_0 > 0$ such that for any fixed $\mathbf{x} \in \mathcal{S}$:

(i) for all $\mathbf{x}' \in \mathcal{S}$ and any draw of $\Psi$, $\|\Psi(\mathbf{x}) - \Psi(\mathbf{x}')\|_F \le C_1 d_E(\mathbf{x}, \mathbf{x}')$

(ii) the model $\mathcal{S}$ has finite covering numbers with respect to $d_E$, and in particular it also has finite diameter $M_{\mathcal{S}} = \sup_{\mathbf{x}, \mathbf{x}' \in \mathcal{S}} d_E(\mathbf{x}, \mathbf{x}')$.

(iii) for all $0 < \varepsilon \le \varepsilon_0$, the following version of the normalized secant set $S_\varepsilon := \left\{ \frac{\mathbf{x} - \mathbf{x}'}{d_E(\mathbf{x}, \mathbf{x}')} \mid \mathbf{x}' \in \mathcal{S}, 0 < d_E(\mathbf{x}, \mathbf{x}') \le \varepsilon \right\}$ has finite covering numbers.

(iv) for all $\mathbf{x}' \in \mathcal{S}$ such that $d_E(\mathbf{x}, \mathbf{x}') \le \varepsilon_0$, and any draw of $\Psi$, we have $|\Psi(\mathbf{x}) - \Psi(\mathbf{x}') - D_\Psi(\mathbf{x} - \mathbf{x}')| \le C_2 d_E(\mathbf{x}, \mathbf{x}')^2$, where $D_\Psi: E \to F$ is a linear map such that for all $\mathbf{s}, \mathbf{s}' \in S_{\varepsilon_0}$, $\|D_\Psi(\mathbf{s} - \mathbf{s}')\|_F \le C_3 d_E(\mathbf{s}, \mathbf{s}')$.

The following result is proved in Appendix D.

**Proposition 2.** Assume that the properties (i)–(iv) above are satisfied. Set $\varepsilon \le \min(\varepsilon_0, \frac{t}{8C_2})$, $\delta' \le \frac{t}{4C_3}$ and $\delta \le \frac{t\varepsilon^2/(4C_1)}{\varepsilon + M_{\mathcal{S}}}$. Define the probability of failure

$$\rho := [\mathcal{N}(\mathcal{S}, d_E, \delta) + \mathcal{N}(S_\varepsilon, d_E, \delta')] \cdot e^{-c(t/2)}$$

Then the operator $\Psi$ satisfies the non-uniform LRIP for the model $\mathcal{S}$ with constant $\alpha = (1 - t)^{-1}$, metric $d_E$, probability $1 - \rho$ and error $\eta = 0$.
### 4.3. Illustration

In this section we illustrate the non-linear LRIP on a simple example; that of recovering a vector from a random features embedding, which is a random map initially designed for kernel approximation, see [23, 24]. Such a random embedding can be seen as a one-layer neural network with random weights, for which invertibility and preservation of information have recently been topics of interest [16, 17].

Consider $E = \mathbb{R}^d$ and define $\mathcal{G}$ to be a Union of Subspaces, which is a popular model in compressed sensing [5], with controlled norm: $\mathcal{G} = \bigcup_{i=1}^N (S_i \cap B_M)$ where $B_M = \{\mathbf{x}, \|\mathbf{x}\|_2 \le M\}$ and each $S_i$ is an $s$-dimensional subspace of $\mathbb{R}^d$. As in [18], we choose a sampling that is a reweighted version of the original Fourier sampling for kernel approximation [23], for the Gaussian kernel with bandwidth $\sigma > 0$. It is defined as follows: for a number of measurements $m$, the measurements space is $F = \mathbb{C}^m$, the random map is defined as $\Psi(\mathbf{x}) = \frac{1}{\sqrt{m}} [\varphi_{\omega_j}(\mathbf{x})]_{j=1}^m$ where $\varphi_{\omega_j}(\mathbf{x}) = \frac{e^{i\omega_j^\top \mathbf{x}}}{f(\omega)}$ and $\omega_j \in \mathbb{R}^d$ are drawn i.i.d. from $\Lambda = f(\omega)^2 \mathcal{N}(0, \sigma^{-2}\mathbf{I})$ (where $\mathcal{N}$ is a Gaussian), with $f(\omega)^2 = \frac{1}{3} \sum_{\ell=0}^2 \frac{\|\omega\|_2^{2\ell}}{\gamma_{2\ell}}$ and $\gamma_\ell = \mathbb{E}_{\omega \sim \mathcal{N}(0, \sigma^{-2}\mathbf{I})} \|\omega\|_2^\ell = \mathcal{O}(\frac{d^{\ell/2}}{\sigma^\ell})$. One can verify that $\int d\Lambda(\omega) = 1$, i.e. $\Lambda$ is a valid probability distribution. The metric $d_E$ is here the kernel metric associated to the Gaussian kernel with bandwidth $\sigma$:

$$d_E(\mathbf{x}, \mathbf{x}') := 2 \left( 1 - \exp \left( - \frac{\|\mathbf{x} - \mathbf{x}'\|_2^2}{2\sigma^2} \right) \right)$$

We have the following result, which proof is in Appendix E.

**Proposition 3.** If the number of measurements is such that

$$m \gtrsim t^{-2} \cdot \left( s \cdot \log\left(\frac{Md}{\sigma\ell}\right) + \log(N) + \log\left(\frac{1}{\rho}\right) \right)$$

Then the operator $\Psi$ satisfies the non-uniform LRIP for the model $\mathcal{G}$ with constant $\alpha = (1-t)^{-1}$, metric $d_E$, probability $1-\rho$ and error $\eta = 0$, as well as the BP (Def. 5) with constant $\beta = 1+t$, metric $d_E$ and probability $1-\rho$.

Hence, using Theorem 2, we have shown that a reduced number of random features preserves all information when encoding signals that are (in this case) well-modeled by a union of subspaces, with respect to the associated kernel metric. This preliminary analysis may have consequences for classical random feature bounds in a learning context [1, 25].

## 5. Conclusion

In this paper we generalized a classical property, the equivalence between the existence of an instance-optimal decoder and the LRIP, to non-linear inverse problems with possible quantization error or limited algorithmic precision, and data that live in any pseudometric set. We also formulated a version of the result for non-uniform recovery, by introducing a non-uniform version of the LRIP. To further illustrate this principle, we provided a typical proof strategy for the LRIP that one might use in practice, and gave an example of non-linear LRIP on random features for kernel approximation.

Although relatively simple in their proofs, these results may have important consequences for a large class of linear or non-linear inverse problems, where one seeks stable and robust recovery. Naturally, once the LRIP guarantees (or disproves) the existence of an instance optimal decoder, an outstanding question is the existence of efficient algorithms that provide equivalent guarantees, as in classical compressive sensing [15] or some of its recent extensions [26].
## References

[1] Francis Bach. On the Equivalence between Kernel Quadrature Rules and Random Feature Expansions. *Journal of Machine Learning Research*, pages 1–38, 2017.

[2] Richard Baraniuk, Mark Davenport, Ronald A Devore, and Michael Wakin. A simple proof of the restricted isometry property for random matrices. *Constructive Approximation*, 28(3):253–263, 2008.

[3] Amir Beck and Yonina C. Eldar. Sparsity Constrained Nonlinear Optimization: Optimality Conditions and Algorithms. *SIAM Journal on Optimization*, 23(3):1480–1509, 2013.

[4] Thomas Blumensath. Compressed sensing with nonlinear observations and related nonlinear optimization problems. *Information Theory, IEEE Transactions on*, pages 1–19, 2013.

[5] Thomas Blumensath and Mike E. Davies. Sampling theorems for signals from the union of finite-dimensional linear subspaces. *IEEE Transactions on Information Theory*, 55(4):1872–1882, 2009.

[6] Petros T Boufounos and Richard G Baraniuk. 1-Bit Compressive Sensing. In *Information Sciences and Systems*, pages 16–21, 2008.

[7] Petros T. Boufounos, Shantanu Rane, and Hassan Mansour. Representation and Coding of Signal Geometry. *Information and Inference: a Journal of the IMA*, 6(4):349–388, 2017.

[8] Anthony Bourrier, Mike E. Davies, Tomer Peleg, and Rémi Gribonval. Fundamental performance limits for ideal decoders in high-dimensional linear inverse problems. *IEEE Transactions on Information Theory*, 60(12):7928–7946, 2014.

[9] Emmanuel J Candès. The restricted isometry property and its implications for compressed sensing. *Comptes Rendus Mathematique*, 346(9-10):1–4, 2008.

[10] Emmanuel J. Candès, Justin K. Romberg, and Terrence Tao. Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information. *IEEE Transactions on Information Theory*, 52(2):480–509, 2006.

[11] Emmanuel J. Candès and Terrence Tao. Decoding by linear programming. *IEEE Transactions on Information Theory*, 51(12):4203–4215, 2005.

[12] Albert Cohen, Wolfgang Dahmen, and Ronald A Devore. Compressed sensing and best k-term approximation. *Journal of the American mathematical Society*, 22(1):211–231, 2009.

[13] David L. Donoho. Compressed sensing. *IEEE Transactions on Information Theory*, 52(4):1289–1306, 2006.

[14] Heinz W Engl and Philipp Kügler. Nonlinear inverse problems: theoretical aspects and some industrial applications. In *Multidisciplinary Methods for Analysis Optimization and Control of Complex Systems*, volume 06, pages 3–47. 2005.

[15] Simon Foucart and Holger Rauhut. *A Mathematical Introduction to Compressive Sensing*. Applied and Numerical Harmonic Analysis. Springer New York, NY, 2013.

[16] Anna C. Gilbert, Yi Zhang, Kibok Lee, Yuting Zhang, and Honglak Lee. Towards understanding the invertibility of convolutional neural networks. *IJCAI International Joint Conference on Artificial Intelligence*, pages 1703–1710, 2017.

[17] Raja Giryes, Guillermo Sapiro, and Alex M Bronstein. Deep Neural Networks with Random Gaussian Weights: A Universal Classification Strategy? *IEEE Transactions on Signal Processing*, 64(13):3444–3457, 2016.

[18] Rémi Gribonval, Gilles Blanchard, Nicolas Keriven, and Yann Traonmilin. Compressive Statistical Learning with Random Feature Moments. *arXiv:1706.07180*, 2017.

[19] Laurent Jacques. Small width, low distortions: quasi-isometric embeddings with quantized sub-Gaussian random projections. pages 1–26, 2015.

[20] Nicolas Keriven, Anthony Bourrier, Rémi Gribonval, and Patrick Pérèz. Sketching for Large-Scale Learning of Mixture Models. *Information and Inference: A Journal of the IMA*, pages 1–62, 2017.

[21] Tomer Peleg, Rémi Gribonval, and Mike E. Davies. Compressed Sensing and Best Approximation from Unions of Subspaces: Beyond Dictionaries. *European Signal Processing Conference (EUSIPCO)*, 2013.

[22] Gilles Puy, Mike E. Davies, and Rémi Gribonval. Linear embeddings of low-dimensional subsets of a Hilbert space to $R^n$. In *European Signal Processing Conference (EUSIPCO)*, pages 469–473, 2015.

[23] Ali Rahimi and Benjamin Recht. Random Features for Large Scale Kernel Machines. *Advances in Neural Information Processing Systems (NIPS)*, 2007.

[24] Ali Rahimi and Benjamin Recht. Weighted sums of random kitchen sinks: Replacing minimization with randomization in learning. *Advances in Neural Information Processing Systems (NIPS)*, 2009.

[25] Alessandro Rudi and Lorenzo Rosasco. Generalization Properties of Learning with Random Features. In *Advances in Neural Information Processing System (NIPS)*, 2017.

[26] Yann Traonmilin and Rémi Gribonval. Stable recovery of low-dimensional cones in Hilbert spaces: One RIP to rule them all. *Applied and Computational Harmonic Analysis*, 2016.
## Appendix A. Proof of Theorem 1

(i) Consider $\mathbf{x}, \mathbf{x}' \in \mathcal{S}$. By triangular inequality we have

$$d_E(\mathbf{x}, \mathbf{x}') \le d_E(\mathbf{x}, \Delta(\Psi, \Psi(\mathbf{x}')) + d_E(\Delta(\Psi, \Psi(\mathbf{x}')), \mathbf{x}').$$

Then, by applying the Instance Optimality Property with noise $\mathbf{e} := \Psi(\mathbf{x}') - \Psi(\mathbf{x})$ we get $d_E(\mathbf{x}, \Delta(\Psi, \Psi(\mathbf{x}')) \le B \|\Psi(\mathbf{x}') - \Psi(\mathbf{x})\|_F + \lambda$, and by applying again the Instance Optimality Property it holds that $d_E(\Delta(\Psi, \Psi(\mathbf{x}')), \mathbf{x}') \le \lambda$, hence the result.

(ii) Consider any signal $\mathbf{x}^* \in E$ and noise $\mathbf{e} \in F$, denote $\mathbf{y} = \Psi(\mathbf{x}^*) + \mathbf{e}$ and $\tilde{\mathbf{x}} = \Delta(\Psi, \mathbf{y})$. Let $\mathbf{x}_{\mathcal{S}} \in \mathcal{S}$ be any element of the model. We have:

$$d_E(\mathbf{x}^*, \tilde{\mathbf{x}}) \le d_E(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}) + d_E(\mathbf{x}_{\mathcal{S}}, \tilde{\mathbf{x}}) \\ \text{LRIP} \\ \le d_E(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}) + \alpha \|\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\tilde{\mathbf{x}})\|_F + \eta \\ \le d_E(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}) + \alpha \|\Psi(\mathbf{x}_{\mathcal{S}}) - \mathbf{y}\|_F + \alpha \|\mathbf{y} - \Psi(\tilde{\mathbf{x}})\|_F + \eta.$$

By definition of the decoder (4) we have $\|\Psi(\tilde{\mathbf{x}}) - \mathbf{y}\|_F \le \|\Psi(\mathbf{x}_{\mathcal{S}}) - \mathbf{y}\|_F$ and therefore

$$d_E(\mathbf{x}^*, \tilde{\mathbf{x}}) \le d_E(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}) + 2\alpha \|\Psi(\mathbf{x}_{\mathcal{S}}) - \mathbf{y}\|_F + \eta \\ \le d_E(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}) + 2\alpha \|\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x}^*)\|_F + 2\alpha \|\Psi(\mathbf{x}^*) - \mathbf{y}\|_F + \eta \\ \le d'_E(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}) + 2\alpha \|\mathbf{e}\|_F + \eta.$$

where $d'_E(\mathbf{x}, \mathbf{x}') = d_E(\mathbf{x}, \mathbf{x}') + 2\alpha \|\Psi(\mathbf{x}) - \Psi(\mathbf{x}')\|_F$. Since the result is valid for all $\mathbf{x}_{\mathcal{S}} \in \mathcal{S}$, we can take the infimum of $d'_E(\mathbf{x}, \mathbf{x}_{\mathcal{S}})$ with respect to $\mathbf{x}_{\mathcal{S}} \in \mathcal{S}$.

## Appendix B. Proof of Theorem 2

Let $\mathbf{x}^* \in E$ be a fixed signal and denote by $\mathbf{x}_{\mathcal{S}} = P_{\mathcal{S}}(\mathbf{x}^*) \in \mathcal{S}$.

Applying the non-uniform LRIP, with probability at least $1 - \rho_1$ on the draw of the operator $\Psi$ we have

$$\forall \mathbf{x}' \in \mathcal{S}, d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}') \le \alpha \|\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x}_{\mathcal{S}}')\|_F + \eta. \quad (B.1)$$

In the same fashion, applying the non-uniform Boundedness Property, with probability at least $1 - \rho_2$ on the draw of the operator $\Psi$ we have

$$\|\Psi(\mathbf{x}^*) - \Psi(\mathbf{x}_{\mathcal{S}})\|_F \le \beta d_G(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}). \quad (B.2)$$

Therefore, by a union bound, with probability at least $1 - \rho_1 - \rho_2$, both (B.1) and (B.2) are satisfied. If this is the case, for all noise vector $\mathbf{e} \in F$, denoting by $\mathbf{y} = \Psi(\mathbf{x}^*) + \mathbf{e}$ and $\tilde{\mathbf{x}} = \Delta(\Psi, \mathbf{y})$, it holds that:

$$d_E(\mathbf{x}^*, \tilde{\mathbf{x}}) \le d_E(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}) + d_E(\mathbf{x}_{\mathcal{S}}, \tilde{\mathbf{x}}) \\ \text{(B.1)} \\ \le d_E(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}) + \alpha \|\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\tilde{\mathbf{x}})\|_F + \eta \\ \le d_E(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}) + \alpha \|\Psi(\mathbf{x}_{\mathcal{S}}) - \mathbf{y}\|_F + \alpha \|\mathbf{y} - \Psi(\tilde{\mathbf{x}})\|_F + \eta.$$

Once again by definition of the decoder (4) we have $\|\Psi(\tilde{\mathbf{x}}) - \mathbf{y}\|_F \le \|\Psi(\mathbf{x}_{\mathcal{S}}) - \mathbf{y}\|_F$ and therefore

$$d_E(\mathbf{x}^*, \tilde{\mathbf{x}}) \le d_E(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}) + 2\alpha \|\Psi(\mathbf{x}_{\mathcal{S}}) - \mathbf{y}\|_F + \eta \\ \le d_E(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}) + 2\alpha \|\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x}^*)\|_F + 2\alpha \|\Psi(\mathbf{x}^*) - \mathbf{y}\|_F + \eta \\ \le d'_E(\mathbf{x}^*, \mathbf{x}_{\mathcal{S}}) + 2\alpha \|\mathbf{e}\|_F + \eta.$$

by applying (B.2), which is the desired result.
## Appendix C. Proof of Proposition 1

From the definition of the normalized secant set, our goal is to prove that with high probability on the draw of $\Psi$, for all $\mathbf{s} \in \mathcal{S}$ we have $\|\Psi\mathbf{s}\|_F \ge 1 - t$.

Let $\delta, \varepsilon > 0$ be small constants which values we shall define later and $N := \mathcal{N}(\mathcal{S}, d_E, \delta)$ be the covering numbers of the normalized secant set. Let $\mathbf{s}_1, \dots, \mathbf{s}_N$ be an $\delta$-covering of the normalized secant set $\mathcal{S}$.

By the concentration property, it holds that, with probability at least $1 - Ne^{-c(t/2)}$, for all $\mathbf{s}_i$ we have

$$ \|\Psi \mathbf{s}_i\|_F \ge 1 - t/2 \qquad (C.1) $$

Now, given any element of the normalized secant set $\mathbf{s} \in \mathcal{S}$, one can find an element of the covering $\mathbf{s}_i$ such that $\|\mathbf{s} - \mathbf{s}_i\|_E \le \delta$. Assuming (C.1) holds, we have

$$ \|\Psi \mathbf{s}\|_F \ge \|\Psi \mathbf{s}_i\|_F - \|\Psi (\mathbf{s} - \mathbf{s}_i)\|_F \ge 1 - \frac{t}{2} - C\delta $$

and therefore by choosing $\delta = t/(2C)$ we obtain the desired result.

## Appendix D. Proof of Proposition 2

Fix any $\mathbf{x}_{\mathcal{S}} \in \mathcal{S}$. Let $\delta, \delta', \varepsilon > 0$ be small constants which values we shall define later, such that $\varepsilon \le \varepsilon_0$. Define $\mathcal{S}_{\varepsilon} = \mathcal{S} \setminus \{\mathbf{x} \in \mathcal{S}, d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}) \le \varepsilon\}$ the model with a ball around $\mathbf{x}_{\mathcal{S}}$ removed, and $N := \mathcal{N}(\mathcal{S}_{\varepsilon}, d_E, \delta)$, $N' := \mathcal{N}(\mathcal{S}_{\varepsilon}, d_E, \delta')$. Let $\mathbf{x}_1, \dots, \mathbf{x}_N$ and $\mathbf{s}_1, \dots, \mathbf{s}_N'$ be, respectively, a $\delta$-covering of $\mathcal{S}_{\varepsilon}$ and a $\delta'$-covering of $\mathcal{S}_{\varepsilon}$, with $\mathbf{x}'_i \in \mathcal{S}$ defined such that $\mathbf{s}_i = \frac{\mathbf{x}_{\mathcal{S}} - \mathbf{x}'_i}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}'_i)}$ and $0 < d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}'_i) \le \varepsilon$.

By the concentration property, it holds that, with probability at least $1 - (N + N')e^{-c(t/2)}$, for all $\tilde{\mathbf{x}} = \mathbf{x}_i$ or $\tilde{\mathbf{x}} = \mathbf{x}'_i$ for all indices $i$, we have

$$ \|\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\tilde{\mathbf{x}})\|_F \ge (1 - t/2) \cdot d_E(\mathbf{x}_{\mathcal{S}}, \tilde{\mathbf{x}}). \qquad (D.1) $$

Our goal is to extend this property to any element $\tilde{\mathbf{x}} = \mathbf{x}$ in the model $\mathcal{S}$.

Let $\mathbf{x} \in \mathcal{S}$ be any element of the model. We distinguish two cases. If $\mathbf{x} \in \mathcal{S}_{\varepsilon}$, i.e. $d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}) \ge \varepsilon$, we consider $\mathbf{x}_i \in \mathcal{S}_{\varepsilon}$ an element of the covering of $\mathcal{S}_{\varepsilon}$ such that $d_E(\mathbf{x}, \mathbf{x}_i) \le \delta$. We have:

$$ \begin{aligned} \left\| \frac{\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x})}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x})} - \frac{\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x}_i)}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}_i)} \right\|_F &\le \left\| \frac{\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x})}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x})} - \frac{\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x}_i)}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x})} \right\|_F + \left\| \frac{\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x}_i)}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x})} - \frac{\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x}_i)}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}_i)} \right\|_F \\ &= \frac{\|\Psi(\mathbf{x}) - \Psi(\mathbf{x}_i)\|_F}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x})} + \|\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x}_i)\|_F \left| \frac{1}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x})} - \frac{1}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}_i)} \right| \\ &\le \frac{C_1 \delta}{\varepsilon} + C_1 M_{\mathcal{S}} \frac{d_E(\mathbf{x}, \mathbf{x}_i)}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}) d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}_i)} \\ &\le \frac{C_1 \delta}{\varepsilon} \left( 1 + \frac{M_{\mathcal{S}}}{\varepsilon} \right) \end{aligned} $$

and therefore

$$ \left\| \frac{\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x})}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x})} \right\|_F \ge 1 - t/2 - \frac{C_1 \delta}{\varepsilon} \left( 1 + \frac{M_{\mathcal{S}}}{\varepsilon} \right) \qquad (D.2) $$

Now, when $d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}) \le \varepsilon$, we define $\mathbf{z} := \frac{\mathbf{x}_{\mathcal{S}} - \mathbf{x}}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x})}$ and note that $\mathbf{z} \in \mathcal{S}_{\varepsilon}$. We approximate it by an element of the covering of the normalized secant set $\mathbf{s}_i = \frac{\mathbf{x}_{\mathcal{S}} - \mathbf{x}_i'}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}_i)} \in \mathcal{S}_{\varepsilon}$ (meaning that $d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}_i') \le \varepsilon$) that verifies $d_E(\mathbf{z}, \mathbf{s}_i) \le \delta'$. Then, we have

$$ \begin{aligned} \left\| \frac{\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x})}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x})} - \frac{\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x}_i')}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}_i')} \right\|_F &\le \left\| \frac{\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x})}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x})} - D_{\Psi}\mathbf{z} \right\|_F + \|D_{\Psi}(\mathbf{z} - \mathbf{s}_i)\|_F + \left\| D_{\Psi}\mathbf{s}_i - \frac{\Psi(\mathbf{x}_{\mathcal{S}}) - \Psi(\mathbf{x}_i)}{d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}_i)} \right\|_F \\ &\stackrel{(iv)}{\le} C_2 d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}) + C_3 \delta' + C_2 d_E(\mathbf{x}_{\mathcal{S}}, \mathbf{x}_i') \le 2C_2\varepsilon + C_3\delta' \end{aligned} $$
and therefore

$$ \left\| \frac{\Psi(\mathbf{x}_{\mathfrak{S}}) - \Psi(\mathbf{x})}{d_E(\mathbf{x}_{\mathfrak{S}}, \mathbf{x})} \right\|_F \geq 1 - t/2 - (2C_2\varepsilon + C_3\delta') \quad (D.3) $$

To conclude, we set $\varepsilon \leq \min(\varepsilon_0, \frac{t}{8C_2})$, $\delta' \leq \frac{t}{4C_3}$, $\delta \leq \frac{t\varepsilon^2/(4C_1)}{\varepsilon+M_{\mathfrak{S}}}$ to obtain the desired result.

## Appendix E. Proof of Proposition 3

*Concentration property.* The concentration result is based on the fact that $\|\Psi(\mathbf{x}) - \Psi(\mathbf{x})\|_2^2 \approx E_{\omega} |\varphi_{\omega}(\mathbf{x}) - \varphi_{\omega}(\mathbf{x}')|^2 = d_E(\mathbf{x}, \mathbf{x}')^2$ by definition of the random features. Using simple function studies and a Bernstein concentration inequality with a control on moments of all orders, it is possible to show (see [18], eq. (160) then Prop. 6.11) that the concentration result (8) is valid with $c(t) \propto \frac{mt^2}{1+t}$ (we do not reproduce the detailed proof here for brevity).

Since this Bernstein inequality is in fact valid for all vectors (not necessarily in the model), as a consequence the Boundedness Property is also satisfied with constant $\beta = 1+t$, metric $d_E$ and probability $1 - e^{-c(t)}$.

We now check hypotheses (i)–(iv) in Proposition 2. We are going to repeatedly use the fact that for $\mathbf{x}, \mathbf{x}' \in \mathfrak{S}$, since $\|\mathbf{x}\|_2 \leq M$, we have

$$ \ell \|\mathbf{x} - \mathbf{x}'\|_2 \leq d_E(\mathbf{x}, \mathbf{x}') \leq L \|\mathbf{x} - \mathbf{x}'\|_2 \quad (E.1) $$

where $\ell := 2 \cdot \frac{1-\exp(-M^2/(2\sigma^2))}{M}$ and $L = \frac{M}{\sigma^2}$.

(i) Using a first-order Taylor expansion we have $\|\Psi(\mathbf{x}) - \Psi(\mathbf{x}')\|_2 \leq (\sum_{j=1}^m \frac{\|\omega_j\|_2^2}{mf(\omega_j)^2})^{\frac{1}{2}} \|\mathbf{x} - \mathbf{x}'\|_2 \lesssim \sqrt{\gamma_2} \|\mathbf{x} - \mathbf{x}'\|_2$ hence by (E.1) hypothesis (i) is valid with $C_1 \propto \frac{\sqrt{\gamma_2}}{\ell}$.

(ii) It is immediate that $M_{\mathfrak{S}} = ML$. Then, using the well-known fact that for each $S_j$ we have $\mathcal{N}(S_j \cap B_M, \|\cdot\|_2, \delta) = \mathcal{O}((\frac{M}{\delta})^s)$, by a union bound and (E.1) we have $\mathcal{N}(\mathfrak{S}, d_E, \delta) \leq \mathcal{O}(N(\frac{LM}{\delta})^s)$.

(iii) Similar to the model, for all $\varepsilon$ the normalized secant set is included in a union of subspaces $S_{\varepsilon} \subset \bigcup_{i,j=1}^N S_{ij} \cap B_{M/\ell}$, with norm controlled by $M/\ell$ by (E.1), and where $S_{ij} = S_i + S_j$ (sum of subspaces). Hence, since $\mathcal{N}(S_{ij} \cap B_{M/\ell}, \|\cdot\|_2, \delta) = \mathcal{O}((\frac{M}{\ell\delta})^{2s})$, by a union bound and (E.1) we have $\mathcal{N}(S_{\varepsilon}, d_E, \delta) \leq \mathcal{O}(N^2(\frac{LM}{\ell\delta})^{2s})$.

(iv) Finally, by a Taylor expansion we have

$$ \varphi_{\omega}(\mathbf{x}) - \varphi_{\omega}(\mathbf{x}') = \frac{i\omega^{\top}}{f(\omega)}(\mathbf{x} - \mathbf{x}')e^{i\omega^{\top}\mathbf{x}} - \frac{1}{2} \cdot \frac{(\omega^{\top}(\mathbf{x} - \mathbf{x}'))^2}{f(\omega)}e^{i\omega^{\top}\mathbf{x}''} $$

where $\mathbf{x}'' = \theta\mathbf{x} + (1-\theta)\mathbf{x}'$. Hence hypothesis (iv) is satisfied with $C_2 = C_1$ and $C_3 = \frac{\sqrt{\gamma_4}}{\ell^2}$, which concludes the proof.