We develop a novel technical approach to approximate vector fields on a surface. Given a set number of samples and some predefined vector fields as constraints, we could approximate all other vector fields on the surface (i.e. values at each face/triangle).