@@ -0,0 +1,979 @@
+from .. import hyperfine%0Aimport numpy as np%0Afrom astropy import units as u%0Afrom astropy import constants%0A%0Adef test_hyperfine():%0A    # Regression test for issue 125%0A%0A    vtau_model = hyperfine.hyperfinemodel(%5B'one'%5D, %7B'one':0.0%7D,%0A                                          %7B'one':100e9%7D,%0A                                          %7B'one':1.0%7D,%0A                                          %7B'one':1.0%7D)%0A%0A%0A    xarr = %5B100%5D*u.GHz%0A    Tex = 10.0%0A    Tbackground = 2.73%0A    tau = 0.1%0A    result = vtau_model.hyperfine(xarr, Tex=Tex,%0A                                  tau=tau, xoff_v=0.0, width=1.0,%0A                                  Tbackground=Tbackground)%0A%0A    hoverk = (constants.h.cgs/constants.k_B.cgs).value%0A    T0 = hoverk * xarr.to(u.Hz).value%0A    oneminusetotheminustau = (1.0-np.exp(-np.array(tau)))%0A    term1 = (np.exp(T0/Tex) - 1)**-1%0A    term2 = (np.exp(T0/Tbackground) - 1)**-1%0A    correctresult = oneminusetotheminustau*T0*(term1 - term2)%0A%0A    assert result == correctresult%0A

add the test I left out