arxiv: v1 [math.na] 27 Nov 2014
|
|
- Victor Hudson
- 5 years ago
- Views:
Transcription
1 arxiv:4.7589v [math.na] 7 Nov 4 Bivariate Lagrange interpolation at the node points of non-degenerate Lissajous curves Wolfgang Erb Christian Kaethner Mandy Ahlborg Thorsten M. Buzug September 6, 4 Abstract Motivated by an application in Magnetic Particle Imaging, we study bivariate Lagrange interpolation at the node points of Lissajous curves. The resulting theory is a generalization of the polynomial interpolation theory developed for a node set known as Padua points. With appropriately defined polynomial spaces, we will show that the node points of non-degenerate Lissajous curves allow unique interpolation and can be used for quadrature rules in the bivariate setting. An explicit formula for the Lagrange polynomials allows to compute the interpolating polynomial with a simple algorithmic scheme. Compared to the already established schemes of the Padua and Xu points, the numerical results for the proposed scheme show similar approximation errors and a similar growth of the Lebesgue constant. Keywords bivariate Lagrange interpolation quadrature formulas Chebyshev polynomials Lissajous curves Introduction A challenging task for multivariate polynomial interpolation is the construction of a suitable set of node points. Depending on the application, these node points should provide a series of favorable properties including a unique interpolation in given polynomial spaces, a slow growth of the Lebesgue constant and simple algorithmic schemes that compute the interpolating polynomial. The construction of suitable point sets for multivariate interpolation has a W. Erb, C. Kaethner, M. Ahlborg, T.M. Buzug Universität zu Lübeck Ratzeburger Allee Lübeck erb@math.uni-luebeck.de kaethner@imt.uni-luebeck.de ahlborg@imt.uni-luebeck.de buzug@imt.uni-luebeck.de
2 W. Erb et al. long-standing history. For an overview, we refer to the survey articles [3, 4] and the references therein. Examples of remarkable constructions in the bivariate setting are the point sets introduced by Morrow and Patterson [], Xu [4], as well as some generalizations of them [7]. A modification of the Morrow-Patterson points, introduced as Padua points [7], is particularly interesting for the purposes of this article. In some applications, the given data points are lying on subtrajectories of the euclidean space. In this case, aside from the above mentioned favorable properties, it is mandatory that the node points are part of these trajectories. Lissajous curves are particularly interesting examples for us, as they are used as a sampling path in a young medical imaging technology called Magnetic Particle Imaging (MPI) [5]. In MPI, the distribution of superparamagnetic iron oxide nanoparticles is reconstructed by measuring the magnetic response of the particles. The measurement process is based on the combination of various magnetic fields that generate and move a magnetic field free point through a region of interest. Although different trajectories are possible, this movement is typically performed in form of a Lissajous curve [9]. The reconstruction of the particle density from the data on the Lissajous trajectory is currently done in a very straight forward way, either by solving a system of linear equations based on a pre-measured system matrix or directly from the measurement data [6]. By using multivariate polynomial interpolation on the nodes of the sampling path, i.e. the Lissajous curve, we hope to obtain a further improvement in the reconstruction process. Of the node points mentioned above, the Padua points, as described in [3], are the ones with the strongest relation to Lissajous curves. They can be characterized as the node points of a particular degenerate Lissajous figure. Moreover, they satisfy a series of remarkable properties: they can be described as an affine variety of a polynomial ideal [5], they form a particular Chebyshev lattice [] and they allow a unique interpolation in the space Π n of bivariate polynomials of degree n [3]. Furthermore, a simple formula for the Lagrange polynomials is available and the Lebesgue constants are growing slowly as O ( log n ) [3]. The aim of this article is to develop, similar to the Padua points, an interpolation theory for node points on Lissajous curves. To this end, we extend the generating curve approach as presented in [3] to particular families of Lissajous curves in [, ]. In this article, we will focus on the node points of non-degenerate Lissajous curves, which are important for the application in MPI [8]. Not all of the above mentioned properties of the Padua points will be carried over to the node points of Lissajous figures. However, the resulting theory will have some interesting resemblences, not only to the theory of the Padua points, but also to the Xu points. We start our investigation by characterizing the node points Lisa n,p of nondegenerate Lissajous curves. Based on the node points Lisa n,p, we will derive suitable quadrature formulas for integration with product Chebyshev weight functions. Next, we will provide the main theoretical results on bivariate in-
3 Lagrange interpolation at Lissajous nodes 3 terpolation based on the points Lisa n,p. We will show that the points Lisa n,p allow unique interpolation in a properly defined space Π L n,p of bivariate polynomials. Further, we will derive a formula for the fundamental polynomials of Lagrange interpolation. This explicit formula allows to compute the interpolating polynomial with a simple algorithmic scheme similar to the one of the Padua points [8]. We conclude this article with some numerical tests for the new bivariate interpolating schemes. Compared to the established interpolating schemes of the Padua and Xu points, the novel interpolation schemes show similar approximation errors and a similar growth of the Lebesgue constant. The node points of non-degenerate Lissajous curves In this article, we consider π-periodic Lissajous curves of the form ( ) γ n,p : R R, γ n,p (t) = sin(nt), sin((n + p)t), () where n and p are positive integers such that n and n + p are relatively prime. Based on the calculations in [] (see also []), the Lissajous curve γ n,p is non-degenerate if and only if p is odd. In this case, γ n,p : [, π) R is an immersed plane curve with precisely n(n + p) n p self-intersection points. In the following, we will always assume that p is odd and sample the Lissajous curve γ n,p along the equidistant points t k := πk, k =,...,. In this way, we get the following set of Lissajous node points: Lisa n,p := γ n,p (t k ) : k =,...,. () To characterize the set Lisa n,p, we divide γ n,p (t k ) for the even and odd integers k. For this decomposition, we use the fact that n and n + p are relatively prime. Then, if n is odd, every odd integer k can be written as k = (i + )n + j(n + p) with i, j Z. If k is even, we can write k = in + (j +)(n+p) with i, j Z. If n is even, the same holds with the roles of n and n+p switched. In this way, we get the decomposition Lisa n,p = Lisa b n,p Lisa w n,p with the sets ) π ) ( (i + )n + j(n + p) Lisa b n,p := γ n,p ( in + (j + )(n + p) Lisa w n,p := γ n,p π : i, j Z : i, j Z, (3). (4) Two examples of Lissajous curves γ n,p with the corresponding node points Lisa b n,p and Lisa w n,p are illustrated in Figure. To get a compact representation
4 4 W. Erb et al (a) Lissajous figure γ,, Lisa, = (b) Lissajous figure γ,3, Lisa,3 = 7. Fig. Illustration of non-degenerate Lissajous curves γ n,p. The node points Lisa n,p of γ n,p are arranged on two different grids (black, white) corresponing to the sets Lisa b n,p and Lisa w n,p. of Lisa b n,p and Lisa w n,p, we use the following notation for the Gauß-Lobatto points: ( ) kπ zk n := cos, n N, k Z. (5) n Then, evaluating the points (3) and (4) explicitly for the Lissajous curve (), we get the following characterization: Lisa b n,p = ( ) i+j( ) z (n+p) (i+)p, i =,..., n + p zn jp :, (6) j =,..., n Lisa w n,p = ( ) i+j( ) z (n+p) ip, z(j+)p n i =,..., n + p :. (7) j =,..., n Since p is assumed to be odd and relatively prime to n, p is relatively prime to n as well as to (n + p). Therefore, by rearranging the points, we can drop the number p in the lower indices of the Gauß-Lobatto points in (6) and (7). Due to the point symmetry of the Lissajous curve γ n,p, the term ( ) i+j which preceeds the points in (6) and (7) can also be dropped by further rearrangement. This leads to the following simple characterization of the point sets Lisa b n,p and Lisa w n,p: ( Lisa b n,p = z (n+p) i + ( Lisa w n,p = z (n+p) i ), zn j : ), zn j + : i =,..., n + p j =,..., n, (8). (9) i =,..., n + p j =,..., n With this characterization, we can also divide the points Lisa n,p into the sets Lisa int n,p and Lisa out n,p denoting the points lying in the interior and on the
5 Lagrange interpolation at Lissajous nodes 5 boundary of the square [, ] respectively. We have ( ) Lisa int n,p := z (n+p) i +, i zn j : =,..., n + p j =,..., n ( ) z (n+p) i, i zn j + : =,..., n + p j, =,..., n ( ) Lisa out n,p := z (n+p) i +, ± : i =,..., n + p ( ) ±, zj n + : j =,..., n. From the representation of the Lisa points in (8) and (9), it is possible to count the number of points in the different sets. They are listed in Table. Table Cardinality of the different Lisa sets. Set Number of elements Lisa n,p n(n + p) + n + p Lisa b n,p (n + )(n + p) Lisa w n,p n(n + p + ) Lisa int n,p Lisa out n,p n(n + p) n p 4n + p From the representation in (3) and (4) and its identification in (6) and (7), we can deduce that ( ) ( ) (i + )n + j(n + p) (i + )n j(n + p) γ n,p π = γ n,p π, ) ) ( in + (j + )(n + p) γ n,p π ( in + (j + )(n + p) = γ n,p π holds for all i, j Z. Moreover, in (3) and (4) the boundary points are represented by j nz and i (n + p)z, respectively. Thus, for interior points in Lisa b n,p Lisa int n,p, i.e. all points in (3) satisfying j nz, there exist at least two different k, k in () that represent the same point. The same holds for all interior points in the second set Lisa w n,p. Therefore, all points in Lisa int n,p are self-intersection points of the Lissajous curve γ n,p. Since Lisa int n,p = n(n+p) n p corresponds to the total number of self-intersection points of a non-degenerate Lissajous curve (see []), we can conclude that Lisa int n,p is precisely the set of all self-intersection points of the Lissajous curve γ n,p. Finally, since Lisa int n,p + Lisa out n,p =, we can also conclude that there are exactly two different k, k that represent the same point in Lisa int n,p and that every point in Lisa out n,p is described by exactly one k in ().
6 6 W. Erb et al. In order to identify the different integers k in () that describe the same point A Lisa n,p, we introduce for k, k Z the equivalence relation k Lisan,p k γ n,p (t k ) = γ n,p (t k ). We say that k Z belongs to the equivalence class [A], A Lisa n,p, if γ n,p (t k ) = A. Therefore, by the above argumentation, there is exactly one k in the equivalence class [A] if A Lisa out n,p and exactly two if A Lisa int n,p. Remark There are some remarkable relations between the Lisa, Padua and Xu points. In formal terms, if p = in the characterization (8) and (9) of the Lisa points, the points Lisa n, correspond with the even Xu points XU n as defined in [4]. Moreover, if p = in (8) and (9), we obtain the even Padua points PD n of the second family (see [8] and (), () in Section 6) with a slight adjustment in the range of the indices. A further comparison of these three point sets in terms of numerical simulations is given in the last section of this article. Finally we would like to add that the Lisa points, similarly to the Padua points, can be considered as two-dimensional Chebyshev lattices of rank (see []). 3 Quadrature formulas based on the Lissajous node points In this section, we study quadrature rules for bivariate integration defined by point evaluations at the points Lisa n,p. As underlying polynomial spaces in R, we consider Π n = spant i (x)t j (y) : i + j n, where T i (x) denotes the Chebyshev polynomial T i (x) = cos(i arccos x) of the first kind. It is well-known (cf. [4]) that T i (x)t j (y) : i + j n is an orthogonal basis of Π n with respect to the inner product f, g := π f(x, y)g(x, y) dxdy. () x y The corresponding orthonormal basis is given by ˆT i (x) ˆT j (y) : i + j n, where, if i =, ˆT i (x) = Ti (x), if i. Using the trajectory γ n,p, it is possible to reduce a double integral of the form used in () into a single integral for a large class of bivariate polynomials. Lemma For all polynomials P Π 8n+4p with P, T (n+p) (x)t n (y) =, the following formula holds: π P (x, y) x dxdy = y π π P (γ n,p (t))dt. ()
7 Lagrange interpolation at Lissajous nodes 7 Proof We check () for all basis polynomials T i (x)t j (y) in the space Π 8n+4p. For the left hand side of () we get the value if (i, j) = (, ) and otherwise. For the right hand side of () we get also if (i, j) = (, ). For (i, j) (, ) we get for P (x, y) = T i (x)t j (y) the expression π π P (γ n,p (t))dt = π = π π π T i (sin(nt))t j (sin((n + p)t))dt ( cos int i π ) ( cos j(n + p)t j π ) dt. We now determine for which indices (i, j) this integral is different from. This is only the case if in = j(n+p) and i j is even. Since the numbers n and n+p are relatively prime, this can only be the case if i = k(n+p), j = nk and k N is an even number. We see that the smallest possible value for k is k = and the second smallest is k = 4. Furthermore, the sum of the respective indices is given by i+j = (n+p)k. Therefore, we can conclude that for all indices (i, j) satisfying i + j =,..., 4n + p and i + j = 4n + p +,..., 8n + 4p the right hand side of () vanishes. If i + j = 4n + p, the above integral is nonzero only if i = (n + p) and j = n. To get a quadrature formula supported on the points Lisa n,p, we define a suitable polynomial subspace with the index set Γ Q n,p N given by Γ Q n,p := Π Q n,p = spant i (x)t j (y) : (i, j) Γ Q n,p (i, j) : i+j 4n 4p m= (i, j) : i+j = 4n+m, j < n(4p m) p Note that the particular index (i, j) = ((n+p), n) is not included in Γ Q n,p and that Lemma is applicable for all polynomials P Π Q n,p. An example of the index set Γ Q n,p is shown in Figure. Clearly, the polynomial space Π Q n,p satisfies Π 4n Π Q n,p Π 4n+4p and the dimension of Π Q n,p can be computed as dim Π Q n,p = Γ Q n,p = 8n(n + p) + 4n (p ) = 4( Lisa n,p n p) (p ). For points A Lisa n,p, we define the quadrature weights w A :=,, if A Lisaout n,p, if A Lisaint n,p. Then, we get the following quadrature rule based on the node set Lisa n,p :.
8 8 W. Erb et al. j i Fig. Illustration of the index set Γ Q, with black and white bullets. We have Γ Q, = 56 black and white bullets. The black bullets correspond to indices describing the polynomial space Π 4n. The black cross is not contained in Γ Q,. It corresponds to the special index (i, j) = (6, 4) appearing in Lemma. Theorem For all P Π Q n,p the quadrature formula is exact. π P (x, y) x y dxdy = A Lisa n,p w A P (A) () Proof For all trigonometric π-periodic polynomials q of degree less than, the following composite trapezoidal quadrature rule is exact: π π q(t)dt = 4n(n+p) q (t k ). Since Π Q n,p Π 8n+4p and Π Q n,p T (n+p) (x)t n (y), we have by Lemma the identity π P (x, y) x k= dxdy = y π π P (γ n,p (t))dt. Thus, if we show that for P Π Q n,p the trigonometric polynomial P (γ n,p (t)) is of degree less than, we immediately get the quadrature formula π P (x, y) x 4n(n+p) dxdy = y k= = A Lisa n,p w A P (A). P (γ n,p (t k ))
9 Lagrange interpolation at Lissajous nodes 9 To finish the proof we consider the representation of the polynomial P Π Q n,p in the orthogonal basis T i (x)t j (y) : (i, j) Γ Q n,p and get P (γ n,p (t)) = (i,j) Γ Q n,p = (i,j) Γ Q n,p a ij T i (sin nt)t j (sin(n + p)t) a ij cos ( ) ( ) int i π cos j(n + p)t j π for some coefficients a ij R. In order for the trigonometric polynomials in this formula to have a degree less than, the indices (i, j) have to satisfy the condition (i + j)n + jp <. In the case that i + j < 4n, we have (i + j)n + jp (i + j)n + 4np < and the condition above is satisfied. In the case that i + j = 4n + m with m 4p, we have (i + j)n + jp = 4n + mn + jp < and the condition above is satisfied for all j with j < n(4p m) p. By definition, this condition is exactly satisfied for all indices (i, j) Γn,p Q and therefore for all polynomials P Πn,p. Q Remark Lemma and Theorem are generalizations of corresponding results proven in [3] for the Padua points. An analogous formula also exists for the Xu points (see [, 4]). Furthermore, the cardinality of the Xu points XU n is known to be minimal for exact integration of bivariate polynomials in Π 4n with respect to a product Chebyshev weight function (see [,4]). Since Lisa n,p > XU n = n(n + ), this is not the case for the Lisa points. On the other hand, as illustrated in Figure, the space Π Q n,p, for which () is exact, shows a remarkable asymmetry. As for multivariate interpolation, the construction of suitable nodes for cubature rules has a long history. For an overview, we refer to the survey article []. 4 Interpolation on the Lissajous node points Given the quadrature formulas of the last section, we now investigate bivariate interpolation at the points Lisa n,p. The corresponding interpolation problem can be formulated as follows: for given function values f A R, A Lisa n,p, we want to find a unique bivariate interpolating polynomial L n,p f such that L n,p f(a) = f A for all A Lisa n,p. (3) To set this problem correctly, we have to fix an underlying interpolation space. This space is linked to Π Q n,p and defined as Π L n,p := spant i (x)t j (y) : (i, j) Γ L n,p
10 W. Erb et al. on the index set Γ L n,p := (i, j) : i+j n p m= (i, j) : i + j = n + m, j < n(p m). p j 4 3 j i (a) Index set Γ, L with Γ, L = 7. (b) Index set Γ,3 L with Γ,3 L = 7. i Fig. 3 Illustration of the index sets Γ L n,p for n =, p = and n =, p = 3. The black bullets correspond to indices describing the polynomial space Π n. Examples of sets Γn,p L with different values of p are given in Figure 3. The reproducing kernel Kn,p L : R R R of the polynomial space Πn,p L is given as Kn,p(x L A, y A ; x B, y B ) = ˆT i (x A ) ˆT i (x B ) ˆT j (y A ) ˆT j (y B ). (i,j) Γ L n,p It is straightforward to check that the kernel K L n,p has the reproducing property P, K L n,p(x, y; ) = P (x, y) for all polynomials P Π L n,p. We have Π n Π L n,p Π (n+p). The dimension of the polynomial space Π L n,p is given as dim(π L n,p) = Γ L n,p = (n + )(n + ) + p m= n(p m) = n + n(p + ) + p = Lisa n,p. Therefore, the dimension dim(π L n,p) of the polynomial space Π L n,p corresponds precisely to the number of distinct points in Lisa n,p. Soon, we will deduce a formula for the fundamental polynomials of Lagrange interpolation with respect to the points in Lisa n,p and show that the p
11 Lagrange interpolation at Lissajous nodes interpolation problem (3) has a unique solution. To this end, we investigate an isomorphism between the polynomial space Π L n,p and the subspace Π trig,l n(n+p) := q Π trig n(n+p) : q(t k) = q(t k ) for all k, k with k Lisan,p k (4) of π-periodic trigonometric polynomials Π trig n(n+p) := n(n+p) q(t) = m= n(n+p) a m cos(mt) + m= b m sin(mt) : a m, b m R. Theorem The operator E γ : Π L n,p Π trig,l n(n+p), E γp (t) = P (γ n,p (t)), t [, π], defines an isometric isomorphism from the space ( Πn,p, L, ) onto the space Π trig,l n(n+p) equipped with the inner product q, q = π q (t)q (t)dt. π Proof The system ˆTi (x) ˆT j (y) : (i, j) Γn,p L forms an orthonormal basis of the space Π L n,p. The image e ij (t) := E γ (ˆTi (x) ˆT ) j (y) (t), of this basis under the linear operator E γ is given by (i, j) Γ L n,p,, if (i, j) = (, ), ( ) cos int i π, if j =, i < (n + p), e ij (t) = ( ) cos j(n + p)t j π, if i =, j n, cos ( ) ( ) int i π cos j(n + p)t j π, otherwise. (5) For (i, j) Γn,p, L j n, the functions e ij (t) are trigonometric polynomials of degree less than n(n + p). The only trigonometric polynomial of exact degree n(n + p) is precisely e,n. By the definition of the operator E γ, the values E γ P (t k ) and E γ P (t k ), t k t k coincide if γ n,p (t k ) = γ n,p (t k ) is a selfintersection point of γ n,p. This is precisely encoded in the constraints given in (4). We can conclude that E γ maps Πn,p L into the space Π trig,l n(n+p). For polyonomials P, P Πn,p, L the product polynomial P P is an element of the space Π 8n+4p and satisfies P P, T (n+p) (x)t n (y) =. Therefore, by Lemma, the set e ij : (i, j) Γn,p L is an orthonormal system in Π trig,l n(n+p), and thus, E γ is an isometric embedding from Πn,p L into Π trig,l n(n+p).
12 W. Erb et al. Now, if we can show that the dimensions of Πn,p L and Π trig,l n(n+p) coincide, the proof is finished. To this end, we consider in Π trig n(n+p) the Dirichlet kernel D n(n+p) (t) := n(n+p) + cos(n(n + p)t) + cos(kt) k= It is well known that the trigonometric polynomials = sin(n(n + p)t) cos t sin t. D k n(n+p) (t) := D n(n+p) (t t k ), k =,...,, are precisely the Lagrange polynomials in the space Π trig n(n+p) with respect to the points t k, k =,...,, i.e. D k n(n+p) (t k ) = δ k,k, k, k. In general, the polynomials Dn(n+p) k do not lie in the subspace Πtrig,L n(n+p). However, we can define a basis for Π trig,l n(n+p) by using the linear combinations l A (t) := Dn(n+p) k (t), A Lisa n,p. (6) k=,...,4n(n+p): k [A] Clearly, the polynomials l A are elements of Π trig,l n(n+p), and l A(t k ) is equal to one if k [A] and zero if k / [A]. Also, by the orthogonality of the functions Dn(n+p) k, we have l A, l B = if A B. Therefore, the system l A : A Lisa n,p forms an orthogonal basis of Π trig,l n(n+p) and dim(πtrig,l n(n+p) ) = Lisa n,p. This corresponds exactly with the dimension of the space Πn,p. L Theorem 3 For A = (x A, y A ) Lisa n,p, the polynomials L A := E γ l A have the representation L A (x, y) = w A (Kn,p(x, L y; x A, y A ) ) ˆT n (y) ˆT n (y A ) (7) and are the fundamental polynomials of Lagrange interpolation in the space Π L n,p on the point set Lisa n,p, i.e. L A (B) = δ A,B, A, B Lisa n,p. The interpolation problem (3) has a unique solution in Π L n,p and the interpolating polynomial L n,p f is given by L n,p f(x, y) = A Lisa n,p f A L A (x, y).
13 Lagrange interpolation at Lissajous nodes 3 Proof From the definition (6) of the trigonometric polynomials l A and the mapping E γ it follows immediately that the polynomials L A = E γ l A satisfy L A (B) = δ A,B for B Lisa n,p. Moreover, since the trigonometric polynomials l A : A Lisa n,p form an orthogonal basis of the space Π trig,l n(n+p), Theorem implies that the polynomials L A : A Lisa n,p form an orthogonal basis of Lagrange polynomials for the space Πn,p L as well. It remains to prove (7). To this end, we compute the decomposition of the polynomials l A in the basis e ij given in (5) and use the inverse of the operator E γ to obtain (7). The proof will be given only for A Lisa b n,p having the representation ( ) (r + )n + s(n + p) ( ) A = (x A, y A ) = γ n,p π = ( ) r+s z (n+p) (r+)p, zn sp. We first suppose that A Lisa b n,p is an interor point such that the two points k, k [A] [, ], k k that represent the same A are given as k = (r + )n + s(n + p) k = (r + )n s(n + p) mod, mod. Using simple trigonometric transformations, the basis function l A can be written as l A (t) = Dn(n+p) k (t) + Dk n(n+p) (t) = + cos((r + )nπ) cos(n(n+p)t) n(n+p) + m= cos ( smπ n ) ( ( ) cos (r+)mπ (n+p) cos(mt) + sin ( ) sin(mt)) (r+)mπ (n+p). Now, using the explicit expression (5) of the basis polynomials e ij and comparing the coefficients in the decomposition of l A, we get the following formula for the inner product l A, e ij = π π l A (t)e ij (t)dt, (i, j) Γn,p: L, if (i, j) = (, ), l A, e ij = = 4n(n+p) 4n(n+p) ( ) (r+s)j 4n(n+p) ( ) (r+s)i 4n(n+p) 4( ) (r+s)(i+j) 4n(n+p), if i =, j = n, cos ( j spπ ) n, if i =, j < n, ( ) cos i (r+)pπ (n+p), if i, j =, ( ) cos i (r+)pπ (n+p) cos ( j spπ ) n, otherwise, ˆT n (y A ), if i =, j = n, 4n(n+p) ˆT i (x A ) ˆT j (y A ), if (i, j) Γn,p L \ (, n).
14 4 W. Erb et al. Therefore, l A (t) can be decomposed as l A (t) = ˆT n (y A ) e,n(t) + (i,j) Γ L n,p j n ˆT i (x A ) ˆT j (y A ) e ij (t). Now, using the inverse mapping E γ together with the definition of the reproducing kernel Kn,p, L we can conclude: L A (x, y) = E γ l A (x, y) = w A (Kn,p(x, L y; x A, y A ) ) ˆT n (y) ˆT n (y A ). If A Lisa b n,p is a point on the boundary of the square [, ], the number k can be represented as k = (r + )n and the basis function l A is given as l A (t) = Dn(n+p) k (t) = + cos((r + )nπ) cos(n(n + p)t) n(n+p) ( ) ( ) r + r + + cos (n + p) mπ cos(mt) + sin (n + p) mπ sin(mt). m= Now, similar calculations to the above yield (7) with the half sized weight function w A = 4n(n+p). Finally, for all points A Lisaw n,p, (7) can be obtained by analogous calculations using the representation (4) instead of (3). Remark 3 (7) has a remarkable resemblence to the Lagrange polynomials of the Padua points. For the Padua points, the analog statement of Theorem 3 can be proved very elegantly by using ideal theory (cf. [5]). This approach was, however, not successful for the more general Lissajous nodes. Here, we had to use the isomorphism E γ and Theorem instead. 5 A simple scheme for the computation of the interpolation polynomial In view of Theorem 3, the solution to the interpolation problem (3) in Πn,p L is given as L n,p f(x, y) = ( w A f A Kn,p(x, L y; x A, y A ) ) ˆT n (y) ˆT n (y A ). A Lisa n,p The representation of the polynomial L n,p f(x, y) in the orthonormal Chebyshev basis ˆT i (x) ˆT j (y) : (i, j) Γ L n,p can now be written as L n,p f(x, y) = (i,j) Γ L n,p c i,j ˆTi (x) ˆT j (y)
15 Lagrange interpolation at Lissajous nodes 5 with the coefficients c i,j = L n,p f, ˆT i (x) ˆT j (y) given by w A f A ˆTi (x A ) ˆT j (y A ), if (i, j) Γn,p\(, L n), A Lisa n,p c i,j = w A f A ˆTn (y A ), if (i, j) = (, n). A Lisa n,p Using a matrix formulation, this identity can be written more compactly. We introduce the coefficient matrix C n,p = (c ij ) R (n+p) n+ by c ij = Next, we define the diagonal matrix Ln,p f, ˆT i (x) ˆT j (y), if (i, j) Γ L n,p,, otherwise. D f (Lisa n,p ) = diag (w A f A, A Lisa n,p ) R Lisan,p Lisan,p. Further, for a general finite set X R of points, we introduce the matrices ˆT (x B ) T x (X ) =. R (n+p) X, ˆT (n+p) (x B ) T y (X ) = B X ˆT (y B ). R n+ X. ˆT n (y B ) B X Finally, we define the mask M n,p = (m ij ) R (n+p) n+ by, if (i, j) Γn,p L \ (, n), m i,j = /, if (i, j) = (, n),, if (i, j) / Γn,p. L Now, the coefficient matrix C n,p of the interpolating polynomial can be computed as C n,p = ( T x (Lisa n,p )D f (Lisa n,p )T y (Lisa n,p ) T ) M n,p, (8) where denotes pointwise multiplication of the matrix entries. For an arbitrary point B R, the evaluation L n,p f(b) of the interpolation polynomial L n,p f at B is then given by L n,p f(b) = T x (B) T C n,p T y (B). (9)
16 6 W. Erb et al. Remark 4 The matrix formulation in (8) and (9) is almost identical to the formulation of the interpolating scheme of the Padua points given in [8]. This is due to the similarity in the representation (7) of the Lagrange polynomials. The main difference between the schemes lies in the form of the mask M n,p. The mask M n,p for the Lisa points has an asymmetric structure determined by the index set Γ L n,p, whereas the matrix is an upper left triangular matrix for the Padua points. Two examples of such a structure are given in Figure 3. 6 Numerical Simulations Based on the results derived in the last sections, we perform numerical simulations on the behaviour of the Lisa points (LS) in comparison to some already established point sets. Unless explicitly mentioned, we assume p = for all numerical simulations of the Lisa points. For the comparison point sets, our focus is on the Xu points (XU) [4] and the Padua points (PD) [7]. Based on the Chebychev-Lobatto points given by (5), and in correspondance to (8) and (9), the odd Xu points XU n+ are defined as the union of the sets XU b n+ = (z n+ i, z n+ j ) : i n, j n, XU w n+ = (z n+ i+, zn+ j+ ) : i n, j n, with the cardinality XU n+ = (n + ). In turn, the even Padua points PD n (nd family) are defined as the union of the sets PD b n = (z n+ i+, zn j ) : i n, j n, () PD w n = (z n+ i, zj+) n : i n, j n. () The cardinality can be calculated as PD n = (n+)(n+). The distributions of the Lisa, Xu and Padua points are shown for small degrees of n in Figure 4. The point sets are introduced in such a way that an equally chosen n results in a similar cardinality. The stability of the mapping f L n,p f is evaluated by means of the growth of the Lebesgue constant. Here, we calculate the values of the Lebesgue constant Λ LS n,p = max B [,] A Lisa n,p L A (B) of the Lisa points up to a degree of n = 6. We compare them with the leastsquares fitting of the Lebesgue constant for the Padua and the Xu points. As shown in [7], it holds for the Padua points that Λ PD n = ( π log(n + ) +.) and as presented for the Xu points in [] that Λ XU n+ = ( π log(n + )). Figure 5(a) indicates that the asymptotic growth of Λ LS n, corresponds to the order O ( log n ) of the Lebesgue constant Λ PD n. In Figure 5(b) it is shown, how a variation of the parameter p of the Lisa points changes the growth of the Lebesgue constant. Here, we consider p =, 3, 5, 7 and excluded each entry for n and p not being relatively prime. In total, these numerical evaluations
17 Lagrange interpolation at Lissajous nodes 7.5 LS XU PD.5 LS XU PD (a) Point sets for n =..5.5 (b) Point sets for n = 5. Fig. 4 Visualizations of the Lisa (LS), Xu (XU) and Padua (PD) point sets. suggest the conjecture that the Lebesgue constant of the Lisa points is of the same order O ( log n ) as the Lebesgue constant of the Padua and Xu points (see [4]). Lebesgue constant 5 5 LS XU PD Lebesgue constant 5 p = p = 3 p = 5 p = degree n (a) Lisa n,, XU n+ and PD n point sets degree n (b) Lisa point sets for p, 3, 5, 7. Fig. 5 Lebesgue constants up to a degree of n = 6 for the Lisa points in comparison to the least-squares fitting of the Lebesgue constant of the Xu and Padua points. For a further evaluation of the Lisa points, we perform numerical interpolations with the Xu, Padua and Lisa points on the Franke-Renka-Brown test set [,3]. In order to simulate the Xu points as well as the Padua points, the numerical algorithms presented in [6, 9] are used. The maximum interpolation errors are computed on a uniform grid of points defined in a region Ω = [, ] [, ]. As mentioned above, the degree n is defined to result in a similar total number of points, i.e. a similar cardinality. For our simulations we take n 5,,, 3. The results are shown in Table 4. It can be seen that the maximum interpolation error of all three point sets shows a similar
18 8 W. Erb et al. behaviour in terms of degree n, with respect to the chosen test function. In terms of the Lisa point sets, we evaluated the behaviour of p in addition to the aforementioned comparisons. We can state that the influence of varying p, with respect to the maximum interpolation error and the nodes used for the evaluation, is almost negligible. Table Interpolation errors for the points Lisa n,. n # F F F 3 F 4 F 5 F 6 F 7 F 8 F 9 F E- 4 E- E-3 6 E-5 E- 3 E-5 8 E- E- E+ 4 E- 4 7 E-3 7 E-3 E-6 E- E-5 E-8 E-5 4 E-3 4 E- 9 E- 88 E-6 E-4 4 E- 5 E-5 E-3 E-4 5 E-4 E-7 5 E-6 4 E E- 7 E-6 3 E-4 E-4 4 E-5 3 E-4 E-3 E-3 9 E- 3 E- Table 3 Interpolation errors for the points XU n+. n # F F F 3 F 4 F 5 F 6 F 7 F 8 F 9 F E- 3 E- E-3 6 E-5 E- 3 E-4 6 E- 3 E- 3 E+ 6 E- 4 5 E-3 6 E-3 E-6 E- E-5 E-8 E-5 5 E-3 4 E- E- 88 E-6 E-4 5 E- 3 E-5 E-3 5 E-5 3 E-4 E-7 5 E-6 4 E E- 7 E-6 E-4 5 E-5 3 E-5 9 E-5 4 E-4 5 E-4 9 E- E- Table 4 Interpolation errors for the Padua points PD n. n # F F F 3 F 4 F 5 F 6 F 7 F 8 F 9 F E- 4 E- E-3 6 E-5 E- 3 E-5 9 E- E- 4 E+ 5 E- 3 6 E-3 7 E-3 3 E-6 E- E-5 E-8 E-5 6 E-3 7 E- E- 86 E-6 E-4 7 E- 3 E-5 E-3 4 E-5 E-4 E-7 7 E-6 4 E E- 7 E-6 E-4 6 E-5 4 E-5 E-4 5 E-4 6 E-4 E- E- Acknowledgements The authors gratefully acknowledge the financial support of the German Federal Ministry of Education and Research (BMBF, grant number 3N9), the German Research Foundation (DFG, grant number BU 436/9- and ER 777/-), the European Union and the State Schleswig-Holstein (EFRE, grant number --4).
19 Lagrange interpolation at Lissajous nodes 9 References. Bogle, M.G.V., Hearst, J.E., Jones, V.F.R., Stoilov, L.: Lissajous knots. J. Knot Theory Ramifications 3(), 4 (994). Bos, L., Caliari, M., De Marchi, S., Vianello, M.: A Numerical Study of the Xu Polynomial Interpolation Formula in Two Variables. Computing 76(3), 3 34 (6) 3. Bos, L., Caliari, M., De Marchi, S., Vianello, M., Xu, Y.: Bivariate Lagrange interpolation at the Padua points: the generating curve approach. J. Approx. Theory 43(), 5 5 (6) 4. Bos, L., De Marchi, S., Vianello, M.: On the Lebesgue constant for the Xu interpolation formula. J. Approx. Theory 4(), 34 4 (6) 5. Bos, L., De Marchi, S., Vianello, M., Xu, Y.: Bivariate Lagrange interpolation at the Padua points: The ideal theory approach. Numer. Math. 8(), (7) 6. Caliari, M., De Marchi, S., Sommariva, A., Vianello, M.: PaduaDM: fast interpolation and cubature at the Padua points in Matlab/Octave. Numer. Algorithms 56(), 45 6 () 7. Caliari, M., De Marchi, S., Vianello, M.: Bivariate polynomial interpolation on the square at new nodal sets. Appl. Math. Comput. 65(), 6 74 (5) 8. Caliari, M., De Marchi, S., Vianello, M.: Bivariate Lagrange interpolation at the Padua points: Computational aspects. J. Comput. Appl. Math. (), 84 9 (8) 9. Caliari, M., Vianello, M., De Marchi, S., Montagna, R.: Hyperd: a numerical code for hyperinterpolation on rectangles. Appl. Math. Comput. 83(), (6). Cools, R.: Constructing cubature formulae: the science behind the art. Acta Numerica 6, 54 (997). Cools, R., Poppe, K.: Chebyshev lattices, a unifying framework for cubature with Chebyshev weight function. BIT 5(), (). Franke, R.: Scattered data interpolation: Tests of some methods. Math. Comp. 38(57), 8 (98) 3. Gasca, M., Sauer, T.: On the history of multivariate polynomial interpolation. J. Comput. Appl. Math. (-), 3 35 () 4. Gasca, M., Sauer, T.: Polynomial interpolation in several variables. Adv. Comput. Math. (4), () 5. Gleich, B., Weizenecker, J.: Tomographic imaging using the nonlinear response of magnetic particles. Nature 435(746), 4 7 (5) 6. Grüttner, M., Knopp, T., Franke, J., Heidenreich, M., Rahmer, J., Halkola, A., Kaethner, C., Borgert, J., Buzug, T.M.: On the formulation of the image reconstruction problem in magnetic particle imaging. Biomed. Tech. / Biomed. Eng. 58(6), (3) 7. Harris, L.A.: Bivariate Lagrange interpolation at the Geronimus nodes. Contemp. Math. 59, (3) 8. Kaethner, C., Ahlborg, M., Bringout, G., Weber, M., Buzug, T.M.: Axially elongated field-free point data acquisition in magnetic particle imaging. IEEE Trans. Med. Imag. (4). Accepted for publication 9. Knopp, T., Biederer, S., Sattel, T.F., Weizenecker, J., Gleich, B., Borgert, J., Buzug, T.M.: Trajectory analysis for magnetic particle imaging. Phys. Med. Biol. 54(), (9). Lamm, C.: There are infinitely many Lissajous knots. Manuscr. Math. 93(), 9 37 (997). Möller, H.M.: Kubaturformeln mit minimaler Knotenzahl. Numer. Math. 5, 85 (976). Morrow, C.R., Patterson, T.N.L.: Construction of algebraic cubature rules using polynomial ideal theory. SIAM J. Numer. Anal. 5, (978) 3. Renka, R.J., Brown, R.: Algorithm 79: Accuracy tests of acm algorithms for interpolation of scattered data in the plane. ACM Trans. Math. Softw. 5(), (999) 4. Xu, Y.: Lagrange interpolation on Chebyshev points of two variables. J. Approx. Theory 87(), 38 (996)
Trivariate polynomial approximation on Lissajous curves 1
Trivariate polynomial approximation on Lissajous curves 1 Stefano De Marchi DMV2015, Minisymposium on Mathematical Methods for MPI 22 September 2015 1 Joint work with Len Bos (Verona) and Marco Vianello
More informationA numerical study of the Xu polynomial interpolation formula in two variables
A numerical study of the Xu polynomial interpolation formula in two variables Len Bos Dept. of Mathematics and Statistics, University of Calgary (Canada) Marco Caliari, Stefano De Marchi Dept. of Computer
More informationBivariate Lagrange interpolation at the Padua points: The generating curve approach
Journal of Approximation Theory 43 (6) 5 5 www.elsevier.com/locate/jat Bivariate Lagrange interpolation at the Padua points: The generating curve approach Len Bos a, Marco Caliari b, Stefano De Marchi
More informationBIVARIATE LAGRANGE INTERPOLATION AT THE PADUA POINTS: THE IDEAL THEORY APPROACH
BIVARIATE LAGRANGE INTERPOLATION AT THE PADUA POINTS: THE IDEAL THEORY APPROACH LEN BOS, STEFANO DE MARCHI, MARCO VIANELLO, AND YUAN XU Abstract The Padua points are a family of points on the square [,
More informationStability of Kernel Based Interpolation
Stability of Kernel Based Interpolation Stefano De Marchi Department of Computer Science, University of Verona (Italy) Robert Schaback Institut für Numerische und Angewandte Mathematik, University of Göttingen
More informationInterpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials
Interpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials Lawrence A. Harris Abstract. We extend the definition of Geronimus nodes to include pairs of real numbers where
More informationOn the Lebesgue constant of subperiodic trigonometric interpolation
On the Lebesgue constant of subperiodic trigonometric interpolation Gaspare Da Fies and Marco Vianello November 4, 202 Abstract We solve a recent conjecture, proving that the Lebesgue constant of Chebyshev-like
More informationPadua2DM: fast interpolation and cubature at the Padua points in Matlab/Octave
Numerical Algorithms manuscript No. (will be inserted by the editor) Padua2DM: fast interpolation and cubature at the Padua points in Matlab/Octave Marco Caliari Stefano De Marchi Alvise Sommariva Marco
More informationPolynomial approximation and cubature at approximate Fekete and Leja points of the cylinder
Polynomial approximation and cubature at approximate Fekete and Leja points of the cylinder Stefano De Marchi, Martina Marchioro and Alvise Sommariva May 18, 2018 arxiv:1110.5513v1 [math.na] 25 Oct 2011
More informationDubiner distance and stability of Lebesgue constants
Dubiner distance and stability of Lebesgue constants Marco Vianello 1 March 15, 2018 Abstract Some elementary inequalities, based on the non-elementary notion of Dubiner distance, show that norming sets
More informationOn Multivariate Newton Interpolation at Discrete Leja Points
On Multivariate Newton Interpolation at Discrete Leja Points L. Bos 1, S. De Marchi 2, A. Sommariva 2, M. Vianello 2 September 25, 2011 Abstract The basic LU factorization with row pivoting, applied to
More informationON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT
ON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT Received: 31 July, 2008 Accepted: 06 February, 2009 Communicated by: SIMON J SMITH Department of Mathematics and Statistics La Trobe University,
More informationConstructing optimal polynomial meshes on planar starlike domains
Constructing optimal polynomial meshes on planar starlike domains F. Piazzon and M. Vianello 1 Dept. of Mathematics, University of Padova (Italy) October 13, 2014 Abstract We construct polynomial norming
More informationDepartment of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2014
Department of Mathematics, University of California, Berkeley YOUR 1 OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2014 1. Please write your 1- or 2-digit exam number on
More informationPolynomial Meshes: Computation and Approximation
Proceedings of the 15th International Conference on Computational and Mathematical Methods in Science and Engineering, CMMSE 215 6 1 July, 215. Polynomial Meshes: Computation and Approximation S. De Marchi
More informationOptimal Polynomial Admissible Meshes on the Closure of C 1,1 Bounded Domains
Optimal Polynomial Admissible Meshes on the Closure of C 1,1 Bounded Domains Constructive Theory of Functions Sozopol, June 9-15, 2013 F. Piazzon, joint work with M. Vianello Department of Mathematics.
More informationQuadrature Formula for Computed Tomography
Quadrature Formula for Computed Tomography orislav ojanov, Guergana Petrova August 13, 009 Abstract We give a bivariate analog of the Micchelli-Rivlin quadrature for computing the integral of a function
More informationOrthogonal Polynomials, Quadratures & Sparse-Grid Methods for Probability Integrals
1/33 Orthogonal Polynomials, Quadratures & Sparse-Grid Methods for Probability Integrals Dr. Abebe Geletu May, 2010 Technische Universität Ilmenau, Institut für Automatisierungs- und Systemtechnik Fachgebiet
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 35, pp. 148-163, 2009. Copyright 2009,. ISSN 1068-9613. ETNA SPHERICAL QUADRATURE FORMULAS WITH EQUALLY SPACED NODES ON LATITUDINAL CIRCLES DANIELA
More informationDISCRETE SUBGROUPS, LATTICES, AND UNITS.
DISCRETE SUBGROUPS, LATTICES, AND UNITS. IAN KIMING 1. Discrete subgroups of real vector spaces and lattices. Definitions: A lattice in a real vector space V of dimension d is a subgroup of form: Zv 1
More informationReconstruction of sparse Legendre and Gegenbauer expansions
Reconstruction of sparse Legendre and Gegenbauer expansions Daniel Potts Manfred Tasche We present a new deterministic algorithm for the reconstruction of sparse Legendre expansions from a small number
More informationStrongly Regular Decompositions of the Complete Graph
Journal of Algebraic Combinatorics, 17, 181 201, 2003 c 2003 Kluwer Academic Publishers. Manufactured in The Netherlands. Strongly Regular Decompositions of the Complete Graph EDWIN R. VAN DAM Edwin.vanDam@uvt.nl
More informationMath 215B: Solutions 1
Math 15B: Solutions 1 Due Thursday, January 18, 018 (1) Let π : X X be a covering space. Let Φ be a smooth structure on X. Prove that there is a smooth structure Φ on X so that π : ( X, Φ) (X, Φ) is an
More informationSYMMETRIC SUBGROUP ACTIONS ON ISOTROPIC GRASSMANNIANS
1 SYMMETRIC SUBGROUP ACTIONS ON ISOTROPIC GRASSMANNIANS HUAJUN HUANG AND HONGYU HE Abstract. Let G be the group preserving a nondegenerate sesquilinear form B on a vector space V, and H a symmetric subgroup
More informationA orthonormal basis for Radial Basis Function approximation
A orthonormal basis for Radial Basis Function approximation 9th ISAAC Congress Krakow, August 5-9, 2013 Gabriele Santin, joint work with S. De Marchi Department of Mathematics. Doctoral School in Mathematical
More informationApplication of modified Leja sequences to polynomial interpolation
Special Issue for the years of the Padua points, Volume 8 5 Pages 66 74 Application of modified Leja sequences to polynomial interpolation L. P. Bos a M. Caliari a Abstract We discuss several applications
More informationSparse polynomial interpolation in Chebyshev bases
Sparse polynomial interpolation in Chebyshev bases Daniel Potts Manfred Tasche We study the problem of reconstructing a sparse polynomial in a basis of Chebyshev polynomials (Chebyshev basis in short)
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationarxiv: v2 [math.na] 27 Dec 2016
An algorithm for constructing Equiangular vectors Azim rivaz a,, Danial Sadeghi a a Department of Mathematics, Shahid Bahonar University of Kerman, Kerman 76169-14111, IRAN arxiv:1412.7552v2 [math.na]
More informationMATHEMATICAL METHODS AND APPLIED COMPUTING
Numerical Approximation to Multivariate Functions Using Fluctuationlessness Theorem with a Trigonometric Basis Function to Deal with Highly Oscillatory Functions N.A. BAYKARA Marmara University Department
More informationMagnetic Nanoparticle Imaging for Biomedical Applications
Magnetic Nanoparticle Imaging for Biomedical Applications by Nurrniza Binti Othman DOCTORAL THESISUBMITTED TO THE GRADUATE SCHOOL OF INFORMATION SCIENCE AND ELECTRICAL ENGINEERING PARTIAL FULFILLMENT OF
More informationOn Riesz-Fischer sequences and lower frame bounds
On Riesz-Fischer sequences and lower frame bounds P. Casazza, O. Christensen, S. Li, A. Lindner Abstract We investigate the consequences of the lower frame condition and the lower Riesz basis condition
More informationSolutions to Problems for Math. H90 Issued 27 Aug. 2007
Problem 1: Given is an ellipse E neither a circle nor degenerate (i.e. a straight line segment). Let be the largest of the areas of triangles inscribed in E. How many inscribed triangles have maximal area?
More informationDiscrete Orthogonal Polynomials on Equidistant Nodes
International Mathematical Forum, 2, 2007, no. 21, 1007-1020 Discrete Orthogonal Polynomials on Equidistant Nodes Alfredo Eisinberg and Giuseppe Fedele Dip. di Elettronica, Informatica e Sistemistica Università
More informationVandermonde Matrices for Intersection Points of Curves
This preprint is in final form. It appeared in: Jaén J. Approx. (009) 1, 67 81 Vandermonde Matrices for Intersection Points of Curves Hakop Hakopian, Kurt Jetter and Georg Zimmermann Abstract We reconsider
More informationCodewords of small weight in the (dual) code of points and k-spaces of P G(n, q)
Codewords of small weight in the (dual) code of points and k-spaces of P G(n, q) M. Lavrauw L. Storme G. Van de Voorde October 4, 2007 Abstract In this paper, we study the p-ary linear code C k (n, q),
More informationMath 120 HW 9 Solutions
Math 120 HW 9 Solutions June 8, 2018 Question 1 Write down a ring homomorphism (no proof required) f from R = Z[ 11] = {a + b 11 a, b Z} to S = Z/35Z. The main difficulty is to find an element x Z/35Z
More informationProblems in Linear Algebra and Representation Theory
Problems in Linear Algebra and Representation Theory (Most of these were provided by Victor Ginzburg) The problems appearing below have varying level of difficulty. They are not listed in any specific
More informationStability constants for kernel-based interpolation processes
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report 59 Stability constants for kernel-based interpolation processes Stefano De Marchi Robert Schaback Dipartimento
More informationMath 396. An application of Gram-Schmidt to prove connectedness
Math 396. An application of Gram-Schmidt to prove connectedness 1. Motivation and background Let V be an n-dimensional vector space over R, and define GL(V ) to be the set of invertible linear maps V V
More informationChapter 2: Linear Independence and Bases
MATH20300: Linear Algebra 2 (2016 Chapter 2: Linear Independence and Bases 1 Linear Combinations and Spans Example 11 Consider the vector v (1, 1 R 2 What is the smallest subspace of (the real vector space
More informationOn the Lebesgue constant of barycentric rational interpolation at equidistant nodes
On the Lebesgue constant of barycentric rational interpolation at equidistant nodes by Len Bos, Stefano De Marchi, Kai Hormann and Georges Klein Report No. 0- May 0 Université de Fribourg (Suisse Département
More informationPriority Program 1253
Deutsche Forschungsgemeinschaft Priority Program 1253 Optimization with Partial Differential Equations Klaus Deckelnick and Michael Hinze A note on the approximation of elliptic control problems with bang-bang
More informationLinear Algebra and Robot Modeling
Linear Algebra and Robot Modeling Nathan Ratliff Abstract Linear algebra is fundamental to robot modeling, control, and optimization. This document reviews some of the basic kinematic equations and uses
More informationALGEBRAIC GEOMETRY I - FINAL PROJECT
ALGEBRAIC GEOMETRY I - FINAL PROJECT ADAM KAYE Abstract This paper begins with a description of the Schubert varieties of a Grassmannian variety Gr(k, n) over C Following the technique of Ryan [3] for
More informationLinear Algebra Notes. Lecture Notes, University of Toronto, Fall 2016
Linear Algebra Notes Lecture Notes, University of Toronto, Fall 2016 (Ctd ) 11 Isomorphisms 1 Linear maps Definition 11 An invertible linear map T : V W is called a linear isomorphism from V to W Etymology:
More informationScattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions
Chapter 3 Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions 3.1 Scattered Data Interpolation with Polynomial Precision Sometimes the assumption on the
More informationMATH JORDAN FORM
MATH 53 JORDAN FORM Let A,, A k be square matrices of size n,, n k, respectively with entries in a field F We define the matrix A A k of size n = n + + n k as the block matrix A 0 0 0 0 A 0 0 0 0 A k It
More informationLinear Algebra using Dirac Notation: Pt. 2
Linear Algebra using Dirac Notation: Pt. 2 PHYS 476Q - Southern Illinois University February 6, 2018 PHYS 476Q - Southern Illinois University Linear Algebra using Dirac Notation: Pt. 2 February 6, 2018
More informationKatholieke Universiteit Leuven Department of Computer Science
Extensions of Fibonacci lattice rules Ronald Cools Dirk Nuyens Report TW 545, August 2009 Katholieke Universiteit Leuven Department of Computer Science Celestijnenlaan 200A B-3001 Heverlee (Belgium Extensions
More informationMATH 225 Summer 2005 Linear Algebra II Solutions to Assignment 1 Due: Wednesday July 13, 2005
MATH 225 Summer 25 Linear Algebra II Solutions to Assignment 1 Due: Wednesday July 13, 25 Department of Mathematical and Statistical Sciences University of Alberta Question 1. [p 224. #2] The set of all
More informationImproved Lebesgue constants on the triangle
Journal of Computational Physics 207 (2005) 625 638 www.elsevier.com/locate/jcp Improved Lebesgue constants on the triangle Wilhelm Heinrichs Universität Duisburg-Essen, Ingenieurmathematik (FB 10), Universitätsstrasse
More informationAlgebra Workshops 10 and 11
Algebra Workshops 1 and 11 Suggestion: For Workshop 1 please do questions 2,3 and 14. For the other questions, it s best to wait till the material is covered in lectures. Bilinear and Quadratic Forms on
More informationLecture Notes Introduction to Cluster Algebra
Lecture Notes Introduction to Cluster Algebra Ivan C.H. Ip Update: May 16, 2017 5 Review of Root Systems In this section, let us have a brief introduction to root system and finite Lie type classification
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 2 Approximation Theory. Antony Jameson
Advanced Computational Fluid Dynamics AA5A Lecture Approximation Theory Antony Jameson Winter Quarter, 6, Stanford, CA Last revised on January 7, 6 Contents Approximation Theory. Least Squares Approximation
More informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More informationBiorthogonal Spline Type Wavelets
PERGAMON Computers and Mathematics with Applications 0 (00 1 0 www.elsevier.com/locate/camwa Biorthogonal Spline Type Wavelets Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan
More informationOn surface-knots with triple point number at most three
On surface-knots with triple point number at most three A. Al Kharusi and T. Yashiro Abstract arxiv:1711.04838v1 [math.at] 13 Nov 2017 In this paper, we show that if a diagram of a surface-knot F has at
More informationwhich arises when we compute the orthogonal projection of a vector y in a subspace with an orthogonal basis. Hence assume that P y = A ij = x j, x i
MODULE 6 Topics: Gram-Schmidt orthogonalization process We begin by observing that if the vectors {x j } N are mutually orthogonal in an inner product space V then they are necessarily linearly independent.
More informationUniqueness of the Solutions of Some Completion Problems
Uniqueness of the Solutions of Some Completion Problems Chi-Kwong Li and Tom Milligan Abstract We determine the conditions for uniqueness of the solutions of several completion problems including the positive
More informationOn the usage of lines in GC n sets
On the usage of lines in GC n sets Hakop Hakopian, Vahagn Vardanyan arxiv:1807.08182v3 [math.co] 16 Aug 2018 Abstract A planar node set X, with X = ( ) n+2 2 is called GCn set if each node possesses fundamental
More informationarxiv: v1 [math.gr] 8 Nov 2008
SUBSPACES OF 7 7 SKEW-SYMMETRIC MATRICES RELATED TO THE GROUP G 2 arxiv:0811.1298v1 [math.gr] 8 Nov 2008 ROD GOW Abstract. Let K be a field of characteristic different from 2 and let C be an octonion algebra
More informationimplies that if we fix a basis v of V and let M and M be the associated invertible symmetric matrices computing, and, then M = (L L)M and the
Math 395. Geometric approach to signature For the amusement of the reader who knows a tiny bit about groups (enough to know the meaning of a transitive group action on a set), we now provide an alternative
More informationNumerical integration formulas of degree two
Applied Numerical Mathematics 58 (2008) 1515 1520 www.elsevier.com/locate/apnum Numerical integration formulas of degree two ongbin Xiu epartment of Mathematics, Purdue University, West Lafayette, IN 47907,
More informationEGYPTIAN FRACTIONS WITH EACH DENOMINATOR HAVING THREE DISTINCT PRIME DIVISORS
#A5 INTEGERS 5 (205) EGYPTIAN FRACTIONS WITH EACH DENOMINATOR HAVING THREE DISTINCT PRIME DIVISORS Steve Butler Department of Mathematics, Iowa State University, Ames, Iowa butler@iastate.edu Paul Erdős
More informationOn the efficiency of the Peaceman-Rachford ADI-dG method for wave-type methods
On the efficiency of the Peaceman-Rachford ADI-dG method for wave-type methods Marlis Hochbruck, Jonas Köhler CRC Preprint 2017/34 (revised), March 2018 KARLSRUHE INSTITUTE OF TECHNOLOGY KIT The Research
More informationORTHONORMAL SAMPLING FUNCTIONS
ORTHONORMAL SAMPLING FUNCTIONS N. KAIBLINGER AND W. R. MADYCH Abstract. We investigate functions φ(x) whose translates {φ(x k)}, where k runs through the integer lattice Z, provide a system of orthonormal
More informationA new stable basis for radial basis function interpolation
A new stable basis for radial basis function interpolation Stefano De Marchi and Gabriele Santin Department of Mathematics University of Padua (Italy) Abstract It is well-known that radial basis function
More informationEXERCISES IN MODULAR FORMS I (MATH 726) (2) Prove that a lattice L is integral if and only if its Gram matrix has integer coefficients.
EXERCISES IN MODULAR FORMS I (MATH 726) EYAL GOREN, MCGILL UNIVERSITY, FALL 2007 (1) We define a (full) lattice L in R n to be a discrete subgroup of R n that contains a basis for R n. Prove that L is
More informationMATH 223A NOTES 2011 LIE ALGEBRAS 35
MATH 3A NOTES 011 LIE ALGEBRAS 35 9. Abstract root systems We now attempt to reconstruct the Lie algebra based only on the information given by the set of roots Φ which is embedded in Euclidean space E.
More informationWe simply compute: for v = x i e i, bilinearity of B implies that Q B (v) = B(v, v) is given by xi x j B(e i, e j ) =
Math 395. Quadratic spaces over R 1. Algebraic preliminaries Let V be a vector space over a field F. Recall that a quadratic form on V is a map Q : V F such that Q(cv) = c 2 Q(v) for all v V and c F, and
More informationNumerical Methods I: Polynomial Interpolation
1/31 Numerical Methods I: Polynomial Interpolation Georg Stadler Courant Institute, NYU stadler@cims.nyu.edu November 16, 2017 lassical polynomial interpolation Given f i := f(t i ), i =0,...,n, we would
More informationMATH 8253 ALGEBRAIC GEOMETRY WEEK 12
MATH 8253 ALGEBRAIC GEOMETRY WEEK 2 CİHAN BAHRAN 3.2.. Let Y be a Noetherian scheme. Show that any Y -scheme X of finite type is Noetherian. Moreover, if Y is of finite dimension, then so is X. Write f
More informationQUATERNIONS AND ROTATIONS
QUATERNIONS AND ROTATIONS SVANTE JANSON 1. Introduction The purpose of this note is to show some well-known relations between quaternions and the Lie groups SO(3) and SO(4) (rotations in R 3 and R 4 )
More informationRational Bézier Patch Differentiation using the Rational Forward Difference Operator
Rational Bézier Patch Differentiation using the Rational Forward Difference Operator Xianming Chen, Richard F. Riesenfeld, Elaine Cohen School of Computing, University of Utah Abstract This paper introduces
More informationQuantum wires, orthogonal polynomials and Diophantine approximation
Quantum wires, orthogonal polynomials and Diophantine approximation Introduction Quantum Mechanics (QM) is a linear theory Main idea behind Quantum Information (QI): use the superposition principle of
More informationEXAMPLES OF REFINABLE COMPONENTWISE POLYNOMIALS
EXAMPLES OF REFINABLE COMPONENTWISE POLYNOMIALS NING BI, BIN HAN, AND ZUOWEI SHEN Abstract. This short note presents four examples of compactly supported symmetric refinable componentwise polynomial functions:
More informationarxiv: v1 [math.co] 3 Nov 2014
SPARSE MATRICES DESCRIBING ITERATIONS OF INTEGER-VALUED FUNCTIONS BERND C. KELLNER arxiv:1411.0590v1 [math.co] 3 Nov 014 Abstract. We consider iterations of integer-valued functions φ, which have no fixed
More informationChapter 8 Integral Operators
Chapter 8 Integral Operators In our development of metrics, norms, inner products, and operator theory in Chapters 1 7 we only tangentially considered topics that involved the use of Lebesgue measure,
More informationSpectra of Multiplication Operators as a Numerical Tool. B. Vioreanu and V. Rokhlin Technical Report YALEU/DCS/TR-1443 March 3, 2011
We introduce a numerical procedure for the construction of interpolation and quadrature formulae on bounded convex regions in the plane. The construction is based on the behavior of spectra of certain
More informationarxiv: v2 [math.fa] 27 Sep 2016
Decomposition of Integral Self-Affine Multi-Tiles Xiaoye Fu and Jean-Pierre Gabardo arxiv:43.335v2 [math.fa] 27 Sep 26 Abstract. In this paper we propose a method to decompose an integral self-affine Z
More information2.3. VECTOR SPACES 25
2.3. VECTOR SPACES 25 2.3 Vector Spaces MATH 294 FALL 982 PRELIM # 3a 2.3. Let C[, ] denote the space of continuous functions defined on the interval [,] (i.e. f(x) is a member of C[, ] if f(x) is continuous
More informationOn interpolation by radial polynomials C. de Boor Happy 60th and beyond, Charlie!
On interpolation by radial polynomials C. de Boor Happy 60th and beyond, Charlie! Abstract A lemma of Micchelli s, concerning radial polynomials and weighted sums of point evaluations, is shown to hold
More informationB 1 = {B(x, r) x = (x 1, x 2 ) H, 0 < r < x 2 }. (a) Show that B = B 1 B 2 is a basis for a topology on X.
Math 6342/7350: Topology and Geometry Sample Preliminary Exam Questions 1. For each of the following topological spaces X i, determine whether X i and X i X i are homeomorphic. (a) X 1 = [0, 1] (b) X 2
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 9: Conditionally Positive Definite Radial Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH
More informationOn the representability of the bi-uniform matroid
On the representability of the bi-uniform matroid Simeon Ball, Carles Padró, Zsuzsa Weiner and Chaoping Xing August 1, 2012 Abstract Every bi-uniform matroid is representable over all sufficiently large
More informationA trigonometric orthogonality with respect to a nonnegative Borel measure
Filomat 6:4 01), 689 696 DOI 10.98/FIL104689M Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat A trigonometric orthogonality with
More informationMath 61CM - Solutions to homework 6
Math 61CM - Solutions to homework 6 Cédric De Groote November 5 th, 2018 Problem 1: (i) Give an example of a metric space X such that not all Cauchy sequences in X are convergent. (ii) Let X be a metric
More informationNodal O(h 4 )-superconvergence of piecewise trilinear FE approximations
Preprint, Institute of Mathematics, AS CR, Prague. 2007-12-12 INSTITTE of MATHEMATICS Academy of Sciences Czech Republic Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations Antti Hannukainen
More informationPositive Denite Matrix. Ya Yan Lu 1. Department of Mathematics. City University of Hong Kong. Kowloon, Hong Kong. Abstract
Computing the Logarithm of a Symmetric Positive Denite Matrix Ya Yan Lu Department of Mathematics City University of Hong Kong Kowloon, Hong Kong Abstract A numerical method for computing the logarithm
More informationTHE CLASSIFICATION OF PLANAR MONOMIALS OVER FIELDS OF PRIME SQUARE ORDER
THE CLASSIFICATION OF PLANAR MONOMIALS OVER FIELDS OF PRIME SQUARE ORDER ROBERT S COULTER Abstract Planar functions were introduced by Dembowski and Ostrom in [3] to describe affine planes possessing collineation
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationbe the set of complex valued 2π-periodic functions f on R such that
. Fourier series. Definition.. Given a real number P, we say a complex valued function f on R is P -periodic if f(x + P ) f(x) for all x R. We let be the set of complex valued -periodic functions f on
More informationExercises Chapter II.
Page 64 Exercises Chapter II. 5. Let A = (1, 2) and B = ( 2, 6). Sketch vectors of the form X = c 1 A + c 2 B for various values of c 1 and c 2. Which vectors in R 2 can be written in this manner? B y
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationApproximate Fekete points and discrete Leja points based on equal area partitions of the unit sphere
Approximate Fekete points and discrete Leja points based on equal area partitions of the unit sphere Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at Computational
More informationPERIODS IMPLYING ALMOST ALL PERIODS FOR TREE MAPS. A. M. Blokh. Department of Mathematics, Wesleyan University Middletown, CT , USA
PERIODS IMPLYING ALMOST ALL PERIODS FOR TREE MAPS A. M. Blokh Department of Mathematics, Wesleyan University Middletown, CT 06459-0128, USA August 1991, revised May 1992 Abstract. Let X be a compact tree,
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationPADÉ INTERPOLATION TABLE AND BIORTHOGONAL RATIONAL FUNCTIONS
ELLIPTIC INTEGRABLE SYSTEMS PADÉ INTERPOLATION TABLE AND BIORTHOGONAL RATIONAL FUNCTIONS A.S. ZHEDANOV Abstract. We study recurrence relations and biorthogonality properties for polynomials and rational
More information