arxiv: v4 [math.na] 20 Oct 2016

Transcription

1 Bézier form of dua bivariate Bernstein poynomias arxiv: v4 [mathna] 20 Oct 2016 Stanisław Lewanowicz a,, Paweł Keer b, Paweł Woźny a a Institute of Computer Science, University of Wrocław, u Joiot-Curie 15, Wrocław, Poand b Facuty of Mathematics and Information Science, Warsaw University of Technoogy, u Koszyowa 75, Warszawa, Poand Abstract Dua Bernstein poynomias of one or two variabes have proved to be very usefu in obtaining Bézier form of the L 2 -soution of the probem of best poynomia approximation of Bézier curve or surface In this connection, the Bézier coefficients of dua Bernstein poynomias are to be evauated at a reasonabe cost In this paper, a set of recurrence reations satisfied by the Bézier coefficients of dua bivariate Bernstein poynomias is derived and an efficient agorithm for evauation of these coefficients is proposed Appications of this resut to some approximation probems of Computer Aided Geometric Design (CAGD) are discussed Key words Dua bivariate Bernstein basis Bézier coefficients Bivariate Jacobi poynomias Bivariate Hahn poynomias Mathematics Subject Cassification (2010) 41A10 41A63 65D17 33C45 1 Introduction and preiminaries The Bernstein-Bézier curves and surfaces have become the standard in the CAGD context due to their favourabe geometric properties (cf [9]) Degree reduction, which consists in approximating a Bézier curve or surface by another one of a ower degree, and approximation of a rationa Bézier curve or surface by a poynomia one have many important appications in geometric modeing, such as data exchange, data compression, and data comparison There have been many papers reevant to this probem (see, eg, [1,2,4,7,10,14,19 26,28,29]) The above approximation probems, with or without constraints, are studied for different choices of the error norm Notice that the Bernstein poynomias do not form an orthogona base, so obtaining the Bézier form of the best L 2 -norm soution requires some effort In a frequenty used approach, the main too appied was transformation between the Bernstein and orthogona poynomia bases Such methods are not ony expensive, but aso may be i-conditioned (cf a remar in [20]) Recenty [28,29], a nove approach to the best L 2 approximation probem, using the dua basis associated with the univariate or bivariate Bernstein basis, was proposed The new methods do not use basis transformation matrices expicity, hence they do not share the abovementioned imitation High efficiency of the methods was obtained thans to the appication of recursive properties of the dua Bernstein poynomias Corresponding author Emai addresses: StanisawLewanowicz@csuniwrocp (Stanisław Lewanowicz), PaweKeer@minipwedup (Paweł Keer), PaweWozny@csuniwrocp (Paweł Woźny)

2 In the present paper, we derive the recurrence reations satisfied by the Bézier coefficients of the dua bivariate Bernstein poynomias and prove in detai a ow-cost agorithm for numerica evauation of these coefficients; see Section 2 This agorithm (without proof) was successfuy used in [15] as a part of the method for poynomia approximation of rationa trianguar Bézier surfaces 1 ; see the discussion in Section 3 Severa usefu properties of the dua bivariate Bernstein basis poynomias, expoited in the proposed method, are obtained in Appendix A In Appendix B, some resuts on the bivariate Jacobi and Hahn orthogona poynomias are coected Beow, we introduce notation and definitions that are used in the paper For y := (y 1,y 2,,y d ) R d, we denote y := y 1 +y 2 ++y d, y := ( y 2 1 +y2 2 ++y2 d )1 2 For n N and c := (c 1, c 2, c 3 ) N 3 such that c < n, we define the foowing sets (cf Figure 1): Θ n := { = ( 1, 2 ) N 2 : 0 n}, Ω c n := { = ( 1, 2 ) N 2 : 1 c 1, 2 c 2, n c 3 }, (11) Γ c n := Θ n \Ω c n c = (0,0,2) c = (0,1,2) c = (2,1,3) Figure 1: Exampes of sets (11) (n = 11) Points of the set Ω c n are mared by white discs, whie the points of the set Γ c n by bac discs Obviousy, Θ n = Ω c n Γ c n Remar 11 In appications to approximation probems reated to trianguar (rationa or poynomia) Bézier patches, the set Θ n corresponds to the set of contro points of a patch, whie its subset Γ c n to the boundary points, where some constraints are to be imposed See Section 3 Let T be the standard triange in R 2, T := {(x 1, x 2 ) : x 1, x 2 0, x 1 +x 2 1} (12) 1 Notice that for a specific choice of parameters, the rationa Bézier surface reduces to a poynomia Bézier surface, so that the above approximation probem, with an additiona assumption, is actuay the poynomia degree reduction probem discussed in [29] In the method proposed there, dua bivariate Bernstein poynomias aso pay the basic roe It shoud be stressed, however, that the present approach is substantiay different and much more efficient than the one proposed in [29] 2

3 For n N and := ( 1, 2 ) Θ n, we denote ( ) n n! := 1! 2!(n )! The shifted factoria is defined for any a C by (a) 0 := 1; (a) := a(a+1) (a+ 1), 1 The Bernstein poynomia basis in Π 2 n, n N, is given by (see, eg, [8], or [9, 184]) ( ) n B n (x) := x 1 1 x 2 2 (1 x )n, := ( 1, 2 ) Θ n, x := (x 1,x 2 ) (13) The (unconstrained) dua bivariate Bernstein basis poynomias [16], are defined so that D n ( ;α) Π2 n, Θ n, (14) D n, Bn α = δ,,, Θ n Here, δ, equas 1 if =, and 0 otherwise, whie the inner product is defined by f,g α := w α (x)f(x)g(x)dx, (15) where the weight function w α is given by T w α (x) := A α x α 1 1 xα 2 2 (1 x )α 3, α := (α 1, α 2, α 3 ), α i > 1, (16) and A α is a normaisation factor (see (B4)) In the seque, we use the notation e (α,n) for the connection coefficients between bivariate Bernstein and dua bivariate Bernstein bases such that D n (x;α) = e (α,n)bn (x), Θ n (17) Θ n Investigating the properties of these coefficients is one of the main goas of the paper The fast recursive scheme for computing the coefficients e (α,n),, Θ n, is formuated in the next section, whie the wide mathematica bacground is given in Section A1 For n N and c := (c 1, c 2, c 3 ) N 3 such that c < n, define the constrained bivariate poynomia space Π c,2 n := { } P Π 2 n : P(x) = x c 1 1 xc 2 2 (1 x )c3 Q(x), where Q Π 2 n c It can be easiy seen that the constrained set of poynomias B n with the range of restricted to Ω c n forms a basis in this space We define the constrained dua bivariate Bernstein basis poynomias, D (n,c) ( ;α) Πn c,2, Ωc n, (18) so that D (n,c), B n α = δ, for, Ω c n, 3

4 where the notation of (15) is used For c = (0,0,0), basis (18) reduces to the unconstrained basis (14) in Π 2 n Obviousy, for any Q Πc,2 n, we have Q(x) = (x) (19) Ω c n Q,D (n,c) α Bn Let E (α,c,n), Ωc n, denote the Bézier coefficients of the constrained bivariate dua Bernstein poynomia D (n,c) (x;α), Ω c n: D (n,c) (x;α) = E (α,c,n)bn (x) (110) Ω c n According to Lemma A7, the foowing formua reates the above coefficients with the unconstrained connection coefficients defined in (17): E (α,c,n) := U V (n)v (n) e c c (α+2c,n c ),, Ω c n, (111) where c := (c 1,c 2 ), and U,V (n), Ω c n, are constants defined in (A19) Hence, for fast evauation of the coefficients (111), it woud be sufficient to have an efficient agorithm for computing the quantities e (α,n) for the fu range, Θ n, with arbitrary parameters α and n In Section 2, we describe such an agorithm that has the computationa compexity order proportiona to the number of these quantities 2 Computing the Bézier coefficients e (α,n) The coefficients e e (α,n),, Θ n, can be arranged in a boc trianguar matrix E (0,n) [ E = E ] E (0,n 1) E (1,n 1) =, (21) Θ n E (0,1) E (1,1) E (n 1,1) E (0,0) E (1,0) E (n 1,0) E (n,0) each boc E, Θ n, being a trianguar matrix E = [ e ] Θ n = e (0,n) e (0,n 1) e (1,n 1) e (0,1) e (1,1) e (n 1,1) e (0,0) e (1,0) e (n 1,0) e (n,0) The proposed agorithm is based on some recurrences satisfied by the eements of the matrix (21) In the seque, we assume that e = 0 if Θ n or Θ n The first recurrence is the foowing (cf Lemma A3): e +v 2 = ( [σ 1 () σ 1 ()]e σ 2()e v 2 +σ 0 ()e +v 2 +σ 2 ()e v 2 )/σ 0 (), (22) 4

5 where v 2 := (0,1), and where for t := (t 1,t 2 ), we define σ 0 (t) := ( t n)(t 2 +α 2 +1), σ 2 (t) := t 2 ( t α 3 n 1), σ 1 (t) := σ 0 (t)+σ 2 (t) Observe that recurrence (22) reates three consecutive bocs of a coumn of the matrix (21), shown in the foowing diagram: E +v 2 E E v 2 According to the convention, the boc E ( 1, 1) has ony zero eements Thus, we can compute a the bocs of this coumn, provided that we have computed the boc E ( 1,0) using another method Now, to initiaise the computation of the first coumn, we need a method to compute the corner boc E (0,0) in the tabe (21) We can use the foowing formua (cf Coroary A5): e (0,0) (α,n) = ( 1) 1 n ( α +3) n 1 Ci h i ( 2 ;α 2,α 3,n 1 ), := ( 1, 2 ) Θ n, (23) n!(α 1 +2) 1 i=0 where the coefficients C i are given by (A17), and we use the notation h i (t;a,b,m) for the univariate Hahn poynomias (cf (B6)) As noticed in Remar B1, we can efficienty evauate the sum in (23) using the Censhaw s agorithm, at the cost of O(n 1 ) operations To compute the remaining part of the ast row of the matrix (21), ie, E (1,0) E (2,0) E (n 1,0) E (n,0), we use the second recurrence (cf Lemma A3): e +v 1 = ( [τ 1 () τ 1 ()]e τ 2()e v 1 +τ 0 ()e +v 1 +τ 2 ()e v 1 )/τ 0 (), (24) where v 1 := (1,0), and for t := (t 1,t 2 ), the coefficients τ j (t) are given by τ 0 (t) := ( t n)(t 1 +α 1 +1), τ 2 (t) := t 1 ( t α 3 n 1), τ 1 (t) := τ 0 (t)+τ 2 (t) The recurrence (24) reates each three consecutive bocs of a row of the matrix (21), shown in the foowing diagram: E v 1 E E +v 1 Again, by the convention, the bocs E ( 1, 2) have ony zero eements The foowing agorithm can be appied to compute the compete set of the coefficients e arranged in the boc matrix E (cf (21)) The agorithm requires O(n 4 ) operations, ie, its cost is proportiona to the tota number of the coefficients Agorithm 21 (Computing the tabe E) Step 1 For 1 = 0,1,,n 1, 2 = 0,1,,n 1, using the Censhaw agorithm, compute e (0,0) ( 1, 2 ) defined by (23), (A17), and put e ( 1, 2 ) (0,0) := e (0,0) ( 1, 2 ) 5

6 Step 2 For 1 = 0,1,,n 1, 1 o for 2 = 0,1,,n 1 1, 1 = 1, 1 +1,,n, 2 = 0,1,,n 1, compute e ( 1, 2 +1) ( 1, 2 ) using the recurrence (22), and put e ( 1, 2 ) ( 1, 2 +1) := e( 1, 2 +1) ( 1, 2 ) ; 2 o for 1 = 1 +1, 1 +2,,n, 2 = 0,1,,n 1, compute e ( 1+1,0) ( 1, 2 ) using the recurrence (24), and put e ( 1, 2 ) ( 1 +1,0) := e( 1+1,0) ( 1, 2 ) Output: The set of the coefficients e (α,n) for, Θ n Remar 22 (i) In Agorithm 21, we made use of the symmetry property e (α,n) = e (α,n) (cf (A4)) (ii) When α 2 = α 3, the cost of competing the tabe E can be reduced significanty using (A5) First, the ower part of the tabe (21), containing the bocs E ( 1, 1 2 (n 1) ) E ( 1,1) E ( 1,0) with 1 = 0,1,,n, is computed using the propery adapted Agorithm 21 In the second step, the upper part of the tabe is computed, using the formua e (α,n) = eˆˆ (α,n), where ˆ := ( 1,n ) and ˆ := ( 1,n ) (iii) Simiar effect of the cost reduction can be obtained if α 1 = α 2, or α 1 = α 3, or α 1 = α 2 = α 3 (cf (A7) (A9)) (iv) Observe that the compexity order of Agorithm 21 equas O(n 4 ), ie, is proportiona to the tota number of the coefficients e (α,c,n) for, Θ n 3 Appications Rationa Bézier curves and surfaces, being the most natura generaization of poynomia Bézier curves and surfaces, are important toos in geometric modeing However, they are sometimes inconvenient in practica appications For this reason, severa agorithms for approximating a rationa Bézier geometric form by a poynomia one have been proposed (see, eg, [10, 18, 21,24]) The foowing constrained approximation probem was recenty considered in [15] 6

7 Probem 31 Given the Bézier coefficients r and positive weights ω, Θ n, of the rationa function R n of degree n, ω r B n (x) Θ R n (x) := n ω B n (x) = r Q n (x), x T, (31) Θ n Θ n where Q n (x) := ω B n (x) ω i Bi n (x), (32) i Θ n find a poynomia of degree m, of the form P m (x) := p B m (x), Θ m x T, (33) with the coefficients p satisfying the conditions p = g for Γ c m, (34) g being prescribed numbers, and c := (c 1,c 2,c 3 ) N 3 being a given parameter vector with c < m, such that the distance between R n and P m, d(r n,p m ) := w α (x)[r n (x) P m (x)] 2 dx, (35) reaches the minimum Ω c m T It has been shown that the soution of Probem 31 is the poynomia (33) with the coefficients given by (cf [15, Thm 22]) p = ( ) m E (α,c,m)( ) u v, Ω c m, (36) where u := ( )( ) n n+m 1 r h I h,, h h+ h Θ n 1 v := ( α +3) 2m h Γ c m ( ) ( m 3 (α i +1) hi + h i )g h i=1 with h 3 := m h, 3 := m, and I h, := T w α (x)q n h (x)bm (x)dx (37) The symbo E (α,c,m) has the meaning given in (110) (see aso (111)) 7

8 Impementation of formua (36) demands 1 o evauation of a the coefficients E (α,c,m) with, Ω c m and 2o computing a the integras I h, with h Θ n, Ω c m The first tas is accompished in two steps (i) We compute a the the coefficients e (µ,c,m),, Ωc m, with M := m c and µ := α+2c by Agorithm 21 and then (ii) use the resut in formua (111) (with n repaced by m) to compute E (α,c,m),, Ωc m Remar that the striing simpicity of Agorithm 21 was obtained thans to using properties of dua bivariate Bernstein poynomias, investigated in Appendix A As for the second tas, observe that in genera, the integras (37) cannot be evauated exacty, and their number, equa to (n+m c )(n+m c +1)/2, may ead to the impression that computation of the coefficients (36) is time-expensive in practice The probem has been defeated in [15], where we have proposed a very effective agorithm which aows to compute the compete set of integras I h,, h Θ n, Ω c m, in a time about twice as ong as the time required to compute a singe integra of this type (the agorithm is an extension of the adaptive, high precision method of [11], and it uses some specia properties of the integras (37)) In this way, we have obtained a very efficient method for computing the soution to the Probem Case of equa weights Notice that in the particuar case where a the weights are equa, ω i = ω, i Θ n, we have Q n (x) = Bn (x), Θ n Consequenty, the rationa function (31) reduces to a poynomia of degree n, so that Probem 31 (with additiona assumption m < n) is actuay the constrained poynomia degree reduction probem that has been discussed [29] In the cited paper, we have deveoped a method to evauate the contro points of the reduced surface in terms of quantities D (m,c), B n, a Ωc m, Θ n Computing these coefficients was the most cumbersome part of the method The main part of the above coection was computed efficienty thans to using some recurrence reations However, a quite arge portion of the coefficients demanded some compex computations By appying the approach of this paper, we may sove the degree reduction probem much faster NB In the case of equa weights, evauation of the integras (37) is a simpe tas as we have (see, eg, [29]) I h, = ( n h )( m ) (α1 +1) h1 + 1 (α 2 +1) h2 + 2 (α 3 +1) n+m h ( α +3) n+m Appendix A: Dua bivariate Bernstein bases A1 Unconstrained dua bivariate Bernstein poynomias In this section, we prove some important properties of the coefficients e (α,n) introduced in (17) Obviousy, we have (cf (19)) e (α,n) = Dn,Dn α,, Θ n (A1) However, we need some aternative formuas that woud be more usefu in obtaining recurrence reations satisfied by these quantities 8

9 In the seque, we adopt the foowing convention: given t := (t 1,t 2 ) Θ n, we use the notation t := (t 2,t 1 ), ˆt := (t 1,t 3 ), ť := (t 3,t 1 ), t := (t 2,t 3 ) and t := (t 3,t 2 ), where t 3 := n t Lemma A1 Dua bivariate Bernstein poynomias have the foowing representation: D n (x;α) = Θ n e (α,n)bn (x), Θ n, where with e (α,n) := 0 i q n C q,i (α,n) := C 2 q,i(α,n)h q,i ( ;α,n)h q,i ( ;α,n) ( ) q 1 i q!( n) q λ q,i (A2) Here H q,i (t;α,n) are the bivariate Hahn poynomias defined by (B5), and λ q,i is the constant given by (B3) Proof In [16], it has been proved that the dua bivariate Bernstein poynomias have the foowing representation: D n (x;α) = b α q,i(n,)pq,i(x), α Θ n, (A3) where 0 i q n b α q,i(n,) := ( 1) i C q,i (α,n)h q,i ( ;α,n) Here Pq,i α, q N; 0 i q, are the two-variabe Jacobi poynomias ([13]; see aso [6, p 86] or (B1)), which are orthonorma on T, ie, Pm, α, Pα n, α equas 1 if (m,) = (n,), and 0 otherwise Using expansion (A3) for both dua poynomias in (A1), we obtain e (α,n) = b α q,i(n,)b α q,i(n,), 0 i q n which is equivaent to (A2) Remar A2 Some symmetry properties of the coefficients e (α,n) shoud be noticed 1 The foowing equation resuts easiy from (A1): e (α,n) = e (α,n) (A4) 2 Let α 2 = α 3 The foowing equaity hods: e (α,n) = eˆˆ (α,n) (A5) This equation can be easiy verified using (A2), definition (B5) of the bivariate Hahn poynomias, and the identity h i (n t ;α 2,α 2,n t 1 ) = ( 1) i h i (t 2 ;α 2,α 2,n t 1 ) (cf (B6)) 9

10 3 The foowing equation hods: e (α,n) = e ( α,n), (A6) where α := (α 2,α 1,α 3 ) By (A1), equation (A6) is equivaent to the equation D n ( ;α),dn ( ;α) α = D ñ ( ; α),dñ ( ; α) which can be easiy verified using the definition (15) of the inner product, α 4 Let α 1 = α 2 By (A6), we obtain the equation α e (α,n) = e (α,n) (A7) 5 Let α 1 = α 3 Then e (α,n) = e (α,n) (A8) Transforming e (α,n), using consecutivey (A6), (A5), and again (A6), the resut foows 6 Let α 1 = α 2 = α 3 By (A5), (A7), and (A8), e (α,n) = e (α,n) = (α,n) = e eˆˆ (α,n) = eǩľ (α,n) = e (α,n) (A9) Lemma A3 The coefficients e e (α,n) defined in (17), and given by (A2), satisfy the foowing bivariate recurrence reations: σ 0 ()e +v 2 σ 1 ()e +σ 2()e v 2 = σ 0 ()e +v 2 σ 1 ()e +σ 2()e v 2, (A10) τ 0 ()e +v 1 τ 1 ()e +τ 2()e v 1 = τ 0 ()e +v 1 τ 1 ()e +τ 2()e v 1, (A11) where v 1 := (1,0), v 2 := (0,1), and for t := (t 1,t 2 ), we define with α := (α 2,α 1,α 3 ), and σ i (t) := ϕ i (α,t), τ i (t) := ϕ i ( α, t), i = 0,1,2, (A12) ϕ 0 (α,t) := ( t n)(t 2 +α 2 +1), ϕ 2 (α,t) := t 2 ( t α 3 n 1), ϕ 1 (α,t) := ϕ 0 (α,t)+ϕ 2 (α,t) (A13) Proof First, we prove the recurrence (A10) By (B5), we have H q,i (t ;α,n) = h i (t 2 ;α 2,α 3,n t 1 )h q i (n t 1 i;α 2 +α 3 +2i+1,α 1,n i) Let D n t 1 t 2 be the difference operator defined according to (B9) Then, by (B8), D n t 1 t 2 h i (t 2 ;α 2,α 3,n t 1 ) = i(i+α 2 +α 3 +1)h i (t 2 ;α 2,α 3,n t 1 ) 10

11 Hence, we obtain D n t 1 t 2 H q,i (t ;α,n) = D n t 1 t 2 h i (t 2 ;α 2,α 3,n t 1 ) h q i (n t 1 i;α 2 +α 3 +2i+1,α 1,n i) = i(i+α 2 +α 3 +1)H q,i (t ;α,n) Consequenty, having in mind the form (A2), we obtain the equation 2 e (α,n) = Dn 1e (α,n) D n 1 which can be simpified to the form (A10) To prove (A11), et us define the difference operator R α t by R α t F(t) = ϕ 0 (α,t)f(t+v 2 ) ϕ 1 (α,t)f(t)+ϕ 2 (α,t)f(t v 2 ), where we use the notation of (A13) The recurrence (A10) can be written as R α e (α,n) = Rα e (α,n) Substituting, and α in pace of, and α, respectivey, gives 2 R α e ( α,n) = R α e ( α,n) (A14) Now, notice that by (A6) we have R α e ( α,n) = ϕ 0 ( α, )e +v 1 ϕ 1 ( α, )e +ϕ 2( α, )e v 1, R α e ( α,n) = ϕ 0 ( α, )e +v 1 ϕ 1 ( α, )e +ϕ 2( α, )e v 1 Using this in (A14), equation (A11) foows Lemma A4 For = ( 1,0), = ( 1, 2 ) Θ n, the foowing formua hods: e (α,n) = G α n max( 1, 1 ) i=0 C i h i ( 2 ;α 2,α 3,n 1 ), (A15) where G α := Γ(η +1)Γ(α 1 +1)/Γ( α +3), η := α 2 +α 3 +1, symbo h i (t;a,b,m) is defined in (B6), and C 0 := c 1, 1 (n,η,α 1 ), C i := (2i+η)( 1 n) i (η +1) i 1 (A16) i!( n) 2 i (α c 1, 3 +1) 1 (n i,2i+η,α 1 ), i 1, i the symbo c h,j (m,α,β) denoting the jth Bézier coefficient of the univariate dua Bernstein poynomia Dh m (x;α,β) (cf [17]) Proof We give a setch of the proof By [16, Thm 33], we have for := ( 1, 2 ) Θ n, D n (x;α) = G α n 1 i=0 f i (n,)(1 x 1 ) i R (α 3,α 2 ) i (x 2 /(1 x 1 ))D n i 1 (x 1 ;2i+η,α 1 ), 11

12 where R (α 3,α 2 ) i are the univariate Jacobi poynomias (cf (B2)), Dj N dua Bernstein poynomias (see, eg, [17]), and (t;µ,ν) are univariate f 0 (n,) := 1, f i (n,) := ( 1) i (2i+η)(η +1) i 1 h i ( 2 ;α 2,α 3,n 1 ), i 1 ( n) i (α 2 +1) i (α 3 +1) i We use the above formua in (A1), then reduce the integration over the triange T to evauating two one-dimensiona integras, and appy the orthogonaity property of the poynomias R (α 3,α 2 ) i (cf [12, 18]) Letting := ( 1,0), the formua (A15) foows Coroary A5 The foowing formua hods: where e (0,0) (α,n) = ( 1) 1 n ( α +3) n 1 Ci h i ( 2 ;α 2,α 3,n 1 ), := ( 1, 2 ) Θ n, n!(α 1 +2) 1 C 0 := (α 1 +2) n (α 2 +α 3 +2) n 1, i=0 Ci := ( 1)i(2i+α 2 +α 3 +1)(α 1 +2) n i ( α +n+3) i, i 1 i!(α 3 +1) i (α 2 +α 3 +i+1) n 1 +1 (A17) Proof The resut foows by putting 1 = 0 in (A15), (A16), and using the expicit form for c 0,1 (n i,2i+α 2 +α 3 +1,α 1 ), given in [17, Eq (211)] A2 Constrained dua bivariate Bernstein poynomias The constrained dua bivariate Bernstein basis poynomias (18) can be expressed in terms of the unconstrained dua bivariate Bernstein poynomias (14) with shifted degree and parameters Namey, we have the foowing resut Lemma A6 ([29]) For Ω c n, the foowing formua hods: where c := (c 1,c 2 ), and D (n,c) (x;α) = U V (n)x c 1 1 xc 2 2 (1 x )c 3 D n c c (x;α+2c), (A18) 3 U := ( α +3) 2 c (α i +1) 1 2c i, V (n) := i=1 ( n c c )( n ) 1 (A19) Lemma A7 The constrained dua bivariate Bernstein poynomias have the Bézier representation D (n,c) (x;α) = E (α,c,n)bn (x), Ω c n where E (α,c,n) := U V (n)v (n) e c c (α+2c,n c ), c := (c 1,c 2 ), notation used being that of (A2) and (A19) 12

13 Proof By Lemma A1, we have D n c c (x;α+2c) = e c Θ n c Putting this resut in (A18) and using the equation x c 1 1 xc 2 2 (1 x )c3 B n c (x) = Ω c n (α+2c,n c )B n c (x) ( n c )( ) n 1 +c B+c n (x), we obtain D (n,c) (x;α) =U V (n) ( )( ) n c n 1 +c e c (α+2c,n c )B+c n (x) Θ n =U V (n) ( )( ) n c n 1 c e c c (α+2c,n c )B n (x) Hence, the emma is proved Appendix B: Bivariate orthogona poynomias The notation ( a1,,a ) r z rf s := b 1,,b s =0 (a 1 ) (a r )!(b 1 ) (b s ) z is used for the generaised hypergeometric series (see, eg, [3, 21]); here r, s Z +, z, a 1,,a r, b 1,, b s C, and (c) is the shifted factoria Reca that two-variabe Jacobi poynomias Pn, α (x), n = 0,1,, = 0,1,,n, are defined by ([13]; see aso [6, p 86]) Pn, α (x) :=λ 1 n, R(2+α 2+α 3 +1,α 1 ) n (x 1 )(1 x 1 ) R (α 3,α 2 ) (x 2 /(1 x 1 )), (B1) where x := (x 1,x 2 ), α := (α 1, α 2, α 3 ), α i > 1, R m (µ,ν) (t) := (µ+1) ( m m,m+µ+ν +1 2F 1 m! µ+1 ) 1 t (B2) is the mth shifted Jacobi poynomia in one variabe [12, 18], and λ 2 n, [ λ α ] 2 (α 1 +1) n (α 2 +1) (α 3 +1) ( +η) n+1 n, :=!(n )!(2 +η)(2n/σ +1)(σ) n+ (B3) with η := α 2 +α 3 +1, σ := α +2 Poynomias (B1) form the orthonorma set with respect to the inner product f,g α := w α (x)f(x)g(x)dx, T 13

14 where T := {(x 1, x 2 ) : x 1, x 2 0, 1 x 1 x 2 0}, and w α (x) := A α x α 1 1 xα 2 2 (1 x 1 x 2 ) α 3 with 3 A α := Γ( α +3) [Γ(α i +1)] 1 (B4) Bivariate Hahn poynomias are defined by (see, eg, [27]) i=1 H q,i (t;α,n) := h i (t 1 ;α 2,α 3,t 1 +t 2 )h q i (t 1 +t 2 i;α 2 +α 3 +2i+1,α 1,n i), (B5) where 0 i q n, n N, α =: (α 1,α 2 α 3 ), α i > 1, i = 1,2,3, t := (t 1,t 2 ) R 2, and h (t;a,b,m) := (a+1) ( M) 3 F 2 (,+a+b+1, t a+1, M ) 1 (B6) are the univariate Hahn poynomias (see, eg, [12, 15]) The atter poynomias satisfy the recurrence reation h +1 (t) = A (t,m)h (t)+b (M)h 1 (t), 0; h 0 (t) 1; h 1 (t) 0, where h (t) h (t;a,b,m), A (t,m) := C (2 +s 1) 2 t D E, B (M) := D E 1, (B7) with s := a+b+1, C := (2 +s+1)/[( +s)(2 +s 1)], D := C ( +M +s)( +b), and E := (+a+1)(m ) Moreover, the poynomias (B6) satisfy the difference equation D M t h j (t;a,b,m) = j(j +a+b+1)h j (t;a,b,m), (B8) where the difference operator D M t is given by D M t F(t) := U(t;a,b,M)F(t +1) V(t;a,b,M)F(t) +W(t;a,b,M)F(t 1) (B9) with U(t;a,b,M) := (t M)(t+a+1), W(t;a,b,M) := t(t b M 1) and V(t;a,b,M) := U(t;a,b,M)+W(t;a,b,M) Remar B1 A inear combination s N (t) := N i=0 γ ih i (t;a,b,m) can be summed up using the foowing Censhaw s agorithm (see, eg, [5, Thm 3211]) Compute the sequence V 0,V 1,, V n+2 from V i := γ i +A i (t;m)v i+1 +B i+1 (M)V i+2, i = N,N 1,,0, with V N+1 = V N+2 = 0, where the coefficients A i (t;m) and B i (M) are defined by (B7) Then, s N (t) = V 0 References [1] YJ Ahn, Using Jacobi poynomias for degree reduction of Bézier curves with C - constraints, Comput Aided Geom Design 20 (2003)

15 [2] YJ Ahn, B-G Lee, Y Par, J Yoo, Constrained poynomia degree reduction in the L 2 -norm equas best weighted Eucidean approximation of Bézier coefficients, Comput Aided Geom Design 21 (2004) [3] GE Andrews, R Asey, R Roy, Specia Functions, Cambridge Univ Press, Cambridge, 1999 [4] G Chen, G Wang, Optima muti-degree reduction of Bézier curves with constraints of endpoints continuity, Comput Aided Geom Design 19 (2002) [5] G Dahquist, A Björc, Numerica Methods in Scientific Computing, vo I, SIAM, Phiadephia, 2008 [6] CF Dun, Y Xu, Orthogona Poynomias of Severa Variabes, Cambridge Univ Press, Cambridge, 2001 [7] M Ec, Least squares degree reduction of Bézier curves, Computer-Aided Design 27 (1995) [8] G Farin, Trianguar Bernstein-Bézier patches, Comput Aided Geom Design 3 (1986) [9] G Farin, Curves and Surfaces for Computer-Aided Geometric Design A Practica Guide, fifth ed, Academic Press, Boston, 2002 [10] QQ Hu, An iterative agorithm for poynomia approximation of rationa trianguar Bézier surfaces, App Math Comput 219 (2013) [11] P Keer, A method for indefinite integration of osciatory and singuar functions, Numer Agor 46 (2007) [12] R Koeoe, RF Swarttouw, The Asey scheme of hypergeometric orthogona poynomias and its q-anaogue, Rep 98-17, Fac Techn Math Informatics, Deft Univ of Technoogy, Deft, 1998 [13] TH Koornwinder, Two-variabe anaogues of the cassica orthogona poynomias, in: RA Asey (ed), Theory and Appication of Specia Functions, Academic Press, New Yor, 1975, pp [14] B-G Lee, Y Par, J Yoo, Appication of Legendre-Bernstein basis transformations to degree eevation and degree reduction, Comput Aided Geom Design 19 (2002) [15] S Lewanowicz, P Keer, P Woźny, Constrained approximation of rationa trianguar Bézier surfaces by poynomia trianguar Bézier surfaces Numer Agor 75 (2017), in press; DOI /s [16] S Lewanowicz, P Woźny, Connections between two-variabe Bernstein and Jacobi poynomias on the triange, J Comput App Math 197 (2006) [17] S Lewanowicz, P Woźny, Bézier representation of the constrained dua Bernstein poynomias, App Math Comput 218 (2011) [18] S Lewanowicz, P Woźny, P Keer, Poynomia approximation of rationa Bézier curves with constraints, Numer Agor 59 (2012) [19] L Lu, G Wang, Optima degree reduction of Bézier curves with G 2 -continuity, Comput Aided Geom Design 23 (2006)

16 [20] L Lu, G Wang, Appication of Chebyshev II Bernstein basis transformations to degree reduction of Bézier curves, J Comput App Math 221 (2008) [21] A Rababah, Distances with rationa trianguar Bézier surfaces, App Math Comp 160 (2005) [22] A Rababah, B-G Lee, J Yoo, A simpe matrix form for degree reduction of Bézier curves using Chebyshev Bernstein basis transformations, App Math Comput 181 (2006) [23] A Rababah, M A-Natour, The weighted dua functionas for the univariate Bernstein basis, App Math Comput 186 (2006) [24] R Sharma, Conversion of a rationa poynomia in trianguar Bernstein-Bézier form to poynomia in trianguar Bernstein-Bézier form, Internat J Comput App Math 8 (2013) [25] H Sunwoo, Matrix representation for muti-degree reduction of Bézier curves, Comput Aided Geom Design 22 (2005) [26] H Sunwoo, N Lee, A unified matrix representation for degree reduction of Bézier curves, Comput Aided Geom Design 21 (2004) [27] MV Tratni, Some mutivariate orthogona poynomias of the Asey tabeau discrete famiies, J Math Physics 32 (1991) [28] P Woźny, S Lewanowicz, Muti-degree reduction of Bézier curves with constraints, using dua Bernstein basis poynomias, Comput Aided Geom Design 26 (2009) [29] P Woźny, S Lewanowicz, Constrained muti-degree reduction of trianguar Bézier surfaces using dua Bernstein poynomias, J Comput App Math 235 (2010)