4th order diffusion tensor interpolationwith divergence and curl constrained Bézier patches

We present a tensor field interpolation method based on tensor-valued Bezier patches. The control points of the patch are determined by imposing physical constraints on the interpolated field by constraining the divergence and curl of the tensor field. The method generalizes to Cartesian tensors of all orders. Solving for the control points requires the solution of a sparse linear system. Results are presented for order 2 and 4 tensors and comparisons are made with linear interpolation and a previously proposed subdivision scheme.