def test_mesh_to_geom():
if not HAVE_PYTAPS:
raise SkipTest
mats = {
0: Material({'H1': 1.0, 'K39': 1.0}, density=42.0),
1: Material({'H1': 0.1, 'O16': 1.0}, density=43.0),
2: Material({'He4': 42.0}, density=44.0),
3: Material({'Tm171': 171.0}, density=45.0),
4: Material({'C12': 1.0}, density=47.0),
5: Material({'1002': 1.0}, density=5.0),
}
m = Mesh(structured_coords=[[0, 1, 2, 3], [0, 1, 2], [0, 1]], mats=mats,
structured=True)
geom = mcnp.mesh_to_geom(m)
exp_geom = (
"Generated from PyNE Mesh\n"
"1 1 42.0 1 -2 5 -6 8 -9\n"
"2 2 43.0 1 -2 6 -7 8 -9\n"
"3 3 44.0 2 -3 5 -6 8 -9\n"
"4 4 45.0 2 -3 6 -7 8 -9\n"
"5 5 47.0 3 -4 5 -6 8 -9\n"
"6 6 5.0 3 -4 6 -7 8 -9\n"
"7 0 -1:4:-5:7:-8:9\n"
"\n"
"1 px 0.0\n"
"2 px 1.0\n"
"3 px 2.0\n"
"4 px 3.0\n"
"5 py 0.0\n"
"6 py 1.0\n"
"7 py 2.0\n"
"8 pz 0.0\n"
"9 pz 1.0\n"
"\n"
"C density = 42.0\n"
"m1\n"
" 1001 -5.0000E-01\n"
" 19039 -5.0000E-01\n"
"C density = 43.0\n"
"m2\n"
" 1001 -9.0909E-02\n"
" 8016 -9.0909E-01\n"
"C density = 44.0\n"
"m3\n"
" 2004 -1.0000E+00\n"
"C density = 45.0\n"
"m4\n"
" 69171 -1.0000E+00\n"
"C density = 47.0\n"
"m5\n"
" 6012 -1.0000E+00\n"
"C density = 5.0\n"
"m6\n"
" 1002 -1.0000E+00\n")
assert_equal(geom, exp_geom)
def test_mesh_to_geom():
if not HAVE_PYTAPS:
raise SkipTest
mats = {
0: Material({'H1': 1.0, 'K39': 1.0}, density=42.0),
1: Material({'H1': 0.1, 'O16': 1.0}, density=43.0),
2: Material({'He4': 42.0}, density=44.0),
3: Material({'Tm171': 171.0}, density=45.0),
4: Material({'C12': 1.0}, density=47.0),
5: Material({'1002': 1.0}, density=5.0),
}
m = Mesh(structured_coords=[[0, 1, 2, 3], [0, 1, 2], [0, 1]], mats=mats,
structured=True)
geom = mcnp.mesh_to_geom(m)
exp_geom = (
"Generated from PyNE Mesh\n"
"1 1 42.0 1 -2 5 -6 8 -9\n"
"2 2 43.0 1 -2 6 -7 8 -9\n"
"3 3 44.0 2 -3 5 -6 8 -9\n"
"4 4 45.0 2 -3 6 -7 8 -9\n"
"5 5 47.0 3 -4 5 -6 8 -9\n"
"6 6 5.0 3 -4 6 -7 8 -9\n"
"7 0 -1:4:-5:7:-8:9\n"
"\n"
"1 px 0.0\n"
"2 px 1.0\n"
"3 px 2.0\n"
"4 px 3.0\n"
"5 py 0.0\n"
"6 py 1.0\n"
"7 py 2.0\n"
"8 pz 0.0\n"
"9 pz 1.0\n"
"\n"
"C density = 42.0\n"
"m1\n"
" 1001 -5.0000e-01\n"
" 19039 -5.0000e-01\n"
"C density = 43.0\n"
"m2\n"
" 1001 -9.0909e-02\n"
" 8016 -9.0909e-01\n"
"C density = 44.0\n"
"m3\n"
" 2004 -1.0000e+00\n"
"C density = 45.0\n"
"m4\n"
" 69171 -1.0000e+00\n"
"C density = 47.0\n"
"m5\n"
" 6012 -1.0000e+00\n"
"C density = 5.0\n"
"m6\n"
" 1002 -1.0000e+00\n")
assert_equal(geom, exp_geom)
def test_mesh_to_geom():
if not HAVE_PYMOAB:
raise SkipTest
mats = {
0: Material({"H1": 1.0, "K39": 1.0}, density=42.0),
1: Material({"H1": 0.1, "O16": 1.0}, density=43.0),
2: Material({"He4": 42.0}, density=44.0),
3: Material({"Tm171": 171.0}, density=45.0),
4: Material({"C12": 1.0}, density=47.0),
5: Material({"1002": 1.0}, density=5.0),
}
m = Mesh(
structured_coords=[[0, 1, 2, 3], [0, 1, 2], [0, 1]], mats=mats, structured=True
)
geom = mcnp.mesh_to_geom(m)
exp_geom = (
"Generated from PyNE Mesh\n"
"1 1 42.0 1 -2 5 -6 8 -9\n"
"2 2 43.0 1 -2 6 -7 8 -9\n"
"3 3 44.0 2 -3 5 -6 8 -9\n"
"4 4 45.0 2 -3 6 -7 8 -9\n"
"5 5 47.0 3 -4 5 -6 8 -9\n"
"6 6 5.0 3 -4 6 -7 8 -9\n"
"7 0 -1:4:-5:7:-8:9\n"
"\n"
"1 px 0.0\n"
"2 px 1.0\n"
"3 px 2.0\n"
"4 px 3.0\n"
"5 py 0.0\n"
"6 py 1.0\n"
"7 py 2.0\n"
"8 pz 0.0\n"
"9 pz 1.0\n"
"\n"
"C name: 0\n"
"C density = 42.00000\n"
"m1\n"
" 1001 -2.1000e+01\n"
" 19039 -2.1000e+01\n"
"C name: 1\n"
"C density = 43.00000\n"
"m2\n"
" 1001 -3.9091e+00\n"
" 8016 -3.9091e+01\n"
"C name: 2\n"
"C density = 44.00000\n"
"m3\n"
" 2004 -4.4000e+01\n"
"C name: 3\n"
"C density = 45.00000\n"
"m4\n"
" 69171 -4.5000e+01\n"
"C name: 4\n"
"C density = 47.00000\n"
"m5\n"
" 6012 -4.7000e+01\n"
"C name: 5\n"
"C density = 5.00000\n"
"m6\n"
" 1002 -5.0000e+00\n"
)
assert_equal(geom, exp_geom)