Improve Pandas float handling

This commit is contained in:
Sterling Harper 2016-02-07 16:06:33 -05:00
parent 5501d1b42a
commit 5dfc9e91fc
6 changed files with 549 additions and 589 deletions

View file

@ -519,7 +519,7 @@
" License: http://mit-crpg.github.io/openmc/license.html\n",
" Version: 0.7.1\n",
" Git SHA1: 34381b40a9445a727e360873aaa6ef892af1cb6a\n",
" Date/Time: 2016-02-07 14:10:52\n",
" Date/Time: 2016-02-07 15:58:16\n",
" MPI Processes: 1\n",
"\n",
" ===========================================================================\n",
@ -605,20 +605,20 @@
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 8.7800E-01 seconds\n",
" Reading cross sections = 1.9200E-01 seconds\n",
" Total time in simulation = 1.3677E+01 seconds\n",
" Time in transport only = 1.3579E+01 seconds\n",
" Time in inactive batches = 1.7900E+00 seconds\n",
" Time in active batches = 1.1887E+01 seconds\n",
" Time synchronizing fission bank = 6.0000E-03 seconds\n",
" Sampling source sites = 2.0000E-03 seconds\n",
" Total time for initialization = 3.2100E-01 seconds\n",
" Reading cross sections = 7.4000E-02 seconds\n",
" Total time in simulation = 8.3830E+00 seconds\n",
" Time in transport only = 8.3670E+00 seconds\n",
" Time in inactive batches = 1.0330E+00 seconds\n",
" Time in active batches = 7.3500E+00 seconds\n",
" Time synchronizing fission bank = 4.0000E-03 seconds\n",
" Sampling source sites = 1.0000E-03 seconds\n",
" SEND/RECV source sites = 3.0000E-03 seconds\n",
" Time accumulating tallies = 1.0000E-03 seconds\n",
" Time accumulating tallies = 0.0000E+00 seconds\n",
" Total time for finalization = 1.0000E-03 seconds\n",
" Total time elapsed = 1.4578E+01 seconds\n",
" Calculation Rate (inactive) = 13966.5 neutrons/second\n",
" Calculation Rate (active) = 8412.55 neutrons/second\n",
" Total time elapsed = 8.7140E+00 seconds\n",
" Calculation Rate (inactive) = 24201.4 neutrons/second\n",
" Calculation Rate (active) = 13605.4 neutrons/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
@ -909,8 +909,8 @@
" <tr>\n",
" <th>0</th>\n",
" <td> 1</td>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0.000000</td>\n",
" <td> 0.000001</td>\n",
" <td> total</td>\n",
" <td> (((total / flux) - (absorption / flux)) - (sca...</td>\n",
" <td> 4.884981e-15</td>\n",
@ -919,8 +919,8 @@
" <tr>\n",
" <th>1</th>\n",
" <td> 1</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 2.00e+01</td>\n",
" <td> 0.000001</td>\n",
" <td> 20.000000</td>\n",
" <td> total</td>\n",
" <td> (((total / flux) - (absorption / flux)) - (sca...</td>\n",
" <td> 1.221245e-15</td>\n",
@ -931,13 +931,13 @@
"</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",
" 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",
"\n",
" score mean std. dev. \n",
"0 (((total / flux) - (absorption / flux)) - (sca... 4.884981e-15 0.011274 \n",
"1 (((total / flux) - (absorption / flux)) - (sca... 1.221245e-15 0.001802 "
" score mean std. dev. \n",
"0 (((total / flux) - (absorption / flux)) - (sca... 4.88e-15 1.13e-02 \n",
"1 (((total / flux) - (absorption / flux)) - (sca... 1.22e-15 1.80e-03 "
]
},
"execution_count": 23,
@ -988,8 +988,8 @@
" <tr>\n",
" <th>0</th>\n",
" <td> 1</td>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0.000000</td>\n",
" <td> 0.000001</td>\n",
" <td> total</td>\n",
" <td> ((absorption / flux) / (total / flux))</td>\n",
" <td> 0.076219</td>\n",
@ -998,8 +998,8 @@
" <tr>\n",
" <th>1</th>\n",
" <td> 1</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 2.00e+01</td>\n",
" <td> 0.000001</td>\n",
" <td> 20.000000</td>\n",
" <td> total</td>\n",
" <td> ((absorption / flux) / (total / flux))</td>\n",
" <td> 0.019319</td>\n",
@ -1010,13 +1010,13 @@
"</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",
" 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",
"\n",
" score mean std. dev. \n",
"0 ((absorption / flux) / (total / flux)) 0.076219 0.000651 \n",
"1 ((absorption / flux) / (total / flux)) 0.019319 0.000086 "
" score mean std. dev. \n",
"0 ((absorption / flux) / (total / flux)) 7.62e-02 6.51e-04 \n",
"1 ((absorption / flux) / (total / flux)) 1.93e-02 8.65e-05 "
]
},
"execution_count": 24,
@ -1060,8 +1060,8 @@
" <tr>\n",
" <th>0</th>\n",
" <td> 1</td>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0.000000</td>\n",
" <td> 0.000001</td>\n",
" <td> total</td>\n",
" <td> ((scatter / flux) / (total / flux))</td>\n",
" <td> 0.923781</td>\n",
@ -1070,8 +1070,8 @@
" <tr>\n",
" <th>1</th>\n",
" <td> 1</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 2.00e+01</td>\n",
" <td> 0.000001</td>\n",
" <td> 20.000000</td>\n",
" <td> total</td>\n",
" <td> ((scatter / flux) / (total / flux))</td>\n",
" <td> 0.980681</td>\n",
@ -1082,13 +1082,13 @@
"</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",
" 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",
"\n",
" score mean std. dev. \n",
"0 ((scatter / flux) / (total / flux)) 0.923781 0.007714 \n",
"1 ((scatter / flux) / (total / flux)) 0.980681 0.002617 "
" score mean std. dev. \n",
"0 ((scatter / flux) / (total / flux)) 9.24e-01 7.71e-03 \n",
"1 ((scatter / flux) / (total / flux)) 9.81e-01 2.62e-03 "
]
},
"execution_count": 25,
@ -1139,8 +1139,8 @@
" <tr>\n",
" <th>0</th>\n",
" <td> 1</td>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0.000000</td>\n",
" <td> 0.000001</td>\n",
" <td> total</td>\n",
" <td> (((absorption / flux) / (total / flux)) + ((sc...</td>\n",
" <td> 1</td>\n",
@ -1149,8 +1149,8 @@
" <tr>\n",
" <th>1</th>\n",
" <td> 1</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 2.00e+01</td>\n",
" <td> 0.000001</td>\n",
" <td> 20.000000</td>\n",
" <td> total</td>\n",
" <td> (((absorption / flux) / (total / flux)) + ((sc...</td>\n",
" <td> 1</td>\n",
@ -1161,13 +1161,13 @@
"</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",
" 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",
"\n",
" score mean std. dev. \n",
"0 (((absorption / flux) / (total / flux)) + ((sc... 1 0.007741 \n",
"1 (((absorption / flux) / (total / flux)) + ((sc... 1 0.002619 "
" score mean std. dev. \n",
"0 (((absorption / flux) / (total / flux)) + ((sc... 1.00e+00 7.74e-03 \n",
"1 (((absorption / flux) / (total / flux)) + ((sc... 1.00e+00 2.62e-03 "
]
},
"execution_count": 26,
@ -1182,15 +1182,6 @@
"# 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": []
}
],
"metadata": {

File diff suppressed because one or more lines are too long

View file

@ -366,7 +366,7 @@
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAD6AgMAAAD1grKuAAAABGdBTUEAALGPC/xhBQAAACBjSFJN\nAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEX///9yEhLpgJFNv8Tq\nQYT7AAAAAWJLR0QAiAUdSAAAAAd0SU1FB+ACBxMQBoBGcLYAAALKSURBVGje7dpLcqQwDAbgHHE2\nYeEj+D4cwQucBUfo+3CEXoSp8OhuhF70T4qpKXmdr21LogK2Pj7A8QmNP+HDhw8fPnz48Kf6VH9G\n+66vy+je8k19jnf8C5dXIPv86ms56lPdjvaYbyodx3ze+XLE76cXFiD4zPji99z0/AJ4n1lfvJ6f\nnl0A6x+578efMSg1wPr172/jPO5yFXM+Ef78gdblM+WPHyguP//t1/g6pA0wfln+ho/fwgYYn19C\n/xwDvwHGc9OvC+hs37DTrwuwfWanXxdQTC9Mvyygs3wjTL8uwPJpn/tNDbSGz7T0SBEWw4vLXzbQ\n6b6RoveIoO6TvPxlA63qs7z8ZQPF9F+SH22vbX8OQKf5Rtv+EgDNJ3X58wZaxWd1+fMGiuFvir8b\nvjp8J/tGy/6jAmRvhW8fwL3vVT+o3grfPoB7r/IpALI3tz8FoJN84/NV873hB8UnM3xzANtf8nb4\ndwmg3grfFEDJO8JPE0i9Ff4pAYL3pI8mkHor/HMCeO9JH00g9SafEsh7T/ppARBvp48UwJnelT5S\nACd7O31TAlnvKx9SQCd7B58KgPO+8iMFuPWe9E8F8BveWX7bAjzX9y4//Jve+fhsH6Ctv7n8PTzj\nvY/v9gEOHz58+PBX+6v/f/wPvnd54f3j6venE/yl769Xv7+j3x/o98/V32/o9+fl389Xnx+g5x/o\n+Qt6/oOeP6HnX+j5G3z+h54/ouefV5/foufP6Pk3ev4On/+j9w/o/Qd6/4Le/6D3T/D9V67Y/ZsV\nQBq+s+8f0ftP+P41axXguP9NWgDuu/Cdfv+N3r/D9/9TAID+A7T/Ae2/gPs/0P4TtP8F7r9J3AIO\n9P+g/Udw/9Oygbf7r9D+L7j/DO1/Q/vv4P4/tP8Q7n9E+y/h/k+0/xTuf4X7b+H+X7T/+BPuf3aM\n8OHDhw8fPnz4w/4vzcvgeY10sY0AAAAldEVYdGRhdGU6Y3JlYXRlADIwMTYtMDItMDdUMTQ6MTY6\nMDYtMDU6MDCStPm5AAAAJXRFWHRkYXRlOm1vZGlmeQAyMDE2LTAyLTA3VDE0OjE2OjA2LTA1OjAw\n4+lBBQAAAABJRU5ErkJggg==\n",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAD6AgMAAAD1grKuAAAABGdBTUEAALGPC/xhBQAAACBjSFJN\nAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEX///9yEhLpgJFNv8Tq\nQYT7AAAAAWJLR0QAiAUdSAAAAAd0SU1FB+ACBxUFD8qiUrQAAALKSURBVGje7dpLcqQwDAbgHHE2\nYeEj+D4cwQucBUfo+3CEXoSp8OhuhF70T4qpKXmdr21LogK2Pj7A8QmNP+HDhw8fPnz48Kf6VH9G\n+66vy+je8k19jnf8C5dXIPv86ms56lPdjvaYbyodx3ze+XLE76cXFiD4zPji99z0/AJ4n1lfvJ6f\nnl0A6x+578efMSg1wPr172/jPO5yFXM+Ef78gdblM+WPHyguP//t1/g6pA0wfln+ho/fwgYYn19C\n/xwDvwHGc9OvC+hs37DTrwuwfWanXxdQTC9Mvyygs3wjTL8uwPJpn/tNDbSGz7T0SBEWw4vLXzbQ\n6b6RoveIoO6TvPxlA63qs7z8ZQPF9F+SH22vbX8OQKf5Rtv+EgDNJ3X58wZaxWd1+fMGiuFvir8b\nvjp8J/tGy/6jAmRvhW8fwL3vVT+o3grfPoB7r/IpALI3tz8FoJN84/NV873hB8UnM3xzANtf8nb4\ndwmg3grfFEDJO8JPE0i9Ff4pAYL3pI8mkHor/HMCeO9JH00g9SafEsh7T/ppARBvp48UwJnelT5S\nACd7O31TAlnvKx9SQCd7B58KgPO+8iMFuPWe9E8F8BveWX7bAjzX9y4//Jve+fhsH6Ctv7n8PTzj\nvY/v9gEOHz58+PBX+6v/f/wPvnd54f3j6venE/yl769Xv7+j3x/o98/V32/o9+fl389Xnx+g5x/o\n+Qt6/oOeP6HnX+j5G3z+h54/ouefV5/foufP6Pk3ev4On/+j9w/o/Qd6/4Le/6D3T/D9V67Y/ZsV\nQBq+s+8f0ftP+P41axXguP9NWgDuu/Cdfv+N3r/D9/9TAID+A7T/Ae2/gPs/0P4TtP8F7r9J3AIO\n9P+g/Udw/9Oygbf7r9D+L7j/DO1/Q/vv4P4/tP8Q7n9E+y/h/k+0/xTuf4X7b+H+X7T/+BPuf3aM\n8OHDhw8fPnz4w/4vzcvgeY10sY0AAAAldEVYdGRhdGU6Y3JlYXRlADIwMTYtMDItMDdUMTY6MDU6\nMTUtMDU6MDAlEzIyAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDE2LTAyLTA3VDE2OjA1OjE1LTA1OjAw\nVE6KjgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<IPython.core.display.Image object>"
]
@ -577,7 +577,7 @@
" License: http://mit-crpg.github.io/openmc/license.html\n",
" Version: 0.7.1\n",
" Git SHA1: 34381b40a9445a727e360873aaa6ef892af1cb6a\n",
" Date/Time: 2016-02-07 14:16:08\n",
" Date/Time: 2016-02-07 16:05:17\n",
" MPI Processes: 1\n",
"\n",
" ===========================================================================\n",
@ -634,20 +634,20 @@
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 4.0900E-01 seconds\n",
" Reading cross sections = 1.0100E-01 seconds\n",
" Total time in simulation = 7.8100E+00 seconds\n",
" Time in transport only = 7.7980E+00 seconds\n",
" Time in inactive batches = 1.3850E+00 seconds\n",
" Time in active batches = 6.4250E+00 seconds\n",
" Time synchronizing fission bank = 1.0000E-03 seconds\n",
" Sampling source sites = 0.0000E+00 seconds\n",
" SEND/RECV source sites = 1.0000E-03 seconds\n",
" Total time for initialization = 3.4700E-01 seconds\n",
" Reading cross sections = 9.1000E-02 seconds\n",
" Total time in simulation = 7.3920E+00 seconds\n",
" Time in transport only = 7.3820E+00 seconds\n",
" Time in inactive batches = 1.0930E+00 seconds\n",
" Time in active batches = 6.2990E+00 seconds\n",
" Time synchronizing fission bank = 2.0000E-03 seconds\n",
" Sampling source sites = 1.0000E-03 seconds\n",
" SEND/RECV source sites = 0.0000E+00 seconds\n",
" Time accumulating tallies = 0.0000E+00 seconds\n",
" Total time for finalization = 1.0000E-03 seconds\n",
" Total time elapsed = 8.2360E+00 seconds\n",
" Calculation Rate (inactive) = 9025.27 neutrons/second\n",
" Calculation Rate (active) = 5836.58 neutrons/second\n",
" Total time for finalization = 2.0000E-03 seconds\n",
" Total time elapsed = 7.7510E+00 seconds\n",
" Calculation Rate (inactive) = 11436.4 neutrons/second\n",
" Calculation Rate (active) = 5953.33 neutrons/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
@ -769,8 +769,8 @@
"</div>"
],
"text/plain": [
" nuclide score mean std. dev.\n",
"0 total (nu-fission / absorption) 1.040166 0.009069"
" nuclide score mean std. dev.\n",
"0 total (nu-fission / absorption) 1.04e+00 9.07e-03"
]
},
"execution_count": 26,
@ -821,8 +821,8 @@
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0</td>\n",
" <td> 0.000001</td>\n",
" <td> total</td>\n",
" <td> absorption</td>\n",
" <td> 0.694707</td>\n",
@ -833,8 +833,8 @@
"</div>"
],
"text/plain": [
" energy low [MeV] energy high [MeV] nuclide score mean std. dev.\n",
"0 0.00e+00 6.25e-07 total absorption 0.694707 0.006699"
" energy low [MeV] energy high [MeV] nuclide score mean std. dev.\n",
"0 0.00e+00 6.25e-07 total absorption 6.95e-01 6.70e-03"
]
},
"execution_count": 27,
@ -883,8 +883,8 @@
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0</td>\n",
" <td> 0.000001</td>\n",
" <td> total</td>\n",
" <td> nu-fission</td>\n",
" <td> 1.201216</td>\n",
@ -895,8 +895,8 @@
"</div>"
],
"text/plain": [
" energy low [MeV] energy high [MeV] nuclide score mean std. dev.\n",
"0 0.00e+00 6.25e-07 total nu-fission 1.201216 0.012288"
" energy low [MeV] energy high [MeV] nuclide score mean std. dev.\n",
"0 0.00e+00 6.25e-07 total nu-fission 1.20e+00 1.23e-02"
]
},
"execution_count": 28,
@ -947,8 +947,8 @@
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0</td>\n",
" <td> 0.000001</td>\n",
" <td> 10000</td>\n",
" <td> total</td>\n",
" <td> absorption</td>\n",
@ -960,11 +960,11 @@
"</div>"
],
"text/plain": [
" energy low [MeV] energy high [MeV] cell nuclide score mean \\\n",
"0 0.00e+00 6.25e-07 10000 total absorption 0.74925 \n",
" energy low [MeV] energy high [MeV] cell nuclide score mean \\\n",
"0 0.00e+00 6.25e-07 10000 total absorption 7.49e-01 \n",
"\n",
" std. dev. \n",
"0 0.008257 "
"0 8.26e-03 "
]
},
"execution_count": 29,
@ -1013,8 +1013,8 @@
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0</td>\n",
" <td> 0.000001</td>\n",
" <td> 10000</td>\n",
" <td> total</td>\n",
" <td> (nu-fission / absorption)</td>\n",
@ -1026,11 +1026,11 @@
"</div>"
],
"text/plain": [
" energy low [MeV] energy high [MeV] cell nuclide \\\n",
"0 0.00e+00 6.25e-07 10000 total \n",
" energy low [MeV] energy high [MeV] cell nuclide \\\n",
"0 0.00e+00 6.25e-07 10000 total \n",
"\n",
" score mean std. dev. \n",
"0 (nu-fission / absorption) 1.663616 0.018624 "
" score mean std. dev. \n",
"0 (nu-fission / absorption) 1.66e+00 1.86e-02 "
]
},
"execution_count": 30,
@ -1078,8 +1078,8 @@
" <tbody>\n",
" <tr>\n",
" <th>0</th>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0</td>\n",
" <td> 0.000001</td>\n",
" <td> 10000</td>\n",
" <td> total</td>\n",
" <td> (((absorption * nu-fission) * absorption) * (n...</td>\n",
@ -1091,11 +1091,11 @@
"</div>"
],
"text/plain": [
" energy low [MeV] energy high [MeV] cell nuclide \\\n",
"0 0.00e+00 6.25e-07 10000 total \n",
" energy low [MeV] energy high [MeV] cell nuclide \\\n",
"0 0.00e+00 6.25e-07 10000 total \n",
"\n",
" score mean std. dev. \n",
"0 (((absorption * nu-fission) * absorption) * (n... 1.040166 0.021928 "
" score mean std. dev. \n",
"0 (((absorption * nu-fission) * absorption) * (n... 1.04e+00 2.19e-02 "
]
},
"execution_count": 31,
@ -1161,8 +1161,8 @@
" <tr>\n",
" <th>0</th>\n",
" <td> 10000</td>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0.000000</td>\n",
" <td> 0.000001</td>\n",
" <td> (U-238 / total)</td>\n",
" <td> (nu-fission / flux)</td>\n",
" <td> 0.000001</td>\n",
@ -1171,8 +1171,8 @@
" <tr>\n",
" <th>1</th>\n",
" <td> 10000</td>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0.000000</td>\n",
" <td> 0.000001</td>\n",
" <td> (U-238 / total)</td>\n",
" <td> (scatter / flux)</td>\n",
" <td> 0.209989</td>\n",
@ -1181,8 +1181,8 @@
" <tr>\n",
" <th>2</th>\n",
" <td> 10000</td>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0.000000</td>\n",
" <td> 0.000001</td>\n",
" <td> (U-235 / total)</td>\n",
" <td> (nu-fission / flux)</td>\n",
" <td> 0.356420</td>\n",
@ -1191,8 +1191,8 @@
" <tr>\n",
" <th>3</th>\n",
" <td> 10000</td>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0.000000</td>\n",
" <td> 0.000001</td>\n",
" <td> (U-235 / total)</td>\n",
" <td> (scatter / flux)</td>\n",
" <td> 0.005555</td>\n",
@ -1201,8 +1201,8 @@
" <tr>\n",
" <th>4</th>\n",
" <td> 10000</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 2.00e+01</td>\n",
" <td> 0.000001</td>\n",
" <td> 20.000000</td>\n",
" <td> (U-238 / total)</td>\n",
" <td> (nu-fission / flux)</td>\n",
" <td> 0.007155</td>\n",
@ -1211,8 +1211,8 @@
" <tr>\n",
" <th>5</th>\n",
" <td> 10000</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 2.00e+01</td>\n",
" <td> 0.000001</td>\n",
" <td> 20.000000</td>\n",
" <td> (U-238 / total)</td>\n",
" <td> (scatter / flux)</td>\n",
" <td> 0.227770</td>\n",
@ -1221,8 +1221,8 @@
" <tr>\n",
" <th>6</th>\n",
" <td> 10000</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 2.00e+01</td>\n",
" <td> 0.000001</td>\n",
" <td> 20.000000</td>\n",
" <td> (U-235 / total)</td>\n",
" <td> (nu-fission / flux)</td>\n",
" <td> 0.008067</td>\n",
@ -1231,8 +1231,8 @@
" <tr>\n",
" <th>7</th>\n",
" <td> 10000</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 2.00e+01</td>\n",
" <td> 0.000001</td>\n",
" <td> 20.000000</td>\n",
" <td> (U-235 / total)</td>\n",
" <td> (scatter / flux)</td>\n",
" <td> 0.003367</td>\n",
@ -1243,25 +1243,25 @@
"</div>"
],
"text/plain": [
" cell energy low [MeV] energy high [MeV] nuclide \\\n",
"0 10000 0.00e+00 6.25e-07 (U-238 / total) \n",
"1 10000 0.00e+00 6.25e-07 (U-238 / total) \n",
"2 10000 0.00e+00 6.25e-07 (U-235 / total) \n",
"3 10000 0.00e+00 6.25e-07 (U-235 / total) \n",
"4 10000 6.25e-07 2.00e+01 (U-238 / total) \n",
"5 10000 6.25e-07 2.00e+01 (U-238 / total) \n",
"6 10000 6.25e-07 2.00e+01 (U-235 / total) \n",
"7 10000 6.25e-07 2.00e+01 (U-235 / total) \n",
" cell energy low [MeV] energy high [MeV] nuclide \\\n",
"0 10000 0.00e+00 6.25e-07 (U-238 / total) \n",
"1 10000 0.00e+00 6.25e-07 (U-238 / total) \n",
"2 10000 0.00e+00 6.25e-07 (U-235 / total) \n",
"3 10000 0.00e+00 6.25e-07 (U-235 / total) \n",
"4 10000 6.25e-07 2.00e+01 (U-238 / total) \n",
"5 10000 6.25e-07 2.00e+01 (U-238 / total) \n",
"6 10000 6.25e-07 2.00e+01 (U-235 / total) \n",
"7 10000 6.25e-07 2.00e+01 (U-235 / total) \n",
"\n",
" score mean std. dev. \n",
"0 (nu-fission / flux) 0.000001 7.377419e-09 \n",
"1 (scatter / flux) 0.209989 2.303838e-03 \n",
"2 (nu-fission / flux) 0.356420 3.951669e-03 \n",
"3 (scatter / flux) 0.005555 6.101004e-05 \n",
"4 (nu-fission / flux) 0.007155 8.053460e-05 \n",
"5 (scatter / flux) 0.227770 1.079289e-03 \n",
"6 (nu-fission / flux) 0.008067 5.254797e-05 \n",
"7 (scatter / flux) 0.003367 1.647058e-05 "
" score mean std. dev. \n",
"0 (nu-fission / flux) 6.66e-07 7.38e-09 \n",
"1 (scatter / flux) 2.10e-01 2.30e-03 \n",
"2 (nu-fission / flux) 3.56e-01 3.95e-03 \n",
"3 (scatter / flux) 5.56e-03 6.10e-05 \n",
"4 (nu-fission / flux) 7.15e-03 8.05e-05 \n",
"5 (scatter / flux) 2.28e-01 1.08e-03 \n",
"6 (nu-fission / flux) 8.07e-03 5.25e-05 \n",
"7 (scatter / flux) 3.37e-03 1.65e-05 "
]
},
"execution_count": 33,
@ -1396,8 +1396,8 @@
" <tr>\n",
" <th>0</th>\n",
" <td> 10000</td>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0.000000</td>\n",
" <td> 0.000001</td>\n",
" <td> U-238</td>\n",
" <td> nu-fission</td>\n",
" <td> 0.000002</td>\n",
@ -1406,8 +1406,8 @@
" <tr>\n",
" <th>1</th>\n",
" <td> 10000</td>\n",
" <td> 0.00e+00</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 0.000000</td>\n",
" <td> 0.000001</td>\n",
" <td> U-235</td>\n",
" <td> nu-fission</td>\n",
" <td> 0.868553</td>\n",
@ -1416,8 +1416,8 @@
" <tr>\n",
" <th>2</th>\n",
" <td> 10000</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 2.00e+01</td>\n",
" <td> 0.000001</td>\n",
" <td> 20.000000</td>\n",
" <td> U-238</td>\n",
" <td> nu-fission</td>\n",
" <td> 0.082149</td>\n",
@ -1426,8 +1426,8 @@
" <tr>\n",
" <th>3</th>\n",
" <td> 10000</td>\n",
" <td> 6.25e-07</td>\n",
" <td> 2.00e+01</td>\n",
" <td> 0.000001</td>\n",
" <td> 20.000000</td>\n",
" <td> U-235</td>\n",
" <td> nu-fission</td>\n",
" <td> 0.092618</td>\n",
@ -1438,17 +1438,17 @@
"</div>"
],
"text/plain": [
" cell energy low [MeV] energy high [MeV] nuclide score mean \\\n",
"0 10000 0.00e+00 6.25e-07 U-238 nu-fission 0.000002 \n",
"1 10000 0.00e+00 6.25e-07 U-235 nu-fission 0.868553 \n",
"2 10000 6.25e-07 2.00e+01 U-238 nu-fission 0.082149 \n",
"3 10000 6.25e-07 2.00e+01 U-235 nu-fission 0.092618 \n",
" cell energy low [MeV] energy high [MeV] nuclide score mean \\\n",
"0 10000 0.00e+00 6.25e-07 U-238 nu-fission 1.62e-06 \n",
"1 10000 0.00e+00 6.25e-07 U-235 nu-fission 8.69e-01 \n",
"2 10000 6.25e-07 2.00e+01 U-238 nu-fission 8.21e-02 \n",
"3 10000 6.25e-07 2.00e+01 U-235 nu-fission 9.26e-02 \n",
"\n",
" std. dev. \n",
"0 1.283958e-08 \n",
"1 6.880390e-03 \n",
"2 8.837250e-04 \n",
"3 5.195308e-04 "
" std. dev. \n",
"0 1.28e-08 \n",
"1 6.88e-03 \n",
"2 8.84e-04 \n",
"3 5.20e-04 "
]
},
"execution_count": 37,
@ -1490,8 +1490,8 @@
" <tr>\n",
" <th>0</th>\n",
" <td> 10002</td>\n",
" <td> 1.00e-08</td>\n",
" <td> 1.08e-07</td>\n",
" <td> 1.000000e-08</td>\n",
" <td> 0.000000</td>\n",
" <td> H-1</td>\n",
" <td> scatter</td>\n",
" <td> 4.619398</td>\n",
@ -1500,8 +1500,8 @@
" <tr>\n",
" <th>1</th>\n",
" <td> 10002</td>\n",
" <td> 1.08e-07</td>\n",
" <td> 1.17e-06</td>\n",
" <td> 1.080060e-07</td>\n",
" <td> 0.000001</td>\n",
" <td> H-1</td>\n",
" <td> scatter</td>\n",
" <td> 2.030757</td>\n",
@ -1510,8 +1510,8 @@
" <tr>\n",
" <th>2</th>\n",
" <td> 10002</td>\n",
" <td> 1.17e-06</td>\n",
" <td> 1.26e-05</td>\n",
" <td> 1.166529e-06</td>\n",
" <td> 0.000013</td>\n",
" <td> H-1</td>\n",
" <td> scatter</td>\n",
" <td> 1.658488</td>\n",
@ -1520,8 +1520,8 @@
" <tr>\n",
" <th>3</th>\n",
" <td> 10002</td>\n",
" <td> 1.26e-05</td>\n",
" <td> 1.36e-04</td>\n",
" <td> 1.259921e-05</td>\n",
" <td> 0.000136</td>\n",
" <td> H-1</td>\n",
" <td> scatter</td>\n",
" <td> 1.853002</td>\n",
@ -1530,8 +1530,8 @@
" <tr>\n",
" <th>4</th>\n",
" <td> 10002</td>\n",
" <td> 1.36e-04</td>\n",
" <td> 1.47e-03</td>\n",
" <td> 1.360790e-04</td>\n",
" <td> 0.001470</td>\n",
" <td> H-1</td>\n",
" <td> scatter</td>\n",
" <td> 2.050773</td>\n",
@ -1540,8 +1540,8 @@
" <tr>\n",
" <th>5</th>\n",
" <td> 10002</td>\n",
" <td> 1.47e-03</td>\n",
" <td> 1.59e-02</td>\n",
" <td> 1.469734e-03</td>\n",
" <td> 0.015874</td>\n",
" <td> H-1</td>\n",
" <td> scatter</td>\n",
" <td> 2.131759</td>\n",
@ -1550,8 +1550,8 @@
" <tr>\n",
" <th>6</th>\n",
" <td> 10002</td>\n",
" <td> 1.59e-02</td>\n",
" <td> 1.71e-01</td>\n",
" <td> 1.587401e-02</td>\n",
" <td> 0.171449</td>\n",
" <td> H-1</td>\n",
" <td> scatter</td>\n",
" <td> 2.213710</td>\n",
@ -1560,8 +1560,8 @@
" <tr>\n",
" <th>7</th>\n",
" <td> 10002</td>\n",
" <td> 1.71e-01</td>\n",
" <td> 1.85e+00</td>\n",
" <td> 1.714488e-01</td>\n",
" <td> 1.851749</td>\n",
" <td> H-1</td>\n",
" <td> scatter</td>\n",
" <td> 2.011925</td>\n",
@ -1570,8 +1570,8 @@
" <tr>\n",
" <th>8</th>\n",
" <td> 10002</td>\n",
" <td> 1.85e+00</td>\n",
" <td> 2.00e+01</td>\n",
" <td> 1.851749e+00</td>\n",
" <td> 20.000000</td>\n",
" <td> H-1</td>\n",
" <td> scatter</td>\n",
" <td> 0.371280</td>\n",
@ -1582,27 +1582,27 @@
"</div>"
],
"text/plain": [
" cell energy low [MeV] energy high [MeV] nuclide score mean \\\n",
"0 10002 1.00e-08 1.08e-07 H-1 scatter 4.619398 \n",
"1 10002 1.08e-07 1.17e-06 H-1 scatter 2.030757 \n",
"2 10002 1.17e-06 1.26e-05 H-1 scatter 1.658488 \n",
"3 10002 1.26e-05 1.36e-04 H-1 scatter 1.853002 \n",
"4 10002 1.36e-04 1.47e-03 H-1 scatter 2.050773 \n",
"5 10002 1.47e-03 1.59e-02 H-1 scatter 2.131759 \n",
"6 10002 1.59e-02 1.71e-01 H-1 scatter 2.213710 \n",
"7 10002 1.71e-01 1.85e+00 H-1 scatter 2.011925 \n",
"8 10002 1.85e+00 2.00e+01 H-1 scatter 0.371280 \n",
" cell energy low [MeV] energy high [MeV] nuclide score mean \\\n",
"0 10002 1.00e-08 1.08e-07 H-1 scatter 4.62e+00 \n",
"1 10002 1.08e-07 1.17e-06 H-1 scatter 2.03e+00 \n",
"2 10002 1.17e-06 1.26e-05 H-1 scatter 1.66e+00 \n",
"3 10002 1.26e-05 1.36e-04 H-1 scatter 1.85e+00 \n",
"4 10002 1.36e-04 1.47e-03 H-1 scatter 2.05e+00 \n",
"5 10002 1.47e-03 1.59e-02 H-1 scatter 2.13e+00 \n",
"6 10002 1.59e-02 1.71e-01 H-1 scatter 2.21e+00 \n",
"7 10002 1.71e-01 1.85e+00 H-1 scatter 2.01e+00 \n",
"8 10002 1.85e+00 2.00e+01 H-1 scatter 3.71e-01 \n",
"\n",
" std. dev. \n",
"0 0.040124 \n",
"1 0.011239 \n",
"2 0.009777 \n",
"3 0.007378 \n",
"4 0.012484 \n",
"5 0.007821 \n",
"6 0.015159 \n",
"7 0.009406 \n",
"8 0.003949 "
"0 4.01e-02 \n",
"1 1.12e-02 \n",
"2 9.78e-03 \n",
"3 7.38e-03 \n",
"4 1.25e-02 \n",
"5 7.82e-03 \n",
"6 1.52e-02 \n",
"7 9.41e-03 \n",
"8 3.95e-03 "
]
},
"execution_count": 38,

View file

@ -430,7 +430,7 @@ class CrossFilter(object):
filter_index = left_index * self.right_filter.num_bins + right_index
return filter_index
def get_pandas_dataframe(self, datasize, summary=None, **kwargs):
def get_pandas_dataframe(self, datasize, summary=None):
"""Builds a Pandas DataFrame for the CrossFilter's bins.
This method constructs a Pandas DataFrame object for the CrossFilter
@ -454,14 +454,6 @@ class CrossFilter(object):
column with a geometric "path" to each distribcell instance.
NOTE: This option requires the OpenCG Python package.
Keyword arguments
-----------------
energy_fmt : None or string
If a format string is provided, energy and energyout filter bins
will be converted from floats to strings using the given format. If
None is provided, the values will be left as floats. The default is
'{:.2e}'.
Returns
-------
pandas.DataFrame
@ -480,15 +472,12 @@ class CrossFilter(object):
# If left and right filters are identical, do not combine bins
if self.left_filter == self.right_filter:
df = self.left_filter.get_pandas_dataframe(datasize, summary,
**kwargs)
df = self.left_filter.get_pandas_dataframe(datasize, summary)
# If left and right filters are different, combine their bins
else:
left_df = self.left_filter.get_pandas_dataframe(datasize, summary,
**kwargs)
right_df = self.right_filter.get_pandas_dataframe(datasize, summary,
**kwargs)
left_df = self.left_filter.get_pandas_dataframe(datasize, summary)
right_df = self.right_filter.get_pandas_dataframe(datasize, summary)
left_df = left_df.astype(str)
right_df = right_df.astype(str)
df = '(' + left_df + ' ' + self.binary_op + ' ' + right_df + ')'
@ -842,7 +831,7 @@ class AggregateFilter(object):
else:
return 0
def get_pandas_dataframe(self, datasize, summary=None, **kwargs):
def get_pandas_dataframe(self, datasize, summary=None):
"""Builds a Pandas DataFrame for the AggregateFilter's bins.
This method constructs a Pandas DataFrame object for the AggregateFilter

View file

@ -462,7 +462,7 @@ class Filter(object):
return filter_bin
def get_pandas_dataframe(self, data_size, summary=None, **kwargs):
def get_pandas_dataframe(self, data_size, summary=None):
"""Builds a Pandas DataFrame for the Filter's bins.
This method constructs a Pandas DataFrame object for the filter with
@ -484,14 +484,6 @@ class Filter(object):
column with a geometric "path" to each distribcell instance.
NOTE: This option requires the OpenCG Python package.
Keyword arguments
-----------------
energy_fmt : None or string
If a format string is provided, energy and energyout filter bins
will be converted from floats to strings using the given format. If
None is provided, the values will be left as floats. The default is
'{:.2e}'.
Returns
-------
pandas.DataFrame
@ -735,12 +727,6 @@ class Filter(object):
lo_bins = np.tile(lo_bins, tile_factor)
hi_bins = np.tile(hi_bins, tile_factor)
# Format the energy values, if necessary.
energy_fmt = kwargs.setdefault('energy_fmt', '{:.2e}')
if energy_fmt is not None:
lo_bins = [energy_fmt.format(E) for E in lo_bins]
hi_bins = [energy_fmt.format(E) for E in hi_bins]
# Add the new energy columns to the DataFrame.
df.loc[:, self.type + ' low [MeV]'] = lo_bins
df.loc[:, self.type + ' high [MeV]'] = hi_bins

View file

@ -2,6 +2,7 @@ from __future__ import division
from collections import Iterable, defaultdict
import copy
from functools import partial
import os
import pickle
import itertools
@ -1244,7 +1245,7 @@ class Tally(object):
return data
def get_pandas_dataframe(self, filters=True, nuclides=True,
scores=True, summary=None, **kwargs):
scores=True, summary=None, float_format='{:.2e}'):
"""Build a Pandas DataFrame for the Tally data.
This method constructs a Pandas DataFrame object for the Tally data
@ -1268,14 +1269,9 @@ class Tally(object):
information in the Summary object is embedded into a Multi-index
column with a geometric "path" to each distribcell intance.
NOTE: This option requires the OpenCG Python package.
Keyword arguments
-----------------
energy_fmt : None or string
If a format string is provided, energy and energyout filter bins
will be converted from floats to strings using the given format. If
None is provided, the values will be left as floats. The default is
'{:.2e}'.
float_format : string
All floats in the DataFrame will be formatted using the given
format string before printing.
Returns
-------
@ -1324,8 +1320,7 @@ class Tally(object):
# Append each Filter's DataFrame to the overall DataFrame
for self_filter in self.filters:
filter_df = self_filter.get_pandas_dataframe(data_size, summary,
**kwargs)
filter_df = self_filter.get_pandas_dataframe(data_size, summary)
df = pd.concat([df, filter_df], axis=1)
# Include DataFrame column for nuclides if user requested it
@ -1376,6 +1371,10 @@ class Tally(object):
# Create and set a MultiIndex for the DataFrame's columns
df.columns = pd.MultiIndex.from_tuples(columns)
# Modify the df.to_string method so that it prints formatted strings.
# Credit to http://stackoverflow.com/users/3657742/chrisb for this trick
df.to_string = partial(df.to_string, float_format=float_format.format)
return df
def get_reshaped_data(self, value='mean'):