In this paper, we propose the generalized B divided difference, with which the $(k−1)$th derivative of B-spline curves of order $k$ can be obtained directly without the need to compute the first $(k−2)$ derivatives as before. Based on the generalized B divided difference, the necessary and sufficient condition for degree-reducible B-spline curves is presented. Algorithms for degree reduction of B-spline curves are proposed using the constrained optimization methods.

Given a set of mixtures, blind source separation attempts to retrieve the source signals without or with very little information of the mixing process. We present a geometric approach for blind separation of nonnegative linear mixtures termed <i>facet component analysis</i>. The approach is based on facet identification of the underlying cone structure of the data. Earlier works focus on recovering the cone by locating its vertices (vertex component analysis) based on a mutual sparsity condition which requires each source signal to possess a stand-alone peak in its spectrum. We formulate alternative conditions so that enough data points fall on the facets of a cone instead of accumulating around the vertices. To find a regime of unique solvability, we make use of both geometric and density properties of the data points and develop an efficient facet identification method by combining data classification and linear

We study error propagation through implicit geometric problems by linearizing and estimating the linearization error. The method is particularly useful for quadratic constraints, which turns out to be no big restriction for many geometric problems in applications.

In this paper, two explicit conversion formulae between triangular and rectangular Bézier patches are derived. Using the formulae, one triangular Bézier patch of degree n can be converted into one rectangular Bézier patch of degree $n \times n$. And one rectangular Bézier patch of degree m × n can be converted into two triangular Bézier patches of degree $m + n$ . Besides, two stable recursive algorithms corresponding to the two conversion formulae are given. Using the algorithms, when converting triangular Bézier patches to rectangular Bézier patches, we can computer the relations between the control points of the two types of patches for any $n\ge2$ based on the relationships for $n=1$. When converting rectangular Bézier patches to triangular Bézier patches, we can computer the relations between the control points of the two types of patches for any $m\ge 2, n \subset N^+ $ and $n \ge 2, m \subset N^+$ based on the relationships for $m=n=1$.

Delaunay meshes (DM) are a special type of triangle mesh where the local Delaunay condition holds everywhere. We present an efficient algorithm to convert an arbitrary manifold triangle mesh M into a Delaunay mesh. We show that the constructed DM has O(Kn) vertices, where n is the number of vertices in M and K is a model-dependent constant. We also develop a novel algorithm to simplify Delaunay meshes, allowing a smooth choice of detail levels. Our methods are conceptually simple, theoretically sound and easy to implement. The DM construction algorithm also scales well due to its O(nK logK) time complexity. Delaunay meshes have many favorable geometric and numerical properties. For example, a DM has exactly the same geometry as the input mesh, and it can be encoded by any mesh data structure. Moreover, the empty geodesic circumcircle property implies that the commonly used cotangent Laplace-Beltrami operator has non-negative weights. Therefore, the existing digital geometry processing algorithms can benefit the numerical stability of DM without changing any codes. We observe that DMs can improve the accuracy of the heat method for computing geodesic distances. Also, popular parameterization techniques, such as discrete harmonic mapping, produce more stable results on the DMs than on the input meshes.