Chordal cubic spline interpolation is fourth order accurate
|
|
- Charlotte Wilson
- 6 years ago
- Views:
Transcription
1 Chordal cubic spline interpolation is fourth order accurate Michael S. Floater Abstract: It is well known that complete cubic spline interpolation of functions with four continuous derivatives is fourth order accurate. In this paper we show that this kind of interpolation, when used to construct parametric spline curves through sequences of points in any space dimension, is again fourth order accurate if the parameter intervals are chosen by chord length. We also show how such chordal spline interpolants can be used to approximate the arc length derivatives of a curve and its length. Keywords: curve parameterization, arc length, spline interpolation, approximation order. 1. Introduction This paper continues a study, started in [ 4 ], of the effect of parameterization on the rate of convergence of parametric curve interpolation. While results for polynomial interpolation of arbitrary degree were obtained in [ 4 ], we begin here with spline interpolation, and study cubic spline interpolation based on the chord length parameterization. While this method for curve fitting has been in widespread use for a long time, and was suggested as early as 1967 by Ahlberg, Nilson, and Walsh [ 1 ], little seems to be known about its approximation order. Roughly speaking, we obtain results analogous to those found in [ 4 ] for cubic polynomials. The main point is that chordal cubic spline interpolation has fourth order accuracy, which provides a sound justification for choosing the chord length parameterization rather than some other parameterization such as uniform or centripetal [ 6 ]: the latter two were found to lead at best to second order accuracy in numerical tests on polynomial interpolation in [ 4 ]. Let f : [a, b] lr d, d 2, be a regular parametric curve, parameterized with respect to arc length, i.e., f is a continuously differentiable function such that f (s = 1, for all s [a, b], where is the Euclidean norm in lr d. We will study cubic spline interpolation with derivative end conditions. Thus suppose that corresponding to some parameter values a = s 0 < s 1 < < s n = b, we are given the points f(s 0,..., f(s n and the (unit derivatives f (s 0 and f (s n. In many applications the values s i are unknown, and so we are forced to choose some appropriate parameter values t 0 < t 1 < < t n, in order to fit a cubic spline interpolant σ : [t 0, t n ] lr d, by which we mean the unique C 2 curve whose restriction to each interval [t i, t i+1 ], i = 0, 1,..., n 1, is a cubic polynomial curve, and such that σ(t i = f(s i, i = 0, 1,..., n, σ (t i = f (1.1 (s i, i = 0, n. 1
2 Our main result is to show that under the chord length parameterization t i+1 t i = f(s i+1 f(s i, i = 0, 1,..., n 1, (1.2 if f C 4 [a, b] then the spline σ converges to the curve f at the rate of h 4 s as h s 0, where h s := max (s i+1 s i. 0 i n 1 We later show that the derivatives of f are approximated by the arc-length derivatives of σ to the same accuracy as for functions. We also show that the length of f is approximated by the length of σ with fourth order accuracy. At the end of the paper we discuss how the derivatives of f are even approximated directly by the derivatives of σ and that the parameterization of σ converges to arc length. However, the convergence orders of the latter approximations are not optimal and can be raised by one by improving the parameterization. 2. Error between curves Our starting point is an estimate for the difference between the length s i = s i+1 s i of the curve piece f [si,s i+1 ] and the chord length f(s i+1 f(s i. It was established in [ 4 ] that if f C 2 [a, b] then for i = 0, 1,..., n 1, 0 s i f(s i+1 f(s i 1 12 ( s i 3 f 2 [s i,s i+1 ], where A = sup x A (x, and so with the values t i of (1.2, s i = O ( ( s i 3 as h s 0. (2.1 This implies, among other things, that for h s small enough, s i /2 s i, and therefore, h s /2 h t h s. Hence, as h s 0, we can freely interchange expressions of the form O(( s i k with O(( k and O(h k s with O(h k t. To measure the distance between the curves σ and f we use the metric [ 7 ] d P (f, σ := inf φ f φ σ [t 0,t n ], (2.2 the infimum taken over continuously differentiable functions φ : [t 0, t n ] [s 0, s n ] with φ (t > 0, t 0 t t n, and φ(t 0 = s 0 and φ(t n = s n. We will prove Theorem 2.1. Suppose σ is the cubic spline interpolant (1.1 based on the chord length parameterization (1.2. If f C 4 [a, b] then d P (f, σ = O ( h 4 s as h s 0. The main idea of the proof is to study the reparameterization g(t := f(φ(t, t [t 0, t n ], (2.3 2
3 where φ : [t 0, t n ] lr is the C 2 cubic spline function satisfying φ(t i = s i, i = 0, 1,..., n, φ (t i = 1, i = 0, n. (2.4 It turns out that when h s is small enough, φ is monotonically increasing, in which case g is well-defined, and also that g σ [t0,t n ] = O ( h 4 s as h s 0. (2.5 To establish both properties, we begin by deriving some bounds on the derivatives of φ. Using the fact that the parameter values t i are chord length, (2.1 implies s i = O ( ( 3 as h s 0, and, dividing by, we deduce that for small enough h s there is some constant C > 0 such that [t i, t i+1 ]φ 1 C( 2, i = 0, 1,..., n 1. (2.6 Lemma 2.2. If (2.6 holds then, for i = 0, 1,..., n 1, φ 1 (ti,t i+1 4C h t, φ (ti,t i+1 6Ch t, φ (ti,t i+1 12C h t. (2.7 Proof: Setting φ i = φ(t i, φ i = φ (t i, and φ i = φ (t i, it is convenient to represent φ in the interval [t i, t i+1 ], i = 0, 1,..., n 1, in terms of its values and second derivatives at the end points, giving φ(t = (1 uφ i + uφ i+1 ( 2 where u = (t t i /. Since 6 ( (2u 3u 2 + u 3 φ i + (u u 3 φ i+1, (2.8 φ (t = [t i, t i+1 ]φ t ( i (2 6u + 3u 2 φ i + (1 3u 2 φ i+1, (2.9 6 the continuity of φ at the nodes t 1,..., t n 1 leads to the familiar linear system in the second derivatives [ 8 ], ([ 1 ], p. 10, φ i 1 + 2φ i φ i+1 = 6[t i 1, t i, t i+1 ]φ, i = 1,..., n 1, (2.10 while the end conditions are captured in the two additional equations 2φ 0 + φ 1 = 6[t 1, t 0, t 1 ]φ, φ n 1 + 2φ n = 6[t n 1, t n, t n+1 ]φ, (2.11 in which t 1 = t 0, hence [t 1, t 0 ]φ = φ (t 0, and similarly t n+1 = t n. The diagonal dominance of the matrix of this linear system implies that φ i 6 max 0 j n [t j 1, t j, t j+1 ]φ, 0 i n. 3
4 Then by (2.6 and the end conditions φ 0 = φ n = 1, [t j 1, t j, t j+1 ]φ [t j, t j+1 ]φ 1 + [t j 1, t j ]φ 1 t j 1 + t j C ( t j ( t j 2 t j 1 + t j Ch t, for 0 j n, hence max 0 i n φ i 6Ch t, (2.12 and since φ is linear in each interval [t i, t i+1 ], this proves the second inequality in (2.7. Next consider the first derivative of φ in [t i, t i+1 ]. Writing equation (2.9 as φ (t 1 = ([t i, t i+1 ]φ 1 t ( i (1 u 2 (2φ i + φ 6 i+1 inequalities (2.6 and (2.12 imply + 2u(1 u(φ i+1 φ i u 2 (φ i + 2φ i+1 φ (t 1 C( 2 + t ( i (1 u 2 + 2u(1 u + u 2 max 2 0 j n φ j 4C h t. Using (2.12 and the linearity of φ in [t i, t i+1 ] we also deduce φ (ti,t i+1 = 1 φ i+1 φ i 1 ( φ i + φ i+1 12C h t., Proof of Theorem 2.1: The equation for φ in (2.7 shows that φ in (2.4 is increasing when h s is small enough and it remains to establish (2.5. Due to (2.4, σ(t i = g(t i, i = 0, 1,..., n, σ (t i = g (t i, i = 0, n, (2.13 and since g C 2 [t 0, t n ], we can bound the error between σ and g using the well-known error analysis for cubic spline interpolation developed by Birkhoff and de Boor [ 2 ], Sharma and Meir [ 8 ], Hall and Meyer [ 5 ], and others, much of which is summarized by de Boor in Chapter V of [ 3 ]. The following error estimate, which derives from the best approximation property of the complete spline interpolation (2.13 is sufficient for our needs: g σ [t0,t n ] 4dist (g, S 0 1, where S1 0 is the space of piecewise linear spline curves [t 0, t n ] lr d over the partition t 0 < t 1 < < t n. This is a straightforward vector-valued generalization of the estimate of ([ 3 ], p. 68. If τ S1 0 is the piecewise linear interpolant to g at the t i we obtain (via the Hermite-Genocchi formula g σ [t0,t n ] 4 g τ [t0,t n ] max ( 2 g (4 (ti,t i+1. ( i n 1
5 Consider then the (local fourth derivative of g in terms of the derivatives of f and φ. The chain rule (with g and φ evaulated at t and f evaluated at φ(t gives g = φ f, g = (φ 2 f + φ f, g = (φ 3 f + 3φ φ f + φ f, g (4 = (φ 4 f (4 + 6(φ 2 φ f + (3(φ 2 + 4φ φ f + φ (4 f. (2.15 Since φ (4 = 0 in each interval (t i, t i+1 and since, by Lemma 2.2, ( φ (ti,t i+1 = O(1, φ (ti,t i+1 = O(1, φ ht (ti,t i+1 = O, as h s 0, we deduce from the expression for g (4 in (2.15 that ( g (4 ht (ti,t i+1 = O, i = 0, 1,..., n 1. (2.16 Hence (2.14 implies that Since ([ 3 ], Chap. V g σ [t0,t n ] = O(h 2 t. (2.17 g σ [t0,t n ] 1 8 h2 t g σ [t0,t n ], this gives d P (f, σ g σ [t0,t n ] = O(h 4 t = O(h 4 s. 3. Error between arc-length derivatives We have seen that the spline curve σ approximates the reparameterization g of the curve f to fourth order for f C 4 [a, b]. Similarly, the derivatives of σ approximate the derivatives of g. Indeed, from inequality (2.17, g (k σ (k [t0,t n ] = O(h 4 k s as h s 0, k = 0, 1, 2. (3.1 The important question though is how to approximate the derivatives of f since these are fixed (arc length derivatives. One approach is to use the arc length derivatives of σ, and we get the same accuracy as for functional interpolation. To see this, first observe that the derivatives of f, being arc length, are the same as the arc length derivatives of g. Thus it is enough to look at the error between the arc length derivatives of σ and the arc length derivatives of g. Denoting arc length derivatives by dots, the first derivatives are ġ(t = g (t g (t, σ(t = σ (t σ (t. 5
6 It follows that ġ σ = g σ g σ + ( σ g g σ. Now since g = φ, equations (2.7 and, with k = 1, (3.1 imply that both g and σ are bounded away from zero as h s 0. Thus from (3.1, f φ σ [t0,t n ] = ġ σ [t0,t n ] = O(h 3 s. Similarly for second derivatives, using the fact that and similarly for σ, we find from (3.1 that g = 1 g 2 g g g g 4 g, f φ σ [t0,t n ] = g σ [t0,t n ] = O(h 2 s. (3.2 Specializing to the interpolation points gives estimates that are independent of φ: Theorem 3.1. If f C 4 [a, b], then f (s i σ(t i = O(h 3 s, f (s i σ(t i = O(h 2 s, as h s 0, i = 0, 1,..., n. Note that the second equation implies that the curvature of σ at t i approximates the curvature of f at s i to order O(h 2 s. We even get approximation in the arc length third derivatives if the global mesh ratio β s := h s / min 0 i n 1 s i is bounded. Let p 1 be the linear polynomial interpolating e at t = t i, t i+1, where e := g σ. Then e (t p 1(t 2 e (4 (ti,t i+1 = 2 g (4 (ti,t i+1, t i < t < t i+1, and since p 1(t = (e (t i+1 e (t i /, it follows from (2.16 and (2.17 that ( h e (ti,t i+1 p 1 (ti,t i g (4 2 (ti,t i+1 = O t. A similar approach to that used to derive (3.2 then gives ( h f φ σ (ti,t i+1 = g 2 σ (ti,t i+1 = O s, s i and so f (s i σ (t i ± = O(β s h s as h s 0, i = 0, 1,..., n. 6
7 4. Error in curve length A well known way of estimating the length L(f of the curve f given only the points f(s 0,..., f(s n is to compute the length n 1 L(P = f(s i+1 f(s i i=0 of the polygon P passing through them. If f C 2 [a, b], the error between L(f and L(P is of order O(h 2 s. In fact Proposition 3.1 of [ 4 ] gives the explicit bound L(P L(f 1 n 1 12 f 2 [a,b] ( s i f 2 [a,b] L(fh2 s. i=0 We next show that the length of the spline interpolant σ gives a higher order estimate. The main point is that the length of f is the same as the length of g, and so it is sufficient to compare the lengths of σ and g. Theorem 4.1. If f C 4 [a, b], then Proof: From (3.1 we have L(σ L(f = O(h 4 s as h s 0. e = O(h 4 s, e = O(h 3 s, (4.1 as h s 0, where e(t := g(t σ(t. Since both σ and g are bounded away from zero for small enough h s, the bound on e in (4.1 implies that L(σ L(g = tn where we have used the identity t 0 tn σ (t g (t dt = 2 t 0 σ g = 2e g + e e σ + g. But since e(t 0 = e(t n = 0, integration by parts implies tn t 0 e (t g(t tn σ (t + g (t dt = e(t d ( t 0 dt e (t g (t σ (t + g (t dt + O( h 6 s, g (t σ (t + g (t dt. Since by (2.15 and (3.1, g, g, σ, and σ are bounded as h s 0, so too is d ( g (t dt σ (t + g (t and the result follows from the bound on e in (4.1. 7
8 5. Error between derivatives We showed in Section 3 that the arc length derivatives of σ approximate the (arc length derivatives of f for f C 4 [a, b]. Interestingly, though, it turns out that the derivatives of σ themselves also approximate the derivatives of f, albeit with the loss of one order of accuracy. We can use the triangle inequality to get the estimate f (k φ σ (k [t0,t n ] f (k φ g (k [t0,t n ] + g (k σ (k [t0,t n ], and the first term on the right-hand side can be estimated from (2.15 and Lemma 2.2. Since the first two equations of (2.15 imply f (φ(t g (t = (1 φ (tf (φ(t, f (φ(t g (t = (1 (φ (t 2 f (φ(t φ (tf(φ(t, (5.1 the bounds on φ in Lemma 2.2 imply and we deduce from (3.1 that Theorem 5.1. If f C 4 [a, b] then f (k φ g (k [t0,t n ] = O(h 3 k s, k = 1, 2, (5.2 f (k φ σ (k [t0,t n ] = O(h 3 k s, k = 1, 2. (5.3 f (s i σ (t i = O(h 2 s, f (s i σ (t i = O(h s, as h s 0, i = 0, 1,..., n. Notice that the case k = 1 in (5.3 implies that σ 1 [t0,t n ] = σ f φ [t0,t n ] = O(h 2 s, (5.4 which means that the parameterization of σ converges to arc length. This makes chordal cubic spline interpolation an attractive method for reparameterizing a given parametric curve, giving an approximate curve whose parameterization is close to arc length. This is certainly not a new idea, but equation (5.4 seems to be the first mathematical justification for this method. 6. Improving the parameterization A natural question raised by the results of the last section is whether there are parameterizations which are even better than the chordal one in the sense that the approximation orders of Theorem 5.1 and (5.4 are one higher. The answer is yes if we use the approach of iterative improvement, studied in [ 4 ]. Suppose we can construct a parameterization ˆt 0 < ˆt 1 < < ˆt n, where each ˆt i is viewed as a function of f(s 0,..., f(s n and f (s 0, f (s n, such that ˆt i s i = O ( ( s i 3 h s, i = 0, 1,..., n 1, (6.1 8
9 and suppose ˆσ : [ˆt 0, ˆt n ] lr d is the cubic spline interpolant resulting from the replacement of the t i in (1.1 by the ˆt i. Then all the bounds on φ in Lemma 2.2 are raised by one power of h s and from (5.1 we can replace (5.2 by and therefore (5.3 by Thus and f (k ˆφ ĝ (k [ˆt 0,ˆt n ] = O(h4 k s, k = 1, 2, f (k ˆφ ˆσ (k [ˆt 0,ˆt n ] = O(h4 k s, k = 1, 2. f (k (s i ˆσ (k 3 (ˆt i = O(h 4 k s, k = 1, 2, i = 0, 1,..., n, ˆσ 1 [ˆt 0,ˆt n ] = O(h3 s. There are at least two ways of finding a parameterization satisfying (6.1. One is to let ˆt i+1 = ˆt i + L(p i [ti,t i+1 ], i = 0, 1,..., n 1, where p i is the (chordal cubic polynomial such that p i (t j = f(s j, j = i 1, i, i + 1, i + 2; a Lagrange interpolant for 1 i n 2 and a Hermite one for i = 0, n 1. Then by Theorem 4.1 of [ 4 ], ˆt i s i = O ( ( s i 3 (s i+1 s i 1 (s i+2 s i, from which (6.1 follows. In practice it would be sufficient [ 4 ] to approximate L(p i [ti,t i+1 ] using for example Simpson s rule or the 2-point Gauss rule to approximate the integral of p i over the interval [t i, t i+1 ]. Another way is to use σ. If we set then ˆt i+1 = ˆt i + L(σ [ti,t i+1 ], i = 0, 1,..., n 1, ˆt i s i = O ( ( s i 3 h 2 s, (6.2 and so (6.1 is satisfied. To see this, we can use the fact that ˆt i s i = L(σ [ti,t i+1 ] L(g [ti,t i+1 ]. By treating each interval [t i, t i+1 ] in (2.17 separately, we obtain the estimates e [ti,t i+1 ] = O(( 2 h 2 t, e [ti,t i+1 ] = O( h 2 t, and (6.2 follows from the same method of proof as for Theorem 4.1. Acknowledgement. I wish to thank the referees for their very thorough reading of the original manuscript and for their help in improving the paper. 9
10 References 1. J. H. Ahlberg, E. N. Nilson, and J. L. Walsh, The theory of splines and their applications, Academic Press, New York, G. Birkhoff and C. de Boor, Error bounds for spline interpolation, J. Math. Mech. 13 (1964, C. de Boor, A practical guide to splines, Springer-Verlag, M. S. Floater, Arc length estimation and the convergence of parametric polynomial interpolation, preprint. 5. C. A. Hall and W. W. Meyer, Optimal error bounds for cubic spline interpolation, J. Approx. Theory 16 (1976, E. T. Y. Lee, Choosing nodes in parametric curve interpolation, Computer-Aided Design 21 (1989, T. Lyche and K. Mørken, A metric for parametric approximation, in Curves and Surfaces, P. J. Laurent, A. Le Méhauté, and L. L. Schumaker (eds., A. K. Peters, Wellesley, (1994, A. Sharma and A. Meir, Degree of approximation of spline interpolation, J. Math. Mech. 15 (1966, Michael S. Floater Department of Informatics Centre of Mathematics for Applications University of Oslo P.B. 1053, Blindern 0316 Oslo, NORWAY michaelf@ifi.uio.no 10
Two chain rules for divided differences
Two chain rules for divided differences and Faà di Bruno s formula Michael S Floater and Tom Lyche Abstract In this paper we derive two formulas for divided differences of a function of a function Both
More informationExtrapolation Methods for Approximating Arc Length and Surface Area
Extrapolation Methods for Approximating Arc Length and Surface Area Michael S. Floater, Atgeirr F. Rasmussen and Ulrich Reif March 2, 27 Abstract A well-known method of estimating the length of a parametric
More informationCurvature variation minimizing cubic Hermite interpolants
Curvature variation minimizing cubic Hermite interpolants Gašper Jaklič a,b, Emil Žagar,a a FMF and IMFM, University of Ljubljana, Jadranska 19, Ljubljana, Slovenia b PINT, University of Primorska, Muzejski
More informationError formulas for divided difference expansions and numerical differentiation
Error formulas for divided difference expansions and numerical differentiation Michael S. Floater Abstract: We derive an expression for the remainder in divided difference expansions and use it to give
More informationApproximation of Circular Arcs by Parametric Polynomial Curves
Approximation of Circular Arcs by Parametric Polynomial Curves Gašper Jaklič Jernej Kozak Marjeta Krajnc Emil Žagar September 19, 005 Abstract In this paper the approximation of circular arcs by parametric
More informationA chain rule for multivariate divided differences
A chain rule for multivariate divided differences Michael S. Floater Abstract In this paper we derive a formula for divided differences of composite functions of several variables with respect to rectangular
More informationAn O(h 2n ) Hermite approximation for conic sections
An O(h 2n ) Hermite approximation for conic sections Michael Floater SINTEF P.O. Box 124, Blindern 0314 Oslo, NORWAY November 1994, Revised March 1996 Abstract. Given a segment of a conic section in the
More informationGeometric Lagrange Interpolation by Planar Cubic Pythagorean-hodograph Curves
Geometric Lagrange Interpolation by Planar Cubic Pythagorean-hodograph Curves Gašper Jaklič a,c, Jernej Kozak a,b, Marjeta Krajnc b, Vito Vitrih c, Emil Žagar a,b, a FMF, University of Ljubljana, Jadranska
More informationCubic Splines MATH 375. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Cubic Splines
Cubic Splines MATH 375 J. Robert Buchanan Department of Mathematics Fall 2006 Introduction Given data {(x 0, f(x 0 )), (x 1, f(x 1 )),...,(x n, f(x n ))} which we wish to interpolate using a polynomial...
More informationBarycentric rational interpolation with no poles and high rates of approximation
Barycentric rational interpolation with no poles and high rates of approximation Michael S. Floater Kai Hormann Abstract It is well known that rational interpolation sometimes gives better approximations
More informationConvergence rates of derivatives of a family of barycentric rational interpolants
Convergence rates of derivatives of a family of barycentric rational interpolants J.-P. Berrut, M. S. Floater and G. Klein University of Fribourg (Switzerland) CMA / IFI, University of Oslo jean-paul.berrut@unifr.ch
More informationNumerische Mathematik
Numer. Math. (997) 76: 479 488 Numerische Mathematik c Springer-Verlag 997 Electronic Edition Exponential decay of C cubic splines vanishing at two symmetric points in each knot interval Sang Dong Kim,,
More informationNonlinear Stationary Subdivision
Nonlinear Stationary Subdivision Michael S. Floater SINTEF P. O. Box 4 Blindern, 034 Oslo, Norway E-mail: michael.floater@math.sintef.no Charles A. Micchelli IBM Corporation T.J. Watson Research Center
More informationC 1 Interpolation with Cumulative Chord Cubics
Fundamenta Informaticae XXI (2001) 1001 16 1001 IOS Press C 1 Interpolation with Cumulative Chord Cubics Ryszard Kozera School of Computer Science and Software Engineering The University of Western Australia
More informationCubic Splines. Antony Jameson. Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305
Cubic Splines Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305 1 References on splines 1. J. H. Ahlberg, E. N. Nilson, J. H. Walsh. Theory of
More informationSelf-Influencing Interpolation in Groundwater Flow
Self-Influencing Interpolation in Groundwater Flow Carolyn Atwood Whitman College Walla Walla, WA Robert Hildebrand University of Puget Sound Tacoma, WA Andrew Homan Ohio Northern University Ada, OH July
More informationOn Parametric Polynomial Circle Approximation
Numerical Algorithms manuscript No. will be inserted by the editor On Parametric Polynomial Circle Approximation Gašper Jaklič Jernej Kozak Received: date / Accepted: date Abstract In the paper, the uniform
More informationNumerical Methods I: Interpolation (cont ed)
1/20 Numerical Methods I: Interpolation (cont ed) Georg Stadler Courant Institute, NYU stadler@cims.nyu.edu November 30, 2017 Interpolation Things you should know 2/20 I Lagrange vs. Hermite interpolation
More informationHigh order parametric polynomial approximation of conic sections
High order parametric polynomial approximation of conic sections Gašper Jaklič a,b,c, Jernej Kozak a,b, Marjeta Krajnc a,b, Vito Vitrih c, Emil Žagar a,b, a FMF, University of Ljubljana, Jadranska 19,
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 informationEXISTENCE OF SET-INTERPOLATING AND ENERGY- MINIMIZING CURVES
EXISTENCE OF SET-INTERPOLATING AND ENERGY- MINIMIZING CURVES JOHANNES WALLNER Abstract. We consider existence of curves c : [0, 1] R n which minimize an energy of the form c (k) p (k = 1, 2,..., 1 < p
More information290 J.M. Carnicer, J.M. Pe~na basis (u 1 ; : : : ; u n ) consisting of minimally supported elements, yet also has a basis (v 1 ; : : : ; v n ) which f
Numer. Math. 67: 289{301 (1994) Numerische Mathematik c Springer-Verlag 1994 Electronic Edition Least supported bases and local linear independence J.M. Carnicer, J.M. Pe~na? Departamento de Matematica
More informationGeometric Interpolation by Planar Cubic Polynomials
1 / 20 Geometric Interpolation by Planar Cubic Polynomials Jernej Kozak, Marjeta Krajnc Faculty of Mathematics and Physics University of Ljubljana Institute of Mathematics, Physics and Mechanics Avignon,
More informationMotion design with Euler-Rodrigues frames of quintic Pythagorean-hodograph curves
Motion design with Euler-Rodrigues frames of quintic Pythagorean-hodograph curves Marjeta Krajnc a,b,, Vito Vitrih c,d a FMF, University of Ljubljana, Jadranska 9, Ljubljana, Slovenia b IMFM, Jadranska
More informationSplines which are piecewise solutions of polyharmonic equation
Splines which are piecewise solutions of polyharmonic equation Ognyan Kounchev March 25, 2006 Abstract This paper appeared in Proceedings of the Conference Curves and Surfaces, Chamonix, 1993 1 Introduction
More informationA generalized mean value property for polyharmonic functions
A generalized mean value property for polyharmonic functions Michael S. Floater December 1, 015 Abstract A well known property of a harmonic function in a ball is that its value at the centreequalsthemeanofitsvaluesontheboundary.
More informationMultiresolution analysis by infinitely differentiable compactly supported functions. N. Dyn A. Ron. September 1992 ABSTRACT
Multiresolution analysis by infinitely differentiable compactly supported functions N. Dyn A. Ron School of of Mathematical Sciences Tel-Aviv University Tel-Aviv, Israel Computer Sciences Department University
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationNodal bases for the serendipity family of finite elements
Foundations of Computational Mathematics manuscript No. (will be inserted by the editor) Nodal bases for the serendipity family of finite elements Michael S. Floater Andrew Gillette Received: date / Accepted:
More informationThe Existence of the Riemann Integral
The Existence of the Riemann Integral James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 18, 2018 Outline The Darboux Integral Upper
More informationArsène Pérard-Gayot (Slides by Piotr Danilewski)
Computer Graphics - Splines - Arsène Pérard-Gayot (Slides by Piotr Danilewski) CURVES Curves Explicit y = f x f: R R γ = x, f x y = 1 x 2 Implicit F x, y = 0 F: R 2 R γ = x, y : F x, y = 0 x 2 + y 2 =
More information3.1 Interpolation and the Lagrange Polynomial
MATH 4073 Chapter 3 Interpolation and Polynomial Approximation Fall 2003 1 Consider a sample x x 0 x 1 x n y y 0 y 1 y n. Can we get a function out of discrete data above that gives a reasonable estimate
More informationEXPLICIT ERROR BOUND FOR QUADRATIC SPLINE APPROXIMATION OF CUBIC SPLINE
J. KSIAM Vol.13, No.4, 257 265, 2009 EXPLICIT ERROR BOUND FOR QUADRATIC SPLINE APPROXIMATION OF CUBIC SPLINE YEON SOO KIM 1 AND YOUNG JOON AHN 2 1 DEPT OF MATHEMATICS, AJOU UNIVERSITY, SUWON, 442 749,
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 informationA sharp upper bound on the approximation order of smooth bivariate pp functions C. de Boor and R.Q. Jia
A sharp upper bound on the approximation order of smooth bivariate pp functions C. de Boor and R.Q. Jia Introduction It is the purpose of this note to show that the approximation order from the space Π
More information1 Definition of the Riemann integral
MAT337H1, Introduction to Real Analysis: notes on Riemann integration 1 Definition of the Riemann integral Definition 1.1. Let [a, b] R be a closed interval. A partition P of [a, b] is a finite set of
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
4 Interpolation 4.1 Polynomial interpolation Problem: LetP n (I), n ln, I := [a,b] lr, be the linear space of polynomials of degree n on I, P n (I) := { p n : I lr p n (x) = n i=0 a i x i, a i lr, 0 i
More informationHermite Interpolation with Euclidean Pythagorean Hodograph Curves
Hermite Interpolation with Euclidean Pythagorean Hodograph Curves Zbyněk Šír Faculty of Mathematics and Physics, Charles University in Prague Sokolovská 83, 86 75 Praha 8 zbynek.sir@mff.cuni.cz Abstract.
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 informationCONTROL POLYGONS FOR CUBIC CURVES
On-Line Geometric Modeling Notes CONTROL POLYGONS FOR CUBIC CURVES Kenneth I. Joy Visualization and Graphics Research Group Department of Computer Science University of California, Davis Overview B-Spline
More informationCubic Splines; Bézier Curves
Cubic Splines; Bézier Curves 1 Cubic Splines piecewise approximation with cubic polynomials conditions on the coefficients of the splines 2 Bézier Curves computer-aided design and manufacturing MCS 471
More informationA RIGOROUS PROOF OF THE ARC LENGTH AND LINE INTEGRAL FORMULA USING THE RIEMANN INTEGRAL
A RGOROUS ROOF OF THE ARC LENGTH AND LNE NTEGRAL FORMULA USNG THE REMANN NTEGRAL ZACHARY DESTEFANO Abstract. n this paper, provide the rigorous mathematical definiton of an arc in general Euclidean Space.
More informationCHAPTER 10 Shape Preserving Properties of B-splines
CHAPTER 10 Shape Preserving Properties of B-splines In earlier chapters we have seen a number of examples of the close relationship between a spline function and its B-spline coefficients This is especially
More informationCurves. Hakan Bilen University of Edinburgh. Computer Graphics Fall Some slides are courtesy of Steve Marschner and Taku Komura
Curves Hakan Bilen University of Edinburgh Computer Graphics Fall 2017 Some slides are courtesy of Steve Marschner and Taku Komura How to create a virtual world? To compose scenes We need to define objects
More informationIntroduction to Computer Graphics. Modeling (1) April 13, 2017 Kenshi Takayama
Introduction to Computer Graphics Modeling (1) April 13, 2017 Kenshi Takayama Parametric curves X & Y coordinates defined by parameter t ( time) Example: Cycloid x t = t sin t y t = 1 cos t Tangent (aka.
More informationFinal Year M.Sc., Degree Examinations
QP CODE 569 Page No Final Year MSc, Degree Examinations September / October 5 (Directorate of Distance Education) MATHEMATICS Paper PM 5: DPB 5: COMPLEX ANALYSIS Time: 3hrs] [Max Marks: 7/8 Instructions
More informationCourse 212: Academic Year Section 1: Metric Spaces
Course 212: Academic Year 1991-2 Section 1: Metric Spaces D. R. Wilkins Contents 1 Metric Spaces 3 1.1 Distance Functions and Metric Spaces............. 3 1.2 Convergence and Continuity in Metric Spaces.........
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 informationNumerical Methods for Differential Equations Mathematical and Computational Tools
Numerical Methods for Differential Equations Mathematical and Computational Tools Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 Part 1. Vector norms, matrix norms and logarithmic
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationDIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION
Meshless Methods in Science and Engineering - An International Conference Porto, 22 DIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION Robert Schaback Institut für Numerische und Angewandte Mathematik (NAM)
More informationKeyframing. CS 448D: Character Animation Prof. Vladlen Koltun Stanford University
Keyframing CS 448D: Character Animation Prof. Vladlen Koltun Stanford University Keyframing in traditional animation Master animator draws key frames Apprentice fills in the in-between frames Keyframing
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation
More informationG-SPLINE INTERPOLATION FOR APPROXIMATING THE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS USING LINEAR MULTI- STEP METHODS
Journal of Al-Nahrain University Vol.0(), December, 00, pp.8- Science G-SPLINE INTERPOLATION FOR APPROXIMATING THE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS USING LINEAR MULTI- STEP METHODS Osama H.
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationResearch Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method
Applied Mathematics Volume 2012, Article ID 605741, 10 pages doi:10.1155/2012/605741 Research Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method M. Zarebnia
More informationWavelets and Image Compression. Bradley J. Lucier
Wavelets and Image Compression Bradley J. Lucier Abstract. In this paper we present certain results about the compression of images using wavelets. We concentrate on the simplest case of the Haar decomposition
More informationBackward divided difference: Representation of numerical data by a polynomial curve
2017; 2(2): 01-06 ISSN: 2456-1452 Maths 2017; 2(2): 01-06 2017 Stats & Maths www.mathsjournal.com Received: 01-01-2017 Accepted: 02-02-2017 Biswajit Das Research Scholar, Assam Down Town University, Panikhaiti,
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationHerz (cf. [H], and also [BS]) proved that the reverse inequality is also true, that is,
REARRANGEMENT OF HARDY-LITTLEWOOD MAXIMAL FUNCTIONS IN LORENTZ SPACES. Jesús Bastero*, Mario Milman and Francisco J. Ruiz** Abstract. For the classical Hardy-Littlewood maximal function M f, a well known
More informationApproximation of Circular Arcs by Parametric Polynomials
Approximation of Circular Arcs by Parametric Polynomials Emil Žagar Lecture on Geometric Modelling at Charles University in Prague December 6th 2017 1 / 44 Outline Introduction Standard Reprezentations
More informationAn interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes
An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationGAKUTO International Series
GAKUTO International Series Mathematical Sciences and Applications, Vol.28(2008) Proceedings of Fourth JSIAM-SIMMAI Seminar on Industrial and Applied Mathematics, pp.139-148 A COMPUTATIONAL APPROACH TO
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationThe Infinity Norm of a Certain Type of Symmetric Circulant Matrix
MATHEMATICS OF COMPUTATION, VOLUME 31, NUMBER 139 JULY 1977, PAGES 733-737 The Infinity Norm of a Certain Type of Symmetric Circulant Matrix By W. D. Hoskins and D. S. Meek Abstract. An attainable bound
More informationLecture 1 INF-MAT : Chapter 2. Examples of Linear Systems
Lecture 1 INF-MAT 4350 2010: Chapter 2. Examples of Linear Systems Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo August 26, 2010 Notation The set of natural
More informationNonlinear Means in Geometric Modeling
Nonlinear Means in Geometric Modeling Michael S. Floater SINTEF P. O. Box 124 Blindern, 0314 Oslo, Norway E-mail: Michael.Floater@math.sintef.no Charles A. Micchelli IBM Corporation T.J. Watson Research
More informationCOMPOSITION SEMIGROUPS ON BMOA AND H AUSTIN ANDERSON, MIRJANA JOVOVIC, AND WAYNE SMITH
COMPOSITION SEMIGROUPS ON BMOA AND H AUSTIN ANDERSON, MIRJANA JOVOVIC, AND WAYNE SMITH Abstract. We study [ϕ t, X], the maximal space of strong continuity for a semigroup of composition operators induced
More informationLocal semiconvexity of Kantorovich potentials on non-compact manifolds
Local semiconvexity of Kantorovich potentials on non-compact manifolds Alessio Figalli, Nicola Gigli Abstract We prove that any Kantorovich potential for the cost function c = d / on a Riemannian manifold
More informationSome Properties of the Augmented Lagrangian in Cone Constrained Optimization
MATHEMATICS OF OPERATIONS RESEARCH Vol. 29, No. 3, August 2004, pp. 479 491 issn 0364-765X eissn 1526-5471 04 2903 0479 informs doi 10.1287/moor.1040.0103 2004 INFORMS Some Properties of the Augmented
More informationA Bernstein-like operator for a mixed algebraic-trigonometric space
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, 21-25 septiembre 2009 (pp. 1 7) A Bernstein-like operator for a mixed algebraic-trigonometric space
More informationSelected solutions for Homework 9
Math 424 B / 574 B Due Wednesday, Dec 09 Autumn 2015 Selected solutions for Homework 9 This file includes solutions only to those problems we did not have time to cover in class. We need the following
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationThe Way of Analysis. Robert S. Strichartz. Jones and Bartlett Publishers. Mathematics Department Cornell University Ithaca, New York
The Way of Analysis Robert S. Strichartz Mathematics Department Cornell University Ithaca, New York Jones and Bartlett Publishers Boston London Contents Preface xiii 1 Preliminaries 1 1.1 The Logic of
More informationELASTIC SPLINES I: EXISTENCE. Albert Borbély & Michael J. Johnson Dept. of Mathematics, Faculty of Science Kuwait University.
ELASTIC SPLINES I: EXISTENCE Albert Borbély & Michael J. Johnson Dept. of Mathematics, Faculty of Science Kuwait University January 21, 214 Abstract. Given interpolation points P 1, P 2,..., P m in the
More informationInterpolation. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Interpolation 1 / 34
Interpolation Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Interpolation 1 / 34 Outline 1 Introduction 2 Lagrange interpolation 3 Piecewise polynomial
More informationA Proof of Markov s Theorem for Polynomials on Banach spaces
A Proof of Markov s Theorem for Polynomials on Banach spaces Lawrence A. Harris Department of Mathematics, University of Kentucky Lexington, Kentucky 40506-007 larry@ms.uky.edu Dedicated to my teachers
More informationarxiv: v2 [math.ds] 9 Jun 2013
SHAPES OF POLYNOMIAL JULIA SETS KATHRYN A. LINDSEY arxiv:209.043v2 [math.ds] 9 Jun 203 Abstract. Any Jordan curve in the complex plane can be approximated arbitrarily well in the Hausdorff topology by
More informationOn the convexity of C 1 surfaces associated with some quadrilateral finite elements
Advances in Computational Mathematics 13 (2000) 271 292 271 On the convexity of C 1 surfaces associated with some quadrilateral finite elements J. Lorente-Pardo a, P. Sablonnière b and M.C. Serrano-Pérez
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationMETRIC HEIGHTS ON AN ABELIAN GROUP
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 44, Number 6, 2014 METRIC HEIGHTS ON AN ABELIAN GROUP CHARLES L. SAMUELS ABSTRACT. Suppose mα) denotes the Mahler measure of the non-zero algebraic number α.
More informationTrajectory Estimation for Exponential Parameterization and Different Samplings
Trajectory Estimation for Exponential Parameterization and Different Samplings Ryszard Kozera 1, Lyle Noakes 2, and Piotr Szmielew 1 1 Warsaw University of Life Sciences - SGGW Faculty of Applied Informatics
More informationLower Bounds for Approximation by MLP Neural Networks
Lower Bounds for Approximation by MLP Neural Networks Vitaly Maiorov and Allan Pinkus Abstract. The degree of approximation by a single hidden layer MLP model with n units in the hidden layer is bounded
More informationReport on Numerical Approximations of FDE s with Method of Steps
Report on Numerical Approximations of FDE s with Method of Steps Simon Lacoste-Julien May 30, 2001 Abstract This is a summary of a meeting I had with Hans Vangheluwe on Thursday May 24 about numerical
More informationExamination paper for TMA4215 Numerical Mathematics
Department of Mathematical Sciences Examination paper for TMA425 Numerical Mathematics Academic contact during examination: Trond Kvamsdal Phone: 93058702 Examination date: 6th of December 207 Examination
More informationBarycentric coordinates for Lagrange interpolation over lattices on a simplex
Barycentric coordinates for Lagrange interpolation over lattices on a simplex Gašper Jaklič gasper.jaklic@fmf.uni-lj.si, Jernej Kozak jernej.kozak@fmf.uni-lj.si, Marjeta Krajnc marjetka.krajnc@fmf.uni-lj.si,
More informationMULTIVARIATE BIRKHOFF-LAGRANGE INTERPOLATION SCHEMES AND CARTESIAN SETS OF NODES. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXXIII, 2(2004), pp. 217 221 217 MULTIVARIATE BIRKHOFF-LAGRANGE INTERPOLATION SCHEMES AND CARTESIAN SETS OF NODES N. CRAINIC Abstract. In this paper we study the relevance
More informationThe Sherrington-Kirkpatrick model
Stat 36 Stochastic Processes on Graphs The Sherrington-Kirkpatrick model Andrea Montanari Lecture - 4/-4/00 The Sherrington-Kirkpatrick (SK) model was introduced by David Sherrington and Scott Kirkpatrick
More information2 The De Casteljau algorithm revisited
A new geometric algorithm to generate spline curves Rui C. Rodrigues Departamento de Física e Matemática Instituto Superior de Engenharia 3030-199 Coimbra, Portugal ruicr@isec.pt F. Silva Leite Departamento
More informationOn a max norm bound for the least squares spline approximant. Carl de Boor University of Wisconsin-Madison, MRC, Madison, USA. 0.
in Approximation and Function Spaces Z Ciesielski (ed) North Holland (Amsterdam), 1981, pp 163 175 On a max norm bound for the least squares spline approximant Carl de Boor University of Wisconsin-Madison,
More informationCURRENT MATERIAL: Vector Calculus.
Math 275, section 002 (Ultman) Fall 2011 FINAL EXAM REVIEW The final exam will be held on Wednesday 14 December from 10:30am 12:30pm in our regular classroom. You will be allowed both sides of an 8.5 11
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationA NOTE ON Q-ORDER OF CONVERGENCE
BIT 0006-3835/01/4102-0422 $16.00 2001, Vol. 41, No. 2, pp. 422 429 c Swets & Zeitlinger A NOTE ON Q-ORDER OF CONVERGENCE L. O. JAY Department of Mathematics, The University of Iowa, 14 MacLean Hall Iowa
More informationGEOMETRIC MODELLING WITH BETA-FUNCTION B-SPLINES, I: PARAMETRIC CURVES
International Journal of Pure and Applied Mathematics Volume 65 No. 3 2010, 339-360 GEOMETRIC MODELLING WITH BETA-FUNCTION B-SPLINES, I: PARAMETRIC CURVES Arne Lakså 1, Børre Bang 2, Lubomir T. Dechevsky
More informationM2R IVR, October 12th Mathematical tools 1 - Session 2
Mathematical tools 1 Session 2 Franck HÉTROY M2R IVR, October 12th 2006 First session reminder Basic definitions Motivation: interpolate or approximate an ordered list of 2D points P i n Definition: spline
More informationSzemerédi s Lemma for the Analyst
Szemerédi s Lemma for the Analyst László Lovász and Balázs Szegedy Microsoft Research April 25 Microsoft Research Technical Report # MSR-TR-25-9 Abstract Szemerédi s Regularity Lemma is a fundamental tool
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationDelay Differential Equations with Constant Lags
Delay Differential Equations with Constant Lags L.F. Shampine Mathematics Department Southern Methodist University Dallas, TX 75275 shampine@smu.edu S. Thompson Department of Mathematics & Statistics Radford
More informationMetrics on the space of shapes
Metrics on the space of shapes IPAM, July 4 David Mumford Division of Applied Math Brown University Collaborators Peter Michor Eitan Sharon What is the space of shapes? S = set of all smooth connected
More information