def multi_init(n, Nexp, NMC, **kwargs):
energy = kwargs.pop("energy", False)
llh = MultiPointSourceLLH(nsource=25,
nsource_bounds=(-Nexp / 2.,
Nexp / 2.) if not energy else
(0., Nexp / 2.),
seed=np.random.randint(2**32),
**kwargs)
for i in xrange(n):
llh.add_sample(str(i), init(Nexp, NMC, energy=energy))
return llh
