[AstroPy] astropy.cosmology

Juande Santander Vela juandesant at gmail.com
Fri Dec 6 16:15:01 EST 2013


I could do that for a couple of them (Enrique's input might be useful here as use cases), but they should return Quantities, so I should familiarize myself with how to do that.

-- 
Juande Santander Vela

> On 6 Dec 2013, at 17:19, Erik Tollerud <erik.tollerud at gmail.com> wrote:
> 
> Right now I think Juande's solution is indeed the best way right now.
> I think the reason for only having the "forward" version of these is
> that it's not clear exactly which of all the cosmology functions
> should have inverses.
> 
> But we could add inverse versions of everything, basically implemented
> the way Juande is suggesting here.  Then they could later be optimized
> when someone wants to make them faster. If anyone wants to issue a PR
> doing that, I think it might be useful...
> 
>> On Thu, Dec 5, 2013 at 12:06 PM, Juande Santander Vela <jdsant at iaa.es> wrote:
>> Hi, Enrique,
>> 
>> I don’t find the age function in astropy.cosmology, but the loopback_time, but I think the answer stands.
>> 
>> If you want to obtain the inverse function, assuming loopback_time is continuous, and it should, you can use scipy.optimize methods to find out the redshift. For instance, the redshift corresponding to a 0.2Gyears lookback time, you can do:
>> 
>> from  scipy.optimize import newton
>> import astropy.cosmology as cosmo
>> 
>> newton(lambda x: cosmo.lookback_time(x).value-0.2, 0)
>> 
>> or, in general:
>> 
>> from  scipy.optimize import newton
>> import astropy.cosmology as cosmo
>> 
>> def redshift_from_lookback_time(gyears):
>>    return newton(lambda x: cosmo.lookback_time(x).value-gyears, 0)
>> 
>> Hope that helps.
>> 
>>> El 05/12/2013, a las 11:01, ES.Enrique Perez <eperez at iaa.es> escribió:
>>> 
>>> Geetings.
>>> 
>>> I may have overlooked the answer in the docs, but really haven't been able to find it:
>>> 
>>> astropy.cosmology age(some_redshift) provides the age corresponding to a given redshfit,
>>> but what is the reverse function that provides the redshift corresponding to a given age?
>>> 
>>> Thanks,
>>> Enrique
>>> _______________________________________________
>>> AstroPy mailing list
>>> AstroPy at scipy.org
>>> http://mail.scipy.org/mailman/listinfo/astropy
>> 
>> --
>> Juande Santander Vela
>> VIA-SKA Project Manager
>> Instituto de Astrofísica de Andalucía (IAA-CSIC)
>> Glorieta de la Astronomía s/n, E-18008, Granada, Spain
>> 
>> Henry Ford: La mayoría de la gente gasta más tiempo y energías en intentar esquivar los problemas que en tratar de resolverlos.
>> 
>> _______________________________________________
>> AstroPy mailing list
>> AstroPy at scipy.org
>> http://mail.scipy.org/mailman/listinfo/astropy
> 
> 
> 
> -- 
> Erik



More information about the AstroPy mailing list