Finite Differences on Sparse Grids for Continuous-Time Heterogeneous Agent Models
We present a finite difference method working on sparse grids to solve higher dimensional heterogeneous agent models. If one wants to solve the arising Hamilton–Jacobi–Bellman equation on a standard full grid, one faces the problem that the number of grid points grows exponentially with the number o...
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
MDPI AG
2025-01-01
|
| Series: | Algorithms |
| Subjects: | |
| Online Access: | https://www.mdpi.com/1999-4893/18/1/40 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Summary: | We present a finite difference method working on sparse grids to solve higher dimensional heterogeneous agent models. If one wants to solve the arising Hamilton–Jacobi–Bellman equation on a standard full grid, one faces the problem that the number of grid points grows exponentially with the number of dimensions. Discretizations on sparse grids only involve <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi mathvariant="script">O</mi><mo>(</mo><mi>N</mi><msup><mrow><mo>(</mo><mo form="prefix">log</mo><mi>N</mi><mo>)</mo></mrow><mrow><mi>d</mi><mo>−</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow></semantics></math></inline-formula> degrees of freedom in comparison to the <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi mathvariant="script">O</mi><mo>(</mo><msup><mi>N</mi><mi>d</mi></msup><mo>)</mo></mrow></semantics></math></inline-formula> degrees of freedom of conventional methods, where <i>N</i> denotes the number of grid points in one coordinate direction and <i>d</i> is the dimension of the problem. While one can show convergence for the used finite difference method on full grids by using the theory introduced by Barles and Souganidis, we explain why one cannot simply use their results for sparse grids. Our numerical studies show that our method converges to the full grid solution for a two-dimensional model. We analyze the convergence behavior for higher dimensional models and experiment with different sparse grid adaptivity types. |
|---|---|
| ISSN: | 1999-4893 |