Constraining cosmological parameters from N-body simulations with Bayesian Neural Networks