H1 = Isotope(ace_path + 'H_001_293.6K.ace')
# -------------------------------
# parameters
# -------------------------------
tau = 1.0e-3
T = 293.5
N = 2000 # # of neutrons want to sample
E = pre_computing(tau)
#--------------------------------------------------------
# case 1, 1000 / 1
#--------------------------------------------------------
mat = Material(E)
mat.addIsotope(U235, 1.0)
mat.addIsotope(H1, 1000.0)
mat.normalize_numdensity()
mc = mc_sampling(E, mat, T, N)
E_for_flux, flux, deviation = mc.solver()
np.savetxt('../results/E_flux.txt', E_for_flux)
np.savetxt('../results/flux1000to1.txt', flux)
np.savetxt('../results/deviation1000to1.txt', deviation)
#--------------------------------------------------------
# case 2, 100 / 1
#--------------------------------------------------------
mat = Material(E)
print 'length mat total ', len(mat.macro_total)
