limited by the number of galaxies
We need a plan to deal with large numbers of galaxies
The pipeline slows at Sim Gal Spec Clean when creating enough memory to put all the reconstructed/simulated spectra into :
Ngal = len(id_fiber)
Nwave = len(tspec_v0[0,:])
flux = np.zeros((Nwave, Ngal)) #!!! The program takes a long time just to instantiate this 2d array
first fix: reduce the number of galaxies for testing purposes
second fix: append to empty array, instead of instantiate full array at beginning; this will probably only kick the can down the road.
longer-term fix: how do we do this for millions of galaxies by Friday?