Efficiently transforming from values of a function on a sparse grid to basis coefficients
In many contexts it is necessary to determine coefficients of a basis expansion of a function f(x_1, …, x_D) from values of the function at points on a sparse grid. Knowing the coefficients, one has an interpolant or a surrogate. For example, such coefficients are used in uncertainty quantification. In this chapter, we present an efficient method for computing the coefficients. It uses basis functions that, like the familiar piecewise linear hierarchical functions, are zero at points in previous levels. They are linear combinations of any, e.g. global, nested basis functions φ_i_k^(k)(x_k). Most importantly, the transformation from function values to basis coefficients is done, exploiting the nesting, by evaluating sums sequentially. When the number of functions in level ℓ_k equals ℓ_k (i.e. when the level index is increased by one, only one point (function) is added) and the basis function indices satisfy ‖𝐢-1‖_1 ≤ b, the cost of the transformation scales as 𝒪(D [b/D+1 + 1] N_sparse), where N_sparse is the number of points on the sparse grid. We compare the cost of doing the transformation with sequential sums to the cost of other methods in the literature.
READ FULL TEXT