
import os
import smooth1d




#(0) Simulate one dataset, multuple parameters:
nIterations = 3
JJ          = [5, 6, 7, 8, 9, 10, 15, 20, 25, 30, 35, 45, 50] #sample sizes
noise_amps  = [0.01, 0.20]
sim         = smooth1d.sim.Simulator()
sim.set_alpha(0.05)
sim.set_filter('None')
# sim.set_dataset_name('Vaughan1982')
# sim.set_dataset_name('Challis1999a')
# sim.set_dataset_name('Challis1999b')
# sim.set_dataset_name('Challis1999c')
# sim.set_dataset_name('Challis1999d')
sim.set_dataset_name('Challis1999e')
sim.set_metadata_labels(['sample_size','noise_amp'], types=[int,float])
sim.set_results_directory( os.path.dirname(__file__) )
sim.set_seed()
### simulate:
for J in JJ:
	for amp in noise_amps:
		print('Sample size = %d, Noise amplitude = %.2f' %(J,amp))
		sim.set_sample_size( J )
		sim.set_noise_amp( amp )
		sim.simulate(nIterations, metadata=[J, amp])
		# sim.save()









