astrophysics

At the moment, when reading a FITS file with BUNIT set to adu with CCDData.read, one gets the following info message:

INFO:astropy:using the unit adu passed to the FITS reader instead of the unit adu in the FITS file.

it might make sense to avoid emitting this message if the units match? (cc @mwcraig)

We could look into using the diskcache package to cache MAST search operations:

https://pypi.org/project/diskcache/

The package is a disk-based key-value store which makes it easy to specify an expiration time for each entry. This would allow us to e.g. cache search queries for 24 hours, which would significantly speed up searches across different Python processes, whilst largely avoiding

The Poisson solve absolute tolerance sets the accuracy of the Poisson solve. Given the value of abs_tol, the error level targeted in the Poisson solve is abs_tol * 4 * pi * G * rho_max, where rho_max is the maximum density on the domain. For Cartesian simulations the default value is 1.e-11, and for non-Cartesian the default value is 1.e-10.

At high resolution, and/or for complicated mass distr

Today one student has tried to evolve a group of orbits using the MWPotential2014 and found that some bound orbits have positive energy obtained from .E() function, which looks incosistent. We have investigated the potential calculation of MWPotential2014 and found that the potential zero for PowerSphericalPotentialwCutoff is chosen at (R=0, z=0) while the potential zero for other two compon

Description
The Monte Carlo Configuration Parameter no_of_packets has currently some inappropriate description in the TARDIS Documentation.

To Reproduce
A Better description that tells the actual use of the input parameter.

🎯 Describe what you'd like

The possibility to plot any variable (radiance, absorbance...) on a double x-axis plot would be really convenient!

💡 Possible implementations

In spectrum.plot the argument wunit could be used.
For exam

The following seminar work of Hernquist details a light profile called the Hernquist profile:

https://ui.adsabs.harvard.edu/abs/1990ApJ...356..359H/abstract

The intensity calculation for a spherical model is given by equations (32-35):

It would be good to implement this light

As per comment here: https://github.com/StingraySoftware/stingray/pull/610/files#r782942772

It was noted in the review for our JOSS paper that the keyword argument noise for the function redshifts_from_comoving_density has a misleading name. See the review comment here. They suggested renaming to fixed_N. My own suggestion would be shot_noise (or perhaps poisson). Since this would require a break