Ψec: A Local Spectral Exterior Calculus
Christian Lessig1 Otto-von-Guericke-Universit¨atMagdeburg
Abstract
We introduce Ψec(Rn), a discretization of Cartan’s exterior calculus of differ- ential forms using wavelets. Our construction consists of differential r-form wavelets with flexible directional localization that provide tight frames for the spaces Ωr(Rn) of forms in R2 and R3. By construction, the wavelets satisfy the de Rahm co-chain complex, the Hodge decomposition, and that the k- dimensional integral of an r-form is an (r − k)-form. They also verify Stokes’ theorem for differential forms, with the most efficient finite dimensional approx- imation attained using directionally localized, curvelet- or ridgelet-like forms. The construction of Ψec(Rn) builds on the geometric simplicity of the exterior calculus in the Fourier domain. We establish this structure by extending existing results on the Fourier transform of differential forms to a frequency description of the exterior calculus, including, for example, a Plancherel theorem for forms and a description of the symbols of all important operators. Keywords: exterior calculus, wavelets, structure preserving discretization 2010 MSC: 42C15, 42C40, 65T60, 58A10
1. Introduction
The exterior calculus, first introduced by Cartan [1], provides a formulation of scalar and vector-valued functions that encodes their differential relation- ships in a coordinate-invariant language. The central objects of the calculus are differential forms, which provide a covariant formalization of vector fields, and the exterior derivative, which is a first order differential operator generalizing
arXiv:1811.12269v3 [cs.NA] 5 Oct 2020 gradient, curl, and divergence. The importance of the exterior calculus for the numerical solution of partial differential equations, even in Rn and with a flat metric, has been shown in various applications. One of the first ones was electromagnetic theory where already in the 1960s it was realized [2] that the electric and magnetic fields have to be discretized differently to obtain satisfactory numerics. Later, it was found that this could be understood by the electric and magnetic fields being differential forms of different degrees [3, 4]. Motivated by these and analogous
Preprint submitted to Applied and Computational Harmonic Analysis October 7, 2020 results in other fields, e.g. in geophysical fluid dynamics [5, 6, 7], various discrete realizations of Cartan’s exterior calculus have been proposed, e.g. [8, 9, 10, 11, 12, 13]. These have been applied to applications such as elasticity [10], fluid mechanics [14, 15], climate simulation [7, 11, 16], magnetohydrodynamics [17], electromagnetic theory [3], and geometry processing [18]. In the present work, we propose Ψec(Rn), an alternative discretization of Cartan’s exterior calculus that is defined over a family of differential form wavelets with flexible directional localization. Towards its construction, we extend existing results from the theoretical physics literature and develop a de- scription of the exterior calculus in the Fourier domain. This shows that under the Fourier transform it becomes a chain complex whose structure is encoded in a simple, geometric way in spherical coordinates in frequency space. The exterior derivative, for instance, acts only on the radial component of a dif- ferential form in frequency space. We therefore discretize the exterior calculus in the Fourier domain using spherical coordinates and attain space-frequency localization by using polar wavelet windows [19, 20] that respect the spherical coordinate structure. The construction yields differential form wavelets well lo- calized in space and frequency that obey by construction the de Rahm complex and the distinction between exact and co-exact forms, i.e. the Hodge-Helmholtz decomposition. At the same time, the definition in spherical coordinates also allows for a flexible directional localization, between fully isotropic and curvelet- and ridgelet-like. In contrast to most existing discretizations in the literature, the “discrete” differential form wavelets in our work are bona fide forms in the sense of the continuous theory. Hence, all operations from the continuous set- ting are well defined. A central question therefore becomes efficient numerical evaluation and we will at least partially address it in the present work. n r,ν,n The properties of Ψec(R ) and its differential form wavelets ψs,a , where n = 2, 3 is the dimension of the space, 0 ≤ r ≤ n is the degree of the form, ν ∈ {d, δ} denotes if the wavelet is exact or co-exact, and s = (js, ks, ts) describes level js, translation ks and orientation ts, can be summarized as follows (see Sec. 5 for the precise technical statements):
r,ν,n ˙ |ν| r n • The ψs,a form tight frames for the homogeneous Sobolev spaces L2 (Ων , R ). • The differential form wavelets satisfy the exterior co-chain complex, i.e. r,δ,n r+1,d,n r,d,n dψs,a = ψs,a and dψs,a = 0. Furthermore, the exterior derivative d preserves the localization described by s.
• Stokes’ theorem for differential forms α ∈ Ωr(Rn), Z Z α = dα, ∂M M
for M ⊂ Rn a manifold, holds for differential form wavelets as Z Z X r−1,δ X r,d αs ψs ∧ χ∂M = αs ψs ∧ χM (1) n n s R s R
2 r n where the αs are the frame coefficients of α ∈ Ω (R ) in the differential
form wavelet frame, and χM and χ∂M are the characteristic functions of M and its boundary ∂M, respectively. When dim(M) = n, Stokes theorem therefore contains the approximation problem for cartoon-like functions (right hand side of Eq. 1). This has been studied extensively in the literature on curvelets, shearlets and related constructions, e.g. [21, 22, 23, 24], and these provide for it quasi-optimal approximation rates.
Moreover, χ∂M is exactly the wavefront set associated with χM .
k,ν,n • The Hodge dual ? ψs,a has a simple and practical description. r,ν,n • The fiber integral of ψs,a along an arbitrary p-dimensional linear sub- r−p,ν,n 0 manifold is again a differential form wavelet ψs0,a and with s preserv- ing the localization of s.
• The Laplace-de Rahm operator ∆ : Ωk → Ωk has a closed form represen- r,ν,n tation and Galerkin projection for the differential form wavelets ψs,a . r,ν,n • The ψs,a have closed form expressions in the spatial domain. This clar- ifies, e.g. how the angular localization in frequency space is carried over to the spatial domain. Multiplicative operators in the exterior calculus, such as the wedge product or the Lie derivative, have currently no natural expression in our calculus. In numerical calculations, these can be determined using the transform method [25, 26], i.e. with the multiplication being evaluated pointwise in space, which is facilitated by the closed form expressions that are available there. A thorough investigation of this has, however, be left for future work. The remainder of the work is structured as follows. After discussing re- lated work in the next section, we present in Sec. 3 some background on the polar wavelets that we use as localization windows in the frequency domain. Afterwards we develop in Sec. 4 the Fourier transform of differential forms. In Sec. 5 we introduce our differential form wavelets and study how the exterior calculus, including the exterior derivative, the wedge product, Hodge dual, and Stokes theorem, can be expressed using these. In Sec. 6 we summarize our work and discuss directions for future work. Our notation and conventions are summarized in Appendix A.
2. Related Work
Our work builds on ideas from what are usually two distinct fields: firstly, wavelet theory and in particular wavelets constructed in polar and spherical coordinates, and, secondly, exterior calculus and its discretizations. We will discuss relevant literature in these fields in order.
3 Polar Wavelets. The idea to combine space-frequency localization with a de- scription in polar coordinates in the Fourier domain goes back at least to work by Fefferman in the theory of partial differential equations, in particular on the second dyadic decomposition, cf. [27, Ch. IX]. It was re-introduced multiple times, e.g. in work on steerable wavelets [28, 29] and curvelets [30, 21, 22] where the motivation was the analysis of directional features in images. A framework that encompasses the latter two approaches was proposed by Unser and co-workers [31, 32, 19, 33] and this defines the polar wavelets that are the scalar building block for the construction of Ψec(Rn). We use the term ‘polar wavelet’, or short ‘polarlet’, since it encapsulates what distinguishes the con- struction from other multi-dimensional wavelets and what is the key to their utility in our work.
Curl and Divergence Free Wavelets. Various constructions for curl and diver- gence free wavelets have been proposed over the years, e.g. [34, 35, 36, 37, 38, 39], and these are related to our differential form wavelets by the musical iso- morphism. Based on seminal work by Lemair´e-Rieusset[34], various authors, e.g. [40, 41, 42, 43, 44] constructed coupled scalar wavelet bases for 1-dimensional domains that satisfy a differential relationship so that the derivative of one yields (up to a constant) the other or a corresponding integral relationship holds. These constructions can be extended to curl or divergence free wavelets using the coordinate representations of the differential operators of vector calculus. Our differential form wavelets satisfy a similar coupling under differentiation. However, our construction is inherently multi-dimensional and we exploit the simple structure of the differential operators in the frequency domain. The latter has also been used in the work by Bostan, Unser and Ward [39] although diver- gence freedom is there only enforced numerically. A precursor to the presented work is those in [45]. There, polar wavelets are used as scalar window functions for the construction of polar divergence free wavelets and it is exploited that a vector field is divergence free when its Fourier transform is tangential to the frequency sphere.
Discretizations of Exterior Calculus. Discretizations of Cartan’s exterior cal- culus try to preserve important structures of the continuous theory, e.g. the de Rahm co-chain complex or Stokes’ theorem, to improve the accuracy and robustness of numerical calculations. Different precursors those good numerical performance can, at least in hindsight, be explained through their adherence to the exterior calculus, exist. For example, spectral discretizations in geophysical fluid dynamics, first introduced in the 1950s and 1960s, yield energy and en- strophy conserving time integrators [5] by implicitly using what one could call a spectral exterior calculus [46]. In the 1970s, Arakawa and Lamb [6, 7] introduced the now eponymous grid-based discretizations for geophysical fluid dynamics to extend energy and enstrophy conservation to finite difference methods, again without an explicit connection to Cartan’s calculus but implicitly respecting its structure. A few years later, N´ed´elec[47, 48], introduced consistent discretiza- tions for the spaces H(curl) and H(div) that respect important properties of the
4 classical differential operators of vector calculus, and hence also of the exterior derivative. The first work that, to our knowledge, made an explicit connection between exterior calculus and numerics was those by Bossavit [3] in compu- tational electromagnetism who realized that the use of differential forms and the calculus defined on them enables one to systematically understand and ex- tend earlier results that provided good numerical performance for problems in computational electromagnetism [2]. Most existing discretizations of exterior calculus are finite element-based, such as finite element exterior calculus developed by Arnold and collabora- tors [10, 49], discrete exterior calculus by Hirani, Desbrun and co-workers [50, 9], mimetic discretizations by Bochev and Hyman [8], and the TRiSK scheme by Thuburn, Ringler, Klemp, and Skamarock [51, 11]. These works all have in common that they construct a discrete analogue of Cartan’s exterior calculus that satisfies important structural properties of the continuous theory, e.g. in that it is a co-chain complex for a discrete exterior derivative or that a dis- crete Stokes’ theorem holds. Our construction, in contrast, is best seen as a multi-resolution projection of the continuous theory that preserve its intrinsic structure. In particular, our differential form wavelets remain genuine forms in the sense of the continuous theory. All operators from there are thus naturally defined. Compared to finite element type discretizations where one seeks dis- crete definitions mimicking the original theory, the questions shifts in our work thus to how these can be computed efficiently numerically or how this can be ensured in the construction of the wavelets. Wavelets and exterior calculus have so far only been combined in two works and again within a finite element framework. Dubos, Kevlahan and collabora- tors [16, 52, 53] extended the TRiSK scheme to a multi-resolution formulation using second generation wavelets and demonstrated that this provides an effi- cient discretization for geophysical fluid dynamics. In work contemporaneous to ours, Budninskiy, Owhadi, and Desbrun [54] constructed operator-adapted wavelets for Discrete Exterior Calculus using the efficient algorithms proposed in [55]. In contrast to Ψec(Rn) whose construction relies on the structure of exterior calculus in the Fourier domain, the above works are defined purely in the spatial domain. The extension of finite element-based discretizations of exterior calculus to spectral convergence rates has been considered by Rufat, Mason and Des- brun [12]. These authors accomplished it by introducing interpolation and histapolation maps to relate discrete exterior calculus and continuous forms. Gross and Atzberger presented a similar approach using hyper-interpolation for radial manifolds, i.e. manifolds that can be described as a height field over S2. A discretizations of exterior calculus for subdivision surfaces was introduced by de Goes et al. [13]. It shares the multi-resolution structure with our approach and, at least in the limit case, also yields continuous differential forms. Another recent discretization of exterior calculus is those by Berry and Gi- annakis [56] for manifold learning problems. They employ the eigenfunctions of the Laplace-de Rahm operator as bases for their forms. This is analogous to the spectral bases, formed by spherical harmonics, first used in geophysical
5 fluid dynamics and that, as mentioned earlier, essentially form a spectral exte- rior calculus. We also build on a spectral formulation but localize it in space and frequency and, furthermore, exploit the additional structure afforded in the Euclidean setting.
3. Polar Wavelets
Polar wavelets are a family of wavelets defined in polar or spherical coor- dinates in the Fourier domain. Their construction uses a compactly supported radial window, hˆ(|ξ|), which controls the overall frequency localization, and an angular one,γ ˆ(ξ¯), which controls the directionality, where ξ¯ = ξ/|ξ|. The mother wavelet is thus given by ψˆ(ξ) = ψˆ(ξ¯) hˆ(|ξ|) with the whole family of functions being generated by dilation, translation and rotation. In two dimensions, the angular window can be described with a Fourier series in the polar angle θˆ. A polar wavelet is there hence given by
j Nj 2 X ˆ j ψˆ (ξ) ≡ ψˆ (ξ) = βj,t einθ hˆ(2−j|ξ|) e−ihξ,2 ki (2) s jkt 2π n n=−Nj where the multi-index s = (j, k, t) describes level j, translation k, and orienta- tion t and we will write (js, ks, ts) whenever confusion might arise. The coeffi- j,t j,t cients βn control the angular localization and in the simplest case βn = δn0 for all j, t and one has isotropic, bump-like wavelet functions. Conversely, when the j,t support of the βn are the entire integers in n, then one can describe compactly supported angular windows. The above formulation then also encompasses sec- ond generation curvelets [22] and provides more generally a framework to realize α-molecule-like constructions [57, 58, 59]. A useful property of polar wavelets is that their spatial representation, given by inverse Fourier transform F−1 of Eq. 2, can be computed in closed form. Using the Fourier transform in polar coordinates, cf. Appendix A.4, one ob- tains [20]
Nj 2j X ψ (x) ≡ ψ (x) = in βj,t einθx h (|2jx − k|) (3) s jkt 2π n n n=−Nj ˆ ˆ where hn(|x|) is the Hankel transform of h(|ξ|) of order n. When h(|ξ|) is the window proposed for the steerable pyramid [60], then hn(|x|) also has a closed form expression [20]. Note also that the angular localization described by the j,t βn is essentially invariant under the Fourier transform and only modified by the factor of in that implements a rotation by π/2. For applications, the wavelets in Eq. 3 are augmented using scaling functions φj,k. They are isotropic and defined using a radial windowg ˆ(|ξ|), i.e. φj,0 = −1 F (ˆg), that is chosen such that the φj,k represent a signal’s low frequency part. To simplify notation, we follow the usual convention and write ψ−1,k ≡ φ0,k. With the scaling functions, the following result holds [19, 45].
6 Tj ×2Nj +1 Proposition 1. Let Uj ∈ R be the matrix formed by the angular local- j,t j −int(2π/Tj ) ization coefficients βn = βn e for the Tj different orientations, and 2Nj +1×2Nj +1 let Dj ∈ R be a diagonal matrix. When the Calder´onadmissibility 2 P∞ ˆ −j 2 2 H condition gˆ(|ξ|) + j=0 h(2 |ξ|) = 1, ∀ξ ∈ Rb , is satisfied and Uj Uj = Dj 2 with tr(Dj) = 1 for all levels j, then any f ∈ L2(R ) has the representation
∞ Tj X X X f(x) = f(y), ψj,k,t(y) ψj,k,t(x) 2 j=−1 k∈Z t=1 with frame functions
2j ψ (x) = ψ R (2jx − k) j,k,t 2π 2πt/Tj
where ψ(x) is given by Eq. 3 and R2πt/Tj is the rotation by 2πt/Tj. Although the above frame is redundant, since it is Parseval tight it still affords most of the conveniences of an orthonormal basis. In the following, we will often write S for the index set of all scales, translations, and rotations. In three dimensions, polar wavelets are defined analogous to Eq. 2 by
3j/2 Lj l 2 X X j ψˆ (ξ) ≡ ψˆ (ξ) = κj,t y (ξ¯) hˆ(2−j|ξ|) e−ihξ,2 ki (4) s j,k,t (2π)3/2 lm lm l=0 m=−l where the directional localization windowγ ˆ(ξ¯) is now expressed using the spher- j,t ical harmonics ylm and the κlm are the expansion coefficients; see Appendix A.3 for a brief summary of harmonic analysis of the sphere. The wavelets have again a closed form expression in the spatial domain,
Lj l 23j/2 X X ψ (x) ≡ ψ (x) = ilκj,t y (¯x) h (|2jx − k|). (5) s j,k,t (2π)3/2 lm lm l l=0 m=−l which can be obtained using the Rayleigh formula, cf. Appendix A.4. Further- more, the wavelets can again be augmented using scaling functions to represent a signal’s low frequency part. Then the following result holds [33, 45] (see again Appendix A.3 for notation).
2 (Lj +1) Proposition 2. Let wj,t ∈ R be the vector formed by the rotated angular j,t Pl m0 localization coefficients κlm = m0=−l W (λj,t)l,m κl,m0 for a localization win- m dow centered at λj,t, where W (λj,t)l,m is the Wigner-D matrix implementing 2 2 rotation in the spherical harmonics domain, and let Glm ∈ R(Lj +1) ×(Lj +1) be the matrix formed by the spherical harmonics product coefficients for fixed 2 P∞ ˆ −j 2 3 (l, m). When the Calder´oncondition gˆ(|ξ|) + j=0 h(2 |ξ|) = 1, ∀ξ ∈ Rb PTj lm 3 is satisfied and δl,0δm,0 = t=0 wj,t G wj,t for all j, t then any f(x) ∈ L2(R )
7 has the representation
∞ Tj X X X f(x) = f(y) , ψj,k,t(y) ψj,k,t(x) 3 j=−1 k∈Z t=1 with frame functions
23j/2 ψ (x) = ψ R (2jx − k), j,k,t (2π)3/2 λj,t
for ψ(x) defined in Eq. 5 and Rλj,t the rotation from the North pole to λj,t. See Appendix B for a proof of the admissibility conditions. We refer to the original works [19, 33] and [20] for a more detailed discussion of polar wavelets.
4. Differential Forms and their Fourier Transform
In this section, we first recall basic facts about differential forms and the exterior calculus defined on them. Our principle reference for this material will be the book by Marsden, Ratiu, and Abraham [61] and we will also use the notation and conventions from there. Afterwards, starting from existing results in theoretical physics [62, 63, 64], we will define and study the Fourier transform of the exterior calculus.
4.1. Tangent and Cotangent Bundle The setting of our discussion in the following will be Rn with n = {2, 3}. Although the tangent bundle T Rn can be identified with Rn × Rn, for our work it is of importance to consistently distinguish Rn and the tangent and cotangent bundles, T Rn and T ∗Rn, respectively. We will denote a generic basis for the tangent space as
∂ ∂ span , ··· , = T n ∂x1 ∂xn xR
i n with the ∂/∂x in general being dependent on the fiber TxR they span (such as in spherical coordinates). The ∂/∂xi induce a dual basis dxi for the cotangent ∗ n i j bundle T R through the fiberwise duality condition dx ∂/∂x = δij. We will also typically use the Einstein summation convention over repeated indices and hence employ “upstairs” indices for components of contra-variant vectors n ∗ n in TxR and “downstairs” ones for covariant vectors in Tx R . For example, a vector field X ∈ X(Rn) is in coordinates given by Xi ∂/∂xi and a co-vector field i by αi dx (summation implied).
8 4.2. Differential Forms and the Exterior Calculus A differential r-form α ∈ Ωr(Rn) on Rn, with 0 ≤ r ≤ n, is a covariant, anti- symmetric tensor field of rank r, i.e. an anti-symmetric, multi-linear functional on the rth power of the tangent bundle. When r = 0 this corresponds to functions, i.e. Ω0(Rn) =∼ F(Rn), and for r = n one has densities, e.g. f(x) dx. All other degrees r correspond to “vector-valued” fields. For example, a 1-form is a co-vector field, i.e. a section of the cotangent bundle T ∗Rn, and since it is a functional on T Rn it naturally pairs with the tangent vector field of a curve. A 2-form, correspondingly, pairs with the tangent vectors of a 2-manifold and in R3 it can hence be thought of as an infinitesimal area. The result of the pairing between an r-form and the tangent vectors of an r-manifold can be integrated and one thus also says that an r-form is to be integrated over an r-manifold [9, 65]. A principal motivation for studying differential forms is that they are closed under the exterior derivative d,
r n r+1 n d : Ω (R ) → Ω (R ) (6) which provides a coordinate invariant derivative for tensor fields. The exterior derivative furthermore satisfies d ◦ d = 0. The differential forms under d hence form a co-chain complex known as the de Rahm complex,
d 0 n d d n n d 0 −−→ Ω (R ) −−→· · · −−→ Ω (R ) −−→ 0.
When a differential form α ∈ Ωr(Rn) satisfies dα = 0 then α is closed and when there exists a β such that β = dα then α is exact. The exterior calculus also provides a product on differential forms known as the wedge product
∧ :Ωr × Ωl → Ωr+l. (7)
It turns the de Rahm complex into a graded algebra. i n ∗ n With the basis {dx }i=1 for the cotangent bundle T R , the coordinate expression for a 1-form α ∈ Ω1(R3), i.e. a section of T ∗Rn, is 1 n α(x) = α1(x) dx + ··· + αn(x) dx . The basis functions for higher degree forms are generated using the wedge prod- uct in Eq. 7, which ensures anti-symmetry. In R3, for example, one has 0 3 Ω (R ) 3 f = f(x) 1 3 1 2 3 Ω (R ) 3 α = α1(x) dx + α2(x) dx + α3(x) dx 2 3 2 3 3 1 1 2 Ω (R ) 3 β = β1(x) dx ∧ dx + β2(x) dx ∧ dx + β3(x) dx ∧ dx 3 3 1 2 3 Ω (R ) 3 γ = γ(x) dx ∧ dx ∧ dx .
For an arbitrary form α ∈ Ωr(Rn) the coordinate expression is
X j1 jr α(x) = αj1,··· ,jr (x) dx ∧ · · · ∧ dx
j1<··· 9 with the summation being over all ordered sequences j1 < ··· < jr, which again ensures the anti-symmetry of the form. Except for our discussions of Stokes’ theorem at the end of this section and in Sec. 5.2.1, we will assume that the coordinate functions αj1,··· ,jr (x) are smooth, or at least sufficiently smooth for the operation under consideration, and real- valued. Differential forms with coordinate functions that are distributions in the sense of Schwartz were introduced by de Rahm [66] and are known as ‘currents’. We will denote the space of such (n − r)-forms as S0(Ωn−r, Rn) and they are dual to the space S(Ωr, Rn) of differential r-forms whose coordinate functions are test functions. The duality pairing hh·, ·ii between S0(Ωr, Rn) and S(Ωr, Rn) is, as in the scalar case, analogous to the L2-inner product for differential forms, which will be introduced in the next paragraph. Hilbert Space Structure and Hodge Dual. So far we did not require a metric and this is an important feature of the exterior calculus. When a metric is available, as is the case in Rn, the exterior calculus can be equipped with additional structure. We will introduce this structure next. 0 n Let a non-degenerate, symmetric, covariant tensor field g ∈ T2 (R ) of rank −1 2 n 2 be a given metric and g ∈ T0 (R ) be the associated co-metric acting on ij co-vectors; the components will be denoted as gij and g , respectively. Unless mentioned otherwise, we will work with the canonical metric on Rn (which, however, is not the identity in spherical coordinates). The metric g and its dual g−1 can be used to identify vectors and co-vectors by “raising” and “lowering” indices. Let X,Y ∈ X(Rn) and α, β ∈ X∗(Rn) =∼ Ω1(Rn). Then [ i i [ i 1 n X (Y ) ≡ g(X,Y ) = gijX Y ⇔ (X )j = gij X ∈ Ω (R ) ] ij ] j ij n α (β) ≡ g(α, β) = g αi βj ⇔ (α ) = g αi ∈ X(R ). The flat operator [ and the sharp operator ] introduced above are together known as musical isomorphism. With these, the classical differential operators from vector calculus can be related to the exterior derivative. Let f ∈ F(R3) and u ∈ X(R3) be a vector field. Then ∇f = (df)] (9a) ∇ × u = (?(du[))]. (9b) ∇ · u = ?d(?u[). (9c) The metric induces, next to the musical isomorphism, also an L2-inner prod- uct hh·, ·ii for differential r-forms. Let α, β ∈ Ωr(Rn). Then Z j1,··· ,jr 1 n hhα, βii = αj1,··· ,jr β dx ∧ · · · ∧ dx (10a) n R where the contra-variant tensor components βj1,··· ,jr are obtained by index rais- ing with the metric, j1,··· ,jr j1,i1 jr ,ir β = g ··· g βi1,··· ,ir . (10b) 10 r n The space L2(Ω , R ) is defined in the usual way as all differential r-forms with finite L2 norm. Unless stated otherwise (e.g. for Stokes’ theorem), we will r n assume that all forms we work with are in L2(Ω , R ). Eq. 10 can equivalently be written using the Hodge dual ?β as Z hhα, βii = α ∧ ?β n R and the equality with Eq. 10 defines ?β. Since the integral in the last equation has to be over a volume form, i.e. one of maximum degree r = n, the Hodge dual has to satisfy r n n−r n ? :Ω (R ) → Ω (R ) and it also holds that ? ? α = (−1)n(n−r)α. By linearity of the Hodge dual, which follows from those of the contraction in Eq. 10b and the wedge product, for the coordinate expression of ?β it suffices to know those for the form basis functions. For R3, for example, these are ?1 = dx1 ∧ dx2 ∧ dx3 (11a) ?dx1 = dx2 ∧ dx3 (11b) ?dx2 = −dx1 ∧ dx3 (11c) ?dx3 = dx1 ∧ dx2. (11d) 4.3. The Fourier Transform of the Exterior Algebra Motivation. Before introducing the formal definition of the Fourier transform of differential forms, it is instructive to understand the duality that underlies it. The classical result that integration becomes under the Fourier transform pointwise evaluation, i.e. Z f(x) dx = fˆ(0), (12) n R 1 n n n suggests that the volume form f(x) dx = f1···ndx ∧ · · · ∧ dx ∈ Ω (R ) should map to a 0-form in frequency space, since for 0-forms “integration” amounts to pointwise evaluation, cf. [65]. This is also required by covariance because if the volume form f(x) dx would map to another volume form (or to a form of degree other than zero) under the Fourier transform then it would “pick up” a Jacobian term under coordinate transformations. But fˆ(0) is already the integrated value and hence invariant under changes of coordinates (e.g. the mass of an object or the energy of a field). This observation for volume forms can be consistently extended to differential forms of lower degree. For example, for a plane in R3, which from the point of view of the exterior algebra is a 2-form, the Fourier transform is known to be in the direction of the plane’s normal, i.e. a 1-form in this direction in frequency space. For the x1-x2 plane we have, for example, F δ1,2(x) = δ3(ξ) 11 3 where δ1,2(x) is the distribution representing the x1-x2 plane in R and δ3(ξ) its Fourier transform in the ξ3 direction. With the above considerations, the duality between differential forms in space and frequency is expected to be where we denote by Ωbn−r(Rbn) the exterior algebra in frequency space Rbn. The spatial coordinate x = xi ∂/∂xi ∈ Rn pairs with the covector ξ = i i n ix ξi ξi dx ∈ Rb in the complex exponential e of the Fourier transform. When we require that the usual correspondence between differential operators ∂/∂xi i and their Fourier transform ξidx (no summation) is to hold (in an appropriate sense, see Proposition 4 below) we also have to expect that the exterior algebra in frequency space is naturally defined over ∂/∂ξi ≡ ∂/∂xi where we use the notation ∂/∂ξi to be explicit if we are considering the primary or dual domain with respect to the Fourier transform. Differential forms in frequency space are hence defined over form basis functions ∂/∂ξj1 ∧ · · · ∧ ∂/∂ξjr and this provides the second cornerstone for the definition of the Fourier transform of differential forms. The Fourier Transform of Differential Forms. With the above motivation, we define the Fourier transform of differential forms, following an approach intro- duced in mathematical physics [63, 67, 64] based on the theory of super vector spaces. Definition 1. Let α ∈ Ωr(Rn) be a differential form of degree r, X j1 jr α(x) = αj1,··· ,jr (x) dx ∧ · · · ∧ dx . j1<··· n−r n Then the Fourier transform αb = F(α) ∈ Ωb (Rb ) of α is the (n − r)-form on frequency space Rbn given by Z 1 X j j −i xpξ i dxq ∂ F(α)(ξ) = α (x) dx 1 ∧ · · · ∧ dx r ∧ e p ∧ e ∂ξq n/2 j1,··· ,jr (2π) n j1<··· 12 From a functional analytic perspective, the usual requirements for the exis- tence of the Fourier transform apply with respect to the coordinate functions α (x) and αj1,··· ,jr (ξ). j1,··· ,jr b dxq ∂ Remark 1. In the literature, e.g. [63, 64], the exponential e ∂ξq of the form basis functions is typically defined with a complex unit in the exponent, which is considered there as only serving “aesthetic” purposes to visually closer match the scalar case. For real-valued functions this is, however, not true and the use of idxq ∂ e ∂ξq leads to inconsistencies in the inverse Fourier transform. In particular, there are no differential r-form wavelets that satisfy the de Rahm complex and that are real-valued in space. We hence use the definition above without the complex unit. To evaluate the exponential in Eq. 13 one uses that the exterior algebra has a Z2-graded structure, with the even grade given by the differential forms of even degree and the odd one by differential forms of odd degree.2 The Z2-grading is respected in the tensor product by using the following multiplication rule (for homogenous elements):