lesgourg / class_public

Public repository of the Cosmic Linear Anisotropy Solving System (master for the most recent version of the standard code; GW_CLASS to include Cosmic Gravitational Wave Background anisotropies; classnet branch for acceleration with neutral networks; ExoCLASS branch for exotic energy injection; class_matter branch for FFTlog)
230 stars 285 forks source link

Accuracy settings for perturbations #398

Open DanielBThomas opened 3 years ago

DanielBThomas commented 3 years ago

Hi

I've been trying to understand which accuracy settings to adjust to get sub% accuracy on the perturbations at late times (particularly the density perturbations and metric potentials).

It seems that not many of the accuracy settings affect these; the primary ones I've found seem to be evolver (=0 seems to be better) and tol_perturb_integration (smaller is better).

In addition, perturb_integration_stepsize seems to be affecting the output without really making it any more accurate (as measured by checking the internal consistency of the output perturbations, e.g. comparing both sides of the poisson equation).

I guess my questions are: Have I missed any important accuracy parameters for late time perturbation evolution? Has anyone else looked at this and could they recommend settings for this?

Thanks Dan