Multiple StellarPopulation instances

Because parameter state and dirtiness is specific to each instance, but the computed model is held by the module level driver object, having multiple StellarPopulation instances can cause a model to not be recomputed even though it should be (example code below) This might also cause problems with shared memory multiprocessing.

Some documentation advising only a single sps instance per program might be a good idea.

from numpy.testing import assert_allclose
import fsps
sps = fsps.StellarPopulation(zcontinuous=1, vactoair_flag=False)
pop = fsps.StellarPopulation(zcontinuous=1, vactoair_flag=False)
w1, s1 = sps.get_spectrum(tage=10.0)
w, s = pop.get_spectrum(tage=1.0)
w2, s2 = sps.get_spectrum(tage=10.0)
assert_allclose(s, s2) # this should not pass, but does
assert_allclose(s1, s2) # this should not fail, but does

Inizia con StellarPopulation.get_spectrum e con lo stato del driver a livello di modulo descritto nell’issue. Riproduci la sequenza fornita con due istanze, quindi traccia il modo in cui vengono condivisi lo stato di modifica dei parametri e il ricalcolo del modello. Il lavoro è completo quando ogni istanza restituisce risultati coerenti con i propri parametri, incluse le due assertions nell’esempio.