crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 22:26:08 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

got mesh domain working for mgxs

Browse source
This commit is contained in:
Sam Shaner 2016-06-28 13:32:27 -04:00
parent d30d010aea
commit 687e81b0dd
8 changed files with 1017 additions and 381 deletions
Download patch file Download diff file Expand all files Collapse all files

403
docs/source/pythonapi/examples/mgxs-part-i.ipynb
View file

@ -183,7 +183,7 @@
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": true
"collapsed": false
},
"outputs": [],
"source": [
@ -208,7 +208,7 @@
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": true
"collapsed": false
},
"outputs": [],
"source": [
@ -229,7 +229,7 @@
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": true
"collapsed": false
},
"outputs": [],
"source": [
@ -324,9 +324,9 @@
"outputs": [],
"source": [
"# OpenMC simulation parameters\n",
"batches = 50\n",
"batches = 20\n",
"inactive = 10\n",
"particles = 2500\n",
"particles = 1000\n",
"\n",
"# Instantiate a Settings object\n",
"settings_file = openmc.Settings()\n",
@ -395,10 +395,20 @@
},
"outputs": [],
"source": [
"# Instantiate a tally Mesh\n",
"mesh = openmc.Mesh(name='mesh')\n",
"mesh.type = 'regular'\n",
"mesh.dimension = [1, 1]\n",
"mesh.lower_left = [-0.63, -0.63]\n",
"mesh.upper_right = [+0.63, +0.63]\n",
"\n",
"# Instantiate a few different sections\n",
"total = mgxs.TotalXS(domain=cell, groups=groups)\n",
"absorption = mgxs.AbsorptionXS(domain=cell, groups=groups)\n",
"scattering = mgxs.ScatterXS(domain=cell, groups=groups)"
"total = mgxs.TotalXS(domain=mesh, groups=groups)\n",
"absorption = mgxs.AbsorptionXS(domain=mesh, groups=groups)\n",
"scattering = mgxs.ScatterXS(domain=mesh, groups=groups)\n",
"#total = mgxs.TotalXS(domain=cell, groups=groups)\n",
"#absorption = mgxs.AbsorptionXS(domain=cell, groups=groups)\n",
"#scattering = mgxs.ScatterXS(domain=cell, groups=groups)"
]
},
{
@ -419,24 +429,22 @@
"data": {
"text/plain": [
"OrderedDict([('flux', Tally\n",
"\tID =\t10000\n",
"\tName =\t\n",
"\tFilters =\t\n",
" \t\tcell\t[1]\n",
" \t\tenergy\t[ 0.00000000e+00 6.25000000e-07 2.00000000e+01]\n",
"\tNuclides =\ttotal \n",
"\tScores =\t['flux']\n",
"\tEstimator =\ttracklength\n",
"), ('absorption', Tally\n",
"\tID =\t10001\n",
"\tName =\t\n",
"\tFilters =\t\n",
" \t\tcell\t[1]\n",
" \t\tenergy\t[ 0.00000000e+00 6.25000000e-07 2.00000000e+01]\n",
"\tNuclides =\ttotal \n",
"\tScores =\t['absorption']\n",
"\tEstimator =\ttracklength\n",
")])"
" \tID =\t10000\n",
" \tName =\t\n",
" \tFilters =\t\n",
" \t\tmesh\t[10000]\n",
" \t\tenergy\t[ 0.00000000e+00 6.25000000e-07 2.00000000e+01]\n",
" \tNuclides =\ttotal \n",
" \tScores =\t['flux']\n",
" \tEstimator =\ttracklength), ('absorption', Tally\n",
" \tID =\t10001\n",
" \tName =\t\n",
" \tFilters =\t\n",
" \t\tmesh\t[10000]\n",
" \t\tenergy\t[ 0.00000000e+00 6.25000000e-07 2.00000000e+01]\n",
" \tNuclides =\ttotal \n",
" \tScores =\t['absorption']\n",
" \tEstimator =\ttracklength)])"
]
},
"execution_count": 13,
@ -513,8 +521,8 @@
" Copyright: 2011-2016 Massachusetts Institute of Technology\n",
" License: http://openmc.readthedocs.io/en/latest/license.html\n",
" Version: 0.7.1\n",
" Git SHA1: 19feb55e6d5e8350398627f39fb55ee8e2e63011\n",
" Date/Time: 2016-05-13 10:19:16\n",
" Git SHA1: d30d010aea2c1cba90c993a9c501be9d5d81921d\n",
" Date/Time: 2016-06-28 13:31:05\n",
" MPI Processes: 1\n",
"\n",
" ===========================================================================\n",
@ -541,57 +549,27 @@
"\n",
" Bat./Gen. k Average k \n",
" ========= ======== ==================== \n",
" 1/1 1.11184 \n",
" 2/1 1.15820 \n",
" 3/1 1.18468 \n",
" 4/1 1.17492 \n",
" 5/1 1.19645 \n",
" 6/1 1.18436 \n",
" 7/1 1.14070 \n",
" 8/1 1.15150 \n",
" 9/1 1.19202 \n",
" 10/1 1.17677 \n",
" 11/1 1.20272 \n",
" 12/1 1.21366 1.20819 +/- 0.00547\n",
" 13/1 1.15906 1.19181 +/- 0.01668\n",
" 14/1 1.14687 1.18058 +/- 0.01629\n",
" 15/1 1.14570 1.17360 +/- 0.01442\n",
" 16/1 1.13480 1.16713 +/- 0.01343\n",
" 17/1 1.17680 1.16852 +/- 0.01144\n",
" 18/1 1.16866 1.16853 +/- 0.00990\n",
" 19/1 1.19253 1.17120 +/- 0.00913\n",
" 20/1 1.18124 1.17220 +/- 0.00823\n",
" 21/1 1.19206 1.17401 +/- 0.00766\n",
" 22/1 1.17681 1.17424 +/- 0.00700\n",
" 23/1 1.17634 1.17440 +/- 0.00644\n",
" 24/1 1.13659 1.17170 +/- 0.00654\n",
" 25/1 1.17144 1.17169 +/- 0.00609\n",
" 26/1 1.20649 1.17386 +/- 0.00610\n",
" 27/1 1.11238 1.17024 +/- 0.00678\n",
" 28/1 1.18911 1.17129 +/- 0.00647\n",
" 29/1 1.14681 1.17000 +/- 0.00626\n",
" 30/1 1.12152 1.16758 +/- 0.00641\n",
" 31/1 1.12729 1.16566 +/- 0.00639\n",
" 32/1 1.15399 1.16513 +/- 0.00612\n",
" 33/1 1.13547 1.16384 +/- 0.00599\n",
" 34/1 1.17723 1.16440 +/- 0.00576\n",
" 35/1 1.09296 1.16154 +/- 0.00622\n",
" 36/1 1.19621 1.16287 +/- 0.00612\n",
" 37/1 1.12560 1.16149 +/- 0.00605\n",
" 38/1 1.17872 1.16211 +/- 0.00586\n",
" 39/1 1.17721 1.16263 +/- 0.00568\n",
" 40/1 1.13724 1.16178 +/- 0.00555\n",
" 41/1 1.18526 1.16254 +/- 0.00542\n",
" 42/1 1.13779 1.16177 +/- 0.00531\n",
" 43/1 1.15066 1.16143 +/- 0.00516\n",
" 44/1 1.12174 1.16026 +/- 0.00514\n",
" 45/1 1.17479 1.16068 +/- 0.00501\n",
" 46/1 1.14146 1.16014 +/- 0.00489\n",
" 47/1 1.20464 1.16135 +/- 0.00491\n",
" 48/1 1.15119 1.16108 +/- 0.00479\n",
" 49/1 1.17938 1.16155 +/- 0.00468\n",
" 50/1 1.15798 1.16146 +/- 0.00457\n",
" Creating state point statepoint.50.h5...\n",
" 1/1 1.07993 \n",
" 2/1 1.23691 \n",
" 3/1 1.17407 \n",
" 4/1 1.08258 \n",
" 5/1 1.21012 \n",
" 6/1 1.23825 \n",
" 7/1 1.18016 \n",
" 8/1 1.20989 \n",
" 9/1 1.15475 \n",
" 10/1 1.13462 \n",
" 11/1 1.12749 \n",
" 12/1 1.09502 1.11125 +/- 0.01623\n",
" 13/1 1.24722 1.15657 +/- 0.04628\n",
" 14/1 1.14700 1.15418 +/- 0.03281\n",
" 15/1 1.13889 1.15112 +/- 0.02560\n",
" 16/1 1.19008 1.15762 +/- 0.02189\n",
" 17/1 1.11320 1.15127 +/- 0.01956\n",
" 18/1 1.12093 1.14748 +/- 0.01736\n",
" 19/1 1.13809 1.14643 +/- 0.01534\n",
" 20/1 1.09886 1.14168 +/- 0.01452\n",
" Creating state point statepoint.20.h5...\n",
"\n",
" ===========================================================================\n",
" ======================> SIMULATION FINISHED <======================\n",
@ -600,27 +578,27 @@
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 4.2300E-01 seconds\n",
" Reading cross sections = 9.3000E-02 seconds\n",
" Total time in simulation = 1.6549E+01 seconds\n",
" Time in transport only = 1.6535E+01 seconds\n",
" Time in inactive batches = 2.3650E+00 seconds\n",
" Time in active batches = 1.4184E+01 seconds\n",
" Time synchronizing fission bank = 5.0000E-03 seconds\n",
" Sampling source sites = 3.0000E-03 seconds\n",
" Total time for initialization = 8.1200E-01 seconds\n",
" Reading cross sections = 2.4200E-01 seconds\n",
" Total time in simulation = 2.8910E+00 seconds\n",
" Time in transport only = 2.8710E+00 seconds\n",
" Time in inactive batches = 8.5700E-01 seconds\n",
" Time in active batches = 2.0340E+00 seconds\n",
" Time synchronizing fission bank = 0.0000E+00 seconds\n",
" Sampling source sites = 0.0000E+00 seconds\n",
" SEND/RECV source sites = 0.0000E+00 seconds\n",
" Time accumulating tallies = 0.0000E+00 seconds\n",
" Total time for finalization = 0.0000E+00 seconds\n",
" Total time elapsed = 1.6981E+01 seconds\n",
" Calculation Rate (inactive) = 10570.8 neutrons/second\n",
" Calculation Rate (active) = 7050.20 neutrons/second\n",
" Total time for finalization = 1.0000E-03 seconds\n",
" Total time elapsed = 3.7090E+00 seconds\n",
" Calculation Rate (inactive) = 11668.6 neutrons/second\n",
" Calculation Rate (active) = 4916.42 neutrons/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" k-effective (Collision) = 1.15984 +/- 0.00411\n",
" k-effective (Track-length) = 1.16146 +/- 0.00457\n",
" k-effective (Absorption) = 1.16177 +/- 0.00380\n",
" Combined k-effective = 1.16105 +/- 0.00364\n",
" k-effective (Collision) = 1.14223 +/- 0.01330\n",
" k-effective (Track-length) = 1.14168 +/- 0.01452\n",
" k-effective (Absorption) = 1.16108 +/- 0.00876\n",
" Combined k-effective = 1.15408 +/- 0.00566\n",
" Leakage Fraction = 0.00000 +/- 0.00000\n",
"\n"
]
@ -664,7 +642,7 @@
"outputs": [],
"source": [
"# Load the last statepoint file\n",
"sp = openmc.StatePoint('statepoint.50.h5')"
"sp = openmc.StatePoint('statepoint.20.h5')"
]
},
{
@ -729,11 +707,11 @@
"text": [
"Multi-Group XS\n",
"\tReaction Type =\ttotal\n",
"\tDomain Type =\tcell\n",
"\tDomain ID =\t1\n",
"\tDomain Type =\tmesh\n",
"\tDomain ID =\t10000\n",
"\tCross Sections [cm^-1]:\n",
" Group 1 [6.25e-07 - 20.0 MeV]:\t6.81e-01 +/- 2.69e-01%\n",
" Group 2 [0.0 - 6.25e-07 MeV]:\t1.40e+00 +/- 5.93e-01%\n",
" Group 1 [6.25e-07 - 20.0 MeV]:\t6.79e-01 +/- 7.64e-01%\n",
" Group 2 [0.0 - 6.25e-07 MeV]:\t1.40e+00 +/- 1.79e+00%\n",
"\n",
"\n",
"\n"
@ -764,40 +742,55 @@
"<div>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <tr>\n",
" <th></th>\n",
" <th>cell</th>\n",
" <th colspan=\"3\" halign=\"left\">mesh 10000</th>\n",
" <th>group in</th>\n",
" <th>nuclide</th>\n",
" <th>mean</th>\n",
" <th>std. dev.</th>\n",
" </tr>\n",
" <tr>\n",
" <th></th>\n",
" <th>x</th>\n",
" <th>y</th>\n",
" <th>z</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>total</td>\n",
" <td>0.667787</td>\n",
" <td>0.001802</td>\n",
" <td>0.666261</td>\n",
" <td>0.005072</td>\n",
" </tr>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>2</td>\n",
" <td>total</td>\n",
" <td>1.292013</td>\n",
" <td>0.007642</td>\n",
" <td>1.292451</td>\n",
" <td>0.022969</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" cell group in nuclide mean std. dev.\n",
"1 1 1 total 0.667787 0.001802\n",
"0 1 2 total 1.292013 0.007642"
" mesh 10000 group in nuclide mean std. dev.\n",
" x y z \n",
"1 1 1 1 1 total 0.666261 0.005072\n",
"0 1 1 1 2 total 1.292451 0.022969"
]
},
"execution_count": 19,
@ -821,11 +814,11 @@
"cell_type": "code",
"execution_count": 20,
"metadata": {
"collapsed": true
"collapsed": false
},
"outputs": [],
"source": [
"absorption.export_xs_data(filename='absorption-xs', format='excel')"
"#absorption.export_xs_data(filename='absorption-xs', format='excel')"
]
},
{
@ -875,9 +868,9 @@
"<div>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <tr>\n",
" <th></th>\n",
" <th>cell</th>\n",
" <th colspan=\"3\" halign=\"left\">mesh 10000</th>\n",
" <th>energy low [MeV]</th>\n",
" <th>energy high [MeV]</th>\n",
" <th>nuclide</th>\n",
@ -885,40 +878,58 @@
" <th>mean</th>\n",
" <th>std. dev.</th>\n",
" </tr>\n",
" <tr>\n",
" <th></th>\n",
" <th>x</th>\n",
" <th>y</th>\n",
" <th>z</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>0.000000e+00</td>\n",
" <td>6.250000e-07</td>\n",
" <td>total</td>\n",
" <td>(((total / flux) - (absorption / flux)) - (sca...</td>\n",
" <td>-3.774758e-15</td>\n",
" <td>0.011292</td>\n",
" <td>-1.554312e-15</td>\n",
" <td>0.033991</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>6.250000e-07</td>\n",
" <td>2.000000e+01</td>\n",
" <td>total</td>\n",
" <td>(((total / flux) - (absorption / flux)) - (sca...</td>\n",
" <td>1.443290e-15</td>\n",
" <td>0.002570</td>\n",
" <td>-2.886580e-15</td>\n",
" <td>0.007258</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" cell energy low [MeV] energy high [MeV] nuclide \\\n",
"0 1 0.00e+00 6.25e-07 total \n",
"1 1 6.25e-07 2.00e+01 total \n",
" mesh 10000 energy low [MeV] energy high [MeV] nuclide \\\n",
" x y z \n",
"0 1 1 1 0.00e+00 6.25e-07 total \n",
"1 1 1 1 6.25e-07 2.00e+01 total \n",
"\n",
" score mean std. dev. \n",
"0 (((total / flux) - (absorption / flux)) - (sca... -3.77e-15 1.13e-02 \n",
"1 (((total / flux) - (absorption / flux)) - (sca... 1.44e-15 2.57e-03 "
" score mean std. dev. \n",
" \n",
"0 (((total / flux) - (absorption / flux)) - (sca... -1.55e-15 3.40e-02 \n",
"1 (((total / flux) - (absorption / flux)) - (sca... -2.89e-15 7.26e-03 "
]
},
"execution_count": 22,
@ -954,9 +965,9 @@
"<div>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <tr>\n",
" <th></th>\n",
" <th>cell</th>\n",
" <th colspan=\"3\" halign=\"left\">mesh 10000</th>\n",
" <th>energy low [MeV]</th>\n",
" <th>energy high [MeV]</th>\n",
" <th>nuclide</th>\n",
@ -964,40 +975,58 @@
" <th>mean</th>\n",
" <th>std. dev.</th>\n",
" </tr>\n",
" <tr>\n",
" <th></th>\n",
" <th>x</th>\n",
" <th>y</th>\n",
" <th>z</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>0.000000e+00</td>\n",
" <td>6.250000e-07</td>\n",
" <td>total</td>\n",
" <td>((absorption / flux) / (total / flux))</td>\n",
" <td>0.076115</td>\n",
" <td>0.000649</td>\n",
" <td>0.076144</td>\n",
" <td>0.001984</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>6.250000e-07</td>\n",
" <td>2.000000e+01</td>\n",
" <td>total</td>\n",
" <td>((absorption / flux) / (total / flux))</td>\n",
" <td>0.019263</td>\n",
" <td>0.000095</td>\n",
" <td>0.019224</td>\n",
" <td>0.000329</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" cell energy low [MeV] energy high [MeV] nuclide \\\n",
"0 1 0.00e+00 6.25e-07 total \n",
"1 1 6.25e-07 2.00e+01 total \n",
" mesh 10000 energy low [MeV] energy high [MeV] nuclide \\\n",
" x y z \n",
"0 1 1 1 0.00e+00 6.25e-07 total \n",
"1 1 1 1 6.25e-07 2.00e+01 total \n",
"\n",
" score mean std. dev. \n",
"0 ((absorption / flux) / (total / flux)) 7.61e-02 6.49e-04 \n",
"1 ((absorption / flux) / (total / flux)) 1.93e-02 9.46e-05 "
" score mean std. dev. \n",
" \n",
"0 ((absorption / flux) / (total / flux)) 7.61e-02 1.98e-03 \n",
"1 ((absorption / flux) / (total / flux)) 1.92e-02 3.29e-04 "
]
},
"execution_count": 23,
@ -1026,9 +1055,9 @@
"<div>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <tr>\n",
" <th></th>\n",
" <th>cell</th>\n",
" <th colspan=\"3\" halign=\"left\">mesh 10000</th>\n",
" <th>energy low [MeV]</th>\n",
" <th>energy high [MeV]</th>\n",
" <th>nuclide</th>\n",
@ -1036,40 +1065,58 @@
" <th>mean</th>\n",
" <th>std. dev.</th>\n",
" </tr>\n",
" <tr>\n",
" <th></th>\n",
" <th>x</th>\n",
" <th>y</th>\n",
" <th>z</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>0.000000e+00</td>\n",
" <td>6.250000e-07</td>\n",
" <td>total</td>\n",
" <td>((scatter / flux) / (total / flux))</td>\n",
" <td>0.923885</td>\n",
" <td>0.007736</td>\n",
" <td>0.923856</td>\n",
" <td>0.023272</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>6.250000e-07</td>\n",
" <td>2.000000e+01</td>\n",
" <td>total</td>\n",
" <td>((scatter / flux) / (total / flux))</td>\n",
" <td>0.980737</td>\n",
" <td>0.003737</td>\n",
" <td>0.980776</td>\n",
" <td>0.010576</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" cell energy low [MeV] energy high [MeV] nuclide \\\n",
"0 1 0.00e+00 6.25e-07 total \n",
"1 1 6.25e-07 2.00e+01 total \n",
" mesh 10000 energy low [MeV] energy high [MeV] nuclide \\\n",
" x y z \n",
"0 1 1 1 0.00e+00 6.25e-07 total \n",
"1 1 1 1 6.25e-07 2.00e+01 total \n",
"\n",
" score mean std. dev. \n",
"0 ((scatter / flux) / (total / flux)) 9.24e-01 7.74e-03 \n",
"1 ((scatter / flux) / (total / flux)) 9.81e-01 3.74e-03 "
" score mean std. dev. \n",
" \n",
"0 ((scatter / flux) / (total / flux)) 9.24e-01 2.33e-02 \n",
"1 ((scatter / flux) / (total / flux)) 9.81e-01 1.06e-02 "
]
},
"execution_count": 24,
@ -1105,9 +1152,9 @@
"<div>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <tr>\n",
" <th></th>\n",
" <th>cell</th>\n",
" <th colspan=\"3\" halign=\"left\">mesh 10000</th>\n",
" <th>energy low [MeV]</th>\n",
" <th>energy high [MeV]</th>\n",
" <th>nuclide</th>\n",
@ -1115,40 +1162,58 @@
" <th>mean</th>\n",
" <th>std. dev.</th>\n",
" </tr>\n",
" <tr>\n",
" <th></th>\n",
" <th>x</th>\n",
" <th>y</th>\n",
" <th>z</th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" <th></th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>0.000000e+00</td>\n",
" <td>6.250000e-07</td>\n",
" <td>total</td>\n",
" <td>(((absorption / flux) / (total / flux)) + ((sc...</td>\n",
" <td>1.0</td>\n",
" <td>0.007763</td>\n",
" <td>1</td>\n",
" <td>0.023356</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>1</td>\n",
" <td>6.250000e-07</td>\n",
" <td>2.000000e+01</td>\n",
" <td>total</td>\n",
" <td>(((absorption / flux) / (total / flux)) + ((sc...</td>\n",
" <td>1.0</td>\n",
" <td>0.003739</td>\n",
" <td>1</td>\n",
" <td>0.010581</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>"
],
"text/plain": [
" cell energy low [MeV] energy high [MeV] nuclide \\\n",
"0 1 0.00e+00 6.25e-07 total \n",
"1 1 6.25e-07 2.00e+01 total \n",
" mesh 10000 energy low [MeV] energy high [MeV] nuclide \\\n",
" x y z \n",
"0 1 1 1 0.00e+00 6.25e-07 total \n",
"1 1 1 1 6.25e-07 2.00e+01 total \n",
"\n",
" score mean std. dev. \n",
"0 (((absorption / flux) / (total / flux)) + ((sc... 1.00e+00 7.76e-03 \n",
"1 (((absorption / flux) / (total / flux)) + ((sc... 1.00e+00 3.74e-03 "
" score mean std. dev. \n",
" \n",
"0 (((absorption / flux) / (total / flux)) + ((sc... 1.00e+00 2.34e-02 \n",
"1 (((absorption / flux) / (total / flux)) + ((sc... 1.00e+00 1.06e-02 "
]
},
"execution_count": 25,
@ -1163,6 +1228,24 @@
"# The scattering-to-total ratio is a derived tally which can generate Pandas DataFrames for inspection\n",
"sum_ratio.get_pandas_dataframe()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
@ -1181,7 +1264,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.6"
"version": "2.7.11"
}
},
"nbformat": 4,

239
docs/source/pythonapi/examples/mgxs-part-iii.ipynb
View file

@ -32,7 +32,7 @@
"name": "stderr",
"output_type": "stream",
"text": [
"/home/wboyd/anaconda2/lib/python2.7/site-packages/matplotlib/__init__.py:1350: UserWarning: This call to matplotlib.use() has no effect\n",
"/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/matplotlib/__init__.py:1350: UserWarning: This call to matplotlib.use() has no effect\n",
"because the backend has already been chosen;\n",
"matplotlib.use() must be called *before* pylab, matplotlib.pyplot,\n",
"or matplotlib.backends is imported for the first time.\n",
@ -94,7 +94,7 @@
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": true
"collapsed": false
},
"outputs": [],
"source": [
@ -458,7 +458,7 @@
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAD6AgMAAAD1grKuAAAABGdBTUEAALGPC/xhBQAAACBjSFJN\nAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEX///9yEhLpgJFNv8Tq\nQYT7AAAAAWJLR0QAiAUdSAAAAAd0SU1FB+AFDhAdBviGIzgAAAWFSURBVGje7Zs7cttADIZ9CSvX\ncrP0iCxUqbBc8Ag6xR6BhV2EvYvwFD4CCx1ABT1jMdgndpegRQnOrCbjpPlGESISC4A/gd27e8H5\n83CX3b4+iKJrRHkS4vkghMPBonRYWGwtfgD2YN+dRDUOoh6lACw0Noi9w2fESuEoAR/uVuMolX03\n9oXGT7F3eFL2iEfhUX1f4cPdL/ishs+68ai+udE4xPhexbjX2FfjGNoPj/DPNX4Tsd+EODr8FvsV\ndf1Hd9P2VvCi4+s/aXvrf+upAD+1/9GV1mkOH5X9vV6THtfvACslcaUCbESL61drBPtdI8SrFMWr\nELsXCkuFDYW75gbiP7d9Cf7bAYI/aCwUShrBvh30+lWQkzVgZ/HD4OixNCgcQpJ3BxU/Ln91elKo\nM5VEE38QtJ+Yv6cQ9xjKNYayyl8TypP8DfJnQ2H/b/N3ye9P83cT33SQv/sQh9gV7zZ/0dNj5HQa\nC5vVzv9+/WFN2w8KVaZ2BwL1+pv4g0x1QRfjq0dB4Q3kT277oP6VNL6gKxNU9a8zK+WLbi/Wwpdi\nhbboKqyxFOulHMj6v4W/AXbmUeAxrv9J/CqEBXaRKsXaodD4nsYvkT/G6H1D4SR/iPy1Roj9JsQ5\ne18/7EUHv1+Fvx/Xj5V9Ugb5K8TW4TZEEdcvoz/up0VTe9qsVIppKVX6a7D6y9ZvwEKjrtQxPtv6\nfXII9vCxKOGaIeAIfEF8IvAG8ie3vRK9rRQl+PPpSctbhfpTUCpviH+kxsZgpT91+snoX1l49KK3\niUQvICRy5aUw6l8leoVwoo3Uv1rKreF/UFLY6d9QP4L9Wf2r7EP9GOSfcsjZ56f60kz+XmVPXv+R\nuP49ff0T/53Rv6n/7m2lvXT9Wqd/VUz8hvh5M/ED6ILmt4mfHYZSaePnTWpsf/SvqV9O6dLYYClL\nEetnoH/LBLFoBvrX189uTv8++kot5vTvQD4/9jP690g9P/4z/bvo/XVG/xYoZZx+8fr3MxAtsf7t\nUOkG2JqsTtCIpgCt/qX1226KqZS7gfzJbe+c9jLrtIZ8lXD+s4umlW6AKIVrlML2/cXjgPFjlJqI\nRC+Fj0bVJe+vSh56pSdR6YkQ1ygF10Wqf0FeLta/iKn9Mv1L24ti2e+7W4n1b3T/W+L+t9H9T/Sv\nVboUmqJJon1/hZq8LnzRDlDrX1u0xRT1+6vEpomMmyYkqi95vIH8yW1PN+122KkLcNLKi/WTF01z\n/cNASrWE/l3ev6T17zX909z9X27/euK/Rf3zWP+Waf9eEv37KkWJ+rfDl6ZglNDa+cEBhwYDvkoN\nP/rX69814NaI3imq0l7OYDy/qSdDGwr7r+Y3VbzoKZr6XX2lfxfOb87qXzr+b1j/Xlp/nP6dn98M\ncdH7cn7zjPObKsYWS3Eb9w8n85smHtqQuPuZ30T2dlIT6F9xFl+n8xslegL9a4c2KRr9W4rp/GYq\numiM9Nec/j2v/yj9u1h//hv9e93vc++f63/u+rPjL3f+5Lbn1j9m/eXWf+7zh/v8+2b9e/Hzn6s/\nuPqHrb8g71n6L3f+5Lbnvn8w33+4718/+5d47//c/gO7/5E7/nPbc/tv3P4fs//I7X9y+6/fqH+v\n6j9z+9/c/ju3/8+eP+TOn9z23PkXc/7Gnf9x5483q38Xzn+582fu/Js9fy8kb/6fO39y23P3n3S8\n/S/c/Tfc/T83uX/pgv1XE/9duP+Lu/+Mvf8td/znti8kb/8ld/9nx9t/Sjw/Ltr/yt1/+337f6/b\nf0zoB3nJ/ucVc/81d/83e/957vzJbc89/8A8f8E9/5HE78XnT/4H/cs5f8Q9/8Q9f8U+/5U7f3Lb\nc88fdrzzjyvm+cuf/Uu887/c88fs88954/8vO4SjPC+2QRIAAAAldEVYdGRhdGU6Y3JlYXRlADIw\nMTYtMDUtMTRUMTI6Mjk6MDYtMDQ6MDAGVIh3AAAAJXRFWHRkYXRlOm1vZGlmeQAyMDE2LTA1LTE0\nVDEyOjI5OjA2LTA0OjAwdwkwywAAAABJRU5ErkJggg==\n",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAD6AgMAAAD1grKuAAAABGdBTUEAALGPC/xhBQAAACBjSFJN\nAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEX///9yEhLpgJFNv8Tq\nQYT7AAAAAWJLR0QAiAUdSAAAAAd0SU1FB+AGEQw6AwjPAt8AAAWFSURBVGje7Zs7cttADIZ9CSvX\ncrP0iCxUqbBc8Ag6xR6BhV2EvYvwFD4CCx1ABT1jMdgndpegRQnOrCbjpPlGESISC4A/gd27e8H5\n83CX3b4+iKJrRHkS4vkghMPBonRYWGwtfgD2YN+dRDUOoh6lACw0Noi9w2fESuEoAR/uVuMolX03\n9oXGT7F3eFL2iEfhUX1f4cPdL/ishs+68ai+udE4xPhexbjX2FfjGNoPj/DPNX4Tsd+EODr8FvsV\ndf1Hd9P2VvCi4+s/aXvrf+upAD+1/9GV1mkOH5X9vV6THtfvACslcaUCbESL61drBPtdI8SrFMWr\nELsXCkuFDYW75gbiP7d9Cf7bAYI/aCwUShrBvh30+lWQkzVgZ/HD4OixNCgcQpJ3BxU/Ln91elKo\nM5VEE38QtJ+Yv6cQ9xjKNYayyl8TypP8DfJnQ2H/b/N3ye9P83cT33SQv/sQh9gV7zZ/0dNj5HQa\nC5vVzv9+/WFN2w8KVaZ2BwL1+pv4g0x1QRfjq0dB4Q3kT277oP6VNL6gKxNU9a8zK+WLbi/Wwpdi\nhbboKqyxFOulHMj6v4W/AXbmUeAxrv9J/CqEBXaRKsXaodD4nsYvkT/G6H1D4SR/iPy1Roj9JsQ5\ne18/7EUHv1+Fvx/Xj5V9Ugb5K8TW4TZEEdcvoz/up0VTe9qsVIppKVX6a7D6y9ZvwEKjrtQxPtv6\nfXII9vCxKOGaIeAIfEF8IvAG8ie3vRK9rRQl+PPpSctbhfpTUCpviH+kxsZgpT91+snoX1l49KK3\niUQvICRy5aUw6l8leoVwoo3Uv1rKreF/UFLY6d9QP4L9Wf2r7EP9GOSfcsjZ56f60kz+XmVPXv+R\nuP49ff0T/53Rv6n/7m2lvXT9Wqd/VUz8hvh5M/ED6ILmt4mfHYZSaePnTWpsf/SvqV9O6dLYYClL\nEetnoH/LBLFoBvrX189uTv8++kot5vTvQD4/9jP690g9P/4z/bvo/XVG/xYoZZx+8fr3MxAtsf7t\nUOkG2JqsTtCIpgCt/qX1226KqZS7gfzJbe+c9jLrtIZ8lXD+s4umlW6AKIVrlML2/cXjgPFjlJqI\nRC+Fj0bVJe+vSh56pSdR6YkQ1ygF10Wqf0FeLta/iKn9Mv1L24ti2e+7W4n1b3T/W+L+t9H9T/Sv\nVboUmqJJon1/hZq8LnzRDlDrX1u0xRT1+6vEpomMmyYkqi95vIH8yW1PN+122KkLcNLKi/WTF01z\n/cNASrWE/l3ev6T17zX909z9X27/euK/Rf3zWP+Waf9eEv37KkWJ+rfDl6ZglNDa+cEBhwYDvkoN\nP/rX69814NaI3imq0l7OYDy/qSdDGwr7r+Y3VbzoKZr6XX2lfxfOb87qXzr+b1j/Xlp/nP6dn98M\ncdH7cn7zjPObKsYWS3Eb9w8n85smHtqQuPuZ30T2dlIT6F9xFl+n8xslegL9a4c2KRr9W4rp/GYq\numiM9Nec/j2v/yj9u1h//hv9e93vc++f63/u+rPjL3f+5Lbn1j9m/eXWf+7zh/v8+2b9e/Hzn6s/\nuPqHrb8g71n6L3f+5Lbnvn8w33+4718/+5d47//c/gO7/5E7/nPbc/tv3P4fs//I7X9y+6/fqH+v\n6j9z+9/c/ju3/8+eP+TOn9z23PkXc/7Gnf9x5483q38Xzn+582fu/Js9fy8kb/6fO39y23P3n3S8\n/S/c/Tfc/T83uX/pgv1XE/9duP+Lu/+Mvf8td/znti8kb/8ld/9nx9t/Sjw/Ltr/yt1/+337f6/b\nf0zoB3nJ/ucVc/81d/83e/957vzJbc89/8A8f8E9/5HE78XnT/4H/cs5f8Q9/8Q9f8U+/5U7f3Lb\nc88fdrzzjyvm+cuf/Uu887/c88fs88954/8vO4SjPC+2QRIAAAAldEVYdGRhdGU6Y3JlYXRlADIw\nMTYtMDYtMTdUMTI6NTg6MDMtMDQ6MDAYC+s2AAAAJXRFWHRkYXRlOm1vZGlmeQAyMDE2LTA2LTE3\nVDEyOjU4OjAzLTA0OjAwaVZTigAAAABJRU5ErkJggg==\n",
"text/plain": [
"<IPython.core.display.Image object>"
]
@ -557,14 +557,14 @@
},
{
"cell_type": "code",
"execution_count": 19,
"execution_count": 20,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"# Specify multi-group cross section types to compute\n",
"mgxs_lib.mgxs_types = ['transport', 'nu-fission', 'fission', 'nu-scatter matrix', 'chi']"
"mgxs_lib.mgxs_types = ['transport', 'nu-fission', 'fission', 'nu-scatter matrix', 'chi', 'chi-delayed']"
]
},
{
@ -578,7 +578,7 @@
},
{
"cell_type": "code",
"execution_count": 20,
"execution_count": 21,
"metadata": {
"collapsed": true
},
@ -600,7 +600,7 @@
},
{
"cell_type": "code",
"execution_count": 21,
"execution_count": 22,
"metadata": {
"collapsed": true
},
@ -619,7 +619,7 @@
},
{
"cell_type": "code",
"execution_count": 22,
"execution_count": 23,
"metadata": {
"collapsed": true
},
@ -640,7 +640,7 @@
},
{
"cell_type": "code",
"execution_count": 23,
"execution_count": 24,
"metadata": {
"collapsed": true
},
@ -660,7 +660,7 @@
},
{
"cell_type": "code",
"execution_count": 24,
"execution_count": 25,
"metadata": {
"collapsed": false
},
@ -688,7 +688,7 @@
},
{
"cell_type": "code",
"execution_count": 25,
"execution_count": 26,
"metadata": {
"collapsed": true
},
@ -700,7 +700,7 @@
},
{
"cell_type": "code",
"execution_count": 26,
"execution_count": 27,
"metadata": {
"collapsed": false
},
@ -725,8 +725,8 @@
" Copyright: 2011-2016 Massachusetts Institute of Technology\n",
" License: http://openmc.readthedocs.io/en/latest/license.html\n",
" Version: 0.7.1\n",
" Git SHA1: 47ef320ad517612376e181ec6a6bc42ca0db98ce\n",
" Date/Time: 2016-05-14 12:29:07\n",
" Git SHA1: d30d010aea2c1cba90c993a9c501be9d5d81921d\n",
" Date/Time: 2016-06-17 12:58:23\n",
" MPI Processes: 1\n",
"\n",
" ===========================================================================\n",
@ -778,32 +778,32 @@
" 22/1 1.04175 1.02516 +/- 0.00588\n",
" 23/1 1.01909 1.02469 +/- 0.00543\n",
" 24/1 1.07119 1.02801 +/- 0.00603\n",
" 25/1 0.97445 1.02444 +/- 0.00665\n",
" 26/1 1.04737 1.02588 +/- 0.00638\n",
" 27/1 1.04656 1.02709 +/- 0.00612\n",
" 28/1 1.03464 1.02751 +/- 0.00578\n",
" 29/1 1.02528 1.02739 +/- 0.00547\n",
" 30/1 1.02799 1.02742 +/- 0.00519\n",
" 31/1 1.05846 1.02890 +/- 0.00516\n",
" 32/1 1.03811 1.02932 +/- 0.00493\n",
" 33/1 1.00894 1.02843 +/- 0.00480\n",
" 34/1 1.02049 1.02810 +/- 0.00460\n",
" 35/1 1.00690 1.02726 +/- 0.00450\n",
" 36/1 1.03129 1.02741 +/- 0.00432\n",
" 37/1 0.98864 1.02597 +/- 0.00440\n",
" 38/1 1.00017 1.02505 +/- 0.00434\n",
" 39/1 1.03635 1.02544 +/- 0.00421\n",
" 40/1 1.07090 1.02696 +/- 0.00434\n",
" 41/1 1.03141 1.02710 +/- 0.00420\n",
" 42/1 1.02624 1.02707 +/- 0.00406\n",
" 43/1 1.02668 1.02706 +/- 0.00394\n",
" 44/1 1.05940 1.02801 +/- 0.00394\n",
" 45/1 1.01149 1.02754 +/- 0.00385\n",
" 46/1 1.06958 1.02871 +/- 0.00392\n",
" 47/1 1.02674 1.02866 +/- 0.00381\n",
" 48/1 1.02542 1.02857 +/- 0.00371\n",
" 49/1 1.03516 1.02874 +/- 0.00362\n",
" 50/1 1.06818 1.02973 +/- 0.00366\n",
" 25/1 0.97414 1.02442 +/- 0.00666\n",
" 26/1 1.04709 1.02584 +/- 0.00639\n",
" 27/1 1.05872 1.02777 +/- 0.00631\n",
" 28/1 1.03930 1.02841 +/- 0.00598\n",
" 29/1 1.01488 1.02770 +/- 0.00570\n",
" 30/1 1.04513 1.02857 +/- 0.00548\n",
" 31/1 0.99538 1.02699 +/- 0.00545\n",
" 32/1 1.00106 1.02581 +/- 0.00532\n",
" 33/1 0.99389 1.02442 +/- 0.00527\n",
" 34/1 0.99938 1.02338 +/- 0.00516\n",
" 35/1 1.02161 1.02331 +/- 0.00495\n",
" 36/1 1.04084 1.02398 +/- 0.00480\n",
" 37/1 0.98801 1.02265 +/- 0.00481\n",
" 38/1 1.01348 1.02232 +/- 0.00464\n",
" 39/1 1.06693 1.02386 +/- 0.00474\n",
" 40/1 1.07729 1.02564 +/- 0.00491\n",
" 41/1 1.03191 1.02585 +/- 0.00475\n",
" 42/1 1.05209 1.02667 +/- 0.00468\n",
" 43/1 1.02997 1.02677 +/- 0.00453\n",
" 44/1 1.07288 1.02812 +/- 0.00460\n",
" 45/1 1.01268 1.02768 +/- 0.00449\n",
" 46/1 1.03759 1.02796 +/- 0.00437\n",
" 47/1 1.02620 1.02791 +/- 0.00425\n",
" 48/1 1.02509 1.02783 +/- 0.00414\n",
" 49/1 1.01043 1.02739 +/- 0.00406\n",
" 50/1 1.01457 1.02707 +/- 0.00397\n",
" Creating state point statepoint.50.h5...\n",
"\n",
" ===========================================================================\n",
@ -813,27 +813,27 @@
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 5.7700E-01 seconds\n",
" Reading cross sections = 1.3400E-01 seconds\n",
" Total time in simulation = 8.0461E+01 seconds\n",
" Time in transport only = 8.0422E+01 seconds\n",
" Time in inactive batches = 6.4060E+00 seconds\n",
" Time in active batches = 7.4055E+01 seconds\n",
" Time synchronizing fission bank = 6.0000E-03 seconds\n",
" Sampling source sites = 2.0000E-03 seconds\n",
" SEND/RECV source sites = 3.0000E-03 seconds\n",
" Time accumulating tallies = 0.0000E+00 seconds\n",
" Total time for initialization = 8.4000E-01 seconds\n",
" Reading cross sections = 2.5400E-01 seconds\n",
" Total time in simulation = 7.0993E+01 seconds\n",
" Time in transport only = 7.0952E+01 seconds\n",
" Time in inactive batches = 5.0170E+00 seconds\n",
" Time in active batches = 6.5976E+01 seconds\n",
" Time synchronizing fission bank = 3.0000E-03 seconds\n",
" Sampling source sites = 1.0000E-03 seconds\n",
" SEND/RECV source sites = 1.0000E-03 seconds\n",
" Time accumulating tallies = 1.0000E-03 seconds\n",
" Total time for finalization = 0.0000E+00 seconds\n",
" Total time elapsed = 8.1067E+01 seconds\n",
" Calculation Rate (inactive) = 3902.59 neutrons/second\n",
" Calculation Rate (active) = 1350.35 neutrons/second\n",
" Total time elapsed = 7.1853E+01 seconds\n",
" Calculation Rate (inactive) = 4983.06 neutrons/second\n",
" Calculation Rate (active) = 1515.70 neutrons/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" k-effective (Collision) = 1.02763 +/- 0.00343\n",
" k-effective (Track-length) = 1.02973 +/- 0.00366\n",
" k-effective (Absorption) = 1.02732 +/- 0.00319\n",
" Combined k-effective = 1.02826 +/- 0.00259\n",
" k-effective (Collision) = 1.02489 +/- 0.00308\n",
" k-effective (Track-length) = 1.02707 +/- 0.00397\n",
" k-effective (Absorption) = 1.02637 +/- 0.00325\n",
" Combined k-effective = 1.02581 +/- 0.00264\n",
" Leakage Fraction = 0.00000 +/- 0.00000\n",
"\n"
]
@ -844,7 +844,7 @@
"0"
]
},
"execution_count": 26,
"execution_count": 27,
"metadata": {},
"output_type": "execute_result"
}
@ -870,7 +870,7 @@
},
{
"cell_type": "code",
"execution_count": 27,
"execution_count": 28,
"metadata": {
"collapsed": false
},
@ -889,7 +889,7 @@
},
{
"cell_type": "code",
"execution_count": 28,
"execution_count": 29,
"metadata": {
"collapsed": false
},
@ -924,7 +924,7 @@
},
{
"cell_type": "code",
"execution_count": 29,
"execution_count": 44,
"metadata": {
"collapsed": false
},
@ -944,19 +944,11 @@
},
{
"cell_type": "code",
"execution_count": 30,
"execution_count": 45,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/home/wboyd/Documents/NSE-CRPG-Codes/openmc/openmc/tallies.py:1988: RuntimeWarning: invalid value encountered in true_divide\n",
" self_rel_err = data['self']['std. dev.'] / data['self']['mean']\n"
]
},
{
"data": {
"text/html": [
@ -978,16 +970,16 @@
" <td>10000</td>\n",
" <td>1</td>\n",
" <td>U-235</td>\n",
" <td>8.055246e-03</td>\n",
" <td>2.857567e-05</td>\n",
" <td>8.046809e-03</td>\n",
" <td>2.697198e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>4</th>\n",
" <td>10000</td>\n",
" <td>1</td>\n",
" <td>U-238</td>\n",
" <td>7.339215e-03</td>\n",
" <td>4.349466e-05</td>\n",
" <td>7.366624e-03</td>\n",
" <td>4.255197e-05</td>\n",
" </tr>\n",
" <tr>\n",
" <th>5</th>\n",
@ -1002,16 +994,16 @@
" <td>10000</td>\n",
" <td>2</td>\n",
" <td>U-235</td>\n",
" <td>3.615565e-01</td>\n",
" <td>2.050486e-03</td>\n",
" <td>3.614917e-01</td>\n",
" <td>2.135233e-03</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>10000</td>\n",
" <td>2</td>\n",
" <td>U-238</td>\n",
" <td>6.742638e-07</td>\n",
" <td>3.795256e-09</td>\n",
" <td>6.741607e-07</td>\n",
" <td>3.924924e-09</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
@ -1027,15 +1019,15 @@
],
"text/plain": [
" cell group in nuclide mean std. dev.\n",
"3 10000 1 U-235 8.055246e-03 2.857567e-05\n",
"4 10000 1 U-238 7.339215e-03 4.349466e-05\n",
"3 10000 1 U-235 8.046809e-03 2.697198e-05\n",
"4 10000 1 U-238 7.366624e-03 4.255197e-05\n",
"5 10000 1 O-16 0.000000e+00 0.000000e+00\n",
"0 10000 2 U-235 3.615565e-01 2.050486e-03\n",
"1 10000 2 U-238 6.742638e-07 3.795256e-09\n",
"0 10000 2 U-235 3.614917e-01 2.135233e-03\n",
"1 10000 2 U-238 6.741607e-07 3.924924e-09\n",
"2 10000 2 O-16 0.000000e+00 0.000000e+00"
]
},
"execution_count": 30,
"execution_count": 45,
"metadata": {},
"output_type": "execute_result"
}
@ -1054,7 +1046,7 @@
},
{
"cell_type": "code",
"execution_count": 31,
"execution_count": 32,
"metadata": {
"collapsed": false
},
@ -1069,13 +1061,13 @@
"\tDomain ID =\t10000\n",
"\tNuclide =\tU-235\n",
"\tCross Sections [cm^-1]:\n",
" Group 1 [6.25e-07 - 20.0 MeV]:\t8.06e-03 +/- 3.55e-01%\n",
" Group 2 [0.0 - 6.25e-07 MeV]:\t3.62e-01 +/- 5.67e-01%\n",
" Group 1 [6.25e-07 - 20.0 MeV]:\t8.05e-03 +/- 3.35e-01%\n",
" Group 2 [0.0 - 6.25e-07 MeV]:\t3.61e-01 +/- 5.91e-01%\n",
"\n",
"\tNuclide =\tU-238\n",
"\tCross Sections [cm^-1]:\n",
" Group 1 [6.25e-07 - 20.0 MeV]:\t7.34e-03 +/- 5.93e-01%\n",
" Group 2 [0.0 - 6.25e-07 MeV]:\t6.74e-07 +/- 5.63e-01%\n",
" Group 1 [6.25e-07 - 20.0 MeV]:\t7.37e-03 +/- 5.78e-01%\n",
" Group 2 [0.0 - 6.25e-07 MeV]:\t6.74e-07 +/- 5.82e-01%\n",
"\n",
"\tNuclide =\tO-16\n",
"\tCross Sections [cm^-1]:\n",
@ -1100,7 +1092,7 @@
},
{
"cell_type": "code",
"execution_count": 32,
"execution_count": 33,
"metadata": {
"collapsed": false
},
@ -1119,7 +1111,7 @@
},
{
"cell_type": "code",
"execution_count": 33,
"execution_count": 34,
"metadata": {
"collapsed": true
},
@ -1131,7 +1123,7 @@
},
{
"cell_type": "code",
"execution_count": 34,
"execution_count": 35,
"metadata": {
"collapsed": true
},
@ -1150,7 +1142,7 @@
},
{
"cell_type": "code",
"execution_count": 35,
"execution_count": 36,
"metadata": {
"collapsed": true
},
@ -1165,7 +1157,7 @@
},
{
"cell_type": "code",
"execution_count": 36,
"execution_count": 41,
"metadata": {
"collapsed": false
},
@ -1191,16 +1183,16 @@
" <td>10000</td>\n",
" <td>1</td>\n",
" <td>U-235</td>\n",
" <td>0.074860</td>\n",
" <td>0.000303</td>\n",
" <td>0.074734</td>\n",
" <td>0.000325</td>\n",
" </tr>\n",
" <tr>\n",
" <th>1</th>\n",
" <td>10000</td>\n",
" <td>1</td>\n",
" <td>U-238</td>\n",
" <td>0.005952</td>\n",
" <td>0.000035</td>\n",
" <td>0.005977</td>\n",
" <td>0.000034</td>\n",
" </tr>\n",
" <tr>\n",
" <th>2</th>\n",
@ -1216,12 +1208,12 @@
],
"text/plain": [
" cell group in nuclide mean std. dev.\n",
"0 10000 1 U-235 0.074860 0.000303\n",
"1 10000 1 U-238 0.005952 0.000035\n",
"0 10000 1 U-235 0.074734 0.000325\n",
"1 10000 1 U-238 0.005977 0.000034\n",
"2 10000 1 O-16 0.000000 0.000000"
]
},
"execution_count": 36,
"execution_count": 41,
"metadata": {},
"output_type": "execute_result"
}
@ -1250,11 +1242,28 @@
},
{
"cell_type": "code",
"execution_count": 37,
"execution_count": 38,
"metadata": {
"collapsed": false
},
"outputs": [],
"outputs": [
{
"ename": "AttributeError",
"evalue": "'RectLattice' object has no attribute 'dimension'",
"output_type": "error",
"traceback": [
"\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
"\u001b[0;31mAttributeError\u001b[0m Traceback (most recent call last)",
"\u001b[0;32m<ipython-input-38-c2821928e70b>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m 1\u001b[0m \u001b[0;31m# Create an OpenMOC Geometry from the OpenCG Geometry\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 2\u001b[0;31m \u001b[0mopenmoc_geometry\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mget_openmoc_geometry\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mmgxs_lib\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mopencg_geometry\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
"\u001b[0;32m/Users/sam/.local/lib/python2.7/site-packages/openmc-0.7.1-py2.7.egg/openmc/mgxs/library.pyc\u001b[0m in \u001b[0;36mopencg_geometry\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m 161\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_opencg_geometry\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 162\u001b[0m \u001b[0;32mfrom\u001b[0m \u001b[0mopenmc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mopencg_compatible\u001b[0m \u001b[0;32mimport\u001b[0m \u001b[0mget_opencg_geometry\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 163\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_opencg_geometry\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mget_opencg_geometry\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_openmc_geometry\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 164\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_opencg_geometry\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 165\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/sam/.local/lib/python2.7/site-packages/openmc-0.7.1-py2.7.egg/openmc/opencg_compatible.pyc\u001b[0m in \u001b[0;36mget_opencg_geometry\u001b[0;34m(openmc_geometry)\u001b[0m\n\u001b[1;32m 997\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 998\u001b[0m \u001b[0mopenmc_root_universe\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mopenmc_geometry\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mroot_universe\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 999\u001b[0;31m \u001b[0mopencg_root_universe\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mget_opencg_universe\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mopenmc_root_universe\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 1000\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 1001\u001b[0m \u001b[0mopencg_geometry\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mopencg\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mGeometry\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/sam/.local/lib/python2.7/site-packages/openmc-0.7.1-py2.7.egg/openmc/opencg_compatible.pyc\u001b[0m in \u001b[0;36mget_opencg_universe\u001b[0;34m(openmc_universe)\u001b[0m\n\u001b[1;32m 746\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 747\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mcell_id\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mopenmc_cell\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mopenmc_cells\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mitems\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 748\u001b[0;31m \u001b[0mopencg_cell\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mget_opencg_cell\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mopenmc_cell\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 749\u001b[0m \u001b[0mopencg_universe\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0madd_cell\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mopencg_cell\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 750\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/sam/.local/lib/python2.7/site-packages/openmc-0.7.1-py2.7.egg/openmc/opencg_compatible.pyc\u001b[0m in \u001b[0;36mget_opencg_cell\u001b[0;34m(openmc_cell)\u001b[0m\n\u001b[1;32m 454\u001b[0m \u001b[0mopencg_cell\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfill\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mget_opencg_universe\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfill\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 455\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 456\u001b[0;31m \u001b[0mopencg_cell\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfill\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mget_opencg_lattice\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfill\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 457\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 458\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mopenmc_cell\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrotation\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/sam/.local/lib/python2.7/site-packages/openmc-0.7.1-py2.7.egg/openmc/opencg_compatible.pyc\u001b[0m in \u001b[0;36mget_opencg_lattice\u001b[0;34m(openmc_lattice)\u001b[0m\n\u001b[1;32m 831\u001b[0m \u001b[0;31m# Create an OpenCG Lattice to represent this OpenMC Lattice\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 832\u001b[0m \u001b[0mname\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mopenmc_lattice\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 833\u001b[0;31m \u001b[0mdimension\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mopenmc_lattice\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdimension\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 834\u001b[0m \u001b[0mpitch\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mopenmc_lattice\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpitch\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 835\u001b[0m \u001b[0mlower_left\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mopenmc_lattice\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlower_left\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;31mAttributeError\u001b[0m: 'RectLattice' object has no attribute 'dimension'"
]
}
],
"source": [
"# Create an OpenMOC Geometry from the OpenCG Geometry\n",
"openmoc_geometry = get_openmoc_geometry(mgxs_lib.opencg_geometry)"
@ -1269,11 +1278,23 @@
},
{
"cell_type": "code",
"execution_count": 38,
"execution_count": 39,
"metadata": {
"collapsed": false
},
"outputs": [],
"outputs": [
{
"ename": "NameError",
"evalue": "name 'openmoc_geometry' is not defined",
"output_type": "error",
"traceback": [
"\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
"\u001b[0;31mNameError\u001b[0m Traceback (most recent call last)",
"\u001b[0;32m<ipython-input-39-feb6ca9af456>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m 1\u001b[0m \u001b[0;31m# Load the library into the OpenMOC geometry\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 2\u001b[0;31m \u001b[0mmaterials\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mload_openmc_mgxs_lib\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mmgxs_lib\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mopenmoc_geometry\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
"\u001b[0;31mNameError\u001b[0m: name 'openmoc_geometry' is not defined"
]
}
],
"source": [
"# Load the library into the OpenMOC geometry\n",
"materials = load_openmc_mgxs_lib(mgxs_lib, openmoc_geometry)"
@ -1596,21 +1617,21 @@
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"display_name": "Python 2",
"language": "python",
"name": "python3"
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.5.1"
"pygments_lexer": "ipython2",
"version": "2.7.11"
}
},
"nbformat": 4,

10
docs/source/pythonapi/examples/mgxs-part-iv.ipynb
View file

@ -1439,21 +1439,21 @@
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"display_name": "Python 2",
"language": "python",
"name": "python3"
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.5.1"
"pygments_lexer": "ipython2",
"version": "2.7.11"
}
},
"nbformat": 4,

1
examples/xml/basic/materials.xml
View file

@ -12,7 +12,6 @@
<density value="1.0" units="g/cc" />
<nuclide name="H-1" ao="2.0" />
<nuclide name="O-16" ao="1.0" />
<sab name="HH2O" xs="71t" />
</material>
</materials>

2
openmc/__init__.py
View file

@ -9,8 +9,8 @@ from openmc.plots import *
from openmc.settings import *
from openmc.surface import *
from openmc.universe import *
from openmc.mgxs_library import *
from openmc.mesh import *
from openmc.mgxs_library import *
from openmc.filter import *
from openmc.trigger import *
from openmc.tallies import *

2
openmc/filter.py
View file

@ -564,7 +564,7 @@ class Filter(object):
# Initialize dictionary to build Pandas Multi-index column
filter_dict = {}
# Append Mesh ID as outermost index of mult-index
# Append Mesh ID as outermost index of multi-index
mesh_key = 'mesh {0}'.format(self.mesh.id)
# Find mesh dimensions - use 3D indices for simplicity

38
openmc/mgxs/library.py
View file

@ -55,9 +55,10 @@ class Library(object):
If true, computes cross sections for each nuclide in each domain
mgxs_types : Iterable of str
The types of cross sections in the library (e.g., ['total', 'scatter'])
domain_type : {'material', 'cell', 'distribcell', 'universe'}
domain_type : {'material', 'cell', 'distribcell', 'universe', 'mesh'}
Domain type for spatial homogenization
domains : Iterable of openmc.Material, openmc.Cell or openmc.Universe
domains : Iterable of openmc.Material, openmc.Cell, openmc.Universe, or
openmc.Mesh
The spatial domain(s) for which MGXS in the Library are computed
correction : {'P0', None}
Apply the P0 correction to scattering matrices if set to 'P0'
@ -188,6 +189,10 @@ class Library(object):
return self.openmc_geometry.get_all_material_cells()
elif self.domain_type == 'universe':
return self.openmc_geometry.get_all_universes()
# FIXME: Change to get tuples of all domain cells
elif self.domain_type == 'mesh':
raise ValueError('Unable to get all domains for a mesh domain ' +
'type. The domains must be set to [openmc.Mesh]')
else:
raise ValueError('Unable to get domains without a domain type')
else:
@ -253,11 +258,21 @@ class Library(object):
@by_nuclide.setter
def by_nuclide(self, by_nuclide):
cv.check_type('by_nuclide', by_nuclide, bool)
if by_nuclide == True and self.domain_type == 'mesh':
raise ValueError('Unable to create MGXS library by nuclide with ' +
'mesh domain')
self._by_nuclide = by_nuclide
@domain_type.setter
def domain_type(self, domain_type):
cv.check_value('domain type', domain_type, openmc.mgxs.DOMAIN_TYPES)
if by_nuclide == True and domain_type == 'mesh':
raise ValueError('Unable to create MGXS library by nuclide with ' +
'mesh domain')
self._domain_type = domain_type
@domains.setter
@ -278,6 +293,9 @@ class Library(object):
elif self.domain_type == 'universe':
cv.check_iterable_type('domain', domains, openmc.Universe)
all_domains = self.openmc_geometry.get_all_universes()
elif self.domain_type == 'mesh':
cv.check_iterable_type('domain', domains, openmc.Mesh)
all_domains = domains
else:
msg = 'Unable to set domains with ' \
'domain type "{}"'.format(self.domain_type)
@ -458,9 +476,13 @@ class Library(object):
Parameters
----------
domain : Material or Cell or Universe or Integral
The material, cell, or universe object of interest (or its ID)
mgxs_type : {'total', 'transport', 'nu-transport', 'absorption', 'capture', 'fission', 'nu-fission', 'kappa-fission', 'scatter', 'nu-scatter', 'scatter matrix', 'nu-scatter matrix', 'multiplicity matrix', 'nu-fission matrix', chi'}
domain : Material or Cell or Universe or Mesh or Integral
The material, cell, universe, or mesh object of interest (or its ID)
mgxs_type : {'total', 'transport', 'nu-transport', 'absorption',
'capture', 'fission', 'nu-fission', 'kappa-fission',
'scatter', 'nu-scatter', 'scatter matrix',
'nu-scatter matrix', 'multiplicity matrix',
'nu-fission matrix', chi'}
The type of multi-group cross section object to return
Returns
@ -482,6 +504,8 @@ class Library(object):
cv.check_type('domain', domain, (openmc.Cell, Integral))
elif self.domain_type == 'universe':
cv.check_type('domain', domain, (openmc.Universe, Integral))
elif self.domain_type == 'mesh':
cv.check_type('domain', domain, (openmc.Mesh, Integral))
# Check that requested domain is included in library
if isinstance(domain, Integral):
@ -761,7 +785,7 @@ class Library(object):
Parameters
----------
domain : openmc.Material or openmc.Cell or openmc.Universe
domain : openmc.Material or openmc.Cell or openmc.Universe or openmc.Mesh
The domain for spatial homogenization
xsdata_name : str
Name to apply to the "xsdata" entry produced by this method
@ -809,7 +833,7 @@ class Library(object):
"""
cv.check_type('domain', domain, (openmc.Material, openmc.Cell,
openmc.Cell))
openmc.Cell, openmc.Mesh))
cv.check_type('xsdata_name', xsdata_name, basestring)
cv.check_type('nuclide', nuclide, basestring)
cv.check_value('xs_type', xs_type, ['macro', 'micro'])

703
openmc/mgxs/mgxs.py
View file

File diff suppressed because it is too large Load diff