Density estimation with adaptive sparse grids for large data sets
Algorithm 847: Spinterp: piecewise multilinear hierarchical sparse grid interpolation in MATLAB