def testGetTable(self): """ model.get_table() single case """ wave, flux = model.get_table(teff=6874, logg=4.21, ebv=0.0) self.assertEqual(len(wave), len(flux)) self.assertEqual(len(wave), 1221) self.assertAlmostEqual(wave[0], 90.9, delta=0.1) self.assertAlmostEqual(wave[-1], 1600000.0, delta=0.1) self.assertAlmostEqual(wave[400], 4370.0, delta=0.1) self.assertAlmostEqual(wave[800], 15925.0, delta=0.1) self.assertAlmostEqual(flux[0], 0.0, delta=0.0001) self.assertAlmostEqual(flux[-1], 0.020199, delta=0.0001) self.assertAlmostEqual(flux[400], 24828659.5845, delta=0.0001) self.assertAlmostEqual(flux[800], 1435461.60457, delta=0.0001)
def testGetTableBinary(self): """ model.get_table() multiple case """ wave,flux = model.get_table(teff=25000, logg=5.12, ebv=0.001, teff2=33240, logg2=5.86, ebv2=0.001) self.assertEqual(len(wave), len(flux)) self.assertEqual(len(wave), 123104) self.assertAlmostEqual(wave[0], 1000, delta=0.1) self.assertAlmostEqual(wave[-1], 24999.8, delta=0.1) self.assertAlmostEqual(wave[40000], 8671.20117188, delta=0.001) self.assertAlmostEqual(wave[80000], 16509.4003906, delta=0.001) self.assertAlmostEqual(flux[0], 170415154318.0, delta=0.5) self.assertAlmostEqual(flux[-1], 2457826.26898, delta=0.0001) self.assertAlmostEqual(flux[40000], 141915936.111, delta=0.001) self.assertAlmostEqual(flux[80000], 12450102.801, delta=0.001)