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

Submitted by:  Mohammad Akhlaghi <makhlaghi>
Submitted on:  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>
Project Administrator


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)

    Do you think this task is very important?
    If so, you can add your encouragement to it.
    This task has 0 encouragements so far.

    Only logged-in users can vote.




    No changes have been made to this item

    Back to the top

    Powered by Savane 3.7