## Tabulation of thin plate splines on a very fine two-dimensional grid

##### View/Open

##### Author

##### Date

1992##### Permanent Link

http://hdl.handle.net/10092/11618A thin plate spline approximation has the form [FORMULA] where {⋋j ∈ R, : j = 1,2, ... , n} and {xj ∈ R² : j = 1, 2, ... , n} are parameters and where p is a linear polynomial. There exist several applications that require s to be tabulated at all the lattice points of a very fine square grid. For example, 10⁸ grid points and n = 500 can occur, and then the direct evaluation of s at every grid point would be impracticable. Fortunately each thin plate spline term is smooth away from its centre xj, so it is possible to apply a scheme that subtabulates by finite differences provided that special attention is given to those terms whose centres are close to the current x. Thus the total work is bounded by a small constant multiple of the number of grid points plus a constant multiple of n ∊ -¹ /³llog hl, where ∊ is a given tolerance on the calculated values of s(x) and where h is the mesh size of the fine grid. We will find that the exponent -1/3 is due to the order of the differences that are employed. An algorithm for this calculation is described and discussed and some numerical results are presented. The errors of the subtabulation procedures are studied in an appendix.

##### Subjects

Field of Research::01 - Mathematical Sciences::0101 - Pure Mathematics::010101 - Algebra and Number Theory##### Collections

- Engineering: Reports [684]