∎

11institutetext:  Juan Guillermo Garrido 22institutetext: Departamento de
Ingeniería Matemática, Universidad de Chile, Santiago, Chile.22email:
<EMAIL_ADDRESS>
Emilio Vilches 33institutetext: Instituto de Ciencias de la Ingeniería,
Universidad de O’Higgins, Rancagua, Chile. 33email<EMAIL_ADDRESS>
# Catching-up Algorithm with Approximate Projections for Moreau’s Sweeping
Processes ††thanks:

Juan Guillermo Garrido and Emilio Vilches

(Received: date / Accepted: date)

###### Abstract

In this paper, we develop an enhanced version of the catching-up algorithm for
sweeping processes through an appropriate concept of approximate projections.
We establish some properties of this notion of approximate projection. Then,
under suitable assumptions, we show the convergence of the enhanced catching-
up algorithm for prox-regular, subsmooth, and merely closed sets. Finally, we
discuss efficient numerical methods for obtaining approximate projections. Our
results recover classical existence results in the literature and provide new
insight into the numerical simulation of sweeping processes.

###### Keywords:

Sweeping process differential inclusions normal cone approximate projections

###### MSC:

34A60 49J52 34G25 49J53 93D30

##  1 Introduction

Given a Hilbert space $\mathcal{H}$, Moreau’s sweeping process is a first-
order differential inclusion involving the normal cone to a family of closed
moving sets $(C(t))_{t\in[0,T]}$. In its simplest form, it can be written as

$\left\\{\begin{aligned} \dot{x}(t)&\in-N\left(C(t);x(t)\right)&\textrm{ a.e.
}t\in[0,T],\\\ x(0)&=x_{0}\in C(0),\end{aligned}\right.$ (SP)

where $N(C(t);\cdot)$ denotes an appropriate normal cone to the sets
$(C(t))_{t\in[0,T]}$. Since its introduction by J.J. Moreau in MO1 ; MO2 , the
sweeping process has allowed the development of various applications in
contact mechanics, electrical circuits, and crowd motion, among others (see,
e.g., Brogliato-M ; Acary-Bon-Bro-2011 ; Maury-Venel ). Furthermore, so far,
we have a well-consolidated existence theory for moving sets in the
considerable class of prox-regular sets.

The most prominent (and constructive) method for finding a solution to the
sweeping process is the so-called _catching-up algorithm_. Developed by J.J.
Moreau in MO2  for convex moving sets, it consists in taking a time
discretization $\\{t_{k}^{n}\\}_{k=0}^{n}$ of the interval $[0,T]$ and
defining a piecewise linear and continuous function
$x_{n}\colon[0,T]\to\mathcal{H}$ with nodes

$x_{k+1}^{n}:=\operatorname{proj}_{C(t_{k+1}^{n})}(x_{k}^{n})\textrm{ for all
}k\in\\{0,\ldots,n-1\\}.$

Moreover, under general assumptions, it could be proved that the sequence
$(x_{n})_{n}$ converges to the unique solution of (SP) (see, e.g., MR2159846
).

The applicability, from the numerical point of view, of the catching-up
algorithm is based on the possibility of calculating an exact formula for the
projection to the moving sets. However, for the majority of sets, the
projection onto a closed set is impossible to obtain exactly, and only
numerical approximations of it are possible. Since there are still no
guarantees on the convergence of the catching-up algorithm with approximate
projections, in this paper, we develop a theoretical framework for the
numerical approximation of the solutions of the sweeping process using an
appropriate concept of approximate projection that is consistent with the
numerical methods for the computation of the projection on a closed set.

Regarding numerical approximations of sweeping processes, we are aware of the
paper MR2800713 , where the author proposes an implementable numerical method
for the particular case of the intersection of the complement of convex sets,
which is used to study crowd motion. Our approach follows a different path and
is based on numerical optimization methods to find an approximate projection
in the following sense: given a closed set $C\subset\mathcal{H}$,
$\varepsilon>0$ and $x\in\mathcal{H}$, we say that $\bar{x}\in C$ is an
_approximate projection_ of $C$ at $x\in\mathcal{H}$ if

$\|x-\bar{x}\|^{2}<\inf_{y\in C}\|x-y\|^{2}+\varepsilon.$

We observe that the set of approximate projections is always nonempty and can
be obtained through numerical methods for optimization. Hence, in this paper,
we study the properties of approximate projections and propose a general
numerical method for the sweeping process based on approximate projections. We
prove that this algorithm converges in three general cases: (i) prox-regular
moving sets (without compactness assumptions), (ii) ball-compact subsmooth
moving sets, and (iii) general ball-compact fixed closed sets. Hence, our
results cover a wide range of existence results for the sweeping process. It
is worth emphasizing that our method generalizes the catching-up algorithm and
provides important insights into the numerical simulation of sweeping
processes.

The paper is organized as follows. Section 2 provides the mathematical tools
needed for the presentation of the paper and also develops the theoretical
properties of approximate projections. Section 3 is devoted to presenting the
proposed algorithm and its main properties. Then, in Section 4, we prove the
convergence of the algorithm when the moving set has uniformly prox-regular
values (without compactness assumptions). Next, in Section 5, we provide the
convergence of the proposed algorithm for ball-compact subsmooth moving sets.
Section 6 shows the convergence for a fixed ball-compact set. Finally, in
Section 7, we discuss numerical aspects for obtaining approximate projections.
The paper ends with concluding remarks.

##  2 Preliminaries

From now on $\mathcal{H}$ stands for a Hilbert space, whose norm, denoted by
$\|\cdot\|$, is induced by an inner product $\langle\cdot,\cdot\rangle$. The
closed (resp. open) ball centered at $x$ with radius $r>0$ is denoted by
$\mathbb{B}[x,r]$ (resp. $\mathbb{B}(x,r)$), and the closed unit ball is
denoted by $\mathbb{B}$. The notation $\mathcal{H}_{w}$ stands for
$\mathcal{H}$ equipped with the weak topology and $x_{n}\rightharpoonup x$
denotes the weak convergence of a sequence $(x_{n})_{n}$ to $x$. For a given
set $S\subset\mathcal{H}$, the _support_ and the _distance function_ of $S$ of
at $x\in\mathcal{H}$ are defined, respectively, as

$\sigma(x,S):=\sup_{z\in S}\langle x,z\rangle\textrm{ and
}d_{S}(x):=\inf_{z\in S}\|x-z\|.$

Given $\rho\in]0,+\infty]$ and $\gamma<1$ positive, the _$\rho$ -enlargement_
and the _$\gamma\rho$ -enlargement_ of $S$ are defined, respectively, as

$U_{\rho}(S)=\\{x\in\mathcal{H}:d_{S}(x)<\rho\\}\textrm{ and
}U_{\rho}^{\gamma}(S):=\\{x\in\mathcal{H}:d_{S}(x)<\gamma\rho\\}.$

Given $A,B\subset\mathcal{H}$ two sets, we define the _excess_ of $A$ over $B$
as the quantity $e(A,B):=\sup_{x\in A}d_{B}(x)$. From this, we define the
_Hausdorff distance_ between $A$ and $B$ as

$d_{H}(A,B):=\max\\{e(A,B),e(B,A)\\}.$

Further properties about Hausdorff distance can be found in (MR2378491, , Sec.
3.16).
A vector $h\in\mathcal{H}$ belongs to the Clarke tangent cone $T(S;x)$ (see
Clarke1983 ); when for every sequence $(x_{n})_{n}$ in $S$ converging to $x$
and every sequence of positive numbers $(t_{n})_{n}$ converging to $0$, there
exists a sequence $(h_{n})_{n}$ in $\mathcal{H}$ converging to $h$ such that
$x_{n}+t_{n}h_{n}\in S$ for all $n\in\mathbb{N}$. This cone is closed and
convex, and its negative polar $N(S;x)$ is the Clarke normal cone to $S$ at
$x\in S$, that is,

$N\left(S;x\right):=\left\\{v\in\mathcal{H}:\left\langle v,h\right\rangle\leq
0\quad\textrm{for all }h\in T(S;x)\right\\}.$

As usual, $N(S;x)=\emptyset$ if $x\notin S$. Through that normal cone, the
_Clarke subdifferential_ of a function
$f\colon\mathcal{H}\to\mathbb{R}\cup\\{+\infty\\}$ is defined by

$\partial f(x):=\left\\{v\in\mathcal{H}:(v,-1)\in
N\left(\operatorname{epi}f,(x,f(x))\right)\right\\},$

where
$\operatorname{epi}f:=\left\\{(y,r)\in\mathcal{H}\times\mathbb{R}:f(y)\leq
r\right\\}$ is the epigraph of $f$. When the function $f$ is finite and
locally Lipschitzian around $x$, the Clarke subdifferential is characterized
(see Clarke1998 ) in the following simple and amenable way

$\partial f(x)=\left\\{v\in\mathcal{H}:\left\langle v,h\right\rangle\leq
f^{\circ}(x;h)\textrm{ for all }h\in\mathcal{H}\right\\},$

where

$f^{\circ}(x;h):=\limsup_{(t,y)\to(0^{+},x)}t^{-1}\left[f(y+th)-f(y)\right],$

is the _generalized directional derivative_ of the locally Lipschitzian
function $f$ at $x$ in the direction $h\in\mathcal{H}$. The function
$f^{\circ}(x;\cdot)$ is in fact the support of $\partial f(x)$. That
characterization easily yields that the Clarke subdifferential of any locally
Lipschitzian function is a set-valued map with nonempty and convex values
satisfying the important property of upper semicontinuity from $\mathcal{H}$
into $\mathcal{H}_{w}$.

Let $f\colon\mathcal{H}\to\mathbb{R}\cup\\{+\infty\\}$ be an lsc (lower
semicontinuous) function and $x\in\operatorname{dom}f$. We say that

  * (i)

An element $\zeta$ belongs to the _proximal subdifferential_ of $f$ at $x$,
denoted by $\partial_{P}f(x)$, if there exist two nonnegative numbers $\sigma$
and $\eta$ such that

$f(y)\geq f(x)+\left\langle\zeta,y-x\right\rangle-\sigma\|y-x\|^{2}\textrm{
for all }y\in\mathbb{B}(x;\eta).$

  * (ii)

An element $\zeta\in\mathcal{H}$ belongs to the _Fréchet subdifferential_ of
$f$ at $x$, denoted by $\partial_{F}f(x)$, if

$\liminf_{h\to 0}\frac{f(x+h)-f(x)-\langle\zeta,h\rangle}{\|h\|}\geq 0.$

  * (iii)

An element $\zeta$ belongs to the _limiting subdifferential_ of $f$ at $x$,
denoted by $\partial_{L}f(x)$, if there exist sequences $(\zeta_{n})$ and
$(x_{n})$ such that $\zeta_{n}\in\partial_{P}f(x_{n})$ for all
$n\in\mathbb{N}$ and $x_{n}\to x$, $\zeta_{n}\rightharpoonup\zeta$, and
$f(x_{n})\to f(x)$.

Through these concepts, we can define the proximal, Fréchet, and limiting
normal cone of a given set $S\subset\mathcal{H}$ at $x\in S$, respectively, as

$N^{P}\left(S;x\right):=\partial_{P}I_{S}(x),\,N^{F}(C;x):=\partial_{F}I_{C}(x)\textrm{
and }N^{L}(S;x)=\partial_{L}I_{S}(x),$

where $I_{S}$ is the indicator function of $S\subset\mathcal{H}$ (recall that
$I_{S}(x)=0$ if $x\in S$ and $I_{S}(x)=+\infty$ if $x\notin S$). It is well-
known that

$N^{P}(S;x)\cap\mathbb{B}=\partial_{P}d_{S}(x)\quad\textrm{ for all }x\in S.$
(1)

The equality (see Clarke1998 )

$\displaystyle N\left(S;x\right)$
$\displaystyle=\overline{\text{co}}^{\ast}N^{L}(S;x)=\operatorname{cl}^{*}\left(\mathbb{R}_{+}\partial
d_{S}(x)\right)$ $\displaystyle\textrm{ for }x\in S,$

gives an expression of the Clarke normal cone in terms of the distance
function.

Now, we recall the concept of uniformly prox-regular sets. Introduced by
Federer in MR110078  and later developed by Rockafellar, Poliquin, and
Thibault in MR1694378 . The prox-regularity generalizes and unifies convexity
and nonconvex bodies with $C^{2}$ boundary. We refer to MR2768810 ;
Thibault-2023-II  for a survey.

###### Definition 1

Let $S$ be a closed subset of $\mathcal{H}$ and $\rho\in]0,+\infty]$. The set
$S$ is called $\rho$-uniformly prox-regular if for all $x\in S$ and $\zeta\in
N^{P}(S;x)$ one has

$\langle\zeta,x^{\prime}-x\rangle\leq\frac{\|\zeta\|}{2\rho}\|x^{\prime}-x\|^{2}\textrm{
for all }x^{\prime}\in S.$

It is important to emphasize that convex sets are $\rho$-uniformly prox-
regular for any $\rho>0$. The following proposition provides a
characterization of uniformly prox-regular sets (see, e.g., MR2768810 ).

###### Proposition 1

Let $S\subset\mathcal{H}$ be a closed set and $\rho\in]0,+\infty]$. The
following assertions are equivalent:

  1. (a)

$S$ is $\rho$-uniformly prox-regular.

  2. (b)

For any positive $\gamma<1$ the mapping $\operatorname{proj}_{S}$ is well-
defined on $U_{\rho}^{\gamma}(S)$ and Lipschitz continuous on
$U_{\rho}^{\gamma}(S)$ with $(1-\gamma)^{-1}$ as a Lipschitz constant, i.e.,

$\left\|\operatorname{proj}_{S}\left(u_{1}\right)-\operatorname{proj}_{S}\left(u_{2}\right)\right\|\leq(1-\gamma)^{-1}\left\|u_{1}-u_{2}\right\|\quad$

for all $u_{1},u_{2}\in U_{\rho}^{\gamma}(S)$.

  3. (c)

For any $x_{i}\in S,v_{i}\in N^{P}\left(S;x_{i}\right)\cap\mathbb{B}$, with
$i=1,2$, one has

$\left\langle
v_{1}-v_{2},x_{1}-x_{2}\right\rangle\geq-\frac{1}{\rho}\left\|x_{1}-x_{2}\right\|^{2},$

that is, the set-valued mapping $N^{P}(S;\cdot)\cap\mathbb{B}$ is
$1/\rho$-hypomonotone.

  4. (d)

For all $\gamma\in]0,1[$, for all $x_{1},x_{2}\in U_{\rho}^{\gamma}(S)$, for
all $\xi\in\partial_{P}d_{S}(x)$, one has

$\langle\xi,x_{2}-x_{1}\rangle\leq\frac{1}{2\rho(1-\gamma)^{2}}\|x_{1}-x_{2}\|^{2}+d_{S}(x_{2})-d_{S}(x_{1}).$

Next, we recall the class of subsmooth sets that includes the concepts of
convex and uniformly prox-regular sets (see MR2115366  and also
(Thibault-2023-II, , Chapter 8) for a survey).

###### Definition 2

Let $S$ be a closed subset of $\mathcal{H}$. We say that $S$ is subsmooth at
$x_{0}\in S$, if for every $\varepsilon>0$ there exists $\delta>0$ such that

$\left\langle\xi_{2}-\xi_{1},x_{2}-x_{1}\right\rangle\geq-\varepsilon\left\|x_{2}-x_{1}\right\|,$
(2)

whenever $x_{1},x_{2}\in\mathbb{B}\left[x_{0},\delta\right]\cap S$ and
$\xi_{i}\in N\left(S;x_{i}\right)\cap\mathbb{B}$ for $i\in\\{1,2\\}$. The set
$S$ is said _subsmooth_ if it is subsmooth at each point of $S$. We further
say that $S$ is uniformly subsmooth, if for every $\varepsilon>0$ there exists
$\delta>0$, such that (2) holds for all $x_{1},x_{2}\in S$ satisfying
$\left\|x_{1}-x_{2}\right\|\leq\delta$ and all $\xi_{i}\in
N\left(S;x_{i}\right)\cap\mathbb{B}$ for $i\in\\{1,2\\}$.
Let $(S(t))_{t\in I}$ be a family of closed sets of $\mathcal{H}$ indexed by a
nonemptyset $I$. The family is called equi-uniformly subsmooth, if for all
$\varepsilon>0$, there exists $\delta>0$ such that for all $t\in I$, the
inequality (2) holds for all $x_{1},x_{2}\in S(t)$ satisfying
$\|x_{1}-x_{2}\|\leq\delta$ and all $\xi_{i}\in N(S(t);x_{i})\cap\mathbb{B}$
with $i\in\\{1,2\\}$.

Given an interval $\mathcal{I}$, a set-valued map
$F\colon\mathcal{I}\rightrightarrows\mathcal{H}$ is said to be measurable if
for all open set $U$ of $\mathcal{H}$, the inverse image
$F^{-1}(U)=\\{t\in\mathcal{I}:F(t)\cap U\neq\emptyset\\}$ is a Lebesgue
measurable set. When $F$ takes nonempty and closed values and $\mathcal{H}$ is
separable, this notion is equivalent to the
$\mathcal{L}\otimes\mathcal{B}(\mathcal{H})$-measurability of the graph
$\operatorname{gph}F:=\\{(t,x)\in\mathcal{I}\times\mathcal{H}:x\in F(t)\\}$
(see, e.g., (MR2527754, , Theorem 6.2.20)).

Given a set-valued map $F\colon\mathcal{H}\rightrightarrows\mathcal{H}$, we
say $F$ is upper semicontinuous from $\mathcal{H}$ into $\mathcal{H}_{w}$ if
for all weakly closed set $C$ of $\mathcal{H}$, the inverse image $F^{-1}(C)$
is a closed set of $\mathcal{H}$. It is known (see, e.g., see (MR2527754, ,
Proposition 6.1.15 (c))) that if $F$ is upper semicontinuous, then the map
$x\mapsto\sigma(\xi,F(x))$ is upper semicontinuous for all
$\xi\in\mathcal{H}$. When $F$ takes convex and weakly compact values, these
two properties are equivalent (see (MR2527754, , Proposition 6.1.17)).
A set $S\subset\mathcal{H}$ is said ball compact if the set $S\cap
r\mathbb{B}$ is compact for all $r>0$. The _projection_ onto
$S\subset\mathcal{H}$ is the (possibly empty) set

$\operatorname{Proj}_{S}(x):=\left\\{z\in S:d_{S}(x)=\|x-z\|\right\\}.$

When the projection set is a singleton, we denote it as
$\operatorname{proj}_{S}(x)$. For $\varepsilon>0$, we define the set of
_approximate projections_ :

$\operatorname{proj}_{S}^{\varepsilon}(x):=\left\\{z\in
S:\|x-z\|^{2}<d_{S}^{2}(x)+\varepsilon\right\\}.$

By definition, the above set is nonempty and open. Moreover, it satisfies
similar properties as the projection map (see Proposition 2 below). The
approximate projections have been considered several times in variational
analysis. In particular, they were used to characterize the subdifferential of
the Asplund function of a given set. Indeed, let $S\subset\mathcal{H}$ and
consider the Asplund function of the set $S$

$\varphi_{S}(x):=\frac{1}{2}\|x\|^{2}-\frac{1}{2}d_{S}^{2}(x)\quad
x\in\mathcal{H}.$

Then, the following formula holds (see, e.g., (MR2848527, , p. 467)):

$\partial\varphi_{S}(x)=\bigcap_{\varepsilon>0}\overline{\text{co}}(\operatorname{proj}_{S}^{\varepsilon}(x)).$

We recall that for any set $S\subset\mathcal{H}$ and $x\in\mathcal{H}$, where
$\operatorname{Proj}_{S}(x)\neq\emptyset$, the following formula holds:

$x-z\in d_{S}(x)\partial_{P}d_{S}(z)\textrm{ for all
}z\in\operatorname{Proj}_{S}(x).$

The next result provides an approximate version of the above formula for any
closed set $S\subset\mathcal{H}$.

###### Lemma 1

Let $S\subset\mathcal{H}$ be a closed set, $x\in\mathcal{H}$, and
$\varepsilon>0$. For each $z\in\operatorname{proj}_{S}^{\varepsilon}(x)$ there
is $v\in\operatorname{proj}_{S}^{\varepsilon}(x)$ such that
$\|z-v\|<2\sqrt{\varepsilon}$ and

$x-z\in(4\sqrt{\varepsilon}+d_{S}(x))\partial_{P}d_{S}(v)+3\sqrt{\varepsilon}\mathbb{B}.$

###### Proof

Fix $\varepsilon>0$, $x\in\mathcal{H}$ and
$z\in\operatorname{proj}_{S}^{\varepsilon}(x)$. According to the Borwein-
Preiss Variational Principle (MR902782, , Theorem 2.6) applied to $y\mapsto
g(y):=\|x-y\|^{2}+I_{S}(y)$, there exists
$v\in\operatorname{proj}_{S}^{\varepsilon}(x)$ such that
$\|z-v\|<2\sqrt{\varepsilon}$ and
$0\in\partial_{P}g(v)+2\sqrt{\varepsilon}\mathbb{B}$. Then, by the sum rule
for the proximal subdifferential (see, e.g., (Clarke1998, , Proposition
2.11)), we obtain that

$x-v\in N^{P}(S;v)+\sqrt{\varepsilon}\mathbb{B},$

which implies that $x-z\in N^{P}(S;v)+3\sqrt{\varepsilon}\mathbb{B}.$ Next,
since $\|x-z\|\leq d_{S}(x)+\sqrt{\varepsilon}$, we obtain that

$x-z\in
N^{P}(S;v)\cap(4\sqrt{\varepsilon}+d_{S}(x))\mathbb{B}+3\sqrt{\varepsilon}\mathbb{B}.$

Finally, the result follows from formula (1) and the above inclusion.

The following proposition displays some properties of approximation
projections for uniformly prox-regular sets.

###### Proposition 2

Let $S\subset\mathcal{H}$ be a $\rho$-uniformly prox-regular set. Then, one
has:

  1. (a)

Let $(x_{n})_{n}$ be a sequence converging to $x\in U_{\rho}(S)$. Then for any
$(z_{n})_{n}$ and any sequence of positive numbers $(\varepsilon_{n})_{n}$
converging to $0$ with
$z_{n}\in\operatorname{proj}_{S}^{\varepsilon_{n}}(x_{n})$ for all
$n\in\mathbb{N}$, we have that $z_{n}\to\operatorname{proj}_{S}(x)$.

  2. (b)

Let $\gamma\in]0,1[$ and $\varepsilon\in]0,\varepsilon_{0}]$ where
$\varepsilon_{0}$ is such that

$\gamma+4\sqrt{\varepsilon_{0}}\left(1+\gamma+\frac{1}{\rho}(1+4\sqrt{\varepsilon_{0}})\right)=1.$

Then, for all $z_{i}\in\operatorname{proj}_{S}^{\varepsilon}(x_{i})$ with
$x_{i}\in U_{\rho}^{\gamma}(S)$ for $i\in\\{1,2\\}$, we have

$(1-\digamma)\|z_{1}-z_{2}\|^{2}\leq\sqrt{\varepsilon}\|x_{1}-x_{2}\|^{2}+M\sqrt{\varepsilon}+\langle
x_{1}-x_{2},z_{1}-z_{2}\rangle,$

where
$\digamma:=\frac{\alpha}{\rho}+4\sqrt{\varepsilon}\left(1+\frac{\alpha}{\rho}+\frac{1}{\rho}(1+\sqrt{\varepsilon})\right)$
with $\alpha:=\max\\{d_{S}(x_{1}),d_{S}(x_{2})\\}$ and $M$ is a nonnegative
constant only dependent on $\varepsilon,\rho,\gamma$.

###### Proof

$(a)$ We observe that for all $n\in\mathbb{N}$

$\displaystyle\|z_{n}\|$ $\displaystyle\leq\|z_{n}-x_{n}\|+\|x_{n}\|\leq
d_{C}(x_{n})+\sqrt{\varepsilon_{n}}+\|x_{n}\|.$

Hence, since $\varepsilon_{n}\to 0$ and $x_{n}\to x$, we obtain $(y_{n})_{n}$
is bounded. On the other hand, since $x\in U_{\rho}(S)$, we obtain
$\operatorname{proj}_{S}(x)$ is well-defined and

$\displaystyle\|z_{n}-\operatorname{proj}_{S}(x)\|^{2}$
$\displaystyle=\|z_{n}-x_{n}\|^{2}-\|x_{n}-\operatorname{proj}_{S}(x)\|^{2}$
$\displaystyle+2\langle
x-\operatorname{proj}_{S}(x),z_{n}-\operatorname{proj}_{S}(x)\rangle+2\langle
z_{n}-\operatorname{proj}_{S}(x),x_{n}-x\rangle$ $\displaystyle\leq
d_{S}^{2}(x_{n})+\varepsilon_{n}-\|x_{n}-\operatorname{proj}_{S}(x)\|^{2}$
$\displaystyle+2\langle
x-\operatorname{proj}_{S}(x),z_{n}-\operatorname{proj}_{S}(x)\rangle+2\langle
z_{n}-\operatorname{proj}_{S}(x),x_{n}-x\rangle$
$\displaystyle\leq\varepsilon_{n}+2\langle
x-\operatorname{proj}_{S}(x),z_{n}-\operatorname{proj}_{S}(x)\rangle$
$\displaystyle+2\langle z_{n}-\operatorname{proj}_{S}(x),x_{n}-x\rangle$

where we have used $z_{n}\in\operatorname{proj}_{S}^{\varepsilon_{n}}(x_{n})$
and that $d_{S}^{2}(x_{n})\leq\|x_{n}-\operatorname{proj}_{S}(x)\|^{2}$.
Moreover, since $x-\operatorname{proj}_{S}(x)\in
N^{P}(S;\operatorname{proj}_{S}(x))$ and $S$ is $\rho$-uniformly prox-regular,
we obtain that

$2\langle
x-\operatorname{proj}_{S}(x),z_{n}-\operatorname{proj}_{s}(x)\rangle\leq\frac{d_{S}(x)}{\rho}\|z_{n}-\operatorname{proj}_{S}(x)\|^{2}.$

Therefore, by using the above inequality and rearranging terms, we obtain that

$\|z_{n}-\operatorname{proj}_{S}(x)\|^{2}\leq\frac{\rho}{\rho-
d_{S}(x)}\left(\varepsilon_{n}+2\langle
z_{n}-\operatorname{proj}_{S}(x),x_{n}-x\rangle\right).$

Finally, since $x_{n}\to x$ and $(z_{n})_{n}$ is bounded, we conclude that
$z_{n}\to\operatorname{proj}_{S}(x)$.
$(b)$ By virtue of Lemma 1, for $i\in\\{1,2\\}$ there exists
$v_{i},b_{i}\in\mathcal{H}$ such that

$b_{i}\in\mathbb{B},v_{i}\in\operatorname{proj}_{S}^{\varepsilon}(x_{i}),\|z_{i}-v_{i}\|<2\sqrt{\varepsilon}\textrm{
and
}\frac{x_{i}-z_{i}-3\sqrt{\varepsilon}b_{i}}{4\sqrt{\varepsilon}+d_{S}(x_{i})}\in\partial_{P}d_{S}(v_{i}).$

The hypomonotonicity of proximal normal cone (see Proposition 1 (b)) implies
that

$\left\langle\zeta_{1}-\zeta_{2},v_{1}-v_{2}\right\rangle\geq\frac{-1}{\rho}\|v_{1}-v_{2}\|^{2},$

where
$\zeta_{i}:=\frac{x_{i}-z_{i}-3\sqrt{\varepsilon}b_{i}}{4\sqrt{\varepsilon}+\alpha}$
for $i\in\\{1,2\\}$ and $\alpha:=\max\\{d_{S}(x_{1}),d_{S}(x_{2})\\}$. On the
one hand, we have

$\|v_{1}-v_{2}\|\leq\|v_{1}-z_{1}\|+\|z_{1}-z_{2}\|+\|z_{2}-v_{2}\|\leq
4\sqrt{\varepsilon}+\|z_{1}-z_{2}\|,$

and for all $z\in\mathcal{H}$ and $i\in\\{1,2\\}$

$|\langle
z,v_{i}-z_{i}\rangle|\leq\frac{\sqrt{\varepsilon}\|z\|^{2}}{2}+\frac{\|v_{i}-z_{i}\|^{2}}{2\sqrt{\varepsilon}}\leq\frac{\sqrt{\varepsilon}\|z\|^{2}}{2}+2\sqrt{\varepsilon}.$

On the other hand,

$\displaystyle\langle(x_{1}-z_{1}-3\sqrt{\varepsilon}b_{1})-(x_{2}-z_{2}-3\sqrt{\varepsilon}b_{2}),v_{1}-v_{2}\rangle$
$\displaystyle=$ $\displaystyle\ 3\sqrt{\varepsilon}\langle
b_{2}-b_{1},v_{1}-v_{2}\rangle+\langle(x_{1}-x_{2})-(z_{1}-z_{2}),v_{1}-v_{2}\rangle$
$\displaystyle=$ $\displaystyle\ 3\sqrt{\varepsilon}\langle
b_{2}-b_{1},v_{1}-v_{2}\rangle+\langle x_{1}-x_{2},v_{1}-z_{1}\rangle+\langle
x_{1}-x_{2},z_{1}-z_{2}\rangle$ $\displaystyle+\langle
x_{1}-x_{2},z_{2}-v_{2}\rangle-\langle
z_{1}-z_{2},v_{1}-z_{1}\rangle-\|z_{1}-z_{2}\|^{2}-\langle
z_{1}-z_{2},z_{2}-v_{2}\rangle$ $\displaystyle\leq$ $\displaystyle\
6\sqrt{\varepsilon}(4\sqrt{\varepsilon}+\|z_{1}-z_{2}\|)+\sqrt{\varepsilon}\|x_{1}-x_{2}\|^{2}+8\sqrt{\varepsilon}+\langle
x_{1}-x_{2},z_{1}-z_{2}\rangle$
$\displaystyle-(1-\sqrt{\varepsilon})\|z_{1}-z_{2}\|^{2}$ $\displaystyle\leq$
$\displaystyle\
24\varepsilon+11\sqrt{\varepsilon}+\sqrt{\varepsilon}\|x_{1}-x_{2}\|^{2}+\langle
x_{1}-x_{2},z_{1}-z_{2}\rangle-(1-4\sqrt{\varepsilon})\|z_{1}-z_{2}\|^{2}.$

It follows that

$\displaystyle\left[1-\frac{\alpha}{\rho}-4\sqrt{\varepsilon}(1+\frac{1}{\rho}(1+4\sqrt{\varepsilon}+\alpha))\right]\|z_{1}-z_{2}\|^{2}$
$\displaystyle\leq$
$\displaystyle\sqrt{\varepsilon}\|x_{1}-x_{2}\|^{2}+\langle
x_{1}-x_{2},z_{1}-z_{2}\rangle+4(4\varepsilon+\sqrt{\varepsilon})(4\frac{\sqrt{\varepsilon}}{\rho}+\gamma)+24\varepsilon+11\sqrt{\varepsilon}$

which proves the desired inequality.

The following result provides a stability result for a family of equi-
uniformly subsmooth sets. We refer to see (MR3574145, , Lemma 2.7) for a
similar result.

###### Lemma 2

Let $\mathcal{C}=\\{C_{n}\\}_{n\in\mathbb{N}}\cup\\{C\\}$ be a family of
nonempty, closed and equi-uniformly subsmooth sets. Assume that

$\lim_{n\to\infty}d_{C_{n}}(x)=0,\text{ for all }x\in C.$

Then, for any sequence $\alpha_{n}\to\alpha\in\mathbb{R}$ and any sequence
$(y_{n})$ converging to $y$ with $y_{n}\in C_{n}$ and $y\in C$, one has

$\limsup_{n\to\infty}\sigma(\xi,\alpha_{n}\partial
d_{C_{n}}(y_{n}))\leq\sigma(\xi,\alpha\partial d_{C}(y))\text{ for all
}\xi\in\mathcal{H}.$

###### Proof

Fix $\xi\in\mathcal{H}$. Since $\partial d_{S}(x)\subset\mathbb{B}$ for all
$x\in\mathcal{H}$, we observe that

$\beta:=\limsup_{n\to\infty}\sigma(\xi,\alpha_{n}\partial
d_{C_{n}}(y_{n}))<+\infty.$

Let us consider a subsequence $(n_{k})_{k}$ such that

$\beta=\lim_{k\to\infty}\sigma(\xi,\alpha_{n_{k}}\partial
d_{C_{n_{k}}}(y_{n_{k}})).$

Given that $\partial d_{C_{n_{k}}}(y_{n_{k}})$ is weakly compact for all
$k\in\mathbb{N}$, there is $v_{n_{k}}\in\partial d_{C_{n_{k}}}(y_{n_{k}})$
such that

$\sigma(\xi,\alpha_{n_{k}}\partial
d_{C_{n_{k}}}(y_{n_{k}}))=\langle\xi,\alpha_{n_{k}}v_{n_{k}}\rangle\textrm{
for all }k\in\mathbb{N}.$

Moreover, the sequence $(v_{n_{k}})$ is bounded. Hence, without loss of
generality, we can assume that $v_{n_{k}}\rightharpoonup v\in\mathbb{B}$. It
follows that $\beta=\langle\xi,\alpha v\rangle$. By equi-uniformly
subsmoothness of $\mathcal{C}$, for any $\varepsilon>0$, there is $\delta>0$
such that for all $D\in\mathcal{C}$ and $x_{1},x_{2}\in D$ with
$\|x_{1}-x_{2}\|<\delta$, one has

$\langle\zeta_{1}-\zeta_{2},x_{1}-x_{2}\rangle\geq-\varepsilon\|x_{1}-x_{2}\|,$
(3)

whenever $\zeta_{i}\in N(D;x_{i})\cap\mathbb{B}$ for $i\in\\{1,2\\}$. Next,
let $y^{\prime}\in C$ such that $\|y-y^{\prime}\|<\delta/2$. Then, since
$d_{C_{n_{k}}}(y^{\prime})$ converges to $0$, there is a sequence
$(y_{n_{k}}^{\prime})$ converging to $y^{\prime}$ with $y_{n_{k}}^{\prime}\in
C_{n_{k}}$ for all $k\in\mathbb{N}$. Hence, there is $k_{0}\in\mathbb{N}$ such
that $\|y_{n_{k}}^{\prime}-y^{\prime}\|<\delta/2$ for all $k\geq k_{0}$. On
the other hand, since $y_{n}\to y$, then there is
$k_{0}^{\prime}\in\mathbb{N}$ such that $\|y_{n_{k}}-y\|<\delta/2$ for all
$k\geq k_{0}^{\prime}$. Hence, if
$k\geq\max\\{k_{0},k_{0}^{\prime}\\}=:\widehat{k}$ we have
$\|y_{n_{k}}-y_{n_{k}}^{\prime}\|<\delta$. Therefore, it follows from the fact
that $0\in\partial d_{C_{n_{k}}}(y_{n_{k}}^{\prime})$ and inequality (3) that

$\langle
v_{n_{k}},y_{n_{k}}-y_{n_{k}}^{\prime}\rangle\geq-\varepsilon\|y_{n_{k}}-y_{n_{k}}^{\prime}\|\text{
for all }k\geq\widehat{k}.$

By taking $k\to\infty$, we obtain that

$\langle v,y-y^{\prime}\rangle\geq-\varepsilon\|y-y^{\prime}\|\textrm{ for all
}y^{\prime}\in C\cap\mathbb{B}(y,\delta/2),$

which implies that $v\in N^{F}(C;y)$. Then, by (MR2986672, , Lemma 4.21),

$v\in N^{F}(C;y)\cap\mathbb{B}=\partial_{F}d_{C}(y)\subset\partial d_{C}(y).$

Finally, we have proved that

$\beta=\langle\xi,\alpha v\rangle\leq\sigma(\xi,\alpha\partial d_{C}(y)),$

which ends the proof.

The following lemma is a convergence theorem for a set-valued map from a
topological space into a Hilbert space.

###### Lemma 3

Let $(E,\tau)$ be a topological space and $\mathcal{G}\colon
E\rightrightarrows\mathcal{H}$ be a set-valued map with nonempty, closed, and
convex values. Consider sequences $(x_{n})_{n}\subset E$,
$(y_{n})_{n}\subset\mathcal{H}$ and
$(\varepsilon_{n})_{n}\subset\mathbb{R}_{+}$ such that

  1. (i)

$x_{n}\to x$ (in $E$), $y_{n}\rightharpoonup y$ (weakly in $\mathcal{H}$) and
$\varepsilon_{n}\to 0$;

  2. (ii)

For all $n\in\mathbb{N}$,
$y_{n}\in\text{co}(\mathcal{G}(x_{k})+\varepsilon_{k}\mathbb{B}:k\geq n)$;

  3. (iii)

$\displaystyle\limsup_{n\to\infty}\sigma(\xi,\mathcal{G}(x_{n}))\leq\sigma(\xi,\mathcal{G}(x))$
for all $\xi\in\mathcal{H}$.

Then, $y\in\mathcal{G}(x)$.

###### Proof

Assume by contradiction that $y\notin\mathcal{G}(x)$. By virtue of Hahn-Banach
Theorem there exists $\xi\in\mathcal{H}\setminus\\{0\\}$, $\delta>0$ and
$\alpha\in\mathbb{R}$ such that

$\langle\xi,y^{\prime}\rangle+\delta\leq\alpha\leq\langle\xi,y\rangle,\
\forall y^{\prime}\in\mathcal{G}(x).$

Then, it follows that $\sigma(\xi,\mathcal{G}(x))\leq\alpha-\delta$. Besides,
according to (ii) we have for all $n\in\mathbb{N}$, there is a finite set
$J_{n}\subset\mathbb{N}$ such that for all $m\in J_{n}$, $m\geq n$ and

$y_{n}=\sum_{j\in J_{n}}\alpha_{j}(y_{j}^{\prime}+\varepsilon_{j}v_{j})$

where for all $j\in J_{n}$, $\alpha_{j}\geq 0$, $v_{j}\in\mathbb{B}$,
$y_{j}^{\prime}\in\mathcal{G}(x_{j})$ and $\sum_{j\in J_{n}}\alpha_{j}=1$.
Also, there exists $N\in\mathbb{N}$ such that for all $n\geq N$,
$\varepsilon_{n}<\frac{\delta}{2\|\xi\|}$. Thus, for $n\geq N$

$\displaystyle\langle\xi,y_{n}\rangle$ $\displaystyle=\sum_{j\in
J_{n}}\alpha_{j}\langle\xi,y_{j}^{\prime}+\varepsilon_{j}v_{j}\rangle$
$\displaystyle\leq\sum_{j\in J_{n}}\alpha_{j}\sup_{k\geq
n}\sigma(\xi,\mathcal{G}(x_{k}))+\sum_{j\in
J_{n}}\alpha_{j}\varepsilon_{j}\langle\xi,v_{j}\rangle$
$\displaystyle\leq\sup_{k\geq
n}\sigma(\xi,\mathcal{G}(x_{k}))+\|\xi\|\sum_{j\in
J_{n}}\alpha_{j}\frac{\delta}{2\|\xi\|}\leq\sup_{k\geq
n}\sigma(\xi,\mathcal{G}(x_{k}))+\frac{\delta}{2}.$

Therefore, as $y_{n}\rightharpoonup y$, letting $n\to\infty$ in the last
inequality we obtain that

$\displaystyle\langle\xi,y\rangle$
$\displaystyle\leq\limsup_{n\to\infty}\sigma(\xi,\mathcal{G}(x_{n}))+\frac{\delta}{2}\leq\sigma(\xi,\mathcal{G}(x))+\frac{\delta}{2}$

Therefore,
$\langle\xi,y\rangle\leq\alpha-\delta/2\leq\langle\xi,y\rangle-\delta/2$,
which is a contradiction. The proof is then complete.

##  3 Catching-up algorithm with errors for sweeping processes

In this section, we propose a numerical method for the existence of solutions
for the sweeping process:

$\left\\{\begin{aligned}
\dot{x}(t)&\in-N\left(C(t);x(t)\right)+F(t,x(t))&\textrm{ a.e. }t\in[0,T],\\\
x(0)&=x_{0}\in C(0),\end{aligned}\right.$ (4)

where $C\colon[0,T]\rightrightarrows\mathcal{H}$ is a set-valued map with
closed values in a Hilbert space $\mathcal{H}$, $N\left(C(t);x\right)$ stands
for the Clarke normal cone to $C(t)$ at $x$, and
$F\colon[0,T]\times\mathcal{H}\rightrightarrows\mathcal{H}$ is a given set-
valued map with nonempty closed and convex values. Our algorithm is based on
the catching-up algorithm, except that we do not ask for an exact calculation
of the projections.

The proposed algorithm is given as follows. For $n\in\mathbb{N}^{\ast}$, let
$(t_{k}^{n})_{k=0}^{n}$ be a uniform partition of $[0,T]$ with uniform time
step $\mu_{n}:=T/n$. Let $(\varepsilon_{n})_{n}$ be a sequence of positive
numbers such that $\varepsilon_{n}/\mu_{n}^{2}\to 0$. We consider a sequence
of piecewise continuous linear approximations $(x_{n})_{n}$ defined as
$x_{n}(0)=x_{0}$ and for any $k\in\\{0,\ldots,n-1\\}$ and
$t\in]t_{k}^{n},t_{k+1}^{n}]$

$x_{n}(t)=x_{k}^{n}+\frac{t-t_{k}^{n}}{\mu_{n}}\left(x_{k+1}^{n}-x_{k}^{n}-\int_{t_{k}^{n}}^{t_{k+1}^{n}}f(s,x_{k}^{n})ds\right)+\int_{t_{k}^{n}}^{t}f(s,x_{k}^{n})ds,$
(5)

where $x_{0}^{n}=x_{0}$ and

$x_{k+1}^{n}\in\operatorname{proj}_{C(t_{k+1}^{n})}^{\varepsilon_{n}}\left(x_{k}^{n}+\int_{t_{k}^{n}}^{t_{k+1}^{n}}f(s,x_{k}^{n})ds\right)\
\text{for }k\in\\{0,1,...,n-1\\}.$ (6)

Here $f(t,x)$ denotes any selection of $F(t,x)$ such that $f(\cdot,x)$ is
measurable for all $x\in\mathcal{H}$. For simplicity, we consider
$f(t,x)\in\operatorname{proj}_{F(t,x)}^{\gamma}(0)$ for some $\gamma>0$. In
Proposition 3, we prove that it is possible to obtain such measurable
selection under mild assumptions.

The above algorithm is called catching-up algorithm with approximate
projections because the projection is not necessarily exactly calculated. We
will prove that the above algorithm converges for several families of
algorithms as long as the inclusion (6) is verified.

Let us consider functions $\delta_{n}(\cdot)$ and $\theta_{n}(\cdot)$ defined
as

$\delta_{n}(t)=\begin{cases}t_{k}^{n}&\textrm{ if
}t\in[t_{k}^{n},t_{k+1}^{n}[\\\ t_{n-1}^{n}&\textrm{ if }t=T,\end{cases}\text{
and }\theta_{n}(t)=\begin{cases}t_{k+1}^{n}&\textrm{ if
}t\in[t_{k}^{n},t_{k+1}^{n}[\\\ T&\textrm{ if }t=T.\end{cases}$

In what follows, we show useful properties satisfied for the above algorithm,
which will help us to prove the existence of the sweeping process (4) in three
cases:

  1. (i)

The set-valued map $t\rightrightarrows C(t)$ takes uniformly prox-regular
values.

  2. (ii)

The set-valued map $t\rightrightarrows C(t)$ takes subsmooth and ball-compact
values.

  3. (iii)

$C(t)\equiv C$ in $[0,T]$ and $C$ is ball-compact.

Throughout this section,
$F\colon[0,T]\times\mathcal{H}\rightrightarrows\mathcal{H}$ will be a set-
valued map with nonempty, closed, and convex values. Moreover, we will
consider the following conditions:

  1. $(\mathcal{H}_{1}^{F})$

For all $t\in[0,T]$, $F(t,\cdot)$ is upper semicontinuous from $\mathcal{H}$
into $\mathcal{H}_{w}$.

  2. $(\mathcal{H}_{2}^{F})$

There exists $h\colon\mathcal{H}\to\mathbb{R}^{+}$ Lipschitz continuous (with
constant $L_{h}>0$) such that

$\displaystyle d\left(0,F(t,x)\right):=\inf\\{\|w\|:w\in F(t,x)\\}\leq h(x),$

for all $x\in\mathcal{H}$ and a.e. $t\in[0,T]$.

  3. $(\mathcal{H}_{3}^{F})$

There is $\gamma>0$ such that the set-valued map
$(t,x)\rightrightarrows\operatorname{proj}_{F(t,x)}^{\gamma}(0)$ has a
selection $f\colon[0,T]\times\mathcal{H}\to\mathcal{H}$ such that $f(\cdot,x)$
is measurable for all $x\in\mathcal{H}$.

The following proposition provides conditions for the feasibility of
hypothesis $(\mathcal{H}_{3}^{F})$.

###### Proposition 3

Let us assume that $\mathcal{H}$ is a separable Hilbert space. Moreover we
suppose $F(\cdot,x)$ is measurable for all $x\in\mathcal{H}$, then
$(\mathcal{H}_{3}^{F})$ holds for all $\gamma>0$.

###### Proof

Let $\gamma>0$ and fix $x\in\mathcal{H}$. Since the set-valued map
$F(\cdot,x)$ is measurable, the map $t\mapsto d(0,F(t,x))$ is a measurable
function. Let us define the set-valued map $\mathcal{F}_{x}\colon
t\rightrightarrows\operatorname{proj}_{F(t,x)}^{\gamma}(0)$. Then,

$\displaystyle\operatorname{gph}\mathcal{F}_{x}$
$\displaystyle=\\{(t,y)\in[0,T]\times\mathcal{H}:y\in\operatorname{proj}_{F(t,x)}^{\gamma}(0)\\}$
$\displaystyle=\\{(t,y)\in[0,T]\times\mathcal{H}:\|y\|^{2}<d(0,F(t,x))^{2}+\gamma\text{
and }y\in F(t,x)\\}$
$\displaystyle=\operatorname{gph}F(\cdot,x)\cap\\{(t,y)\in[0,T]\times\mathcal{H}:\|y\|^{2}<d(0,F(t,x))^{2}+\gamma\\}.$

Hence, $\operatorname{gph}\mathcal{F}_{x}$ is a measurable set. Consequently,
$\mathcal{F}_{x}$ has a measurable selection (see (MR2527754, , Theorem
6.3.20)). Denoting by $t\mapsto f(t,x)$ such selection, we obtain the result.

Now, we establish the main properties of the proposed algorithm.

###### Theorem 3.1

Assume, in addition to $(\mathcal{H}_{1}^{F})$, $(\mathcal{H}_{2}^{F})$ and
$(\mathcal{H}_{3}^{F})$, that $C\colon[0,T]\rightrightarrows\mathcal{H}$ is a
set-valued map with nonempty and closed values such that

$d_{H}(C(t),C(s))\leq L_{C}|t-s|\text{ for all }t,s\in[0,T].$ (7)

Then, the sequence of functions
$(x_{n}\colon[0,T]\to\mathcal{H})_{n\in\mathbb{N}}$ generated by the numerical
scheme (5) and (6) satisfies the following properties:

  1. $(a)$

There are nonnegative constants $K_{1},K_{2},K_{3},K_{4},K_{5}$ such that for
all $n\in\mathbb{N}$ and $t\in[0,T]$:

    1. (i)

$d_{C(\theta_{n}(t))}(x_{n}(\delta_{n}(t))+\int_{\delta_{n}(t)}^{\theta_{n}(t)}f(s,x_{n}(\delta_{n}(t)))ds)\leq(L_{C}+h(x(\delta_{n}(t)))+\sqrt{\gamma})\mu_{n}.$

    2. (ii)

$\|x_{n}(\theta_{n}(t))-x_{0}\|\leq K_{1}.$

    3. (iii)

$\|x_{n}(t)\|\leq K_{2}.$

    4. (iv)

$\|x_{n}(\theta_{n}(t))-x_{n}(\delta_{n}(t))\|\leq
K_{3}\mu_{n}+\sqrt{\varepsilon_{n}}.$

    5. (v)

$\|x_{n}(t)-x_{n}(\theta_{n}(t))\|\leq K_{4}\mu_{n}+2\sqrt{\varepsilon_{n}}$.

  2. $(b)$

There exists $K_{5}>0$ such that for all $t\in[0,T]$ and $m,n\in\mathbb{N}$ we
have

$d_{C(\theta_{n}(t))}(x_{m}(t))\leq
K_{5}\mu_{m}+L_{C}\mu_{n}+2\sqrt{\varepsilon_{m}}.$

  3. $(c)$

There exists $K_{6}>0$ such that for all $n\in\mathbb{N}$ and almost all
$t\in[0,T]$, $\|\dot{x}_{n}(t)\|\leq K_{6}$.

  4. $(d)$

For all $n\in\mathbb{N}$ and $k\in\\{0,1,...,n-1\\}$, there is $v_{k+1}^{n}\in
C(t_{k+1}^{n})$ such that for all $t\in]t_{k}^{n},t_{k+1}^{n}[$:

$\dot{x}_{n}(t)\in-\frac{\lambda_{n}(t)}{\mu_{n}}\partial_{P}d_{C(\theta_{n}(t))}(v_{k+1}^{n})+f(t,x_{n}(\delta_{n}(t)))+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\mathbb{B},$
(8)

where
$\lambda_{n}(t)=4\sqrt{\varepsilon_{n}}+(L_{C}+h(x(\delta_{n}(t)))+\sqrt{\gamma})\mu_{n}$.
Moreover, $\|v_{k+1}^{n}-x_{n}(\theta_{n}(t))\|<2\sqrt{\varepsilon_{n}}$.

###### Proof

$(a)$: Set $\mu_{n}:=T/n$ and let $(\varepsilon_{n})$ be a sequence of
nonnegative numbers such that $\varepsilon_{n}/\mu_{n}^{2}\to 0$. We define
$\mathfrak{c}:=\sup_{n\in\mathbb{N}}\frac{\sqrt{\varepsilon_{n}}}{\mu_{n}}$.
We denote by $L_{h}$ the Lipschitz constant of $h$. For all $t\in[0,T]$ and
$n\in\mathbb{N}$, we define
$\tau_{n}(t):=x_{n}(\delta_{n}(t))+\int_{\delta_{n}(t)}^{\theta_{n}(t)}f(s,x_{n}(\delta_{n}(t)))ds$.
Since
$f(t,x_{n}(\delta_{n}(t)))\in\operatorname{proj}_{F(t,x_{n}(\delta_{n}(t)))}^{\gamma}(0)$
we obtain that

$\displaystyle d_{C(\theta_{n}(t))}(\tau_{n}(t))$ $\displaystyle\leq
d_{C(\theta_{n}(t))}(x_{n}(\delta_{n}(t)))+\left\|\int_{\delta_{n}(t)}^{\theta_{n}(t)}f(s,x_{n}(\delta_{n}(t)))ds\right\|$
$\displaystyle\leq
L_{C}\mu_{n}+\int_{\delta_{n}(t)}^{\theta_{n}(t)}\|f(s,x_{n}(\delta_{n}(t)))\|ds$
$\displaystyle\leq
L_{C}\mu_{n}+\int_{\delta_{n}(t)}^{\theta_{n}(t)}(h(x_{n}(\delta_{n}(t)))+\sqrt{\gamma})ds$
$\displaystyle\leq(L_{C}+h(x_{n}(\delta_{n}(t)))+\sqrt{\gamma})\mu_{n},$

which proves $(i)$. Moreover, since
$x_{n}(\theta_{n}(t))\in\operatorname{proj}_{C(\theta_{n}(t))}^{\varepsilon_{n}}(\tau_{n}(t))$,
we get that

$\displaystyle\|x_{n}(\theta_{n}(t))-\tau_{n}(t)\|$ $\displaystyle\leq
d_{C(\theta_{n}(t))}(\tau_{n}(t))+\sqrt{\varepsilon_{n}}$ (9)
$\displaystyle\leq(L_{C}+h(x_{n}(\delta_{n}(t)))+\sqrt{\gamma})\mu_{n}+\sqrt{\varepsilon_{n}},$

which yields

$\displaystyle\|x_{n}(\theta_{n}(t))-x_{n}(\delta_{n}(t))\|$
$\displaystyle\leq(L_{C}+2h(x_{n}(\delta_{n}(t)))+2\sqrt{\gamma})\mu_{n}+\sqrt{\varepsilon_{n}}$
(10)
$\displaystyle\leq(L_{C}+2h(x_{0})+2\sqrt{\gamma}+2L_{h}\|x_{n}(\delta_{n}(t))-x_{0}\|)\mu_{n}$
$\displaystyle+\sqrt{\varepsilon_{n}}.$

Hence, for all $t\in[0,T]$

$\displaystyle\|x_{n}(\theta_{n}(t))-x_{0}\|$
$\displaystyle\leq(1+2L_{h}\mu_{n})\|x_{n}(\delta_{n}(t))-x_{0}\|$
$\displaystyle+(L_{C}+2h(x_{0})+2\sqrt{\gamma})\mu_{n}+\sqrt{\varepsilon_{n}}.$

The above inequality means that for all $k\in\\{0,1,...,n-1\\}$:

$\displaystyle\|x_{k+1}^{n}-x_{0}\|$
$\displaystyle\leq(1+2L_{h}\mu_{n})\|x_{k}^{n}-x_{0}\|+(L_{C}+2h(x_{0})+2\sqrt{\gamma})\mu_{n}+\sqrt{\varepsilon_{n}}.$

Then, by (Clarke1998, , p. 183), we obtain that for all $k\in\\{0,...,n-1\\}$

$\displaystyle\|x_{k+1}^{n}-x_{0}\|$
$\displaystyle\leq(k+1)((L_{C}+2h(x_{0})+2\sqrt{\gamma})\mu_{n}+\sqrt{\varepsilon_{n}})\exp(2L_{h}(k+1)\mu_{n})$
(11) $\displaystyle\leq
T(L_{C}+2h(x_{0})+\sqrt{\gamma}+\mathfrak{c})\exp(2L_{h}T)=:K_{1}.$

which proves $(ii)$.
$(iii)$: By definition of $x_{n}$, for $t\in]t_{k}^{n},t_{k+1}^{n}]$ and
$k\in\\{0,1...,n-1\\}$, using (5)

$\displaystyle\|x_{n}(t)\|$
$\displaystyle\leq\|x_{k}^{n}\|+\|x_{k+1}^{n}-\tau_{n}(t)\|+\int_{t_{k}^{n}}^{t}\|f(s,x_{k}^{n})\|ds$
$\displaystyle\leq
K_{1}+\|x_{0}\|+(L_{C}+\sqrt{\gamma}+h(x_{k}^{n}))\mu_{n}+\sqrt{\varepsilon_{n}}+(h(x_{k}^{n})+\sqrt{\gamma})\mu_{n},$

where we have used (9). Moreover, it is clear that for $k\in\\{0,...,n\\}$

$\displaystyle h(x_{k}^{n})\leq h(x_{0})+L_{h}\|x_{k}^{n}-x_{0}\|\leq
h(x_{0})+L_{h}K_{1}.$

Therefore, for all $t\in[0,T]$

$\displaystyle\|x_{n}(t)\|$ $\displaystyle\leq
K_{1}+\|x_{0}\|+(L_{C}+2(h(x_{0})+L_{h}K_{1}+\sqrt{\gamma}))\mu_{n}+\sqrt{\varepsilon_{n}}$
$\displaystyle\leq
K_{1}+\|x_{0}\|+T(L_{C}+2(h(x_{0})+L_{h}K_{1}+\sqrt{\gamma})+\mathfrak{c})=:K_{2},$

which proves $(iii)$.
$(iv)$: From (10) and (11) it is easy to see that there exists $K_{3}>0$ such
that for all $n\in\mathbb{N}$ and $t\in[0,T]$:
$\|x_{n}(\theta_{n}(t))-x_{n}(\delta_{n}(t))\|\leq
K_{3}\mu_{n}+\sqrt{\varepsilon_{n}}$.
$(v)$: To conclude this part, we consider $t\in]t_{k}^{n},t_{k+1}^{n}]$ for
some $k\in\\{0,1,...,n-1\\}$. Then $x_{n}(\theta_{n}(t))=x_{k+1}^{n}$ and also

$\displaystyle\|x_{n}(\theta_{n}(t))-x_{n}(t)\|\leq$
$\displaystyle\|x_{k+1}^{n}-x_{k}^{n}\|+\|x_{k+1}^{n}-\tau_{n}(t)\|+\int_{t_{k}^{n}}^{t}\|f(s,x_{k}^{n})\|ds$
$\displaystyle\leq$ $\displaystyle
K_{3}\mu_{n}+\sqrt{\varepsilon_{n}}+(L_{C}+\sqrt{\gamma}+h(x_{0})+L_{h}K_{1})\mu_{n}+\sqrt{\varepsilon_{n}}$
$\displaystyle+\mu_{n}(h(x_{k}^{n})+\sqrt{\gamma})$ $\displaystyle\leq$
$\displaystyle(\underbrace{K_{3}+L_{C}+2(h(x_{0})+L_{h}K_{1})+2\sqrt{\gamma}}_{=:K_{4}})\mu_{n}+2\sqrt{\varepsilon_{n}},$

and we conclude this first part.
$(b)$: Let $m,n\in\mathbb{N}$ and $t\in[0,T]$, then

$\displaystyle d_{C(\theta_{n}(t))}(x_{m}(t))$ $\displaystyle\leq
d_{C(\theta_{n}(t))}(x_{m}(\theta_{m}(t)))+\|x_{m}(\theta_{m}(t))-x_{m}(t)\|$
$\displaystyle\leq
d_{H}(C(\theta_{n}(t)),C(\theta_{m}(t)))+K_{4}\mu_{m}+2\sqrt{\varepsilon_{m}}$
$\displaystyle\leq
L_{C}|\theta_{n}(t)-\theta_{m}(t)|+K_{4}\mu_{m}+2\sqrt{\varepsilon_{m}}$
$\displaystyle\leq
L_{C}(\mu_{n}+\mu_{m})+K_{4}\mu_{m}+2\sqrt{\varepsilon_{m}}$

where we have used (v). Hence, by setting $K_{5}:=K_{4}+L_{C}$ we prove (b).
$(c)$: Let $n\in\mathbb{N}$, $k\in\\{0,1,...,n-1\\}$ and
$t\in]t_{k}^{n},t_{k+1}^{n}]$. Then,

$\displaystyle\|\dot{x}_{n}(t)\|$
$\displaystyle=\left\|\frac{1}{\mu_{n}}\left(x_{k+1}^{n}-x_{k}^{n}-\int_{t_{k}^{n}}^{t_{k+1}^{n}}f(s,x_{k}^{n})ds\right)+f(t,x_{k}^{n})\right\|$
$\displaystyle\leq\frac{1}{\mu_{n}}\|x_{n}(\theta_{n}(t))-\tau_{n}(t)\|+\|f(t,x_{k}^{n})\|$
$\displaystyle\leq\frac{1}{\mu_{n}}((L_{C}+h(x_{k}^{n})+\sqrt{\gamma})\mu_{n}+\sqrt{\varepsilon_{n}})+h(x_{k}^{n})+\sqrt{\gamma}$
$\displaystyle\leq\frac{\sqrt{\varepsilon_{n}}}{\mu_{n}}+L_{C}+2(h(x_{0})+L_{h}K_{1}+\sqrt{\gamma})$
$\displaystyle\leq\mathfrak{c}+L_{C}+2(h(x_{0})+L_{h}K_{1}+\sqrt{\gamma})=:K_{6},$

which proves $(c)$.
$(d)$: Fix $k\in\\{0,1,...,n-1\\}$ and $t\in]t_{k}^{n},t_{k+1}^{n}[$. Then,
$x_{k+1}^{n}\in\operatorname{proj}_{C(t_{k+1}^{n})}^{\varepsilon_{n}}(\tau_{n}(t))$.
Hence, by Lemma 1, there exists $v_{k+1}^{n}\in C(t_{k+1}^{n})$ such that
$\|x_{k+1}-v_{k+1}^{n}\|<2\sqrt{\varepsilon_{n}}$ and

$\displaystyle\tau_{n}(t)-x_{k+1}^{n}\in\alpha_{n}(t)\partial_{P}d_{C(t_{k+1}^{n})}(v_{k+1}^{n})+3\sqrt{\varepsilon_{n}}\mathbb{B},\
\forall t\in]t_{k}^{n},t_{k+1}^{n}[,$

where
$\alpha_{n}(t)=4\sqrt{\varepsilon_{n}}+d_{C(\theta_{n}(t))}(\tau_{n}(t))$. By
virtue of $(i)$,

$\alpha_{n}(t)\leq
4\sqrt{\varepsilon_{n}}+(L_{C}+h(x(\delta_{n}(t)))+\sqrt{\gamma})\mu_{n}=:\lambda_{n}(t).$

Then, for all $t\in]t_{k}^{n},t_{k+1}^{n}[$

$-\mu_{n}(\dot{x}_{n}(t)-f(t,x_{k}^{n}))\in\lambda_{n}(t)\partial_{P}d_{C(t_{k+1}^{n})}(v_{k+1}^{n})+3\sqrt{\varepsilon_{n}}\mathbb{B},$

which implies that $t\in]t_{k}^{n},t_{k+1}^{n}[$

$\dot{x}_{n}(t)\in-\frac{\lambda_{n}(t)}{\mu_{n}}\partial_{P}d_{C(t_{k+1}^{n})}(v_{k+1}^{n})+f(t,x_{k}^{n})+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\mathbb{B}.$

##  4 Prox-Regular Case

In this section, we will study the algorithm under the assumption of uniform
prox-regularity of the moving set. The classic catching-up algorithm in this
framework was studied MR2159846 , where the existence of solutions for (4) was
established for a set-valued map $F$ taking values in a fixed compact set.

###### Theorem 4.1

Suppose, in addition to the assumptions of Theorem 3.1, that $C(t)$ is
$\rho$-uniformly prox-regular for all $t\in[0,T]$, and there exists a
nonnegative integrable function $k$ such that for all $t\in[0,T]$ and
$x,x^{\prime}\in\mathcal{H}$

$\langle y-y^{\prime},x-x^{\prime}\rangle\leq k(t)\|x-x^{\prime}\|^{2},\
\forall y\in F(t,x),\forall y^{\prime}\in F(t,x^{\prime}).$ (12)

Then, the sequence of functions $(x_{n})_{n}$ generated by the algorithm (5)
and (6) converges uniformly to an absolutely continuous function $x$, which is
the unique solution of (4).

###### Proof

Consider $m,n\in\mathbb{N}$ with $m\geq n$ big enough such that for all
$t\in[0,T]$, $d_{C(\theta_{n}(t))}(x_{m}(t))<\rho$, this can be guaranteed by
Theorem $(b)$. Then, for a.e. $t\in[0,T]$

$\frac{d}{dt}\left(\frac{1}{2}\|x_{n}(t)-x_{m}(t)\|^{2}\right)=\langle\dot{x}_{n}(t)-\dot{x}_{m}(t),x_{n}(t)-x_{m}(t)\rangle.$

Let $t\in[0,T]$ where the above equality holds. Let $k,j\in\\{0,1,...,n-1\\}$
such that $t\in]t_{k}^{n},t_{k+1}^{n}]$ and $t\in]t_{j}^{m},t_{j+1}^{m}]$. On
the one hand, we have that

$\displaystyle\langle\dot{x}_{n}(t)-\dot{x}_{m}(t),x_{n}(t)-x_{m}(t)\rangle=$
$\displaystyle\
\langle\dot{x}_{n}(t)-\dot{x}_{m}(t),x_{n}(t)-x_{k+1}^{n}\rangle$ (13)
$\displaystyle+\langle\dot{x}_{n}(t)-\dot{x}_{m}(t),x_{k+1}^{n}-v_{k+1}^{n}\rangle$
$\displaystyle+\langle\dot{x}_{n}(t)-\dot{x}_{m}(t),v_{k+1}^{n}-v_{j+1}^{m}\rangle$
$\displaystyle+\langle\dot{x}_{n}(t)-\dot{x}_{m}(t),v_{j+1}^{m}-x_{j+1}^{m}\rangle$
$\displaystyle+\langle\dot{x}_{n}(t)-\dot{x}_{m}(t),x_{j+1}^{m}-x_{m}(t)\rangle$
$\displaystyle\leq$ $\displaystyle\
2K_{6}(K_{4}(\mu_{n}+\mu_{m})+4(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}}))$
$\displaystyle+\langle\dot{x}_{n}(t)-\dot{x}_{m}(t),v_{k+1}^{n}-v_{j+1}^{m}\rangle$

where $v_{k+1}^{n}\in C(t_{k+1}^{n})$ and $v_{j+1}^{m}\in C(t_{j+1}^{m})$ are
the given in Theorem $(d)$. We can see that

$\displaystyle\max\\{d_{C(t_{k+1}^{n})}(v_{j+1}^{m}),d_{C(t_{j+1}^{m})}(v_{k+1}^{n})\\}$
$\displaystyle\leq d_{H}(C(t_{j+1}^{m}),C(t_{k+1}^{n}))$ $\displaystyle\leq
L_{C}|t_{j+1}^{m}-t_{k+1}^{n}|\leq L_{C}(\mu_{n}+\mu_{m}).$

From now, $m,n\in\mathbb{N}$ are big enough such that
$L_{C}(\mu_{n}+\mu_{m})<\frac{\rho}{2}$. Moreover, as $h$ is
$L_{h}$-Lipschitz, we have that for all $p\in\mathbb{N}$,
$i\in\\{0,1,...,p\\}$ and $t\in[0,T]$

$\|f(t,x_{i}^{p})\|\leq h(x_{i}^{p})+\sqrt{\gamma}\leq
h(x_{0})+L_{h}K_{1}+\sqrt{\gamma}=:\alpha.$

From the other hand, using (8) and Proposition 1 we have that

$\displaystyle\frac{1}{\digamma}\max\\{\left\langle\zeta_{n}-\dot{x}_{n}(t),v_{j+1}^{m}-v_{k+1}^{n}\right\rangle,\left\langle\zeta_{m}-\dot{x}_{m}(t),v_{k+1}^{n}-v_{j+1}^{m}\right\rangle\\}$
$\displaystyle\leq$
$\displaystyle\frac{2}{\rho}\|v_{k+1}^{n}-v_{j+1}^{m}\|^{2}+L_{C}(\mu_{n}+\mu_{m})$

where $\xi_{n},\xi_{m}\in\mathbb{B}$,
$\digamma:=\sup\\{\frac{\lambda_{\ell}(t)}{\mu_{\ell}}:t\in[0,T],\ell\in\mathbb{N}\\}$
and
$\zeta_{i}:=f(t,x_{i}(\delta_{i}(t)))+\frac{3\sqrt{\varepsilon_{i}}}{\mu_{i}}\xi_{i}$
for $i\in\\{n,m\\}$. Therefore, we have that

$\displaystyle\langle\dot{x}_{n}(t)-\dot{x}_{m}(t),v_{k+1}^{n}-v_{j+1}^{m}\rangle$
$\displaystyle\hskip
59.75078pt=\langle\dot{x}_{n}(t)-\zeta_{n},v_{k+1}^{n}-v_{j+1}^{m}\rangle+\langle\zeta_{n}-\zeta_{m},v_{k+1}^{n}-v_{j+1}^{m}\rangle$
$\displaystyle\hskip
59.75078pt+\langle\zeta_{m}-\dot{x}_{m}(t),v_{k+1}^{n}-v_{j+1}^{m}\rangle$
$\displaystyle\hskip 59.75078pt\leq
2\digamma(\frac{2}{\rho}\|v_{k+1}^{n}-v_{j+1}^{m}\|^{2}+L_{C}(\mu_{n}+\mu_{m}))+\langle\zeta_{n}-\zeta_{m},v_{k+1}^{n}-v_{j+1}^{m}\rangle$
$\displaystyle\hskip
59.75078pt\leq\frac{4\digamma}{\rho}(\|x_{n}(t)-x_{m}(t)\|+3(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}})+K_{4}(\mu_{n}+\mu_{m}))^{2}$
$\displaystyle\hskip 59.75078pt+2\digamma
L_{C}(\mu_{n}+\mu_{m})+\langle\zeta_{n}-\zeta_{m},v_{k+1}^{n}-v_{j+1}^{m}\rangle.$

Moreover, using Theorem 3.1 and property (12),

$\displaystyle\langle\zeta_{n}-\zeta_{m},v_{k+1}^{n}-v_{j+1}^{m}\rangle$
$\displaystyle\hskip 59.75078pt=\langle
f(t,x_{n}(\delta_{n}(t)))-f(t,x_{m}(\delta_{m}(t))),x_{n}(\delta_{n}(t))-x_{m}(\delta_{m}(t))\rangle$
$\displaystyle\hskip 59.75078pt+\langle
f(t,x_{n}(\delta_{n}(t)))-f(t,x_{m}(\delta_{m}(t))),v_{k+1}^{n}-x_{k+1}^{n}\rangle$
$\displaystyle\hskip 59.75078pt+\langle
f(t,x_{n}(\delta_{n}(t)))-f(t,x_{m}(\delta_{m}(t))),x_{k+1}^{n}-x_{k}^{n}\rangle$
$\displaystyle\hskip 59.75078pt+\langle
f(t,x_{n}(\delta_{n}(t)))-f(t,x_{m}(\delta_{m}(t))),x_{j}^{m}-x_{j+1}^{m}\rangle$
$\displaystyle\hskip 59.75078pt+\langle
f(t,x_{n}(\delta_{n}(t)))-f(t,x_{m}(\delta_{m}(t))),x_{j+1}^{m}-v_{j+1}^{m}\rangle$
$\displaystyle\hskip
59.75078pt+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\langle\xi_{n},v_{k+1}^{n}-v_{j+1}^{m}\rangle+\frac{3\sqrt{\varepsilon_{m}}}{\mu_{m}}\langle\xi_{m},v_{j+1}^{m}-v_{k+1}^{n}\rangle$
$\displaystyle\hskip 59.75078pt\leq\
k(t)\|x_{n}(\delta_{n}(t))-x_{m}(\delta_{m}(t))\|^{2}$ $\displaystyle\hskip
59.75078pt+2\alpha(3(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}})+K_{3}(\mu_{n}+\mu_{m}))$
$\displaystyle\hskip
59.75078pt+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\|v_{k+1}^{n}-v_{j+1}^{m}\|+\frac{3\sqrt{\varepsilon_{m}}}{\mu_{m}}\|v_{j+1}^{m}-v_{k+1}^{n}\|$
$\displaystyle\hskip 59.75078pt\leq
k(t)(\|x_{n}(t)-x_{m}(t)\|+3(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}})+(K_{3}+K_{4})(\mu_{n}+\mu_{m}))^{2}$
$\displaystyle\hskip
59.75078pt+2\alpha(3(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}})+K_{3}(\mu_{n}+\mu_{m}))$
$\displaystyle\hskip
59.75078pt+6\left(\frac{\sqrt{\varepsilon_{n}}}{\mu_{n}}+\frac{\sqrt{\varepsilon_{m}}}{\mu_{m}}\right)(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}}+K_{2}).$

These two inequalities and (13) yield

$\displaystyle\frac{d}{dt}\|x_{n}(t)-x_{m}(t)\|^{2}$ $\displaystyle\hskip
28.45274pt\leq\
4\left(\frac{4\digamma}{\rho}+k(t)\right)\|x_{n}(t)-x_{m}(t)\|^{2}+4\alpha(3(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}})+K_{3}(\mu_{n}+\mu_{m}))$
$\displaystyle\hskip 28.45274pt+4\digamma
L_{C}(\mu_{n}+\mu_{m})+12\left(\frac{\sqrt{\varepsilon_{n}}}{\mu_{n}}+\frac{\sqrt{\varepsilon_{m}}}{\mu_{m}}\right)(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}}+K_{2})$
$\displaystyle\hskip
28.45274pt+\frac{16\digamma}{\rho}(3(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}})+K_{4}(\mu_{n}+\mu_{m}))^{2}$
$\displaystyle\hskip
28.45274pt+4k(t)(3(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}})+(K_{3}+K_{4})(\mu_{n}+\mu_{m}))^{2}.$

Hence, using Gronwall’s inequality, we have for all $t\in[0,T]$ and $n,m$ big
enough:

$\|x_{n}(t)-x_{m}(t)\|^{2}\leq
A_{m,n}\exp\left(\frac{16\digamma}{\rho}T+4\int_{0}^{T}k(s)ds\right),$ (14)

where,

$\displaystyle A_{m,n}$ $\displaystyle=\ 4\alpha
T(3(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}})+K_{3}(\mu_{n}+\mu_{m}))$
$\displaystyle+4T\digamma
L_{C}(\mu_{n}+\mu_{m})+12T\left(\frac{\sqrt{\varepsilon_{n}}}{\mu_{n}}+\frac{\sqrt{\varepsilon_{m}}}{\mu_{m}}\right)(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}}+K_{2})$
$\displaystyle+\frac{16T\digamma}{\rho}(3(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}})+K_{4}(\mu_{n}+\mu_{m}))^{2}$
$\displaystyle+4\|k\|_{1}(3(\sqrt{\varepsilon_{n}}+\sqrt{\varepsilon_{m}})+(K_{3}+K_{4})(\mu_{n}+\mu_{m}))^{2}.$

Since $A_{m,n}$ goes to $0$ when $m,n\to\infty$, it shows that
$(x_{n})_{n\in\mathbb{N}}$ is a Cauchy sequence in the space of continuous
functions with the uniform convergence. Therefore, it converges uniformly to
some continuous function $x\colon[0,T]\to\mathcal{H}$. It remains to check
that $x$ is absolutely continuous, and it is the unique solution of (4). First
of all, by Theorem $(c)$ and (MR3626639, , Lemma 2.2), $x$ is absolutely
continuous and there is a subsequence of $(\dot{x}_{n})$ which converges
weakly in $L^{1}([0,T];\mathcal{H})$ to $\dot{x}$. So, without relabeling, we
have $\dot{x}_{n}\rightharpoonup\dot{x}$ in $L^{1}([0,T];\mathcal{H})$. On the
other hand, using Theorem $(d)$ and defining $v_{n}(t):=v_{k+1}^{n}$ for
$t\in]t_{k}^{n},t_{k+1}^{n}]$ we have

$\displaystyle\dot{x}_{n}(t)$
$\displaystyle\in-\frac{\lambda_{n}(t)}{\mu_{n}}\partial_{P}d_{C(\theta_{n}(t))}(v_{n}(t))+f(t,x_{n}(\delta_{n}(t)))+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\mathbb{B}$
$\displaystyle\in-\kappa_{1}\partial
d_{C(\theta_{n}(t))}(v_{n}(t))+\kappa_{2}\mathbb{B}\cap
F(t,x_{n}(\delta_{n}(t)))+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\mathbb{B}.$

where, by Theorem 3.1, $\kappa_{1}$ and $\kappa_{2}$ are nonnegative numbers
which do not depend of $n\in\mathbb{N}$ and $t\in[0,T]$. We also have
$v_{n}\to x$, $\theta_{n}\to\text{Id}_{[0,T]}$ and
$\delta_{n}\to\text{Id}_{[0,T]}$ uniformly. Theorem $(b)$ ensures that
$x(t)\in C(t)$ for all $t\in[0,T]$. By Mazur’s Lemma, there is a sequence
$(y_{j})$ such that for all $n$, $y_{n}\in\text{co}(\dot{x}_{k}:k\geq n)$ and
$(y_{n})$ converges strongly to $\dot{x}$ in $L^{1}([0,T];\mathcal{H})$. That
is to say

$y_{n}(t)\in\text{co}\left(-\kappa_{1}\partial
d_{C(\theta_{k}(t))}(v_{k}(t))+\kappa_{2}\mathbb{B}\cap
F(t,x_{k}(\delta_{k}(t)))+\frac{3\sqrt{\varepsilon_{k}}}{\mu_{k}}\mathbb{B}:k\geq
n\right).$

Hence, there exists $(y_{n_{j}})$ which converges to $\dot{x}$ almost
everywhere in $[0,T]$. Then, by virtue of Lemma 2, $(\mathcal{H}_{1}^{F})$ and
Lemma 3, we obtain that

$\dot{x}(t)\in-\kappa_{1}\partial d_{C(t)}(x(t))+\kappa_{2}\mathbb{B}\cap
F(t,x(t))\textrm{ for a.e. }t\in[0,T].$

Finally, since $\partial d_{C(t)}(x(t))\subset N(C(t);x(t))$ for all
$t\in[0,T]$, we have $x$ is the solution of (4).

###### Remark 1

The property required for $F$ in (12) is a classical monotonicity assumption
in the theory of existence of solutions for differential inclusions (see,
e.g., (MR1189795, , Theorem 10.5)).

######  Remark 2 (Rate of convergence)

In the precedent proof, we have established the following estimation:

$\|x_{n}(t)-x_{m}(t)\|^{2}\leq
A_{m,n}\exp\left(\frac{16\digamma}{\rho}T+4\int_{0}^{T}k(s)ds\right)$

for $m,n$ such that $\mu_{n}+\mu_{m}<\frac{\rho}{2L_{C}}$. Hence, by letting
$m\to\infty$, we obtain that

$\|x_{n}(t)-x(t)\|^{2}\leq
A_{n}\exp\left(\frac{16\digamma}{\rho}T+4\int_{0}^{T}k(s)ds\right)\textrm{ for
all }n>\frac{2L_{C}T}{\rho},$

where

$A_{n}:=\lim_{m\to\infty}A_{m,n}\leq
D(\sqrt{\varepsilon_{n}}+\mu_{n}+\frac{\sqrt{\varepsilon_{n}}}{\mu_{n}}),$

where $D$ is a nonnegative constant. Hence, the above estimation provides a
rate of convergence for our scheme.

##  5 Subsmooth case

In this section, we study the sweeping process (4) in a more general setting
than the uniformly prox-regular case. We now assume $(C(t))_{t\in[0,T]}$ is a
equi-uniformly subsmooth family. The classical catching-up algorithm was
studied in MR3574145  under this framework. In this case, we make the
assumption about the ball compactness of the moving sets. We will see that our
algorithm allows us to prove the existence of a solution, but we only ensure
that a subsequence converges to this solution, which is expected due to the
lack of uniqueness of solutions in this case.

###### Theorem 5.1

Suppose, in addition to assumptions of theorem 3.1, that the family
$(C(t))_{t\in[0,T]}$ is equi-uniformly subsmooth and the set $C(t)$ are ball-
compact for all $t\in[0,T]$. Then, the sequence of continuous functions
$(x_{n})_{n}$ generated by algorithm (5) and (6) converges uniformly (up to a
subsequence) to an absolutely continuous function $x$, which is a solution of
(4).

###### Proof

From Theorem $(d)$ we have for all $n\in\mathbb{N}$ and
$k\in\\{0,\ldots,n-1\\}$, there is $v_{k+1}^{n}\in C(t_{k+1}^{n})$ such that
$\|v_{k+1}^{n}-x_{k+1}^{n}\|<2\sqrt{\varepsilon_{n}}$ and for all
$t\in]t_{k}^{n},t_{k+1}^{n}]$:

$\dot{x}_{n}(t)\in-\frac{\lambda_{n}(t)}{\mu_{n}}\partial_{P}d_{C(\theta_{n}(t))}(v_{k+1}^{n})+f(t,x_{n}(\delta_{n}(t)))+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\mathbb{B}.$

where
$\lambda_{n}(t)=4\sqrt{\varepsilon_{n}}+(L_{C}+h(x(\delta_{n}(t)))+\sqrt{\gamma})\mu_{n}$.
As $h$ is $L_{h}$-Lipschitz it follows that

$\lambda_{n}(t)\leq(4\mathfrak{c}+L_{C}+h(x_{0})+\sqrt{\gamma}+L_{h}K_{1})\mu_{n}.$

Defining $v_{n}(t):=v_{k+1}^{n}$ on $]t_{k}^{n},t_{k+1}^{n}]$, then for all
$n\in\mathbb{N}$ and almost all $t\in[0,T]$

$\displaystyle\dot{x}_{n}(t)$
$\displaystyle\in-M\partial_{P}d_{C(\theta_{n}(t))}(v_{n}(t))+f(t,x_{n}(\delta_{n}(t)))+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\mathbb{B}$
(15) $\displaystyle\in-M\partial
d_{C(\theta_{n}(t))}(v_{n}(t))+M\mathbb{B}\cap
F(t,x_{n}(\delta_{n}(t)))+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\mathbb{B}.$

where $M:=4\mathfrak{c}+L_{C}+h(x_{0})+L_{h}K_{1}+\sqrt{\gamma}$. Moreover, by
Theorem $(b)$, we have

$\displaystyle d_{C(t)}(x_{n}(t))\leq
d_{C(\theta_{n}(t))}(x_{n}(t))+L_{C}\mu_{n}\leq(K_{5}+2L_{C})\mu_{n}+2\sqrt{\varepsilon_{n}}.$
(16)

for all $t\in[0,T]$.
Next, fix $t\in[0,T]$ and define $K(t):=\\{x_{n}(t):n\in\mathbb{N}\\}$. We
claim that $K(t)$ is relatively compact. Indeed, let $x_{m}(t)\in K(t)$ and
take $y_{m}(t)\in\operatorname{Proj}_{C(t)}(x_{m}(t))$ (the projection exists
due to the ball compactness of $C(t)$ and the boundedness of $K(t)$).
Moreover, according to (16) and Theorem $(a)$,

$\displaystyle\|y_{n}(t)\|$ $\displaystyle\leq
d_{C(t)}(x_{n}(t))+\|x_{n}(t)\|\leq(K_{5}+2L_{C})\mu_{n}+2\sqrt{\varepsilon_{n}}+K_{2}.$

This entails that $y_{n}(t)\in C(t)\cap R\,\mathbb{B}$ for all
$n\in\mathbb{N}$ for some $R>0$. Thus, by the ball compactness of $C(t)$,
there exists a subsequence $(y_{m_{k}}(t))_{m_{k}}$ of $(y_{m}(t))_{m}$
converging to some $y(t)$ as $k\to+\infty$. Then,

$\displaystyle\|x_{m_{k}}(t)-y(t)\|$ $\displaystyle\leq
d_{C(t)}(x_{m_{k}}(t))+\|y_{m_{k}}(t)-y(t)\|$
$\displaystyle\leq(K_{5}+2L_{C})\mu_{m_{k}}+2\sqrt{\varepsilon_{m_{k}}}+\|y_{m_{k}}(t)-y(t)\|,$

which implies that $K(t)$ is relatively compact. Moreover, it is not difficult
to see by Theorem $(c)$ that $K:=(x_{n})$ is equicontinuous. Therefore, by
virtue of Theorem $(c)$, Arzela-Ascoli’s and Dunford-Pettis’s Theorems, we
obtain the existence of a Lipschitz function $x$ and a subsequence
$(x_{j})_{j}$ of $(x_{n})_{n}$ such that

  1. (i)

$(x_{j})$ converges uniformly to $x$ on $[0,T]$.

  2. (ii)

$\dot{x}_{j}\rightharpoonup\dot{x}$ in $L^{1}\left([0,T];\mathcal{H}\right)$.

  3. (iii)

$x_{j}(\theta_{j}(t))\to x(t)$ for all $t\in[0,T]$.

  4. (iv)

$x_{j}(\delta_{j}(t))\to x(t)$ for all $t\in[0,T]$.

  5. (v)

$v_{j}(t)\to x(t)$ for all $t\in[0,T]$.

From (16) it is clear that $x(t)\in C(t)$ for all $t\in[0,T]$. By Mazur’s
Lemma, there is a sequence $(y_{j})$ such that for all $j$,
$y_{j}\in\text{co}(\dot{x}_{k}:k\geq j)$ and $(y_{j})$ converges strongly to
$\dot{x}$ in $L^{1}([0,T];\mathcal{H})$. That is to say

$y_{j}(t)\in\text{co}\left(-M\partial
d_{C(\theta_{n}(t))}(v_{n}(t))+M\mathbb{B}\cap
F(t,x_{n}(\delta_{n}(t)))+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\mathbb{B}:n\geq
j\right).$

On the other hand, there exists $(y_{n_{j}})$ which converges to $\dot{x}$
almost everywhere in $[0,T]$. Then, using Lemma 2, Lemma 3 and
$(\mathcal{H}_{1}^{F})$, we have

$\dot{x}(t)\in-M\partial d_{C(t)}(x(t))+M\mathbb{B}\cap F(t,x(t))\
\text{a.e.}$

Finally, since $\partial d_{C(t)}(x(t))\subset N(C(t);x(t))$ for all
$t\in[0,T]$, it follows that $x$ is the solution of (4).

##  6 Fixed set

In this section, we consider a closed and nonempty set $C\subset\mathcal{H}$,
and we look for a solution of the particular case of (4) given by

$\left\\{\begin{aligned}
\dot{x}(t)&\in-N\left(C;x(t)\right)+F(t,x(t))&\textrm{ a.e. }t\in[0,T],\\\
x(0)&=x_{0}\in C,\end{aligned}\right.$ (17)

where $F\colon[0,T]\times\mathcal{H}\rightrightarrows\mathcal{H}$ is a set-
valued map defined as above. The existence of a solution using classical
catching up was done in MR3956966 . Now we use similar ideas to get the
existence of a solution using our proposed algorithm. We emphasize that in
this case, no regularity of the set $C$ is required.

###### Theorem 6.1

Let $C\subset\mathcal{H}$ be a ball-compact set and
$F\colon[0,T]\times\mathcal{H}\rightrightarrows\mathcal{H}$ be a set-valued
map satisfying $(\mathcal{H}_{1}^{F})$, $(\mathcal{H}_{2}^{F})$ and
$(\mathcal{H}_{3}^{F})$. Then, for any $x_{0}\in S$, the sequence of functions
$(x_{n})_{n}$ generated by the algorithm (6) converges uniformly (up to a
subsequence) to a Lipschitz solution $x$ of the sweeping process (17) such
that

$\displaystyle\|\dot{x}(t)\|$ $\displaystyle\leq 2(h(x(t))+\sqrt{\gamma})$
$\displaystyle\textrm{ a.e. }t\in[0,T].$

###### Proof

We are going to use the properties of Theorem 3.1, where now we have
$L_{C}=0$. First of all, from Theorem $(d)$ we have for all $n\in\mathbb{N}$
and $k\in\\{0,1,\ldots,n-1\\}$, there is $v_{k+1}^{n}\in C$ such that
$\|v_{k+1}^{n}-x_{k+1}^{n}\|<2\sqrt{\varepsilon_{n}}$ and for all
$t\in]t_{k}^{n},t_{k+1}^{n}]$:

$\dot{x}_{n}(t)\in-\frac{\lambda_{n}(t)}{\mu_{n}}\partial_{P}d_{C}(v_{k+1}^{n})+f(t,x_{n}(\delta_{n}(t)))+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\mathbb{B}.$

where
$\lambda_{n}(t)=4\sqrt{\varepsilon_{n}}+(h(x(\delta_{n}(t)))+\sqrt{\gamma})\mu_{n}$.
Defining $v_{n}(t):=v_{k+1}^{n}$ on $]t_{k}^{n},t_{k+1}^{n}]$, we get that for
all $n\in\mathbb{N}$ and a.e. $t\in[0,T]$

$\displaystyle\dot{x}_{n}(t)$
$\displaystyle\in-\frac{\lambda_{n}(t)}{\mu_{n}}\partial_{P}d_{C}(v_{n}(t))+f(t,x_{n}(\delta_{n}(t)))+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\mathbb{B}$
$\displaystyle\in-\frac{\lambda_{n}(t)}{\mu_{n}}\partial
d_{C}(v_{n}(t))+(h(t,x_{n}(\delta_{n}(t)))+\sqrt{\gamma})\mathbb{B}\cap
F(t,x_{n}(\delta_{n}(t)))+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\mathbb{B}.$

Moreover, by Theorem $(b)$, we have

$d_{C}(x_{n}(t))\leq K_{5}\mu_{n}+2\sqrt{\varepsilon_{n}}\textrm{ for all
}t\in[0,T].$

Next, fix $t\in[0,T]$ and define $K(t):=\\{x_{n}(t):n\in\mathbb{N}\\}$. We
claim that $K(t)$ is relatively compact. Indeed, let $x_{m}(t)\in K(t)$ and
take $y_{m}(t)\in\operatorname{Proj}_{C}(x_{m}(t))$ (the projection exists due
to the ball compactness of $C$ and the boundedness of $K(t)$). Moreover,
according to the above inequality and Theorem $(a)$,

$\|y_{n}(t)\|\leq d_{C}(x_{n}(t))+\|x_{n}(t)\|\leq
K_{5}\mu_{n}+2\sqrt{\varepsilon_{n}}+K_{2},$

which entails that $y_{n}(t)\in C\cap R\,\mathbb{B}$ for all $n\in\mathbb{N}$
for some $R>0$. Thus, by the ball-compactness of $C$, there exists a
subsequence $(y_{m_{k}}(t))_{m_{k}}$ of $(y_{m}(t))_{m}$ converging to some
$y(t)$ as $k\to+\infty$. Then,

$\displaystyle\|x_{m_{k}}(t)-y(t)\|$ $\displaystyle\leq
d_{C}(x_{m_{k}}(t))+\|y_{m_{k}}(t)-y(t)\|$ $\displaystyle\leq
K_{5}\mu_{m_{k}}+2\sqrt{\varepsilon_{m_{k}}}+\|y_{m_{k}}(t)-y(t)\|,$

which implies that $K(t)$ is relatively compact. Moreover, it is not difficult
to see by Theorem $(c)$ that the set $K:=(x_{n})$ is equicontinuous.
Therefore, by virtue of Theorem $(c)$, Arzela-Ascoli’s and Dunford-Pettis’s
Theorems, we obtain the existence of a Lipschitz function $x$ and a
subsequence $(x_{j})_{j}$ of $(x_{n})_{n}$ such that

  1. (i)

$(x_{j})$ converges uniformly to $x$ on $[0,T]$.

  2. (ii)

$\dot{x}_{j}\rightharpoonup\dot{x}$ in $L^{1}\left([0,T];\mathcal{H}\right)$.

  3. (iii)

$x_{j}(\theta_{j}(t))\to x(t)$ for all $t\in[0,T]$.

  4. (iv)

$x_{j}(\delta_{j}(t))\to x(t)$ for all $t\in[0,T]$.

  5. (v)

$v_{j}(t)\to x(t)$ for all $t\in[0,T]$.

  6. (vi)

$x(t)\in C$ for all $t\in[0,T]$.

By Mazur’s Lemma, there is a sequence $(y_{j})$ such that for all $j$,
$y_{j}\in\text{co}(\dot{x}_{k}:k\geq j)$ and $(y_{j})$ converges strongly to
$\dot{x}$ in $L^{1}([0,T];\mathcal{H})$. i.e.,

$y_{j}(t)\in\text{co}\left(-\alpha_{n}\partial
d_{C}(v_{n}(t))+\beta_{n}\mathbb{B}\cap
F(t,x_{n}(\delta_{n}(t)))+\frac{3\sqrt{\varepsilon_{n}}}{\mu_{n}}\mathbb{B}:n\geq
j\right),$

where
$\alpha_{n}:=\frac{4\sqrt{\varepsilon_{n}}}{\mu_{n}}+h(t,x_{n}(\delta_{n}(t)))+\sqrt{\gamma}$
and
$\beta_{n}:=\frac{4\sqrt{\varepsilon_{n}}}{\mu_{n}}+h(t,x_{n}(\delta_{n}(t)))$.
On the other hand, there exists $(y_{n_{j}})$ which converges to $\dot{x}$
almost everywhere in $[0,T]$. Then, using Lemma 2, Lemma 3 and
$(\mathcal{H}_{1}^{F})$, we have

$\dot{x}(t)\in-(h(x(t))+\sqrt{\gamma})\partial
d_{C}(x(t))+(h(x(t))+\sqrt{\gamma})\mathbb{B}\cap F(t,x(t))\text{ for a.e.
}t\in[0,T].$

Finally, since $\partial d_{C}(x(t))\subset N(C;x(t))$ for all $t\in[0,T]$, we
obtain that $x$ is the solution of (17).

##  7 Numerical methods for approximate projections

As stated before, in most cases, finding an explicit formula for the
projection onto a closed set is not possible. Therefore, one must resort to
numerical algorithms to obtain approximate projections. There are several
papers discussing this issue for different notions of approximate projections
(see, e.g., pmlr-v139-usmanova21a ). These algorithms are called _projection
oracles_ and provide an approximate solution $\bar{z}\in\mathcal{H}$ to the
following optimization problem:

$\min_{z\in C}\|x-z\|^{2},$ ($P_{x}$)

where $C$ is a given closed set and $x\in\mathcal{H}$. Whether the approximate
solution $\bar{z}$ belongs to the set $C$ or not depends on the notion of
approximate projection. In our case, to implement our algorithm, we need that
$\bar{z}\in C$. In this line, a well-known projection oracle fulfilling this
property can be obtained via the celebrated Frank-Wolfe algorithm (see, e.g.,
frank1956algorithm ; pmlr-v28-jaggi13 ), where a linear sub-problem of
($P_{x}$) is solved in each iteration. For several types of convex sets, this
method has been successfully developed (see pmlr-v28-jaggi13 ; MR4275646 ;
MR4314104 ). Besides, in ding2018frank , it was shown that an approximate
solution of the linear sub-problem is enough to obtain a projection oracle.

Another important approach to obtaining approximate projections is the use of
the Frank-Wolfe algorithm with separation oracles (see dadush2022simple ).
Roughly speaking, a separation oracle determines whether a given point belongs
to a set and, in the negative case, provides a hyperplane separating the point
from the set (see MR0936633  for more details). For particular sets, it is
easy to get an explicit separation oracle (see (MR0936633, , p. 49)). An
important example is the case of a sublevel set: let
$g\colon\mathcal{H}\to\mathbb{R}$ be a continuous convex function and
$\lambda\in\mathbb{R}$. Then
$[g\leq\lambda]:=\\{x\in\mathcal{H}:g(x)\leq\lambda\\}$ has a separation
oracle described as follows: to verify that any point belongs to
$[g\leq\lambda]$ is straightforward. When a point $x\in\mathcal{H}$ does not
belong to $[g\leq\lambda]$, we can consider any $x^{\ast}\in\partial g(x)$.
Then, for all $y\in[g\leq\lambda]$,

$\langle x^{\ast},x\rangle\geq g(x)-g(y)+\langle x^{\ast},y\rangle>\langle
x^{\ast},y\rangle,$

where we have used that $g(x)>\lambda\geq g(y)$. Hence, the above inequality
shows the existence of the desired hyperplane, which provides a separation
oracle for $[g\leq\lambda]$. Therefore, if $C$ is the sublevel set of some
convex function, we can use the algorithm proposed in dadush2022simple  to get
an approximate solution $\bar{z}\in\operatorname{proj}_{S}^{\varepsilon}(x)$.
Moreover, the sublevel set enables us to consider the case

$C(t,x):=\bigcap_{i=1}^{m}\\{x\in\mathcal{H}:g_{i}(t,x)\leq
0\\}=\\{x\in\mathcal{H}:g(t,x):=\max_{i=1,...,m}g_{i}(t,x)\leq 0\\},$

where for all $t\in[0,T]$, $g_{i}(t,\cdot)\colon\mathcal{H}\to\mathbb{R}$,
$i=1,\ldots,m$ are convex functions. We refer to (MR4027814, , Proposition
5.1) for the proper assumptions on these functions to ensure the Lipschitz
property of the map $t\rightrightarrows C(t)$ holds (7).

##  8 Concluding remarks

In this paper, we have developed an enhanced version of the catching-up
algorithm for sweeping processes through an appropriate concept of approximate
projections. We provide the proposed algorithm’s convergence for three
frameworks: prox-regular, subsmooth, and merely closed sets. Some insights
into numerical procedures to obtain approximate projections were given mainly
in the convex case. Finally, the convergence of our algorithm for other
notions of approximate solutions will be explored in forthcoming works.

## References

  * (1) Acary, V., Bonnefon, O., Brogliato, B.: Nonsmooth modeling and simulation for switched circuits.  Springer (2011)
  * (2) Adly, S., Haddad, T.: Well-posedness of nonconvex degenerate sweeping process via unconstrained evolution problems.  Nonlinear Anal. Hybrid Syst. 36, 100,832, 13 (2020)
  * (3) Aliprantis, C.D., Border, K.C.: Infinite dimensional analysis, 3rd edn.  Springer, Berlin (2006)
  * (4) Aussel, D., Daniilidis, A., Thibault, L.: Subsmooth sets: functional characterizations and related concepts.  Trans. Amer. Math. Soc. 357(4), 1275–1301 (2005)
  * (5) Bomze, I.M., Rinaldi, F., Zeffiro, D.: Frank-Wolfe and friends: a journey into projection-free first-order optimization methods.  4OR 19(3), 313–345 (2021)
  * (6) Borwein, J.M., Preiss, D.: A smooth variational principle with applications to subdifferentiability and to differentiability of convex functions.  Trans. Amer. Math. Soc. 303(2), 517–527 (1987)
  * (7) Bounkhel, M., Thibault, L.: Nonconvex sweeping process and prox-regularity in Hilbert space.  J. Nonlinear Convex Anal. 6(2), 359–374 (2005)
  * (8) Brogliato, B.: Nonsmooth mechanics, 3rd edn.  Commun. Numer. Methods Eng. Springer (2016)
  * (9) Clarke, F.: Optimization and nonsmooth analysis.  Wiley Intersciences, New York (1983)
  * (10) Clarke, F., Ledyaev, Y., Stern, R., Wolenski, P.: Nonsmooth analysis and control theory, _Grad. Texts in Math._ , vol. 178.  Springer-Verlag, New York (1998)
  * (11) Colombo, G., Thibault, L.: Prox-regular sets and applications.  In: Handbook of nonconvex analysis and applications, pp. 99–182. Int. Press, Somerville, MA (2010)
  * (12) Combettes, C.W., Pokutta, S.: Complexity of linear minimization and projection on some sets.  Oper. Res. Lett. 49(4), 565–571 (2021)
  * (13) Dadush, D., Hojny, C., Huiberts, S., Weltge, S.: A simple method for convex optimization in the oracle model.  In: International Conference on Integer Programming and Combinatorial Optimization, pp. 154–167. Springer (2022)
  * (14) Deimling, K.: Multivalued differential equations, _De Gruyter Series in Nonlinear Analysis and Applications_ , vol. 1.  Walter de Gruyter & Co., Berlin (1992)
  * (15) Ding, L., Udell, M.: Frank-wolfe style algorithms for large scale optimization.  Large-Scale and Distributed Optimization pp. 215–245 (2018)
  * (16) Federer, H.: Curvature measures.  Trans. Amer. Math. Soc. 93, 418–491 (1959)
  * (17) Frank, M., Wolfe, P., et al.: An algorithm for quadratic programming.  Naval research logistics quarterly 3(1-2), 95–110 (1956)
  * (18) Grötschel, M., Lovász, L., Schrijver, A.: Geometric algorithms and combinatorial optimization, _Algorithms and Combinatorics: Study and Research Texts_ , vol. 2.  Springer-Verlag, Berlin (1988)
  * (19) Haddad, T., Noel, J., Thibault, L.: Perturbed sweeping process with a subsmooth set depending on the state.  Linear Nonlinear Anal. 2(1), 155–174 (2016)
  * (20) Hiriart-Urruty, J.B., López, M.A., Volle, M.: The $\varepsilon$-strategy in variational analysis: illustration with the closed convexification of a function.  Rev. Mat. Iberoam. 27(2), 449–474 (2011)
  * (21) Jaggi, M.: Revisiting Frank-Wolfe: Projection-free sparse convex optimization.  In: S. Dasgupta, D. McAllester (eds.) Proceedings of the 30th International Conference on Machine Learning, _Proceedings of Machine Learning Research_ , vol. 28(1), pp. 427–435. PMLR, Atlanta, Georgia, USA (2013)
  * (22) Jourani, A., Vilches, E.: Moreau-Yosida regularization of state-dependent sweeping processes with nonregular sets.  J. Optim. Theory Appl. 173(1), 91–116 (2017)
  * (23) Maury, B., Venel, J.: Un modéle de mouvement de foule.  ESAIM Proc. 18, 143–152 (2007)
  * (24) Moreau, J.J.: Rafle par un convexe variable I, expo. 15.  Sém, Anal. Conv. Mont. pp. 1–43 (1971)
  * (25) Moreau, J.J.: Rafle par un convexe variable II, expo. 3.  Sém, Anal. Conv. Mont. pp. 1–36 (1972)
  * (26) Papageorgiou, N.S., Kyritsi-Yiallourou, S.T.: Handbook of applied analysis, _Adv. Mech. Math._ , vol. 19.  Springer, New York (2009)
  * (27) Penot, J.P.: Calculus without derivatives, _Grad. Texts in Math_ , vol. 266\.  Springer, New York (2013)
  * (28) Poliquin, R.A., Rockafellar, R.T., Thibault, L.: Local differentiability of distance functions.  Trans. Amer. Math. Soc. 352(11), 5231–5249 (2000)
  * (29) Thibault, L.: Unilateral Variational Analysis in Banach Spaces.  World Scientific Publishing Co. Pte. Ltd. (2023).  Part II: Special Classes of Functions and Sets
  * (30) Usmanova, I., Kamgarpour, M., Krause, A., Levy, K.: Fast projection onto convex smooth constraints.  In: M. Meila, T. Zhang (eds.) Proceedings of the 38th International Conference on Machine Learning, _Proceedings of Machine Learning Research_ , vol. 139, pp. 10,476–10,486. PMLR (2021)
  * (31) Venel, J.: A numerical scheme for a class of sweeping processes.  Numer. Math. 118(2), 367–400 (2011)
  * (32) Vilches, E.: Existence and Lyapunov pairs for the perturbed sweeping process governed by a fixed set.  Set-Valued Var. Anal. 27(2), 569–583 (2019)