We use Bayesian GPR to model the relationship between the crystal structures and the response: lattice energies evaluated at multiple levels of accuracy and corresponding computational complexity. Each structure is uniquely identified by the vector of the structural descriptor values x ∈ R s , concatenated symmetry functions as described in Section 2.3. In a system with SF r radial and SF a angular symmetry functions, the dimensionality of structures with