search

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)
234 stars 289 forks source link

Beyond a_today. #163

Open rainwoodman opened 7 years ago

rainwoodman commented 7 years ago

Could you recommended a way to patch class allowing going beyond a_today?

Briefly looking at the code, currently the max of z_table is hard coded to a_today -- and thus it appears an additional parameter is necessary.

If we use class to guide a separate universe simulation, a_box may go beyond 1.0. However, (I think) we only want to specify the local cosmology parameter at a=1.0 rather than at the final a_box. (@eelregit could you comment if my understanding of the parameter scaling method is correct?)

This motivation might be wrong -- but going beyond a=1.0 may still be useful in general.

rainwoodman commented 7 years ago

Here is a patch that seem to have worked. It adds a_max to the parameters and replaces a_today with a_max in the integration loop of bg.

https://github.com/nickhand/classylss/pull/13/commits/36b690da99d339625e5b17e38cac9d317af7a760