Functions of bounded Musielak-Orlicz-type deformation and anisotropic Total Generalized Variation
for image-denoising problems


Abstract

In the first part of this paper we introduce the space of bounded deformation fields with generalized Orlicz growth. We establish their main properties, provide a modular representation, and characterize a decomposition of the modular into an absolutely continuous part and a singular part weighted via a recession function. A further analysis in the variable exponent case is also provided. The second part of the paper contains a notion of Musielak-Orlicz anisotropic Total Generalized Variation. We establish a duality representation, and show well-posedness of the corresponding image reconstruction problem.

2020 Mathematics Subject Classification: 26A45, 26B30, 46E35, 46E99, 49J27, 49J45

Keywords and phrases: Functions of bounded deformation, generalized Orlicz spaces, modular, Total Generalized Variation, imaging.

1 Introduction↩︎

Ever since the early 90s, the variational approach has been a cornerstone of the mathematical image reconstruction. Generally speaking, in its simpler formulation, this methodology can be summarized as follows: for a given bounded domain \(\Omega \subset \mathbb{R}^2\), grayscale images are identified with functions \(u : \Omega \to \mathbb{R}\). Starting from a datum \(f\), possibly corrupted by noise, a clean image is obtained through the minimization of an energy functional of the form \[\mathcal{I}(u) := \alpha \|Ku - f\|_H + \beta \, \mathcal{R}(u),\] for \(\alpha,\beta > 0\). The structure of \(\mathcal{I}\) relies on the competition between the fidelity term \(\|Ku - f\|_H\), for a given normed space \(H\) and a forward operator \(K\), intended to keep track of the original datum \(f\), and a regularizing term \(\mathcal{R}(u)\), whose task is to filter out the noise and return a clean image. The strength of both effects is quantified by the two tuning parameters \(\alpha\) and \(\beta\).
In this paper we propose a novel regularizer based on a Musielak-Orlicz generalization of classical Total Generalized Variation functionals, and provide a first study of the associated space of functions with bounded deformation with generalized Orlicz growth.
Depending on specific applications, various types of variational regularizers have been proposed in the literature. A classical example is the ROF model, where \(\mathcal{R}(u)\) is the total variation of \(u\), cf. [1]. A known drawback of such functional is the so-called staircasing effect, namely the tendency of producing artificially piecewise-affine artifacts in the image denoising, see [2], [3]. In order to counteract such problems, several higher-order models have been introduced. A particularly successful one is the second-order Total Generalized Variation, defined as \[TGV^{2}_\alpha(u) := \sup \left\{ \int_\Omega u \, \mathrm{div}^2 \psi \, dx \ : \ \psi \in C^2_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym}), \|\psi\|_\infty \le \alpha_1, \|\mathrm{div}\, \psi\|_\infty \le \alpha_2 \right\},\] for every \(u \in L^1_{\rm loc}(\Omega)\), with \(\alpha = (\alpha_1, \alpha_2) \in \mathbb{R}^+ \times \mathbb{R}^+\). We refer to the seminal paper [4] for a general introduction and higher-order generalizations, and to [5] for an in-depth study of the second-order formulation.
In the last few years, due to an increasing interest in accurate and adaptive image reconstruction, an emerging trend has been considering anisotropic models, allowing to tailor the denoising effects to specific features of different parts of the images. First steps in this direction have already been taken, e.g., in [2], [6], [7]. In particular, in [8] the authors have proposed an anisotropic ROF-type functional with \(\mathcal{R}\) of the form \(\mathcal{R}(u) = \int_\Omega \phi(x, |Du|) \, dx\) with \[\phi(x, t) := \begin{dcases} \tfrac{1}{p(x)} \, t^{p(x)} & \textrm{ if } t \in [0,1], \\ t - 1 + \tfrac{1}{p(x)} & \textrm{ if } t > 1, \end{dcases}\] and with the variable exponent \(p : \Omega \to [c,2]\), with \(c > 1\).
In parallel, the study of variational problems in variable, Orlicz, and generalized Orlicz spaces has boomed in the past ten years. Among the vast literature, we single out the contributions [9][14] and the references therein.
The extension of the model in [8] to incorporate the case \(c = 1\), as well as an Orlicz generalization for regularizers \(\mathcal{R}(u) \simeq \int_\Omega \varphi(x,|Du|)\) for \(\Phi\)-functions \(\varphi\) possibly having linear growth, is the subject of [15], where a thorough mathematical foundation of the space of anisotropic functions with bounded variation with generalized Orlicz growth has also been provided. A further step in this direction has been taken in [16], where the associated Euler-Lagrange equations and gradient flows have also been analyzed (see also [17]). The effect of an anisotropic norm in the Chan-Vese (or anisotropic Mumford-Shah) is the subject of [18]. Anisotropic models, at least in the ROF setting, have indeed proven to allow for more accurate computations, as well as for a staircasing reduction and a sharper texture recognition. We refer to [19][21] for an overview and further references.
Our contribution is twofold. First, we introduce a novel anisotropic functional setting of functions with bounded deformation, namely the space \(BD^\varphi\) of bounded deformation fields with generalized Orlicz growth. After analyzing its basic properties, we establish a link to a modular representation and provide a further study in the special case of variable Lebesgue integrability.
Second, we analyze a Musielak-Orlicz anisotropic Total Generalized Variation class of regularizers. Let \(\Omega\) be a bounded, connected and open set in \(\mathbb{R}^n\), with \(n \ge 2\). Inspired by [4], for a fixed generalized weak \(\Phi-\)function \(\varphi\in \Phi_w(\Omega)\) (see Section 2.3 for the detailed assumptions), we define our anisotropic Total Generalized Variation of order \(2\) with weight \(\alpha = (\alpha_1, \alpha_2) \in \mathbb{R}^+ \times \mathbb{R}^+\) as \[\label{TGV} \begin{align} TGV^{\varphi,2}_\alpha(u):= \sup\bigg\{ &\int_\Omega u \,\mathrm{div}^2 \psi \, dx \ : \ \\ & \;\psi \in C^2_c (\Omega; \mathbb{R}^{n \times n}_{\rm sym}), \|\psi\|_{\varphi^*} \le \alpha_1, \|\mathrm{div}\, \psi\|_{\varphi^*} \le \alpha_2 \bigg\}, \end{align}\tag{1}\] for every \(u \in L^1_{\textrm{loc}}(\Omega)\), where \(\|\cdot\|_{\varphi^*}\) is the norm associated to the Fenchel conjugate function of \(\varphi\), cf. Definition 10, and where the divergence of order \(2\) of \(\xi \in \mathbb{R}^{n \times n}_{\rm sym}\) is defined as \[\mathrm{div}^2 \xi := \sum_{i = 1}^n \frac{\partial^2 \xi_{ii}}{\partial x_i^2} + 2 \sum_{i<j} \frac{\partial^2 \xi_{ij}}{\partial x_i \partial x_j},\] (see also [4] for higher-order notions of divergence). Note that, in principle, \(TGV^{\varphi,2}(u)\) might take infinite values.

A prototypical example of admissible functions \(\varphi\) is given by \(\varphi(t) := \tfrac{1}{p(x)} t^{p(x)}\) (see also [15]), with \(p : \Omega \to [1,+\infty]\), thus possibly describing linear, superlinear, or even infinite growth in different portions of the domain.
After studying the basic properties of \(TGV^{\varphi,2}_\alpha\), we identify a dual representation, and show how classical existence and stability for classical Total Generalized Variation carry over to our Musielak-Orlicz setting.
The paper is organized as follows. In Section 2 we collect some definitions and technical results on \(\Phi\)-functions and generalized Orlicz spaces. Section 3 contains the definition and main properties of the novel space of bounded deformation fields with generalized Orlicz growth \(BD^\varphi\). In particular, in Subsection 3.1 we establish a connection between \(BD^{\varphi}\) and an associated weighted-modular representation, while in Subsection 3.2 we further establish a decomposition of the modular function in terms of the modular of the absolutely continuous part of the symmetric gradient and a singular part weighted by the recession function, cf. Theorem 10. Eventually, in Subsection 3.3 we specify our setting to the variable exponent case. The focus of Section 4 is on our anisotropic Total Generalized Variation. In Subsection 4.1 we establish its main properties, and in Subsection 4.2 we identify its dual reformulation in terms of a notion of anisotropic measure variation. Finally, in Subsection 4.3 we show stability with respect to the data and well-posedness of the minimization problem. In Appendices 5 and 6 we give the proofs of some auxiliary results we need for the proof of Theorem 10.

2 Preliminary results↩︎

In this section, after giving some notation in Subsection 2.1, we collect some definitions and technical results on modulars and norms in Subsection 2.2 and on generalized \(\Phi\)-functions and Orlicz spaces in Subsection 2.3.

2.1 Notation↩︎

The notation \(f \approx g\) means that there exist two constants \(c_1, c_2 > 0\) such that \(c_1 f \le g \le c_2 f\). Given an exponent \(p \in (1,+\infty)\), we denote by \(p'\) its Hölder conjugate exponent, namely \(p':=\frac{p}{p-1}\). If \(p = 1\), then its Hölder conjugate exponent is \(p' = +\infty\) and viceversa. In the sequel \(\Omega\) denotes an open subset of \(\mathbb{R}^n\). We denote by \(L^0(\Omega)\) the set of measurable functions \(f : \Omega \to \mathbb{R}\) and we use the notation \(\mathcal{C}^k_c(\Omega; \mathbb{R}^n) := \{ \sigma \in \mathcal{C}^k(\Omega; \mathbb{R}^n) \ : \ {\rm supp}(\sigma) \textrm{ is compact in } \Omega \}\) with \(k \in \mathbb{N}\cup \{+\infty\}\). We denote with \(|\cdot|\) the Frobenius norm for matrices, i.e. \(|M|^2 := {\rm tr}(M M^\textrm{T})\) where \(M \in \mathbb{R}^{n \times n}\), and with \(:\) the scalar product associated with the Frobenius norm, i.e. \(A : B := {\rm tr}(A B^\textrm{T})\) for \(A,B \in \mathbb{R}^{n \times n}\). With the same notation, namely \(|\cdot|\), we also indicate the standard vector Euclidean norm, and with \(\cdot\) the scalar product in \(\mathbb{R}^n\).
We say that a function \(\omega: [0, +\infty) \to [0, +\infty]\) is a modulus of continuity if it is increasing and such that \(\omega(0) = \lim_{t \to 0^+} \omega(t)=0\). Note that we do not require concavity and we allow extended real values.
We denote by \(\mathbb{R}^{n \times n}_{\rm sym}\) the space of real-valued \(n \times n\) symmetric matrices. By \(\mathcal{M}(\Omega;X)\) we denote the space of all finite Radon measures in \(\Omega\) with values in \(X\), which can be \(\mathbb{R}\), \(\mathbb{R}^n\) or \(\mathbb{R}^{n \times n}_{\rm sym}\) depending on the situation. Given a set \(O \subset X\), we denote by \(I_O\) the indicator function of the set \(O\), namely \[I_O(x) := \begin{dcases} 0 & \textrm{ if } x \in O, \\ + \infty & \textrm{ if } x \notin O, \end{dcases}\] and with \(\chi_O\) the characteristic function of the set \(O\), namely \[\chi_O(x) := \begin{dcases} 1 & \textrm{ if } x \in O, \\ 0 & \textrm{ if } x \notin O. \end{dcases}\] Given a function \(u : \Omega \to \mathbb{R}^m\), we denote by \(Du\) its distributional gradient. We indicate by \(BV(\Omega;\mathbb{R}^m)\) the space of functions \(u : \Omega \to \mathbb{R}^m\) such that \(Du \in \mathcal{M}(\Omega;\mathbb{R}^m)\). If \(m=1\) we denote the latter space simply by \(BV(\Omega)\). When \(n=m\), for any function \(u : \Omega \to \mathbb{R}^n\) we use the notation \(E u\) for the distributional symmetric gradient of \(u\), namely \(E u = \frac{1}{2} (Du + Du ^ \textrm{T})\), where \(Du^\textrm{T}\) stands for the transpose of \(Du\). We write \[BD(\Omega) := \{u \in L^1(\Omega;\mathbb{R}^n) \ : \ E u \in \mathcal{M}(\Omega;\mathbb{R}^{n \times n}_{\rm sym}) \}.\] We will often indicate by \(c\) a generic constant, whose value may change from line to line.

2.2 Modulars and norms↩︎

In order to define the spaces we use in our work, we need to give different notions of norms and modulars that generate our spaces. Note that our terminology differs from the original one in [22], but follows the lines of [15], [23] in order to get a clearer connection between the modulars and their associated norms.

Definition 1. Let \(X\) be a real vector space and let \(\|\cdot\|: X \to [0, +\infty]\) be a function on \(X\). Consider the following conditions:

  • \(\|f\| = 0\) if and only if \(f = 0\),

  • \(\|af\| = |a| \, \|f\|\) for all \(f \in X\) and \(a \in \mathbb{R}\),

  • \(\|f + g\| \le \|f\| + \|g\|\) for all \(f, g \in X\),

  • there exists \(c > 0\) such that \(\|f + g\| \le c \, (\|f\| + \|g\|)\) for all \(f,g \in X\).

We say that \(\|\cdot\| : X \to [0, +\infty]\) is

  • a quasi-seminorm if and only if it satisfies (N2) and (N3’),

  • a seminorm if and only if it satisfies (N2) and (N3),

  • a quasinorm if and only if it satisfies (N1), (N2) and (N3’),

  • a norm if and only if it satisfies (N1), (N2) and (N3).

Note that conditions (N2) and (N3), equivalently (N2) and (N3’), already imply that if \(f = 0_X\) then \(\|f\| = 0\), namely condition (N1) only adds one of the two implications.

Definition 2. Let \(X\) be a real vector space. A function \(\varrho: X \to [0,+\infty]\) is called a quasi-semimodular on \(X\) if and only if

  • \(\varrho(0_X) = 0\),

  • the function \(\lambda \mapsto \varrho(\lambda x)\) is non-decreasing on \([0,+\infty)\) for every \(x \in X\),

  • \(\varrho(-x) = \varrho(x)\) for every \(x \in X\),

  • there exists \(\beta \in (0,1]\) such that for every \(x,y \in X\) and every \(\theta \in [0,1]\) \[\varrho\left(\beta(\theta x + (1-\theta) y) \right) \le \theta \varrho(x) + (1-\theta) \varrho(y).\]

If (iv) holds with \(\beta = 1\), then \(\varrho\) is called semimodular. A quasi-semimodular is called a quasimodular when \(\varrho(x) = 0\) if and only if \(x = 0_X\). A modular is defined analogously from a semimodular.

Definition 3. Let \(X\) be a real vector space. If \(\varrho\) is a quasi-semimodular in \(X\), then the modular space \[X_\varrho:= \{x \in X \ : \ \|x\|_\varrho< +\infty\}\] is defined by the quasi-seminorm \[\|x\|_\varrho:= \inf \left\{\lambda > 0 \ : \ \varrho\left(\frac{x}{\lambda}\right) \le 1 \right\}.\]

2.3 Generalized \(\Phi\)-functions and Orlicz spaces↩︎

In this subsection we give some definitions regarding \(\Phi\)-functions and we then define generalized Orlicz (or Musielak-Orlicz) spaces, following the lines of [11]. Since we are also going to make use of singular measures, the assumptions always need to get adjusted to hold for every point, as in [23].

Definition 4. A function \(f : \Omega \to \mathbb{R}\) is called \(L\)-almost increasing if there exists \(L \ge 1\) such that \(f(s) \le L f(t)\) for every \(s \le t\). If this inequality holds for \(L = 1\), the function is called non-decreasing. We define analogously \(L\)-almost decreasing and non-increasing functions.

Definition 5. We say that a function \(\varphi: \Omega \times [0, +\infty) \to [0, +\infty]\) is a weak \(\Phi\)-function, and we write \(\varphi\in \Phi_w(\Omega)\), if the following conditions hold:

  • the function \(x \mapsto \varphi(x,|f(x)|)\) is measurable for every function \(f \in L^0(\Omega)\),

  • the function \(t \mapsto \varphi(x,t)\) is non-decreasing for every \(x \in \Omega\),

  • \(\varphi(x,0) = \underset{t \to 0^+}{\lim} \varphi(x,t) = 0\) and \(\underset{t \to +\infty}{\lim} \varphi(x,t) = +\infty\) for every \(x \in \Omega\),

  • the function \(t \mapsto \tfrac{\varphi(x, t)}{t}\) is \(L\)-almost increasing on \((0,+\infty)\) for every \(x \in \Omega\) with a constant \(L\) independent of \(x\).

If, moreover, \(\varphi\in \Phi_w(\Omega)\) is convex and left-continuous with respect to \(t\) for every \(x \in \Omega\), then \(\varphi\) is a convex \(\Phi\)-function, and we write \(\varphi\in \Phi_c(\Omega)\). If \(\varphi\) does not depend on \(x\), then we omit the set and write \(\varphi\in \Phi_w\) or \(\varphi\in \Phi_c\).

Note that since we deal with conjugates of functions possibly with linear growth at infinity, it is needed to allow for extended real-valued \(\Phi\)-functions. Now we are in the position of defining generalized Orlicz spaces.

Definition 6. Let \(\varphi\in \Phi_w(\Omega)\) and define \[\varrho_\varphi(f) := \int_\Omega \varphi(x,|f|) \, dx \qquad \textrm{for every } f \in L^0(\Omega).\] Then, the set \[L^\varphi(\Omega) := (L^0(\Omega))_{\varrho_\varphi} = \left\{ f \in L^0(\Omega) \ : \ \varrho_\varphi(\lambda f) < +\infty \quad \textrm{ for some } \lambda > 0 \right\}\] is called a generalized Orlicz space* with quasinorm given by \(\|f\|_\varphi:= \|f\|_{\varrho_\varphi}\). We use the abbreviation \(\|f\|_\varphi:= \left\| |f| \right\|_\varphi\) for vector-valued functions.*

Thanks to [11], it is known that \(\|\cdot\|_\varphi\) is a quasinorm in \(L^\varphi(\Omega)\) if \(\varphi\in \Phi_w(\Omega)\), and a norm if \(\varphi\in \Phi_c(\Omega)\).
Now we need to define some standard conditions on \(\Phi\)-functions that grant them additional properties, which sometimes are inherited by their conjugate functions as well.

Definition 7. Let \(\varphi: \Omega \times [0,+\infty) \to [0,+\infty]\) and let \(p,q > 0\). Then we define the following conditions:

  • There exists \(\beta \in (0,1]\) such that \(\varphi(x,\beta) \le 1 \le \varphi(x, \tfrac{1}{\beta})\) for every \(x \in \Omega\),

  • For every \(K > 0\) there exists \(\beta \in (0,1]\) such that for every \(x,y \in \Omega\) \[\varphi(x,\beta t) \le \varphi(y,t) + 1 \quad \textrm{when } \varphi(y, t) \in \left[0, \tfrac{K}{|x-y|^n}\right],\]

  • For every \(K > 0\) there exists a modulus of continuity \(\omega\) such that for every \(x,y \in \Omega\) \[\varphi\left(x,\frac{t}{1 + \omega(|x-y|)}\right) \le \varphi(y,t) + \omega(|x-y|) \quad \textrm{when } \varphi(y, t) \in \left[0, \frac{K}{|x-y|^n}\right],\]

  • There exists \(L_p \ge 1\) such that \(t \mapsto \tfrac{\varphi(x,t)}{t^p}\) is \(L_p\)-almost increasing on \((0,+\infty)\) for every \(x \in \Omega\),

  • There exists \(L_q \ge 1\) such that \(t \mapsto \tfrac{\varphi(x,t)}{t^q}\) is \(L_q\)-almost decreasing on \((0,+\infty)\) for every \(x \in \Omega\).

We say that (aInc), respectively (aDec), holds if (aInc)\(_p\), respectively (aDec)\(_q\), holds for some \(p > 1\), respectively \(q > 1\).

Remark 1. Note that, by Definition 5, every function \(\varphi\in \Phi_w(\Omega)\) satisfies (aInc)\(_1\).

Note that if \(\varphi\in \Phi_w(\Omega)\), then \(\varphi(\, \cdot \, ,1) \approx 1\) implies (A0). Moreover, if \(\varphi\) satisfies (aDec), then (A0) and \(\varphi(\,\cdot\, ,1) \approx 1\) are equivalent. Assumption (A1) is an almost continuity condition and (aInc) and (aDec) are quantitative versions of the \(\nabla_2\) and \(\Delta_2\) conditions and measure lower and upper growth rates. We define the condition (A1) as in [15]. We point out that it slightly differs from [11], [24], where it is assumed that \[\varphi(x,\beta t) \le \varphi(y,t) \quad \textrm{when } \varphi(y, t) \in \left[1, \tfrac{1}{|B|}\right]\] and \(x\) and \(y\) belong to the ball \(B\). Lastly, the “vanishing (A1)” condition (VA1) is a continuity condition for \(\varphi\), introduced in [25] to prove maximal regularity of minimizers.
In order to state the last property which we need for our analysis, we first need to define the recession function of \(\varphi\).

Definition 8. Let \(\varphi\in \Phi_w(\Omega)\). We call recession function of \(\varphi\) the function \(\varphi^\infty : \Omega \to [0,+\infty]\) defined by \[\varphi^\infty(x) := \limsup_{t \to +\infty} \frac{\varphi(x,t)}{t} \qquad \textrm{for all } x \in \Omega.\]

Note that, since the continuity of \(\varphi\) is not necessarily inherited by \(\varphi^\infty\), this makes the non-autonomous case not trivial. For our analysis, we also need the following weaker version of (VA1) where at least one of the points has to belong to the set \(\{\varphi^\infty < +\infty\}\).

Definition 9. Let \(\varphi\in \Phi_w(\Omega)\). We say that \(\varphi\) satisfies restricted (VA1) if and only if it satisfies (A1) and for every \(K > 0\) there exists \(\omega\) modulus of continuity such that \[\varphi\left(x, \frac{t}{1 + \omega(|x-y|)}\right) \le \varphi(y,t) + \omega(|x-y|) \quad \textrm{when } \varphi(y, t) \in \left[0, \frac{K}{|x-y|^n}\right],\] for every \(x,y \in \Omega\) such that \(\varphi^\infty(x) < +\infty\) or \(\varphi^\infty(y) < +\infty\).

We now need to define the associate space of the generalized Orlicz space \(L^\varphi(\Omega)\). In general, the associate space is a variant of the dual function space which works better at the end-points \(p = 1\) and \(p = \infty\).

Definition 10. Let \(\varphi\in \Phi_w(\Omega)\). We define the associate space \((L^\varphi)'(\Omega)=\{u \in L^0(\Omega): \|u\|_{(L^\varphi)'}<+\infty\}\) with \[\|u\|_{(L^\varphi)'} := \sup_{\|v\|_\varphi\le 1} \int_\Omega uv \, dx.\]

Thanks to [11], it is well know that for weak \(\Phi\)-functions the associate space corresponds to the generalized Orlicz space \(L^{\varphi^*}(\Omega)\), where \[\varphi^*(x, t) := \sup_{s \ge 0} [st - \varphi(x, s)]\] is the (Fenchel) conjugate function of \(\varphi\in \Phi_w(\Omega)\). Moreover, it can be shown that (see [15]) for all \(u \in L^0(\Omega)\), \(\|u\|_{L^{\varphi}}\) is equivalent to \(\sup_{\|v\|_{L^{\varphi^*}} \le 1} \int_\Omega uv \, dx.\)

Let us state some well-known properties of the conjugate function.

Lemma 1. Let \(\varphi\in \Phi_w(\Omega)\). Then,

  • \(\varphi^* \in \Phi_c(\Omega)\), meaning that \(\varphi^*\) is always convex and left-continuous,

  • If \(p,q \in (1,+\infty)\), then \(\varphi\) satisfies (aInc\(_p\)) or (aDec\(_q\)) if and only if \(\varphi^*\) satisfies respectively (aDec\(_{p'}\)) or (aInc\(_{q'}\)),

  • If \(\varphi\in \Phi_c(\Omega)\), then \[\varphi^*\left(x, \tfrac{\varphi(x, t)}{t}\right) \le \varphi(x, t) \quad \textrm{for every } x \in \Omega \textrm{ and } t \in \mathbb{R}^+,\] and \(\varphi^{**} = \varphi\),

  • If \(\varphi\) satisfies (A0) or (A1), then so does \(\varphi^*\).

0◻

Conditions (i), (ii) and (iv) have been proven in [11], respectively in Lemma 2.4.1, Proposition 2.4.9 and Lemmas 3.7.6 and 4.1.7, while condition (iii) has been shown in [10], see also [11].
Lastly, we report the following result regarding some immersions of generalized Orlicz spaces into Lebesgue spaces.

Lemma 2. Let \(\Omega \subset \mathbb{R}^n\) be such that \(|\Omega| < +\infty\) and let \(\varphi\in \Phi_w(\Omega)\). If \(\varphi\) satisfies (A0), then \[L^\infty(\Omega) \hookrightarrow L^\varphi(\Omega) \hookrightarrow L^1(\Omega).\] Moreover, if \(\varphi\) satisfies (aInc)\(_p\), then \(L^\varphi(\Omega) \hookrightarrow L^p(\Omega)\). 0◻

Note that the second continuous inclusion comes from the fact that every weak \(\Phi\)-function satisfies (aInc)\(_1\), see Remark 1. This Lemma was shown in [11]. We point out that condition (A0) is essential both for establishing the above embeddings and to ensure that \(L^\varphi(\Omega)\) is a Banach space, cf. [11]. In particular, of (A0) is violated, the first inclusion in Lemma 2 might fail to hold. We refer to [11] for a counterexample.

3 The space of bounded deformations fields with generalized Orlicz growth↩︎

In [15], given a function \(\varphi\in \Phi_w(\Omega)\), an anisotropic total variation (there called dual norm) is introduced exploiting the conjugate function of \(\varphi\), namely for any \(u \in L^1(\Omega)\) we have \[\label{def:EHH95TV} V_\varphi(u) := \sup \left\{ \int_\Omega u \, \mathrm{div}\, w \, dx \ : \ w \in C^1_c(\Omega; \mathbb{R}^n), \|w\|_{\varphi^*} \le 1 \right\}.\tag{2}\] Moreover, the space of functions of bounded variation with generalized Orlicz growth is defined as \[\label{def:BV95phi} BV^\varphi(\Omega) := \{ u \in L^\varphi(\Omega) \ : \ \|u\|_{BV^\varphi(\Omega)} < +\infty \},\tag{3}\] where \[\label{norm95BVphi} \|u\|_{BV^\varphi(\Omega)} := \|u\|_\varphi+ V_\varphi(u).\tag{4}\] Equivalently, given \(\varphi\in \Phi_w(\Omega)\) in [26] a vectorial anisotropic total variation is introduced, namely for any \(m \in \mathbb{N}^+\) and \(u \in L^1_{\rm loc}(\Omega;\mathbb{R}^m)\) we have \[\label{var95BHHZ} V^m_\varphi(u) := \sup \left \{\sum_{i = 1}^m \int_\Omega u_i \, \mathrm{div}\, w_i \, dx \ : \ w \in \left[C^1_c(\Omega; \mathbb{R}^n)\right]^m, \|w\|_{\varphi^*} \le 1 \right\}.\tag{5}\] One can equivalently write \[V^m_\varphi(u) = \sup \left \{ \int_\Omega u \cdot \mathrm{div}\, w \, dx \ : \ w \in C^1_c(\Omega; \mathbb{R}^{m \times n}), \|w\|_{\varphi^*} \le 1 \right\},\] where the divergence \(\mathrm{div}\, w\) is a vector in \(\mathbb{R}^m\), namely taken a matrix-valued field \(\xi\), it is computed row-wise, i.e. \[\label{divrow} (\mathrm{div}\, \xi)_i := \mathrm{div}\, \xi_i = \sum_{j = 1}^n \frac{\partial \xi_{ij}}{\partial x_j},\tag{6}\] for any \(i \in [1,m] \cap \mathbb{N}\) and any \(\xi : \Omega \to \mathbb{R}^{m \times n}\). Moreover, we say that \(u\) belongs to \(BV^\varphi(\Omega;\mathbb{R}^m)\) if \[\|u\|_{BV^\varphi(\Omega;\mathbb{R}^m)} := \|u\|_\varphi+ V_\varphi^m(u) < +\infty.\] In the case \(m = 1\), 5 reduces to 2 and we denote \(BV^\varphi(\Omega;\mathbb{R})\) simply by \(BV^\varphi(\Omega)\).
Note that in general Harjulehto & Hästö in their works, see e.g. [11], use the symbol \(C^1_0(\Omega; \mathbb{R}^n)\) to denote the space \(C^1_c(\Omega; \mathbb{R}^n)\). For this reason, we gave the definition of the total variation with our notation instead of the one used in [15].

Remark 2. Clearly, when \(\varphi(x,t) = |t|\), \(BV^\varphi(\Omega;\mathbb{R}^m)= BV(\Omega;\mathbb{R}^m)\). A priori, without assuming any other property on \(\varphi\) except for it to be in the space \(\Phi_w(\Omega)\), the spaces \(BV(\Omega)\) and \(BV^\varphi(\Omega)\) are not related to each other. Considering \(\Omega \subset \mathbb{R}^n\) bounded, one gets the continuous inclusion \(BV^\varphi(\Omega) \hookrightarrow BV(\Omega)\) assuming (A0) on the function \(\varphi\), see Lemma 2. The fact that \(\varphi\) satisfies (aInc)\(_1\) is already included in the definition of weak \(\Phi\)-function, see Remark 1.

Remark 3. In [27] and [28], in the simple anisotropic case, (i.e. \(\varphi(x,t)\) with linear behavior in \(t\)) the dual function \(\varphi^0\) is used instead of the conjugate function \(\varphi^*\) in order to define a weighted total variation (they call it generalized total variation). But under the assumptions considered in [27] and [28], the two definitions coincide. Indeed, let us assume that \(\varphi: \Omega \times [0, +\infty) \to [0,+\infty)\) is a continuous function which is convex in the second variable, satisfies (A0) and such that \[\varphi(x,t) = t \, \varphi(x,1) \quad \textrm{ for any } x \in \Omega \textrm{ and } t \in [0, +\infty).\] Assume that there exist \(0 < \lambda \le \Lambda < +\infty\) such that \[\lambda t \le \varphi(x,t) \le \Lambda t \qquad \textrm{ for any } x \in \Omega \textrm{ and for any } t \in [0, +\infty).\] Moreover, let us consider \(\psi : \Omega \times \mathbb{R}^n \to [0, +\infty)\) defined as \(\psi(x,\xi) := \varphi(x, |\xi|)\) for any \(x \in \Omega\) and \(\xi \in \mathbb{R}^n\). Note that \(\psi\) is continuous, convex in the second variable, and such that \[\psi(x, t \xi) = \varphi(x, |t \xi|) = |t|\varphi(x, |\xi|) = |t| \psi(x,\xi)\] and \[\lambda |\xi| \le \psi(x, \xi) \le \Lambda |\xi|,\] for every \(t \in \mathbb{R}\), \(\xi \in \mathbb{R}^n\) and every \(x \in \Omega\).
In [27], considered \(u \in BV(\Omega)\) and under these assumptions on \(\psi\), the anisotropic total variation (there called
generalized total variation of \(u\) with respect to \(\psi\) in \(\Omega\)) has been characterized, see in particular [27], as \[\label{Vpsieq} \overline{V_\psi}(u) := \sup \left\{\int_\Omega u \, \mathrm{div}w \, dx \ : \ w \in C^1_c(\Omega;\mathbb{R}^n), \psi^0(x,w(x)) \le 1 \textrm{ for all } x \in \Omega \right\},\qquad{(1)}\] where \(\psi^0\) is the dual function of \(\psi\), namely \[\psi^0(x, \xi^*) := \sup \{(\xi^*, \xi) \ : \ \xi \in \mathbb{R}^n, \psi(x, \xi) \le 1\}.\] Note that in [27] they use the notation \(C^1_0(\Omega;\mathbb{R}^n)\) to indicate the space that we denote with \(C^1_c(\Omega;\mathbb{R}^n)\). Then, given a function \(u \in BV(\Omega)\), the anisotropic total variation of \(u\) with respect to \(\varphi\) in \(\Omega\) introduced in [15] as 2 , coincides with the notion given by [27], namely ?? .
Indeed, by [27] the right-hand side of ?? coincides with \[\label{Vpsieq2} \sup\left\{\int_\Omega u \, \mathrm{div}\, w \, dx \ : \ w \in C^1_c(\Omega;\mathbb{R}^n), \int_\Omega \psi^*(x, w(x)) \, dx < +\infty \right\}.\qquad{(2)}\] Considering the even extension of \(\varphi(x,\cdot)\), for any \(x \in \Omega\), and since \(\varphi(x,\cdot)\) is convex, in view of our definition of \(\psi\) and Lemma 11, which guarantees that \[\label{psi42varphi42mod} \psi^*(x,\cdot)= \varphi^*(x,|\cdot|) \qquad \textrm{ for a.e. } x \in \Omega,\qquad{(3)}\] we have that ?? is equal to \[\label{Vpsieq3} \sup \left\{\int_\Omega u \, \mathrm{div}\, w \, dx \ : \ w \in C^1_c(\Omega;\mathbb{R}^n), \int_\Omega \varphi^*(x, |w(x)|) \, dx < +\infty \right\}.\qquad{(4)}\] Taking into account [27], namely \[\psi^*(x, \xi^*) = \begin{dcases} 0 & \textrm{ if } \psi^0(x, \xi^*) \le 1, \\ +\infty & \textrm{ if } \psi^0(x, \xi^*) > 1, \end{dcases}\] for every \(x \in \Omega\) and every \(\xi^* \in \mathbb{R}^n\), we have that ?? can be equivalently written as \[\label{Vpsieq4} \sup \left\{\int_\Omega u \, \mathrm{div}\, w \, dx \ : \ w \in C^1_c(\Omega;\mathbb{R}^n), \int_\Omega \varphi^*(x, |w(x)|) \, dx \le 1 \right\},\qquad{(5)}\] where ?? has been exploited again. Now, using the definition \[\|w\|_{\varphi^*} := \inf \left\{\lambda > 0 \ : \ \int_\Omega \varphi^* \left(x, \left|\frac{w(x)}{\lambda} \right| \right) \, dx < 1 \right\},\] we get that ?? is equal to 2 .

In order to treat the minimization problem using \(TGV^{\varphi,2}_\alpha\) as a regularizer, see 1 , we need to introduce a variant of the generalized variation associated to \(\varphi\) defined in 2 . For every \(v \in \mathcal{M}(\Omega;\mathbb{R}^m)\) with \(m \in \mathbb{N}^+\), we define \[\label{V95tilde} \widetilde{V}^m_\varphi(v) := \sup \left\{\int_\Omega \psi \cdot dv(x) \ : \ \psi \in C^1_c(\Omega;\mathbb{R}^m), \|\psi\|_{\varphi^*} \le 1 \right\}.\tag{7}\] If we are working on a set \(A \subsetneq \Omega\), then we will denote it by \(\widetilde{V}^m_\varphi(v;A)\).

Remark 4. Note that for any \(u \in BV(\Omega)\), in order to know that at least that \(Du \in \mathcal{M}(\Omega;\mathbb{R}^n)\), an integration by parts yields \[\widetilde{V}^n_\varphi(Du) = V_\varphi(u),\] where the second term is defined as in 2 , namely the anisotropic total variation introduced in [15].

In the following, we will either use the above notation for \(m = n\) or \(m = n \times n\). Note that in the latter case the scalar product coincides with the one underlying the Frobenius norm in \(\mathbb{R}^{n \times n}\), namely reading for any \(v \in \mathcal{M}(\Omega;\mathbb{R}^{n\times n}_{\rm sym})\) as \[\widetilde{V}^{n\times n}_\varphi(v) = \sup \left\{\int_\Omega \psi : dv(x) \ : \ \psi \in C^1_c(\Omega;\mathbb{R}^{n \times n}_{\rm sym}), \|\psi\|_{\varphi^*} \le 1 \right\}.\] Considering test functions in \(C^1_c(\Omega;\mathbb{R}^{n \times n}_{\rm sym})\) or \(C^1_c(\Omega;\mathbb{R}^{n \times n})\) is equivalent when \(v \in \mathcal{M}(\Omega;\mathbb{R}^{n \times n}_{\rm sym})\), since the scalar product only takes into account the symmetric part of the test function. In general the notation \(\widetilde{V}^{n\times n}_\varphi\) could be used also for measures in \(\mathcal{M}(\Omega;\mathbb{R}^{n \times n})\) without any symmetry property, but that is not the case in this work. We avoid specifying the symmetry in order to lighten the notation.

Remark 5. Considering an open set \(\Omega \subset \mathbb{R}^n\), if \(v\) has the special structure \(v = Du - w\), with \(u,w \in \mathcal{M}(\Omega;\mathbb{R}^n)\), an integration by parts entails \[\widetilde{V}^n_\varphi(Du - w) = \sup \left\{\int_\Omega u \, \mathrm{div}\, \psi \, dx + \int_\Omega \psi \cdot dw(x) \ : \ \psi \in C^1_c(\Omega;\mathbb{R}^n), \|\psi\|_{\varphi^*} \le 1 \right\},\] see also [29]. Note that it is possible to rewrite it like this by [30], indeed having \(u \in \mathcal{M}(\Omega;\mathbb{R}^n)\) and \(Du = v + w \in \mathcal{M}(\Omega;\mathbb{R}^n)\) implies that \(u \in BV(\Omega)\).
In the same way, if \(\Omega\) is also bounded and \(C^1\), by [31], having \(Eu \in \mathcal{M}(\Omega;\mathbb{R}^{n \times n}_{\rm sym})\) implies \(u \in BD(\Omega)\). Namely, in the special case \(v = Eu \in \mathcal{M}(\Omega;\mathbb{R}^{n\times n}_{\rm sym})\), 7 rewrites as \[\label{var95grad95sym} \widetilde{V}^{n \times n}_\varphi(E u) = \sup \left\{\int_\Omega E^* \psi \cdot u \, dx \ : \ \psi \in C^1_c(\Omega;\mathbb{R}^{n \times n}_{\rm sym}), \|\psi\|_{\varphi^*} \le 1 \right\},\qquad{(6)}\] where \(E^* \psi = -\mathrm{div}\, \psi\), and \(\mathrm{div}\) denotes the matrix divergence of order \(1\), computed row-wise (i.e. with respect to columns), cf. 6 .

Remark 6. We point out that, here and in the following, the set \(\Omega\) is only required to have a \(C^1\) boundary in order to apply [31]. The extension of the above result to the case of Lipschitz boundary will be the subject of the forthcoming paper [32].

Thanks to the definition of generalized variation associated to \(\varphi\) in 7 , we can now define the space of bounded deformation fields with generalized Orlicz growth, namely \[\label{BD95phi} BD^\varphi(\Omega) := \left\{ w \in L^\varphi(\Omega;\mathbb{R}^n) \ : \ \|w\|_{BD^\varphi(\Omega)} < +\infty \right\},\tag{8}\] where \[\|w\|_{BD^\varphi(\Omega)} := \|w\|_\varphi+ \widetilde{V}^{n \times n}_\varphi(Ew).\]

Lemma 3. Let \(\Omega\) be a bounded and open set and let \(\varphi\in \Phi_w(\Omega)\). Then \(\widetilde{V}^{n\times n}_\varphi\) is a seminorm.

Proof. We need to check that \(\widetilde{V}^{n \times n}_\varphi\) satisfies conditions (N2) and (N3) in Definition 1. The homogeneity property (N2) is obvious, in fact \(\widetilde{V}^{n \times n}_\varphi(av) = |a|\widetilde{V}^{n \times n}_\varphi(v)\) for every \(v \in \mathcal{M}(\Omega;\mathbb{R}^{n \times n}_{\rm sym})\) and any \(a \in \mathbb{R}\). In order to check (N3), namely the triangle inequality, let us consider \(v,w \in \mathcal{M}(\Omega;\mathbb{R}^{n \times n}_{\rm sym})\). We have \[\int_\Omega \psi : d(v+w) = \int_\Omega \psi : dv + \int_\Omega \psi : dw \le \widetilde{V}^{n \times n}_\varphi(v) + \widetilde{V}^{n \times n}_\varphi(w)\] for any \(\psi \in C^1_c(\Omega;\mathbb{R}^{n \times n}_{\rm sym})\) such that \(\|\psi\|_{\varphi^*} \le 1\). By taking the supremum over such \(\psi\) we have the thesis. ◻

This result entails that if \(\varphi\in \Phi_w(\Omega)\), then \(\|\cdot\|_{BD^\varphi(\Omega)}\) is a quasinorm. In addition, if \(\varphi\) verifies (A0), or if \(\varphi\in \Phi_c(\Omega)\), then \(\|\cdot\|_{BD^\varphi(\Omega)}\) is a norm, see [11].

Remark 7. As for \(BV(\Omega)\) and \(BV^\varphi(\Omega)\), see Remark 2, without assuming any other property on \(\varphi\) except for it to be in the space \(\Phi_w(\Omega)\), the spaces \(BD(\Omega)\) and \(BD^\varphi(\Omega)\) are not related to each other. Considering a set \(\Omega\) bounded, one gets the continuous immersion \(BD^\varphi(\Omega) \hookrightarrow BD(\Omega)\) assuming (A0) on the function \(\varphi\), see Lemma 2.

Remark 8. Since a priori we define the space starting from a function \(\varphi\in \Phi_w(\Omega)\), one may ask whether the space defined starting from \(\varphi^{**}\), always a convex \(\Phi\)-function, see Lemma 1(i), namely \[BD^{\varphi^{**}}(\Omega) := \left\{ w \in L^{\varphi^{**}}(\Omega;\mathbb{R}^n) \ : \ \|w\|_{BD^{\varphi^{**}}(\Omega)} < +\infty \right\},\] with \[\|w\|_{BD^{\varphi^{**}}(\Omega)} := \|w\|_{\varphi^{**}} + \widetilde{V}^{n \times n}_{\varphi^{**}}(E w),\] differs from the space \(BD^\varphi(\Omega)\). Thanks to [11] it is possible to prove that \(BD^\varphi(\Omega) = BD^{\varphi^{**}}(\Omega)\), indeed the norms \(\|\cdot\|_\varphi\) and \(\|\cdot\|_{\varphi^{**}}\) are equivalent and the variations coincide, implying that the norms \(\|\cdot\|_{BD^\varphi(\Omega)}\) and \(\|\cdot\|_{BD^{\varphi^{**}}(\Omega)}\) are equivalent as well. Note that some results needed for this consideration, e.g. [11] or [11], are only stated without the \(x\)-dependence, but they hold true even in the generalized Orlicz setting, while others, e.g. [11] or [11], are already stated with the \(x\)-dependence.

With the next example we want to show that, at least considering some \(\Phi\)-functions for which a Korn-type inequality does not hold, the space \(BV^\varphi(\Omega)\) is a proper subset of \(BD^\varphi(\Omega)\). The example follows the lines of [33], exploiting a construction developed in [34], see also [31].

Example 1. We first remark that 5 , \(BV^\varphi(\Omega;\mathbb{R}^n) \subset BD^\varphi(\Omega)\). We show below that, at least in some situations, the inclusion is strict.
Indeed, without loss of generality we can consider an open set \(\Omega \supset \supset B\), with \(B\) a small ball centered at \(0\), and a function such that \(\varphi(x,t) = |t|\) for a.e. \(x \in B\). Then, the same field \(u\) provided by [33], which has compact support in \(B\) and is in \(BD(\Omega) \setminus BV(\Omega; \mathbb{R}^n)\) provides the desired example of a field in \(BD^\varphi(\Omega) \setminus BV^\varphi(\Omega; \mathbb{R}^n)\).
Note that many \(\Phi\)-functions, e.g. many out of the ones listed in [11], fulfill the assumption of coinciding with \(|t|\) on a ball \(B\), meaning that this example holds true in many contexts.

The validity of Korn-type inequalities in the context of Orlicz spaces has been studied in different contexts, for example [35] where it was shown that the \(\nabla_2\)-condition and the \(\Delta_2\)-condition are sharp conditions for Korn-type inequalities in Orlicz spaces (for the definition of these two conditions see [14] or [36]). In [37] the validity of Korn-type inequalities for Young functions was shown, while in [36] the equivalence of these Korn-type inequalities to two balance conditions, similar to the \(\nabla_2\)- and \(\Delta_2\)-conditions, has been proven. Note that, under the assumptions considered in [37] and [36], their definition of Young functions and our definition of convex \(\Phi\)-functions, when independent of \(x\), are equivalent. In the context of generalized Orlicz spaces (or Musielak-Orlicz spaces) instead, there are results in the case of variable exponents Lebesgue spaces when the exponent is bounded away from \(1\), see [38] and [10], but not much else is known.

3.1 Weighted modular and associated norms↩︎

Analogously to [15], we define the anisotropic modular, or dual modular, associated to the anisotropic total variation \(\widetilde{V}^{n \times n}_\varphi\), defined as in 7 . In particular, we define for any \(v \in \mathcal{M}(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\) \[\label{varrhosymvarphi} \varrho_{\widetilde{V}^{n \times n}_\varphi}(v) := \sup \left\{ \int_\Omega \psi : dv(x) - \int_\Omega\varphi^*(x,|\psi|) \, dx \ : \ \psi \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym}) \right\}.\tag{9}\]

Lemma 4. Let \(\Omega\) be a bounded and open set and let \(\varphi\in \Phi_w(\Omega)\). Then \(\varrho_{\widetilde{V}^{n \times n}_\varphi}\) is a semimodular.

Proof. We need to check that \(\varrho_{\widetilde{V}^{n \times n}_\varphi}\) satisfies the conditions of Definition 2. Conditions (i) and (iii) are obvious, respectively thanks to the properties of \(\varphi^*\) and that if \(\psi\) is an admissible test function, then also \(-\psi\) is. Condition (iv) is also easy to prove using the linearity of the first integral in the definition of \(\varrho_{\widetilde{V}^{n \times n}_\varphi}\) and the fact that the supremum of a sum is smaller than the sum of the suprema.
We are only left to prove condition (ii), namely that the function \(\lambda \mapsto \varrho_{\widetilde{V}^{n \times n}_\varphi}(\lambda v)\) on \([0, +\infty)\) is non-decreasing for any \(v \in \mathcal{M}(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\). Indeed, if \(\lambda = 0\), then \(\varrho_{\widetilde{V}^{n \times n}_\varphi}(\lambda v) = 0\), but we know already that for any other \(\mu > \lambda = 0\) then \(\varrho_{\widetilde{V}^{n \times n}_\varphi}(\mu v) \ge 0\), since we can always choose \(\psi = 0\) as a test function. Let us instead fix \(0 < \lambda < \mu < +\infty\). We point out that since both \(\psi\) and \(-\psi\) are possible test functions, we can also only consider without loss of generality functions in the space \[\overline{T}^+ := \left\{ \psi \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym}) \ : \ \int_\Omega \psi : dv(x) \ge 0 \right\},\] since the second integral in the definition of the modular is always zero or negative. Then we have \[\begin{align} \varrho_{\widetilde{V}^{n \times n}_\varphi}(\lambda v) &= \sup \left\{ \int_\Omega \lambda\psi : dv(x) - \int_\Omega\varphi^*(x,|\psi|) \, dx \ : \ \psi \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym}) \right\} \\ &= \sup \left\{ \lambda \int_\Omega \psi : dv(x) - \int_\Omega\varphi^*(x,|\psi|) \, dx \ : \ \psi \in \overline{T}^+ \right\} \\ &\le \sup \left\{ \mu \int_\Omega \psi : dv(x) - \int_\Omega\varphi^*(x,|\psi|) \, dx \ : \ \psi \in \overline{T}^+ \right\} \\ &= \sup \left\{\int_\Omega \mu\psi : dv(x) - \int_\Omega\varphi^*(x,|\psi|) \, dx \ : \ \psi \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym}) \right\} = \varrho_{\widetilde{V}^{n \times n}_\varphi}(\mu v), \end{align}\] for any \(v \in \mathcal{M}(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\), proving condition (ii).
Since condition (iv) holds with the constant \(\beta = 1\), then \(\varrho_{\widetilde{V}^{n \times n}_\varphi}\) is a semimodular. ◻

Note that, analogously to the anisotropic total variation, see Remark 5, considering a bounded, \(C^1\) and open set \(\Omega \subset \mathbb{R}^n\), in the special case \(v = Eu \in \mathcal{M}(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\) this rewrites as \[\label{varrhosymvarphi95Eu} \varrho_{\widetilde{V}^{n \times n}_\varphi}(Eu) = \sup \left\{ \int_\Omega \left[u \cdot E^*\psi - \varphi^*(x,|\psi|) \right] dx \ : \ \psi \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym}) \right\},\tag{10}\] where \(E^*\psi = - \mathrm{div}\, \psi\).

Remark 9. On the one hand, note that testing 10 with \(\psi = 0\), we get that the supremum in \(\varrho_{\widetilde{V}^{n \times n}_\varphi}\) is always non-negative. On the other hand, if \(\varrho_{\varphi^*}(|\psi|) = +\infty\), then the integral in 10 is equal to \(-\infty\), since the first term is finite. Here we used again [31] yielding that \(Eu \in \mathcal{M}(\Omega;\mathbb{R}^{n\times n}_{\rm sym})\) implies \(u \in BD(\Omega)\).

This means that test functions \(\psi\) with \(\varrho_{\varphi^*}(|\psi|) = +\infty\) can be omitted and \(\varrho_{\widetilde{V}^{n \times n}_\varphi}\) can be equivalently rewritten as \[\varrho_{\widetilde{V}^{n \times n}_\varphi}(Eu) := \sup \left\{ \int_\Omega \left[u \cdot E^*\psi - \varphi^*(x,|\psi|) \right] dx \ : \ \psi \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym}), \varrho_{\varphi^*}(|\psi|) < +\infty \right\}.\] The analogous reasoning can be repeated for the general expression of \(\varrho_{\widetilde{V}^{n \times n}_\varphi}\) in 9 .

In Lemma 4 we showed that \(\varrho_{\widetilde{V}^{n \times n}_\varphi}\) is a semimodular, implying that, thanks to the Luxemburg method, we obtain a seminorm \(\|\cdot\|_{\varrho_{\widetilde{V}^{n \times n}_\varphi}}\). The following result relates the anisotropic variation we introduced in 7 with the seminorm associated to the modular defined in 9 , similarly to [15].

Lemma 5. Let \(\Omega \subset \mathbb{R}^n\) be a bounded, connected, \(C^1\) and open set. Let \(\varphi\in \Phi_w(\Omega)\) and let \(u \in BD^\varphi(\Omega)\). Then \[\|Eu\|_{\varrho_{\widetilde{V}^{n \times n}_\varphi}} \le \widetilde{V}^{n\times n}_\varphi(Eu) \le 2 \, \|Eu\|_{\varrho_{\widetilde{V}^{n \times n}_\varphi}}.\]

Proof. Since \(u \in BD^\varphi(\Omega)\), the case \(\widetilde{V}^{n \times n}_\varphi(E u) = +\infty\) is excluded a priori. Moreover, if \(\widetilde{V}^{n \times n}_\varphi(Eu) = 0\), then \(\varrho_{\widetilde{V}^{n \times n}_\varphi}(\frac{Eu}{\lambda}) = 0\) for every \(\lambda > 0\), yielding \(\|Eu\|_{\varrho_{\widetilde{V}^{n \times n}_\varphi}} = 0\). Lastly, since the claim is homogeneous, the case \(\widetilde{V}^{n \times n}_\varphi(Eu) \in (0,+\infty)\) can always be reduced to \(\widetilde{V}^{n \times n}_\varphi(Eu) = 1\).
By the definition of \(\widetilde{V}^{n \times n}_\varphi\), see ?? , for any \(w \in C^1_c(\Omega;\mathbb{R}^{n \times n}_{\rm sym})\) such that \(\|w\|_{\varphi^*} \le 1\) it follows that \[\int_\Omega u \cdot E^*w \, dx \le \widetilde{V}^{n \times n}_\varphi(Eu) \, \|w\|_{\varphi^*} = \|w\|_{\varphi^*} \le 1 + \varrho_{\varphi^*}(|w|),\] where the last inequality comes from a general property of the Luxemburg norms, see [10]. Thus, by the definition of \(\varrho_{\widetilde{V}^{n \times n}_\varphi}\) given in 9 , \[\varrho_{\widetilde{V}^{n \times n}_\varphi}(Eu) \le \sup_{w \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym})} [1 + \varrho_{\varphi^*}(|w|) - \varrho_{\varphi^*}(|w|)] = 1,\] implying \(\|E u\|_{\varrho_{\widetilde{V}^{n \times n}_\varphi}} \le 1\). This concludes the proof of the first inequality.
We next establish the opposite inequality \(2 \, \|Eu\|_{\varrho_{\widetilde{V}^{n \times n}_\varphi}} \ge 1\), which is equivalent to \[\label{eq:mod95magg} \varrho_{\widetilde{V}^{n \times n}_\varphi} \left(\frac{2Eu}{\lambda} \right) > 1 \qquad \textrm{ for every } \lambda < 1.\tag{11}\] Indeed, on the one hand if \(\|v\|_{\varrho_{\widetilde{V}^{n \times n}_\varphi}} \ge 1\), by Definition 3 we have that \[\inf \left\{ k > 0 : \varrho_{\widetilde{V}^{n \times n}_\varphi}\left( \frac{v}{k} \right) \le 1 \right\} \ge 1,\] entailing that for any \(\lambda < 1\) it holds \(\|v\|_{\varrho_{\widetilde{V}^{n \times n}_\varphi}} > \lambda\). By [11] then this implies that \(\varrho_{\widetilde{V}^{n \times n}_\varphi}(\tfrac{v}{\lambda}) > 1\) for any \(\lambda < 1\). On the other hand, having \(\varrho_{\widetilde{V}^{n \times n}_\varphi}(\tfrac{v}{\lambda}) > 1\) for any \(\lambda < 1\) entails that, by Definition 3, \[\|v\|_{\varrho_{\widetilde{V}^{n \times n}_\varphi}} = \inf \left\{ \lambda > 0 \ : \ \varrho_{\widetilde{V}^{n \times n}_\varphi}\left( \frac{v}{\lambda} \right) \le 1 \right\} \ge 1\] directly from the definition of \(\inf\). Now, by [11] and thanks to the fact that \(\varphi^* \in \Phi_c(\Omega)\), see Lemma 1(iii), \(\varrho_{\varphi^*}(|w|) \le 1\) if and only if \(\|w\|_{\varphi^*} \le 1\), and we can immediately conclude that \[\begin{align} \varrho_{\widetilde{V}^{n \times n}_\varphi} \left(\frac{2Eu}{\lambda} \right) & \ge \sup \left\{ \int_\Omega \left[ \frac{2 u}{\lambda} \cdot E^* w - \varphi^*(x, |w|) \right] dx \ : \ w \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym}), \|w\|_{\varphi^*}\le 1 \right\} \\ & \ge \frac{2}{\lambda} \, \widetilde{V}^{n \times n}_\varphi(Eu) - 1 > 1. \end{align}\] This proves 11 and yields the thesis. ◻

The following result, instead, relates the seminorm \(\widetilde{V}^{n\times n}_\varphi\), see Lemma 3, with the norm of the symmetrized gradient, in particular considering the norm associated to the associate space \((L^{\varphi^*}(\Omega; \mathbb{R}^{n \times n}))'\). Obviously, in the associate space norm we test with functions in \(L^{\varphi^*}(\Omega; \mathbb{R}^{n \times n})\), while in \(\widetilde{V}^{n\times n}_\varphi\) the test functions are smooth. This means that some approximation is needed, but a priori we are not able to use density in \(L^{\varphi^*}(\Omega; \mathbb{R}^{n \times n})\) since \(\varphi^*\) is in general not doubling.
In order to state the result, we need the classical decomposition of the symmetrized gradient, see for instance [39], \[Eu = \mathcal{E}u \, \mathcal{L}^n + E^s u.\] Moreover, let us remark the definition of the space \(LD_{\rm loc}(\Omega)\), namely \[LD_{\rm loc}(\Omega) := \{u \in L^1_{\rm loc}(\Omega; \mathbb{R}^n) \ : \ Eu \in L^1_{\rm loc}(\Omega;\mathbb{R}^{n \times n}_{\rm sym})\}.\] The result, analogous to [15], reads as follows.

Lemma 6. Let \(\Omega \subset \mathbb{R}^n\) be a bounded, connected and open set. Let \(\varphi\in \Phi_w(\Omega)\) and let \(u \in LD_{\rm loc}(\Omega)\). Then \[\widetilde{V}^{n \times n}_\varphi(Eu) \le \| \mathcal{E}u \|_{(L^{\varphi^*}(\Omega; \mathbb{R}^{n \times n}))'}.\] Moreover, if additionally \(C^1_c(\Omega; \mathbb{R}^{n \times n})\) is dense in \(L^{\varphi^*}(\Omega; \mathbb{R}^{n \times n})\), then \[\widetilde{V}^{n \times n}_\varphi(Eu) = \| \mathcal{E}u \|_{(L^{\varphi^*}(\Omega;\mathbb{R}^{n \times n}))'}.\]

Proof. Since \(u \in LD_{\rm loc}(\Omega)\), it follows from the definition of \(\widetilde{V}^{n \times n}_\varphi\) that \[\label{eq:integrationByParts} \widetilde{V}^{n \times n}_\varphi(Eu) = \sup \left\{ \int_\Omega \mathcal{E}u : w \, dx \ : \ w \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym}), \|w\|_{\varphi^*} \le 1 \right\}.\tag{12}\] The definition of the associate space norm (extended in a standard way to the vector valued setting), see Definition 10, implies that, chosen \(w \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\) such that \(\|w\|_{\varphi^*} \le 1\), \[\int_\Omega \mathcal{E}u : w \, dx \le \int_\Omega |\mathcal{E}u| |w| \, dx \le \| \mathcal{E}u\|_{(L^{\varphi^*}(\Omega;\mathbb{R}^{n \times n}))'} \|w\|_{\varphi^*}.\] Taking the supremum over those \(w\), we conclude that \(\widetilde{V}^{n \times n}_\varphi(\mathcal{E}u) \le \|\mathcal{E}u\|_{(L^{\varphi^*}(\Omega;\mathbb{R}^{n \times n}))'}\).
Under the density assumption, namely \(C^1_c(\Omega;\mathbb{R}^{n \times n})\) being dense in \(L^{\varphi^*}(\Omega;\mathbb{R}^{n \times n})\), we next show the opposite inequality, namely \(\|\mathcal{E}u\|_{(L^{\varphi^*}(\Omega;\mathbb{R}^{n \times n}))'} \le \widetilde{V}^{n \times n}_\varphi(\mathcal{E}u)\). Indeed, let \(w \in L^{\varphi^*}(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\) with \(\|w\|_{\varphi^*} = 1\) and let \((w_j)_j \subset C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\) be a sequence such that \(w_j \to w\) in \(L^{\varphi^*}(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\) and pointwise almost everywhere. Since we know that also \(w_j/\|w_j\|_{\varphi^*} \to w\) in \(L^{\varphi^*}(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\), we may assume that \(\|w_j\|_{\varphi^*} = 1\) for every \(j \in \mathbb{N}\). By Fatou’s Lemma, we have \[\liminf_{j \to \infty} \int_\Omega \mathcal{E}u : w_j \, dx \ge \int_\Omega \mathcal{E}u : w \, dx.\] So, from 12 follows that \[\widetilde{V}^{n \times n}_\varphi(\mathcal{E}u) \ge \sup \left\{ \int_\Omega \mathcal{E}u : w \, dx \ : \ w \in L^{\varphi^*}(\Omega; \mathbb{R}^{n \times n}_{\rm sym}), \| w\|_{\varphi^*} \le 1 \right\}.\] Let \(h\in L^{\varphi^*}(\Omega)\). We can set \[w := \begin{dcases} \frac{\mathcal{E} u}{|\mathcal{E} u|} \, h & \textrm{ if } |\mathcal{E} u|\neq 0, \\ 0 & \textrm{ otherwise.} \end{dcases}\] This gives \[\widetilde{V}^{n \times n}_\varphi( \mathcal{E}u) \ge \sup \left\{ \int_\Omega |\mathcal{E}u| \, h \, dx \ : \ h \in L^{\varphi^*}(\Omega), \| h\|_{\varphi^*}\le 1 \right\} = \| \mathcal{E}u \|_{(L^{\varphi^*}(\Omega;\mathbb{R}^{n \times n}))'}.\] Hence \(\widetilde{V}^{n \times n}_\varphi(\mathcal{E}u) = \| \mathcal{E}u \|_{(L^{\varphi^*}(\Omega;\mathbb{R}^{n \times n}))'}\). This concludes the proof. ◻

In the second part of Lemma 6 we assume that \(C^1_c(\Omega; \mathbb{R}^{n \times n})\) is dense in \(L^{\varphi^*}(\Omega; \mathbb{R}^{n \times n})\). If \(\varphi^*\) satisfies (A0) and (aDec), then this holds automatically by [11] (note that this theorem is stated for scalar functions, but the vectorial case can be deduced working components-wise thanks to the equivalence of euclidean norms in spaces of finite dimension). We know already that, by Lemma 1(ii) & (iv), \(\varphi^*\) satisfies these conditions if and only if \(\varphi\) satisfies (A0) and (aInc), meaning that it is not clear whether this density is there or not in the subsets where the function \(\varphi\) is linear.

3.2 Explicit expression of the anisotropic modular↩︎

In this subsection we derive an explicit formula for the dual modular \(\varrho_{\widetilde{V}^{n \times n}_\varphi}\) in terms of the modular \(\varrho_\varphi\) of the absolutely continuous part of the symmetric gradient and its singular part with weight given by the recession function, see Definition 8. In particular, analogously to [15], the main result of the section reads as follows.

Theorem 10. Let \(\Omega \subset \mathbb{R}^n\) be a bounded, connected and open set. Let \(\varphi\in \Phi_c(\Omega) \cap C(\Omega \times [0,+\infty))\) satisfy (A0), (aDec) and restricted (VA1). If \(u \in BD(\Omega)\), then \[\varrho_{\widetilde{V}^{n \times n}_\varphi}(Eu) = \varrho_\varphi(|\mathcal{E}u|) + \int_\Omega \varphi^\infty \, d|E^su|.\] 0◻

The following result is an immediate consequence of Theorem 10, when the \(\Phi\)-function \(\varphi\) does not depend on \(x\), i.e. when it is a classical Orlicz \(\Phi\)-function.

Corollary 1. Let \(\varphi\in \Phi_c\) be independent of \(x\) and satisfy \((aDec)\). If \(u \in BD(\Omega)\), then \[\varrho_{\widetilde{V}^{n \times n}_\varphi}(Eu) = \varrho_\varphi(|\mathcal{E}u|) + \varphi^\infty |E^su|\] in any Borel subset of \(\Omega\) and so

  • \(BD^\varphi(\Omega) = BD(\Omega)\) if \(\varphi^\infty < \infty\),

  • \(BD^\varphi(\Omega) = LD^\varphi(\Omega)\) if \(\varphi^\infty = \infty\).

0◻

We observe that in (2) we cannot say that \(BD^\varphi(\Omega) = W^{1,\varphi}(\Omega)\), see e.g. [40], since the latter equality would follow from a Korn’s inequality. But, as proven in [35], the latter, for classical convex \(\Phi\)-functions holds if and only if \(\varphi\) satisfies both \(\nabla_2\) and \(\Delta_2\) conditions. Now, recall first that \(\varphi\in \Phi_c\) implies that \(\varphi= \varphi^{**}\) and that the \(\nabla_2\) condition for \(\varphi\) is equivalent, by [10], to \(\varphi^*\) satisfying \(\Delta_2\). Furthermore, [10] states that \(\varphi\in \Phi_w\) satisfies \(\nabla_2\) if and only if it satisfies \((aInc)\). By assumption \(\varphi\) satisfies \((aDec)\) but not \((aInc)\), thus \(\varphi\) satisfies \(\Delta_2\) but not \(\nabla_2\). As a particular case, one can consider \(\varphi(t) = t \log(1+t)\), see [10].

In order to prove Theorem 10 we need the following lemma, proven in [15], which shows the importance of the recession function \(\varphi^\infty\), see Definition 8.

Lemma 7. Let \(\Omega \subset \mathbb{R}^n\) be a bounded, connected and open set. Let \(\varphi\in \Phi_c(\Omega)\) satisfy restricted (VA1). If \(\psi \in C(\Omega)\) with \(\varrho_{\varphi^*}(\psi) < +\infty\), then \(|\psi| \le \varphi^\infty\). 0◻

In [12] it was shown that log-Hölder continuity is not sufficient when working in \(BV^{p(\cdot)}\) and here, similarly, the (A1) condition (corresponding to log-Hölder continuity in the general case of convex \(\Phi\)-functions) is not sufficient in the previous result, as in the following ones, in view of [15]. That is why we make use of the restricted (VA1) condition, which corresponds to strong log-Hölder continuity in the case of the variable exponent function, see [15].
In the next two results we consider the singular and absolutely continuous parts of the symmetric gradient separately. Then we combine them to handle the whole function in Theorem 10. Since the proofs are similar to the ones of [15], they can be found in our setting in Appendix 5.

Proposition 11. Let \(\Omega \subset \mathbb{R}^n\) be a bounded, connected and open set. Let \(\varphi\in \Phi_c(\Omega)\) satisfy (A0), (aDec) and restricted (VA1). If \(u \in BD(\Omega)\), then \[\sup_{\psi \in T^\varphi} \int_\Omega \psi : dE^su = \int_\Omega \varphi^\infty \, d|E^su|,\] where \[\label{def95test} T^\varphi:= \left\{ \psi \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym}) \ : \ \varrho_{\varphi^*}(|\psi|) < +\infty \right\}.\qquad{(7)}\] 0◻

Proposition 12. Let \(\Omega \subset \mathbb{R}^n\) be a bounded, connected and open set. Let \(\varphi\in \Phi_c(\Omega) \cap C(\Omega \times [0,+\infty))\) satisfy (A0) and (aDec). If \(u \in BD(\Omega)\), then \[\sup_{\psi \in T^\varphi} \int_\Omega \left[\mathcal{E}u : \psi - \varphi^*(x,|\psi|) \right] dx = \varrho_\varphi(|\mathcal{E}u|).\] 0◻

Note that removing the assumption that \(\varphi\) is continuous in both variables is not trivial.
Having the two previous propositions at hands, we can now prove Theorem 10, exploiting the ideas used in [15].

Proof of Theorem 10. Since \(Eu = \mathcal{E}u + E^su\) and since \(u \in BD(\Omega)\), an integration by parts implies that \[\int_\Omega u \cdot E^*\psi \, dx = \int_\Omega \psi : dEu = \int_\Omega \mathcal{E}u : \psi \, dx + \int_\Omega \psi : dE^su\] for any \(\psi \in T^\varphi\). Hence, the claim follows from Propositions 11 and 12 once we prove that \[\begin{align} &\sup_{\psi \in T^\varphi} \left\{\int_\Omega [\mathcal{E}u : \psi - \varphi^*(x, |\psi|)] \, dx + \int_\Omega \psi : dE^su \right\} \\ = &\sup_{\psi \in T^\varphi} \int_\Omega [\mathcal{E}u : \psi - \varphi^*(x, |\psi|)] \, dx + \sup_{\psi \in T^\varphi} \int_\Omega \psi : dE^su. \end{align}\] The first inequality, namely “\(\le\)”, is obvious, so we only need to prove the opposite one.
Let us first assume that \(\varrho_\varphi(|\mathcal{E}u|) + \int_\Omega \varphi^\infty \, d|E^su| < +\infty\) and let us fix \(\varepsilon> 0\). We can always choose \(\psi_1, \psi_2\in T^\varphi\) such that \[\label{eq95aux:choice951} \sup_{\psi \in T^\varphi} \int_\Omega [\mathcal{E}u : \psi - \varphi^*(x, |\psi|)] \, dx \le \int_\Omega [\mathcal{E}u : \psi_1 - \varphi^*(x, |\psi_1|)] \, dx + \varepsilon< +\infty\tag{13}\] and \[\label{eq95aux:choice952} \sup_{\psi \in T^\varphi} \int_\Omega \psi : dE^su \le \int_\Omega \psi_2 : dE^su + \varepsilon< +\infty.\tag{14}\] Since \(u \in BD(\Omega)\) and \(\psi_1, \psi_2 \in T^\varphi\), we have \(|\mathcal{E}u|\,|\psi_i| \in L^1(\Omega)\) and \(\varrho_{\varphi^*}(|\psi_i|) < +\infty\) for \(i \in \{1,2\}\). Thus, by the absolute continuity of the integral, we can find \(\delta > 0\) such that \[\label{eq95aux:stima951952} \left| \int_{\Omega \setminus \Omega_1} [\mathcal{E}u : \psi_i - \varphi^*(x, |\psi_i|)] \, dx \right| \le \int_{\Omega \setminus \Omega_1} [|\mathcal{E}u| \, |\psi_i| + \varphi^*(x, |\psi_i|)] \, dx \le \varepsilon\tag{15}\] for \(i \in \{1,2\}\) and any \(\Omega_1 \subset \Omega\) with \(|\Omega \setminus \Omega_1|< \delta\) and \[\left|\int_{\Omega \setminus \Omega_2} \psi_2 : dE^su \, \right| \le \int_{\Omega \setminus \Omega_2} \varphi^\infty \, d|E^su| \le \varepsilon\] for any \(\Omega_2 \subset \Omega\) with \(|E^su|(\Omega \setminus \Omega_2) < \delta\).
Since \({\rm supp}\, E^su\) is a set with zero Lebesgue measure, we can find a finite collection of open rectangles \(Q_i \subset \Omega\) such that \(|E^su|(\bigcup_{i \in \mathbb{N}} Q_i) > |E^su|(\Omega) - \delta\) and \(|\bigcup_{i \in \mathbb{N}} 2Q_i| < \delta\). Let us now choose \(\theta \in C^1_c(\Omega)\) such that \(0 \le \theta \le 1\), \(\theta=1\) in \(\Omega_2 := \bigcup_{i \in \mathbb{N}} Q_i\) and \(\theta = 0\) in \(\Omega_1 := \Omega \setminus \bigcup_{i \in \mathbb{N}} 2Q_i\), and let us fix \(\psi_\varepsilon:= \theta \, \psi_2 + (1 - \theta) \psi_1 \in C^1_c(\Omega; \mathbb{R}^{n\times n}_{\rm sym})\). Since \(\psi_\varepsilon\) is a pointwise convex combination, we get \[\varphi^*(\cdot, |\psi_\varepsilon|) \le \varphi^*(\cdot, \max\{|\psi_2|, |\psi_1|\}) \le \varphi^*(\cdot, |\psi_2|) + \varphi^*(\cdot, |\psi_1|).\] This yields that \(\varrho_{\varphi^*}(|\psi_\varepsilon|) \le \varrho_{\varphi^*}(|\psi_2|) + \varrho_{\varphi^*}(|\psi_1|) < +\infty\), yielding \(\psi_\varepsilon\in T^\varphi\). We then get that \[\begin{align} \varrho_{\widetilde{V}^{n \times n}_\varphi}(Eu) &= \sup_{\psi \in T^\varphi} \left\{\int_\Omega [\mathcal{E}u : \psi - \varphi^*(x, |\psi|)] \, dx + \int_\Omega \psi : dE^su \right\} \\ &\ge \int_\Omega [\mathcal{E}u : \psi_\varepsilon- \varphi^*(x, |\psi_\varepsilon|)] \, dx + \int_\Omega \psi_\varepsilon: dE^su \\ &\ge \int_{\Omega_1} [\mathcal{E}u : \psi_1 - \varphi^*(x, |\psi_1|)] \, dx + \int_{\Omega_2} \psi_2 : dE^su - c_\theta \\ &\ge \int_\Omega [\mathcal{E}u : \psi_1 - \varphi^*(x, |\psi_1|)] \, dx + \int_\Omega \psi_2 : dE^su - 5\varepsilon, \end{align}\] where, thanks to 15 , \[c_\theta := \int_{\Omega \setminus \Omega_1} [|\mathcal{E}u| \, |\psi_\varepsilon| + \varphi^*(x, |\psi_\varepsilon|)] \, dx + \int_{\Omega \setminus\Omega_2} |\psi_\varepsilon|\, d|E^su| \le 3\varepsilon\] by the absolute integrability assumptions and by Lemma 7. By the choices of \(\psi_1\) and \(\psi_2\), namely 13 and 14 , we then get \[\varrho_{\widetilde{V}^{n \times n}_\varphi}(Eu) \ge \sup_{\psi \in T^\varphi} \int_\Omega [\mathcal{E}u : \psi - \varphi^*(x, |\psi|)] \, dx + \sup_{\psi \in T^\varphi} \int_\Omega \psi : dE^su - 7 \varepsilon.\] The lower bound follows by letting \(\varepsilon\to 0^+\). This concludes the proof in the case \(\varrho_\varphi(|\mathcal{E}u|) + \int_\Omega \varphi^\infty \, d|E^su| < +\infty\).
If, instead, \(\varrho_\varphi(|\mathcal{E}u|) = +\infty\) and \(\int_\Omega \varphi^\infty \, d|E^su| < +\infty\), we simply estimate \[\begin{align} &\sup_{\psi \in T^\varphi} \left\{ \int_\Omega [\mathcal{E}u : \psi - \varphi^*(x, |\psi|)] \, dx + \int_\Omega \psi : dE^su \right\} \\ \ge &\sup_{\psi \in T^\varphi} \int_\Omega [\mathcal{E}u : \psi - \varphi^*(x, |\psi|)] \, dx - \sup_{\psi \in T^\varphi} \int_\Omega \psi : dE^su \\ = & \;\varrho_\varphi(|\mathcal{E}u|) - \int_\Omega \varphi^\infty \, d|E^su| = +\infty. \end{align}\] We are only left to consider the case \(\int_\Omega \varphi^\infty \, d|E^su| = +\infty\). By the proof of Proposition 11, see Appendix 5, there exists \(\psi_\varepsilon\in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\) such that \[\int_\Omega \psi_\varepsilon: dE^su > \frac{1}{\varepsilon}\] and \(|\psi_\varepsilon| \le \frac{\varphi(\cdot, k)}{k}\) for some \(k = k_\varepsilon\). For any function \(\theta : \Omega \to [0,1]\), we have that \[\mathcal{E}u : (\theta \psi_\varepsilon) - \varphi^*(\cdot, |\theta \psi_\varepsilon|) \ge \mathcal{E}u : (\theta \psi_\varepsilon) - \varphi^* \left(\cdot, \frac{\varphi(\cdot, k)}{k} \right) \ge - \left(\frac{\varphi^+(k)}{k} \, |\mathcal{E}u| + \varphi^+(k) \right),\] see Lemma 1(iii). Since the function on the right-hand side is integrable, we can choose \(\delta_k > 0\) such that its integral over any measurable \(A\) with \(|A| < \delta_k\) is at least \(-1\). Furthermore, since \({\rm supp}\, E^su\) has measure zero, we can choose \(\theta \in C^\infty_c(\Omega)\) as before to have support with Lebesgue measure at most \(\delta_k\) and such that satisfies \[\int_\Omega (\theta \psi_\varepsilon) : dE^su > \frac{1}{2} \int_\Omega \psi_\varepsilon: dE^su > \frac{1}{2\varepsilon}.\] Then, \[\begin{align} &\sup_{\psi \in T^\varphi} \left\{\int_\Omega [\mathcal{E}u : \psi - \varphi^*(x, |\psi|)] \, dx + \int_\Omega \psi : dE^su \right\} \\ \ge & \int_\Omega (\theta \psi_\varepsilon) : dE^su + \int_\Omega [\mathcal{E}u : (\theta \psi_\varepsilon) - \varphi^*(x, |\theta \psi_\varepsilon|)] \, dx \ge \frac{1}{2\varepsilon} - 1. \end{align}\] Letting \(\varepsilon\to +\infty\), the claim follows. This concludes the proof of the theorem. ◻

3.3 The variable exponent case↩︎

In this section we focus on the special \(\Phi\)-function \[\label{nostraphi} \varphi(x,t) := \frac{1}{p(x)} \, t^{p(x)},\tag{16}\] with \(p:\Omega \to [1,+\infty)\) a measurable function. Such a \(p\) is called a variable exponent. Let us recall that, choosing \(\varphi\) as in 16 and considering the set \[\label{set95superlin} Y := \{x \in \Omega \ : \ p(x) = 1\},\tag{17}\] we have that the conjugate \(\varphi^*\) is of the form \[\label{phi42def} \varphi^*(x,t) = \begin{dcases} \tfrac{1}{q(x)} \, t^{q(x)} & \textrm{ if } x \in \Omega \setminus Y, \\ \infty \chi_{(1, \infty)}(t) & \textrm{ if } x \in Y, \end{dcases}\tag{18}\] see [15], where \(q\) is the function defined pointwise in \(\Omega \setminus Y\) by \(\tfrac{1}{q(x)} + \tfrac{1}{p(x)} = 1\) and with the convention that \(\infty \cdot 0 = 0\). From now, given a set \(A \subset \Omega\), we write \[p^+_A := \underset{x \in A}{\textrm{ess-sup} \;} p(x) \qquad \textrm{ and } \qquad p^-_A := \underset{x \in A}{\textrm{ess-inf} \;} p(x),\] respectively \(q^\pm_A\) for the conjugate exponent function \(q\). In the following, we will assume \(p^+_\Omega < +\infty\). We start by recalling the following notion, introduced in [26].

Definition 11. A variable exponent \(p\) is strongly \(\log\)-Hölder continuous in \(Y\)* if \[p(x)-1 \le \frac{\omega(|x-y|)}{\log(e+ \frac{1}{|x - y|})},\] for some \(\omega:[0,\infty)\to [0, \infty)\) with \(\lim_{t \to 0}\omega(t)=0\) and every \(y \in Y\) and \(x\in\Omega\).*

At this point, it is worth observing that the above property is weaker than the ‘classical \(\log\)-Hölder continuity’ as stated in [26], which is, in turn, equivalent to the restricted (VA1) condition, introduced in Definition 9, for the function 16 .
In [26] the following technical result has been proven.

Lemma 8. Suppose that \(p\) is strongly \(\log\)-Hölder continuous in \(Y\). If \(\rho_{q(\cdot)}(w)<\infty\) for some continuous \(w\), then \(|w|\le 1\) in \(Y\).

Remark 13. The explicit definition of the function \(\varphi\) in 16 entails some additional properties other than the explicit form of its conjugate function given by 18 . Indeed, it is immediate to see that \(\varphi\in \Phi_c(\Omega)\), see Definition 5 and [11]. Moreover, arguing as in [11], it is easy to see that \(\varphi\) satisfies (A0), (aInc)\(_{p^-_\Omega}\) and (aDec)\(_{p^+_\Omega}\), see Definition 7.

The first result we prove in this subsection is the equivalent of Theorem 10 in this specific setting. In particular, we show that the same conclusion of Theorem 10 holds under less restrictive hypotheses on our function \(\varphi\). Before doing so, we need to recall the definition of the modular associated to the variable exponent space when \(p^+_\Omega < +\infty\), given in an equivalent way to the one in [12], namely \[\varrho_{L^{p(\cdot)}(\Omega; \mathbb{R}^{n \times n}_{\rm sym})}(u) := \int_{\Omega} \frac{1}{p(x)} |u(x)|^{p(x)} \, dx\] for any \(u \in \mathcal{M}(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\). Now, the theorem reads as follows.

Theorem 14. Let \(\Omega \subset \mathbb{R}^n\) be a bounded, connected and open set. Let \(\varphi\) be defined as in 16 and let \(p : \Omega \to [1, +\infty)\) be continuous in \(\Omega\) and strongly \(\log\)-Hölder continuous in \(Y\), with \(p^+_\Omega < +\infty\). If \(u \in BD(\Omega)\), then \[\label{rho95old} \varrho_{\widetilde{V}^{n \times n}_\varphi}(Eu) = \varrho_{\textrm{old}}(Eu) := |Eu|(Y) + \varrho_{L^{p(\cdot)}(\Omega \setminus Y; \mathbb{R}^{n \times n}_{\rm sym})}(\mathcal{E}u).\qquad{(8)}\]

Proof. First, let us observe that the claim of Lemma 7 is automatically satisfied in the case of \(\varphi\) given by 16 with the only requirement of strong \(\log\)-Hölder continuity in \(Y\) and the assumption \(p^+_\Omega\) is finite. Indeed, on the one hand, in view of 18 and Definition 8 (which is a limit in this specific case), it suffices to prove the inequality in the set \(Y\) given by 17 , since in \(\Omega \setminus Y\) the recession function is \(+\infty\). On the other hand though, whenever \(\psi \in C(\Omega)\) with \(\varrho_{\varphi^*}(|\psi|) < +\infty\), then by 18 we have that \(\varphi^*(x, |\psi|) = 0\) for almost every \(x \in Y\), implying that \(|\psi| \le 1=\varphi^\infty(x)\) in \(Y\), in view of Lemma 8, which proves our claim.
As a second step, let us remark that choosing \(\varphi\) as in 16 , with \(p(\cdot)\) continuous, , Proposition 12 is valid. In fact, by assumption \(\varphi\) is continuous, and by Remark 13, we know that (A0) and (aDec) are satisfied. Hence, we can conclude that \[\sup_{\psi \in T^\varphi} \int_\Omega \left[\mathcal{E}u : \psi - \varphi^*(x,|\psi|) \right] \, dx = \frac{1}{p(x)} \int_\Omega |\mathcal{E}u|^{p(x)} \, dx\] for any \(u \in BD(\Omega)\), where \(T^\varphi\) is defined as in ?? .
We are left to prove that Proposition 11 holds true under these hypotheses and choosing \(\varphi\) as in 16 , namely that \[\sup_{\psi \in T^\varphi} \int_\Omega \psi : dE^su = \int_\Omega \varphi^\infty \, d|E^su|.\] Note that the inequality "\(\le\)" is immediately verified thanks to the fact that Lemma 7 holds under our assumptions. So it suffices to prove that \[\int_\Omega \varphi^\infty d |E^s u|\leq \sup_{\psi \in T^\varphi}\int_\Omega \psi : d E^s u.\] Thanks to the continuity assumption on \(p\) in \(\Omega\), we obtain that 23 and 24 hold. In particular, the second one holds for any \(\beta \in [0,1]\), due to the non decreasing monotonicity of \((\cdot)^{p(x)}\). Observe that \(h_k\) is lower semicontinuous by definition (it is a lower semicontinuous envelope, see [41]) and increasing in \(k\). Consequently the same arguments as in the proof of Proposition 11 lead to the desired conclusion. Moreover, let us remark that for every \(u \in BD(\Omega)\) it holds \[\sup_{\psi \in T^\varphi} \int_\Omega \psi : d E^s u = \int_Y d|E^s u| = |E^s u|(Y),\] since \(\varphi^\infty(x) = 1\) in \(Y\) as observed before, which indeed says that \(E^s u\) is concentrated on \(Y\).
Finally, following the same steps of the proof of Theorem 10, since under our hypotheses the equivalent formulations of Lemma 7, Proposition 11 and Proposition 12 hold, we get the thesis. ◻

Remark 15. As a corollary of Theorem 14, under the same hypotheses we get the coincidence of the spaces \(BD^\varphi\) and \(BD \cap LD^{p(\cdot)}\), with \[LD^{p(\cdot)}(U) := \{u \in L^{p(\cdot)}(U; \mathbb{R}^n) \ : \ Eu \in L^{p(\cdot)}(U;\mathbb{R}^{n \times n}_{\rm sym})\},\] when \(\varphi\) is as in 16 . Alternatively, this latter equality can be proven directly arguing as in [26] even removing the continuity hypotesis on \(p(\cdot)\), replacing it by lower semicontinuity .

4 The imaging problem↩︎

In this section we finally deal with the novel regularizer based on a Musielak-Orlicz generalization of classical Total Generalized Variation functionals that we defined in 1 , namely the anisotropic Total Generalized Variation of order \(2\), with weight \(\alpha = (\alpha_1, \alpha_2) \in \mathbb{R}^+ \times \mathbb{R}^+\), given by \[\begin{align} TGV^{\varphi,2}_\alpha(u):= \sup\bigg\{ &\int_\Omega u \,\mathrm{div}^2 \psi \, dx \ : \ \\ & \;\psi \in C^2_c (\Omega; \mathbb{R}^{n \times n}_{\rm sym}), \|\psi\|_{\varphi^*} \le \alpha_1, \|\mathrm{div}\, \psi\|_{\varphi^*} \le \alpha_2 \bigg\}, \end{align}\] where we consider \(u \in L^1_{\rm loc}(\Omega)\).
In Subsection 4.1 we prove some properties of this functional, while in Subsection 4.2 we give a dual characterization of this functional exploiting the space \(BD^\varphi\) we defined in Section 3. Finally, in Subsection 4.3 we show how the anisotropic Total Generalized Variation inherits classical existence and stability properties.

4.1 Basic properties of \(TGV^{\varphi,2}_\alpha(u)\)↩︎

Following the lines of [4], we define the space of functions of bounded variation of order \(2\) with Orlicz growth and weight \(\alpha\) as \[\label{def:BGV} BGV^{\varphi,2}_\alpha(\Omega) := \{ u \in L^\varphi(\Omega) \ : \ TGV^{\varphi,2}_\alpha(u) < +\infty \},\tag{19}\] endowed with the norm \[\|u\|_{BGV^{\varphi,2}_\alpha} := \|u\|_\varphi+ TGV^{\varphi,2}_\alpha(u),\] where the functional \(TGV^{\varphi,2}_\alpha\) is defined as in 1 .

Remark 16. Note that it would also be possible to define the functionals \(TGV^{\varphi,k}_\alpha\) and the spaces \(BGV^{\varphi,k}_\alpha(\Omega)\) for general \(k \in \mathbb{N}^+\) following the lines of [4]. Moreover, for \(k = 1\) and \(\alpha \in \mathbb{R}^+\), and chosen \(u \in L^1_{\rm loc}(\Omega)\), we would have \[TGV^{\varphi,1}_\alpha(u) = \sup \left\{\int_\Omega u \, \mathrm{div}\psi \, dx \ : \ \psi \in C^1_c (\Omega; \mathbb{R}^n), \|\psi\|_{\varphi^*} \le \alpha \right\} = \alpha V_\varphi(u),\] where the right-hand side is defined as in 2 . Hence, we can indeed speak of a generalization of an anisotropic total variation. Contrary to Remark 4, which was true for \(u \in BV(\Omega)\) in order to know at least that \(Du \in \mathcal{M}(\Omega;\mathbb{R}^n)\), here it is enough to take \(u \in L^1_{{\rm loc}}(\Omega)\).

In the following result we prove some basic properties of the anisotropic Total Generalized Variation.

Proposition 17. Let \(\Omega\) be a bounded, connected and open set in \(\mathbb{R}^n\), with \(n \ge 2\). Let \(\varphi\in \Phi_w(\Omega)\) satisfy (A0) and let \(\alpha = (\alpha_1, \alpha_2) \in \mathbb{R}^+ \times \mathbb{R}^+\). Then, the following statements hold:

  • \(TGV^{\varphi,2}_\alpha\) is a seminorm on the normed space \(BGV^{\varphi,2}_\alpha(\Omega)\).

  • For \(u \in L^1_{\rm loc}(\Omega)\), there holds \(TGV^{\varphi,2}_\alpha(u) = 0\) if and only if \(u\) is a polynomial of degree \(0\) or \(1\).

  • For eevery \(\tilde{\alpha} = (\tilde{\alpha}_1, \tilde{\alpha}_2) \in \mathbb{R}^+ \times \mathbb{R}^+\), the seminorms \(TGV^{\varphi,2}_\alpha\) and \(TGV^{\varphi,2}_{\tilde{\alpha}}\) are equivalent.

  • \(TGV^{\varphi,2}_\alpha\) is rotationally invariant, namely for any matrix \(R \in SO(n)\) and \(u \in BGV^{\varphi,2}_\alpha(\Omega)\) we can define \(\tilde{u}(x) := u(Rx)\) and we have that \(\tilde{u} \in BGV^{\varphi,2}_\alpha(R^\textrm{T}\Omega)\) and \(TGV^{\varphi,2}_\alpha(\tilde{u}) = TGV^{\varphi,2}_\alpha(u)\).

  • For \(r > 0\) and \(u \in BGV^{\varphi,2}_\alpha(\Omega)\), define \(\bar{u}(x) := u(rx)\). Then \(\bar{u} \in BGV^{\varphi,2}_\alpha(r^{-1}\Omega)\) and \[TGV^{\varphi,2}_\alpha(\bar{u}) = r^{-n} \, TGV^{\varphi,2}_{\bar\alpha}(u) \qquad \textrm{ with } \bar\alpha = (\bar\alpha_1, \bar\alpha_2) := (\alpha_1 r^2, \alpha_2 r).\]

Proof. (i). Thanks to the definition of the functional \(TGV^{\varphi,2}_\alpha\) given in 1 , it is obvious to see that the two conditions (N2) and (N3) in Definition 1 hold, implying that \(TGV^{\varphi,2}_\alpha\) is a seminorm. Note that \(BGV^{\varphi,2}_\alpha(\Omega)\) is a normed space because \(\varphi\in \Phi_w(\Omega)\) and satisfies (A0), so \(\|\cdot\|_\varphi\) is a norm, see [11].
(ii). Let us denote the set of the test functions for the functional \(TGV^{\varphi,2}_\alpha\) by \[K^{\varphi,2}_\alpha(\Omega) := \{ \psi \in C^2_c (\Omega; \mathbb{R}^{n \times n}_{\rm sym}) \ : \ \|\psi\|_{\varphi^*} \le \alpha_1, \|\mathrm{div}\, \psi\|_{\varphi^*} \le \alpha_2 \}.\] Let us first assume that \(u\) is a polynomial of degree \(0\) or \(1\), implying that \(\nabla^2 u = 0\). Then, since the test function \(\psi\) is in \(K^{\varphi,2}_\alpha(\Omega)\) and thanks to two integration by parts, it is obvious that \(TGV^{\varphi,2}_\alpha(u) = 0\).
Suppose now that \(TGV^{\varphi,2}_\alpha(u) = 0\) for some \(u \in L^1_{\rm loc}(\Omega)\). For any \(\psi \in C^2_c (\Omega; \mathbb{R}^{n \times n}_{\rm sym})\), one can find a \(\lambda > 0\) such that \[\|\lambda \psi\|_\infty \le \alpha_1 \quad \textrm{ and } \quad \|\mathrm{div}(\lambda \psi)\|_\infty \le \alpha_2.\] Since \(\varphi\) satisfies (A0), meaning that \(\varphi^*\) satisfies (A0) as well, see Lemma 1(iv), and since \(\Omega\) is bounded, thanks to Lemma 2 we have that \(\lambda\psi \in K^{\varphi,2}_\alpha(\Omega)\). Testing in 1 with \(\lambda\psi\), we then get \[\int_\Omega u \,\mathrm{div}^2 \psi \, dx = 0 \qquad \textrm{ for any } \psi \in C^2_c (\Omega; \mathbb{R}^{n \times n}_{\rm sym}),\] implying that \(\nabla^2 u = 0\) in the weak sense. Since \(\Omega\) is connected, this implies that \(u\) has to be a polynomial of degree less than \(2\).
(iii). The equivalence of the seminorms \(TGV^{\varphi,2}_\alpha\) and \(TGV^{\varphi,2}_{\tilde{\alpha}}\) can be proven noticing that, chosen \[c := \frac{\max\{\alpha_1, \alpha_2\}}{\min\{\tilde{\alpha}_1, \tilde{\alpha}_2\}},\] we have that \(c \, K^{\varphi,2}_\alpha(\Omega) \subset K^{\varphi,2}_{\tilde{\alpha}}(\Omega)\), meaning that \(c \, TGV^{\varphi,2}_\alpha(u) \le TGV^{\varphi,2}_{\tilde{\alpha}}(u)\) for any \(u \in L^1_{\rm loc}(\Omega)\). Interchanging the roles of \(\alpha\) and \(\tilde{\alpha}\) gives us the desired equivalence.
(iv). Thanks to [4], we know that given a matrix \(R \in SO(n)\), then \[\psi \in C^2_c (R^\textrm{T}\Omega; \mathbb{R}^{n \times n}_{\rm sym}) \quad \iff \quad \widetilde{\psi} := (\psi \circ R^\textrm{T})R^\textrm{T}\in C^2_c (\Omega; \mathbb{R}^{n \times n}_{\rm sym}),\] with \(\mathrm{div}^l \, \widetilde{\psi} = ((\mathrm{div}^l \, \psi) \circ R^\textrm{T})R^\textrm{T}\), with \(l = 0,1\) where \(\mathrm{div}^0 = id\). Hence, for \(l = 0,1\), \[\|\mathrm{div}^l \, \widetilde{\psi}\|_{\varphi^*} = \| ((\mathrm{div}^l \, \psi) \circ R^\textrm{T})R^\textrm{T}\|_{\varphi^*} = \| (\mathrm{div}^l \, \psi) \circ R^\textrm{T}\|_{\varphi^*} = \| \mathrm{div}^l \, \psi \|_{\varphi^*},\] implying that \(\psi \in K^{\varphi,2}_\alpha(R^\textrm{T}\Omega)\) if and only if \(\widetilde{\psi} \in K^{\varphi,2}_\alpha(\Omega)\). Lastly, for any \(\psi \in K^{\varphi,2}_\alpha(R^\textrm{T}\Omega)\) we have \[\int_{R^\textrm{T}\Omega} u(Rx) \,\mathrm{div}^2 \psi(x) \, dx = \int_\Omega u(x) \,\mathrm{div}^2 \psi(R^\textrm{T}x) \, dx = \int_\Omega u(x) \, \mathrm{div}^2 \widetilde{\psi}(x) \, dx,\] yielding \(TGV^{\varphi,2}_\alpha(u \circ R) = TGV^{\varphi,2}_\alpha(u)\).
(v). Note that \[\mathrm{div}(v \circ r^{-1} \mathbb{I}_n) = r^{-1} (\mathrm{div}\, v) \circ r^{-1} \mathbb{I}_n,\] implying that, chosen \(\psi \in C^l_c(r^{-1}\Omega;\mathbb{R}^{2 \times 2}_{\rm sym})\) and \(\bar\psi := r^2 \, \psi \circ r^{-1} \mathbb{I}_n\), we have \[\| \mathrm{div}^l \bar\psi \|_{\varphi^*} = r^{2-l} \| \mathrm{div}^l \psi \|_{\varphi^*},\] for \(l = 0,1\). As a consequence, \(\psi \in K^{\varphi,2}_\alpha(r^{-1} \Omega)\) if and only if \(\bar\psi \in K^{\varphi,2}_{\bar\alpha}(\Omega)\) with \(\bar\alpha = (\bar\alpha_1, \bar\alpha_2) := (\alpha_1 r^2, \alpha_2 r)\). Lastly, this yields \[\int_{r^{-1} \Omega} u(rx) \,\mathrm{div}^2 \psi(x) \, dx = r^{-n} \int_\Omega u(x) \,\mathrm{div}^2 \psi(r^{-1} x) \, dx = \int_\Omega u(x) \, \mathrm{div}^2 \bar\psi(x) \, dx,\] implying \(TGV^{\varphi,2}_\alpha(u \circ r \, \mathbb{I}_n) = r^{-n} \, TGV^{\varphi,2}_{\bar\alpha}(u)\). ◻

Note that it is possible to prove that \(TGV^{\varphi,2}_\alpha\) is weakly lower semicontinuous with respect to the weak \(L^p\) convergence, as well as to the weak* convergence in \(L^\infty\).

Lemma 9. Let \(\Omega\) be a bounded set in \(\mathbb{R}^n\) and let \(\varphi\in \Phi_w(\Omega)\). Then, \(TGV^{\varphi,2}_\alpha\) is weakly lower semicontinuous in \(L^p(\Omega)\) for every \(1 \le p < +\infty\), as well as weak* lower semicontinuous in \(L^\infty(\Omega)\) for any \(\alpha \in \mathbb{R}^+ \times \mathbb{R}^+\).

Proof. We only prove the lemma for \(p \in [1, +\infty)\), as the case \(p = +\infty\) is analogous.
Chosen \(\psi \in C^2_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\) with \(\|\psi\|_{\varphi^*} \le \alpha_1\) and \(\|\mathrm{div}\, \psi\|_{\varphi^*} \le \alpha_2\) as in 1 , then \(\mathrm{div}^2 \psi \in L^\infty(\Omega)\). Then, chosen \((u_k)_k \subset L^p(\Omega)\) such that \(u_k \rightharpoonup u\) weakly in \(L^p(\Omega)\) we get \[\int_\Omega u \, \mathrm{div}^2\psi \, dx = \liminf_{k \to +\infty} \int_\Omega u_k \, \mathrm{div}^2\psi \, dx \le \liminf_{k \to +\infty} TGV^{\varphi,2}_\alpha(u_k).\] The thesis follows by taking the supremum over all such \(\psi\). ◻

The next proposition is the equivalent of [4] in our setting, proving that \(BGV^{\varphi,2}_\alpha(\Omega)\) is a Banach space.

Proposition 18. Let \(\Omega\) be a bounded, connected and open set in \(\mathbb{R}^n\), with \(n \ge 2\). Let \(\varphi\in \Phi_w(\Omega)\) satisfy (A0). Then, for any weight \(\alpha = (\alpha_1, \alpha_2) \in \mathbb{R}^+ \times \mathbb{R}^+\) the space \(BGV^{\varphi,2}_\alpha(\Omega)\) is a Banach space if endowed with the norm \(\|\cdot\|_{BGV^{\varphi,2}_\alpha}\).

Proof. Thanks to Proposition 17(iii), we know that all the seminorms \(TGV^{\varphi,2}_\alpha\), respectively norms \(\|\cdot\|_{BGV^{\varphi,2}_\alpha}\), are equivalent for any \(\alpha \in \mathbb{R}^+ \times \mathbb{R}^+\), so we only need to prove that \(BGV^{\varphi,2}_\alpha(\Omega)\) is complete for some \(\alpha \in \mathbb{R}^+ \times \mathbb{R}^+\).
Note that, since \(\Omega\) is bounded and \(\varphi\) satisfies (A0), then by Lemma 2 we have \(L^\varphi(\Omega) \hookrightarrow L^1(\Omega)\). This, together with Lemma 9, implies that \(TGV^{\varphi,2}_\alpha\) is a lower semicontinuous functional with respect to \(L^\varphi(\Omega)\) for any \(\alpha\) and \(\varphi\) chosen as in the hypotheses.
Now let \((u_n)_n\) be a Cauchy sequence in \(BGV^{\varphi,2}_\alpha(\Omega)\). This yields that \((u_n)_n\) is a Cauchy sequence also in \(L^\varphi(\Omega)\), implying the existence of a limit \(u \in L^\varphi(\Omega)\) such that \[TGV^{\varphi,2}_\alpha(u) \le \liminf_{n \to +\infty} TGV^{\varphi,2}_\alpha(u_n),\] thanks to the lower semicontinuity of the functional. As a consequence, \(u \in BGV^{\varphi,2}_\alpha(\Omega)\) and the only thing left to be shown is that \(u\) is the limit in the corresponding norm as well. In order to do so, let us choose \(\varepsilon> 0\) and \(n \in \mathbb{N}\) such that for any \(\bar{n} \ge n\) it holds that \(TGV^{\varphi,2}_\alpha(u_n - u_{\bar{n}}) \le \varepsilon\). Thanks again to the lower semicontinuity of the functional (with respect to the \(L^1\) weak convergence), letting \(\bar{n} \to +\infty\) we then get \[TGV^{\varphi,2}_\alpha(u_n - u) \le \liminf_{\bar{n} \to +\infty} TGV^{\varphi,2}_\alpha(u_n - u_{\bar{n}}) \le \varepsilon,\] implying that \(u_n \to u\) in \(BGV^{\varphi,2}_\alpha(\Omega)\). ◻

4.2 Dual formulation↩︎

The first step in order to treat minimization problems in image denoising with this new regularizer is to prove an equivalent result to [5], namely proving that 1 can be rewritten equivalently as a minimization problem which can be interpreted as an optimal balancing between the first and second derivative of \(u\) in terms of “sparse” penalization, see [5]. In [5] this was done via the Radon norm (i.e. the norm in \(\mathcal{M}\) associated to the Radon measure \(Du\)), while here we need to use the variant of the generalized variation associated to \(\varphi\) defined in 7 . Given the anisotropic nature of such regularizers, we expect them to better capture fine texture details in the images. The result, providing a dual characterization of \(TGV^{\varphi,2}_\alpha\), reads as follows.

Theorem 19. Let \(\Omega \subset \mathbb{R}^n\) be a bounded, connected, \(C^1\) and open set. Let \(\varphi\in \Phi_w(\Omega)\) satisfy (A0) and let \(\alpha = (\alpha_1, \alpha_2) \in \mathbb{R}^+ \times \mathbb{R}^+\). Then for any \(u \in L^1(\Omega)\) we have \[\label{eq:dualTGV} TGV^{\varphi, 2}_\alpha(u) = \min_{w \in BD(\Omega)} \left\{ \alpha_2 \widetilde{V}^n_\varphi(Du - w) + \alpha_1 \widetilde{V}^{n \times n}_\varphi(E w)\right\},\qquad{(9)}\] where \(TGV^{\varphi, 2}_\alpha\) is defined as in 1 and where we have adopted the convention that the above right-hand side is infinite if \(\widetilde{V}^{n \times n}_\varphi(E w) = +\infty\). If, additionally, \(u \in BV^\varphi(\Omega)\), defined as in 3 , then \[\label{eq:dualTGV95phi} TGV^{\varphi, 2}_\alpha(u) = \min_{w \in BD^\varphi(\Omega)} \left\{ \alpha_2 \widetilde{V}^n_\varphi(Du - w) + \alpha_1 \widetilde{V}^{n \times n}_\varphi(E w)\right\},\qquad{(10)}\] where the space \(BD^\varphi(\Omega)\) is defined as in 8 .

Proof. Let us define \(X_2 := C^2_0(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\), \(X_1 := C^1_0(\Omega; \mathbb{R}^n)\) and \(\Lambda := \mathrm{div}\in \mathcal{L}(X_2,X_1)\). We stress that \[\label{chiusure} X_2 = \overline{C^2_c(\Omega;\mathbb{R}^{n \times n}_{\rm sym})}^{\|\cdot\|_{C^2}} \qquad \textrm{ and } \qquad X_1 = \overline{C^1_c(\Omega;\mathbb{R}^n) }^{\|\cdot\|_{C^1}}.\tag{20}\] Chosen \(u \in L^1(\Omega)\), we define \[F_1(M) := I_{\{ \|\cdot\|_{\varphi*} \le \alpha_2\}}(M), \quad \textrm{ for any } M \in X_2\] and \[F_2^u(v) := I_{\{\| \cdot\|_{\varphi^*} \le \alpha_1\}}(v) - \int_\Omega u\,\mathrm{div}v \, dx, \quad \textrm{ for any } v \in X_1.\] Since \(\varphi\) satisfies (A0), then by Lemma 1(iv) also \(\varphi^*\) satisfies (A0) and, by Lemma 2, we have that \(F_1\) and \(F^u_2\) are continuous in the \(C^2\)- and \(C^1\)-topologies, respectively. Then, thanks to a density argument and to 20 , we have that \[TGV^{\varphi, 2}_\alpha(u) = -\inf_{M \in X_2} \left[ F_1(M) + F_2^u(\Lambda M) \right].\] Let us remark that we can write the space \(X_1\) equivalently as \[\label{unione95domini} X_1 = \bigcup_{\lambda \ge 0} \lambda \left(\mathrm{dom}F_2^u - \Lambda \, \mathrm{dom}F_1\right).\tag{21}\] Indeed, if \(z \in \mathrm{dom}F_2^u\) then, trivially, \(z \in X_1\) implies \(z \in \bigcup_{\lambda \ge 0}\lambda\left(\mathrm{dom}F_2^u - \Lambda \, \mathrm{dom}F_1\right)\) and for some \(\lambda \ge 0\) it is always verified that \(\tfrac{z}{\lambda} \in \mathrm{dom}F^u_2\). This implies that \(X_1 \subset \bigcup_{\lambda \ge 0}\lambda\left(\mathrm{dom}F_2^u - \Lambda \, \mathrm{dom}F_1\right)\). On the other hand, the converse is true since \(X_1\) is a closed vector space (and \(\mathrm{dom}F^u_2\) and \(\mathrm{dom}F_1\) are subsets of \(X_1\)).
Now, since \(X_1\) is a closed set, it means that also the right-hand side of 21 is closed, which allows us to use [42], a Fenchel-Rockafellar duality result, getting \[TGV^{\varphi, 2}_\alpha(u) = -\inf_{M \in X_2} \left[F_1(M) + F_2^u(\Lambda M) \right] = \min_{w \in X_1^*} \left[ (F_2^u)^*(w) + F_1^*(-\Lambda^* w) \right].\] Now, exploiting the definition of Fenchel conjugate, we get \[(F_2^u)^*(w) = \alpha_2 \sup_{v \in X_1} \left\{ \langle w , v \rangle_Y + \int_\Omega u \, \mathrm{div}v \, dx \ : \ \|v\|_{\varphi^*} \le 1 \right\},\] for every \(w\in {\mathcal{D}}'(\Omega;\mathbb{R}^n)\), and \[F_1^*(-\Lambda^* w) = \alpha_1 \sup_{M \in X_2} \left\{ \langle -\Lambda^* w, M \rangle_X \ : \ \|M\|_{\varphi^*} \le 1 \right\},\] for every \(w\in {\mathcal{D}}'(\Omega;\mathbb{R}^n)\) with \(-\Lambda^* w \in {\mathcal{D}}'(\Omega;\mathbb{R}^{n \times n})\), where \(\langle\cdot,\cdot\rangle_{X_1}\) and \(\langle\cdot,\cdot\rangle_{X_2}\) denote, respectively, the duality products between \(X_1\) and \(X_1^*\) and between \(X_2\) and \(X_2^*\). Now, we observe that if \(F_1^*(-\Lambda^* w) = +\infty\), then by Lemma 2 we find \(\|-\Lambda^* w\|_\mathcal{M}= +\infty\). Moreover, since \(-\Lambda^* = E\), we infer that for every \(w\in {\mathcal{D}}'(\Omega;\mathbb{R}^n)\) with \(-\Lambda^* w \in {\mathcal{D}}'(\Omega;\mathbb{R}^{n \times n})\) we have \[F_1^*(-\Lambda^* w) = \begin{dcases} \alpha_1 \widetilde{V}_\varphi^{n \times n}(E w) & \textrm{if } E w\in \mathcal{M}(\Omega;\mathbb{R}^{n \times n}) \\ +\infty & \textrm{if } E w \in {\mathcal{D}}'(\Omega;\mathbb{R}^{n\times n}) \setminus \mathcal{M}(\Omega;\mathbb{R}^{n \times n}). \end{dcases}\] Analogously, by observing that, since \(u\in L^1(\Omega)\), \[(F_2^u)^*(w) = \alpha_2 \sup_{v \in X_1} \left\{ \langle Du-w , v\rangle_{X_1} \ : \ \|v\|_{\varphi^*} \le 1 \right\},\] for every \(w\in {\mathcal{D}}'(\Omega;\mathbb{R}^n)\), the same argument as for \(F_1^*\) yields \[(F_2^u)^*(w) = \begin{dcases} \alpha_2 \widetilde{V}_\varphi^n (Du-w) &\textrm{if } Du - w \in \mathcal{M}(\Omega;\mathbb{R}^n) \\ +\infty & \textrm{if } Du-w \in {\mathcal{D}}'(\Omega;\mathbb{R}^n) \setminus \mathcal{M}(\Omega;\mathbb{R}^n), \end{dcases}\] again for every \(w\in {\mathcal{D}}'(\Omega;\mathbb{R}^n)\). Thus, to summarize, by classical properties of \(BD\) and \(BV\), see [31] and [30], we have that the following characterization holds: \[\label{eq95min95inf} TGV^{\varphi, 2}_\alpha(u) = \begin{dcases} \min_{w \in BD(\Omega)} \left\{ \alpha_2 \widetilde{V}_\varphi^n (Du-w) + \alpha_1 \widetilde{V}_\varphi^{n \times n}(E w) \right\} & \textrm{if } u\in BV(\Omega;\mathbb{R}^n), \\ +\infty & \textrm{otherwise in } L^1(\Omega;\mathbb{R}^n). \end{dcases}\tag{22}\] To conclude the proof we observe that if \(u \in BV^\varphi(\Omega)\), then by [11] the above infimum problem can be considered on the subspace of \(BD(\Omega)\) where additionally \(\|w\|_{\varphi^{**}} < +\infty\), namely \(BD^\varphi(\Omega)\), see Remark 8. ◻

Note that in the proof we exploited the condition (A0) on the function \(\varphi\) in order to have the inclusions \(BV^\varphi(\Omega) \hookrightarrow BV(\Omega)\) and \(BD^\varphi(\Omega) \hookrightarrow BD(\Omega)\), see Remarks 2 and 7.

The following result reduces the functional-analytic setting to the space \(BV^\varphi(\Omega)\), see 3 , by establishing equivalence of the two spaces \(BGV^{\varphi,2}_\alpha\) and \(BV^\varphi\).

Corollary 2. Under the same assumptions of Theorem 19, we have that \(TGV^{\varphi,2}_\alpha(u) < +\infty\) if and only if \(u\in BV^\varphi(\Omega)\).

Proof. By Remark 4, choosing \(w = 0\) as a test function in the characterization of \(TGV^{\varphi,2}_\alpha(u)\), namely ?? , we directly infer \[TGV^{\varphi,2}_\alpha(u) \le \alpha_2 \widetilde{V}_\varphi^n(Du) = \alpha_2 V_\varphi(u),\] (see 2 ), which, in turn, yields \[\|u\|_\varphi+ TGV^{\varphi,2}_\alpha(u) \le \max\{1,\alpha_2 \} \|u\|_{BV^\varphi(\Omega)}.\] Conversely, since \(\widetilde{V}^n_\varphi\) satisfies the triangle inequality, letting \(w_u\) be the minimum in ?? , we deduce \[\alpha_2 V_{\varphi}(u) = \alpha_2 \widetilde{V}_\varphi^n(Du) \le \alpha_2 \widetilde{V}_\varphi^n(Du - w_u) + \alpha_2 \widetilde{V}_\varphi^n(w_u) \le TGV^{\varphi,2}_\alpha(u) + \alpha_2 \widetilde{V}_\varphi^n(w_u).\] ◻

Remark 20. Thanks to Proposition 17, in particular thesis (i), and to Corollary 2, we have that \(TGV^{\varphi,2}_\alpha\) is a seminorm also on the space \(BV^\varphi(\Omega)\).

Remark 21. Differently from the classical \(TGV\)-setting, it is not clear whether it is possible to show that there exists a constant \(C > 0\) such that \[\frac{1}{C} \|u\|_{BV^\varphi(\Omega)} \le \|u\|_\varphi+ TGV^{\varphi,2}_\alpha(u)\] or not. Roughly speaking, this is due to the absence of a Musielak-Orlicz equivalent of Sobolev-Korn inequality. We refer the interested reader to the discussion and the references in Section 3.

We also stress that, as it follows from the analysis in Proposition 22, \[\frac{1}{C} (\|u\|_{\varphi}+ |Du|(\Omega))\leq \|u\|_{\varphi}+ TGV^{\varphi,2}_\alpha(u)\] for every \(u \in BV^\varphi(\Omega)\).

Thanks to what we have proven in Subsections 3.2 and 3.3, we expect also the variation \(\widetilde{V}^m_\varphi\) and, as a consequence, the total generalized variation \(TGV^{\varphi,2}_\alpha\) to have a decomposition depending on the absolutely continuous part and the singular part of the gradient of the functions \(u \in BV^\varphi(\Omega)\). We leave this open question for future works.

4.3 Existence and stability of solutions to the minimization problem↩︎

In this section we show how our notion of total generalized variation inherits classical existence and stability properties. Throughout this section we will always tacitly assume that \(TGV^{\varphi,2}_\alpha\) is extended to \(+\infty\) outside of its domain. We first prove that a Poincaré-Wirtinger inequality also holds in our setting, entailing the coercivity needed to prove the existence of solutions to minimum problems in image denoising.
To this aim, from now on we denote by \(\mathcal{P}^1(\Omega)\) the space of affine functions, as well as by \(\Pi:L^1(\Omega)\to \mathcal{P}^1(\Omega)\) a linear projection.

Proposition 22. Let \(\Omega \subset \mathbb{R}^n\) be a bounded, connected, open set with \(C^1\) boundary. Let \(\varphi\in \Phi_w(\Omega)\) satisfy (A0) and let \(\alpha = (\alpha_1, \alpha_2) \in \mathbb{R}^+ \times \mathbb{R}^+\). Let \(1 < p \le \tfrac{n}{n-1}\). Then, there exists a constant \(c > 0\) such that \[\|u\|_{L^p(\Omega)} \le c \,TGV^{\varphi,2}_\alpha(u),\] for every \(u \in {\rm ker}\,\Pi \cap L^p(\Omega)\).

Proof. Because \(\varphi\) satisfies (A0), by Remark 2 we have \(BV^\varphi(\Omega) \hookrightarrow BV(\Omega)\) and so \[|Du - w|(\Omega) \le \widetilde{V}_\varphi^n(Du - w).\] Analogously, thanks to Remark 7 we get \(BD^\varphi(\Omega) \hookrightarrow BD(\Omega)\), entailing \[|Ew|(\Omega) \le \widetilde{V}_\varphi^{n\times n}(Ew),\] for every \(w \in BD(\Omega)\). In view of [5], we infer the existence of a constant \(c > 0\) such that \[\begin{align} \|u\|_{L^p(\Omega)} &\le c \, TGV^{2}_\alpha(u) \le c \left[ \alpha_2 |Du - w|(\Omega) + \alpha_1 |Ew|(\Omega) \right] \\ &\le c \left[ \alpha_2 \widetilde{V}_\varphi^n(Du - w) + \alpha_1 \widetilde{V}_\varphi^{n \times n}(Ew) \right], \end{align}\] for every \(w \in BD(\Omega)\). The thesis follows then by Theorem 19. ◻

In view of the Lemma 9, we deduce existence of minimizers whenever \(TGV^{\varphi,2}_\alpha\) is augmented by suitable fidelity terms. In what follows, given two Hilbert spaces \(X\) and \(H\), we denote by \(\mathcal{L}(X;H)\) the set of linear and continuous operators between \(X\) and \(H\).

Theorem 23. Let \(\Omega \subset \mathbb{R}^n\) be a bounded, connected, open set with \(C^1\) boundary. Let \(\varphi\in \Phi_w(\Omega)\) satisfy (A0) and let \(\alpha = (\alpha_1, \alpha_2) \in \mathbb{R}^+ \times \mathbb{R}^+\). Let \(1 < p \le \tfrac{n}{n-1}\). Let further \(H\) be an Hilbert space, and let \(K \in \mathcal{L}(L^p(\Omega);H)\) be injective on \(\mathcal{P}^1(\Omega)\). Then, for every \(f \in H\), the minimum problem \[\label{eq:minK} \min_{u \in L^p(\Omega)} \frac{1}{2} \|Ku-f\|_H + TGV^{\varphi,2}_\alpha(u)\qquad{(11)}\] has a solution.

Proof. Arguing exactly as in the proof of [5], in view of Propositions 17 and 22, we find that any minimizing sequence for ?? is uniformly bounded in \(L^p(\Omega)\) for every \(1\leq p \le \tfrac{n}{n-1}\). The thesis follows then by Lemma 9. ◻

We conclude this section by showing stability of minimizers of ?? with respect to the fidelity datum \(f\).

Theorem 24. Let \(\Omega \subset \mathbb{R}^n\) be a bounded, connected, open set with \(C^1\) boundary. Let \(\varphi\in \Phi_w(\Omega)\) satisfy (A0) and let \(\alpha = (\alpha_1, \alpha_2) \in \mathbb{R}^+ \times \mathbb{R}^+\). Let \(1 < p \le \tfrac{n}{n-1}\). Let \(H\) be an Hilbert space and let \(K \in \mathcal{L}(L^p(\Omega);H)\) be injective on \(\mathcal{P}^1(\Omega)\). Let \((f_k)_k \subset H\) be such that \(f_k \to f\) strongly in \(H\). Then, setting \[\mathcal{F}_k(u) := \frac{1}{2} \|Ku - f_k\|_H + TGV^{\varphi,2}_\alpha(u), \qquad \mathcal{F}_0(u) := \frac{1}{2} \|Ku - f\|_H + TGV^{\varphi,2}_\alpha(u),\] for every \(u\in L^p(\Omega)\), we have that \((\mathcal{F}_k)_k\) \(\Gamma\)-converges to \(\mathcal{F}_0\) with respect to the weak \(L^p\)-topology.

Proof. The limsup inequality is immediately verified by considering for each \(u\in L^p(\Omega)\) the constant recovery sequence \(\bar{u}_k = u\) for every \(k \in \mathbb{N}\). To show the liminf inequality, let \(u \in L^p(\Omega)\) and let \((u_k)_k \subset L^p(\Omega)\) be such that \(u_k \rightharpoonup u\) weakly in \(L^p(\Omega)\). The fact that \[\mathcal{F}_0(u) \le \liminf_{k \to +\infty} \mathcal{F}_k(u_k)\] is a direct consequence of Lemma 9. ◻

As a consequence of Theorem 24 we infer convergence of the associated minimizers.

Corollary 3. Under the same hypotheses of Theorem 24, we have

  • \(\min \mathcal{F}_k \to \min \mathcal{F}_0\).

  • Letting \(u_k \in \textrm{argmin }\mathcal{F}_k\) for every \(k \in \mathbb{N}\), we have that \((u_k)_k\) is precompact in \(L^{\frac{n}{n-1}}(\Omega)\), and that every limiting point \(u\) satisfies \(u \in \textrm{argmin }\mathcal{F}_0\).

Proof. The result follows from the Fundamental Theorem of Gamma-convergence, cf. [43], by combining Theorem 24 and Proposition 22. ◻

5 Proofs of the results in Subsection 3.2↩︎

In Proposition 11 and Proposition 12 we consider respectively the singular and absolutely continuous parts of the symmetric gradient separately and we then combine them to handle the whole function in Theorem 10.
In order to prove Proposition 11, we need the following result, analogous to [15], which holds a characterization of the singular part of the symmetric gradient of our function.

Lemma 10. Let \(\Omega \subset \mathbb{R}^n\) be a bounded, connected and open set. If \(u \in BD(\Omega)\), then \[|E^s u|(\Omega) = \sup \left\{ \int_\Omega \psi : dE^su \ : \ \psi \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym}), \, |\psi| \le 1 \right\}.\]

Proof. By the definitions of total variation of a measure and of weak derivative, together with the fact that \(Eu = \mathcal{E}u + E^s u\), we see that \[\begin{align} |Eu|(\Omega) &= \sup_{\psi \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym}), |\psi| \le 1} \left( \int_\Omega \psi : d\mathcal{E}u + \int_\Omega \psi : dE^su \right) \\ &\le \sup_{\psi \in L^\infty(\Omega; \mathbb{R}^{n \times n}_{\rm sym}), |\psi| \le 1} \int_\Omega \psi : d \mathcal{E}u + \sup_{\psi \in L^\infty(\Omega; \mathbb{R}^{n \times n}_{\rm sym}), |\psi| \le 1} \int_\Omega \psi : dE^su \\ &= |\mathcal{E}u|(\Omega) + |E^su|(\Omega). \end{align}\] Since \(|\mathcal{E}u|(\Omega) + |E^su|(\Omega) = |Eu|(\Omega)\) as \(\mathcal{E}\) and \(E^s\) are mutually singular, each inequality has to be an equality, and so the claim follows. ◻

We can now give the proof of Proposition 11, following the lines of [15].

Proof of Proposition 11. By the definition of the total variation of a measure, Lemma 7 and the definition of \(T^\varphi\), see ?? , we have \[\sup_{\psi \in T^\varphi} \int_\Omega \psi : dE^su \le \sup_{\psi \in T^\varphi} \int_\Omega |\psi| \, d|E^su| \le \int_\Omega \varphi'_\infty \, d|E^su|.\] In order to tackle the opposite inequality, let us define \(h_k : \Omega \to [0,+\infty]\) by \[h_k(x) := \lim_{r \to 0^+} \inf_{y \in B(x,r)} \frac{\varphi(y, k)}{k} \qquad \textrm{ for any } x \in \Omega.\] Notice that \(h_k\) is lower semicontinuous with \[\label{leqhk} h_k \le \tfrac{\varphi(\cdot, k)}{k} \le \varphi^\infty.\tag{23}\] From the first inequality it follows that \(\varphi^*(\cdot, h_k) \le \varphi(\cdot, k)\), implying that \(\varrho_{\varphi^*}(h_k) \le \varrho_\varphi(k) < +\infty\) since \(\varphi\) satisfies (A0) and (aDec), and \(\Omega\) is bounded.
Let us now remark that \(h_k \to \varphi^\infty\) as \(k \to +\infty\). Indeed, if \(\varphi^\infty(x) = +\infty\), then we can use (A1) in all sufficiently small balls, since \(\varphi^+(k) < +\infty\), to conclude that \[\label{geqhk} h_k(x) = \lim_{r \to 0^+} \inf_{y \in B(x,r)} \frac{\varphi(y, k)}{k} \ge \frac{\varphi(x, \beta k) - 1}{k} \to \beta \varphi^\infty(x) = +\infty\tag{24}\] as \(k \to +\infty\). If instead \(\varphi^\infty(x)< +\infty\), then we use the same inequality but now with \(\beta := \tfrac{1}{1 + \omega(r)}\), where \(\omega\) is the modulus of continuity coming from the restricted (VA1) condition, and we obtain the desired convergence as \(\beta \to 1^-\).
Note that \(h_k\) is increasing in \(k\) since \(\varphi\) is convex. Then, by monotone convergence, \[\int_\Omega \varphi^\infty \, d|E^su| = \lim_{k \to \infty} \int_\Omega h_k \, d|E^su|.\] Let us fix \(\varepsilon> 0\) and assume that \(\int_\Omega \varphi^\infty \, d|E^su| < +\infty\). Then, we can find \(h = h_k\) and \(K > 0\) such that \[\begin{align} \int_\Omega \varphi^\infty \, d|E^su| &\le \int_\Omega h \, d|E^su| + \varepsilon \le \sum_{j = 1}^{K^2} \int_\Omega \frac{1}{K} \, \chi_{\{h > \tfrac{j}{K}\}} \, d|E^su| + 2\varepsilon\\ &= \frac{1}{K} \sum_{j = 1}^{K^2} |E^su|\left(\left\{h > \frac{j}{K}\right\}\right) + 2\varepsilon. \end{align}\] Since \(h\) is lower semicontinuous, then \(\{h > \tfrac{j}{K}\}\) is open, and hence by Lemma 10 we can choose \(\psi_j \in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\) with \(|\psi_j| \le 1\) such that \[\int_\Omega \varphi^\infty \, d|E^su| \le \frac{1}{K} \sum_{j = 1}^{K^2} \int_{\{h > \tfrac{j}{K}\}} \psi_j : dE^su + 3\varepsilon = \int_\Omega \left(\sum_{j = 1}^{K^2} \frac{1}{K} \psi_j \right) : dE^su + 3\varepsilon.\] Defining \(\psi_\varepsilon:= \sum_{j = 1}^{K^2} \tfrac{1}{K} \psi_j\), we can notice that \(\psi_\varepsilon\in C^1_c(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\) and \[|\psi_\varepsilon| \le \frac{1}{K} \sum_{j = 1}^{K^2} |\psi_j| \le \frac{1}{K} \sum_{j = 1}^{K^2} \chi_{\{h > \frac{j}{K}\}} \le h,\] giving \(\varrho_{\varphi^*}(|\psi_\varepsilon|) < +\infty\). Therefore \(\psi_\varepsilon\in T^\varphi\) and \[\int_\Omega \varphi^\infty \, d|E^su| \le \int_\Omega \psi_\varepsilon: dE^su + 3\varepsilon \le \sup_{\psi \in T^\varphi} \int_\Omega \psi : dE^su +3\varepsilon.\] The upper bound then follows letting \(\varepsilon\to 0^+\). If instead \(\int_\Omega \varphi^\infty \, d|E^su| = +\infty\), then a similar argument gives \(\tfrac{1}{3\varepsilon} \le \sup_{\psi \in T^\varphi} \int_\Omega \psi : dE^su\) and the claim again follows. ◻

Lastly, we now give the proof of Proposition 12, following the lines of [15].

Proof of Proposition 12. The upper bound follows directly from Young’s inequality, indeed \[\sup_{\psi \in T^\varphi} \int_\Omega [\mathcal{E}u : \psi - \varphi^*(x, |\psi|)] \, dx \le \int_\Omega \varphi(x, |\mathcal{E}u|) \, dx = \varrho_\varphi(|\mathcal{E}u|).\] In order to prove the lower bound let us choose \(g_i \in C(\Omega; \mathbb{R}^{n \times n}_{\rm sym}) \cap L^\varphi(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\) such that \(g_i \to \mathcal{E}u\) pointwise a.e and in \(L^1(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\). Then Fatou’s Lemma and \(L^1\)-convergence yield \[\int_\Omega \varphi(x, |\mathcal{E}u|) \, dx \le \liminf_{i \to +\infty} \int_\Omega \varphi(x, |g_i|) \, dx \quad \textrm{ and } \quad \lim_{i \to +\infty} \int_\Omega g_i : \psi \, dx = \int_\Omega \mathcal{E}u : \psi \, dx\] for any \(\psi \in T^\varphi\). Thus, it suffices to show that \[\int_\Omega \varphi(x, |g|) \, dx \le \sup_{\psi \in T^\varphi} \int_\Omega [g : \psi - \varphi^*(x, |\psi|)] \, dx\] for any \(g \in C(\Omega; \mathbb{R}^{n \times n}_{\rm sym}) \cap L^\varphi(\Omega; \mathbb{R}^{n \times n}_{\rm sym})\). Moreover, replacing \(\psi\) by \(\frac{g}{\varepsilon+ |g|} |\psi|\) and letting \(\varepsilon\to 0^+\), we see that \(g : \frac{g}{\varepsilon+ |g|} |\psi| \to |g| \, |w|\) pointwise. Thus, by monotone convergence, the vector-values of \(g\) and \(\psi\) are irrelevant and we only need to show that \[\label{eq:thesis952} \int_\Omega \varphi(x, |g|) \, dx \le \sup_{\psi \in C^1_c(\Omega)} \int_\Omega [|g \psi| - \varphi^*(x, |\psi|)] \, dx\tag{25}\] for any \(g\in C(\Omega) \cap L^\varphi(\Omega)\). We can also exclude test functions such that \(\varrho_{\varphi^*}(\psi) = +\infty\), see Remark 9.
Now let us denote with \(\varphi'\) the left-derivative of \(\varphi\) with respect to the second variable. Then \(\varphi'\) is left-continuous and \(\varphi(x,s) \ge \varphi(x,s_0) + \varphi'(x,s_0)(s-s_0)\) by convexity. For any \(\varepsilon> 0\), let us define \[\label{final95choice} \psi_\varepsilon(x,t) := \int_{-\infty}^{+\infty} \varphi(x, \max\{\tau,0\}) \, \zeta_\varepsilon(t-\tau) \, d\tau = (\varphi*_t \zeta_\varepsilon)(x, t),\tag{26}\] where \(\zeta_\varepsilon\) is a mollifier in \(\mathbb{R}\) with support in \([0,\varepsilon]\) and where \(*_t\) denotes the convolution only in the second variable. Since \(\varphi\) and \(\varphi'\) are increasing in the second variable and left-continuous, we have that \(\psi_\varepsilon\nearrow \varphi\) and \(\psi_\varepsilon' \nearrow \varphi'\) as \(\varepsilon\to 0^+\). Moreover, \(\psi_\varepsilon' = \varphi*_t \zeta_\varepsilon'\) is continuous in both \(x\), since \(\varphi\) is, and \(t\), as a convolution with a smooth function.
Let \(v_i \in C_c(\Omega)\) with \(0 \le v_i \le 1\) and \(v_i \nearrow 1\) as \(i \to +\infty\). By uniform continuity in \({\rm supp}\, v_i\), we can choose \(\delta = \delta_{\varepsilon, i} > 0\) such that, for any \(\varepsilon> 0\) and \(i \in \mathbb{N}\), \[\psi_\varepsilon'(x, |g(x)| \, v_i(x)) - \varepsilon\le \psi_\varepsilon'(y, |g(y)| \, v_i(y))\] for all \(x \in B(y,\delta)\) and \(y \in \Omega\). Then \[\psi_{\varepsilon, i} := \max\{\psi_\varepsilon'(\,\cdot\,, |g| \, v_i) - \varepsilon, 0\} *_x \eta_\delta \le \psi_\varepsilon'(\cdot, |g|) \le \varphi'(\cdot, |g|),\] where \(\eta_\delta\) is a mollifier in \(\mathbb{R}\) with support in \([0,\delta]\) and where \(*_x\) denotes the convolution only in the first variable. We have that \(\psi_{\varepsilon,i} \to \varphi'(\cdot, |g|)\), so we conclude by Fatou’s Lemma that \[\int_\Omega |g| \, \varphi'(x, |g|) \, dx \le \liminf_{i \to +\infty, \, \varepsilon\to 0} \int_\Omega |g \, \psi_{\varepsilon, i}| \, dx.\] Since \(\varphi\) satisfies (A0) and (aDec), we see that \[\varphi^*(x, |\psi_{\varepsilon,i}|) \le \varphi^*(x, \varphi'(x, |g|)) \le |g| \, \varphi'(x, |g|) \le c \, \varphi(x, |g|),\] where \(c > 0\) is a constant. Since \(g \in L^\varphi(\Omega)\) and \(\varphi\) satisfies (aDec), then \(\varrho_\varphi(g) < +\infty\). Thus, dominated convergence with upper bound \(c \, \varphi(\cdot, g)\) yields \[\int_\Omega \varphi^*(x, \varphi'(x, |g|)) \, dx = \lim_{i \to +\infty, \, \varepsilon\to 0} \int_\Omega \varphi^*(x, |\psi_{\varepsilon,i}|) \, dx.\] Since \(\psi_\varepsilon\) is a valid test function and \(\varphi'(\cdot, |g|) < +\infty\) a.e. in \(\Omega\), the last equality, together with “Young’s equality”, see [44], implies that \[\begin{align} \sup_{\psi \in C^1_c(\Omega), \, \varrho_{\varphi^*}(\psi) < +\infty} &\int_\Omega [|g \psi| - \varphi^*(x, |\psi|)] \, dx \ge \liminf_{i \to +\infty, \, \varepsilon\to 0} \int_\Omega [|g| \, |\psi_{\varepsilon,i}| - \varphi^*(x, |\psi_{\varepsilon,i}|)]\, dx \\ \ge &\int_\Omega [|g| \, \varphi'(x, |g|) - \varphi^*(x, \varphi'(x, |g|))] \, dx = \int_\Omega \varphi(x, |g|) \, dx, \end{align}\] proving 25 and completing the proof of the lower bound. ◻

Remark 25. The natural idea would be to choose \(\psi := \varphi'(x, |g|)\) in the supremum in 25 , instead of the family of test functions that we chose in 26 , where \(\varphi'\) is the left-derivative of \(\varphi\) with respect to the second variable. However, this function is not in general smooth and we cannot use regular approximation by smooth functions since \(\varphi^*\) is not doubling.

6 A Convex Analysis lemma↩︎

The following result is presented for readers convenience and it is due to [45], see Lemma 2.33 therein. It has been used in Section 3.

Lemma 11. Let \(X\) be a Banach space with norm \(\|\cdot \|_X\), let \(\phi : \mathbb{R}\to \mathbb{R}\) be an even and convex function, and define \(f : X \to \mathbb{R}\) as \(f(x) := \phi(\|x\|_X)\). Then \[f^*(x^*) = \phi^*(\|x^*\|_{X^*})\] where \(X^*\) denotes the dual of \(X\) with norm \(\|\cdot\|_{X^*}\).

Proof. By definition \[f^*(x^*) = \sup_{x \in X} [\langle x^*, x \rangle - f(x)] = \sup_{x \in X} [\langle x^*, x \rangle - \phi(\|x\|_{X})],\] which can be equivalently written as \[\begin{align} \label{24613Amar} f^*(x^*) &= \sup_{t \ge 0} \sup_{\|x\|_X = t} [\langle x^*, x \rangle - \phi(t)] \\ &= \sup_{t \ge 0} \left\{ \left[ \sup_{\|x\|_X = t} \langle x^*, x \rangle \right] - \phi(t) \right\} = \sup_{t \ge 0} [t \|x^*\|_{X^*} - \phi(t)]. \end{align}\tag{27}\] On the other hand, by definition \[\label{24614Amar} \phi^*(\|x^*\|_{X^*}) = \sup_{t \in \mathbb{R}} [t \|x^*\|_{X^*} - \phi(t)],\tag{28}\] hence it suffices to prove that 27 coincides with 28 , i.e. that it suffices to take the supremum on nonnegative real numbers. Using the fact that \(\phi\) is even, for any \(t > 0\) we have \[-t \|x^*\|_{X^*} - \phi(-t) = -t \|x^*\|_{X^*} - \phi(t) \le t \|x^*\|_{X^*} - \phi(t),\] hence \[\sup_{t \in \mathbb{R}} [t\|x^*\|_{X^*} - \phi(t)] = \sup_{t \ge 0} [t\|x^*\|_{X^*} - \phi(t)],\] which concludes the proof. ◻

Acknowledgements↩︎

The research of E. Davoli and S. Riccó was funded by the Austrian Science Fund (FWF) projects 10.55776/F65, 10.55776/Y1292, 10.55776/P35359, and 10.55776/F100800. For open-access purposes, the authors have applied a CC BY public copyright license to any author-accepted manuscript version arising from this submission. G. Bertazzoni and E. Zappale are members of INdAM GNAMPA, whose support is gratefully acknowledged. They are also very grateful to the Institute of Analysis and Scientific Computing at the TU Wien for its kind support and hospitality during the preparation of this work. E. Zappale’s research has also been supported by PRIN 2022: Mathematical Modelling of Heterogeneous Systems (MMHS) - Next Generation EU CUP B53D23009360006.

References↩︎

[1]
L. I. Rudin, S. Osher, E. Fatemi:Nonlinear total variation based noise removal algorithms, Phys. D: Nonlinear Phenomena 60, 259–268, (1992).
[2]
A. Chambolle, P.L. Lions:Image recovery via total variation minimization and related problems, Numer. Math. 76, 167–188, (1997).
[3]
K. Jalalzai:Some remarks on the staircasing phenomenon in total variation-based image denoising, J. Math. Imaging Vis. 54, 256–268, (2016).
[4]
K. Bredies, K. Kunisch, T. Pock:Total generalized variation, SIAM J. Imaging Sci. 3(3), 492–526, (2010).
[5]
K. Bredies, T. Valkonen:Inverse problems with second-order total generalized variation constraints, Proceedings of SampTA 2011 - 9th International Conference on Sampling Theory and Applications, Singapore, (2011).
[6]
P. Blomgren, T. F. Chan, P. Mulet, C. Wong:Total variation image restoration: Numerical methods and extensions, in Proceedings of the IEEE International Conference on Image Processing, Vol. III, IEEE, Los Alamitos, CA, 384–387, (1997).
[7]
D. M. Strong, T. F. Chan:Spatially and scale adaptive total variation based regularization and anisotropic diffusion in image processing, Technical Report CAM96-46, University of California, Los Angeles, CA, (1996).
[8]
Y. Chen, S. Levine, M. Rao:Variable exponent, linear growth functionals in image restoration, SIAM J. Appl. Math. 66(4), 1383–1406, (2006).
[9]
I. Chlebicka, P. Gwiazda, A. Świerczewska-Gwiazda, A. Wróblewska-Kamińska:Partial Differential Equations in Anisotropic Musielak–Orlicz Spaces, Springer, Cham, (2021).
[10]
L. Diening, P. Harjulehto, P. Hästö, M. Růžička:Lebesgue and Sobolev spaces with variable exponents, in: Lecture Notes in Mathematics, vol. 2017, Springer, Heidelberg, (2011).
[11]
P. Harjulehto, P. Hästö:Orlicz spaces and generalized Orlicz spaces, Lecture Notes in Mathematics, vol. 2236, Springer, (2019).
[12]
P. Harjulehto, P. Hästö, V. Latvala:Minimizers of the variable exponent, non-uniformly convex Dirichlet energy, J. Math. Pures Appl. 89, 174–197, (2008).
[13]
P. Harjulehto, P. Hästö, V. Latvala, O. Toivanen:Critical variable exponent functionals in image restoration, Appl. Math. Lett. 26, 56–60, (2013).
[14]
M. M. Rao, Z. D. Ren:Theory of Orlicz spaces, CRC Press, (1991).
[15]
M. Eleuteri, P. Harjulehto, P. Hästö:Bounded variation spaces with generalized Orlicz growth related to image denoising, Mathematische Zeitschrift 310, 26, (2025).
[16]
W. Górny, M. Łasica, A. Matsoukas:Euler-Lagrange equations for variable-growth total variation, Nonlinear Anal. TMA 263, 113984, (2026).
[17]
W. Górny, M. Łasica, A. Matsoukas:Adaptive double-phase Rudin–Osher–Fatemi denoising model, https://arxiv.org/abs/2510.04382, (2025).
[18]
J. S. Moll, V. Pallardó-Julià:Anisotropic Chan–Vese segmentation, Nonlinear Analysis: Real World Applications 73, 103908, (2023).
[19]
H. Chen, C. Wang, Y. Song, Z. Li:Split Bregmanized anisotropic total variation model for image deblurring, J. Vis. Commun. Image Represent. 31, 282–293 (2015).
[20]
S. Esedo\(\bar{\textsc{g}}\)lu, S. J. Osher:Decomposition of images by the anisotropic Rudin-Osher-Fatemi model, Commun. Pure Appl. Math. 57, 1609–1626, (2004).
[21]
M. Łasica, S. Moll, P.B. Mucha:Total variation denoising in \(l^1\) anisotropy, SIAM J. Imaging Sci. 10, 1691–1723, (2017).
[22]
J. Musielak:Orlicz spaces and modular spaces, Lecture Notes in Mathematics, vol. 1034, Springer, Berlin, (1983).
[23]
P. Hästö, J. Juusti, H. Rafeiro:Riesz spaces with generalized Orlicz growth, https://arxiv.org/abs/2204.14128, (2022).
[24]
P. Hästö:The maximal operator on generalized Orlicz spaces, J. Funct. Anal. 269, 4038–4048, (2015).
[25]
P. Hästö, J. Ok:Maximal regularity for non-autonomous differential equations, J. Eur. Math. Soc. (JEMS) 24(4), 1285–1334, (2022).
[26]
G. Bertazzoni, P. Harjulehto, P. Hästo, E. Zappale:Relaxation of quasi-convex functionals with variable exponent growth, https://arxiv.org/abs/2510.04672, (2025).
[27]
M. Amar, G. Bellettini:Approximation by \(\Gamma\)-convergence of a total variation with discontinuous coefficients, Asympt. Anal. 10(3), 225–243, (1995).
[28]
J. S. Moll:The anisotropic total variation flow, Math. Ann. 332, 177–218, (2005).
[29]
E. Davoli, R. Ferreira, I. Fonseca, J.A. Iglesias:Dyadic partition-based training schemes for TV/TGV denoising, J. Math. Imaging Vis. 66, 1070–1108, (2024).
[30]
L. Ambrosio, N. Fusco, D. Pallara:Functions of bounded variation and free discontinuity problems, Oxford Mathematical Monographs, New York, (2000).
[31]
R. Temam:Mathematical problems in plasticity. Modern Applied Mathematics Series, Gauthier-Villars, (1985).
[32]
J. Deutsch, S. Riccò:A regularity theorem for \(BV^{\mathcal A}(\Omega)\), in preparation.
[33]
L. Ambrosio, A. Coscia, G. Dal Maso:Fine properties of functions with bounded deformation, Arch. Rational Mech. Anal. 139, 201–238, (1997).
[34]
D. Ornstein:A non-inequality for differential operators in the \(L_1\) norm, Arch. Rational Mech. Anal. 11, 40–49, (1962).
[35]
D. Breit, L. Diening:Sharp Conditions for Korn inequalities in Orlicz spaces, J. Math. Fluid Mech. 14, 565–573, (2012).
[36]
D. Breit, A. Cianchi, L. Diening:Trace-free Korn inequalities in Orlicz spaces, SIAM J. Math. Anal. 49(4), 2496–2526, (2017).
[37]
A. Cianchi:Korn type inequalities in Orlicz spaces, J. Funct. Anal. 267, 2313–2352, (2014).
[38]
L. Diening, M. Růžička:Calderon–Zygmund operators on generalized Lebesgue spaces \(L^{p(\cdot)}\) and problems related to fluid dynamics, J. Reine Angew. Math. 563, 197–220, (2003).
[39]
G. Bellettini, A. Coscia, G. Dal Maso:Compactness and lower semicontinuity properties in \(BD(\Omega)\), Math. Z. 228, 337–351, (1998).
[40]
D. Baruah, P. Harjulehto, P. Hästö:Capacities in Generalized Orlicz Spaces, J. Func. Spaces, 8459874, (2018).
[41]
G. Dal Maso:An introduction to \(\Gamma\)-convergence, Progress in Nonlinear Differential Equations and Their Applications, Volume 8, (1993).
[42]
H. Attouch, H. Brezis:Duality for the sum of convex functions in general Banach spaces, Asp. Math. Appl. 34, 125–133, (1986).
[43]
A. Braides:Gamma-convergence for beginners, Oxford Lecture Series in Mathematics and Its Applications, (2002).
[44]
C. P. Niculescu, L.-E. Persson:Convex functions and their applications: A contemporary approach(2nd ed.), CMS Books in Mathematics, Springer, Cham, (2018).
[45]
M. Amar:Semicontinuitá e rilassamento. Appunti del corso di dottorato in Modelli e Metodi Matematici per la tecnologia e la societá, a.a. 2003-2004.