arXiv:2405.15735v3 Announce Type: replace Abstract: We introduce a numerical method to approximate arbitrary differential operators on vector fields in weak form, given point cloud data sampled randomly from a $d$-dimensional manifold embedded in $\mathbb{R}^n$. This method generalizes the local linear mesh method to the local curved mesh method, thus allowing for the estimation of differential operators with nontrivial Christoffel symbols, such as the Bochner and Hodge Laplacians. In particular, we leverage the potentially small intrinsic dimension of the manifold $(d \ll n)$ to construct local parameterizations that incorporate both local meshes and higher-order curvature information. The former is constructed using low-dimensional meshes obtained from local data projected onto the tangent spaces, while the latter is obtained by fitting local polynomials using generalized moving least squares. Theoretically, we prove spectral convergence of the proposed method for the estimation of the Bochner Laplacian. We provide numerical results supporting the theoretical convergence rates for the Bochner and Hodge Laplacians on simple manifolds.