A neural multilevel method for high-dimensional parametric PDEs