arxiv: v3 [math.na] 1 Dec 2015
|
|
- Alaina Lucas
- 6 years ago
- Views:
Transcription
1 arxiv: v3 [math.na] Dec 5 Constrained approximation of rational triangular Bézier surfaces by polynomial triangular Bézier surfaces Stanisław Lewanowicz a,, Paweł Keller b, Paweł Woźny a a Institute of Computer Science, University of Wrocław, ul. Joliot-Curie 5, Wrocław, Poland b Faculty of Mathematics and Information Science, Warsaw University of Technology, ul. Koszyowa 75, -66 Warszawa, Poland Abstract. We propose a novel approach to the problem of polynomial approximation of rational Bézier triangular patches with prescribed boundary control points. The method is very efficient thans to using recursive properties of the bivariate dual Bernstein polynomials and applying a smart algorithm for evaluating a collection of two-dimensional integrals. Some illustrative examples are given. Key words: Rational triangular Bézier surface; Polynomial approximation; Bivariate dual Bernstein basis; Two-dimensional integral; Adaptive quadrature. Introduction and preliminaries Rational triangular Bézier surfaces are an important tool in computer graphics. However, they may be sometimes inconvenient in practical applications. The reason is that evaluation of integrals or derivatives of rational expressions is cumbersome. Also, it happens that a rational surface produced in one CAD system is to be imported into another system which can handle only polynomial surfaces. In order to solve the two problems above, different algorithms for approximating the rational surface by polynomial surface are proposed [3, 8, 7, ]. The spectrum of methods contains hybrid algorithm [], progressive iteration approximation [3, 8], least squares approximation and linear programming [7], and approximation by Bézier surfaces with control points obtained by successive degree elevation of the rational Bézier surface [, ]. As a rule, no geometric constraints are imposed, which means a serious drawbac: if we start with a patchwor of smoothly connected rational Bézier triangles and approximate each patch separately, we do not obtain a smooth composite surface. In this paper, we propose a method for solving the problem of the constrained least squares approximation of a rational triangular Bézier patch by a polynomial triangular Bézier patch; see Problem. below. The method is based on the idea of using constrained dual bivariate Bernstein polynomials. Using a fast recursive scheme of evaluation of Bézier form coefficients of the bivariate dual Bernstein polynomials, and applying a swift adaptive scheme of numerical computation of a collection of double integrals involving rational functions resulted in high efficiency of the method. Corresponding author addresses: Stanislaw.Lewanowicz@ii.uni.wroc.pl (Stanisław Lewanowicz), Pawel.Keller@mini.pw.edu.pl (Paweł Keller), Pawel.Wozny@ii.uni.wroc.pl (Paweł Woźny)
2 The outline of the paper is as follows. Section brings a complete solution to the approximation problem. Some comments on the algorithmic implementation of the method are given in Section 3; some technical details of the implementation are presented in Appendix A. In Section 4, several examples are given to show the efficiency of the method. In Appendix B, some basic information on the Hahn orthogonal polynomials is recalled. We end this section by introducing some notation. For y := (y,y,...,y d ) R d, we denote y := y +y +...+y d, and y := ( y +y +...+y d For n N and c := (c, c, c 3 ) N 3 such that c < n, we define the following sets (cf. Figure ): Θ n := { = (, ) N : n}, Ω c n := { = (, ) N : c, c, n c 3 }, (.) Γ c n := Θ n \Ω c n. ) c = (,,) c = (,,) c = (,,3) Figure : Examples of sets (.) (n = ). Points of the set Ω c n are mared by white discs, while the points of the set Γ c n by blac discs. Obviously, Θ n = Ω c n Γ c n. Throughout this paper, the symbolπ n denotes the space of all polynomials of two variables, of total degree at most n. Let T be the standard triangle in R, T := {(x, x ) : x, x, x +x }. (.) For n N, and := (, ) Θ n, we denote, ( ) n n! :=!!(n )!. The shifted factorial is defined for any a C by (a) := ; (a) := a(a+) (a+ ),. The Bernstein polynomial basis in Π n, n N, is given by (see, e.g., [5], or [6, 7.3]), ( ) n B n (x) := x x ( x )n, := (, ) Θ n, x := (x,x ). (.3)
3 The (unconstrained) bivariate dual Bernstein basis polynomials [], are defined so that D n ( ;α) Π n, Θ n, (.4) D n, Bn l α = δ,l,,l Θ n. Here δ,l equals if = l, and otherwise, while the inner product is defined by f,g α := w α (x)f(x)g(x)dx, (.5) where the weight function w α is given by T w α (x) := A α x α xα ( x )α 3, α := (α, α, α 3 ), α i >, (.6) with A α := Γ( α +3)/[Γ(α +)Γ(α +)Γ(α 3 +)]. For n N and c := (c, c, c 3 ) N 3 such that c < n, define the constrained bivariate polynomial space Π c, n := { } P Π n : P(x) = xc xc ( x )c3 Q(x), where Q Π n c. It can be easily seen that the constrained set {B n} Ω c of degree n bivariate Bernstein polynomials forms a basis in this space. We define constrained dual bivariate Bernstein basis n polynomials, D (n,c) ( ;α) Πn c,, Ω c n, (.7) so that D (n,c), B n l α = δ,l for, l Ω c n, (.8) where the notation of (.5) is used. For c = (,,), basis (.7) reduces to the unconstrained basis (.4) in Π n. Notice that the solution of the least squares approximation problem in the space Π n (c,) result. can be given in terms of the polynomials D (n,c). Namely, we have the following Lemma. Let F be a function defined on the standard triangle T (cf. (.)). The polynomial S n Π n (c,), which gives the minimum value of the norm F S n L := F S n,f S n α, is given by S n = Ω c n F,D (n,c) α Bn. (.9) Proof. Π n (c,) : Obviously, S n has the following representation in the Bernstein basis of the space S n = Ω c n S n,d (n,c) α Bn. 3
4 On the other hand, a classical characterization of the best approximation polynomial S n is that F S n,q α = holds for any polynomial Q Π n (c,) (see, e.g. [4, Thm 4.5.]). In particular, for Q = D (n,c), we obtain F,D (n,c) α = S n,d (n,c), α Ωc n. Hence, the formula (.9) follows. The coefficients El (α,c,n) in the Bézier form of the dual Bernstein polynomials, D (n,c) (x;α) = El (α,c,n)bn l (x), Ωc n, (.) l Ω c n play important role in the proposed method. Using the duality property (.8), we obtain the following expression for the coefficients of the above expansion: El (α,c,n) = D (n,c),d (n,c) l. (.) α In a recent paper [], an efficient algorithm was obtained for evaluating all these coefficients for, l Ω c n, with the computational complexity O(n 4 ), i.e., proportional to the total number of these coefficients. See Section 3. for details. Polynomial approximation of Bézier triangular surfaces with constraints In this paper, we consider the following approximation problem. Problem. Let R n be a rational triangular Bézier surface of degree n, ω r B R n (x) := Q n (x) n(x) ω(x) = Θ n ω B n (x), x T, Θ n with the control points r R 3 and positive weights ω R, Θ n. Find a Bézier triangular surface of degree m, P m (x) := p B m (x), x T, Θ m with the control points p R 3, satisfying the conditions p = g for Γ c m, (.) g R 3 being prescribed control points, and c := (c,c,c 3 ) N 3 being a given parameter vector with c < m, such that the distance between the surfaces R n and P m, d(r n,p m ) := w α (x) R n (x) P m (x) dx, (.) reaches the minimum. T 4
5 Remar. Remember that continuity conditions for any two adjacent triangular Bézier patches are given in terms of several rows of the control net "parallel" to the control polygon of their common boundary (see, e.g., [6, Section 7]). Therefore, constraints (.) are natural, in a sense (cf. Fig. ). In Section 4, we give several examples of practical usage of this approach. Clearly, the Bézier triangular patch being the solution of Problem. can be obtained in a componentwise way. Hence it is sufficient to give a method for solving the above problem in the case where R n and P m are scalar functions, and g are numbers. All the details of the proposed method are given in the following theorem. Theorem.3 Given the coefficients r and positive weights ω, Θ n, of the rational function ω r B R n (x) := Q n (x) n(x) ω(x) = Θ n ω B n (x), (.3) Θ n the coefficients p of the degree m polynomial minimising the error with the constraints are given by where p = P m (x) := p B m (x), (.4) Θ m R n P m L := R n P m,r n P m α, (.5) l Ω c m h Θ n p = g for Γ c m, (.6) ( ) m El l (α,c,m)( ) u l v l, Ω c m, (.7) u l := ( )( ) n n+m ω h r h I h+l, h h+l v l := ( α +3) m with h 3 := m h, l 3 := m l, and I j := T h Γ c m The symbol El (α,c,m) has the meaning given in (.). ( ) ( m 3 (α i +) hi +l h i )g h, i= w α (x) Bn+m j (x) dx, j Ω c n+m ω(x). (.8) 5
6 Proof. Observe that R n P m L = W S m L where W := R n T m, T m := g B m, S m := p B m, Γ c m Ω c m the notation being that of (.). Thus, we want S m to be the best approximation polynomial for the function W in the space Π c, m. Its Bézier coefficients are given by p = α = ) El ( R (α,c,m) n,bl m α T m,bl m α, Ω c m, W,D (m,c) l Ω c m where we have used Lemma.. We obtain R n, Bl m α = B n ω h r h h ω, Bm l h Θ n α = ( )( )( ) n m n+m ω h r h h l h+l ω, Bn+m h+l h Θ n = ( )( )( ) n m n+m ω h r h I h+l, h l h+l h Θ n where we use the notation (.8). Further, using equations (.3) and (.5), we obtain T m,bl m α = g h Bh m,bm l α h Γ c m = h Γ c m g h ( m h )( ) m (α +) h +l (α +) h +l (α 3 +) m h l. l ( α +3) m α Hence, the formula (.7) follows. Remar.4 In general, the integrals (.8) cannot be evaluated exactly. In Section 3., we show that they can be efficiently computed numerically up to high precision using an extension of the method of [9]. In the special case where all the weights ω i, i Θ n, are equal, the rational function (.3) reduces to a polynomial of degree n, so that the problem is actually the constrained polynomial degree reduction problem (see, e.g., [8]). Evaluation of the integrals is then a simple tas. 3 Implementation of the method In this section, we discuss some computational details of the polynomial approximation of the rational Bézier function, described in Section (see Theorem.3). 6
7 3. Computing the coefficients E l (α,c,m) We have to compute all the coefficients El (α,c,m) with, l Ωc m. It has been shown [] that they can be given in terms of e l (µ,m) := DM,DM l µ,,l Θ M, with M := m c and µ := α + c, where D M DM ( ;µ) are the unconstrained dual Bernstein polynomial of total degree M (cf. (.4)). See Eq. (3.) for details. Obviously, e l (µ,m) = el (µ,m). The following algoritm is based on the recurrence relations satisfied by e l e l (µ,m), obtained in the paper cited above. Algorithm 3. (Computing the coefficients E l (α,c,m)) Step Let M := m c, µ := α+c. Step For l =,,...,M, l =,,...,M l, compute and e l := ( )l ( µ +3) M M!(µ +) l M l i= C i h i (l ;µ,µ 3,M l ), (3.) where = (,), l = (l,l ), h i (t;a,b,n) are the Hahn polynomials (cf. (B.)), Ci := (µ +) M ( µ µ +) M l, i =, ( ) i(i+ µ µ +)(µ +) M i ( µ +M +3) i, i ; i!(µ 3 +) i ( µ µ +i+) M l + next put e l := e l. Step 3 For =,,...,M, o for =,,...,M, l =, +,...,M, l =,,...,M l, compute e +v l := ( [σ () σ (l)]e l σ ()e v l +σ (l)e l+v +σ (l)e l v )/σ (), where = (, ), l = (l,l ), v := (,), and where for t := (t,t ) we define σ (t) := ( t M)(t +µ +), σ (t) := t ( t µ 3 M ), σ (t) := σ (t)+σ (t), next put e l +v := e +v l ; o for l = +, +,...,M, l =,,...,M l, compute e +v l := ( [τ () τ (l)]e l τ ()e v l +τ (l)e l+v +τ (l)e l v )/τ (), 7
8 where = (,), l = (l,l ), v := (,), and for t := (t,t ) the coefficients τ j (t) are given by τ (t) := ( t M)(t +µ +), τ (t) := t ( t µ 3 M ), τ (t) := τ (t)+τ (t); next put e l +v := e +v l. Step 4 For, l Ω c m, compute where c := (c,c ), and E l (α,c,m) := U V V l e c l c, (3.) 3 U := ( α +3) c (α i +) c i, V h := i= ( m c h c )( m h As noticed in Remar B., the sum in (3.) can be evaluated efficiently using the Clenshaw s algorithm, at the cost of O(M l ) operations. ). 3. Computing the integrals I j The most computationally expensive part of the proposed method is the evaluation of the collection of integrals (.8). For example, for n + m =, if c = (,,), there are 76 two-dimensional integrals to be computed. It is obvious that using any standard quadrature would completely ruin the efficiency of the algorithm. Moreover, if any of the parameters α i (i =,,3) in (.6) is smaller than and the corresponding constrain parameter c i equals zero, then the integrands in (.8) are singular functions, and standard quadratures may fail to deliver any approximations to the integrals. Therefore, for evaluating the complete set of integrals (.8), we introduce a special scheme which is based on the general method [9] for approximating singular integrals. The proposed numerical quadrature is of the automatic type, which means that the required number of nodes is adaptively selected, depending on the complexity of the rational Bézier function, so that the requested accuracy of the approximation is always achieved. Most importantly, the algorithm is extremely effective in the considered application. In the example given at the beginning of this subsection (n+m = ), the time required to compute the whole collection of 76 integrals is only twice longer than the time needed to approximate a single separate integral of a similar type. First, we shall write the integral (.8) in a different form which is better suited for fast numerical evaluation. Observe that bivariate Bernstein polynomials (.3) can be expressed in terms of univariate Bernstein polynomials. Namely, we have Bj N (x) = BN j (x )B N j j (x /( x )), j := (j,j ), x := (x,x ), where Bi M (t) := ( M) i t i ( t) M i, i M, are univariate Bernstein polynomials. Further, the bivariate weight function w α (see (.6)) can be expressed as w α (x) = A α v α +α 3,α (x )v α3,α (x /( x )), Based on the Maple implementation of the algorithm. If the collection consists of 99 integrals (n+m = 4), the computation time increases by only 5% (compared to the case of 76 integrals). The detailed report from the efficiency test can be found at the end of Appendix B. 8
9 where v α,β (t) := ( t) α t β is the univariate Jacobi weight function. Hence, the integral (.8) can be written as x I j = w α (x) BN j (x) ω(x) dx dx ( = A α v α +α 3 +,α (s)bj N (s) = A α ( N j ) where we denoted N := n+m, and ( v a,b (t) ) v α3,α (t) BN j j (t) ω (s,t) dt ds ) v c,d (s) ω (s,t) ds dt, (3.3) a a(j) := α 3 +N j, b b(j ) := α +j, c c(j ) := α +α 3 +N j +, d d(j ) := α +j, ω (s,t) := ω(s,( s)t) = i= j= } (3.4) n n i w i (t)bi n (s), w i (t) = ω i,j Bj n i (t). (3.5) Note that the computation of values of the integrand is now much more effective, because the coefficients w i of the function ω ( i n) in (3.5) do not depend of the inner integration variable s. The main idea is, however, to compute the values of ω only once (at a properly selected set of quadrature nodes), and obtain a tool for fast computation of the integrals (3.3) for different values of a, b, c, and d, i.e. for different values of j. For arbitrary fixed t [,], define the function It is easy to see that we can write with where we use the notation ψ t (s) := ω (s,t). (3.6) I j = A α ( N j J(α,β,f) := ) J(a,b,Φ), Φ(t) := J(c,d,ψ t ), (3.7) ( x) α x β f(x)dx. The functions ψ t and Φ are analytic in a closed complex region containing the interval [,] (it is proved in Appendix B). This implies that (cf. [6, Chapter 3]) they can be accurately and efficiently approximated by polynomials given in terms of the (shifted) Chebyshev polynomials of the first ind, M t ψ t (s) S Mt (s) := γ [t] i T i (s ), Φ(t) ŜM(t) := i= M ˆγ l T i (t ), l= s,t, (3.8) 9
10 where M may depend on j, and the prime denotes a sum with the first term halved. Once the above expansions are computed (this can be done in a time proportional to M t log(m t ) and M log(m)), the integrals J(,, ) can be easily evaluated using the following algorithm that was proved in [4]. Algorithm 3. (Computing the integral J(α,β;S), S being a polynomial) Given numbers α, β >, let r := β α, u := α+β+. Let S M be a polynomial defined by S M (x) = M γ i T i (x ). i= Compute the sequence d i, i M+, by d M+ = d M :=, d i := rd i +(i u)d i+ γ i, i = M,M,...,. i+u Output: J(α,β;S M ) = B ( γ rd +ud ), where B := Γ(α+)Γ(β +)/Γ(α+β +). By the repeated use of the above very fast scheme, we may efficiently approximate the whole set of integrals I j for j Ω c n+m. The remaining technical details of the adaptive implementation of the proposed quadrature and the complete formulation of the integration algorithm are presented in Appendix A. 3.3 Main algorithm The method presented in this paper is summarized in the following algorithm. Algorithm 3.3 (Polynomial approximation of the rational Bézier triangular surface) Given the coefficients r and positive weights ω, Θ n, of the rational function (.3), the coefficients p of the degree m polynomial (.4), minimising the error (.5), with the constraints (.6), can be computed in the following way. Step Compute the table {El (α,c,m)},l Ω c by Algorithm 3.. m Step Compute the table {I j } j Ω c n+m by Algorithm A.. Step 3 For Γ c m, put p := g. Step 4 For Ω c m, compute p by (.7). Output: Set of the coefficients p, Θ m. 4 Examples In this section, we present some examples of approximation of rational triangular Bézier patches by triangular Bézier patches. No theoretical justification is nown for the best choice of the vector parameter α in the distance functional (.) if we use the error function (x) := R n (x) P m (x) (4.)
11 to measure the quality of the approximation. On the base of numerical experiments, we claim that α = (,, ) usually leads to slightly better results than the ones obtained for other "natural" choices of parameters, including the usually preferred α = (,, ) (meaning w α (x) = ). The computations were performed in 6-decimal-digit arithmetic. In the implementation of Algorithm A., we have assumed ε = 5 6 in (A.3), and used the initial values M = M = Example Let R 6 be the degree 6 rational triangular Bézier patch [8, Example ], R 6 (x) := ω r B 6 (x) Θ 6 ω B 6 (x), x T, (4.) Θ 6 T being the standard triangle (.), and the control points r and the associated weights ω being listed in Table. We let α = (,, ), c = (,,), and constructed the degree 5 Table : Control points r (upper entries) and the associated weights ω (lower entries) of the surface (4.), with = (, ) Θ 6. \ (6,,) (5,,3) (4,.5,3.5) (3,.,4) (.5,.5,) (.4,.4,) (,,) (5.,,3) (4.5,,3) (3,.6,4) (,.9,3) (.,,) (.4,.8,.6) (4.5,,5) (4,.,4) (3,,3) (,.,) (.8,.5,.5) (4,3,6) (.5,.5,5) (.5,.8,4) (,,3) (3.5, 3.5, 4) (.5, 3, 5) (.5, 3.5, 3) (3,4.,) (,4,) (,5,) best approximating polynomial patch P 5 (x) := p B 5 (x), x T, Θ 5 under the restriction p = g for Γ c 5, where Γ c 5 := { = (, ) : =, or =, or = 5}, and the set of points g, Γ c 5, is obtained in the following way. As well nown, the boundary of the patch (4.) is formed by three degree 6 rational Bézier curves. The least squares degree
12 5 polynomial approximation to each of these rational curves, with the endpoints preservation, is constructed using an extension of the method of [4], described in [3] (the input data: m = 5, α = β =, = l =, notation used being that of [3]). Now, the set of points g is the appropriate collection of all control points of the three resulting Bézier curves. We have repeated the computations for α = (,,) (with α = β =, in [3]), obtaining slightly worse results. The maximum errors max x T (x) (cf. (4.)) of the obtained results (see Fig. ) are about 5% less than those reported in [8, Table ] Figure : Constrained degree 5 polynomial approximation of the degree 6 rational triangular Bézier surface, with c = (,,). Upper part: Rational surface R 6 (x) and the approximating surface P 5 (x) with α = (,, ). Lower part: The error (x) plots corresponding to α = (,,) and α = (,, ), respectively. The maximum errors are.6 and.3, respectively. Notice that the original surface and the approximating surface agree at the corner points. 4. Example Let R be the composite rational surface, R R R 5 (y), y := ( x,x x ), x T R, (x) := R Y 5 (z), z := (x x, x ), x T Y, (4.3) where for C {R,Y}, R C 5 (w) := ω C rc B5 (w) Θ 5 ω C B5 (w), w T, (4.4) Θ 5
13 T being the standard triangle (.), and T R :={x = (x,x ) : x x, x }, T Y :={x = (x,x ) : x x, x }. The control points r C and the associated weights ωc of the rational patches (4.4) can be found at the webpage pwo/programs.html. The surface (4.3) is shown in Fig. 3 (the left plot). Now, we show how to obtain the degree m polynomial approximations of the rational subpatches, which form a C -continuous composite surface. o Let P Y m be the triangular Bézier patch of degree m approximating the rational patch R Y 5 without constraints, i.e., for c = (,,). Let py be the control points of the patch PY m. o We approximate the rational patch R R 5 by the triangular Bézier patch PR m of degree m, with constraints of the type c = (,,), where the points g Γ c m are chosen so that the C -continuity is obtained (cf. [6, Section 7]): g (,i) := p Y (i,), g (,i) := p Y (i+,) +(py (i+,) py (i,) ), i =,,...,m, i =,,...,m. The results, obtained for m = 5 and m = 6, with α = (,, ), are shown in Fig. 3. It can be observed that approximation of the rational composite surface (4.3) by two jointed polynomial patches of degree m = 5 (the middle plot) resulted in some visible differences. Increasing the degree of the approximating polynomials to m = 6 (the right plot) already gave a very satisfactory result Figure 3: The composite rational Bézier surface (4.3) (left) and the C -continuous composite polynomial surfaces of degree (5,5) (middle) and (6,6) (right). 3
14 Appendix A: The adaptive algorithm for computing the integrals I j We start with proving that the functions ψ t (3.6), t [,], and Φ (3.7) are analytic in a closed complex region containing the interval [,]. The assertion is clearly true in the case ofψ t (z) = ω (z,t), as the bivariate polynomial ω has no roots in [,] [,]. Similarly, for any s [,], the function z ω (s,z) is analytic in a rectangular region [ σ,+σ] [ σ,σ], where σ > does not depend on s. Thus, if s [, ], then ω (s,z) dz = C for any closed contour C [ σ,+σ] [ σ,σ]. Consequently, if α,β >, then ( ) ( ) ( s) α s β ω (s,z) ds dz = ( s) α s β ω (s,z) dz ds =. C C Therefore, by Morera s theorem (see, e.g., [, Chapter.3]), the function Φ(z) = J(α,β,ψ z ) is also analytic in [ σ, + σ] [ σ, σ]. The polynomials S Mt and ŜM in (3.8), which approximates the functions ψ t and Φ, are determined to satisfy the interpolation conditions S M (s j ) = ω (s j,t ), j M, M, Ŝ M (t ) = J(c,d;S M ), where, for simplicity, we denote M M t, and the interpolation nodes are given by s j = + In such a case, the coefficients γ [t ] i γ [t ] i = δ i,m M ˆγ l = δ l,m M M = cos jπ M, t = + and ˆγ l in (3.8) are given by M j= cos π M. ω (s j,t ) cos ijπ M, i M, J(c,d;S M )cos lπ M, l M, (A.) (A.) where δ j, is the Kronecer delta, the double prime means that the first and the last term of the sum are to be halved. The sets of coefficients (A.) can be very efficiently computed by means of the FFT with only O ( M log(m ) ) and O ( M log(m) ) arithmetic operations (cf. [7] or [4, Section 5.]; the authors recall that the FFT is not only fast, but also resistant to round-off errors). The presented approach is very convenient from the practical point of view because if the accuracy of the approximation (3.8) is not satisfactory, then we may double the value of M (or M) and reuse the previously computed results. The expansions (3.8) are accepted if M i=m 3 γ [t ] i max {, max i 3 γ[t ] i } 6ε and 4 M i=m 3 ˆγ i max {, max i 3 ˆγ i } 56ε, (A.3)
15 where ε is the computation precision. Here is the complete algorithm for efficient approximation of the whole set of integrals I j for j Ω c n+m. The functions (parameters) a, b, c, and d are defined in (3.4). Algorithm A. (Numerical computation of the set of integrals I j, j Ω c n+m) Let M := M, where M is an arbitrary integer greater than 7. Phase I. For {,,...,M} do the following Steps 6: Step. Computet according to (A.), and compute w i (t ) in (3.5) for i {,,...,n}. Step. Let M := M, where M is an arbitrary integer greater than 7. Step 3. Compute the values ω (s j,t ) for j {,,...,M }, where s j is given by (A.). Step 4. Using the FFT, compute the coefficients γ [t ] i ( i M ) defined in (A.). Step 5. If the first condition of (A.3) is not satisfied, then set M := M, compute the additional values ω (s j,t ) for j {,3,5,...,M }, and go to Step 4. Step 6. Compute the set of quantities W[t,j ] := J (c(j ),d(j );S M ), by applying Algorithm 3., for j {c,c +,...,N c c 3 }, where N = n+m. Phase II. For j {c,c +,...,N c 3 c } perform the following Steps 7 9: Step 7. Compute the coefficients ˆγ l ( l M) defined in (A.), by means of the FFT, using the stored values W[t,j ], M, in place of J ( c(j ),d(j );S M ). Step 8. If the second condition of (A.3) is not satisfied, then set M := M, and repeat Steps 6 for {,3,5,...,M }. Step 9. For j {c,c +,...,N c 3 j }, compute the integrals ( ) N ( ) I j I (j,j ) := A α J a(j),b(j );ŜM j using Algorithm 3.. Output: Set of the integrals I j for j Ω c n+m. Remar A. In Steps 4 and 7 of the above algorithm the coefficients γ [t ] i ( i M ) or ˆγ l ( l M) are recalculated each time the value of M or M is doubled. Such a procedure is advised if we use a system (lie, e.g., Maple or Matlab) equipped with a fast built-in FFT subroutine. If we are to program the FFT summation algorithm by ourselves, it should rather be done in such a way that practically all results computed for a previous value of M or M are reused (cf., e.g., [7]). In Table we present the results of the efficiency test, where the proposed quadrature (implemented in Maple) is compared to the Maple built-in integration subroutine. We have used the Bézier surface form Example 4. (n = 6), and set the parameters m and c to several different values, to obtain collections of integrals of different sizes (equal to Ω c n+m ). The experiment was performed in the 64-bit version of Maple 6 on the computer equipped with the 3.7GHz 5
16 Table : Comparison of the computation times of the Maple library function and the proposed adaptive quadrature (Algorithm A.) in the case of several collections of integrals (.8). The number of integrals which are to be computed equals Ω c n+m. Ω c computation time (in seconds) n+m Maple library function the proposed method FAILURE.89 i7 processor. All parameters α i in (.6) were set to (the efficiency of the proposed method does not depend on α, but the Maple built-in integration subroutine wors most efficiently with this selection). We have to eep in mind that Maple is an interpretative programming language with a pretty slow code interpreter. Therefore, the 4.7 times longer computation time of our quadrature, compared to the computation time of the Maple library function, in the case of -element collection of integrals is in fact an excellent result. The last collection of 99 integrals (n+m = 4) was too difficult to be computed by the Maple built-in subroutine (in 4-decimal digit arithmetic, assumed during this test). Appendix B: Hahn orthogonal polynomials The notation ( ) a,...,a r rf s z := b,...,b s = (a ) (a r )!(b ) (b s ) z is used for the generalized hypergeometric series (see, e.g., [,.]); here r, s Z +, z, a i,b j C, and (c) is the shifted factorial. The Hahn polynomials (see, e.g., [,.5]) h l (t) h l (t;a,b,m) := (a+) l ( M) l 3 F ( l,l+a+b+, t a+, M where l =,,...,M, a, b >, and M N, satisfy the recurrence relation ), (B.) h l+ (t) = A l (t,m)h l (t)+b l (M)h l (t), l ; h (t) ; h (t), (B.) with the coefficients A l (t,m) := C l (l +s ) t D l E l, B l (M) := D l E l, (B.3) where s := a+b+, C l := (l+s+)/[(l +s)(l+s )], D l := C l l(l+m +s)(l+b), and E l := (l+a+)(m l). Remar B. A linear combination of Hahn polynomials, s N (t) := N i= γ ih i (t;a,b,m), can be summed using the following Clenshaw s algorithm (see, e.g., [4, Thm 3..]). Compute the 6
17 sequence V,V,...,V n+ from V i := γ i +A i (t;m)v i+ +B i+ (M)V i+, i = N,N,...,, with V N+ = V N+ =, where the coefficients A i (t;m) and B i (M) are defined by (B.3). Then s N (t) = V. References [] L.V. Ahlfors, Complex Analysis, 3rd Ed., McGraw-Hill, 979. [] G.E. Andrews, R. Asey, R. Roy, Special Functions, Cambridge Univ. Press, Cambridge, 999. [3] J. Chen, G.J. Wang, Progressive-iterative approximation for triangular Bézier surfaces, Comp. Aided-Design 43 () [4] G. Dahlquist, A. Björc, Numerical Methods in Scientific Computing, vol. I, SIAM, Philadelphia, 8. [5] G. Farin, Triangular Bernstein-Bézier patches, Comput. Aided Geom. Design 3 (986) [6] G. Farin, Curves and Surfaces for Computer-Aided Geometric Design. A Practical Guide, 5th ed., Academic Press, Boston,. [7] W.M. Gentleman, Implementing Clenshaw-Curtis quadrature II. Computing the cosine transformation, Comm. ACM 5 (97) [8] Q.Q. Hu, An iterative algorithm for polynomial approximation of rational triangular Bézier surfaces, Appl. Math. Comput. 9 (3) [9] P. Keller, A method for indefinite integration of oscillatory and singular functions, Numer. Algor. 46 (7) 9 5. [] R. Koeoe, R.F. Swarttouw, The Asey scheme of hypergeometric orthogonal polynomials and its q-analogue, Rep. 98-7, Fac. Techn. Math. Informatics, Delft Univ. of Technology, Delft, 998. [] S. Lewanowicz, P. Keller, P. Woźny, Bézier form of dual bivariate Bernstein polynomials, arxiv:5.846 [math.na] (5). [] S. Lewanowicz, P. Woźny, Connections between two-variable Bernstein and Jacobi polynomials on the triangle, J. Comput. Appl. Math. 97 (6) [3] S. Lewanowicz, P. Woźny, P. Keller, Weighted polynomial approximation of rational Bézier curves, arxiv: [math.na] (5). [4] S. Lewanowicz, P. Woźny, P. Keller, Polynomial approximation of rational Bézier curves with constraints, Numer. Algor. 59 () [5] A. Rababah, Distances with rational triangular Bézier surfaces, Appl. Math. Comp. 6 (5) [6] T.J. Rivlin, Chebyshev Polynomials: From Approximation Theory to Algebra and Number Theory, nd ed., Wiley, New Yor, 99. [7] R. Sharma, Conversion of a rational polynomial in triangular Bernstein-Bézier form to polynomial in triangular Bernstein-Bézier form, Internat. J. Comput. Appl. Math. 8 (3) [8] P. Woźny, S. Lewanowicz, Constrained multi-degree reduction of triangular Bézier surfaces using dual Bernstein polynomials, J. Comput. Appl. Math. 35 () [9] P. Woźny, S. Lewanowicz, Multi-degree reduction of Bézier curves with constraints, using dual Bernstein basis polynomials, Comput. Aided Geom. Design 6 (9) [] H.X. Xu, G.J. Wang, Approximating rational triangular Bézier surfaces by polynomial triangular Bézier surfaces, J. Comput. Appl. Math. 8 (9) [] H.X. Xu, G.J. Wang, New algorithm of polynomial triangular B-B surfaces approximation to rational triangular B-B surfaces, J. Inform. Comput. Sci. 7 () [] L. Zhang, G.J. Wang, An effective algorithm to approximate rational triangular B-B surfaces using polynomial forms, Chinese J. Computers 9 (6) 5 6 (in Chinese). 7
Journal of Computational and Applied Mathematics
Journal of Computational and Applied Mathematics 233 2010) 1554 1561 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: wwwelseviercom/locate/cam
More informationWeighted G 1 -Multi-Degree Reduction of Bézier Curves
Vol. 7, No. 2, 216 Weighted G 1 -Multi-Degree Reduction of Bézier Curves Abedallah Rababah Department of Mathematics, Jordan University of Science and Technology Irbid 2211 Jordan Salisu Ibrahim Department
More informationConstrained Ultraspherical-Weighted Orthogonal Polynomials on Triangle
International Journal of Mathematical Analysis Vol. 9, 15, no., 61-7 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/1.1988/ijma.15.411339 Constrained Ultraspherical-Weighted Orthogonal Polynomials on Triangle
More informationarxiv:math/ v1 [math.na] 12 Jul 2004
arxiv:math/0407177v1 [math.na] 12 Jul 2004 On improving the accuracy of Horner s and Goertzel s algorithms Alica Smoktunowicz and Iwona Wróbel Faculty of Mathematics and Information Science, Warsaw University
More informationScientific Computing
2301678 Scientific Computing Chapter 2 Interpolation and Approximation Paisan Nakmahachalasint Paisan.N@chula.ac.th Chapter 2 Interpolation and Approximation p. 1/66 Contents 1. Polynomial interpolation
More informationGauss-Lobatto to Bernstein polynomials transformation
Gauss-Lobatto to Bernstein polynomials transformation Loredana Coluccio, Alfredo Eisinberg, Giuseppe Fedele Dip. Elettronica Informatica e Sistemistica, Università degli Studi della Calabria, 87036, Rende
More informationThe Essentials of CAGD
The Essentials of CAGD Chapter 4: Bézier Curves: Cubic and Beyond Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book www.farinhansford.com/books/essentials-cagd c 2000
More informationComparison of Clenshaw-Curtis and Gauss Quadrature
WDS'11 Proceedings of Contributed Papers, Part I, 67 71, 211. ISB 978-8-7378-184-2 MATFYZPRESS Comparison of Clenshaw-Curtis and Gauss Quadrature M. ovelinková Charles University, Faculty of Mathematics
More informationKey Words: cardinal B-spline, coefficients, moments, rectangular rule, interpolating quadratic spline, hat function, cubic B-spline.
M a t h e m a t i c a B a l a n i c a New Series Vol. 24, 2, Fasc.3-4 Splines in Numerical Integration Zlato Udovičić We gave a short review of several results which are related to the role of splines
More informationReading. w Foley, Section 11.2 Optional
Parametric Curves w Foley, Section.2 Optional Reading w Bartels, Beatty, and Barsky. An Introduction to Splines for use in Computer Graphics and Geometric Modeling, 987. w Farin. Curves and Surfaces for
More informationON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT
ON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT Received: 31 July, 2008 Accepted: 06 February, 2009 Communicated by: SIMON J SMITH Department of Mathematics and Statistics La Trobe University,
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More informationA matrix over a field F is a rectangular array of elements from F. The symbol
Chapter MATRICES Matrix arithmetic A matrix over a field F is a rectangular array of elements from F The symbol M m n (F ) denotes the collection of all m n matrices over F Matrices will usually be denoted
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationarxiv: v4 [math.na] 20 Oct 2016
Bézier form of dua bivariate Bernstein poynomias arxiv:151008246v4 [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
More informationA Reverse Technique for Lumping High Dimensional Model Representation Method
A Reverse Technique for Lumping High Dimensional Model Representation Method M. ALPER TUNGA Bahçeşehir University Department of Software Engineering Beşiktaş, 34349, İstanbul TURKEY TÜRKİYE) alper.tunga@bahcesehir.edu.tr
More informationarxiv: v1 [math.cv] 18 Aug 2015
arxiv:508.04376v [math.cv] 8 Aug 205 Saddle-point integration of C bump functions Steven G. Johnson, MIT Applied Mathematics Created November 23, 2006; updated August 9, 205 Abstract This technical note
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 informationarxiv: v1 [cs.lg] 26 Jul 2017
Updating Singular Value Decomposition for Rank One Matrix Perturbation Ratnik Gandhi, Amoli Rajgor School of Engineering & Applied Science, Ahmedabad University, Ahmedabad-380009, India arxiv:70708369v
More information8 A pseudo-spectral solution to the Stokes Problem
8 A pseudo-spectral solution to the Stokes Problem 8.1 The Method 8.1.1 Generalities We are interested in setting up a pseudo-spectral method for the following Stokes Problem u σu p = f in Ω u = 0 in Ω,
More informationNew series expansions of the Gauss hypergeometric function
New series expansions of the Gauss hypergeometric function José L. López and Nico. M. Temme 2 Departamento de Matemática e Informática, Universidad Pública de Navarra, 36-Pamplona, Spain. e-mail: jl.lopez@unavarra.es.
More informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationLectures 9-10: Polynomial and piecewise polynomial interpolation
Lectures 9-1: Polynomial and piecewise polynomial interpolation Let f be a function, which is only known at the nodes x 1, x,, x n, ie, all we know about the function f are its values y j = f(x j ), j
More informationResearch Article Approximate Implicitization Using Linear Algebra
Journal of Applied Mathematics Volume 2012, Article ID 293746, 25 pages doi:10.1155/2012/293746 Research Article Approximate Implicitization Using Linear Algebra Oliver J. D. Barrowclough and Tor Dokken
More informationSome approximation theorems in Math 522. I. Approximations of continuous functions by smooth functions
Some approximation theorems in Math 522 I. Approximations of continuous functions by smooth functions The goal is to approximate continuous functions vanishing at ± by smooth functions. In a previous homewor
More informationProbabilistic analysis of the asymmetric digital search trees
Int. J. Nonlinear Anal. Appl. 6 2015 No. 2, 161-173 ISSN: 2008-6822 electronic http://dx.doi.org/10.22075/ijnaa.2015.266 Probabilistic analysis of the asymmetric digital search trees R. Kazemi a,, M. Q.
More informationLinear Recurrence Relations for Sums of Products of Two Terms
Linear Recurrence Relations for Sums of Products of Two Terms Yan-Ping Mu College of Science, Tianjin University of Technology Tianjin 300384, P.R. China yanping.mu@gmail.com Submitted: Dec 27, 2010; Accepted:
More information1.1. The analytical denition. Denition. The Bernstein polynomials of degree n are dened analytically:
DEGREE REDUCTION OF BÉZIER CURVES DAVE MORGAN Abstract. This paper opens with a description of Bézier curves. Then, techniques for the degree reduction of Bézier curves, along with a discussion of error
More informationAsymptotics of Integrals of. Hermite Polynomials
Applied Mathematical Sciences, Vol. 4, 010, no. 61, 04-056 Asymptotics of Integrals of Hermite Polynomials R. B. Paris Division of Complex Systems University of Abertay Dundee Dundee DD1 1HG, UK R.Paris@abertay.ac.uk
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND First Printing, 99 Chapter LINEAR EQUATIONS Introduction to linear equations A linear equation in n unknowns x,
More informationOutline. 1 Interpolation. 2 Polynomial Interpolation. 3 Piecewise Polynomial Interpolation
Outline Interpolation 1 Interpolation 2 3 Michael T. Heath Scientific Computing 2 / 56 Interpolation Motivation Choosing Interpolant Existence and Uniqueness Basic interpolation problem: for given data
More informationPascal Eigenspaces and Invariant Sequences of the First or Second Kind
Pascal Eigenspaces and Invariant Sequences of the First or Second Kind I-Pyo Kim a,, Michael J Tsatsomeros b a Department of Mathematics Education, Daegu University, Gyeongbu, 38453, Republic of Korea
More informationGeneralized Akiyama-Tanigawa Algorithm for Hypersums of Powers of Integers
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 16 (2013, Article 1332 Generalized Aiyama-Tanigawa Algorithm for Hypersums of Powers of Integers José Luis Cereceda Distrito Telefónica, Edificio Este
More informationIB Mathematics HL Year 2 Unit 11: Completion of Algebra (Core Topic 1)
IB Mathematics HL Year Unit : Completion of Algebra (Core Topic ) Homewor for Unit Ex C:, 3, 4, 7; Ex D: 5, 8, 4; Ex E.: 4, 5, 9, 0, Ex E.3: (a), (b), 3, 7. Now consider these: Lesson 73 Sequences and
More informationa 11 x 1 + a 12 x a 1n x n = b 1 a 21 x 1 + a 22 x a 2n x n = b 2.
Chapter 1 LINEAR EQUATIONS 11 Introduction to linear equations A linear equation in n unknowns x 1, x,, x n is an equation of the form a 1 x 1 + a x + + a n x n = b, where a 1, a,, a n, b are given real
More informationTwo finite forms of Watson s quintuple product identity and matrix inversion
Two finite forms of Watson s uintuple product identity and matrix inversion X. Ma Department of Mathematics SuZhou University, SuZhou 215006, P.R.China Submitted: Jan 24, 2006; Accepted: May 27, 2006;
More informationGENERALIZED GAUSS RADAU AND GAUSS LOBATTO FORMULAE
BIT 0006-3835/98/3804-0101 $12.00 2000, Vol., No., pp. 1 14 c Swets & Zeitlinger GENERALIZED GAUSS RADAU AND GAUSS LOBATTO FORMULAE WALTER GAUTSCHI 1 1 Department of Computer Sciences, Purdue University,
More informationGenerating Functions for the q-bernstein Bases
Generating Functions for the q-bernstein Bases Ron Goldman Department of Computer Science, Rice University Houston, Texas 77251 USA Plamen Simeonov Department of Computer and Mathematical Sciences University
More information= i 0. a i q i. (1 aq i ).
SIEVED PARTITIO FUCTIOS AD Q-BIOMIAL COEFFICIETS Fran Garvan* and Dennis Stanton** Abstract. The q-binomial coefficient is a polynomial in q. Given an integer t and a residue class r modulo t, a sieved
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 informationarxiv: v1 [cs.it] 12 Jun 2016
New Permutation Trinomials From Niho Exponents over Finite Fields with Even Characteristic arxiv:606.03768v [cs.it] 2 Jun 206 Nian Li and Tor Helleseth Abstract In this paper, a class of permutation trinomials
More informationarxiv: v1 [math.na] 22 Jun 2010
On computing Bézier curves by Pascal matrix methods arxiv:1006.4327v1 [math.na] 22 Jun 2010 Licio Hernanes Bezerra a, Leonardo Koller Sacht b a Universidade Federal de Santa Catarina, Departamento de Matemática,
More informationInterpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials
Interpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials Lawrence A. Harris Abstract. We extend the definition of Geronimus nodes to include pairs of real numbers where
More informationMATHEMATICAL ENGINEERING TECHNICAL REPORTS
MATHEMATICAL ENGINEERING TECHNICAL REPORTS Combinatorial Relaxation Algorithm for the Entire Sequence of the Maximum Degree of Minors in Mixed Polynomial Matrices Shun SATO (Communicated by Taayasu MATSUO)
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 35, pp. 148-163, 2009. Copyright 2009,. ISSN 1068-9613. ETNA SPHERICAL QUADRATURE FORMULAS WITH EQUALLY SPACED NODES ON LATITUDINAL CIRCLES DANIELA
More informationSpectra of Multiplication Operators as a Numerical Tool. B. Vioreanu and V. Rokhlin Technical Report YALEU/DCS/TR-1443 March 3, 2011
We introduce a numerical procedure for the construction of interpolation and quadrature formulae on bounded convex regions in the plane. The construction is based on the behavior of spectra of certain
More informationA Combinatorial Interpretation of the Numbers 6 (2n)! /n! (n + 2)!
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 8 (2005), Article 05.2.3 A Combinatorial Interpretation of the Numbers 6 (2n)! /n! (n + 2)! Ira M. Gessel 1 and Guoce Xin Department of Mathematics Brandeis
More informationA SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 3, Fall 2009 A SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS SERGEY KHASHIN ABSTRACT. A new approach based on the use of new
More informationOn the decay of elements of inverse triangular Toeplitz matrix
On the decay of elements of inverse triangular Toeplitz matrix Neville Ford, D. V. Savostyanov, N. L. Zamarashkin August 03, 203 arxiv:308.0724v [math.na] 3 Aug 203 Abstract We consider half-infinite triangular
More informationMACMAHON S PARTITION ANALYSIS IX: k-gon PARTITIONS
MACMAHON S PARTITION ANALYSIS IX: -GON PARTITIONS GEORGE E. ANDREWS, PETER PAULE, AND AXEL RIESE Dedicated to George Szeeres on the occasion of his 90th birthday Abstract. MacMahon devoted a significant
More informationA gradient recovery method based on an oblique projection and boundary modification
ANZIAM J. 58 (CTAC2016) pp.c34 C45, 2017 C34 A gradient recovery method based on an oblique projection and boundary modification M. Ilyas 1 B. P. Lamichhane 2 M. H. Meylan 3 (Received 24 January 2017;
More informationarxiv: v1 [cs.sc] 29 Jul 2009
An Explicit Construction of Gauss-Jordan Elimination Matrix arxiv:09075038v1 [cssc] 29 Jul 2009 Yi Li Laboratory of Computer Reasoning and Trustworthy Computing, University of Electronic Science and Technology
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 informationSolutions of the Harmonic Oscillator Equation in a B-Polynomial Basis
Adv. Studies Theor. Phys., Vol. 3, 2009, no. 12, 451-460 Solutions of the Harmonic Oscillator Equation in a B-Polynomial Basis Muhammad I. Bhatti Department of Physics and Geology University of Texas-Pan
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 7 Interpolation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationMATH 167: APPLIED LINEAR ALGEBRA Chapter 3
MATH 167: APPLIED LINEAR ALGEBRA Chapter 3 Jesús De Loera, UC Davis February 18, 2012 Orthogonal Vectors and Subspaces (3.1). In real life vector spaces come with additional METRIC properties!! We have
More informationNumerical Solution of Nonlocal Parabolic Partial Differential Equation via Bernstein Polynomial Method
Punjab University Journal of Mathematics (ISSN 116-2526) Vol.48(1)(216) pp. 47-53 Numerical Solution of Nonlocal Parabolic Partial Differential Equation via Bernstein Polynomial Method Kobra Karimi Department
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 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 informationCOMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE
BIT 6-85//41-11 $16., Vol. 4, No. 1, pp. 11 118 c Swets & Zeitlinger COMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE WALTER GAUTSCHI Department of Computer Sciences,
More informationBrief Review of Vector Algebra
APPENDIX Brief Review of Vector Algebra A.0 Introduction Vector algebra is used extensively in computational mechanics. The student must thus understand the concepts associated with this subject. The current
More informationConstruction of Smooth Fractal Surfaces Using Hermite Fractal Interpolation Functions. P. Bouboulis, L. Dalla and M. Kostaki-Kosta
BULLETIN OF THE GREEK MATHEMATICAL SOCIETY Volume 54 27 (79 95) Construction of Smooth Fractal Surfaces Using Hermite Fractal Interpolation Functions P. Bouboulis L. Dalla and M. Kostaki-Kosta Received
More informationTridiagonal test matrices for eigenvalue computations: two-parameter extensions of the Clement matrix
Tridiagonal test matrices for eigenvalue computations: two-parameter extensions of the Clement matrix Roy Oste a,, Joris Van der Jeugt a a Department of Applied Mathematics, Computer Science and Statistics,
More informationNumerische Mathematik
Numer. Math. (1998) 80: 39 59 Numerische Mathematik c Springer-Verlag 1998 Electronic Edition Numerical integration over a disc. A new Gaussian quadrature formula Borislav Bojanov 1,, Guergana Petrova,
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K. R. MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Second Online Version, December 1998 Comments to the author at krm@maths.uq.edu.au Contents 1 LINEAR EQUATIONS
More informationNumerical Analysis: Interpolation Part 1
Numerical Analysis: Interpolation Part 1 Computer Science, Ben-Gurion University (slides based mostly on Prof. Ben-Shahar s notes) 2018/2019, Fall Semester BGU CS Interpolation (ver. 1.00) AY 2018/2019,
More informationOn cardinality of Pareto spectra
Electronic Journal of Linear Algebra Volume 22 Volume 22 (2011) Article 50 2011 On cardinality of Pareto spectra Alberto Seeger Jose Vicente-Perez Follow this and additional works at: http://repository.uwyo.edu/ela
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationMathematics for Engineers. Numerical mathematics
Mathematics for Engineers Numerical mathematics Integers Determine the largest representable integer with the intmax command. intmax ans = int32 2147483647 2147483647+1 ans = 2.1475e+09 Remark The set
More informationConstruction algorithms for plane nets in base b
Construction algorithms for plane nets in base b Gunther Leobacher, Friedrich Pillichshammer and Thomas Schell arxiv:1506.03201v1 [math.nt] 10 Jun 2015 Abstract The class of (0, m, s)-nets in base b has
More informationDirect method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
Journal of Computational and Applied Mathematics 22 (28) 51 57 wwwelseviercom/locate/cam Direct method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
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 informationMORE CONSEQUENCES OF CAUCHY S THEOREM
MOE CONSEQUENCES OF CAUCHY S THEOEM Contents. The Mean Value Property and the Maximum-Modulus Principle 2. Morera s Theorem and some applications 3 3. The Schwarz eflection Principle 6 We have stated Cauchy
More informationBivariate Lagrange interpolation at the Padua points: The generating curve approach
Journal of Approximation Theory 43 (6) 5 5 www.elsevier.com/locate/jat Bivariate Lagrange interpolation at the Padua points: The generating curve approach Len Bos a, Marco Caliari b, Stefano De Marchi
More informationSplit Algorithms and ZW-Factorization for Toeplitz and Toeplitz-plus-Hankel Matrices
Split Algorithms and ZW-Factorization for Toeplitz and Toeplitz-plus-Hanel Matrices Georg Heinig Dept of Math& CompSci Kuwait University POBox 5969 Safat 13060, Kuwait Karla Rost Faultät für Mathemati
More informationSome Formulas for the Principal Matrix pth Root
Int. J. Contemp. Math. Sciences Vol. 9 014 no. 3 141-15 HIKARI Ltd www.m-hiari.com http://dx.doi.org/10.1988/ijcms.014.4110 Some Formulas for the Principal Matrix pth Root R. Ben Taher Y. El Khatabi and
More informationThe generalized order-k Fibonacci Pell sequence by matrix methods
Journal of Computational and Applied Mathematics 09 (007) 33 45 wwwelseviercom/locate/cam The generalized order- Fibonacci Pell sequence by matrix methods Emrah Kilic Mathematics Department, TOBB University
More informationPolynomial Harmonic Decompositions
DOI: 10.1515/auom-2017-0002 An. Şt. Univ. Ovidius Constanţa Vol. 25(1),2017, 25 32 Polynomial Harmonic Decompositions Nicolae Anghel Abstract For real polynomials in two indeterminates a classical polynomial
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationChapter 11 - Sequences and Series
Calculus and Analytic Geometry II Chapter - Sequences and Series. Sequences Definition. A sequence is a list of numbers written in a definite order, We call a n the general term of the sequence. {a, a
More informationMIT Algebraic techniques and semidefinite optimization February 14, Lecture 3
MI 6.97 Algebraic techniques and semidefinite optimization February 4, 6 Lecture 3 Lecturer: Pablo A. Parrilo Scribe: Pablo A. Parrilo In this lecture, we will discuss one of the most important applications
More informationGenerating Functions of Partitions
CHAPTER B Generating Functions of Partitions For a complex sequence {α n n 0,, 2, }, its generating function with a complex variable q is defined by A(q) : α n q n α n [q n ] A(q). When the sequence has
More informationChebyshev finite difference method for a Two Point Boundary Value Problems with Applications to Chemical Reactor Theory
Iranian Journal of Mathematical Chemistry, Vol. 3, o., February 22, pp. 7 IJMC Chebyshev finite difference method for a Two Point Boundary Value Problems with Applications to Chemical Reactor Theory ABBAS
More informationResearch Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
Mathematical Problems in Engineering Volume 1, Article ID 693453, 1 pages doi:11155/1/693453 Research Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
More informationSMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract
SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm Abstract Dedicated to Prof. Michele Caputo on the occasion of his 8th birthday We consider ordinary fractional
More informationChapter 1 Computer Arithmetic
Numerical Analysis (Math 9372) 2017-2016 Chapter 1 Computer Arithmetic 1.1 Introduction Numerical analysis is a way to solve mathematical problems by special procedures which use arithmetic operations
More informationOn the Chebyshev quadrature rule of finite part integrals
Journal of Computational and Applied Mathematics 18 25) 147 159 www.elsevier.com/locate/cam On the Chebyshev quadrature rule of finite part integrals Philsu Kim a,,1, Soyoung Ahn b, U. Jin Choi b a Major
More informationarxiv: v4 [cs.gr] 23 Nov 2015
A de Casteljau Algorithm for Bernstein type Polynomials based on (p, q)-integers arxiv:507.040v4 [cs.gr 3 Nov 05 Khalid Khan, ) D.K. Lobiyal ) and Adem ilicman ) ) School of Computer and System Sciences,
More informationQuadrature Formula for Computed Tomography
Quadrature Formula for Computed Tomography orislav ojanov, Guergana Petrova August 13, 009 Abstract We give a bivariate analog of the Micchelli-Rivlin quadrature for computing the integral of a function
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationMacMahon s Partition Analysis VIII: Plane Partition Diamonds
MacMahon s Partition Analysis VIII: Plane Partition Diamonds George E. Andrews * Department of Mathematics The Pennsylvania State University University Park, PA 6802, USA E-mail: andrews@math.psu.edu Peter
More informationDepartment of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Spring Semester 2015
Department of Mathematics, University of California, Berkeley YOUR OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Spring Semester 205. Please write your - or 2-digit exam number on this
More informationBlock Lanczos Tridiagonalization of Complex Symmetric Matrices
Block Lanczos Tridiagonalization of Complex Symmetric Matrices Sanzheng Qiao, Guohong Liu, Wei Xu Department of Computing and Software, McMaster University, Hamilton, Ontario L8S 4L7 ABSTRACT The classic
More informationMINIMAL NORMAL AND COMMUTING COMPLETIONS
INTERNATIONAL JOURNAL OF INFORMATION AND SYSTEMS SCIENCES Volume 4, Number 1, Pages 5 59 c 8 Institute for Scientific Computing and Information MINIMAL NORMAL AND COMMUTING COMPLETIONS DAVID P KIMSEY AND
More informationAbsolute value equations
Linear Algebra and its Applications 419 (2006) 359 367 www.elsevier.com/locate/laa Absolute value equations O.L. Mangasarian, R.R. Meyer Computer Sciences Department, University of Wisconsin, 1210 West
More informationarxiv: v1 [cs.na] 8 Sep 2014
Degree reduction of disk rational Bézier curves Mao Shi a a College of Mathematics and Information Science of Shaanxi Normal University,Xi an 7162, China arxiv:149.2455v1 [cs.na] 8 Sep 214 Abstract This
More informationarxiv: v3 [math.ca] 16 Nov 2010
A note on parameter derivatives of classical orthogonal polynomials Rados law Szmytowsi arxiv:0901.639v3 [math.ca] 16 Nov 010 Atomic Physics Division, Department of Atomic Physics and Luminescence, Faculty
More information