[bug #52295] Cosmology library integrals crash for high z

[Boud Roukema]

Follow-up Comment #6, bug#52295 (group gnuastro):

Commit d74dbcbe on branch bug52295 looks fine to me to handle the higher
redshifts:

 for z in 0.01 0.1 1 10 20 100 200 1000 1200 2000; do
./bin/cosmiccal/astcosmiccal --quiet --properdistance --oradiation=0.0
--omatter=0.3 --olambda=0.7 --H0=70 --redshift=$z; done |& grep "^[0-9]"
42.730926
418.454488
3303.828806
9440.249626
10742.296524
12598.733839
13051.764583
13660.529297
13703.558079
13805.213886


As a side comment, changing from qng to qag doesn't hurt, but there are more
significant issues for calculational speed. For the extreme slowdown due to
shell overhead, see task 14691: https://savannah.gnu.org/task/?14691 ;
cosmdist is 50000 or so times faster, which will be non-negligible for
LSST/Rubin photo-z conversions.



    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?52295>

_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/