using grid_type = Grid< T, N >using array_type = blitz::Array< T, N >using function_type = Linear_function< T, N >Linear_interpolation(array_type points, const grid_type & grid)explicitNearest neighbour weighted interpolation for 3-d fields.
using grid_type = Grid< T, 3 >using index_type = blitz::TinyVector< int, 3 >using vertex_field = blitz::Array< return_type, 3 >using vec3 = blitz::TinyVector< T, 3 >using return_type = Returnoperator()(const vec3 & x, const index_type & i00) -> return_typeLinear_function_irregular(const vertex_field & points, const grid_type & grid)explicitusing return_type = Returnusing function_type = Functionusing index_type = blitz::TinyVector< int, N >using grid_type = Grid< T, N >using vertex_type = blitz::TinyVector< T, N >