bugGNU 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

Fri 27 Oct 2017 12:18:21 PM UTC, original submission:  

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:

$ astcosmiccal -z200
gsl: qng.c:189: ERROR: failed to reach tolerance with highest-order rule
Default GSL error handler invoked.
Aborted (core dumped)

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.

Mohammad Akhlaghi <makhlaghi>
Group administrator


(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

No files currently attached


Depends on the following items: None found

Items that depend on this one: None found


Carbon-Copy List
  • -email is unavailable- added by makhlaghi (Submitted the item)

    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

    Back to the top

    Powered by Savane 3.12