def test_num_den_to_mesh_stdout():
if not HAVE_PYMOAB:
raise SkipTest
filename = os.path.join(thisdir, "files_test_alara",
"num_density_output.txt")
m = Mesh(structured=True, structured_coords=[[0, 1], [0, 1], [0, 1, 2]])
p = subprocess.Popen(["cat", filename], stdout=subprocess.PIPE)
lines, err = p.communicate()
lines = str_to_unicode(lines)
num_density_to_mesh(lines.split("\n"), "shutdown", m)
# expected composition results:
exp_comp_0 = {
10010000: 5.3390e19,
10020000: 3.0571e17,
10030000: 1.2082e12,
20030000: 7.4323e09,
20040000: 7.1632e02,
}
exp_comp_1 = {
10010000: 4.1240e13,
10020000: 4.7443e11,
10030000: 2.6627e13,
20030000: 8.3547e10,
20040000: 2.6877e19,
}
# actual composition results
act_comp_0 = m.mats[0].to_atom_frac()
act_comp_1 = m.mats[1].to_atom_frac()
assert_equal(len(exp_comp_0), len(act_comp_0))
for key, value in exp_comp_0.items():
assert_almost_equal(value / act_comp_0[key], 1.0, 15)
assert_equal(len(exp_comp_1), len(act_comp_1))
for key, value in exp_comp_1.items():
assert_almost_equal(value / act_comp_1[key], 1.0, 15)
# compare densities
exp_density_0 = 8.96715e-05
exp_density_1 = 1.785214e-04
assert_almost_equal(exp_density_0, m.mats[0].density)
assert_almost_equal(exp_density_1, m.mats[1].density)
def test_num_den_to_mesh_stdout():
if not HAVE_PYTAPS:
raise SkipTest
filename = os.path.join(thisdir, "files_test_alara", "num_density_output.txt")
m = Mesh(structured=True, structured_coords=[[0, 1], [0, 1], [0, 1, 2]])
p = subprocess.Popen(["cat", filename], stdout=subprocess.PIPE)
lines, err = p.communicate()
num_density_to_mesh(lines.split("\n"), "shutdown", m)
# expected composition results:
exp_comp_0 = {
10010000: 5.3390e19,
10020000: 3.0571e17,
10030000: 1.2082e12,
20030000: 7.4323e09,
20040000: 7.1632e02,
}
exp_comp_1 = {
10010000: 4.1240e13,
10020000: 4.7443e11,
10030000: 2.6627e13,
20030000: 8.3547e10,
20040000: 2.6877e19,
}
# actual composition results
act_comp_0 = m.mats[0].to_atom_frac()
act_comp_1 = m.mats[1].to_atom_frac()
assert_equal(len(exp_comp_0), len(act_comp_0))
for key, value in exp_comp_0.iteritems():
assert_almost_equal(value / act_comp_0[key], 1.0, 15)
assert_equal(len(exp_comp_1), len(act_comp_1))
for key, value in exp_comp_1.iteritems():
assert_almost_equal(value / act_comp_1[key], 1.0, 15)
# compare densities
exp_density_0 = 8.96715e-05
exp_density_1 = 1.785214e-04
assert_almost_equal(exp_density_0, m.mats[0].density)
assert_almost_equal(exp_density_1, m.mats[1].density)
def test_num_den_to_mesh_shutdown():
if not HAVE_PYTAPS:
raise SkipTest
filename = os.path.join(thisdir, "files_test_alara",
"num_density_output.txt")
m = Mesh(structured=True, structured_coords=[[0, 1], [0, 1], [0, 1, 2]])
with open(filename) as f:
lines = f.readlines()
num_density_to_mesh(lines, 'shutdown', m)
# expected composition results:
exp_comp_0 = {
10010000: 5.3390e+19,
10020000: 3.0571e+17,
10030000: 1.2082e+12,
20030000: 7.4323e+09,
20040000: 7.1632e+02
}
exp_comp_1 = {
10010000: 4.1240e+13,
10020000: 4.7443e+11,
10030000: 2.6627e+13,
20030000: 8.3547e+10,
20040000: 2.6877e+19
}
# actual composition results
act_comp_0 = m.mats[0].to_atom_frac()
act_comp_1 = m.mats[1].to_atom_frac()
assert_equal(len(exp_comp_0), len(act_comp_0))
for key, value in exp_comp_0.iteritems():
assert_almost_equal(value / act_comp_0[key], 1.0, 15)
assert_equal(len(exp_comp_1), len(act_comp_1))
for key, value in exp_comp_1.iteritems():
assert_almost_equal(value / act_comp_1[key], 1.0, 15)
# compare densities
exp_density_0 = 8.96715E-05
exp_density_1 = 1.785214E-04
assert_almost_equal(exp_density_0, m.mats[0].density)
assert_almost_equal(exp_density_1, m.mats[1].density)
def test_num_den_to_mesh_1_y():
if not HAVE_PYTAPS:
raise SkipTest
filename = os.path.join(thisdir, "files_test_alara",
"num_density_output.txt")
m = Mesh(structured=True, structured_coords=[[0, 1], [0, 1], [0, 1, 2]])
num_density_to_mesh(filename, '1 y', m)
# expected results:
exp_comp_0 = {
10010000: 5.3390e+19,
10020000: 3.0571e+17,
10030000: 1.1424e+12,
20030000: 7.3260e+10,
20040000: 7.1632e+02
}
exp_comp_1 = {
10010000: 4.1240e+13,
10020000: 4.7443e+11,
10030000: 2.5176e+13,
20030000: 1.5343e+12,
20040000: 2.6877e+19
}
# actual results
act_comp_0 = m.mats[0].to_atom_frac()
act_comp_1 = m.mats[1].to_atom_frac()
assert_equal(len(exp_comp_0), len(act_comp_0))
for key, value in exp_comp_0.iteritems():
assert_almost_equal(value / act_comp_0[key], 1.0, 15)
assert_equal(len(exp_comp_1), len(act_comp_1))
for key, value in exp_comp_1.iteritems():
assert_almost_equal(value / act_comp_1[key], 1.0, 15)
# compare densities
exp_density_0 = 8.96715E-05
exp_density_1 = 1.78521E-04
assert_almost_equal(exp_density_0, m.mats[0].density)
assert_almost_equal(exp_density_1, m.mats[1].density)
def test_num_den_to_mesh_1_y():
if not HAVE_PYTAPS:
raise SkipTest
filename = os.path.join(thisdir, "files_test_alara", "num_density_output.txt")
m = Mesh(structured=True, structured_coords=[[0, 1], [0, 1], [0, 1, 2]])
num_density_to_mesh(filename, "1 y", m)
# expected results:
exp_comp_0 = {
10010000: 5.3390e19,
10020000: 3.0571e17,
10030000: 1.1424e12,
20030000: 7.3260e10,
20040000: 7.1632e02,
}
exp_comp_1 = {
10010000: 4.1240e13,
10020000: 4.7443e11,
10030000: 2.5176e13,
20030000: 1.5343e12,
20040000: 2.6877e19,
}
# actual results
act_comp_0 = m.mats[0].to_atom_frac()
act_comp_1 = m.mats[1].to_atom_frac()
assert_equal(len(exp_comp_0), len(act_comp_0))
for key, value in exp_comp_0.iteritems():
assert_almost_equal(value / act_comp_0[key], 1.0, 15)
assert_equal(len(exp_comp_1), len(act_comp_1))
for key, value in exp_comp_1.iteritems():
assert_almost_equal(value / act_comp_1[key], 1.0, 15)
# compare densities
exp_density_0 = 8.96715e-05
exp_density_1 = 1.78521e-04
assert_almost_equal(exp_density_0, m.mats[0].density)
assert_almost_equal(exp_density_1, m.mats[1].density)
def test_num_den_to_mesh_shutdown():
if not HAVE_PYTAPS:
raise SkipTest
filename = os.path.join(thisdir, "files_test_alara",
"num_density_output.txt")
m = Mesh(structured=True, structured_coords=[[0,1],[0,1],[0,1,2]])
with open(filename) as f:
lines = f.readlines()
num_density_to_mesh(lines, 'shutdown', m)
# expected composition results:
exp_comp_0 = {10010000:5.3390e+19,
10020000:3.0571e+17,
10030000:1.2082e+12,
20030000:7.4323e+09,
20040000:7.1632e+02}
exp_comp_1 = {10010000:4.1240e+13,
10020000:4.7443e+11,
10030000:2.6627e+13,
20030000:8.3547e+10,
20040000:2.6877e+19}
# actual composition results
act_comp_0 = m.mats[0].to_atom_frac()
act_comp_1 = m.mats[1].to_atom_frac()
assert_equal(len(exp_comp_0), len(act_comp_0))
for key, value in exp_comp_0.iteritems():
assert_almost_equal(value/act_comp_0[key], 1.0, 15)
assert_equal(len(exp_comp_1), len(act_comp_1))
for key, value in exp_comp_1.iteritems():
assert_almost_equal(value/act_comp_1[key], 1.0, 15)
# compare densities
exp_density_0 = 8.96715E-05
exp_density_1 = 1.785214E-04
assert_almost_equal(exp_density_0, m.mats[0].density)
assert_almost_equal(exp_density_1, m.mats[1].density)
