GNU Astronomy Utilities - Tasks: task #14689, faster ageofuniverse in cosmology...
You are not allowed to post comments on this tracker with your current authentication level.
task #14689: faster ageofuniverse in cosmology library + other cosmo extensions
Submitter: | Boud Roukema <boud> | ||
Submitted: | Tue 03 Oct 2017 02:53:24 PM UTC | ||
Should Start On: | Tue 03 Oct 2017 12:00:00 AM UTC | Should be Finished on: | Wed 03 Oct 2018 12:00:00 AM UTC |
Category: | Libraries | Priority: | 5 - Normal |
Item Group: | Enhancement | Status: | In Progress |
Privacy: | Public | Assigned to: | boud |
Percent Complete: | 20% | Open/Closed: | Open |
Effort: | 0.00 |
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.
Follows 1 latest change.
Date | Changed by | Updated Field | Previous Value | => | Replaced by |
---|---|---|---|---|---|
2017-10-18 | boud | Should be Finished on | 2017-10-03 | 2018-10-03 |
Powered by Savane 3.13-d3ae.
Corresponding source code
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 :
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.