I recently had a bug in my code that obviously was caused by an issue with floating point precision but had me scratching my head how it came about. Here it is:
In [1]:
import numpy as np
from astropy.table import Table
from astropy import cosmology
cosmo = cosmology.WMAP9 …