[gnuastro-devel] [task #14665] Cosmological functions in library

[Boud Roukema]

Follow-up Comment #1, task #14665 (project gnuastro):

See https://bitbucket.org/broukema/gnuastro/branch/cosmo-dev
for a proof-of-concept of a 100x speed up for Universe ages.
It's a bit hacked but enough to get started on 
improving the cosmo functions, I think.

I get the output from commit _dd4443a_ :

ageofuniverse benchmark: old 1.81705  new 0.013959, old/new 130.17
0.815234


I don't think we'll get much of a speed-up in proper distances
compared to https://arxiv.org/abs/1702.0886 (ignore the comments
there about luminosity distance being the main distance; it's
main in terms of observational analysis, but proper distance
is simpler in terms of the FLRW models.)

We would need an analytical integral of (sinh x)^(-2/3), because
d_{radial_comoving} = \int c dt/a(t) - see 
http://ned.ipac.caltech.edu/level5/March02/Sahni/Sahni3_2.html

The best that
https://en.wikipedia.org/wiki/List_of_integrals_of_hyperbolic_functions
suggests
would be to rewrite it in terms of \int (sinh x)^(4/3) , which
wouldn't help much. :P

We can implement what looks best for proper distances after re-reading
1702.0886.

The non-flat cases should be added too.


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/task/?14665>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/