The mixed directional difference–summation algorithm for generating the Bézier net of a trivariate four-direction Box-spline