High-Dimensional Dynamic Stochastic Model Representation

We propose a generic and scalable method for computing global solutions of nonlinear, high-dimensional dynamic stochastic economic models. First, within an MPI–TBB parallel time-iteration framework, we approximate economic policy functions using an adaptive, high-dimensional model representation scheme, combined with adaptive sparse grids. With increasing dimensions, the number of points in this efficiently-chosen combination of low-dimensional grids grows much more slowly than standard tensor product grids, sparse grids, or even adaptive sparse grids. Moreover, the adaptivity within the individual component functions adds an additional layer of sparsity, since grid points are added only where they are most needed — that is to say, in regions of the computational domain with steep gradients or at non-differentiabilities. Second, we introduce a performant vectorization scheme of the interpolation compute kernel. Third, we validate our claims with numerical experiments conducted on “Piz Daint" (Cray XC50) at the Swiss National Super-computing Center. We observe significant speedups over the state-of-the-art techniques, and almost ideal strong scaling up to at least 1, 000 compute nodes. Fourth, to demonstrate the broad applicability of our method, we compute global solutions to two different versions of a dynamic stochastic economic model: a high-dimensional international real business cycle model with capital adjustment costs, and with or without irreversible investment. We solve these models up to 300 continuous state variables globally