Recurrence Relations for Orthogonal Polynomials on Triangular Domains

: In Farouki et al , 2003, Legendre-weighted orthogonal polynomials P n , r ( u , v , w ) , r = 0, 1, . . . , n , n ≥ 0 on the triangular domain T = { ( u , v , w ) : u , v , w ≥ 0, u + v + w = 1 } are constructed, where u , v , w are the barycentric coordinates. Unfortunately, evaluating the explicit formulas requires many operations and is not very practical from an algorithmic point of view. Hence, there is a need for a more efﬁcient alternative. A very convenient method for computing orthogonal polynomials is based on recurrence relations. Such recurrence relations are described in this paper for the triangular orthogonal polynomials, providing a simple and fast algorithm for their evaluation.


Introduction
Orthogonal polynomials are very strong tools in approximation theory and play a fundamental role in finding the least-squares approximation in explicit form.They are well-studied, and a lot of research has been done for the univariate case, see [1].Orthogonal polynomials over a square region can be constructed using the tensor product of univariate orthogonal polynomials, see [2,3].Orthogonal polynomials over triangular domains have to be determined in a different way, see [4].A scheme for orthogonal bivariate polynomials on triangular regions is constructed in [5].The Legendre weight function W(u, v, w) = 1 is used; these Legendre-weighted orthogonal polynomials are given in the Bernstein basis form, in order to take full advantage of the numerical stability property of the later basis.In [5], to find these bivariate orthogonal polynomials on triangular domains, the sums have to be expanded, and many multiplications have to be performed.This leads to a time-consuming and costly method of construction.There is an essential need for an easy method to find these orthogonal polynomials.It is known that the most convenient method to compute orthogonal polynomials is using a recurrence relation.For the case under consideration, such an algorithm is described in this paper in Section 6.
Connections between bivariate Bernstein and Jacobi bases are considered in [6] without considering the recurrence relations for these polynomials.Orthogonal polynomials on the interior of the triangle based on second-order linear partial differential equations is considered in [7].The approach is different when considering the boundary of the triangle which needs applying the directional derivatives.
This paper is organized as follows.The univariate Legendre and Bernstein polynomials are introduced in Section 2. In Sections 3 and 4, bivariate polynomials and bivariate polynomials on triangular domains are introduced.Legendre-weighted orthogonal polynomials are given in Section 5, and their recurrence relations are given in Section 6.We end this paper with conclusions in Section 7.

Univariate Legendre and Bernstein Polynomials
The univariate Legendre polynomials, L n (x), of degree n are orthogonal on the interval [−1, 1] with respect to the weight function W(x) = 1.They are traditionally defined over [−1, 1] to highlight their symmetry, while the Bernstein polynomials, B n i (u), are defined over [0, 1] to highlight their geometric behaviour.Since the final results will be given in Bernstein form, the suitable interval for both polynomials is [0, 1].In this case, the shifted Legendre polynomials For u ∈ [0, 1], the Bernstein polynomial basis is defined as, see [8], For u ∈ [0, 1], the univariate Legendre polynomials satisfy the following three term recurrence relation, see [1]: The polynomials L n satisfy the following orthogonality conditions:

Bivariate Orthogonal Polynomials
Let G ⊆ R 2 be a finite simply connected domain bounded by a curve Γ.

The integrals
are finite and called moments of W. Similar to the univariate case, weight functions are used to define bivariate orthogonal polynomials.A system of polynomials P n,i (x, y), n ≥ 0, i = 0, 1, . . .n is called orthogonal over G with respect to the weight function W(x, y) if all the leading coefficients c n,i of P n,i (x, y) are positive and where δ is the Kronecker delta, see [9,10].Like the univariate case, for any weight function W over a domain G there exists a system of orthogonal polynomials.Up to this normalization, the system of orthogonal polynomials is unique.For n ≥ 1, we define the space L n of polynomials of degree n that are orthogonal to all polynomials of degree < n over the domain G, i.e., and Π n−1 to be the space of all polynomials of degree n − 1 over the domain G.These polynomials are usually represented in a triangular table : P 0,0 (x, y) P 1,0 (x, y), P 1,1 (x, y) P 2,0 (x, y), P 2,1 (x, y), P 2,2 (x, y) . . .

Bivariate Polynomials on Triangular Domains
The bivariate polynomials on triangular domains are used in design and modelling in Computer Aided Geometric Design; they are defined using the barycentric coordinates.Let p k = (x k , y k ), k = 1, 2, 3 be the vertices (non-collinear) of a triangle T and p be a point, all in the plane.It is always possible to express p as a barycentric combination of p k , k = 1, 2, 3 in the form p=up 1 + vp 2 + wp 3 , where u, v, w ≥ 0, u + v + w = 1, see [14] and the references therein.Using Cramer's rule, the barycentric coordinates are given as the ratio of areas of subtriangles of the base triangle in explicit form: Consider the triangular region T that is given by the barycentric coordinates as follows and the triples α = (i, j, k) of non-negative integers with |α| = i + j + k.A system of polynomials {P n,r (u, v, w), r = 0, 1, . . .n} ∞ n=0 is orthogonal over T with respect to the weight function W(u, v, w) = 1 if all the leading coefficients are positive and T W(u, v, w)P n,r (u, v, w)P m,s (u, v, w)dA = δ n,m δ r,s .
In the following section, the Legendre weight function W(u, v, w) = 1 is considered.

Legendre-Weighted Orthogonal Polynomials
Farouki, Goodman, and Sauer determined in [5] closed-form representation of degree-ordered system of Legendre-weighted orthogonal polynomials P n,r (u, v, w), r = 0, 1, . . ., n, n ≥ 0 on the triangular domain T = {(u, v, w : u, v, w ≥ 0, u + v + w = 1}.These orthogonal polynomials are orthogonal with respect to the weight function W(u, v, w) = 1.Since the Bernstein polynomials are stable, it is convenient to write these polynomials in Bernstein form.These polynomials P n,r (u, v, w) are orthogonal to each polynomial of degree ≤ n − 1 and also orthogonal to each polynomial P n,s (u, v, w), r = s.They are given explicitly by: This formula is further simplified in [5] to: where L r , r = 0, 1, . . ., n are the Legendre polynomials of degree r and Traditionally, univariate Legendre polynomials are computed using different methods.The most convenient method is the recurrence relation in Equation (2).The method of construction in [5] is complicated and time-consuming, since many computations are required.Hence, there is a need for a simple evaluation method, which we present in the next section.

Recurrence Relation
Any three consecutive univariate orthogonal polynomials are related by a recurrence relation.In this section, it is shown that the bivariate orthogonal polynomials in Equation ( 5) also possess a recurrence relation.Any bivariate orthogonal polynomial from the kth row is related by a recurrence relation to two orthogonal polynomials from the preceding two rows for all r = n − 1.If r = n − 1, then P n,n−1 (u, v, w) is related by a recurrence relation to two orthogonal polynomials from the k − 2nd and k − 3rd rows.The exact formulas are given in the following theorem.
Theorem 1.The bivariate orthogonal polynomials P n,r (u, v, w), r = 0, 1, . . ., n given by Equation ( 5) are computed using the following recurrence relations: Proof.The construction of P n,r (u, v, w) suggests that the recurrence relation depends on the values of r.Therefore, we have to distinguish between the three cases r = n, r = n − 1, 0 ≤ r ≤ n − 2. In all cases, we begin with Equation ( 6) and do the proper arrangements and implications to get the required recurrence relations as follows: The case r = n: We begin by substituting r = n in Equation ( 6) to get .
Using Equation (4.3.2) in [1], we get the univariate Jacobi polynomials in terms of the Bernstein polynomials in the following form: and thus q n,r (w) satisfies the following recurrence relation: Substituting the last relation in Equation ( 6) gives: This is further simplified to: Substituting α = 0, β = 2r + 1, n = n − r − 1 in Equation ( 11), we get the following two equations: Thus the last equation for P n,r (u, v, w) is simplified to q n−1,r (w) (n − r)(n + r + 1)(2n − 1) q n−2,r (w).
Applying Equation ( 8) twice, we get the required result.This completes the proof of the Theorem.

Conclusions
In this paper, recurrence relations for the Legendre-weighted orthogonal polynomials P n,r (u, v, w), r = 0, 1, . . ., n, n ≥ 0 on the triangular domain T = {(u, v, w) : u, v, w ≥ 0, u + v + w = 1} are presented.Unlike the method of construction in [5], our recursive algorithm enables us to compute these polynomials efficiently in a very simple fashion.
A function W(x, y) ≥ 0 is called a weight function over G if G W(x, y)dxdy exists and is positive.The weight function W induces a function space L 2 (W, G) with the inner product: f , g = G W(x, y) f (x, y)g(x, y)dxdy.