Stanisław Lewanowicz, Paweł Woźny
Explicit formulae for the Bézier coefficients of the constrained dual Bernstein basis polynomials are derived in terms of the Hahn orthogonal polynomials. Using difference properties of the latter polynomials, efficient recursive scheme is obtained to compute these coefficients. Applications of this result to some problems of CAGD is discussed.
Filip Chudy, Paweł Woźny
A new algorithm for computing a point on a polynomial or rational curve in Bézier form is proposed. The method has a geometric interpretation and uses only convex combinations of control points. The new algorithm's computational complexity is linear with respect to the number of control points and its memory complexity is $O(1)$. Some remarks on similar methods for surfaces in rectangular and triangular Bézier form are also given.
Stanisław Lewanowicz, Paweł Keller, Paweł Woźny
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 thanks 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.
Stanisław Lewanowicz, Paweł Keller, Paweł Woźny
Dual Bernstein polynomials of one or two variables have proved to be very useful in obtaining Bézier form of the $L^2$-solution of the problem of best polynomial approximation of Bézier curve or surface. In this connection, the Bézier coefficients of dual Bernstein polynomials are to be evaluated at a reasonable cost. In this paper, a set of recurrence relations satisfied by the Bézier coefficients of dual bivariate Bernstein polynomials is derived and an efficient algorithm for evaluation of these coefficients is proposed. Applications of this result to some approximation problems of Computer Aided Geometric Design (CAGD) are discussed.
Stanisław Lewanowicz, Paweł Woźny, Paweł Keller
We present an efficient method to solve the problem of the constrained least squares approximation of the rational Bézier curve by the Bézier curve. The presented algorithm uses the dual constrained Bernstein basis polynomials, associated with the Jacobi scalar product, and exploits their recursive properties. Examples are given, showing the effectiveness of the algorithm.
Filip Chudy, Paweł Woźny
New differential-recurrence properties of dual Bernstein polynomials are given which follow from relations between dual Bernstein and orthogonal Hahn and Jacobi polynomials. Using these results, a fourth-order differential equation satisfied by dual Bernstein polynomials has been constructed. Also, a fourth-order recurrence relation for these polynomials has been obtained; this result may be useful in the efficient solution of some computational problems.
Paweł Woźny, Filip Chudy
It is well-known that a $d$-dimensional polynomial Bézier curve of degree $n$ can be subdivided into two segments using the famous de Casteljau algorithm in $O(dn^2)$ time. Can this problem be solved more efficiently? In this paper, we show that it is possible to do this in $O(dn\log{n})$ time using the fast Fourier transform and its inverse. Experiments show that the direct application of the new method performs well only for small values of $n$, as the algorithm is numerically unstable. However, a slightly modified version -- which still has $O(dn\log{n})$ computational complexity -- offers good numerical quality, which is confirmed by numerical experiments conducted in \textsf{Python}. Moreover, the new method has a nice property: if a Bézier curve is extended by an additional control point, the subdivision can be updated in $O(d)$ time. A similar idea can be applied to speed up the subdivision of rational Bézier curves and rectangular Bézier surfaces, as well as to compute the derivatives of Bézier curves more efficiently.
Filip Chudy, Paweł Woźny
We present new representations of Gauss--Legendre polynomials and their derivatives in the shifted power basis and in bases related to symmetric orthogonal Jacobi polynomials. Using these representations and certain recurrence relations, we propose efficient $O(n^2+dn)$ methods for evaluating a Gauss--Legendre curve of degree $n$ in $\mathbb E^d$. We also propose algorithms for multipoint evaluation with computational complexity $O(Mdn+dn^2)$, where $M$ is the number of evaluation points.
Przemysław Gospodarczyk, Paweł Woźny
In this paper, we present a new iterative approximate method of solving boundary value problems. The idea is to compute approximate polynomial solutions in the Bernstein form using least squares approximation combined with some properties of dual Bernstein polynomials which guarantee high efficiency of our approach. The method can deal with both linear and nonlinear differential equations. Moreover, not only second order differential equations can be solved but also higher order differential equations. Illustrative examples confirm the versatility of our method.
Przemysław Gospodarczyk, Paweł Woźny
In the paper, we give methods of construction of dual bases for the B-spline basis and truncated power basis. Explicit formulas for the dual B-spline basis are obtained using the Legendre-like orthogonal basis of the polynomial spline space presented in (Wei et al., Comput.-Aided Des. 45 (2013), 85-92) and a connection between orthogonal and dual bases of any space given in (Lewanowicz and Woźny, J. Approx. Theory 138 (2006), 129-150). Construction of the dual truncated power basis is performed in two phases. We start with explicit formulas for the dual power basis of the space of polynomials. Then, we expand this basis using an iterative algorithm proposed in (Woźny, J. Comput. Appl. Math. 260 (2014), 301-311). As a result, we obtain the dual truncated power basis. We also present some applications of the proposed dual polynomial spline bases and illustrative examples.
Przemysław Gospodarczyk, Paweł Woźny
In the paper, we show that the transformations between modified Jacobi and Bernstein bases of the constrained space of polynomials of degree at most $n$ can be performed with the complexity $O(n^2)$. As a result, the algorithm of degree reduction of Bézier curves that was first presented in (Bhrawy et al., J. Comput. Appl. Math. 302 (2016), 369--384), and then corrected in (Lu and Xiang, J. Comput. Appl. Math. 315 (2017), 65--69), can be significantly improved, since the necessary transformations are done in those papers with the complexity $O(n^3)$. The comparison of running times shows that our transformations are also faster in practice.
Paweł Woźny, Przemysław Gospodarczyk, Stanisław Lewanowicz
This paper deals with the merging problem of segments of a composite Bézier curve, with the endpoints continuity constraints. We present a novel method which is based on the idea of using constrained dual Bernstein polynomial basis (P. Woźny, S. Lewanowicz, Comput. Aided Geom. Design 26 (2009), 566--579) to compute the control points of the merged curve. Thanks to using fast schemes of evaluation of certain connections involving Bernstein and dual Bernstein polynomials, the complexity of our algorithm is significantly less than complexity of other merging methods.