@@ -0,0 +1,751 @@
+import numpy as np%0Aimport os.path%0Aimport itertools%0A%0Afrom read_spectrum_fits import QSORecord%0Afrom hdf5_spectrum_container import Hdf5SpectrumContainer, Hdf5SpectrumIterator%0A%0AFILENAME = '/mnt/gastro/yishay/sdss_QSOs/spectra.hdf5'%0AMAX_WAVELENGTH_COUNT = 4992%0A%0A%0Adef return_spectra_2(qso_record_table, mem_map_file=FILENAME):%0A%0A    spectra = Hdf5SpectrumContainer(mem_map_file, True)%0A%0A    # we assume that the order of spectra is the same as in the QSO list%0A    for i, j in itertools.izip(Hdf5SpectrumIterator(spectra), qso_record_table):%0A        qso_rec = QSORecord.from_row(j)%0A        ar_wavelength = i.get_wavelength()%0A        ar_flux = i.get_flux()%0A        if ar_wavelength.size == 0:%0A            continue%0A        yield ar_wavelength, ar_flux, qso_rec%0A

add missing file