Ejemplo n.º 1
0
    def test_normalize_by_int_spectrum(self):
        
        spectrum = Spectrum(x = numpy.linspace(3000.0, 10000.0, num=10000),
                            y = numpy.linspace(1e-13, 1e-11, num=10000))
        norm_spectrum = spectrum.normalize_by_int()

        astlib_spectrum = astsed.SED(wavelength=numpy.linspace(3000.0, 10000.0, num=10000),
                                     flux = numpy.linspace(1e-13, 1e-11, num=10000))
        astlib_spectrum.normalise()


        self.assertEqual(norm_spectrum.y[200], astlib_spectrum.flux[200])
        self.assertAlmostEqual(norm_spectrum.y[200] / norm_spectrum.norm_constant, spectrum.y[200])