期刊
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING
卷 63, 期 3, 页码 455-471出版社
JOHN WILEY & SONS LTD
DOI: 10.1002/nme.1296
关键词
tricubic; interpolation; computational dynamics
The purpose of this paper is to give a local tricubic interpolation scheme in three dimensions that is both C-1 and isotropic. The algorithm is based on a specific 64 x 64 matrix that gives the relationship between the derivatives at the corners of the elements and the coefficients of the tricubic interpolant for this element. In contrast with global interpolation where the interpolated function usually depends on the whole data sit, our tricubic local interpolation only uses data in a neighbourhood of an element. We show that the resulting interpolated function and its three first derivatives are continuous if one uses cubic interpolants. The implementation of the interpolator can be downloaded as a static and dynamic library for most platforms. The major difference between this work and current local interpolation schemes is that we do not separate the problem into three one-dimensional problems. This allows for a much easier and accurate computation of higher derivatives of the extrapolated field. Applications to the computation of Lagrangian coherent structures in ocean data are briefly discussed. Copyright (c) 2005 John Wiley & Sons, Ltd.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据