def test_planck(self):
import astropy.units as u
import astropy.constants as const
wave = np.logspace(-1, 3, 10000) * u.um
I = util.planck(wave, 300, unit=u.Unit('W/(m2 um sr)'))
assert allclose(wave[I.value.argmax()], 2.8977685e3 / 300, rtol=1e-3)
F = (util.davint(wave.value, I.value, wave.value[0], wave.value[-1])
* np.pi)
assert allclose(F, const.sigma_sb.si.value * 300.**4, rtol=1e-5)
def test_davint(self):
x = np.linspace(0, 2 * pi)
y = np.sin(x)
assert allclose(util.davint(x, y, 0, 2 * pi), 0)
