arxiv: v1 [math.ap] 9 Mar 2017

Transcription

1 FRACTIONAL SOBOLEV METRICS ON SPACES OF IMMERSED CURVES MARTIN BAUER, MARTINS BRUVERIS, AND BORIS KOLEV arxiv: v1 [math.ap] 9 Mar 2017 Abstrat. Motivated by appliations in the field of shape analysis, we study reparametrization invariant, frational order Sobolev-type metris on the spae of smooth regular urves Imm(, R d ) and on its Sobolev ompletions I q (, R d ). We prove loal well-posedness of the geodesi equations both on the Banah manifold I q (, R d ) and on the Fréhetmanifold Imm(, R d ) provided the order of the metri is greater or equal to one. In addition we show that the H s -metri indues a strong Riemannian metri on the Banah manifold I s (, R d ) of the same order s, provided s > 3. These investigations an be also interpreted as 2 a generalization of the analysis for right invariant metris on the diffeomorphism group. Contents 1. Introdution 1 Connetions to shape analysis 3 2. Parametrized and unparametrized urves 5 3. Riemannian metris on immersions 5 4. The geodesi equation 8 5. Smoothness of the metri Smoothness of the spray Strong Riemannian metris 18 Appendix A. Smooth urves in S r (Z) 20 Appendix B. Square-root of an operator in Eλ r (Z) 20 Referenes Introdution The interest in Riemannian metris on infinite-dimensional manifolds is fueled by their onnetions to mathematial physis and in partiular fluid dynamis. It was Arnold who disovered in 1966 that the inompressible Euler equation, whih desribes the motion of an ideal fluid, has an interpretation as the geodesi equation on an infinite-dimensional manifold; the manifold in question is the group of volume-preserving diffeomorphisms equipped with the L 2 -metri. Sine then many other PDEs in mathematial physis have been reinterpreted as geodesi equations. Examples inlude Date: Marh 10, Mathematis Subjet Classifiation. 58D05, 35Q35. Key words and phrases. Sobolev metris of frational order. 1

2 2 MARTIN BAUER, MARTINS BRUVERIS, AND BORIS KOLEV Burgers equation, whih is the geodesi equation of the L 2 -metri on the group of all diffeomorphisms of the irle, Diff( ), and the Camassa Holm equation [11], the geodesi equation of the H 1 -metri [23] on the same group. Interestingly, geodesi equations orresponding to frational orders in the Sobolev sale have also found appliations in physis: Wunsh showed that the geodesi equation of the homogenous H 1/2 -metri on Diff( ) is onneted to the Constantin Lax Majda equation [14, 40, 18], whih itself is a simplified model of the vortiity equation. The geometri interpretation of a PDE as the geodesi equation enables one to show loal well-posedness of the PDE. This was done first by Ebin and Marsden [16] for the Euler equation. Using a similar method Constantin and Kolev showed in [13] that the geodesi equation of Sobolev H n -metris on Diff( ) with integer n 1 is loally well-posed. In [17] this was extend by Esher and Kolev to the Sobolev H r -metris of frational order r 1 2. Similar results were shown for the diffeomorphism group of ompat manifolds by Shkoller in [33, 34] and by Preston and Misiolek in [31]. Frational metris on Diff(R d ) have been studied in [5] by Bauer, Esher and Kolev. The loal well-posedness of the geodesi equation for frational order metris on the diffeomorphism group of a general manifold M remains an open problem. In this paper we study the loal well-posedness of a family of PDEs that arise as geodesi equations on the spae Imm(, R d ) of immersed urves. To be preise Imm(, R d ) onsists of smooth, losed urves with nowhere vanishing derivatives. We an regard the diffeomorphism group Diff( ) Imm(, ) as an open subset of the spae of immersions. If we replae on the right hand side by R d we obtain the spae of urves. The PDE tt = D s t, D s t t, D s t D s 1 2 t 2 D 2 s, where D s = 1 θ and = (t, θ) is a time-dependent urve, is the geodesi equation for the L 2 -metri G (h, k) = h, k dθ. This is a (weak) Riemannian metri on Imm(, R d ). The weight in the integral makes the metri invariant under the natural Diff( )-ation and leads to the appearane of ar length derivatives D s in the geodesi equation. This PDE an be seen as a generalization of the geodesi equation on Diff( ), φ tt = 2 φ tφ tx, φ x whih beomes, when written in terms of the Eulerian veloity u = φ t φ, Burgers equation, u t = 2u x u. While the behaviour of Burgers equation is well-known, to our knowledge, nothing is known about the loal well-posedness of the L 2 -geodesi equation on the spae of urves.

3 FRACTIONAL SOBOLEV METRICS ON SPACES OF IMMERSED CURVES 3 The situation improves when we add higher derivatives to the Riemannian metri. By doing so, we get the H n -metri G (h, k) = A h, k dθ, A = () j α j Ds 2j, where α j 0 are onstants. The orresponding geodesi equation is with (A t ) t = D s t, D s A t A t, D s t D s W (, t )D 2 s, W (, t ) = 1 2 t j () k+1 Ds 2j k t, Ds k t. j=1 k=1 Interestingly, the behaviour of the geodesi equation is better understood in the seemingly more ompliated ase where n 1. For n = 1 the PDE is loally but generally not globally well-posed [30], while for n 2 the PDE has solutions that are global in time [9]. Inspired by related work [5, 17] on geodesi equations on the diffeomorphism group we onsider H r -metris of non-integer order r, i.e. G (h, k) = A h, k dθ, with A being, for eah fixed urve, a Fourier multiplier of order 2r; the preise assumptions made on A are desribed in Setion 3. The geodesi equation takes the form (A t ) t = D s t, D s A t A t, D s t D s (w(, t ) + w 0 (, t )D 2 s, where w(, t ) and w 0 (, t ) are expressions that are defined in Theorem 4.1. It is a nonlinear evolution equation of seond order in t and of order 2s in θ; the right hand side is quadrati in t and highly nonlinear in. For non-integer s, both A t and w(, t ) are nonloal funtions of both and t. While not immediately obvious, we show in Example 4.3 that it ontains the geodesi equations for the L 2 -metri and the H n -metris as speial ases. Our ontribution is to show in Corollary 6.5 that the geodesi equation for the H r -metri is loally well-posed in the Sobolev spae H q for q 2r and r 1. Connetions to shape analysis. Throughout the artile we will assume that the operator A defining the metri is equivariant with respet to the ation of the diffeomorphism group Diff( ), i.e., A φ (h φ) = (A h) φ, φ Diff( ). This is equivalent to requiring that the Riemannian metri G is invariant under the ation of Diff( ), meaning G φ (h φ, k φ) = G (h, k), φ Diff( ). This assumption is neessary in order to apply the lass of Sobolev metris in shape analysis. The mathematial analysis of shapes has been the fous of intense researh in reent years [21, 26, 15, 42, 35] and has found appliations in fields suh as image analysis, omputer vision, biomedial imaging and funtional

4 4 MARTIN BAUER, MARTINS BRUVERIS, AND BORIS KOLEV Figure 1. Example of a geodesi between two shapes with respet to an H 2 -metri [2]. data analysis. An important lass of shapes are outlines of planar objets. Mathematially, these shapes an be represented by equivalene lasses of parametrized urves modulo the reparametrization group Diff( ). This yields the following geometri piture, π : Imm(, R d ) Imm(, R d )/Diff( ). A key step in shape analysis is to define an effiiently omputable distane funtion between shapes in order to measure similarity between shapes. This distane funtion an be the geodesi distane indued by a Riemannian metri. Sine it is diffiult to work with the quotient Imm(, R d )/Diff( ) diretly, the standard approah is instead to define a Riemannian metri on Imm(, R d ), that is Diff( )-invariant. Suh a metri an then indue a Riemannian metri on the quotient spae, making π a Riemannian submersion; this metri would be given by the formula G π() (u, u) = inf G (h, h). T π.h=u Hene we will only look at Diff( )-invariant metris on Imm(, R d ) in this paper. The arguably simplest invariant metri is the L 2 -metri, orresponding to the operator A = Id. However, its indued geodesi distane funtion is identially zero, both, on the spae of parametrized urves as well as the quotient spae of unparametrized urves; this was shown by Mihor and Mumford [28, 29], see also [1]. This means that between any two urves there exists a path of arbitrary short length. This property makes the metri illsuited for most appliations in shape analysis, beause a notion of distane between shapes is one of the basi tools there. As a onsequene higher order metris, mostly of integer order, were studied and used suessfully in appliations; see [30, 22, 36, 41, 39, 37]. For an overview on various metris on the shape spae of parametrized and unparametrized urves see [30, 3, 4]. As an example we have inluded an optimal deformation between two unparametrized urves with respet to a seond order Sobolev metri in Figure 1. Frational order metris have been briefly mentioned [27, Setion 3]. However, so far, an analysis of the well-posedness of the orresponding geodesi equation was missing. This is the question we will fous on in this paper.

5 FRACTIONAL SOBOLEV METRICS ON SPACES OF IMMERSED CURVES 5 2. Parametrized and unparametrized urves In this artile we onsider the spae of smooth regular urves with values in R d Imm(, R d ) := { C (, R d ) : } 0. As an open subset of the Fréhet spae C (, R d ), it is a Fréhet manifold. Its tangent spae at any urve is the vetor spae of smooth funtions: T Imm(, R d ) = C (, R d ). The group of smooth diffeomorphisms of the irle Diff( ) := { φ C (, ) : φ > 0 } ats on the spae of regular urves via omposition from the right: Imm(, R d ) Diff( ) Imm(, R d ), (, φ) φ. Taking the quotient with respet to this group ation, we obtain the shape spae of un-parameterized urves B i (, R d ) := Imm(, R d )/Diff( ). Note that the ation of Diff( ) on Imm(, R d ) is not free, and thus that the quotient spae B i (, R d ) is not a manifold, but only an orbifold with finite isotropy groups. A way to overome this diffiulty is to onsider the slightly smaller spae of free immersions Imm f (, R d ), i.e., those immersions upon whih Diff( ) ats freely. This spae is an open and dense subset of Imm(, R d ) and the orresponding quotient spae is again a Fréhet manifold, see [12]. B i,f (, R d ) := Imm f (, R d )/Diff( ) 3. Riemannian metris on immersions Let G be a Riemannian metri on Imm(, R d ). Motivated by appliations in the field of shape analysis we require G to be invariant with respet to the diffeomorphism group Diff( ), i.e., (3.1) G (h, k) = G φ (h φ, k φ), φ Diff( ). This invariane is a neessary assumption for G to indue a Riemannian metri on the shape spae B i,f (, R d ), suh that the projetion map is a Riemannian submersion. We assume that the metri is given in the form (3.2) G (h, k) = A h, k ds, where A : C (, R d ) C (, R d ) is a ontinuous linear operator that depends on the foot point. Assoiated to the metri is a map ˇG : T Imm(, R d ) T * Imm(, R d ). In terms of A we have ˇG = A ds. If G is Diff( )-invariant, then the family of operators A has to be Diff( )-equivariant, A φ (h φ) = A (h) φ, φ Diff( ).

6 6 MARTIN BAUER, MARTINS BRUVERIS, AND BORIS KOLEV Example 3.1 (Integer order Sobolev metris). Motivated by their appliability in shape analysis [43, 38, 27, 30] important examples are Sobolev metris with onstant oeffiients, (3.3) G n (h, k) = α j D j s h, Dsk S j ds, 1 and sale-invariant Sobolev metris, (3.4) Gn (h, k) = α j l 2j 3 D j s h, D j sk ds, with onstants α j 0, for j = 0,..., n; one requires α 0, α n > 0 and alls n the order of the metri. Here, D s = θ denotes differentiation with respet to ar length, ds = dθ integration with respet to ar length and l = S ds the orresponding urve length. 1 Using integration by parts we obtain a formula for the operator A : (3.5) A n = () j α j Ds 2j for metris with onstant oeffiients and (3.6) An = () j α j l 2j 3 D 2j s. for sale-invariant metris. Beause it will be important later, we note that if is a onstant speed urve, then the operator for a metri with onstant oeffiients is A n = () j α j ( 2π l ) 2j 2j θ ; it is a differential operator with onstant oeffiients and the oeffiients depend only on the length of the urve. Until now loal and global well-posedness results have been established only for the Sobolev metris of integer order [30, 7, 9, 6]. The main goal of this artile is to extend these results to metris of frational order; in partiular to metris, for whih the operator A is defined using Fourier multipliers of a ertain lass. Given a urve Imm(, R d ), let ψ Diff( ) be a diffeomorphism suh that ψ has onstant speed. Reparametrization invariane of the metri implies that (3.7) G (h, k) = G ψ (h ψ, k ψ ) ; in other words, G is determined by its behaviour on onstant speed urves. Remark 3.2. The situation is similar to that of right-invariant metris on diffeomorphism groups, whih are determined by their behaviour at the identity diffeomorphism. For urves, the invariane property is weaker and the spae of onstant speed urves is still quite large.

7 FRACTIONAL SOBOLEV METRICS ON SPACES OF IMMERSED CURVES 7 Let us write the invariane property (3.7) in terms of A : by a straightforward alulation we obtain G (h, k) = G ψ = S 1 = whih implies the identity (h ψ, k ψ ) A ψ (h ψ R ψ A ψ A = R ψ A ψ ), k ψ R ψ (h), k R ψ. l 2π dθ ds, The lass of metris on onstant speed urves is large. To make it more manageable we will restrit the possible dependane of A on the urve. Assumption. We assume from this point onwards, that the operator A ψ depends on only through its length l ; in other words A ψ = A(l ), where λ A(λ) is a smooth urve with values in the spae of linear maps, L(C (, R d ), C (, R d )) 1. Then, with this assumption, (3.8) A = R ψ A(l ) R ψ. Requiring this form for the operator A imposes a restrition on A(l ). For α onsider φ α Diff( ), defined as φ α (θ) = θ + α mod 2π. If a urve has onstant speed, then so does φ α and thus l = l φα. Therefore (3.8) implies A(l ) R φα = R φα A(l ), or equivalently after differentiating with respet to α, A(l ) θ = θ A(l ). This means that A(l ) has to be a Fourier multiplier, i.e., A(l ).u(θ) = m Z a(l, m).ˆu(m) exp(imθ), where a(l, ) : Z L(C d ) is alled the symbol of A(l ). We will write A(l ) = a(l, D) or A(l ) = op (a(l, )). For now we make no assumptions on the symbol apart from requiring that the map l A(l ) is smooth. More ontrol on the symbol will be neessary in order to prove well-posedness of the geodesi equation. Example 3.3 (Integer order Sobolev metris). For Sobolev metris with onstant oeffiients the symbol of the operator A n (l ) is ( ) a n (l, m) = 2πm 2j () j α j I d, l 1 This is equivalent to requiring that the map R+ C (, R d ) C (, R d ) given by (λ, h) A(λ)h is smooth.

8 8 MARTIN BAUER, MARTINS BRUVERIS, AND BORIS KOLEV with I d C d d the identity matrix, and for sale-invariant metris it is a n (l, m) = () j α j (2πm) 2j I d, l 3 Remark 3.4. Even though all known situations orrespond to Fourier multipliers of the form A(l ) = a(l, D) with a(l, m) = a(l, m)i d, a multiple of the identity matrix, we will treat in this artile general matrix-valued symbols, beause doing so introdues no additional diffiulties. 4. The geodesi equation Geodesis between two urves 0, 1 Imm(, R d ) are ritial points of the energy funtional E() = 1 1 A t, t ds dt, 2 0 where = (t, θ) is a path in Imm(, R d ) joining 0 and 1. The geodesi equation is obtained by omputing the derivative of E on paths with fixed endpoints. We will use the following notations and t := t := θ v := = D s, A = R ψ A(l ) R ψ A := R ψ A (l ) R ψ, where A (l ) is the derivative of A(l ) with respet to the parameter l and ψ (θ) = 2π l θ 0 dσ. We have the following ommutation rules for D s and R ψ, D s R ψ = 2π l R ψ θ θ R ψ = l 2π R ψ D s, whih imply that beause A(l ) ommutes with θ, that the operator A ommutes with D s, D s A = A D s. This will be used when omputing the geodesi equation. Theorem 4.1. Assume that, for eah λ R +, the operator A(λ) : C (, R d ) C (, R d ) is invertible with a ontinuous inverse. Then the weak Riemannian metri (3.2) on Imm(, R d ) has a geodesi spray whih is given by (4.1) F : (, h) (h, S (h)), T Imm(, R d ) T T Imm(, R d ) where (4.2) S (h) = A {(D,h A )h + D s h, v A h + A h, D s h v + (w + w 0 )D s v}.

9 FRACTIONAL SOBOLEV METRICS ON SPACES OF IMMERSED CURVES 9 with w(, h) = A h, D s h ds and 1 w 0 (, h) = 2π A h, ψ D s h (l A + A )h, h ds. The geodesi equation is (4.3) (A t ) t = D s t, v A t A t, D s t v (w(, t ) + w 0 (, t ))D s v. To ompute these equations, we will first derive a variational formula for the reparametrization funtion ψ. Lemma 4.2. The derivative of the map ψ is given by D,h ψ (θ) = 2π θ D s h, v d s 1 ) D s h, v ds ψ (θ). l l ( 0 Proof. Using the variational formulas for and l, D,h = D s h, v D,h l = D s h, v ds, the lemma follows diretly from ψ (θ) = 2π l θ 0 dσ. Proof of Theorem 4.1. We an write the energy of a path (t, θ) as E() = 1 1 l A(l )( t ψ ), t ψ dθ dt π The variation of the energy in diretion h is then given by 1 l D,h E = A(l )( t ψ ), h t ψ + ( t ψ 0 2π + 1 l 2 2π D (,hl l A(l ) + A (l ) ) ( t ψ S 1 1 = A t, h t + ( t D,h ψ ) ψ ds D,hl (l A + A ) t, t ds dt. ) D,h ψ dθ ), t ψ dθ dt Let us start with the term involving h t. After integrating by parts we get 1 1 A t, h t ds dt = (A t ) t + (A t ) t( ) 0 0, h ds dt 1 = (A t ) t + D s t, v A t, h ds dt. 0 Next we deal with ψ. Using the formula ( 1 = D,h ψ ) ψ D,h ψ ψ,

10 10 MARTIN BAUER, MARTINS BRUVERIS, AND BORIS KOLEV and Lemma 4.2 we obtain a formula for the variation of ψ, ( D,h ψ ) ψ (θ) = 1 ( θ D s h, v d s 1 ( D s h, v ds 0 2π ) ) ψ (θ). Therefore we have A t, ( t D,h ψ ) ψ ds ( θ = A t, D s t D s h, v d s 1 ( ) ) D s h, v d s ψ ds 0 2π S 1 θ = A t, D s t D s h, v d s ds D s h, v ds A t, ψ D s t ds 2π Beause A is symmetri and ommutes with D s, A t, D s t ds = 0, and thus the funtion w(θ) = θ 0 A t, D s t ds, is periodi with D s w = A t, D s t. Next we integrate by parts, A t, ( t D,h ψ ) ψ ds = w D s h, v ds + 1 D s h, v ds A t, ψ D s t ds 2π S 1 = D s (wv), h ds 1 D s v, h ds A t, ψ D s t ds. 2π Finally we onsider the term involving D,h l. Sine D,h l = D s h, v ds = D s v, h ds, we have 1 2 D,hl (l A + A ) t, t ds = = 1 D s v, h ds (l A + A 2 ) t, t ds. Grouping all these expressions together, we an express D,h E as (4.4) D,h E() = where 1 0 w 0 = (A t ) t D s t, v A t D s (wv) w 0 D s v, h ds dt, ( 1 2π A t, ψ D s t + 1 ) 2 (l A + A ) t, t ds.

11 FRACTIONAL SOBOLEV METRICS ON SPACES OF IMMERSED CURVE1 Thus we obtain the geodesi equation (A t ) t = D s t, v A t D s (wv) w 0 D s v. The existene of the geodesi spray follows from (4.4), whih an be rewritten as D,h E() = 1 0 G ( tt + S ( t ), h) dt. Note that in this last step we used the invertibility of A on C (, R d ). Next we will see how the geodesi equation simplifies for Sobolev metris with onstant oeffiients and sale-invariant metris. First note, that using integration by parts we an write 1 2π A t, ψ D s t ds = 1 θ A t, D s t dσ ds l 0 = 1 θ A t, D s t d s ds. l We shall also make use of the following identity, valid for j 1 and h C (, R d ), ( ) 2j Ds 2j 1 h, D s h = D s () k+1 Ds 2j k h, Ds k h. 2 Setting W j = 1 2 we obtain as well as 2j k=1 k=1 () k+1 D 2j k s t, D k s t for j 1, W 0 = 1 2 t 2, θ 0 D 2j s t, D s t d s = W j (θ) W j (0), W j (θ) ds = (1 2j) Ds 2j t, t ds. Example 4.3 (Metris with onstant oeffiients). For metris with onstant oeffiients we have ( ) 2π 2j A n = () j α j Ds 2j, A n (l ) = () j α j 2j l θ, and thus their l -derivatives are (A n ) = 2l () j jα j Ds 2j, Therefore l A n + (A n ) = l A (l ) = 2l () j jα j ( 2π l () j (1 2j)α j D 2j s, ) 2j 2j θ.

12 12 MARTIN BAUER, MARTINS BRUVERIS, AND BORIS KOLEV and w 0 = 1 () j α j [ W j (θ) W j (0) ds + 1 ] (1 2j) Ds 2j t, t ds l 2 = () j α j W j (0). Sine w = θ () j α j Ds 2j t, D s t d s = it follows that 0 w(θ) + w 0 = () j α j (W j (θ) W j (0)), () j α j W j (θ). Thus the geodesi equation has the form (A t ) t = D s t, v A t A t, D s t v () j α j W j D s v. Thus we have regained the formula from [30]. Example 4.4 (Sale-invariant metris). For sale-invariant metris we have A n = () j α j l 2j 3 Ds 2j, and thus their l -derivatives are (A n ) = 3l () j α j l 2j 3 Ds 2j, An (l ) = l 3 A (l ) = 3l 4 () j α j (2π) 2j 2j θ, () j α j (2π) 2j 2j θ. Therefore l A similar alulation as before gives w 0 = w(θ) = and therefore w(θ) + w 0 = () j α j l 2j 3 A n + ( A n ) = 2l A n. ( W j (0) + 1l ( j () j α j l 2j 3 (W j (θ) W j (0)), () j α j l 2j 3 ) ) Ds 2j t, t ds ( ( W j (θ) + 1l j + 1 ) ) Ds 2j t, t ds. 2

13 FRACTIONAL SOBOLEV METRICS ON SPACES OF IMMERSED CURVE3 5. Smoothness of the metri For any q > 3 2 we an onsider the Sobolev ompletion Iq = I q (, R d ), whih is an open set of the Hilbert vetor spae H q (, R d ). It is the aim of this setion to show that the metri G extends to a smooth weak metri on I q, for high enough q. It turns out that the smoothness of the metri redues to the smoothness of the mapping A = R ψ A(l ) R ψ. We will begin our investigations with this question. It is well-known (see for instane [16, Appendix A]) that, given a differential operator A of order r with smooth oeffiients and a diffeomorphism ψ, the onjugate operator A ψ = R ψ A R ψ is again a differential operator of order r, whose oeffiients are polynomial expressions in ψ, its derivatives and 1/ψ ; in partiular, the mapping ψ A ψ = R ψ A R ψ, D q ( ) L(H q (, R d ), H q r (, R d )) is smooth for q > 3/2 and q r 1, where D q ( ) is group of H q - diffeomorphisms of. In [17] Esher and Kolev extended this result to Fourier multipliers by showing that the map ψ A ψ remains smooth when A is a Fourier multiplier of lass S r (Z) (to be defined below). The proof that the mapping A is smooth that we present here is inspired by [17, 5], but in the present ase we have to deal with Fourier multipliers that depend (in a nie way) on a parameter λ. First we introdue the lasses S r (Z) and Sλ r (Z) of Fourier multipliers (for further details see [32] for instane). A Fourier multiplier a(d) ats on a funtion u via (5.1) a(d).u(θ) = m Z a(m).ˆu(m) exp(imθ), and its symbol is a funtion a : Z L(C d ). Definition 5.1. Given r R, the Fourier multiplier a(d) belongs to the lass S r (Z), if a(m) satisfies α a(m) C α m r α, for eah α N, where m := (1 + m 2 ) 1/2. Here we define the differene operator α via a(m) = a(m + 1) a(m), α = α. We equip the spae S r (Z) with the topology indued by the seminorms p α (a(d)) = sup α a(m) m (r α), m Z where is some norm on L(C d ). With this topology S r (Z) is a Fréhet spae. The lass S r (Z) oinides with the one defined in [17]; the equivalene is shown in [32].

14 14 MARTIN BAUER, MARTINS BRUVERIS, AND BORIS KOLEV Example 5.2. Any linear differential operator of order r with onstant oeffiients belongs to S r (Z). Furthermore the operator Λ 2r = op ( m 2r), whih defines the H r -norm via u H r = Λ 2r u u dθ, belongs to S 2r (Z). Remark 5.3. A Fourier multiplier a(d) of lass S r (Z) extends for any q R to a bounded linear operator and the linear embedding is ontinuous. a(d) : H q (, R d ) H q r (, R d ), S r (Z) L(H q (, R d ), H q r (, R d )) Now we introdue the lass of parameter-dependent symbols. Definition 5.4. Given r R, the lass S r λ (Z) = C (R +, S r (Z)) onsists of smooth urves R + S r (Z), λ a(λ, D) of Fourier multipliers. We will all suh a urve a λ-symbol. We will often write a λ (D) for a(λ, D) or, by slight abuse of notation, to denote the whole urve λ a(λ, D). Using the material in Appendix A, we an give an alternative desription of the spae Sλ r (Z) of one-parameter families of symbols. Lemma 5.5. Given r R, a family of smooth urves a(, m) C (R +, L(C d )) with m Z defines an element a λ (D) in the lass Sλ r (Z) if and only if for eah α, β N, (5.2) β λ α a(λ, m) C α,β m r α, holds loally uniformly in λ R +. Proof. Let a λ (D) Sλ r (Z). Sine all derivatives of a smooth urve are loally ) bounded, it follows that for all α, β N, the expression p α ( β λ a λ(d) is bounded loally uniformly in λ and hene a(λ, m) satisfies (5.2) for some onstants C α,β. Conversely, let a family a(λ, m) of smooth urves satisfying the above estimates be given. Define the urves a β (λ, m) = β λa(λ, m). The estimate (5.2) shows that the urves λ a β (λ, ) are loally bounded in S r (Z) and hene by Lemma A.1 the urve λ a λ (D) is an element of Sλ r(z). The following is the main theorem that will be used to show the smoothness of the metri. Set with ψ (θ) = 2π l θ 0 dσ. A = R ψ A(l ) R ψ,

15 FRACTIONAL SOBOLEV METRICS ON SPACES OF IMMERSED CURVE5 Proposition 5.6. Let r 1 and A(λ) = a λ (D) belong to the lass S r λ (Z). Then the map A, I q (, R d ) L(H q (, R d ), H q r (, R d )), is smooth provided q > 3 2 and q r. Proof. It was established in [17, Thm. 3.7] that for A S r (Z) the mapping ψ A ψ := R ψ A R ψ, D q L(H q, H q r ) is smooth. Using the uniform boundedness priniple [24, 5.18] it follows that the map ψ (A A ψ ), D q L(S r (Z), L(H q, H q r )) is smooth 2. Beause the inlusion L(S r (Z), L(H q, H q r )) C (S r (Z), L(H q, H q r ), is bounded [24, 5.3], the following map is smooth ψ (A A ψ ), D q C (S r (Z), L(H q, H q r )) ; via the exponential law [24, 3.12] this is equivalent to the smoothness of the joint map (A, ψ) A ψ, S r (Z) D q L(H q, H q r )). Next we note that the maps l, I q R + and ψ, I q D q are smooth this an be seen from their definitions and sine S r λ (Z) = C (R +, S r (Z)) the omposition A(l ), I q S r (Z) is smooth as well. To onlude the proof we note that A an be written as the omposition (A(l ), ψ ) A(l ) ψ. Corollary 5.7. Let r 1 2 and A(λ) = a λ(d) belong to the lass Sλ 2r(Z). Then the bilinear form G (h, k) = A h, k ds extends smoothly to I q (, R d ) provided q > 3 2 Proof. It is enough to show that and q 2r. ˇG = M A, I q L(H q, H q r ), is smooth, where M denotes pointwise multipliation by. Now M, I q L(H ρ, H ρ ), is smooth for 0 ρ q 1 and the onlusion follows by Proposition 5.6 and the fat that omposition of ontinuous linear mappings between Banah spaes is smooth. 2 When E, F are Fréhet spaes or more generally onvenient vetor spaes, L(E, F ) is the spae of bounded linear maps equipped with the topology of uniform onvergene on bounded sets; see [24, 5.3] for details.

16 16 MARTIN BAUER, MARTINS BRUVERIS, AND BORIS KOLEV 6. Smoothness of the spray In order to prove the existene and smoothness of the spray, we will require moreover an elliptiity ondition on λ-symbols. For the purpose of this artile, we will adopt the following definition. Definition 6.1. An element a λ (D) of the lass Sλ r (Z) is alled loally uniformly ellipti, if a(λ, m) GL(C d ) for all (λ, m) R + Z and a(λ, m) C m r holds loally uniformly in λ. Remark 6.2. A Fourier multiplier a(d) of lass S r (Z) is ellipti, if a(m) GL(C d ) for all m Z and a(m) C m r, holds for all m. Suh an a(d) indues a bounded isomorphism between H q (, R d ) and H q r (, R d ), for all q R. We summarize our onsiderations by introduing the following lass of operators whih will be denoted E r λ (Z). Definition 6.3. A family a λ (D) of Fourier multipliers is an element of the lass Eλ r (Z), if (1) a λ (D) is in the lass Sλ r(z), (2) a(λ, m) is a positive Hermitian matrix for all (λ, m) R + Z and (3) a λ (D) is loally uniformly ellipti. Theorem 6.4. Let r 1, q 2r and A(λ) = a λ (D) belong to Eλ 2r (Z). Then G (h, k) = A h, k ds defines a smooth weak Riemannian metri of order r on I q (, R d ) with a smooth geodesi spray. Proof. It is shown in Corollary 5.7 that G extends smoothly to I q (, R d ). We an write G (h, h) = 2π m Z a(l, m)ĥ(m), ĥ(m), and beause a(l, m) are positive Hermitian matries, G (h, h) = 0 only for h = 0. Thus G is a Riemannian metri. It is a weak Riemannian metri, beause the inner produt G (, ) indues on eah tangent spae T I q (, R d ) the H r -topology, while the manifold I q (, R d ) itself arries the H q -topology and by our assumptions q > r. It remains to show that the geodesi spray of G exists and is smooth. For eah λ R +, the operator A(λ) is an ellipti Fourier multiplier and thus indues a bi-bounded linear isomorphism on C (, R d ). Thus we an apply Theorem 4.1, whih shows that the geodesi spray exists on the spae Imm(, R d ) of smooth immersions. We will show that the map S (h) extends smoothly to the Sobolev ompletion I q (, R d ). Then F (, h) = (h, S (h)) is neessarily the geodesi spray of G on I q (, R d ).

17 FRACTIONAL SOBOLEV METRICS ON SPACES OF IMMERSED CURVE7 We have the following formula for S (h), S (h) = A {(D,h A )h + D s h, v A h + A h, D s h v + (w + w 0 )D s v}. The map A A, U L(H q, H q 2r ) L(H q 2r, H q ), defined on the open subset U of invertible operators is smooth and so is A by Proposition 5.6; therefore so is the omposition A : I q L(H q 2r, H q ). Thus we need to show that (, h) (D,h A )h + D s h, v A h + A h, D s h v + (w + w 0 )D s v, is smooth between I q H q H q 2r. The first term (D,h A )h is the derivative of (, h) A h with respet to the first variable and hene the map (, h) (D,h A )h is smooth between I q H q H q 2r. Ar-length derivation (, u) D s u, I q H ρ H ρ, 0 ρ q, is smooth and pointwise multipliation in C (, R) extends to a ontinuous bilinear mapping H σ H ρ H ρ, σ > 1, 0 ρ σ. 2 Therefore, noting that v = D s and q 1 > 1 2, the expression D sh, v A h + A h, D s h v is a smooth map I q H q H q 2r. Next we use the fat that the antiderivative u ( θ θ 0 ) u(σ) dσ, H ρ H ρ+1, ρ 0, is a bounded linear mapping. Thus (, h) w(, h) maps smoothly I q H q H q 2r+1 and beause q 2r +1 > 1 2 we have that (, h) w(, h)d sv is smooth between I q H q H q 2r. The last term is w 0 (, h). First we note that the map ψ between I q D q is smooth. A term-byterm inspetion shows that the map (, h) f(, h), where f(, h) is the integrand in the definition of w 0 (, h) is a smooth map I q H q L 1 and thus (, h) w 0 (, h) is smooth as well. Thus the extension of G to I q (, R d ) has a smooth spray. As a orollary, using the Cauhy Lipshitz theorem, we obtain the loal existene of geodesis on the Hilbert manifold I q (, R d ). Corollary 6.5. Let A λ = a λ (D) belong to the lass Eλ 2r (Z), where r 1, and let q 2r. Consider the geodesi flow on the tangent bundle T I q indued by the Fourier multiplier A λ. Then, given any ( 0, h 0 ) T I q, there exists a unique non-extendable geodesi (, h) C (J, T I q ) with (0) = 0 and h(0) = h 0 on the maximal interval of existene J, whih is open and ontains 0. As a onsequene of a no-loss-no-gain argument, we also obtain loal wellposedness on the spae of smooth urves.

18 18 MARTIN BAUER, MARTINS BRUVERIS, AND BORIS KOLEV Corollary 6.6. Let A λ = a λ (D) belong to the lass Eλ 2r (Z), where r 1. Consider the geodesi flow on the tangent bundle T Imm(, R d ) indued by the Fourier multiplier A λ. Then, given any ( 0, h 0 ) T Imm(, R d ), there exists a unique non-extendable geodesi (, h) C (J, T Imm(, R d )) with (0) = 0 and h(0) = h 0 on the maximal interval of existene J, whih is open and ontains 0. Proof. Sine the metri is invariant by re-parametrization, it is invariant in partiular by translations of the parameter. A similar observation has been used in [17, Setion 4] to show that loal existene of geodesi still holds in the smooth ategory. The proof is based on a no loss-no gain result in spatial regularity, initially formulated in [16] (see also [8]). The same argument is still true here and leads to the loal existene of the geodesis on the Fréhet manifold Imm(, R d ). 7. Strong Riemannian metris The goal of this setion is to show, that metris of order s > 3 2 (indued by a Fourier multiplier A λ in the lass Eλ 2s (Z)), indue strong smooth Riemannian metris on the Sobolev ompletion I s of the same order as the metri. Let A(λ) be of lass Eλ 2s(Z) with s > 3 2. In Setion 6 it was shown that the orresponding Riemannian metri (7.1) G (h, k) = A h, k ds, an be extended smoothly to I q (, R d ) for q 2s. Now we want to improve this to q s. Of partiular interest is the ase q = s, when the topologies indued by the inner produts G (, ) oinide with the manifold topology. Any positive Hermitian matrix a has a unique positive square root b whih depends smoothly on a. We have, moreover, the stronger result that an operator A(λ) = a λ (D) in the lass Eλ 2s has a square root B(λ) = b λ(d) in the lass Eλ s (see Lemma B.2). With B = R ψ B(l ) R ψ we have the identities (7.2) A(l ) = B(l ) 2, A = B 2 ; furthermore, beause b(λ, m) is a Hermitian matrix, the operator B(l ) is L 2 (dθ)-symmetri and for eah urve the operator B is L 2 (ds)-symmetri. Therefore we an rewrite the metri G in the symmetri form (7.3) G (h, k) = A h, k ds = B h, B k ds. We obtain therefore the following expression for the operator ˇG on T I s : where ˇG = B t M B, B : H s (, R d ) L 2 (, R d ), B t : L 2 (, R d ) H s (, R d ),

19 FRACTIONAL SOBOLEV METRICS ON SPACES OF IMMERSED CURVE9 and B t is the transpose of B. The latter formula an now be used to obtain the following result onerning the smooth extension of this family of inner produts into a strong Riemannian metri on I s, provided s > 3/2. Theorem 7.1. Let s > 3/2 and B λ Eλ s (Z). Then the expression G (h, k) = B h, B k ds, defines a smooth and strong Riemannian metri on I s (, R d ). Remark 7.2. For strong metris on a Lie group the invariane of the metri implies the geodesi and metri ompleteness of the spae, see [20, Lemma 5.2]. This has been used in [5] to show ompleteness of the H s -metri on Diff(R d ), see also [10] for integer orders on the diffeomorphism group of a general manifold. Unfortunately, there is no automati analogue of this result in our situation. To prove the global well-posedness on the spae of regular urves additional assumptions on the dependene of the operator A on the length will be neessary. In future work, we plan to follow this line of researh and use a similar strategy as in [7] for integer orders to obtain this result. Proof. Sine B λ Eλ s (Z) and s > 3/2, the mapping ˇG = B t M B defines, for eah I s, a bounded isomorphism between H s (, R d ) and its dual H s (, R d ). Thus we need only to show that the mapping ˇG, I s L(H s (, R d ), H s (, R d )) is smooth. Sine transposition between Banah spaes is itself a bounded operator it follows that the transpose is smooth iff B t, I s L(L 2 (, R d ), H s (, R d )) B, I s L(H s (, R d ), L 2 (, R d )) is smooth, whih is the ase by Proposition 5.6. Now, the mapping M, I s L(L 2 (, R d ), L 2 (, R d )), is smooth for s > 3/2. Finally, sine omposition of bounded operators between Banah spaes is itself a bounded operator, it follows that the omposition is smooth. B t M B, I s L(H s (, R d ), H s (, R d )), Remark 7.3. A smooth, strong Riemannian metri on a Hilbert manifold has a smooth spray (see [25] for instane). However, formula (4.2) is no longer useful in that ase beause it is not lear that this expression extends to a smooth map from I s to H s (, R d ). Following Lang [25, Proposition 7.2], we introdue the H s inner produt on H s (, R d ): h, k H s := Λ s h, Λ s k dθ

20 20 MARTIN BAUER, MARTINS BRUVERIS, AND BORIS KOLEV so that G (h, k) = P h, k H s, where P := ( Λ s B ) t Λ s M B. Here, Q t is the transpose of a bounded operator Q : H s (, R d ) H s (, R d ). In that ase, there is an alternative expression for the spray, whih is nevertheless equivalent to (4.2). It is given impliitly by the formula S (h), P k H s = 1 2 (D,kP ) h, h H s (D,h P ) h, k H s. Appendix A. Smooth urves in S r (Z) The right framework to work with smooth maps on infinite dimensional (other than Banah spaes) are onvenient vetor spaes (see [19, 24]). Note that every Fréhet spae is a onvenient vetor spae. We will need the following lemma about reognising smooth urves in onvenient vetor spaes. Lemma A.1. [19, ] Let : R E be a urve in a onvenient vetor spae. Let V E be a point-separating subset of bounded linear funtionals suh that the bornology of E has a basis of σ(e, V)-losed sets. Then the following are equivalent: (1) is smooth; (2) There exist loally bounded urves k : R E suh that l : R R is smooth with (l ) (k) = l k, for eah l V. To apply this lemma to the spae S r (Z) of Fourier multipliers we need to hoose a suitable set V of linear funtionals. This is aomplished in the following lemma. Lemma A.2. Let r R and V = {λ ev m : m Z, λ L(C d ) } S r (Z). Then the bornology of S r (Z) has a basis onsisting of σ(s r (Z), V)-losed sets. Proof. We an regard V as a subset of both S r (Z) as well as l (Z, L(C d )). It is shown in [19, ] that the bornology of l (Z, L(C d )) has a basis of σ(l (Z, L(C d )), V)-losed sets. We embed S r (Z) as a losed subspae into ι : S r (Z) ( ) l (Z, L(C d )), (a(m)) m Z m (r α) α a(m) α N m Z, α N. Beause the bornology of the produt is the produt bornology it follows that a basis for the bornology on α N l (Z, L(C d )) is given by σ( α l, W)- losed sets, where W = α N V pr α and pr α denotes the anonial projetion onto the α-th fator. Here we note that on S r (Z) the set ι * (W) is ontained in the linear spae of V and hene the bornology of S r (Z) has a basis onsisting of σ(s r (Z), V)-losed sets. Appendix B. Square-root of an operator in E r λ (Z) Sine a positive definite Hermitian matrix has a unique positive square root, whih depends smoothly on its oeffiients, we an define formally the square root B λ = op ( a(λ, m) 1/2) of an element A λ = a λ (D) in the lass Eλ r (Z). In order to prove this we need the following lemma.

21 FRACTIONAL SOBOLEV METRICS ON SPACES OF IMMERSED CURVES 21 Lemma B.1. [5, Lemma 4.8] Let a, b, x L(C d ) be three matries satisfying bx + xb = a, with b Hermitian and positive definite. Then d x b a, 2 where denotes the Frobenius norm, i.e. x = tr xx *. The following lemma together with its proof is a generalisation of [5, Lemma 4.7] to our situation of one-parameter families of symbols. Lemma B.2. The positive square root of an operator in the lass Eλ r(z) belongs to the lass E r/2 λ (Z). Conversely, the square of an operator in the lass Eλ r (Z) belongs to the lass E2r λ (Z). Proof. We will prove the estimate β λ α b(λ, m) m r/2 α, whih holds loally uniformly in λ, by indution over α + β. If α + β = 0 the statement is b(λ, m) m r/2. Assume that it has been proven for α + β k. Then let α + β = k + 1 and, omitting the arguments (λ, m), we obtain using the produt rule, ( ) ( ) β λ α a = b β λ α b + β λ α b b + ( ) ( ) α β j i j λ i b β j λ α i b, (i,j) X with X = {(i, j) : 0 i α, 0 j β, i + j < α + β}. Then by the indution assumption j λ i b β j α i b m r/2 i m r/2 α+i m r α, λ and hene we obtain via Lemma B.1, β λ α b m r/2 m r α m r/2 α. This ompletes the indution. Referenes [1] M. Bauer, M. Bruveris, P. Harms, and P. W. Mihor. Vanishing geodesi distane for the Riemannian metri with geodesi equation the KdV-equation. Ann. Global Anal. Geom., 41(4): , [2] M. Bauer, M. Bruveris, P. Harms, and J. Møller-Andersen. A numerial framework for Sobolev metris on the spae of urves. SIAM J. Imaging Si., 10(1):47 73, [3] M. Bauer, M. Bruveris, and P. W. Mihor. Overview of the geometries of shape spaes and diffeomorphism groups. J. Math. Imaging Vision, 50(1-2):60 97, [4] M. Bauer, M. Bruveris, and P. W. Mihor. Why use Sobolev metris on the spae of urves. In Riemannian omputing in omputer vision, pages Springer, Cham, [5] M. Bauer, J. Esher, and B. Kolev. Loal and global well-posedness of the frational order EPDiff equation on R d. Journal of Differential Equations, 258(6): , [6] M. Bauer, P. Harms, and P. W. Mihor. Sobolev metris on shape spae of surfaes. J. Geom. Meh., 3(4): , 2011.

22 22 MARTIN BAUER, MARTINS BRUVERIS, AND BORIS KOLEV [7] M. Bruveris. Completeness properties of Sobolev metris on the spae of urves. J. Geom. Meh., 7(2): , [8] M. Bruveris. Regularity of maps between Sobolev spaes. Ann. Glob. Anal. Geom., pages 1 14, Online first. [9] M. Bruveris, P. W. Mihor, and D. Mumford. Geodesi ompleteness for Sobolev metris on the spae of immersed plane urves. Forum Math. Sigma, 2:e19, 38, [10] M. Bruveris and F.-X. Vialard. On ompleteness of groups of diffeomorphisms To appear in J. Eur. Math. So. [11] R. Camassa and D. D. Holm. An integrable shallow water equation with peaked solitons. Phys. Rev. Lett., 71(11): , [12] V. Cervera, F. Masaró, and P. W. Mihor. The ation of the diffeomorphism group on the spae of immersions. Differential Geom. Appl., 1(4): , [13] A. Constantin and B. Kolev. Geodesi flow on the diffeomorphism group of the irle. Comment. Math. Helv., 78(4): , [14] P. Constantin, P. D. Lax, and A. Majda. A simple one-dimensional model for the three-dimensional vortiity equation. Comm. Pure Appl. Math., 38(6): , [15] I. L. Dryden and K. Mardia. Statistial Shape Analysis. John Wiley & Son, [16] D. G. Ebin and J. E. Marsden. Groups of diffeomorphisms and the motion of an inompressible fluid. Ann. of Math. (2), 92: , [17] J. Esher and B. Kolev. Right-invariant Sobolev metris of frational order on the diffeomorphism group of the irle. J. Geom. Meh., 6(3): , [18] J. Esher, B. Kolev, and M. Wunsh. The geometry of a vortiity model equation. Commun. Pure Appl. Anal., 11(4): , July [19] A. Fröliher and A. Kriegl. Linear Spaes and Differentiation Theory. Pure and Applied Mathematis (New York). John Wiley & Sons, Ltd., Chihester, A Wiley- Intersiene Publiation. [20] F. Gay-Balmaz and T. S. Ratiu. The geometry of the universal Teihmüller spae and the Euler-Weil-Petersson equation. Adv. Math., 279: , [21] D. G. Kendall. Shape manifolds, Prorustean metris, and omplex projetive spaes. Bull. London Math. So., 16(2):81 121, [22] E. Klassen, A. Srivastava, M. Mio, and S. Joshi. Analysis of planar shapes using geodesi paths on shape spaes. IEEE T. Pattern Anal., 26(3): , [23] S. Kouranbaeva. The Camassa-Holm equation as a geodesi flow on the diffeomorphism group. J. Math. Phys., 40(2): , [24] A. Kriegl and P. W. Mihor. The Convenient Setting of Global Analysis, volume 53 of Mathematial Surveys and Monographs. Amerian Mathematial Soiety, Providene, RI, [25] S. Lang. Fundamentals of Differential Geometry, volume 191 of Graduate Texts in Mathematis. Springer-Verlag, New York, [26] H. L. Le and D. G. Kendall. The Riemannian struture of Eulidean shape spaes: a novel environment for statistis. Annals of Statistis, 21(3): , [27] A. C. Mennui, A. Yezzi, and G. Sundaramoorthi. Properties of Sobolev-type metris in the spae of urves. Interfaes Free Bound., 10(4): , [28] P. W. Mihor and D. Mumford. Vanishing geodesi distane on spaes of submanifolds and diffeomorphisms. Do. Math., 10: (eletroni), [29] P. W. Mihor and D. Mumford. Riemannian geometries on spaes of plane urves. J. Eur. Math. So., 8:1 48, [30] P. W. Mihor and D. Mumford. An overview of the Riemannian metris on spaes of urves using the Hamiltonian approah. Appl. Comput. Harmon. Anal., 23(1):74 113, [31] G. Misio lek and S. C. Preston. Fredholm properties of Riemannian exponential maps on diffeomorphism groups. Invent. Math., 179(1): , [32] M. Ruzhansky and V. Turunen. Pseudo-differential operators and symmetries, volume 2 of Pseudo-Differential Operators. Theory and Appliations. Birkhäuser Verlag, Basel, Bakground analysis and advaned topis. [33] S. Shkoller. Geometry and urvature of diffeomorphism groups with H 1 metri and mean hydrodynamis. J. Funt. Anal., 160(1): , 1998.

23 FRACTIONAL SOBOLEV METRICS ON SPACES OF IMMERSED CURVES 23 [34] S. Shkoller. Analysis on groups of diffeomorphisms of manifolds with boundary and the averaged motion of a fluid. J. Differential Geom., 55(1): , [35] A. Srivastava and E. Klassen. Funtional and Shape Data Analysis. Springer Series in Statistis, [36] A. Srivastava, E. Klassen, S. H. Joshi, and I. H. Jermyn. Shape analysis of elasti urves in Eulidean spaes. IEEE T. Pattern Anal., 33(7): , [37] G. Sundaramoorthi, A. C. Mennui, S. Soatto, and A. Yezzi. A new geometri metri in the spae of urves, and appliations to traking deforming objets by predition and filtering. SIAM J. Imaging Si., 4(1): , [38] G. Sundaramoorthi, A. Yezzi, and A. C. Mennui. Sobolev ative ontours. International Journal of Computer Vision, 73: , [39] G. Sundaramoorthi, A. Yezzi, and A. C. Mennui. Coarse-to-fine segmentation and traking using Sobolev ative ontours. IEEE Transations on Pattern Analysis and Mahine Intelligene, 30(5): , [40] M. Wunsh. On the geodesi flow on the group of diffeomorphisms of the irle with a frational Sobolev right-invariant metri. J. Nonlinear Math. Phys., 17(1):7 11, [41] L. Younes. Computable elasti distanes between shapes. SIAM J. Appl. Math., 58(2): (eletroni), [42] L. Younes. Shapes and Diffeomorphisms, volume 171 of Applied Mathematial Sienes. Springer-Verlag, Berlin, [43] L. Younes, P. W. Mihor, J. Shah, and D. Mumford. A metri on shape spae with expliit geodesis. Atti Aad. Naz. Linei Cl. Si. Fis. Mat. Natur., 19(1):25 57, Faulty for Mathematis, Florida State University, USA address: bauer@math.fsu.edu Department of Mathematis, Brunel University London, Uxbridge, UB8 3PH, United Kingdom address: martins.bruveris@brunel.a.uk Aix Marseille Université, CNRS, Centrale Marseille, I2M, UMR 7373, Marseille, Frane address: boris.kolev@math.nrs.fr