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)