GNU Astronomy Utilities - Bugs: bug #52295, Cosmology library integrals crash...
You are not allowed to post comments on this tracker with your current authentication level.
bug #52295: Cosmology library integrals crash for high z
Submitter: | Mohammad Akhlaghi <makhlaghi> | ||
Submitted: | Fri 27 Oct 2017 12:18:21 PM UTC | ||
Category: | Libraries | Severity: | 3 - Normal |
Item Group: | Crash | Status: | Confirmed |
Privacy: | Public | Assigned to: | None |
Open/Closed: | Open |
No files currently attached
Depends on the following items: None found
Items that depend on this one: None found
Carbon-Copy List
There are 0 votes so far. Votes easily highlight which items people would like to see resolved in priority, independently of the priority of the item set by tracker managers.
No changes have been made to this item
The cosmology library functions use GSL's qng integration function for their integrations. However, when you run them with a high redshift (for example 200, as shown in the example below with the CosmicCalculator program which uses those library functions), they crash:
As the error suggests, this may be fixed by decreasing the tollerance level, but the question is this: how much can the tolerance level be decreased to give a reasonable result at all redshifts? Or maybe, should we set a varying tolerance level based on redshift?
As an observer, I hardly ever need any calculation above redshift of ~10, but I guess simulations do need it. So another question is this: should we set a limit on the highest redshift (maybe until recombination at z=1100)?
Please share your thoughts on the best way forward.