Approximation of Solution Operators for High-dimensional PDEs