Upcoming Events

Using Adaptive Sparse Grids to Solve High-Dimensional Dynamic Models

November 6, 2013 3:00PM to 4:00PM
Simon Scheidegger, University of Zurich
Building 240, Room 4301
LANS Informal Seminar
We present a flexible and scalable method to compute global solutions of high-dimensional non-smooth dynamic models. Within a time-iteration setup, we interpolate policy functions using an adaptive sparse grid algorithm with piecewise multi-linear (hierarchical) basis functions. As the dimensionality increases, sparse grids grow considerably slower than standard tensor product grids. In addition, the grid scheme we use is automatically refined locally and can thus captures steep gradients or even kinks.

To further increase the maximal problem size we can handle, our implementation is fully hybrid parallel, i.e. using a combination of MPI and OpenMP. This parallelization enables us to efficiently use modern high-performance computing architectures. Our time iteration algorithm scales up nicely to at least 1,200 parallel processes. To demonstrate the performance of our method, we apply it to an international real business cycle model with capital adjustment costs, irreversible investment, and more than twenty continuous state variables.