Hig order approximation to non-smoot multivariate functions Anat Amir David Levin arxiv:164.281v1 [mat.na] 11 Apr 216 April 12, 216 Abstract Approximations of non-smoot multivariate functions return low-order approximations in te vicinities of te singularities. Most prior works solve tis problem for univariate functions. In tis work we introduce a metod for approximating non-smoot multivariate functions of te form f = g + r + were g, r C M+1 R n and te function r + is defined b { r, r r + =, R n., r < Given scattered or uniform data points X R n, we investigate approximation b quasi-interpolation. We design a correction term, suc tat te corrected approximation acieves full approximation order on te entire domain. We also sow tat te correction term is te solution to a Moving Least Squares MLS problem, and as suc can bot be easil computed and is smoot. Last, we prove tat te suggested metod includes a ig-order approximation to te locations of te singularities. 1 Introduction Approximation of non-smoot functions is a complicated problem. Most approximation tools return smoot approximations, tus reling on te smootness of te original function for te approximation to be correct. However, te need to approximate non-smoot functions does exist in man applications. For a ig-order approximation of non-smoot functions, we need to allow our approximation to be non-smoot. Oterwise, in te vicinities of te singularities, we will get a low-order approximation. In tis work we will suggest a metod tat will allow us to properl approximate non-smoot functions of a given model. We will concentrate on functions f : R n R wic ma be modelled as f = g + r + were g, r C M+1 R n and te function r + is defined b r + = { r, r, r <, R n. Israel Scool of Matematical Sciences, Tel Aviv Universit, Ramat Aviv, Tel Aviv 699781, 1
Suc functions are obviousl continuous, but are non-smoot across te persurface Γ r := {z R n : rz = }. As an example for suc functions, consider sock waves, wic are solutions of non-linear perbolic PDEs [1]. Anoter example migt be a signed distance function [11], were te distance is measured from a disconnected set. Our goal is to acieve ig-order approximations of suc functions. To acieve tat we will concentrate on a specific famil of approximation tools. Consider a quasi-interpolation operator Q [14]. Suc an operator receives te values of a function φ : R n R on a set of data points X R n. Te quasi-interpolation operator Q returns an approximation defined b Qφ := q x φx, R n, were {q x } are te quasi-interpolation basis functions, eac is smoot and as finite support. Let be te fill distance of X, satisfing: 1. := min {L : B, L X, R n }. 2. Tere exists N > suc tat R n we ave Denote # X B, n N. R 1 := inf {R > : suppq x Bx, R, x X}. We assume tat te operator Q as a bounded Lebesgue constant { } L 1 := sup q x : R n, and reproduces polnomials in Π M R n := {p : R n R : degp M}, were degp is te total degree of te polnomial p. Ten, te error in te quasi-interpolation, Eφ := φ Qφ satisfies for all φ C M+1 R n and R n were and Eφ C 1 φ C M+1 M+1 C 1 = 1 + L 1 R 1 M+1 D J φ φ C M+1 := max J =M+1 J! wit J a multi-index and te maximum norm. Tat is, te operator Q as full approximation order for smoot functions [14]. However, as explained above, since te function Qφ is alwas smoot, te operator gives low-order approximations in te vicinities of singularities. 2
One example of a quasi-interpolation operator is te MLS approximation [7]. Given a function φ : R n R and a point R n te MLS approximation is defined as Qφ := p were x p := arg min px φx 2. p Π M R n Here is a smoot weigt function wit compact support. Te MLS approximation is especiall important in tis work, since we base muc of our results on our abilit to adapt te MLS approximation to our non-smoot scenario. Our initial goal for tis work was to generalize te work done b Lipman and Levin [8]. In tat work, te autors address te problem of approximation of univariate functions of te form fx = gx + M j=1 j j! x sj +, { were g P C M+1 R and x s j x s + = j, x s. Indeed, suc, x < s functions are continuous but not smoot. In [8], te univarite case is solved b modelling te error terms of te approximation b a quasi-interpolate Q. Tat is, one searces for variables s,, suc tat te errors in te quasiinterpolation approximation of te term r,s x := M j=1 j j! x s j + gives te best Least-squares approximation to te errors of te function f at te data points. It is sown tat b adding te error of te approximation of te new term r,s to te approximation Qf, full approximation order for te function f is acieved. Anoter approac to tis problem was proposed b Harten [5]. Te autor introduces te essentiall non-oscillator ENO and te subcell resolution SR scemes. Te ENO sceme bases te approximation at eac point on onl some of te data points in its vicinit. Tus, disregarding points from te oter side of te singularit wic contaminate te approximation. Te SR sceme locates te singularities b intersecting polnomials from supposedl different sides of te singularities. For an examination of tese metods for univariate functions wit a jump discontinuit in te derivative see [1]. Oter approaces were suggested b Markakis and Barack [9], were te autors revise te Lagrange interpolation formula to approximate univariate discontinuous functions, and b Plaskota et al [12], [13], were te autors suggest using adaptive metods for tis approximation. Batenkov et al [2], [3], [4] address a similar problem, te reconstruction of a piecewise smoot function from its integral measurements. One disadvantage of te metods mentioned above, is tat te do not easil adapt to multivariate singular functions. Te main advantage of te metod we suggest in tis paper is its abilit to deal wit te multivariate case. Our metod enables us to approximate multivariate functions wic ave non-continuous derivatives across smoot persurfaces, Γ r. Indeed, te dimension n of te domain of te function f affects 3
te required number of data points in X and te dimension of Π M R n, but te correction procedure is not affected b te iger dimension. 2 Main results Our goal is to fix te approximation of te function f = g + r +. Remark 1 Te decomposition f = g + r + is not unique. Indeed, f = g + r + = g + r + r +. However, we onl correct te approximation error, for wic, as we sow later, we do ave uniqueness. We will acieve tis, following te main idea in [8], b investigating te error terms Ef = Eg +Er + = Er + + O M+1. }{{} O M+1 Tis work relies eavil on our abilit to determine weter te function r returns a positive or negative value at eac data point x X. For tis we will need te following definition: Definition 1 For a set X R n wit fill-distance and a function f we will sa tat te set P X determines te signs of te nodes in X wit respect to f if 1. x P eiter f x > or f x = O M+1. 2. x X, if f x > ten x P. In section 5 we will introduce an algoritm tat determines te signs of te data points in X wit respect to te function r. For now, let us assume tat a set P wic determines te signs of te data points in X wit respect to te function r is known. Remark 2 Apparentl, once we ave a set P wic determines te signs of te nodes in X wit respect to te function r, we can approximate a point were r as a positive value using onl te data points in P, and a point were r as negative value using onl te data points in X \ P. However, we predict weter r as a positive or negative value onl on te data points X, and not on te entire domain. Specificall, for a point close to te singularit location, Γ r, we can not tell weter we sould approximate based on P or on X \ P. Hence we can not rel onl on te set P to fix te approximation. Definition 2 For R n define p min p min := arg min p Π M R n x ω Π M R n b Eχ r px Efx 2. Here ω : R + R + is a smoot positive weigt function wit finite support, suppω = [, ρ] 4
and χ r is te indicator function defined b { 1, z P χ r z =, z / P, z Rn. Remark 3 Peraps a more natural coice for te polnomial p min been te polnomial minimizing te sum x ω Ep + x Efx 2. would ave However, tis computation is not linear, ence we defined te polnomial p min as we did. Let {p α } α I be a basis of omogeneous polnomials for te linear space Π M R n. And denote { } z D = + z : ω > = B, ρ. In te following we state te main results of tis work. Te proofs of tese teorems are given in section 3. Note tat tese results rel upon te linear independence of te vectors {Eχ r p α x : x D X} α I. In section 4 we discuss tis issue, and suppl a proof for te linear independence of te vectors for an MLS approximation. Wile we did not succeed in proving tis linear independence for a general quasi-interpolation operator, we believe tat tese vectors are indeed independent except for degenerate cases. Teorem 1 Approximation of te singular part For eac R n denote p min = λ α p α. α I Ten we ave te following: 1. Te mappings λ α are smoot. 2. Tere exists a constant C 2 > suc tat if te vectors are linearl independent, ten {Eχ r p α x : x D X} α I p min z rz C 2 g C M+1 + r C M+1 M+1 for all z B, R 1. Remark 4 As will be explained later in te proof of Teorem 1, for te computation of te polnomial p min we ave to find te solution δ to te linear sstem ADA T δ = AD φ were A is te matrix A = E χ r p α x α,x 5,
D is te diagonal matrix wit values x D x,x = ω, and φ is te vector φ x = Efx. If te vectors {Eχ r p α x : x D X} α I are linearl dependent ten te matrix ADA T is singular, and te sstem does not ave a unique solution. In tis case we suggest using te pseudo-inverse of tis matrix, tus guaranteeing tat te solution we coose is te one of minimal norm. For a point R n far enoug from te singularit location, Γ r, te errors satisf Eχ r p α x =, x D X and α I. Tus, te vectors are clearl dependent, and b using te pseudo-inverse we guarantee tat te polnomial p min we coose will be te zero polnomial. Te corrected approximation Denote b G R n te set of all points R n for wic te vectors {Eχ r p α x : x D X} α I are linearl independent. For a point R n follow tese steps: 1. Find p min p min defined b := arg min p Π M R n x ω 2. Define te corrected approximation Eχ r px Efx 2. Qf := Qf + { E p min +, G, oterwise. Tus we get Teorem 2 Corrected approximation errors Tere exists a constant C 3 > suc tat if G ten Êf = f Qf C 3 g C M+1 + r C M+1 M+1. Remark 5 Te complexit of te entire approximation procedure for a point R n is ON 2 + J 2 N + J 3 + J N + 1 + N + 1, were J = dimπ M R n and N := min {n N : R n, n # {x X : q x } D X}. 6
We first ave to compute te quasi-interpolation errors at te data points in te vicinit of, D X, at ON 2. Afterwards, we solve te linear sstem ADA T δ = AD φ at OJ 2 N to build te sstem, and OJ 3 to solve te sstem. Last, we compute te polnomial p min at te data points in te vicinit of and at itself, at OJ N + 1, and set te approximation to be at ON + 1. Qf + Ep min + Teorem 3 Approximation of te singularit location Γ r Assume tat rz, z Γ r. Define r : G R b r = p min. If { C1 := { G : r = } C 2 := { G : r = } ten d H C 1, C 2 = O M+1, were d H is te Hausdorff distance of te two sets., 3 Proofs 3.1 Proof of Teorem 1 Proof 1 Let us define an MLS operator Q for a function φ : R n R n b Qφ = p φ, R and were p φ, := arg min p Π M R n x ω Eχ r px Eχ r φx 2. Write p = α I δ α p α and p φ, = α I δ φ, α p α, ten from te linearit of E we ave δ φ, = arg min δ x ω α δ α Eχ r p α x Eχ r φx 2. Recall tat te polnomial p α is a omogenous polnomial for eac α I, and denote d α := degp α, 7
ten we ma equivalentl solve te problem x δ φ, = arg min ω δ 2 δ α dα E χ r p α x Eχ r φx. α To solve tis problem we follow [7], from wic we know tat if te vectors { } E χ r p α x : x D X α I are linearl independent ten te solution to tis problem is te vector δ φ, defined b δ φ, = ADA T 1 AD φ were A is te matrix A = E χ r p α x α,x D is te diagonal matrix wit values x D x,x = ω,, and φ is te vector φ x = Eχ r φx. Moreover, te operator Q clearl reproduces polnomials in Π M R n and as suc is a quasi-interpolation operator Qφ = q x Eχ r φx wit basis functions of finite support q x := dα p α ADA T 1 AD = α,x α I p α ADA T 1 AD α,x α I R 2 := inf {R : supp q x Bx, R, x X} and a bounded Lebesgue constant { } L 2 := sup q x z : z R n. First note tat A α,x = E χ r p α x L 1 R1 M. 8
Assume tat te polnomials {p α } are eac of te form p α u = u z α, wit α a multi-index. Ten we ave q x z ADA T 1 L 1 R M 1 z x ω, wit Hence, were Note tat { ADA ADA T 1 T 1 v = sup : v } v p min L 2 ADA T 1 L 1 R M 1 K, { } z x K := sup ω : z R n. z = p α z ADA T 1 AD Efx. α,x α I. Also, since te operator Q reproduces polnomials we ave tat p T alor z = p α z ADA T 1 AD Eχ α,x r p T alor x, α I were p T alor is te Talor approximation of te function r at te point. Denote #I = dimπ M R n. Ten, p min Compute, z p T alor z = α I p α z ADA T 1 AD Efx Eχ α,x r p T alor x #I z M q x Efx Eχ r p T alor x #I z M { L 2 max Efx Eχr p T alor x } B,R 2 Efx Eχ r p T alor x = Egx + Er + x Eχ r p T alor x Egx + Er + x Eχ r p T alor x = Egx + Eχ r rx Eχ r p T alor x = Egx + Eχ r r p T alor x C 1 g C M+1 M+1 + 1 + L 1 r C M+1 x M+1 9
Hence for all z B, R 1 we ave Clearl for C 2 = max p min z rz p min z p T alor z + p T alor z rz #IR 1 M L 2 C1 g C M+1 + 1 + L 1 r C M+1R 2 M+1 2M+1 + R 1 M+1 r C M+1 M+1 { #I ADA T 1 L 1 KC 1 R 1 M, #I ADA T 1 1 + L 1 L 1 KR 2 R 1 R 2 M + R M+1 1 we get for all z B, R 1 tat p min z rz C 2 g C M+1 + r C M+1 M+1. Moreover, from te representation p min = α I p α ADA T 1 AD Efx, α,x } {{ } λ α it is clear tat te mappings of te coefficients λ α are smoot. 3.2 Proof of Teorem 2 Proof 2 If G, ten te vectors } {Eχ r p α x : x D} α I are linearl independent and from Teorem 1 we ave p min z rz C 2 g C M+1 + r C M+1 M+1, for all z B, R 1. Consequentl, we get, Er + Ep min + 1 + L 1 C 2 g C M+1 + r C M+1 M+1. Terefore, Êf = f Qf = f Qf E p min + }{{} Ef = Ef E p min + Er + E p min + + Eg 1 + L 1 C 2 g C M+1 + r C M+1 M+1 + C 1 g C M+1 M+1. 1
3.3 Proof of Teorem 3 Proof 3 Pick C 1, ten r = and r. Using Talor s approximation we ave Hence, r ± ɛ r = ±ɛ r 2 + Oɛ 2 r 2. r ɛ r < < r + ɛ r. From Teorem 1 we know tat r is a smoot function and tat r ɛ r < < r + ɛ r. Hence, tere must exist u = + λ r wit λ < ɛ suc tat ru =. But ten, using Teorem 1 and Talor s approximation again we get Hence, and consequentl O M+1 = ru = λ r 2 + Oλ 2 r 2. λ r 2 = O M+1, u = λ r = O M+1. Tat is, tere must exist u wit u = O M+1 and ru =, ence, Similarl we sow tat sup inf d, u = O M+1. u C 2 C 1 sup inf d, u = O M+1, u C 1 C 2 and we ave d H C 1, C 2 = O M+1. 4 Linear independence of {Eχ r p α x : x D X} α I We ave based our correction algoritm at a point R n on te fact tat te vectors {Eχ r p α x : x D X} α I are linearl independent. We ave not succeeded in suppling a proof for tis linear independence for a general quasi-interpolation operator. However, if we use an MLS approximation as our quasi-interpolation operator tese vectors are indeed independent. Let us assume tat our quasi-interpolation operator is an MLS approximation, tat is for a function φ : R n R and a point R n we ave Qφ := p were p := arg min x px φx 2. p Π M R n 11
Here is our weigt function wit compact support [, ϱ]. Tat an MLS approximation is indeed a quasi-interpolation operator is sown in [7]. Denote Υ := min {ɛ > : X B, ɛ is uni-solvent for Π M R n, R n }. We will assume tat 4Υ ϱ and tat B, 2Υ D = B, ρ, R n. Teorem 4 Linear independence Assume tat rz for all z Γ r. Ten for R n, if bot P D and X \P D are non-empt, te vectors are linearl independent. Proof 4 Pick R n. If te vectors {Eχ r p α x : x D X} α I {Eχ r p α x : x D X} α I are linearl dependent ten tere exists a non-zero polnomial p Π M R n suc tat Eχ r px =, x D X. Since bot P D and X \P D are not empt, te intersection Γ r D must be non-empt, ence we ma define Note tat z := arg min u. u Γ r D z d, z + Υ = z Υ, z Hence we ave from B, 2Υ D = B, ρ tat z B z + Υ z, Υ D, and from te uni-solvenc of X B exists u X B z + Υ z + Υ z z, Υ wit pu. Note tat rz d z ± Υ rz, u rz d z ± Υ rz, z + Υ z z }{{} 2Υ Since 4Υ ϱ, it follows tat rz B z + Υ rz, Υ 12 z z, Υ for Π M R n tere X D + d z z + Υ z, u 3Υ }{{} Υ X Bu, ϱ P
and rz B z Υ rz, Υ X Bu, ϱ X \ P. Hence, bot Bu, ϱ P and Bu, ϱ X \ P are uni-solvent for Π M R n. r > ϱ ρ u z Υ r < Γ r Compute, Eχ r pu = χ r pu Qχ r pu = χ r pu p u u, wit u x p u = arg min p Π M R n u x arg min p Π M R n px χ r px 2 = px px 2 + \P u x p 2 x Define K : Π M R n \ {} R + b Kp := \P u x p 2 x u x p 2 x Note tat since P Bu, ϱ and X \P B, ϱ are uni-solvent for Π M R n, ten if p is not te zero polnomial we get u x u x p 2 x > and p 2 x >, \P ence te constant Kp is well-defined. Assume tat p u is not te zero polnomial, and denote K 1 = Kp u, ten. 13
we ave p u = arg min p Π M R n Kp=K 1 = arg min u x p Π M R n Kp=K 1 = arg min p Π M R n Kp=K 1 = arg min p Π M R n Kp=K 1 u x u x u x px px 2 + \P px px 2 + K 1 u x p 2 x [ ] px px 2 + K 1 p 2 x [ 1 + K 1 px px 1 + K 1 Clearl, te minimum is attained for a polnomial p satisfing px = px 1 + K 1, x P. Since P Bu, ϱ is uni-solvent for Π M R n we must ave p u = 1 1 + K 1 p. u x p 2 x ] 2 + K 1 p 2 x 1 + K 1 Terefore tere exists a constant K 2 suc tat p u = K 2 p. Note tat if p u is te zero polnomial a scenario we disregarded earlier we take K 2 =. Ten we ma write, K 2 = arg min c arg min c arg min c arg min c p ṷ p = u x u x u x u x c px px 2 + c 1 px 2 + p 2 x c 2 2 p 2 x = u x p x 2 c u x p 2 x u x \P u x \P u x u x u x u x u x 2 p 2 x p 2 x c px 2 = c px 2 = p 2 x c + p 2 x p 2 x 2 + 14
Here we see tat te minimum is attained for u x p 2 x c =. p 2 x u x Recall tat p is not te zero polnomial. Ten, from te uni-solvenc of P Bu, ϱ and X \ P Bu, ϱ for Π M R n we ave tat < u x p 2 x < u x p 2 x, and Tus, a contradiction to for u D X. < c < 1. Eχ r pu = χ r pu p u u = χ r c pu, Eχ r px =, x D X, 5 Sign determination algoritm Recall tat our metod relies upon our abilit to correctl identif a set P X wic determines te signs of te data points in X wit respect to te function r. Tat is, 1. For all x P eiter rx > or rx = O M+1. 2. If rx > ten x P. To acieve tis goal we rel on an important propert of te quasi-interpolation approximation. In te immediate vicinit of a singularit of a function φ te quasi-interpolation operator Qφ as large errors. On te oter and, if a point is distanced from te singularit b R 1 ten te approximation error at tis point can be bounded b C 1 φ C M+1 M+1. Note tat weter or not a quasi-interpolation operator as significant errors in te proximit of a singularit, is a question tat is ver muc related to te question of te linear independence of te vectors {Eχ r p α x : x D X} α I. Pick z Γ r, and let p T alor be te M-t Talor approximation to te function r centred at z, ten for a point x Dz X we ave Efx = Egx + Er + x = Eχ r p T alor x + O M+1. If te vectors {Eχ r p α x : x Dz X} α I are linearl independent, ten we sould ave a significant error at some data point x Dz X. Specificall, 15
for an MLS quasi-interpolation, we ave explained in section 4 w for a nonzero polnomial te approximation error is significant in te immediate vicinit of te singularit. Tus, we suggest determining te sign of te function r, based on an MLS quasi-interpolation approximation, or an oter quasi-interpolation operator for wic we can guarantee significant errors near te singularit, b performing te following steps: 1. Compute te approximation errors on te data set X of an MLS approximation. 2. Run a cluster algoritm on te approximation errors of X to divide te data set into two subsets, one wit small errors N and te oter wit significant errors S. 3. Define M := N \ Bx, R 1. x S Since te singularit Γ r divides R n into two connected sets, te union Bx, x M must also ave two connected sets. Arbitraril pick one of tese connected sets and denote it b P. 4. For eac data point in te set X Bx, R 1 x S approximate its value using MLS approximation based on eac of te connected components separatel, if te approximation error is smaller wen basing te approximation on te connected component P add tis point to te set P. It is important to note tat since we arbitraril coose te set P, we ma pick te wrong connected set. Tat is, instead of coosing te connected set on wic te function r returns onl positive values, we coose te connected set on wic te function r returns negative values. However, tis coice will ave no effect on te approximation algoritm. Indeed, Er + = Er r = }{{} Er +E r = O + E r +, }{{} O M+1 r + were E is te approximation error of te quasi-interpolation operator. Tus, te actual coice of te connected set is insignificant. 16
6 Numerical results For our test we took our data set X to be a random set of 41 2 data points wit fill distance.33 in te region [.4,.4] 2 R 2. Using te sampled data points, we approximated a function, and compared our approximation to te actual function values. We performed tis comparison on a uniform mes of 6 2 points wit fill distance =.1 in te region [.3,.3] 2 R 2. Note tat we avoided testing points close to te edge, wising to avoid approximation errors resulting from partial neigbourood near te edges. For te quasi-interpolation we used te MLS approximation wit weigt function x ω = e x 2 4 2. We approximated te functions f 1 = g + r 1 + and f 2 = g + r 2 + were gx, = e x+ r 1 x, = x 2 + 2 1 2 5 r 2 x, = x 4 + 4 1 4. 5 In figure 1, one can see a comparison between te errors of te original MLS approximation and our corrected approximation for bot f 1 and f 2, ere we used M = 4. In table 1 and figure 2, tere is a comparison of te errors for varing M, te total degree of te approximating polnomials. Table 1: Comparison of te errors for varing M. M Ef 1 Êf 1 Ef 2 Êf 2 1 1.6 1.22 1.31 1.35 2 1.54 2.42 2.1 1.9 3 1.57 3.2 2.23 2.56 4 1.6 4.26 2.42 4.22 5 1.66 5.6 2.48 5. 6 1.67 6.23 2.47 6.19 In figure 3, we drew te curves r 1 = and r 2 =, on wic te functions f 1 and f 2 respectivel ave singularities, and approximations to tese curves 17
1-3 1-3 4 4 2 2-2 -2-4 -4-6.4-6.4.2.4.2 -.2 -.2 -.4 -.4.2.4.2 -.2 -.2 -.4 -.4 1-4 1-4 4 4 2 2-2 -2-4 -4-6.4-6.4.2.4.2 -.2 -.2 -.4 -.4.2.4.2 -.2 -.2 -.4 -.4 Figure 1: Comparison between te original MLS errors left and te errors of our corrected approximation rigt for f 1 above and f 2 below. -1-1 -2-2 -3-3 -4-4 log E -5 log E -5-6 -6-7 -7-8 Original Error Corrected Error -8 Original Error Corrected Error -9 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 M -9 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 M Figure 2: Te errors of te original MLS approximation and te corrected approximation for varing M values for f 1 left and f 2 rigt. 18
.3 r= p min =.3 r= p min =.2.2.1.1 -.1 -.1 -.2 -.2 -.3 -.3 -.4 -.4 -.3 -.2 -.1.1.2.3 -.4 -.4 -.3 -.2 -.1.1.2.3 Figure 3: Comparison between te curve r 1 = and r 2 = red and te curve p min = blue dotted for f 1 left and f 2 rigt. we derive from our corrected approximations, based on te values of p min. We ave used te MATLAB contour command to draw tese curves. Here again we used M = 4. In figure 4, we outline te Hausdorff distance between te singularit curve r 1 = and our approximation of tis curve, for varing values of fill-distance and constant M = 4. Note tat as we took varing values we did not cange te number of data points, but onl te size of te region in wic we test te procedure. In figure 5, we compare te original approximation errors to te corrected approximation errors for varing values and constant M = 4 for bot f 1 and f 2. We also tested our sign determination algoritm on our data set for te functions f 1 and f 2. In figure 6 we sow te two clusters a k-means algoritm [6] clusters te original approximation errors into. In section 5 we ave denoted te cluster wit small errors b N and te cluster wit significant errors b S. In figure 7 we sow te data points wic lie witin a fixed-distance neigbourood of S, and te data points in N wic do not lie witin tis neigbourood. In section 5 we ave denoted te later set b M. Tese are te data points we know are distanced from te singularit curve and as suc, teir approximation is not influenced b te singularit. Last, as explained in section 5, we test eac of te data points wic are not in M to see weter te are better approximated b te connected component P or te oter connected component. If te better approximation is b te set P, ten we add te data point to te set P. Te final set P is our sign determining set wit respect to te functions r 1 and r 2. In figure 8 we sow our final set P, for bot f 1 and f 2, on wic we estimate te functions r 1 and r 2 respectivel return onl positive values. In bot test cases our estimation was correct for all data points. References [1] F. Arandiga, A. Coen, R. Donat, and N. Dn. Interpolation and approximation of piecewise smoot functions. SIAM Journal on Numerical Analsis, 431:41 57, 25. 19
-2.5-3 -3.5 log distance -4-4.5-5 -5.5-6 -6.5-2.4-2.2-2 -1.8-1.6-1.4-1.2-1 -.8 log Figure 4: Te Hausdorff distance of te curve r 1 = and our approximation of tis curve p min = for varing values of. -1-2 -2-3 -3-4 -4 log E -5-6 log E -5-6 -7-7 -8-9 Original Error Corrected Error -8 Original Error Corrected Error -1-3.5-3 -2.5-2 -1.5-1 -.5 log -9-3.5-3 -2.5-2 -1.5-1 -.5 log Figure 5: Te errors of te original MLS approximation and te corrected approximation for varing values for f 1 left and f 2 rigt. 2
.4.4.3.3.2.2.1.1 -.1 -.1 -.2 -.2 -.3 -.3 -.4 -.4 -.3 -.2 -.1.1.2.3.4 -.4 -.4 -.3 -.2 -.1.1.2.3.4.4.4.3.3.2.2.1.1 -.1 -.1 -.2 -.2 -.3 -.3 -.4 -.4 -.3 -.2 -.1.1.2.3.4 -.4 -.4 -.3 -.2 -.1.1.2.3.4 Figure 6: Te two clusters of data points for f 1 left and f 2 rigt. Te clusters of small approximation errors are above, and tose wit significant approximation errors are below. 21
.4.4.3.3.2.2.1.1 -.1 -.1 -.2 -.2 -.3 -.3 -.4 -.4 -.3 -.2 -.1.1.2.3.4 -.4 -.4 -.3 -.2 -.1.1.2.3.4.4.4.3.3.2.2.1.1 -.1 -.1 -.2 -.2 -.3 -.3 -.4 -.4 -.3 -.2 -.1.1.2.3.4 -.4 -.4 -.3 -.2 -.1.1.2.3.4 Figure 7: Te data points wic we predict are close to te singularit curve above and te data points wic are far from te curve below. On te left te singularit curve is r 1 =, and on te rigt te singularit curve is r 2 =..4.4.3.3.2.2.1.1 -.1 -.1 -.2 -.2 -.3 -.3 -.4 -.4 -.3 -.2 -.1.1.2.3.4 -.4 -.4 -.3 -.2 -.1.1.2.3.4 Figure 8: Our estimation of te data points on wic te function r 1 returns positive values left and on wic te function r 2 returns positive values rigt. 22
[2] D. Batenkov. Complete algebraic reconstruction of piecewise-smoot functions from fourier data. arxiv preprint arxiv:1211.68, 212. [3] D. Batenkov, N. Sarig, and Y. Yomdin. Algebraic reconstruction of piecewise-smoot functions from integral measurements. arxiv preprint arxiv:113.3969, 211. [4] D. Batenkov and Y. Yomdin. Algebraic fourier reconstruction of piecewise smoot functions. Matematics of Computation, 81277:277 318, 212. [5] A. Harten. ENO scemes wit subcell resolution. Journal of Computational Psics, 831:148 184, 1989. [6] A. K. Jain, M. N. Murt, and P. J. Flnn. Data clustering: A review. ACM Comput. Surv., 313:264 323, Sept. 1999. [7] D. Levin. Te approximation power of moving least-squares. Matematics of Computation of te American Matematical Societ, 67224:1517 1531, 1998. [8] Y. Lipman and D. Levin. Approximating piecewise-smoot functions. IMA Journal of Numerical Analsis, 34:1159 1183, 21. [9] C. Markakis and L. Barack. Hig-order difference and pseudospectral metods for discontinuous problems. arxiv preprint arxiv:146.4865, 214. [1] P. Morse and K. Ingard. Teoretical Acoustics. International series in pure and applied psics. Princeton Universit Press, 1986. [11] S. Oser and R. Fedkiw. Level Set Metods and Dnamic Implicit Surfaces. Applied Matematical Sciences. Springer, 23. [12] L. Plaskota and G. W. Wasilkowski. Te power of adaptive algoritms for functions wit singularities. Journal of fixed point teor and applications, 62:227 248, 29. [13] L. Plaskota, G. W. Wasilkowski, and Y. Zao. An adaptive algoritm for weigted approximation of singular functions over r. SIAM Journal on Numerical Analsis, 513:147 1493, 213. [14] H. Wendland. Scattered data approximation, volume 17. Cambridge universit press, 24. 23