Research Article | Open Access

Xuqiong Luo, Qikui Du, "Local Lagrange Interpolations Using Bivariate Splines of Degree Seven on Triangulated Quadrangulations", *Journal of Mathematics*, vol. 2013, Article ID 638254, 12 pages, 2013. https://doi.org/10.1155/2013/638254

# Local Lagrange Interpolations Using Bivariate Splines of Degree Seven on Triangulated Quadrangulations

**Academic Editor:**Josefa Linares-Perez

#### Abstract

A local Lagrange interpolation scheme using bivariate splines of degree seven over a checkerboard triangulated quadrangulation is constructed. The method provides optimal order approximation of smooth functions.

#### 1. Introduction

Suppose that is a regular triangulation of a connected polygonal domain . For integers , we define where is the dimensional space of bivariate polynomials of degree . For integers , we further define the following super spline space: where, as usual, means that all polynomials on triangles sharing the vertex have common derivatives up to order at that vertex.

We consider the following Lagrange interpolation problem.

*Problem 1. *Let be a set of points in the plane, and let be a quadrangulation with vertices at the points of . Find a triangulation of and a set of additional points such that, for every choice of the data , there is a unique spline satisfying

We call a Lagrange interpolation set for and and a Lagrange interpolation pair. Although constructing Lagrange interpolation pairs sounds simple at first glance, it is in fact a complex problem, especially since we want a local and stable method which has linear complexity and provides optimal order approximation. In order to construct Lagrange interpolation pairs, both and must be carefully chosen.

For , the first result on local Lagrange interpolation by splines on triangulations was given by NÃ¼rnberger and Zeilfelder [1], where, by subdividing about half of the triangles with a Clough-Tocher split, a local Lagrange interpolation scheme for cubic splines on triangulations whose interior vertices have degree six was developed. NÃ¼rnberger and Zeilfelder [2] constructed a local Lagrange interpolation set for , where is the refining triangulation of Powell-Sabin type (I). NÃ¼rnberger and Zeilfelder [3] used a coloring algorithm to divide all the triangles in into two kinds: white triangles and black triangles, got a new triangulation through refining all the white triangles by the Clough-Tocher refinement, and then gave a local Lagrange interpolation set for . As to more results, the reader is referred to several survey papers [4â€“6].

For , NÃ¼rnberger et al. [7] constructed a local Lagrange interpolation set for . Liu and Fan [8] constructed a local Lagrange interpolation set for , where is the triangulation by refining some of triangles in with the double Clough-Tocher splits.

For general cases , a local Lagrange interpolation of has been proposed by NÃ¼rnberger et al. in [9], where , , and are taken as the following related values
with .

Let be a quadrangulation of which consists of nondegenerate convex quadrilaterals. By adding one or two diagonals of each quadrilateral, some triangulated quadrangulations can be obtained. NÃ¼rnberger et al. [10] constructed a Lagrange interpolation scheme based on cubic splines on certain triangulations obtained from checkerboard quadrangulations. In [11], they also constructed a local Lagrange interpolation method based on cubic splines on certain triangulations obtained from a separable quadrangulation. Further they [12] described local Lagrange interpolation methods based on cubic splines on triangulations obtained from arbitrary strictly convex quadrangulations by adding one or two diagonals. Their construction makes use of a fast algorithm for coloring quadrangulations, and the overall algorithm has linear complexity while providing optimal order approximation of smooth functions.

However, to the authors' knowledge, the local Lagrange interpolation schemes based on splines on any triangulated quadrangulation have not been developed. In this paper, we shall construct a Lagrange interpolation scheme on over the triangulation obtained by adding the two diagonals of each quadrilateral from checkerboard quadrangulations.

The paper is organized as follows. In Section 2 we introduce some notation and describe the Bernstein-BÃ©zier representation of splines. In Section 3 we introduce the checkerboard triangulations. In Section 4 several lemmas of Lagrange minimal determining sets are established. In Section 5 the main results of construction of the Lagrange interpolation pair and error bounds for the interpolating splines are presented.

#### 2. Preliminaries

Throughout the paper we shall make extensive use of the well-known Bernstein-BÃ©zier representation of splines. Let in with vertices , and the corresponding polynomial piece is written in the form: where are the Bernstein basis polynomials of degree associated with . As usual, we identify the Bernstein-BÃ©zier coefficients with the set of domain points . We write for the union of the sets of domain points associated with the triangles of .

Given an integer , let We have similar definitions at the other vertices of . If is a vertex of a triangulation , we, respectively, define that the ring and the disk of radius around are the set where the union is taken over all triangles attached to .

We recall [13] that supposes that is a subspace of , and then is said to be a determining set for that provides that, for any for all implies that . Note that is a linear functional. The set is called a minimal determining set (MDS) for if there is no smaller determining set. Further, following [7], a basis for a spline space is called a stable local basis provided that constants , and exist depending only on the smallest angle in such that(1)for each , there is a vertex of for which ,(2)for all choices of the coefficient vector , Here is defined to be the set of all triangles surrounding vertex , and is defined to be the union of the , where are vertices of .

It is well known that a spline in is uniquely determined by its Bernstein-BÃ©zier coefficient set . In order to describe smoothness conditions for splines, we recall some notations introduced in [14]. Suppose that and are two adjoining triangles from , which share the oriented edge , and let where and are the Bernstein polynomials of degree on the triangles and , respectively. Given integers , let be linear functional defined on by These are called smoothness functionals of order . According to [14], a spline belongs to for some if and only if

#### 3. Checkerboard Triangulations

*Definition 1 (see [10]). *Suppose that is a quadrangulation consisting of quadrilaterals with largest interior angle less than . Suppose that the quadrilaterals can be colored black and white in such a way that any two quadrilaterals sharing an edge have the opposite color. Then we call a checkerboard quadrangulation. The triangulation which is obtained by drawing in both diagonals of all quadrilaterals will be called a checkerboard triangulation.

Let and denote the sets of black and white quadrilaterals of , respectively. Following [10], throughout this paper, we also assume that all interior vertices of are of degree four. This assumption can ensure that there exists such that for every interior vertex of , there is a unique quadrilateral sharing the vertex . For , let be the set of white quadrilaterals which share edges with black quadrilaterals. Let and for , and let be the total number of vertices of . As shown in Figure 1, a typical checkerboard triangulation is displayed in which the quadrilaterals in the set are shaded grey. It is noted that the other black quadrilaterals have not been colored.

#### 4. Lagrange Minimal Determining Sets

Lemma 2 (see [10]). *The set of all domain points in a triangle is a Lagrange minimal determining set for the space .*

Let be a quadrilateral with vertices in counterclockwise order and be the triangulated quadrangulation of with being the intersection point of the two diagonals of . As shown in Figure 2, let and for , where .

Lemma 3. *The set
**
is a stable Lagrange MDS for . These domain points are marked with in Figure 2.*

*Proof. *Equivalently, we consider the related homogenous interpolation problem. By Lemma 2, all Bernstein-BÃ©zier coefficients of associated to 36 domain points in must be zero. By the smoothness conditions on the edge and the smoothness conditions at vertex , all Bernstein-BÃ©zier coefficients of must be zero except for 14 coefficients (), , , and . Let
then the homogenous Lagrange interpolation conditions at the associated 14 domain points lead to
where

It is easy to see that . Thus all the remaining 14 Bernstein-BÃ©zier coefficients of are zero. Similarly, all the Bernstein-BÃ©zier coefficients of must be also zero by smoothness conditions and given Lagrange interpolation conditions. Further, using smoothness conditions along two edges and and smoothness conditions at vertices and , all Bernstein-BÃ©zier coefficients of are zero except for . Using the interpolation condition at domain point : , that is; since . Therefore all B coefficients of are zero.

The construction in Lemma 3 is stable in the sense that the maximum coefficient of is bounded by , where is a constant depending only on the smallest angle in . Thus, we say that is a stable Lagrange MDS for . The proof is completed.

In the following lemmas, we will consider four cases depending on how many edges of adjoin with the other .

Lemma 4. *Suppose thatâ€‰ consists of two triangulated quadrangulations and sharing the edge , where and . Let and be the points where the two diagonals of and intersect, respectively. Let the set
**
where is the set defined in Lemma 3 for the triangulated quadrangulations of and
**
be a stable Lagrange MDS for . These domain points are marked with in Figure 3.*

*Proof. *Using Lemma 3, we can see that all of the B coefficients of associated with domain points in are uniquely determined by the data. Writing in B form and using , , and smoothness conditions, all B coefficients are determined except for the B coefficients associated with . Because , where

It is easy to see that . Thus, all of the remaining B coefficients of are uniquely determined by the data. Then, using Lemma 3 again, it is easy to see that all of the B coefficients of associated with domain points in are uniquely and stable determined by the data. Therefore, is a stable Lagrange MDS for . The proof is completed.

Lemma 5. *Suppose that consists of three triangulated quadrangulations , , and as in Figure 4, where and . Let the set
**
where is the set defined in Lemma 3 for the triangulated quadrangulations of and
**
be a stable Lagrange MDS for . These domain points are marked with in Figure 4.*

*Proof. *Using Lemma 3, we can see that all of the B coefficients of associated with domain points in and are uniquely determined by the data. Arguing as in Lemma 4, all the B coefficients of associated with domain points are uniquely determined by the data too. Writing in B form and using , and smoothness conditions, all B coefficients are determined except for . The Lagrange interpolation and and smoothness conditions at imply that these coefficients must satisfy the linear system , where
where . The determinant of the matrix is
By the geometric meaning of the and , we have and . Thus . So, all B coefficients of are uniquely determined. By the same principle, all B coefficients of are uniquely determined. Using Lemma 3 again, it is easy to see that all B coefficients of are uniquely and stable determined by the data. Therefore is a stable Lagrange MDS for . The proof is completed.

Lemma 6. *Suppose that consists of three triangulated quadrangulations , , and as in Figure 5, where and . Let the set
**
where is the set defined in Lemma 3 for the triangulated quadrangulations of and
**
be a stable Lagrange MDS for . These domain points are marked with in Figure 5.*

*Proof. *Using Lemma 3 and , , and smoothness conditions, we can see that all of the B coefficients of associated with domain points in and are uniquely determined by the data. Writing in B form and using , , and smoothness conditions, all B coefficients are determined except for the B coefficients of associated. The Lagrange interpolation and and smoothness conditions at and imply that these coefficients must satisfy the linear system , where where and . The determinant of the matrix **G** is
By the geometric meaning of the and , we have and . Thus . So, all B coefficients of are uniquely determined. Applying Lemmas 3 and 4, it is easy to see that all B coefficients of are uniquely and stable determined by the data. Therefore is a stable Lagrange MDS for . The proof is completed.

Lemma 7. *Suppose that consists of three triangulated quadrangulations , , , and as in Figure 6, where and . Let the set
**
where is the set defined in Lemma 3 for the triangulated quadrangulations of and
**
be a stable Lagrange MDS for . These domain points are marked with in Figure 6.*

*Proof. *Using Lemma 3 and , , and smoothness conditions, we can see that all of the B coefficients of associated with domain points in , , and are uniquely determined by the data. Writing in B form and using , , and smoothness conditions, all B-coefficients are determined except for the B-coefficients of associated. The Lagrange interpolation and , smoothness conditions at and imply that these coefficients must satisfy the linear system , where