Abstract: Given a set of data points lying on a smooth manifold, we present methods to interpolate those with piecewise Bézier splines. The spline is composed of Bézier curves (resp. surfaces) patched together such that the spline is continuous and differentiable at any point of its domain.