Sobolev space explained

In mathematics, a Sobolev space is a vector space of functions equipped with a norm that is a combination of Lp-norms of the function together with its derivatives up to a given order. The derivatives are understood in a suitable weak sense to make the space complete, i.e. a Banach space. Intuitively, a Sobolev space is a space of functions possessing sufficiently many derivatives for some application domain, such as partial differential equations, and equipped with a norm that measures both the size and regularity of a function.

Sobolev spaces are named after the Russian mathematician Sergei Sobolev. Their importance comes from the fact that weak solutions of some important partial differential equations exist in appropriate Sobolev spaces, even when there are no strong solutions in spaces of continuous functions with the derivatives understood in the classical sense.

Motivation

In this section and throughout the article

\Omega

is an open subset of

\Rn.

There are many criteria for smoothness of mathematical functions. The most basic criterion may be that of continuity. A stronger notion of smoothness is that of differentiability (because functions that are differentiable are also continuous) and a yet stronger notion of smoothness is that the derivative also be continuous (these functions are said to be of class

C1

- see Differentiability classes). Differentiable functions are important in many areas, and in particular for differential equations. In the twentieth century, however, it was observed that the space

C1

(or

C2

, etc.) was not exactly the right space to study solutions of differential equations. The Sobolev spaces are the modern replacement for these spaces in which to look for solutions of partial differential equations.

Quantities or properties of the underlying model of the differential equation are usually expressed in terms of integral norms. A typical example is measuring the energy of a temperature or velocity distribution by an

L2

-norm. It is therefore important to develop a tool for differentiating Lebesgue space functions.

The integration by parts formula yields that for every

u\inCk(\Omega)

, where

k

is a natural number, and for all infinitely differentiable functions with compact support

\varphi\in

infty
C
c

(\Omega),

\int\OmegauD\alpha\varphidx=(-1)|\alpha|\int\Omega\varphiD\alphaudx,

where

\alpha

is a multi-index of order

|\alpha|=k

and we are using the notation:

D\alphaf=

\partial|f
\partial
\alpha1
x
1
...\partial
\alphan
x
n

.

The left-hand side of this equation still makes sense if we only assume

u

to be locally integrable. If there exists a locally integrable function

v

, such that

\int\OmegauD\alpha\varphidx=(-1)|\alpha|\int\Omega\varphivdx    forall\varphi\in

infty(\Omega),
C
c

then we call

v

the weak

\alpha

-th partial derivative
of

u

. If there exists a weak

\alpha

-th partial derivative of

u

, then it is uniquely defined almost everywhere, and thus it is uniquely determined as an element of a Lebesgue space. On the other hand, if

u\inCk(\Omega)

, then the classical and the weak derivative coincide. Thus, if

v

is a weak

\alpha

-th partial derivative of

u

, we may denote it by

D\alphau:=v

.

For example, the function

u(x)=\begin{cases} 1+x&-1<x<0\\ 10&x=0\\ 1-x&0<x<1\\ 0&else \end{cases}

is not continuous at zero, and not differentiable at −1, 0, or 1. Yet the function

v(x)=\begin{cases} 1&-1<x<0\\ -1&0<x<1\\ 0&else \end{cases}

satisfies the definition for being the weak derivative of

u(x),

which then qualifies as being in the Sobolev space

W1,p

(for any allowed

p

, see definition below).

The Sobolev spaces

Wk,p(\Omega)

combine the concepts of weak differentiability and Lebesgue norms.

Sobolev spaces with integer k

One-dimensional case

In the one-dimensional case the Sobolev space

Wk,p(\R)

for

1\lep\leinfty

is defined as the subset of functions

f

in

Lp(\R)

such that

f

and its weak derivatives up to order

k

have a finite norm. As mentioned above, some care must be taken to define derivatives in the proper sense. In the one-dimensional problem it is enough to assume that the

(k{-}1)

-th derivative

f(k-1)

is differentiable almost everywhere and is equal almost everywhere to the Lebesgue integral of its derivative (this excludes irrelevant examples such as Cantor's function).

With this definition, the Sobolev spaces admit a natural norm,

\|f\|k,p=\left

k
(\sum
i=0

\left\|f(i)\right

p
\|
p
1
p
\right)

=\left

k
(\sum
i=0

\int\left|f(i)(t)\right|pdt\right

1
p
)

.

One can extend this to the case

p=infty

, with the norm then defined using the essential supremum by

\|f\|k,infty=maxi=0,\ldots,k\left\|f(i)\right\|infty=maxi=0,\ldots,k\left(ess\supt\left|f(i)(t)\right|\right).

Equipped with the norm

\|\|k,p,Wk,p

becomes a Banach space. It turns out that it is enough to take only the first and last in the sequence, i.e., the norm defined by

\left\|f(k)\right\|p+\|f\|p

is equivalent to the norm above (i.e. the induced topologies of the norms are the same).

The case

Sobolev spaces with are especially important because of their connection with Fourier series and because they form a Hilbert space. A special notation has arisen to cover this case, since the space is a Hilbert space:

Hk=Wk,2.

The space

Hk

can be defined naturally in terms of Fourier series whose coefficients decay sufficiently rapidly, namely,

Hk(T)= \{f\inL2(T):

infty
\sum
n=-infty

\left(1+n2+n4+...+n2k\right)\left|\widehat{f}(n)\right|2<infty \},

where

\widehat{f}

is the Fourier series of

f,

and

T

denotes the 1-torus. As above, one can use the equivalent norm
2
\|f\|
k,2
infty
=\sum
n=-infty

\left(1+|n|2\right)k\left|\widehat{f}(n)\right|2.

Both representations follow easily from Parseval's theorem and the fact that differentiation is equivalent to multiplying the Fourier coefficient by

in

.

Furthermore, the space

Hk

admits an inner product, like the space

H0=L2.

In fact, the

Hk

inner product is defined in terms of the

L2

inner product:

\langle

u,v\rangle
Hk

=

k
\sum
i=0

\left\langleDiu,Div\right

\rangle
L2

.

The space

Hk

becomes a Hilbert space with this inner product.

Other examples

In one dimension, some other Sobolev spaces permit a simpler description. For example,

W1,1(0,1)

is the space of absolutely continuous functions on (or rather, equivalence classes of functions that are equal almost everywhere to such), while

W1,infty(I)

is the space of bounded Lipschitz functions on, for every interval . However, these properties are lost or not as simple for functions of more than one variable.

All spaces

Wk,infty

are (normed) algebras, i.e. the product of two elements is once again a function of this Sobolev space, which is not the case for

p<infty.

(E.g., functions behaving like |x|−1/3 at the origin are in

L2,

but the product of two such functions is not in

L2

).

Multidimensional case

The transition to multiple dimensions brings more difficulties, starting from the very definition. The requirement that

f(k-1)

be the integral of

f(k)

does not generalize, and the simplest solution is to consider derivatives in the sense of distribution theory.

A formal definition now follows. Let

k\in\N,1\leqslantp\leqslantinfty.

The Sobolev space

Wk,p(\Omega)

is defined to be the set of all functions

f

on

\Omega

such that for every multi-index

\alpha

with

|\alpha|\leqslantk,

the mixed partial derivative

f(\alpha)=

\partial|f
\partial
\alpha1
x
1
...\partial
\alphan
x
n

exists in the weak sense and is in

Lp(\Omega),

i.e.

\left\|f(\alpha)\right

\|
Lp

<infty.

That is, the Sobolev space

Wk,p(\Omega)

is defined as

Wk,p(\Omega)=\left\{u\inLp(\Omega):D\alphau\inLp(\Omega)\forall|\alpha|\leqslantk\right\}.

k

is called the order of the Sobolev space

Wk,p(\Omega).

There are several choices for a norm for

Wk,p(\Omega).

The following two are common and are equivalent in the sense of equivalence of norms:

\|u

\|
Wk,(\Omega)

:=\begin{cases}\left(\sum|\alpha\left\|D\alphau\right

p
\|
Lp(\Omega)
1
p
\right)

&1\leqslantp<infty;\ max|\left\|D\alphau\right

\|
Linfty(\Omega)

&p=infty;\end{cases}

and

\|u

\|'
Wk,(\Omega)

:=\begin{cases}\sum|\left\|D\alphau\right

\|
Lp(\Omega)

&1\leqslantp<infty;\ \sum|\left\|D\alphau\right

\|
Linfty(\Omega)

&p=infty.\end{cases}

With respect to either of these norms,

Wk,p(\Omega)

is a Banach space. For

p<infty,Wk,p(\Omega)

is also a separable space. It is conventional to denote

Wk,2(\Omega)

by

Hk(\Omega)

for it is a Hilbert space with the norm

\|

\|
Wk,(\Omega)
.

Approximation by smooth functions

It is rather hard to work with Sobolev spaces relying only on their definition. It is therefore interesting to know that by the Meyers–Serrin theorem a function

u\inWk,p(\Omega)

can be approximated by smooth functions. This fact often allows us to translate properties of smooth functions to Sobolev functions. If

p

is finite and

\Omega

is open, then there exists for any

u\inWk,p(\Omega)

an approximating sequence of functions

um\inCinfty(\Omega)

such that:

\left\|um-u\right

\|
Wk,p(\Omega)

\to0.

If

\Omega

has Lipschitz boundary, we may even assume that the

um

are the restriction of smooth functions with compact support on all of

\Rn.

Examples

In higher dimensions, it is no longer true that, for example,

W1,1

contains only continuous functions. For example,

|x|-1\inW1,1(B3)

where

B3

is the unit ball in three dimensions. For

k>n/p

, the space

Wk,p(\Omega)

will contain only continuous functions, but for which

k

this is already true depends both on

p

and on the dimension. For example, as can be easily checked using spherical polar coordinates for the function

f:Bn\to\R\cup\{infty\}

defined on the n-dimensional ball we have:

f(x)=|x|-\alpha\inWk,p(Bn)\Longleftrightarrow\alpha<\tfrac{n}{p}-k.

Intuitively, the blow-up of f at 0 "counts for less" when n is large since the unit ball has "more outside and less inside" in higher dimensions.

Absolutely continuous on lines (ACL) characterization of Sobolev functions

Let

1\leqslantp\leqslantinfty.

If a function is in

W1,p(\Omega),

then, possibly after modifying the function on a set of measure zero, the restriction to almost every line parallel to the coordinate directions in

\Rn

is absolutely continuous; what's more, the classical derivative along the lines that are parallel to the coordinate directions are in

Lp(\Omega).

Conversely, if the restriction of

f

to almost every line parallel to the coordinate directions is absolutely continuous, then the pointwise gradient

\nablaf

exists almost everywhere, and

f

is in

W1,p(\Omega)

provided

f,|\nablaf|\inLp(\Omega).

In particular, in this case the weak partial derivatives of

f

and pointwise partial derivatives of

f

agree almost everywhere. The ACL characterization of the Sobolev spaces was established by Otto M. Nikodym (1933); see .

A stronger result holds when

p>n.

A function in

W1,p(\Omega)

is, after modifying on a set of measure zero, Hölder continuous of exponent

\gamma=1-\tfrac{n}{p},

by Morrey's inequality. In particular, if

p=infty

and

\Omega

has Lipschitz boundary, then the function is Lipschitz continuous.

Functions vanishing at the boundary

See also: Trace operator.

The Sobolev space

W1,2(\Omega)

is also denoted by

H1(\Omega).

It is a Hilbert space, with an important subspace
1
H
0(\Omega)
defined to be the closure of the infinitely differentiable functions compactly supported in

\Omega

in

H1(\Omega).

The Sobolev norm defined above reduces here to
\|f\|
H1

=\left(\int\Omega|f|2+|\nablaf|2

12
\right)

.

When

\Omega

has a regular boundary,
1
H
0(\Omega)
can be described as the space of functions in

H1(\Omega)

that vanish at the boundary, in the sense of traces (see below). When

n=1,

if

\Omega=(a,b)

is a bounded interval, then
1
H
0(a,b)
consists of continuous functions on

[a,b]

of the form

f(x)=

x
\int
a

f'(t)dt,    x\in[a,b]

where the generalized derivative

f'

is in

L2(a,b)

and has 0 integral, so that

f(b)=f(a)=0.

When

\Omega

is bounded, the Poincaré inequality states that there is a constant

C=C(\Omega)

such that:

\int\Omega|f|2\leqslantC2\int\Omega|\nablaf|2,    f\in

1
H
0(\Omega).

When

\Omega

is bounded, the injection from
1
H
0(\Omega)
to

L2(\Omega),

is compact. This fact plays a role in the study of the Dirichlet problem, and in the fact that there exists an orthonormal basis of

L2(\Omega)

consisting of eigenvectors of the Laplace operator (with Dirichlet boundary condition).

Traces

See also: Trace operator.

Sobolev spaces are often considered when investigating partial differential equations. It is essential to consider boundary values of Sobolev functions. If

u\inC(\Omega)

, those boundary values are described by the restriction

u|\partial\Omega.

However, it is not clear how to describe values at the boundary for

u\inWk,p(\Omega),

as the n-dimensional measure of the boundary is zero. The following theorem resolves the problem:

Tu is called the trace of u. Roughly speaking, this theorem extends the restriction operator to the Sobolev space

W1,p(\Omega)

for well-behaved Ω. Note that the trace operator T is in general not surjective, but for 1 < p < ∞ it maps continuously onto the Sobolev–Slobodeckij space
1-1,p
p
W

(\partial\Omega).

Intuitively, taking the trace costs 1/p of a derivative. The functions u in W1,p(Ω) with zero trace, i.e. Tu = 0, can be characterized by the equality

1,p
W
0

(\Omega)=\left\{u\inW1,p(\Omega):Tu=0\right\},

where

1,p
W
0

(\Omega):=\left\{u\inW1,p(\Omega):\exists\{um\}

infty\subset
m=1
infty(\Omega),
C
c

suchthatum\tourm{in}W1,p(\Omega)\right\}.

In other words, for Ω bounded with Lipschitz boundary, trace-zero functions in

W1,p(\Omega)

can be approximated by smooth functions with compact support.

Sobolev spaces with non-integer k

Bessel potential spaces

For a natural number k and one can show (by using Fourier multipliers) that the space

Wk,p(\Rn)

can equivalently be defined as

Wk,p(\Rn)=Hk,p(\Rn):= \{f\inLp(\Rn):l{F}-1[(1+|\xi|2)

k
2

l{F}f]\inLp(\Rn) \},

with the norm

\|f\|
Hk,p(\Rn)

:=\left\|l{F}-1[(1+|\xi|2)

k
2

l{F}f]

\right\|
Lp(\Rn)

.

This motivates Sobolev spaces with non-integer order since in the above definition we can replace k by any real number s. The resulting spaces

Hs,p(\Rn):=\left\{f\inlS'(\Rn):l{F}-1\left[(1+|\xi|2

s
2
)

l{F}f\right]\inLp(\Rn)\right\}

are called Bessel potential spaces[1] (named after Friedrich Bessel). They are Banach spaces in general and Hilbert spaces in the special case p = 2.

For

s\geq0,Hs,p(\Omega)

is the set of restrictions of functions from

Hs,p(\Rn)

to Ω equipped with the norm
\|f\|
Hs,p(\Omega)

:=inf\left

\{\|g\|
Hs,p(\Rn)

:g\inHs,p(\Rn),g|\Omega=f\right\}.

Again, Hs,p(Ω) is a Banach space and in the case p = 2 a Hilbert space.

Using extension theorems for Sobolev spaces, it can be shown that also Wk,p(Ω) = Hk,p(Ω) holds in the sense of equivalent norms, if Ω is domain with uniform Ck-boundary, k a natural number and . By the embeddings

Hk+1,p(\Rn)\hookrightarrowHs',p(\Rn)\hookrightarrowHs,p(\Rn)\hookrightarrowHk,p(\Rn),k\leqslants\leqslants'\leqslantk+1

the Bessel potential spaces

Hs,p(\Rn)

form a continuous scale between the Sobolev spaces

Wk,p(\Rn).

From an abstract point of view, the Bessel potential spaces occur as complex interpolation spaces of Sobolev spaces, i.e. in the sense of equivalent norms it holds that

\left[Wk,p(\Rn),Wk+1,p(\Rn)\right]\theta=Hs,p(\Rn),

where:

1\leqslantp\leqslantinfty, 0<\theta<1,s=(1-\theta)k+\theta(k+1)=k+\theta.

Sobolev–Slobodeckij spaces

Another approach to define fractional order Sobolev spaces arises from the idea to generalize the Hölder condition to the Lp-setting. For

1\leqslantp<infty,\theta\in(0,1)

and

f\inLp(\Omega),

the Slobodeckij seminorm (roughly analogous to the Hölder seminorm) is defined by

[f]\theta,:=\left(\int\Omega\int\Omega

|f(x)-f(y)|p
|x-y|\theta

dxdy\right

1
p
)

.

Let be not an integer and set

\theta=s-\lfloors\rfloor\in(0,1)

. Using the same idea as for the Hölder spaces, the Sobolev–Slobodeckij space[2]

Ws,p(\Omega)

is defined as

Ws,p(\Omega):=\left\{f\inW\lfloor(\Omega):\sup|\alpha|[D\alphaf]\theta,<infty\right\}.

It is a Banach space for the norm

\|f\|

Ws,(\Omega)

:=

\|f\|
W\lfloor(\Omega)

+\sup|\alpha|[D\alphaf]\theta,.

If

\Omega

is suitably regular in the sense that there exist certain extension operators, then also the Sobolev–Slobodeckij spaces form a scale of Banach spaces, i.e. one has the continuous injections or embeddings

Wk+1,p(\Omega)\hookrightarrowWs',p(\Omega)\hookrightarrowWs,p(\Omega)\hookrightarrowWk,(\Omega),k\leqslants\leqslants'\leqslantk+1.

There are examples of irregular Ω such that

W1,p(\Omega)

is not even a vector subspace of

Ws,p(\Omega)

for 0 < s < 1 (see Example 9.1 of [3])

From an abstract point of view, the spaces

Ws,p(\Omega)

coincide with the real interpolation spaces of Sobolev spaces, i.e. in the sense of equivalent norms the following holds:

Ws,p(\Omega)=\left(Wk,p(\Omega),Wk+1,p(\Omega)\right)\theta,,k\in\N,s\in(k,k+1),\theta=s-\lfloors\rfloor.

Sobolev–Slobodeckij spaces play an important role in the study of traces of Sobolev functions. They are special cases of Besov spaces.

The constant arising in the characterization of the fractional Sobolev space

Ws,p(\Omega)

can be characterized through the Bourgain-Brezis-Mironescu formula:

\lims(1-s) \int\Omega\int\Omega

|f(x)-f(y)|p
|x-y|s

dxdy=

2
n-1
2
\pi
\Gamma
(p+1
2
)
p\Gamma
(p+n
2
)

\int\Omega\vert\nablaf\vertp;

and the condition

\limsups(1-s) \int\Omega\int\Omega

|f(x)-f(y)|p
|x-y|s

dxdy<infty

characterizes those functions of

Lp(\Omega)

that are in the first-order Sobolev space

W1,p(\Omega)

[4] .

Extension operators

If

\Omega

is a domain whose boundary is not too poorly behaved (e.g., if its boundary is a manifold, or satisfies the more permissive "cone condition") then there is an operator A mapping functions of

\Omega

to functions of

\Rn

such that:
  1. Au(x) = u(x) for almost every x in

\Omega

and

A:Wk,p(\Omega)\toWk,p(\Rn)

is continuous for any 1 ≤ p ≤ ∞ and integer k.

We will call such an operator A an extension operator for

\Omega.

Case of p = 2

Extension operators are the most natural way to define

Hs(\Omega)

for non-integer s (we cannot work directly on

\Omega

since taking Fourier transform is a global operation). We define

Hs(\Omega)

by saying that

u\inHs(\Omega)

if and only if

Au\inHs(\Rn).

Equivalently, complex interpolation yields the same

Hs(\Omega)

spaces so long as

\Omega

has an extension operator. If

\Omega

does not have an extension operator, complex interpolation is the only way to obtain the

Hs(\Omega)

spaces.

As a result, the interpolation inequality still holds.

Extension by zero

Like above, we define

s
H
0(\Omega)
to be the closure in

Hs(\Omega)

of the space
infty
C
c(\Omega)
of infinitely differentiable compactly supported functions. Given the definition of a trace, above, we may state the following

If

u\in

s
H
0(\Omega)
we may define its extension by zero

\tildeu\inL2(\Rn)

in the natural way, namely

\tildeu(x)=\begin{cases}u(x)&x\in\Omega\ 0&else\end{cases}

For its extension by zero,

Ef:=\begin{cases}f&rm{on}\Omega,\ 0&rm{otherwise}\end{cases}

is an element of

Lp(\Rn).

Furthermore,

\|Ef

\|
Lp(\Rn)

=\|f

\|
Lp(\Omega)

.

In the case of the Sobolev space W1,p(Ω) for, extending a function u by zero will not necessarily yield an element of

W1,p(\Rn).

But if Ω is bounded with Lipschitz boundary (e.g. ∂Ω is C1), then for any bounded open set O such that Ω⊂⊂O (i.e. Ω is compactly contained in O), there exists a bounded linear operator

E:W1,p(\Omega)\toW1,p(\Rn),

such that for each

u\inW1,p(\Omega):Eu=u

a.e. on Ω, Eu has compact support within O, and there exists a constant C depending only on p, Ω, O and the dimension n, such that

\|Eu

\|
W1,p(\Rn)

\leqslantC

\|u\|
W1,p(\Omega)

.

We call

Eu

an extension of

u

to

\Rn.

Sobolev embeddings

See main article: Sobolev inequality.

It is a natural question to ask if a Sobolev function is continuous or even continuously differentiable. Roughly speaking, sufficiently many weak derivatives (i.e. large k) result in a classical derivative. This idea is generalized and made precise in the Sobolev embedding theorem.

Write

Wk,p

for the Sobolev space of some compact Riemannian manifold of dimension n. Here k can be any real number, and 1 ≤ p ≤ ∞. (For p = ∞ the Sobolev space

Wk,infty

is defined to be the Hölder space Cn where k = n + α and 0 < α ≤ 1.) The Sobolev embedding theorem states that if

k\geqslantm

and

k-\tfrac{n}{p}\geqslantm-\tfrac{n}{q}

then

Wk,p\subseteqWm,q

and the embedding is continuous. Moreover, if

k>m

and

k-\tfrac{n}{p}>m-\tfrac{n}{q}

then the embedding is completely continuous (this is sometimes called Kondrachov's theorem or the Rellich–Kondrachov theorem). Functions in

Wm,infty

have all derivatives of order less than m continuous, so in particular this gives conditions on Sobolev spaces for various derivatives to be continuous. Informally these embeddings say that to convert an Lp estimate to a boundedness estimate costs 1/p derivatives per dimension.

There are similar variations of the embedding theorem for non-compact manifolds such as

\Rn

. Sobolev embeddings on

\Rn

that are not compact often have a related, but weaker, property of cocompactness.

See also

References

translation of Mat. Sb., 4 (1938) pp. 471–497.

External links

Notes and References

  1. Bessel potential spaces with variable integrability have been independently introduced by Almeida & Samko (A. Almeida and S. Samko, "Characterization of Riesz and Bessel potentials on variable Lebesgue spaces", J. Function Spaces Appl. 4 (2006), no. 2, 113–144) and Gurka, Harjulehto & Nekvinda (P. Gurka, P. Harjulehto and A. Nekvinda: "Bessel potential spaces with variable exponent", Math. Inequal. Appl. 10 (2007), no. 3, 661–676).
  2. In the literature, fractional Sobolev-type spaces are also called Aronszajn spaces, Gagliardo spaces or Slobodeckij spaces, after the names of the mathematicians who introduced them in the 1950s: N. Aronszajn ("Boundary values of functions with finite Dirichlet integral", Techn. Report of Univ. of Kansas 14 (1955), 77–94), E. Gagliardo ("Proprietà di alcune classi di funzioni in più variabili", Ricerche Mat. 7 (1958), 102–137), and L. N. Slobodeckij ("Generalized Sobolev spaces and their applications to boundary value problems of partial differential equations", Leningrad. Gos. Ped. Inst. Učep. Zap. 197 (1958), 54–112).
  3. 2012-07-01 . Hitchhikerʼs guide to the fractional Sobolev spaces. Bulletin des Sciences Mathématiques . en. 136. 5. 521–573 . 10.1016/j.bulsci.2011.12.004. 0007-4497 . free. Di Nezza. Eleonora. Palatucci. Giampiero. Valdinoci. Enrico. 1104.4345.
  4. Book: Bourgain, Jean. Jean Bourgain . Brezis . Haïm . Haïm Brezis. Mironescu . Petru . 2001 . Menaldi . José Luis . Optimal control and partial differential equations. In honour of Professor Alain Bensoussan’s 60th birthday. Proceedings of the conference, Paris, France, December 4, 2000. Amsterdam: IOS Press; Tokyo: Ohmsha . 439-455 . Another look at Sobolev spaces . 978-1-58603-096-4.