diff --git a/docs/source/pythonapi/examples/mgxs-part-ii.ipynb b/docs/source/pythonapi/examples/mgxs-part-ii.ipynb index b41b32241..764bb0d41 100644 --- a/docs/source/pythonapi/examples/mgxs-part-ii.ipynb +++ b/docs/source/pythonapi/examples/mgxs-part-ii.ipynb @@ -34,14 +34,14 @@ "name": "stderr", "output_type": "stream", "text": [ - "/usr/local/lib/python2.7/dist-packages/matplotlib/__init__.py:1318: UserWarning: This call to matplotlib.use() has no effect\n", + "/usr/lib/pymodules/python2.7/matplotlib/__init__.py:1173: 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", "\n", " warnings.warn(_use_error_msg)\n", - "/usr/local/lib/python2.7/dist-packages/ipykernel/__main__.py:11: QAWarning: pyne.rxname is not yet QA compliant.\n", - "/usr/local/lib/python2.7/dist-packages/ipykernel/__main__.py:11: QAWarning: pyne.ace is not yet QA compliant.\n" + "/usr/local/lib/python2.7/dist-packages/IPython/kernel/__main__.py:11: QAWarning: pyne.rxname is not yet QA compliant.\n", + "/usr/local/lib/python2.7/dist-packages/IPython/kernel/__main__.py:11: QAWarning: pyne.ace is not yet QA compliant.\n" ] } ], @@ -388,7 +388,7 @@ }, { "cell_type": "code", - "execution_count": 13, + "execution_count": null, "metadata": { "collapsed": false }, @@ -425,7 +425,7 @@ }, { "cell_type": "code", - "execution_count": 14, + "execution_count": null, "metadata": { "collapsed": false }, @@ -450,8 +450,9 @@ " Copyright: 2011-2015 Massachusetts Institute of Technology\n", " License: http://mit-crpg.github.io/openmc/license.html\n", " Version: 0.7.1\n", - " Git SHA1: ea9fb637f63f9374c7436456141afa850b84acf9\n", - " Date/Time: 2016-01-14 08:03:40\n", + " Git SHA1: 263266f4f8807fd38c6ac282fae259ae73fa1eee\n", + " Date/Time: 2016-01-20 17:57:53\n", + " MPI Processes: 1\n", "\n", " ===========================================================================\n", " ========================> INITIALIZATION <=========================\n", @@ -509,101 +510,8 @@ " 30/1 1.23670 1.22056 +/- 0.00368\n", " 31/1 1.21396 1.22024 +/- 0.00352\n", " 32/1 1.21389 1.21995 +/- 0.00337\n", - " 33/1 1.24649 1.22111 +/- 0.00342\n", - " 34/1 1.23204 1.22156 +/- 0.00330\n", - " 35/1 1.20768 1.22101 +/- 0.00322\n", - " 36/1 1.22271 1.22107 +/- 0.00309\n", - " 37/1 1.21796 1.22096 +/- 0.00298\n", - " 38/1 1.23842 1.22158 +/- 0.00293\n", - " 39/1 1.23080 1.22190 +/- 0.00285\n", - " 40/1 1.23572 1.22236 +/- 0.00279\n", - " 41/1 1.21691 1.22218 +/- 0.00271\n", - " 42/1 1.24616 1.22293 +/- 0.00272\n", - " 43/1 1.21903 1.22282 +/- 0.00264\n", - " 44/1 1.22967 1.22302 +/- 0.00257\n", - " 45/1 1.22053 1.22295 +/- 0.00250\n", - " 46/1 1.24087 1.22344 +/- 0.00248\n", - " 47/1 1.20251 1.22288 +/- 0.00248\n", - " 48/1 1.20331 1.22236 +/- 0.00246\n", - " 49/1 1.22724 1.22249 +/- 0.00240\n", - " 50/1 1.24798 1.22313 +/- 0.00243\n", - " Triggers unsatisfied, max unc./thresh. is 1.32110 for scatter-p1 in tally 10054\n", - " The estimated number of batches is 80\n", - " Creating state point statepoint.050.h5...\n", - " 51/1 1.22253 1.22311 +/- 0.00237\n", - " 52/1 1.24330 1.22359 +/- 0.00236\n", - " 53/1 1.23251 1.22380 +/- 0.00231\n", - " 54/1 1.21133 1.22352 +/- 0.00228\n", - " 55/1 1.24503 1.22399 +/- 0.00228\n", - " 56/1 1.22013 1.22391 +/- 0.00223\n", - " 57/1 1.23877 1.22423 +/- 0.00220\n", - " 58/1 1.23793 1.22451 +/- 0.00218\n", - " 59/1 1.21018 1.22422 +/- 0.00215\n", - " 60/1 1.22417 1.22422 +/- 0.00211\n", - " 61/1 1.23094 1.22435 +/- 0.00207\n", - " 62/1 1.23310 1.22452 +/- 0.00204\n", - " 63/1 1.22488 1.22453 +/- 0.00200\n", - " 64/1 1.22702 1.22457 +/- 0.00196\n", - " 65/1 1.18834 1.22391 +/- 0.00204\n", - " 66/1 1.23112 1.22404 +/- 0.00200\n", - " 67/1 1.21611 1.22390 +/- 0.00197\n", - " 68/1 1.22513 1.22392 +/- 0.00194\n", - " 69/1 1.21741 1.22381 +/- 0.00191\n", - " 70/1 1.22484 1.22383 +/- 0.00188\n", - " 71/1 1.19662 1.22338 +/- 0.00190\n", - " 72/1 1.23315 1.22354 +/- 0.00187\n", - " 73/1 1.22796 1.22361 +/- 0.00185\n", - " 74/1 1.21417 1.22346 +/- 0.00182\n", - " 75/1 1.21020 1.22326 +/- 0.00181\n", - " 76/1 1.23413 1.22343 +/- 0.00179\n", - " 77/1 1.22184 1.22340 +/- 0.00176\n", - " 78/1 1.20309 1.22310 +/- 0.00176\n", - " 79/1 1.23458 1.22327 +/- 0.00174\n", - " 80/1 1.20724 1.22304 +/- 0.00173\n", - " Triggers satisfied for batch 80\n", - " Creating state point statepoint.080.h5...\n", - "\n", - " ===========================================================================\n", - " ======================> SIMULATION FINISHED <======================\n", - " ===========================================================================\n", - "\n", - "\n", - " =======================> TIMING STATISTICS <=======================\n", - "\n", - " Total time for initialization = 3.8300E-01 seconds\n", - " Reading cross sections = 1.0700E-01 seconds\n", - " Total time in simulation = 2.5629E+02 seconds\n", - " Time in transport only = 2.5623E+02 seconds\n", - " Time in inactive batches = 1.5743E+01 seconds\n", - " Time in active batches = 2.4055E+02 seconds\n", - " Time synchronizing fission bank = 3.4000E-02 seconds\n", - " Sampling source sites = 2.1000E-02 seconds\n", - " SEND/RECV source sites = 1.3000E-02 seconds\n", - " Time accumulating tallies = 1.0000E-03 seconds\n", - " Total time for finalization = 7.0000E-03 seconds\n", - " Total time elapsed = 2.5674E+02 seconds\n", - " Calculation Rate (inactive) = 6352.03 neutrons/second\n", - " Calculation Rate (active) = 1662.87 neutrons/second\n", - "\n", - " ============================> RESULTS <============================\n", - "\n", - " k-effective (Collision) = 1.22327 +/- 0.00148\n", - " k-effective (Track-length) = 1.22304 +/- 0.00173\n", - " k-effective (Absorption) = 1.22407 +/- 0.00129\n", - " Combined k-effective = 1.22373 +/- 0.00113\n", - " Leakage Fraction = 0.00000 +/- 0.00000\n", - "\n" + " 33/1 1.24649 1.22111 +/- 0.00342\n" ] - }, - { - "data": { - "text/plain": [ - "0" - ] - }, - "execution_count": 14, - "metadata": {}, - "output_type": "execute_result" } ], "source": [ @@ -628,7 +536,7 @@ }, { "cell_type": "code", - "execution_count": 15, + "execution_count": null, "metadata": { "collapsed": false }, @@ -647,7 +555,7 @@ }, { "cell_type": "code", - "execution_count": 16, + "execution_count": null, "metadata": { "collapsed": true }, @@ -667,7 +575,7 @@ }, { "cell_type": "code", - "execution_count": 17, + "execution_count": null, "metadata": { "collapsed": false }, @@ -702,46 +610,11 @@ }, { "cell_type": "code", - "execution_count": 18, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "Multi-Group XS\n", - "\tReaction Type =\tnu-fission\n", - "\tDomain Type =\tcell\n", - "\tDomain ID =\t10000\n", - "\tNuclide =\tU-235\n", - "\tCross Sections [barns]:\n", - " Group 1 [0.821 - 20.0 MeV]:\t3.31e+00 +/- 1.88e-01%\n", - " Group 2 [0.00553 - 0.821 MeV]:\t3.97e+00 +/- 1.24e-01%\n", - " Group 3 [4e-06 - 0.00553 MeV]:\t5.50e+01 +/- 2.02e-01%\n", - " Group 4 [6.25e-07 - 4e-06 MeV]:\t8.83e+01 +/- 3.56e-01%\n", - " Group 5 [2.8e-07 - 6.25e-07 MeV]:\t2.90e+02 +/- 4.54e-01%\n", - " Group 6 [1.4e-07 - 2.8e-07 MeV]:\t4.49e+02 +/- 4.10e-01%\n", - " Group 7 [5.8e-08 - 1.4e-07 MeV]:\t6.87e+02 +/- 2.56e-01%\n", - " Group 8 [0.0 - 5.8e-08 MeV]:\t1.44e+03 +/- 2.82e-01%\n", - "\n", - "\tNuclide =\tU-238\n", - "\tCross Sections [barns]:\n", - " Group 1 [0.821 - 20.0 MeV]:\t1.06e+00 +/- 2.30e-01%\n", - " Group 2 [0.00553 - 0.821 MeV]:\t1.21e-03 +/- 2.25e-01%\n", - " Group 3 [4e-06 - 0.00553 MeV]:\t5.82e-04 +/- 3.09e+00%\n", - " Group 4 [6.25e-07 - 4e-06 MeV]:\t6.54e-06 +/- 3.27e-01%\n", - " Group 5 [2.8e-07 - 6.25e-07 MeV]:\t1.07e-05 +/- 4.39e-01%\n", - " Group 6 [1.4e-07 - 2.8e-07 MeV]:\t1.55e-05 +/- 4.12e-01%\n", - " Group 7 [5.8e-08 - 1.4e-07 MeV]:\t2.30e-05 +/- 2.57e-01%\n", - " Group 8 [0.0 - 5.8e-08 MeV]:\t4.24e-05 +/- 2.81e-01%\n", - "\n", - "\n", - "\n" - ] - } - ], + "outputs": [], "source": [ "nufission = xs_library[fuel_cell.id]['nu-fission']\n", "nufission.print_xs(xs_type='micro', nuclides=['U-235', 'U-238'])" @@ -756,34 +629,11 @@ }, { "cell_type": "code", - "execution_count": 19, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "Multi-Group XS\n", - "\tReaction Type =\tnu-fission\n", - "\tDomain Type =\tcell\n", - "\tDomain ID =\t10000\n", - "\tCross Sections [cm^-1]:\n", - " Group 1 [0.821 - 20.0 MeV]:\t2.52e-02 +/- 2.19e-01%\n", - " Group 2 [0.00553 - 0.821 MeV]:\t1.51e-03 +/- 1.22e-01%\n", - " Group 3 [4e-06 - 0.00553 MeV]:\t2.06e-02 +/- 2.02e-01%\n", - " Group 4 [6.25e-07 - 4e-06 MeV]:\t3.31e-02 +/- 3.56e-01%\n", - " Group 5 [2.8e-07 - 6.25e-07 MeV]:\t1.09e-01 +/- 4.54e-01%\n", - " Group 6 [1.4e-07 - 2.8e-07 MeV]:\t1.69e-01 +/- 4.10e-01%\n", - " Group 7 [5.8e-08 - 1.4e-07 MeV]:\t2.58e-01 +/- 2.56e-01%\n", - " Group 8 [0.0 - 5.8e-08 MeV]:\t5.40e-01 +/- 2.82e-01%\n", - "\n", - "\n", - "\n" - ] - } - ], + "outputs": [], "source": [ "nufission = xs_library[fuel_cell.id]['nu-fission']\n", "nufission.print_xs(xs_type='macro', nuclides='sum')" @@ -798,151 +648,11 @@ }, { "cell_type": "code", - "execution_count": 20, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stderr", - "output_type": "stream", - "text": [ - "/usr/local/lib/python2.7/dist-packages/numpy/lib/shape_base.py:872: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n", - " return c.reshape(shape_out)\n", - "/usr/local/lib/python2.7/dist-packages/numpy/lib/shape_base.py:872: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n", - " return c.reshape(shape_out)\n" - ] - }, - { - "data": { - "text/html": [ - "
\n", - "\n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - "
cellgroup ingroup outnuclidemeanstd. dev.
1261000211H-10.2340220.003645
1271000211O-161.5603050.006280
1241000212H-11.5880250.002815
1251000212O-160.2851470.001392
1221000213H-10.0107760.000186
1231000213O-160.0000000.000000
1201000214H-10.0000230.000010
1211000214O-160.0000000.000000
1181000215H-10.0000000.000000
1191000215O-160.0000000.000000
\n", - "
" - ], - "text/plain": [ - " cell group in group out nuclide mean std. dev.\n", - "126 10002 1 1 H-1 0.234022 0.003645\n", - "127 10002 1 1 O-16 1.560305 0.006280\n", - "124 10002 1 2 H-1 1.588025 0.002815\n", - "125 10002 1 2 O-16 0.285147 0.001392\n", - "122 10002 1 3 H-1 0.010776 0.000186\n", - "123 10002 1 3 O-16 0.000000 0.000000\n", - "120 10002 1 4 H-1 0.000023 0.000010\n", - "121 10002 1 4 O-16 0.000000 0.000000\n", - "118 10002 1 5 H-1 0.000000 0.000000\n", - "119 10002 1 5 O-16 0.000000 0.000000" - ] - }, - "execution_count": 20, - "metadata": {}, - "output_type": "execute_result" - } - ], + "outputs": [], "source": [ "nuscatter = xs_library[moderator_cell.id]['nu-scatter']\n", "df = nuscatter.get_pandas_dataframe(xs_type='micro')\n", @@ -958,7 +668,7 @@ }, { "cell_type": "code", - "execution_count": 21, + "execution_count": null, "metadata": { "collapsed": true }, @@ -980,143 +690,22 @@ }, { "cell_type": "code", - "execution_count": 22, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "Multi-Group XS\n", - "\tReaction Type =\ttransport\n", - "\tDomain Type =\tcell\n", - "\tDomain ID =\t10000\n", - "\tNuclide =\tU-235\n", - "\tCross Sections [cm^-1]:\n", - " Group 1 [6.25e-07 - 20.0 MeV]:\t7.81e-03 +/- 4.75e-01%\n", - " Group 2 [0.0 - 6.25e-07 MeV]:\t1.82e-01 +/- 1.89e-01%\n", - "\n", - "\tNuclide =\tU-238\n", - "\tCross Sections [cm^-1]:\n", - " Group 1 [6.25e-07 - 20.0 MeV]:\t2.17e-01 +/- 1.31e-01%\n", - " Group 2 [0.0 - 6.25e-07 MeV]:\t2.53e-01 +/- 2.08e-01%\n", - "\n", - "\tNuclide =\tO-16\n", - "\tCross Sections [cm^-1]:\n", - " Group 1 [6.25e-07 - 20.0 MeV]:\t1.45e-01 +/- 1.50e-01%\n", - " Group 2 [0.0 - 6.25e-07 MeV]:\t1.74e-01 +/- 2.66e-01%\n", - "\n", - "\n", - "\n" - ] - } - ], + "outputs": [], "source": [ "condensed_xs.print_xs()" ] }, { "cell_type": "code", - "execution_count": 23, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stderr", - "output_type": "stream", - "text": [ - "/usr/local/lib/python2.7/dist-packages/numpy/lib/shape_base.py:872: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n", - " return c.reshape(shape_out)\n", - "/usr/local/lib/python2.7/dist-packages/numpy/lib/shape_base.py:872: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n", - " return c.reshape(shape_out)\n" - ] - }, - { - "data": { - "text/html": [ - "
\n", - "\n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - "
cellgroup innuclidemeanstd. dev.
3100001U-23520.8281270.098842
4100001U-2389.5822950.012550
5100001O-163.1573580.004725
0100002U-235485.2176490.916465
1100002U-23811.1760810.023196
2100002O-163.7881670.010090
\n", - "
" - ], - "text/plain": [ - " cell group in nuclide mean std. dev.\n", - "3 10000 1 U-235 20.828127 0.098842\n", - "4 10000 1 U-238 9.582295 0.012550\n", - "5 10000 1 O-16 3.157358 0.004725\n", - "0 10000 2 U-235 485.217649 0.916465\n", - "1 10000 2 U-238 11.176081 0.023196\n", - "2 10000 2 O-16 3.788167 0.010090" - ] - }, - "execution_count": 23, - "metadata": {}, - "output_type": "execute_result" - } - ], + "outputs": [], "source": [ "df = condensed_xs.get_pandas_dataframe(xs_type='micro')\n", "df" @@ -1138,7 +727,7 @@ }, { "cell_type": "code", - "execution_count": 24, + "execution_count": null, "metadata": { "collapsed": false }, @@ -1157,7 +746,7 @@ }, { "cell_type": "code", - "execution_count": 25, + "execution_count": null, "metadata": { "collapsed": false }, @@ -1202,187 +791,13 @@ }, { "cell_type": "code", - "execution_count": 26, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "[ NORMAL ] Ray tracing for track segmentation...\n", - "[ NORMAL ] Dumping tracks to file...\n", - "[ NORMAL ] Computing the eigenvalue...\n", - "[ NORMAL ] Iteration 0:\tk_eff = 0.574633\tres = 0.000E+00\n", - "[ NORMAL ] Iteration 1:\tk_eff = 0.679931\tres = 4.254E-01\n", - "[ NORMAL ] Iteration 2:\tk_eff = 0.660910\tres = 1.832E-01\n", - "[ NORMAL ] Iteration 3:\tk_eff = 0.658975\tres = 2.798E-02\n", - "[ NORMAL ] Iteration 4:\tk_eff = 0.642976\tres = 2.927E-03\n", - "[ NORMAL ] Iteration 5:\tk_eff = 0.625710\tres = 2.428E-02\n", - "[ NORMAL ] Iteration 6:\tk_eff = 0.606520\tres = 2.685E-02\n", - "[ NORMAL ] Iteration 7:\tk_eff = 0.587277\tres = 3.067E-02\n", - "[ NORMAL ] Iteration 8:\tk_eff = 0.568777\tres = 3.173E-02\n", - "[ NORMAL ] Iteration 9:\tk_eff = 0.551415\tres = 3.150E-02\n", - "[ NORMAL ] Iteration 10:\tk_eff = 0.535708\tres = 3.052E-02\n", - "[ NORMAL ] Iteration 11:\tk_eff = 0.521916\tres = 2.849E-02\n", - "[ NORMAL ] Iteration 12:\tk_eff = 0.510221\tres = 2.575E-02\n", - "[ NORMAL ] Iteration 13:\tk_eff = 0.500691\tres = 2.241E-02\n", - "[ NORMAL ] Iteration 14:\tk_eff = 0.493392\tres = 1.868E-02\n", - "[ NORMAL ] Iteration 15:\tk_eff = 0.488318\tres = 1.458E-02\n", - "[ NORMAL ] Iteration 16:\tk_eff = 0.485438\tres = 1.028E-02\n", - "[ NORMAL ] Iteration 17:\tk_eff = 0.484705\tres = 5.896E-03\n", - "[ NORMAL ] Iteration 18:\tk_eff = 0.486046\tres = 1.510E-03\n", - "[ NORMAL ] Iteration 19:\tk_eff = 0.489362\tres = 2.766E-03\n", - "[ NORMAL ] Iteration 20:\tk_eff = 0.494546\tres = 6.824E-03\n", - "[ NORMAL ] Iteration 21:\tk_eff = 0.501481\tres = 1.059E-02\n", - "[ NORMAL ] Iteration 22:\tk_eff = 0.510042\tres = 1.402E-02\n", - "[ NORMAL ] Iteration 23:\tk_eff = 0.520095\tres = 1.707E-02\n", - "[ NORMAL ] Iteration 24:\tk_eff = 0.531508\tres = 1.971E-02\n", - "[ NORMAL ] Iteration 25:\tk_eff = 0.544144\tres = 2.194E-02\n", - "[ NORMAL ] Iteration 26:\tk_eff = 0.557872\tres = 2.378E-02\n", - "[ NORMAL ] Iteration 27:\tk_eff = 0.572557\tres = 2.523E-02\n", - "[ NORMAL ] Iteration 28:\tk_eff = 0.588073\tres = 2.632E-02\n", - "[ NORMAL ] Iteration 29:\tk_eff = 0.604293\tres = 2.710E-02\n", - "[ NORMAL ] Iteration 30:\tk_eff = 0.621101\tres = 2.758E-02\n", - "[ NORMAL ] Iteration 31:\tk_eff = 0.638382\tres = 2.781E-02\n", - "[ NORMAL ] Iteration 32:\tk_eff = 0.656032\tres = 2.782E-02\n", - "[ NORMAL ] Iteration 33:\tk_eff = 0.673950\tres = 2.765E-02\n", - "[ NORMAL ] Iteration 34:\tk_eff = 0.692043\tres = 2.731E-02\n", - "[ NORMAL ] Iteration 35:\tk_eff = 0.710227\tres = 2.685E-02\n", - "[ NORMAL ] Iteration 36:\tk_eff = 0.728423\tres = 2.628E-02\n", - "[ NORMAL ] Iteration 37:\tk_eff = 0.746558\tres = 2.562E-02\n", - "[ NORMAL ] Iteration 38:\tk_eff = 0.764569\tres = 2.490E-02\n", - "[ NORMAL ] Iteration 39:\tk_eff = 0.782396\tres = 2.412E-02\n", - "[ NORMAL ] Iteration 40:\tk_eff = 0.799990\tres = 2.332E-02\n", - "[ NORMAL ] Iteration 41:\tk_eff = 0.817301\tres = 2.249E-02\n", - "[ NORMAL ] Iteration 42:\tk_eff = 0.834292\tres = 2.164E-02\n", - "[ NORMAL ] Iteration 43:\tk_eff = 0.850928\tres = 2.079E-02\n", - "[ NORMAL ] Iteration 44:\tk_eff = 0.867178\tres = 1.994E-02\n", - "[ NORMAL ] Iteration 45:\tk_eff = 0.883018\tres = 1.910E-02\n", - "[ NORMAL ] Iteration 46:\tk_eff = 0.898427\tres = 1.827E-02\n", - "[ NORMAL ] Iteration 47:\tk_eff = 0.913390\tres = 1.745E-02\n", - "[ NORMAL ] Iteration 48:\tk_eff = 0.927892\tres = 1.665E-02\n", - "[ NORMAL ] Iteration 49:\tk_eff = 0.941926\tres = 1.588E-02\n", - "[ NORMAL ] Iteration 50:\tk_eff = 0.955484\tres = 1.512E-02\n", - "[ NORMAL ] Iteration 51:\tk_eff = 0.968563\tres = 1.439E-02\n", - "[ NORMAL ] Iteration 52:\tk_eff = 0.981162\tres = 1.369E-02\n", - "[ NORMAL ] Iteration 53:\tk_eff = 0.993283\tres = 1.301E-02\n", - "[ NORMAL ] Iteration 54:\tk_eff = 1.004929\tres = 1.235E-02\n", - "[ NORMAL ] Iteration 55:\tk_eff = 1.016105\tres = 1.172E-02\n", - "[ NORMAL ] Iteration 56:\tk_eff = 1.026818\tres = 1.112E-02\n", - "[ NORMAL ] Iteration 57:\tk_eff = 1.037075\tres = 1.054E-02\n", - "[ NORMAL ] Iteration 58:\tk_eff = 1.046885\tres = 9.989E-03\n", - "[ NORMAL ] Iteration 59:\tk_eff = 1.056259\tres = 9.459E-03\n", - "[ NORMAL ] Iteration 60:\tk_eff = 1.065206\tres = 8.954E-03\n", - "[ NORMAL ] Iteration 61:\tk_eff = 1.073740\tres = 8.471E-03\n", - "[ NORMAL ] Iteration 62:\tk_eff = 1.081872\tres = 8.012E-03\n", - "[ NORMAL ] Iteration 63:\tk_eff = 1.089614\tres = 7.573E-03\n", - "[ NORMAL ] Iteration 64:\tk_eff = 1.096980\tres = 7.156E-03\n", - "[ NORMAL ] Iteration 65:\tk_eff = 1.103981\tres = 6.760E-03\n", - "[ NORMAL ] Iteration 66:\tk_eff = 1.110631\tres = 6.382E-03\n", - "[ NORMAL ] Iteration 67:\tk_eff = 1.116944\tres = 6.024E-03\n", - "[ NORMAL ] Iteration 68:\tk_eff = 1.122932\tres = 5.684E-03\n", - "[ NORMAL ] Iteration 69:\tk_eff = 1.128608\tres = 5.361E-03\n", - "[ NORMAL ] Iteration 70:\tk_eff = 1.133985\tres = 5.055E-03\n", - "[ NORMAL ] Iteration 71:\tk_eff = 1.139076\tres = 4.764E-03\n", - "[ NORMAL ] Iteration 72:\tk_eff = 1.143892\tres = 4.490E-03\n", - "[ NORMAL ] Iteration 73:\tk_eff = 1.148447\tres = 4.228E-03\n", - "[ NORMAL ] Iteration 74:\tk_eff = 1.152753\tres = 3.982E-03\n", - "[ NORMAL ] Iteration 75:\tk_eff = 1.156819\tres = 3.749E-03\n", - "[ NORMAL ] Iteration 76:\tk_eff = 1.160659\tres = 3.528E-03\n", - "[ NORMAL ] Iteration 77:\tk_eff = 1.164283\tres = 3.319E-03\n", - "[ NORMAL ] Iteration 78:\tk_eff = 1.167701\tres = 3.122E-03\n", - "[ NORMAL ] Iteration 79:\tk_eff = 1.170923\tres = 2.936E-03\n", - "[ NORMAL ] Iteration 80:\tk_eff = 1.173960\tres = 2.759E-03\n", - "[ NORMAL ] Iteration 81:\tk_eff = 1.176822\tres = 2.594E-03\n", - "[ NORMAL ] Iteration 82:\tk_eff = 1.179516\tres = 2.437E-03\n", - "[ NORMAL ] Iteration 83:\tk_eff = 1.182052\tres = 2.289E-03\n", - "[ NORMAL ] Iteration 84:\tk_eff = 1.184438\tres = 2.150E-03\n", - "[ NORMAL ] Iteration 85:\tk_eff = 1.186682\tres = 2.018E-03\n", - "[ NORMAL ] Iteration 86:\tk_eff = 1.188792\tres = 1.895E-03\n", - "[ NORMAL ] Iteration 87:\tk_eff = 1.190775\tres = 1.778E-03\n", - "[ NORMAL ] Iteration 88:\tk_eff = 1.192639\tres = 1.668E-03\n", - "[ NORMAL ] Iteration 89:\tk_eff = 1.194389\tres = 1.565E-03\n", - "[ NORMAL ] Iteration 90:\tk_eff = 1.196032\tres = 1.468E-03\n", - "[ NORMAL ] Iteration 91:\tk_eff = 1.197575\tres = 1.375E-03\n", - "[ NORMAL ] Iteration 92:\tk_eff = 1.199023\tres = 1.290E-03\n", - "[ NORMAL ] Iteration 93:\tk_eff = 1.200381\tres = 1.209E-03\n", - "[ NORMAL ] Iteration 94:\tk_eff = 1.201654\tres = 1.133E-03\n", - "[ NORMAL ] Iteration 95:\tk_eff = 1.202848\tres = 1.061E-03\n", - "[ NORMAL ] Iteration 96:\tk_eff = 1.203968\tres = 9.941E-04\n", - "[ NORMAL ] Iteration 97:\tk_eff = 1.205016\tres = 9.308E-04\n", - "[ NORMAL ] Iteration 98:\tk_eff = 1.205999\tres = 8.707E-04\n", - "[ NORMAL ] Iteration 99:\tk_eff = 1.206921\tres = 8.159E-04\n", - "[ NORMAL ] Iteration 100:\tk_eff = 1.207783\tres = 7.635E-04\n", - "[ NORMAL ] Iteration 101:\tk_eff = 1.208590\tres = 7.145E-04\n", - "[ NORMAL ] Iteration 102:\tk_eff = 1.209345\tres = 6.684E-04\n", - "[ NORMAL ] Iteration 103:\tk_eff = 1.210053\tres = 6.250E-04\n", - "[ NORMAL ] Iteration 104:\tk_eff = 1.210714\tres = 5.848E-04\n", - "[ NORMAL ] Iteration 105:\tk_eff = 1.211333\tres = 5.466E-04\n", - "[ NORMAL ] Iteration 106:\tk_eff = 1.211913\tres = 5.115E-04\n", - "[ NORMAL ] Iteration 107:\tk_eff = 1.212453\tres = 4.784E-04\n", - "[ NORMAL ] Iteration 108:\tk_eff = 1.212960\tres = 4.462E-04\n", - "[ NORMAL ] Iteration 109:\tk_eff = 1.213433\tres = 4.175E-04\n", - "[ NORMAL ] Iteration 110:\tk_eff = 1.213875\tres = 3.901E-04\n", - "[ NORMAL ] Iteration 111:\tk_eff = 1.214288\tres = 3.643E-04\n", - "[ NORMAL ] Iteration 112:\tk_eff = 1.214674\tres = 3.404E-04\n", - "[ NORMAL ] Iteration 113:\tk_eff = 1.215036\tres = 3.181E-04\n", - "[ NORMAL ] Iteration 114:\tk_eff = 1.215372\tres = 2.972E-04\n", - "[ NORMAL ] Iteration 115:\tk_eff = 1.215686\tres = 2.772E-04\n", - "[ NORMAL ] Iteration 116:\tk_eff = 1.215980\tres = 2.586E-04\n", - "[ NORMAL ] Iteration 117:\tk_eff = 1.216255\tres = 2.416E-04\n", - "[ NORMAL ] Iteration 118:\tk_eff = 1.216511\tres = 2.258E-04\n", - "[ NORMAL ] Iteration 119:\tk_eff = 1.216749\tres = 2.105E-04\n", - "[ NORMAL ] Iteration 120:\tk_eff = 1.216973\tres = 1.963E-04\n", - "[ NORMAL ] Iteration 121:\tk_eff = 1.217181\tres = 1.834E-04\n", - "[ NORMAL ] Iteration 122:\tk_eff = 1.217375\tres = 1.710E-04\n", - "[ NORMAL ] Iteration 123:\tk_eff = 1.217556\tres = 1.595E-04\n", - "[ NORMAL ] Iteration 124:\tk_eff = 1.217725\tres = 1.488E-04\n", - "[ NORMAL ] Iteration 125:\tk_eff = 1.217883\tres = 1.389E-04\n", - "[ NORMAL ] Iteration 126:\tk_eff = 1.218030\tres = 1.296E-04\n", - "[ NORMAL ] Iteration 127:\tk_eff = 1.218167\tres = 1.208E-04\n", - "[ NORMAL ] Iteration 128:\tk_eff = 1.218294\tres = 1.128E-04\n", - "[ NORMAL ] Iteration 129:\tk_eff = 1.218414\tres = 1.041E-04\n", - "[ NORMAL ] Iteration 130:\tk_eff = 1.218525\tres = 9.793E-05\n", - "[ NORMAL ] Iteration 131:\tk_eff = 1.218628\tres = 9.128E-05\n", - "[ NORMAL ] Iteration 132:\tk_eff = 1.218725\tres = 8.502E-05\n", - "[ NORMAL ] Iteration 133:\tk_eff = 1.218815\tres = 7.903E-05\n", - "[ NORMAL ] Iteration 134:\tk_eff = 1.218899\tres = 7.376E-05\n", - "[ NORMAL ] Iteration 135:\tk_eff = 1.218976\tres = 6.888E-05\n", - "[ NORMAL ] Iteration 136:\tk_eff = 1.219049\tres = 6.362E-05\n", - "[ NORMAL ] Iteration 137:\tk_eff = 1.219117\tres = 5.979E-05\n", - "[ NORMAL ] Iteration 138:\tk_eff = 1.219180\tres = 5.591E-05\n", - "[ NORMAL ] Iteration 139:\tk_eff = 1.219239\tres = 5.181E-05\n", - "[ NORMAL ] Iteration 140:\tk_eff = 1.219294\tres = 4.824E-05\n", - "[ NORMAL ] Iteration 141:\tk_eff = 1.219345\tres = 4.496E-05\n", - "[ NORMAL ] Iteration 142:\tk_eff = 1.219392\tres = 4.197E-05\n", - "[ NORMAL ] Iteration 143:\tk_eff = 1.219437\tres = 3.871E-05\n", - "[ NORMAL ] Iteration 144:\tk_eff = 1.219478\tres = 3.666E-05\n", - "[ NORMAL ] Iteration 145:\tk_eff = 1.219516\tres = 3.373E-05\n", - "[ NORMAL ] Iteration 146:\tk_eff = 1.219552\tres = 3.112E-05\n", - "[ NORMAL ] Iteration 147:\tk_eff = 1.219585\tres = 2.941E-05\n", - "[ NORMAL ] Iteration 148:\tk_eff = 1.219616\tres = 2.711E-05\n", - "[ NORMAL ] Iteration 149:\tk_eff = 1.219645\tres = 2.542E-05\n", - "[ NORMAL ] Iteration 150:\tk_eff = 1.219671\tres = 2.369E-05\n", - "[ NORMAL ] Iteration 151:\tk_eff = 1.219696\tres = 2.207E-05\n", - "[ NORMAL ] Iteration 152:\tk_eff = 1.219720\tres = 2.049E-05\n", - "[ NORMAL ] Iteration 153:\tk_eff = 1.219741\tres = 1.903E-05\n", - "[ NORMAL ] Iteration 154:\tk_eff = 1.219761\tres = 1.770E-05\n", - "[ NORMAL ] Iteration 155:\tk_eff = 1.219780\tres = 1.639E-05\n", - "[ NORMAL ] Iteration 156:\tk_eff = 1.219797\tres = 1.533E-05\n", - "[ NORMAL ] Iteration 157:\tk_eff = 1.219813\tres = 1.396E-05\n", - "[ NORMAL ] Iteration 158:\tk_eff = 1.219828\tres = 1.319E-05\n", - "[ NORMAL ] Iteration 159:\tk_eff = 1.219842\tres = 1.216E-05\n", - "[ NORMAL ] Iteration 160:\tk_eff = 1.219854\tres = 1.152E-05\n", - "[ NORMAL ] Iteration 161:\tk_eff = 1.219867\tres = 1.049E-05\n", - "[ NORMAL ] Iteration 162:\tk_eff = 1.219878\tres = 1.005E-05\n" - ] - } - ], + "outputs": [], "source": [ "# Generate tracks for OpenMOC\n", - "openmoc_geometry.initializeFlatSourceRegions()\n", "track_generator = openmoc.TrackGenerator(openmoc_geometry, num_azim=128, spacing=0.1)\n", "track_generator.generateTracks()\n", "\n", @@ -1400,21 +815,11 @@ }, { "cell_type": "code", - "execution_count": 27, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "openmc keff = 1.223729\n", - "openmoc keff = 1.219878\n", - "bias [pcm]: -385.1\n" - ] - } - ], + "outputs": [], "source": [ "# Print report of keff and bias with OpenMC\n", "openmoc_keff = solver.getKeff()\n", @@ -1435,7 +840,7 @@ }, { "cell_type": "code", - "execution_count": 28, + "execution_count": null, "metadata": { "collapsed": false }, @@ -1475,254 +880,13 @@ }, { "cell_type": "code", - "execution_count": 29, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "[ NORMAL ] Importing ray tracing data from file...\n", - "[ NORMAL ] Computing the eigenvalue...\n", - "[ NORMAL ] Iteration 0:\tk_eff = 0.495594\tres = 0.000E+00\n", - "[ NORMAL ] Iteration 1:\tk_eff = 0.557313\tres = 5.044E-01\n", - "[ NORMAL ] Iteration 2:\tk_eff = 0.518115\tres = 1.245E-01\n", - "[ NORMAL ] Iteration 3:\tk_eff = 0.509017\tres = 7.033E-02\n", - "[ NORMAL ] Iteration 4:\tk_eff = 0.496280\tres = 1.756E-02\n", - "[ NORMAL ] Iteration 5:\tk_eff = 0.488358\tres = 2.502E-02\n", - "[ NORMAL ] Iteration 6:\tk_eff = 0.482660\tres = 1.596E-02\n", - "[ NORMAL ] Iteration 7:\tk_eff = 0.479524\tres = 1.167E-02\n", - "[ NORMAL ] Iteration 8:\tk_eff = 0.478569\tres = 6.497E-03\n", - "[ NORMAL ] Iteration 9:\tk_eff = 0.479591\tres = 1.992E-03\n", - "[ NORMAL ] Iteration 10:\tk_eff = 0.482389\tres = 2.136E-03\n", - "[ NORMAL ] Iteration 11:\tk_eff = 0.486774\tres = 5.834E-03\n", - "[ NORMAL ] Iteration 12:\tk_eff = 0.492576\tres = 9.091E-03\n", - "[ NORMAL ] Iteration 13:\tk_eff = 0.499632\tres = 1.192E-02\n", - "[ NORMAL ] Iteration 14:\tk_eff = 0.507800\tres = 1.433E-02\n", - "[ NORMAL ] Iteration 15:\tk_eff = 0.516944\tres = 1.635E-02\n", - "[ NORMAL ] Iteration 16:\tk_eff = 0.526942\tres = 1.801E-02\n", - "[ NORMAL ] Iteration 17:\tk_eff = 0.537682\tres = 1.934E-02\n", - "[ NORMAL ] Iteration 18:\tk_eff = 0.549061\tres = 2.038E-02\n", - "[ NORMAL ] Iteration 19:\tk_eff = 0.560985\tres = 2.116E-02\n", - "[ NORMAL ] Iteration 20:\tk_eff = 0.573368\tres = 2.172E-02\n", - "[ NORMAL ] Iteration 21:\tk_eff = 0.586133\tres = 2.207E-02\n", - "[ NORMAL ] Iteration 22:\tk_eff = 0.599208\tres = 2.226E-02\n", - "[ NORMAL ] Iteration 23:\tk_eff = 0.612528\tres = 2.231E-02\n", - "[ NORMAL ] Iteration 24:\tk_eff = 0.626035\tres = 2.223E-02\n", - "[ NORMAL ] Iteration 25:\tk_eff = 0.639676\tres = 2.205E-02\n", - "[ NORMAL ] Iteration 26:\tk_eff = 0.653402\tres = 2.179E-02\n", - "[ NORMAL ] Iteration 27:\tk_eff = 0.667171\tres = 2.146E-02\n", - "[ NORMAL ] Iteration 28:\tk_eff = 0.680942\tres = 2.107E-02\n", - "[ NORMAL ] Iteration 29:\tk_eff = 0.694681\tres = 2.064E-02\n", - "[ NORMAL ] Iteration 30:\tk_eff = 0.708356\tres = 2.018E-02\n", - "[ NORMAL ] Iteration 31:\tk_eff = 0.721940\tres = 1.969E-02\n", - "[ NORMAL ] Iteration 32:\tk_eff = 0.735407\tres = 1.918E-02\n", - "[ NORMAL ] Iteration 33:\tk_eff = 0.748734\tres = 1.865E-02\n", - "[ NORMAL ] Iteration 34:\tk_eff = 0.761904\tres = 1.812E-02\n", - "[ NORMAL ] Iteration 35:\tk_eff = 0.774897\tres = 1.759E-02\n", - "[ NORMAL ] Iteration 36:\tk_eff = 0.787700\tres = 1.705E-02\n", - "[ NORMAL ] Iteration 37:\tk_eff = 0.800299\tres = 1.652E-02\n", - "[ NORMAL ] Iteration 38:\tk_eff = 0.812684\tres = 1.599E-02\n", - "[ NORMAL ] Iteration 39:\tk_eff = 0.824844\tres = 1.548E-02\n", - "[ NORMAL ] Iteration 40:\tk_eff = 0.836773\tres = 1.496E-02\n", - "[ NORMAL ] Iteration 41:\tk_eff = 0.848462\tres = 1.446E-02\n", - "[ NORMAL ] Iteration 42:\tk_eff = 0.859908\tres = 1.397E-02\n", - "[ NORMAL ] Iteration 43:\tk_eff = 0.871105\tres = 1.349E-02\n", - "[ NORMAL ] Iteration 44:\tk_eff = 0.882052\tres = 1.302E-02\n", - "[ NORMAL ] Iteration 45:\tk_eff = 0.892745\tres = 1.257E-02\n", - "[ NORMAL ] Iteration 46:\tk_eff = 0.903184\tres = 1.212E-02\n", - "[ NORMAL ] Iteration 47:\tk_eff = 0.913368\tres = 1.169E-02\n", - "[ NORMAL ] Iteration 48:\tk_eff = 0.923296\tres = 1.128E-02\n", - "[ NORMAL ] Iteration 49:\tk_eff = 0.932972\tres = 1.087E-02\n", - "[ NORMAL ] Iteration 50:\tk_eff = 0.942394\tres = 1.048E-02\n", - "[ NORMAL ] Iteration 51:\tk_eff = 0.951567\tres = 1.010E-02\n", - "[ NORMAL ] Iteration 52:\tk_eff = 0.960491\tres = 9.733E-03\n", - "[ NORMAL ] Iteration 53:\tk_eff = 0.969168\tres = 9.378E-03\n", - "[ NORMAL ] Iteration 54:\tk_eff = 0.977604\tres = 9.035E-03\n", - "[ NORMAL ] Iteration 55:\tk_eff = 0.985801\tres = 8.704E-03\n", - "[ NORMAL ] Iteration 56:\tk_eff = 0.993762\tres = 8.384E-03\n", - "[ NORMAL ] Iteration 57:\tk_eff = 1.001491\tres = 8.076E-03\n", - "[ NORMAL ] Iteration 58:\tk_eff = 1.008992\tres = 7.778E-03\n", - "[ NORMAL ] Iteration 59:\tk_eff = 1.016271\tres = 7.490E-03\n", - "[ NORMAL ] Iteration 60:\tk_eff = 1.023329\tres = 7.214E-03\n", - "[ NORMAL ] Iteration 61:\tk_eff = 1.030174\tres = 6.945E-03\n", - "[ NORMAL ] Iteration 62:\tk_eff = 1.036809\tres = 6.689E-03\n", - "[ NORMAL ] Iteration 63:\tk_eff = 1.043238\tres = 6.441E-03\n", - "[ NORMAL ] Iteration 64:\tk_eff = 1.049466\tres = 6.200E-03\n", - "[ NORMAL ] Iteration 65:\tk_eff = 1.055498\tres = 5.970E-03\n", - "[ NORMAL ] Iteration 66:\tk_eff = 1.061338\tres = 5.748E-03\n", - "[ NORMAL ] Iteration 67:\tk_eff = 1.066993\tres = 5.533E-03\n", - "[ NORMAL ] Iteration 68:\tk_eff = 1.072465\tres = 5.328E-03\n", - "[ NORMAL ] Iteration 69:\tk_eff = 1.077760\tres = 5.128E-03\n", - "[ NORMAL ] Iteration 70:\tk_eff = 1.082882\tres = 4.937E-03\n", - "[ NORMAL ] Iteration 71:\tk_eff = 1.087836\tres = 4.753E-03\n", - "[ NORMAL ] Iteration 72:\tk_eff = 1.092627\tres = 4.575E-03\n", - "[ NORMAL ] Iteration 73:\tk_eff = 1.097260\tres = 4.404E-03\n", - "[ NORMAL ] Iteration 74:\tk_eff = 1.101737\tres = 4.240E-03\n", - "[ NORMAL ] Iteration 75:\tk_eff = 1.106065\tres = 4.080E-03\n", - "[ NORMAL ] Iteration 76:\tk_eff = 1.110247\tres = 3.929E-03\n", - "[ NORMAL ] Iteration 77:\tk_eff = 1.114288\tres = 3.781E-03\n", - "[ NORMAL ] Iteration 78:\tk_eff = 1.118191\tres = 3.639E-03\n", - "[ NORMAL ] Iteration 79:\tk_eff = 1.121961\tres = 3.502E-03\n", - "[ NORMAL ] Iteration 80:\tk_eff = 1.125602\tres = 3.372E-03\n", - "[ NORMAL ] Iteration 81:\tk_eff = 1.129119\tres = 3.245E-03\n", - "[ NORMAL ] Iteration 82:\tk_eff = 1.132513\tres = 3.124E-03\n", - "[ NORMAL ] Iteration 83:\tk_eff = 1.135790\tres = 3.006E-03\n", - "[ NORMAL ] Iteration 84:\tk_eff = 1.138954\tres = 2.894E-03\n", - "[ NORMAL ] Iteration 85:\tk_eff = 1.142008\tres = 2.786E-03\n", - "[ NORMAL ] Iteration 86:\tk_eff = 1.144954\tres = 2.681E-03\n", - "[ NORMAL ] Iteration 87:\tk_eff = 1.147797\tres = 2.580E-03\n", - "[ NORMAL ] Iteration 88:\tk_eff = 1.150539\tres = 2.483E-03\n", - "[ NORMAL ] Iteration 89:\tk_eff = 1.153185\tres = 2.389E-03\n", - "[ NORMAL ] Iteration 90:\tk_eff = 1.155738\tres = 2.300E-03\n", - "[ NORMAL ] Iteration 91:\tk_eff = 1.158201\tres = 2.214E-03\n", - "[ NORMAL ] Iteration 92:\tk_eff = 1.160575\tres = 2.130E-03\n", - "[ NORMAL ] Iteration 93:\tk_eff = 1.162865\tres = 2.050E-03\n", - "[ NORMAL ] Iteration 94:\tk_eff = 1.165074\tres = 1.973E-03\n", - "[ NORMAL ] Iteration 95:\tk_eff = 1.167202\tres = 1.899E-03\n", - "[ NORMAL ] Iteration 96:\tk_eff = 1.169255\tres = 1.827E-03\n", - "[ NORMAL ] Iteration 97:\tk_eff = 1.171235\tres = 1.759E-03\n", - "[ NORMAL ] Iteration 98:\tk_eff = 1.173142\tres = 1.693E-03\n", - "[ NORMAL ] Iteration 99:\tk_eff = 1.174981\tres = 1.629E-03\n", - "[ NORMAL ] Iteration 100:\tk_eff = 1.176753\tres = 1.568E-03\n", - "[ NORMAL ] Iteration 101:\tk_eff = 1.178461\tres = 1.508E-03\n", - "[ NORMAL ] Iteration 102:\tk_eff = 1.180108\tres = 1.452E-03\n", - "[ NORMAL ] Iteration 103:\tk_eff = 1.181694\tres = 1.397E-03\n", - "[ NORMAL ] Iteration 104:\tk_eff = 1.183223\tres = 1.344E-03\n", - "[ NORMAL ] Iteration 105:\tk_eff = 1.184696\tres = 1.294E-03\n", - "[ NORMAL ] Iteration 106:\tk_eff = 1.186115\tres = 1.245E-03\n", - "[ NORMAL ] Iteration 107:\tk_eff = 1.187482\tres = 1.198E-03\n", - "[ NORMAL ] Iteration 108:\tk_eff = 1.188800\tres = 1.152E-03\n", - "[ NORMAL ] Iteration 109:\tk_eff = 1.190069\tres = 1.110E-03\n", - "[ NORMAL ] Iteration 110:\tk_eff = 1.191290\tres = 1.067E-03\n", - "[ NORMAL ] Iteration 111:\tk_eff = 1.192469\tres = 1.026E-03\n", - "[ NORMAL ] Iteration 112:\tk_eff = 1.193602\tres = 9.889E-04\n", - "[ NORMAL ] Iteration 113:\tk_eff = 1.194694\tres = 9.503E-04\n", - "[ NORMAL ] Iteration 114:\tk_eff = 1.195747\tres = 9.154E-04\n", - "[ NORMAL ] Iteration 115:\tk_eff = 1.196759\tres = 8.808E-04\n", - "[ NORMAL ] Iteration 116:\tk_eff = 1.197735\tres = 8.464E-04\n", - "[ NORMAL ] Iteration 117:\tk_eff = 1.198674\tres = 8.155E-04\n", - "[ NORMAL ] Iteration 118:\tk_eff = 1.199579\tres = 7.845E-04\n", - "[ NORMAL ] Iteration 119:\tk_eff = 1.200450\tres = 7.550E-04\n", - "[ NORMAL ] Iteration 120:\tk_eff = 1.201289\tres = 7.263E-04\n", - "[ NORMAL ] Iteration 121:\tk_eff = 1.202097\tres = 6.984E-04\n", - "[ NORMAL ] Iteration 122:\tk_eff = 1.202875\tres = 6.728E-04\n", - "[ NORMAL ] Iteration 123:\tk_eff = 1.203623\tres = 6.470E-04\n", - "[ NORMAL ] Iteration 124:\tk_eff = 1.204344\tres = 6.220E-04\n", - "[ NORMAL ] Iteration 125:\tk_eff = 1.205037\tres = 5.994E-04\n", - "[ NORMAL ] Iteration 126:\tk_eff = 1.205706\tres = 5.757E-04\n", - "[ NORMAL ] Iteration 127:\tk_eff = 1.206349\tres = 5.551E-04\n", - "[ NORMAL ] Iteration 128:\tk_eff = 1.206968\tres = 5.328E-04\n", - "[ NORMAL ] Iteration 129:\tk_eff = 1.207565\tres = 5.132E-04\n", - "[ NORMAL ] Iteration 130:\tk_eff = 1.208138\tres = 4.944E-04\n", - "[ NORMAL ] Iteration 131:\tk_eff = 1.208690\tres = 4.749E-04\n", - "[ NORMAL ] Iteration 132:\tk_eff = 1.209221\tres = 4.566E-04\n", - "[ NORMAL ] Iteration 133:\tk_eff = 1.209733\tres = 4.395E-04\n", - "[ NORMAL ] Iteration 134:\tk_eff = 1.210226\tres = 4.232E-04\n", - "[ NORMAL ] Iteration 135:\tk_eff = 1.210700\tres = 4.076E-04\n", - "[ NORMAL ] Iteration 136:\tk_eff = 1.211156\tres = 3.916E-04\n", - "[ NORMAL ] Iteration 137:\tk_eff = 1.211595\tres = 3.767E-04\n", - "[ NORMAL ] Iteration 138:\tk_eff = 1.212017\tres = 3.628E-04\n", - "[ NORMAL ] Iteration 139:\tk_eff = 1.212423\tres = 3.486E-04\n", - "[ NORMAL ] Iteration 140:\tk_eff = 1.212814\tres = 3.350E-04\n", - "[ NORMAL ] Iteration 141:\tk_eff = 1.213191\tres = 3.227E-04\n", - "[ NORMAL ] Iteration 142:\tk_eff = 1.213554\tres = 3.107E-04\n", - "[ NORMAL ] Iteration 143:\tk_eff = 1.213902\tres = 2.989E-04\n", - "[ NORMAL ] Iteration 144:\tk_eff = 1.214238\tres = 2.869E-04\n", - "[ NORMAL ] Iteration 145:\tk_eff = 1.214561\tres = 2.770E-04\n", - "[ NORMAL ] Iteration 146:\tk_eff = 1.214872\tres = 2.659E-04\n", - "[ NORMAL ] Iteration 147:\tk_eff = 1.215171\tres = 2.558E-04\n", - "[ NORMAL ] Iteration 148:\tk_eff = 1.215458\tres = 2.458E-04\n", - "[ NORMAL ] Iteration 149:\tk_eff = 1.215735\tres = 2.361E-04\n", - "[ NORMAL ] Iteration 150:\tk_eff = 1.216002\tres = 2.281E-04\n", - "[ NORMAL ] Iteration 151:\tk_eff = 1.216258\tres = 2.194E-04\n", - "[ NORMAL ] Iteration 152:\tk_eff = 1.216505\tres = 2.110E-04\n", - "[ NORMAL ] Iteration 153:\tk_eff = 1.216742\tres = 2.028E-04\n", - "[ NORMAL ] Iteration 154:\tk_eff = 1.216970\tres = 1.952E-04\n", - "[ NORMAL ] Iteration 155:\tk_eff = 1.217191\tres = 1.876E-04\n", - "[ NORMAL ] Iteration 156:\tk_eff = 1.217402\tres = 1.809E-04\n", - "[ NORMAL ] Iteration 157:\tk_eff = 1.217605\tres = 1.734E-04\n", - "[ NORMAL ] Iteration 158:\tk_eff = 1.217800\tres = 1.673E-04\n", - "[ NORMAL ] Iteration 159:\tk_eff = 1.217988\tres = 1.602E-04\n", - "[ NORMAL ] Iteration 160:\tk_eff = 1.218169\tres = 1.546E-04\n", - "[ NORMAL ] Iteration 161:\tk_eff = 1.218344\tres = 1.485E-04\n", - "[ NORMAL ] Iteration 162:\tk_eff = 1.218511\tres = 1.429E-04\n", - "[ NORMAL ] Iteration 163:\tk_eff = 1.218672\tres = 1.378E-04\n", - "[ NORMAL ] Iteration 164:\tk_eff = 1.218828\tres = 1.323E-04\n", - "[ NORMAL ] Iteration 165:\tk_eff = 1.218977\tres = 1.273E-04\n", - "[ NORMAL ] Iteration 166:\tk_eff = 1.219121\tres = 1.225E-04\n", - "[ NORMAL ] Iteration 167:\tk_eff = 1.219259\tres = 1.182E-04\n", - "[ NORMAL ] Iteration 168:\tk_eff = 1.219392\tres = 1.131E-04\n", - "[ NORMAL ] Iteration 169:\tk_eff = 1.219520\tres = 1.090E-04\n", - "[ NORMAL ] Iteration 170:\tk_eff = 1.219643\tres = 1.050E-04\n", - "[ NORMAL ] Iteration 171:\tk_eff = 1.219761\tres = 1.008E-04\n", - "[ NORMAL ] Iteration 172:\tk_eff = 1.219876\tres = 9.712E-05\n", - "[ NORMAL ] Iteration 173:\tk_eff = 1.219985\tres = 9.360E-05\n", - "[ NORMAL ] Iteration 174:\tk_eff = 1.220090\tres = 8.946E-05\n", - "[ NORMAL ] Iteration 175:\tk_eff = 1.220191\tres = 8.611E-05\n", - "[ NORMAL ] Iteration 176:\tk_eff = 1.220288\tres = 8.299E-05\n", - "[ NORMAL ] Iteration 177:\tk_eff = 1.220382\tres = 7.977E-05\n", - "[ NORMAL ] Iteration 178:\tk_eff = 1.220472\tres = 7.682E-05\n", - "[ NORMAL ] Iteration 179:\tk_eff = 1.220560\tres = 7.421E-05\n", - "[ NORMAL ] Iteration 180:\tk_eff = 1.220643\tres = 7.140E-05\n", - "[ NORMAL ] Iteration 181:\tk_eff = 1.220724\tres = 6.824E-05\n", - "[ NORMAL ] Iteration 182:\tk_eff = 1.220801\tres = 6.591E-05\n", - "[ NORMAL ] Iteration 183:\tk_eff = 1.220875\tres = 6.353E-05\n", - "[ NORMAL ] Iteration 184:\tk_eff = 1.220946\tres = 6.092E-05\n", - "[ NORMAL ] Iteration 185:\tk_eff = 1.221015\tres = 5.826E-05\n", - "[ NORMAL ] Iteration 186:\tk_eff = 1.221081\tres = 5.613E-05\n", - "[ NORMAL ] Iteration 187:\tk_eff = 1.221145\tres = 5.433E-05\n", - "[ NORMAL ] Iteration 188:\tk_eff = 1.221206\tres = 5.230E-05\n", - "[ NORMAL ] Iteration 189:\tk_eff = 1.221266\tres = 5.025E-05\n", - "[ NORMAL ] Iteration 190:\tk_eff = 1.221322\tres = 4.833E-05\n", - "[ NORMAL ] Iteration 191:\tk_eff = 1.221376\tres = 4.642E-05\n", - "[ NORMAL ] Iteration 192:\tk_eff = 1.221429\tres = 4.437E-05\n", - "[ NORMAL ] Iteration 193:\tk_eff = 1.221479\tres = 4.286E-05\n", - "[ NORMAL ] Iteration 194:\tk_eff = 1.221528\tres = 4.123E-05\n", - "[ NORMAL ] Iteration 195:\tk_eff = 1.221574\tres = 3.993E-05\n", - "[ NORMAL ] Iteration 196:\tk_eff = 1.221619\tres = 3.785E-05\n", - "[ NORMAL ] Iteration 197:\tk_eff = 1.221662\tres = 3.669E-05\n", - "[ NORMAL ] Iteration 198:\tk_eff = 1.221704\tres = 3.547E-05\n", - "[ NORMAL ] Iteration 199:\tk_eff = 1.221743\tres = 3.387E-05\n", - "[ NORMAL ] Iteration 200:\tk_eff = 1.221781\tres = 3.243E-05\n", - "[ NORMAL ] Iteration 201:\tk_eff = 1.221818\tres = 3.127E-05\n", - "[ NORMAL ] Iteration 202:\tk_eff = 1.221854\tres = 3.018E-05\n", - "[ NORMAL ] Iteration 203:\tk_eff = 1.221888\tres = 2.924E-05\n", - "[ NORMAL ] Iteration 204:\tk_eff = 1.221921\tres = 2.786E-05\n", - "[ NORMAL ] Iteration 205:\tk_eff = 1.221952\tres = 2.677E-05\n", - "[ NORMAL ] Iteration 206:\tk_eff = 1.221982\tres = 2.579E-05\n", - "[ NORMAL ] Iteration 207:\tk_eff = 1.222011\tres = 2.468E-05\n", - "[ NORMAL ] Iteration 208:\tk_eff = 1.222039\tres = 2.377E-05\n", - "[ NORMAL ] Iteration 209:\tk_eff = 1.222066\tres = 2.288E-05\n", - "[ NORMAL ] Iteration 210:\tk_eff = 1.222093\tres = 2.204E-05\n", - "[ NORMAL ] Iteration 211:\tk_eff = 1.222118\tres = 2.117E-05\n", - "[ NORMAL ] Iteration 212:\tk_eff = 1.222141\tres = 2.044E-05\n", - "[ NORMAL ] Iteration 213:\tk_eff = 1.222165\tres = 1.954E-05\n", - "[ NORMAL ] Iteration 214:\tk_eff = 1.222187\tres = 1.899E-05\n", - "[ NORMAL ] Iteration 215:\tk_eff = 1.222208\tres = 1.797E-05\n", - "[ NORMAL ] Iteration 216:\tk_eff = 1.222229\tres = 1.745E-05\n", - "[ NORMAL ] Iteration 217:\tk_eff = 1.222248\tres = 1.686E-05\n", - "[ NORMAL ] Iteration 218:\tk_eff = 1.222267\tres = 1.623E-05\n", - "[ NORMAL ] Iteration 219:\tk_eff = 1.222286\tres = 1.568E-05\n", - "[ NORMAL ] Iteration 220:\tk_eff = 1.222304\tres = 1.508E-05\n", - "[ NORMAL ] Iteration 221:\tk_eff = 1.222320\tres = 1.458E-05\n", - "[ NORMAL ] Iteration 222:\tk_eff = 1.222337\tres = 1.387E-05\n", - "[ NORMAL ] Iteration 223:\tk_eff = 1.222353\tres = 1.337E-05\n", - "[ NORMAL ] Iteration 224:\tk_eff = 1.222368\tres = 1.286E-05\n", - "[ NORMAL ] Iteration 225:\tk_eff = 1.222382\tres = 1.230E-05\n", - "[ NORMAL ] Iteration 226:\tk_eff = 1.222396\tres = 1.201E-05\n", - "[ NORMAL ] Iteration 227:\tk_eff = 1.222409\tres = 1.131E-05\n", - "[ NORMAL ] Iteration 228:\tk_eff = 1.222422\tres = 1.082E-05\n", - "[ NORMAL ] Iteration 229:\tk_eff = 1.222435\tres = 1.051E-05\n", - "[ NORMAL ] Iteration 230:\tk_eff = 1.222447\tres = 1.015E-05\n" - ] - } - ], + "outputs": [], "source": [ "# Generate tracks for OpenMOC\n", - "openmoc_geometry.initializeFlatSourceRegions()\n", "track_generator = openmoc.TrackGenerator(openmoc_geometry, num_azim=128, spacing=0.1)\n", "track_generator.generateTracks()\n", "\n", @@ -1733,21 +897,11 @@ }, { "cell_type": "code", - "execution_count": 30, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "openmc keff = 1.223729\n", - "openmoc keff = 1.222447\n", - "bias [pcm]: -128.2\n" - ] - } - ], + "outputs": [], "source": [ "# Print report of keff and bias with OpenMC\n", "openmoc_keff = solver.getKeff()\n", @@ -1790,7 +944,7 @@ }, { "cell_type": "code", - "execution_count": 31, + "execution_count": null, "metadata": { "collapsed": false }, @@ -1816,32 +970,11 @@ }, { "cell_type": "code", - "execution_count": 32, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "data": { - "text/plain": [ - "(9.9999999999999994e-12, 20.0)" - ] - }, - "execution_count": 32, - "metadata": {}, - "output_type": "execute_result" - }, - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEeCAYAAABlggnIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4VFX6wPHvTCadoNLFBiq+rh0RxIZgr2vFuq6yFhRF\ndBd7xQayVlBUrKirrgULyE+Fta2o7KLYWHlBBWwogkqA1Cm/P+6dySSZJDNJpub9PE8eJnfu3PfM\nMDnvPefce44nFAphjDHGAHjTXQBjjDGZw5KCMcaYCEsKxhhjIiwpGGOMibCkYIwxJsKSgjHGmAhf\nugtgsouIBIFuqvpr1LbjgfNVdViM/b3ABOAwIAgsAUaq6ioR2Qm4DygBQsCVqvqa+7rbgeOBcJxF\nqnpyE+X5AghEbf6vqp4jIguAfVW1PMH3eCRwgKqOSeR1LRzzdGAkUAwUAO8Bl6rqmvaKEWc5zgAu\nwPnb9wEfAH9L9DOKOt7hwCBVvS4Zn5tJPUsKJtn+AvQH+qtqrYjcCtwOnA48AVyjqq+IyPbAByLS\nRVX9wB7Aiar6YRwxhkYnqTBV7d+aAqvqDGBGa14bi4hcCRwCHKWqv4iID7jLjTGkveLEUY6BwDXA\nAFX93U3Y9+Ik5lNbediBQBdo/8/NpIclBZNsXwDzVbXW/f0jYJT7eFdVDZ/hbw38BgREpBAnkVwi\nIlsBXwEXq+p3TcTwxNoYbtXgnJk/DnR1n3pVVa8VkV5NbD8DOE5VjxSRTXEqzS3cONNU9TYR6QP8\nC3gV2B2nYrxKVZ9tUIZS4ApgF1X9BUBV/SJyCXC0iOQDV+EkwV7ApziJ9E5gP5wW0Dz3/a8TkfNw\nWhw1QBVOq+vLprY3+Eg2xukyLgV+V9WgiFwLbBdV3quAY939lgGjVHWF+1ndDwhOi+9+t1wjgTwR\nWYPz/9Qun5tJHxtTMEmlqh+q6icAIrIRcC3wrPtcQEQ8IvI18AJwq6qGgN44FcflqroL8CHwcjNh\n3hKRBVE/3aKe8wBnA1+r6gBgH2BrEenczHZwurMA/gH8S1V3AvYC/iQiJ7rP9QVeU9XdgcuAiTHK\nti1QoapfN/hcKlX16ahkuRlOa+rPOGfzvYCdgJ1x/k7/7p7Z3wkcrKqDgKnAXiKSF2t7jLLMAuYC\ny0TkIxGZDAxU1XcAROTPwA443UH9gf8DHnJfOwWnC+8POAnsHGA1TnJ4RlWvbufPzaSJJQWTqFjz\nonip36ffiHvG/y7wrqpOCW9X1ZCqboXTUrhCRIap6lJVPUJVl7j73AZsJSJbNHH4oaraP+pnVYPy\n/h9wnIi8inNme4Xbh97UdgCPiJQAe+J0seA+9xhwqHvcWlWd5e6/ALcbpYEgLf+dhYAPVTXo/n4I\ncL+qBtwkORk41H3+OZxutsnAGuARt7XVaHvDIKrqV9U/4SSg23FaUNNE5Bl3lyOAwcB8dzzmAmAb\n97n9cZINqlquqju6ic5D/ZZae31uJk0sKZhErcLpkonW092OiMyKOmM/wt02DHgfeFRVR7nb8kXk\nJBHxAKjqMmAO0F9EdhSR08IHd/fxALW0gqrOxzk7nQr0Af4jIns0tT3qpV4aV3p51HW71kRtDxG7\nG+t/QL6bFCNEpMj9rDZ2N62PETc6Zr77Xk7Dqby/wjnLnt7c9gYxzxSRI1X1J1V9SlVHArsCw0Wk\nqxt3Qji5ArtRN+bhb3CsviJS5r7vhicK7fG5mTTJuKQgIjuLyLsi8qiIDE13eUwj/wdcGK7M3S6h\nP+N0TaCqh0Wdsc8UkT2BF4HTVPWO8EHcbpMbgZPc4/QGhgFv41QUd7v9zwDnAZ+q6o+tKK9HRCbg\nDGi/DFwELAS2EZHxMbb3c+Ojqutwuq7Od8u4AXAaMJs4KzJVrQZuBR4RkR7ucQpxBpqLVXVFjGO9\nDpwrIj63y+h84A0R6Soi3wK/qurdON1MOzW1PUZx/MBEEdk8atu2wFKc8ZzXgbPdyh7gemCa+3gO\nMCLqc/iX+1n5cVoc0e+5zZ+bSZ+MSwrAIGAFzpdtYZrLYhobAxQBX4jIp8A7OH3Kjzex//U4leyt\nUS2IF9znjsGp/BbgXLUyVlU/VtUvgNHADBH5H3AU0OhyVFdz0/yGz2LvBHYRkc+B/wLfAE/hVMwN\ntz8T9VpwrsrZX0Q+wxlYfV5VpzXYp9myqOp4nDGT1933+glOt9JRDcoZdhPwk7vf/3DOsseo6mr3\nuX+JyHxgPHBWU9tjlGMaMAmYKSIqIotwKu5D3K6ph4CZwIci8gXOeMbp7ssvAP7g/p+/B9yiqh/j\nJIc/isjdDd5Hmz83kx6eTJs6W0T+AHyPM9B2lqpeluYiGWNMh5HSS1JFZHecPsthbrN4Ck4ztxon\nAXwN7ILTUvg91eUzxpiOLmXdRyJyKfAgUOhuOhooUNU9gctxroYA59royTj9sJNSVT5jjDGpPRP/\nCuemmCfc3/cGXgNQ1Xkispv7+AOcW++NMcakWMpaCqo6nfqXtZUB0fOtBNwuJWOMMWmSzj77cpzE\nEOaNunknLn5/IBQIpG6g3Ofz4vcnVESL1UFipTqexcquWKmO11KswkJfk5cGpzMpzAWOBJ4TkcHA\nZ4keIBAIUV5e2e4Fa0rnzsUpi2exsitWquNZrOyKlep4LcXq3r2syefSkRTCp/YvAgeKyFz39xFp\nKIsxxpgoKU0K7lQGe7qPQzh3qhpjjMkQNrBrjDEmwpKCMcaYCEsKxhhjIiwpGGOMibCkYIwxJsIm\nnDPGGOCbb77m/vsnU1VVRWVlBYMH78WZZ45M6BjvvPMW22+/I16vh0cffYi//S37Jnm2loIxpsNb\nu3Yt48ZdxZgxY5k06X4eeOAxvvnmK15+udECds16/vlnqKhYR5cuXbMyIYC1FIwxhvfee4cBAway\nySabAuD1ernmmhvx+XxMnnwnn3/+KQAHHngIw4efxM03X09BQQErVqxg9epVXHXVdaxatYolSxZz\n003Xc801N3DTTdfxwAOPcvrpJ9G//wCWLfuGQCDIhAm3o7qIl1+ezrhxtwDwxz8ezCuvvM6KFT8y\nfvwNBIPOFBUXXXQJW2/dL/I8wHXXXcHRRx9P167dGD9+HHl5PkKhENdddxM9evRs82dhScEYk3GG\nDClh0aK8djvettsGePfdiiafX7VqFRtvvEm9bUVFRcyd+29++ulHpk59DL/fz6hRZzFgwG54PB56\n9erNJZdcyYwZL/HKKy8yduwV9Ou3DZdcciU+X13VWlFRwQEHHMKeew5i7NixfPjh+3Tp0rVeLI87\nE9G9997FCSecwt57D2HJksVMmHAjDz30eOR5d28A5s//D9tttyPnnTeazz77hHXr1llSMMbkpuYq\n8GTo1asXixdrvW0//vgDixcvYued+wPg8/nYfvsdWbp0KQDbbCMAdO/eI9KSaEp43x49elJTU9Po\n+fACmMuXL2OXXXYFoF+/bVi58ucYRwvh8Xg44oij+Mc/pvG3v11Ip06ljBx5ftzvtzk2pmCM6fD2\n2msf5s17nx9++B4Av9/PPffcRVlZGZ999klk2xdffMpmm23W5HG8Xm+k6yeap/6pPgUFhaxevQqA\nn35aQXn5GgC22KIvn3zyMQBLlihdu3aNxK6srKS2tpalS78hFArx73+/w8479+fuu6cwdOj+PPnk\nNNqDtRSMMR1eSUkpV101jokTbyYYDFJRUcHeew/h+ONP4ueff+bcc/9CbW0t++9/INtssy1QV9FH\nV/g77LATN998HZdccmWjRBBt223/QFlZGeeccwZ9+vSld2+n6+qCCy7i1ltv4plnnsTv93P55dcC\nMHz4yYwceQa9e29Cr1698Xg8bLvtH7j55uvJz88nGAxy4YV/bZfPwhMKpW49gvZWXe0PddSpby1W\nZsVKdTyLlV2xUh0vjqmzm8xYWd199M036S6BMcbklqxOCvvvn8f11xdSVZXukhhjTG7I6qTwn/8E\nWL7cw0EHlfD551n9VowxJiNkdU3avTs88kgVo0fXcOKJxdxxRwF+f7pLZYwx2SurkwI4N30MH+5n\n9uwK3n8/jyOOKOGrr5oe9TfGGNO0rE8KYZtsEuLZZysZPryWI44o4cEH84lxubAxxphm5ExSAPB6\n4cwza5k1q4Lp0/MZPryY77+3VoMxpnkffzyfffYZyL/+9Ua97aeffhK33DIu5mtmzZrB/fffA8DL\nL0/H7/ezZMliHnvsoZj7f/jhB4wZM4pRo85i9OiR3HLLONavX9e+b6Qd5FRSCNtyyxAzZlQwZEiA\nAw8s4ZlnfGTx7RjGmBTYYos+zJlTlxS+/vorqpq5tNHj8URuUHvyyccIBoP067cNZ5xxVqN9lyxZ\nzJ133sE119zAlCkPMXnyA2y9dT+eeuqJ9n8jbZSzdzT7fDBmTA377+/n/POLmDXLx+23V9O9u2UH\nYzJZ8ZTJlPx9PN52PIsOlnai4pIrqBw1OubzHo+Hrbbqx3fffcv69esoLe3E66/P4qCDDuXnn3/i\nqKMO5uWX689SChAKhZg582VWr17N9ddfxfDhJ/HSSy9EZj8Ne/nlFxg58ly6desW2XbCCadEHp92\n2glsvvkW5OcXMHbsFdxww9VUVFQQCPg5++xR7Lrrbhx//JE8/fR08vPzue++yfTp05devTbmmWee\npKamhl9//ZVjjjkuUrbWysmWQrQddgjyxhsV9OsXZNiwEl59NWfzoDE5ofi+ye2aEAC869dRfN/k\nFvcbOnQ/3nnnLQAWLfofO+ywE86sD9Hd0PW7pI844ii6du3KuHG30NQMEStWrGCzzTYHnIn2Ro8e\nyQUXnMOoUU6roqqqijPOOJvrr7+Zxx57iEGDBnPPPVO58cZbGT/+Ridq1LQZ0Y/XrFnDbbdN4oEH\nHuXpp5/kt99+a/F9NifnkwJAYSFcc00NDz9cxbhxhVxwQRFr1qS7VMaYWCrPG02wtFO7HjNY2onK\n82K3EoBIZX7AAQczZ84bfPLJx5HZUd09mngc22effcLo0SMZPXokH3zwHj169OT7778DoHfvTZg8\n+QHuuOMefvllZeQ1m2++BQDffruMnXd2Zkrt1q07paWl/PbbrzHLC7DLLruSl5dHUVERW265FStW\n/NBi+ZrToU6bd989wJtvrueGGwoZOrSUu+6qYt99A+kuljEmSuWo0U128yRb796bUFVVyfPPP8O5\n546OzJoaCASorKzE5/OxdGnj+XU8Hg/BYF1dstNOuzB58gOR37t06cadd06gT59+dO3qdCF99NF/\n653xe73OOfoWW/Tl008/pl+/bfjll5WsW7eWzp03oKCggFWrfqFXr41ZsmQxffr0BZwWDTitjWXL\nlrLpppu36TPoUEkBoFMnmDixmjff9DNmTBGHHurnmmuqKSlJd8mMMekSPWi8//4H8vrr/8emm27G\nDz98j8fj4fjjT6o3S2n06wB23rk/Y8eOYcSIs2POjiqyLX/961huvvn6yDTYPXr04Kabbg0fKbLv\naaeNYPz4G3j77Teprq7i0kuvIi8vj1NO+TOXXDKGXr02pnPnzpH9169fz0UXjWLt2rWMGHFOveda\n9Vlk4iypItITmKmqA5vbr62zpP7+O1xxRRGffJLHPfdUMmBA8zc2ZNIshxYrs2KlOp7Fyq5YyYr3\n8cfzeeedN7n44ksTipVVs6SKiAe4BFiW7Fgbbgj33VfFFVdUc9ppxYwfX0CMRZGMMSYjRbdw2kvG\nJQXgXOBJIGVzn/7xj37eequChQvzOOSQEr78MhM/FmOMqa9//wFcdNEl7XrMlNZ+IrK7iLzlPvaK\nyP0i8r6IvCUiW7m7HQCMBAaJyHGpKlvPniGeeKKSs86q4dhji7nnnnwCNgZtjOlgUpYURORS4EGg\n0N10NFCgqnsClwO3A6jqcap6HjBPVV9IVfnAmVzvlFP8vPZaBbNn+zj66GKWLrVpMowxHUcqWwpf\nAcdSN8y+N/AagKrOA3aL3llV/5zCstWzxRYhXnyxksMO83PooSVMm5Zv02QYYzqElCUFVZ0ORK92\nUAaUR/0eEJGM6cz3euG882p5+eVKnnwyn5NPLubHH9NdKmOMSa503qdQjpMYwryqmtBk1z6fl86d\ni9u3VA0MHAjvvRdiwgQvgwd7uO22Ek44IfnNhlS8N4uVvfEsVnbFSnW8tsRKZ1KYCxwJPCcig4HP\nEj2A3x9M2XXGY8bAYYcVc8YZHl54IcStt1bRpUvy4uXq9dq5GivV8SxWdsVKdbw47lNo8rl0dNeE\nT7NfBKpEZC7OIPPFaShLQgYMgDlzKujVK8TQoaXMmZOX7iIZY0y7SmlLQVWXAXu6j0PAeamM3x6K\ni+HGG6s55BA/F15YxNChfsaNq6ZT+87fZYwxaZExA7vZZq+9Arz99noCARg6tJQPPrBWgzEm+1lS\naIOyMrjrrmpuvrmKc84p4rrrCmlmoSZjjMl4lhTawcEHB3j77Qq++87DgQeW8Nln9rEaY7KT1V7t\npGvXEA8/XMWYMTWcdFIxt99egN/f8uuMMSaTWFJoRx4PHH+8nzlzKvjwwzwOP7yEJUvsIzbGZA+r\nsZKgd+8Qzz5byUkn1XLkkcVMnZpPMKHb8owxJj0sKSSJxwMjRtQya1YFL72Uz/HHF/Pddza5njEm\ns1lSSLIttwwxY0YFQ4cGOOigEp55xmeT6xljMpYlhRTIy4MLL6zhuecquf/+Ak4/vYiVK63VYIzJ\nPJYUUmiHHYK8/noF22wTZNiwEmbOTOfUU8YY05glhRQrLISrr67h0UcrufHGQs4/v4g1a9JdKmOM\ncVhSSJNBg4K8+eZ6OnVyJtd7+22bJsMYk36WFNKotBRuvbWaO+6o4qKLirj88kLWr093qYwxHZkl\nhQwwbJgzud6aNR7237+U+fPtv8UYkx5W+2SIDTeE++6r4qqrqjn99GKuvdZLTU26S2WM6WgsKWSY\nI4/08+abFSxcCAcfXMLChfZfZIxJHatxMlDPniGefz7I2WfXcPzxxUyaVEAgkO5SGWM6AksKGcrj\ngVNO8fP66xW8+WYeRx1VzNKldsObMSa5LClkuM03DzF9eiVHHunn0ENLeOyxfJsmwxiTNJYUsoDX\nCyNH1vLKK5U89VQ+J51UzE8/WavBGNP+LClkkW22CfLqqxXsumuA/fazaTKMMe3PkkKWyc+Hyy6r\nYdq0Sm64oZAxY4pYty7dpTLG5ApLCllq4EBnmgyvN8SwYaX85z/2X2mMaTurSbJYp05w553VjBtX\nzYgRxUyYUEBtbbpLZYzJZpYUcsBhhzk3vH3ySR5HHFHC11/bILQxpnUsKeSInj1DPP10JSecUMvh\nh5cwbZpdumqMSZwlhRzi8cCZZzqXrj7+eD6nnVbML79Yq8EYE7+MSwoiMkBEHhWRx0SkR7rLk422\n2SbI//1fBdtuG2DYsBLeeMPWajDGxCfjkgJQCFwEvArskeayZK2CAmeFtwcfrOKKK4oYO9bWajDG\ntCzjkoKqvg9sB4wFPklzcbLeHnsEeOut9VRUOGs1LFiQcf/lxpgMktIaQkR2F5G33MdeEblfRN4X\nkbdEZCt3+0BgPnAo8NdUli9Xde4MU6ZUcfnl1Zx6ajF33GGzrhpjYmsyKbiV9mgR2dH9fYyIfCEi\nj4tI50QDicilwIM43UMARwMFqroncDlwu7u9E/AI8HfgH4nGMU07+mg/c+ZU8N57eRxzTDHff2+D\n0MaY+pprKYwHDgTWichewI04ff0fA5NaEesr4FggXBPtDbwGoKrzgN3cx2+p6mmqeqaq/qcVcUwz\nevcO8dxzley/f4CDDiphxgybP8kYU6e5GuFwoL+q1orIGOA5VZ0DzBGRRYkGUtXpItInalMZUB71\ne0BEvKoajPeYPp+Xzp2LEy1Kq6UyXrJjXXMNHHxwkNNPL+KDD+DWW4spLU1auIhc+gzTGc9iZVes\nVMdrS6zmkoJfVcOTJgzDaTmEtcdYRDlOYogcM5GEAOD3Bykvr2yHosSnc+filMVLRaxtt4U5c+Dq\nq0sZPNjDAw9UscMOCf0XJCzXPsN0xbNY2RUr1fFaitW9e1mTzzVXuVeIyBYisj2wLfAGgDvGsKZ1\nRa1nLnCYe8zBwGftcEyToLIyePTRIGPG1DB8eDFTp9qd0MZ0ZM0lhSuBD4B5wDhV/VVERgGzgWvb\nEDNc5bwIVInIXJxB5ovbcEzTRiec4GfWrAqmT8/n1FPtTmhjOqomu49U9W0R6QuUqOpv7uaPgCGq\nurg1wVR1GbCn+zgEnNea45jk6Ns3xIwZFdx6awH771/C3XdXMWyYXbtqTEfS3CWp56tqdVRCCF8l\ntFJEnk5J6UzK5ec7d0Lfc08VF11UxPXXF1JTk+5SZa7lyz0sWtT8EFsoBN98Yy0vkx2a+zYfLCLT\nRWTD8AYRGYbT928TJuS4IUMCvPlmBV9/7eWww2w67qYcfXQJQ4Y0f9nWRx95GTy4U4pKZEzbNHf1\n0VHA34D5InI2zqDwacC5qvpSKgrXkoJuG9E9xWtRdu9Asbrj3kgCCc9CFSztRMUlV1A5anTbCpbh\ngnFcrFVRYQnVZI8mWwqqGlLV24DrgH8BpwK7ZUpCAPDY4sQZy7t+HSV/H9/yjlnOa1NJmRzT7Fda\nRC4G7sQZEH4bmC4i/VJQrriEOlmTPJN51+d+0s6LY1ZyjzUUTBZpsvtIRP4FFAODVfUb4AERGQH8\nW0SuVtWHUlXIptSs+q3D3nySrlihEEyZks+UKQVMmlTF/vs3vjqpe4+Ep8bKWvFU+JYUTDZprqXw\nDrC3mxAAUNVHgX2BUckumMlMHg+cf34tDz9cxV//WsSECQVx9avnqrZ2H9XUgIi1eE3maG5M4YZY\n006oqmKL33R4gwcHmD27gvffz+OMM4roqMM7bW0pVFTAb79ZU8Jkjlad56hqdXsXxGSfHj1CPP98\nJd26hTj88BKWL+94lZt1H5lcY9dOmDYpKIDbb6/mz3+u5bDDSpg7t2OtB+3x1J8o6t//zuPOOwvS\nVBpj2s6SgmkzjwfOPLOW++6r4uyzi9JdnJRq2Aq4664Cxo8vjL2zMVmgxRVWROQM4DagS9TmkKp2\nrFNC06IhQwLMnFkBg9NdktRpmBRidRU1131kXUsm08Sz7NZ1wFBgoTuJnTFN2nLL+l+RmhqniylX\nxZMUmmPTlJtME0/30feq+oUlBNMap5xSzNq16S5F8rS1pRBOCpYcTKaIp6XwkYg8j7PITviqo5Cq\nPp68Yplc8c67+bBV/W3R8yxl+xxJ7dVSCAbjuzvamGSLp6WwIbAO596Eoe7PsOQVyWS7YGn8N2Pl\n2hxJid7MFgo5WaThDYDPPuuzq5hMWrTYUlDVM0SkABB3/y+i1m42ppGKS66g5O/j4577KJvnSGrr\nQHE4GTTsPpowoZDvv/dy8cW2mIVJrRbPa0RkN2AxMA14BFjurqlsTEyVo0azeumP/LKyvN7Pc8+u\np0f3II88XMEvK8vTXcx20dakYGMKJtPE09idBJyoqruqan/gWHebMQkZOjTAzJkBrr66kKlT89Nd\nnHbRMCk0V7kPGtR4MZ7w/iKdeOqpeIb4jEmueJJCqbsMJwCq+iHQse5QMu1m551hxowKpk3LzaTQ\n3D7LltX/c1P1suOOzvhLRYWHefOcpLBwoZfvv7f7Sk16xPPN+01Ejg7/IiLHAKuTVyST6zbfPOTc\n5NaCYBBWr87su7va0n3088+xXzxsWPPLexqTTPG0V88BnhSRhwEP8DXwp6SWyuS8jTaq/3tTazCU\n0InXBl/N0Fcyc7b2tiSFhl1NDedRMiYdWmwpqOpiVR0EbA5soaoD3emzjWmTeC5dLWMdh867iXff\nzcyL+BMZU2hJ+PJUY9KpyaQgIg+6/74lIm8BM4FX3N/fTFUBTe6quOSKuBJDp9A6Jk/OzGv242kp\n2JVFJps01310v/vvuKht9vU27aZy1Ohm72SO7lJauNDLkiVe+vXLrGXe2nNCO+s+MpmguZXXPnIf\nfgD8pqpvA5sCR+Dct2BMypx6ai2PP555Vyy1Z0vBuo9MJojn6qMngeNFZHfgeqAc50a2pBCR/UVk\nqog8KSI7JSuOyS6nnFLL88/7qE7zmn+//w7/+Ediyamt3Ue33VbA77/Ht+8331hiMW0TT1Loq6rX\nAMcBD6vqjcBGLbymLYpV9RycNRwOSmIck0X69g2x3XZBZs1K7w1eL7yQz8UXJ3abTnRS6NGjjOee\ni/0emuo+mjixkLffbvyaQAAqK+tvGzy4E4sW1f1ZL17s5dFHM6+FZTJXPEkhT0S6AUcDr4rIxkBJ\nsgqkqjNFpBS4EHgsWXFMduneozPv/jufc0aW0L1H50Y/Xfv2pnjK5JSXK54J8Bq2FBYujH0lVaLd\nR7fcUsAWW5Q12h7dmrrrrgIuu8zuNTXxiycp/B2YB8xS1c+Bt4EbWxNMRHZ3r2RCRLwicr+IvO9e\n0bSVu70bMBm4VlVXtSaOyQ2ZONtq46mynRp/6dKmK/SGM6C21+D0kiWx/3yjj29XPplExXOfwlOq\nupWqXiQiGwDDVfWZRAOJyKXAg0B4AdujgQJV3RO4HLjd3X470BMYLyLHJRrH5I54L1kNizXb6kMP\n5bNgQfKmjAhXwG+80XS3VlM3qSVy89o55xSzaJG3UXdRc2UypjXimSX1LBF5RER6AAuBZ0Xk5lbE\n+gpnMr3wV3Zv4DUAd26l3dzHp6vq4ap6mqq+0Io4JkfEmm117nvr6NUzwI8/1G1rzpVXFnHXXe13\nj0NTZ95NTYHd1LYePcr44gtvg/08TY43AAwZUsoddzR+LzvvXH9ajDPOKOaII4qbPI4xzYln1G4U\ncADO1BYvA2NwupOuSiSQqk4XkT5Rm8pwrmQKC4iIV1XjvhDd5/PSuXPqvvypjGexYtttN9hqKw/v\nvFPCscc2rm2jj+33O5XuBhvkJRxzzRrYYIPG2wsKnHOa9euL2Xhj8Pmc8YHi4nx8Pg8+X+NyFBfX\nP3UvKnIGfr/6qn4Fn5+fx/nn1x8UbljuQCCfzp3z8Pm8kVgrVtT/TL/7zst33znb8vO9MY+TiGz6\nfmRqrFSy3J96AAAgAElEQVTHa0usuC7lUNVfReQwYLKq+kWkPUauynESQ1hCCQHA7w9SXh5He7qd\ndO5cnLJ4Fqtp552Xx003FbLffpV4vfWX94w+9tdfFwNe1qwJUF5elVCMnj3L+PLLdXTtWj/xVFTk\nA3n07etj+fK1hELFgI+qqlr8fh9+P4CXFSsqKXVP4M8+u/6ZfHV1LVBIbW2A6Ma63+8H6icK5/2U\nRb3WT3l5NZ07F+P3hyKvr3vfZfVeW1tbBHgZPjzEww8n9hmEZdv3IxNjpTpeS7G6d298gUJYPJ2t\nC0VkJs5Ku7NF5Fngv4kWMoa5wGEA7qI9n7XDMU0HcNBBAQoK4NVXmz+n+fxzDz17BlmzpnWd7Oti\nLAgX3RUUCNT134f/DT/ft28ZP//s4b//9fLTT/X/zF5+Ob/RsQCCwfYdDFizpu7xjBl2WaqJTzxJ\nYQQwEdhdVWtwblw7qw0xw38KLwJVIjIXZ3D54jYc03QgHg9cfnk1N91U2OzNbJ9/7mGvvQL8/nti\nlW24so41YBsI1D0OBptfOe3dd/M4/PDG02B/+23sP7vp02PfixCrbA29917jy1z79Surt38wCN9+\na6PQpnlNnmqJyEhVfYC6sYNhIhJ+egBwQ6LBVHUZsKf7OAScl+gxjAHYb78AIgHuvbeApi5E/ewz\nOOKIAPPmJTbDaq27AnnDChnqn80HAnUDzMEgfPcdvPNO3Z/U008ndnZeU9O4wm54OetDDxWw444B\n9tij/vZjjy1h5cq1jV4/fXpdGZ591seFFxbH3M+YsHhaCrFOLex0w6TdzTdXN7msZyAACxZ42Gcf\nP2vXJvZ1ralx/m1YIUP9M3W/3xPZZ4MN4Lff6u/73nvJuft6zJhiBg1qfOxffmn+fZaX25+taVlz\n39oFAKp6fWqKYkxiNtssxNVX18BfGz+3eLGXHj2cVd7WrXMq83iv3w+3FPx+Dw0nBm7YHRMMeth2\n2wB+f+L3B8RzY1kiN59tv33z93TYjWwmHs21FB4IPxCR25vZz5i0OfXU2pjbZ8/2MWxYCJ8Piopg\n/fr4jxnuxnGuJKovuvUQHlMoKHD2TbTSbWtSeO21tg8eB4Oxu8lMxxXvrZ77JbUUxrRSw7Pzmhqo\nqoInn8zn5JOdGrysLJRQF1JzYwrRiSIQcH4KC+te094eeyx5Vw2tXw9/+1sh220X/13jJveld8pJ\nY9rZXnuVUloaYsCAAHvu6aG8vC4pbLxxy6fmixZ5WbXKSSCxkkJ05e90H0FhYYja2sT76+Ppbvro\no+QtQzpoUCm//JK8KUBMdrKkYHLK1KmVrFnjYZ99AoBzR2dZGayN84KbIUPqLiGN1X1UXV3/6qNw\nSyHWvi2ZObPlP79Yg92tFd0VVVNDvYTw+OP5DBvmZ7PNbOCho2vuW7mziES+ktGPgZCqZuZK6qZD\n69+/cS3aqVNi3UdhgUDj10Tf0BYK1bUUWpMUohNMU5I1OHzKKfWnQBg7tohTTqnhrrvSvIqRSbsm\nk4KqWrvS5IRExxTCYnUfVVVFtxQ8+P3OQHZrrj5KtTvvdCYofvppH+++G/tPf/lyDwMHdrJ7GTow\nq/hNzisriz1lRUtiJYXwPQwQvnLHQ2EhLF/upap1UwulzG+/OVlrzJjYE6WFQh6WLg3PpZSyYpkM\nY2MKJud17tx+LYXoaTUCAdyWQoinny7g99/bcQDAla57C3baqRPLlrUik5qsZy0Fk/Na230Ua5A3\nehwgfI1/obts1MqV7d9/lMqkEAzC5MnOLK0VFR4qKpzts2d76q37bHJbiy0FEfEC5wL7u/u/iTOF\ndvufFhnTRt17dK7/O3VL+nFry6+vVwef7CwJWnHJFVSOGg3E6j6C4uJQ5Pf2lsoby/75z/r3RPTp\nU8bKlWs58sg8dtqpiDlzKlJXGJM28aT/icBBOLOjPopzI9sdySyUMYlIZMnORDVc+7lhS8Hvd+5o\nhuSc1bfHXcvGJCKepHAQcJyqvqKqLwHHAYckt1jGxC/RtZwTFb32c8MxhWDQGWSG1t2rYEymiWeg\nOc/dL9yQ9QH29TcZo3LU6Ej3TrTw6lOzZ+fxyCMFPP1086teBYPQq1fdilShGJMBN+w+8vs9HHts\nLc88k5/TSeGzz/J488089tvPJkrKdfEkhX8Ab4vIUzhTZp8MPJ3UUhnTjpq6o3nVKg9duoTwuu3l\niji6zBvepxAIOFc3QW62FE49te7y1Wefzbek0AHEO6ZwI7CF+3OTqt6c1FIZ046auvpou+06MWlS\n3ZrIv//uweNpfmAgeu6jUMjpQiopcX7PxdlGZ8+uO2+Mvg8jFxOgccTTUviPqu4KzEp2YYxJhrKy\nEOvWxb5c9Icf6ravXOlh881DLF/e9KWl0d1H4bmPwlcf5XpFOWtWPoFAFbNm+TjzTFvBLVfFkxR+\nFpEhwDxVtYlRTNbp3j3E6tXOdffhs/qw6CuGfvnFw+abByMDx/WO4V7q+kv0xmOhBmCgeynryvYt\nd0baGP6C80OP1h2i4WW+JrPE0320G/A2UCkiQfcnBxvKJlcVF8Mf/hBkwYLm53BcudLLZpvV3WxQ\nU2jrDCRDw8t8TWZpMSmoandV9boT5PncxzZDqskqAwcG+M9/WkoKHnr2dJoOeXkh5h5wVVIvde3I\noi/zNZmlxaQgIsNEZG7dr7JURPZKcrmMaVeDBsWXFLp3d5LCBhuEmLf3Raxe+iO9egbwEGLFj+V4\nPUE8hPAQ4pGHKygqDLLy53I8hPDl1T3XEX4uu7Qq8viXleVx/ZjMF0/30R3ASABV/RI4FLg7mYUy\npr0NGhRg/vy8yFQU4UHh6DGFNWs8bLhhiIkTqzjwwEBk3/CU2OvXO11RYbW1kJdX93z0+gtdu+b+\nLDATJxamuwgmCeJJCoWq+kX4F1VdhM2uarJMjx4hNtooxMKFzlc+fE/C+vV1Ffm6dR46dQpxxhm1\nbLBBKHKJabjS79evjIqKuv1rapykEEtT243JdPFU7ioitwJP4Ny8dhKwOKmlMiYJhg+vZfLkAh54\noCpyiWr0DWvr1kEndwjB66VRUmiottbTZOWf38GmLKqsrN+KMtkrnpbCmUAnnLuYpwGlwNnJLJQx\nyXDeeTUsWuRl8uSCyIIzsVoK4Jzph7uDvFF/JY88UjdVRk0N+Hyxb3bzNXG61a1bbnYrbbFFWcs7\nmazQYktBVX8Fzk9BWeoRkf2Ak1XVEpBpF506wVNPVXLyycU8/7yPbt2C9bqDopOCz9e4+6i4OMQR\nR9TdoVZbW5cw7r+/knPPdU6Vd9ghwOGH+7n11sZ97pm+ZKcxTSYFEVmgqv1FJNapTSiZl6WKyFbA\nLkBRsmKYjmnTTUO8+moFU6cWsOWWQSZMqKu4o7uPnJaC8zhckZeV1W8V1NR4Ii2C6LUUHnywkv/9\nzwYVTHZqMimoan/335QvuaSqXwN3iMgTqY5tcl/nzjB2bA1+P/ztb0WsXu2ha9dQvZZC9JhC2AYb\n1E8K4auPAPbbz2lBDBgQYKutQnz5ZezY1lIwma65lsKfm3uhqj7emoAisjswQVWHuau6TQF2AqqB\ns9yEYEzS+Xxw6KF+nnoqn1GjaqisrJsGIy+v/jxH4CSTaNXVdQPKXbrAjBkBOnWqirzemGzUXCvg\nMeDvwMHAsBg/CRORS4EHgXCb/WigQFX3BC4nauVEY1Lh/PNrmDo1n9WrPZSU1I0RxOo+Ck+RHVZZ\n6aGgoG7bgQeG2H77oPua+vs++GDzazkYkymaSwq7Ao8AgjPf1zM4Z/IjVHVEK+N9BRwLkdVL9gZe\nA1DVeTjzLEWo6mmtjGNMXLbfPsguuwR58MH8SNcRhLuPnK9pU0mhqqrpS09Dofr9REcdVTdAHZ5V\n1ZhM1NyYwifAJ8AVIjIQOBG4RUT+C/xTVd9KNJiqTheRPlGbyoDoe98DIuJV1biu2/P5vHTunLqL\no1MZz2KlLtb553s444wCunYlsn9pqYe1az107uwlL8+p4Lt2zat3vGDQR0mJJ7ItOl5hYf2kEN6e\nl+fhm28CXH65l2nTUj5cl1SJ/r/G2j8Tvx/ZGK8tseK6M1lV/ysi84F9gAnAaTj3K7RVOU5iCIs7\nIQD4/UHKy1PXLA8v72ixcivWzjvDr7+W0a1bILJ/bW0+lZVeysurCYVKAQ/FxX7Ky6sJf2XXrg3g\n9Xoir4mOt2aND6j7o3S2lxEKhcjLq6SoqBCoW+AnF8TzWXdvYf9M/H5kY7yWYnXv3vR9Jc0mBXcg\neAhwPM6cR58Ck4CZrSloDHOBI4HnRGQw8Fk7HdeYuIXvxP3667oz93jGFJrrPopepcyYbNLc1Uf3\n4wwyLwCeBS5X1faa7zb81/UicGDULKytHaswpk2GDvXzxRd1SaH+NBfO17VxUvBQ0MTJ/gEHxF5y\nJJxgvLnVc2RySHMthXOA1UB/92e8iISfC6nqlq0JqKrLgD3dxyHgvNYcx5j2dM89VVRHrSvo88XX\nUthoo9iDxuEpuJsST1KYO3c9e+3VHr20qbFsmYc+fWwQPds1lxRaVekbk4169KhfmcXqPmp485rT\nUmhdJej11r1u/Pgq9t3Xz3ffeTnxxBLuuquSY47xZ90Ec4sWeenTxxZlzHbNXX20LIXlMCaj5OWF\nIpekhtdc2HDD+gmgupomu49aEt1S6NIlxNZbh9h6a6dCDQY9WZcQoPFluCY7Wc+mMTFEjykEg05l\n16VL8zevxSPWmELIelxMBrGkYEwMeXl1k9yFk8OGG9bfJ3qai0SOC7k5B9Lzz/tYtMiqlGxn/4PG\nxBA9plBY6JzKd+0a/9VH0ZYsWRt5HJ5VNTopNGwpZGvCmDEjn0cf7WCrC+UgSwrGxBCdFPLzYc6c\n9Y0q63hbChtsUPc4P79uFtawXOk+2mijEG+84WPFCg+LF3v55hsPN95YwEcfWTWTTWytZWNiiB5T\nWLPGU+9y1AUL1vHggwVMmVKQ8JhCrJZCtClTKjnkEH/sJzPc2LHVLFniZciQUrp0cRYp+vZbL5Mn\nF9KvX4C5cytaPohJO0vhxsTgjCl4CIVg9WoP3brVVf6bbBJi442dAYfmWgp/+EOAww+vrbctvH9T\n3UfHH++PLPSTbQIBmDixmn/+s4Kff/aw5ZZ1M9YsWWJziWcLaykYE4PPF8Lvh6VLPZSVhShtcA9Z\nOEk0N6bw5psVjVoE4W6jbB03aM733ztvbtddg/zvf+tYvdrDgAF1Ge7bbz315j4ymclaCsbEEB5T\nePddH0OGBBpV4r17O0mhuWmw8/Ia37kcbhU0N9CcraIXJSopgc02C/HFF3Uz4zz0UG5NAJirrKVg\nTAzhMYV3382L2cffo4fTNVJY2OipZoWTSC4mhVjvo0ePEPPmreP33z0cfHAp96W+WCZB1lIwJob8\nfOfMd+5cH/vs03jqhvAdx0VF8dfoc+asZ+pUZ/rU6KSQK11JTSW3vn1D9O8f5Oqrq+ttf+YZHw8/\nbJewZhpLCsbEUFIS4osv8igrC7Hxxo1ru3AyKCqK/5g77RSkZ8/GLYVcmTG14dxQDV14YU2D34u5\n4grnAwyFbLrxTJEjX0dj2ldxMaxf72HHHWNP8BYeYE6kpdCUXGkphBNeov71rzx69ixj883LIpcB\nm/SxpGBMDCUlTgW3ySaxK7rwWEIiLYWmJJIUGk7fnUl22inuRRMBuP56p2lw7rl1s/9ddpmXffct\nYd06GDeu0FoPaWBJwZgYwmMGTZ39hm9CS3Tuo7Do/vd4k8Ltt7e9hrzssuqWd2olny+xhDVqVC0/\n/riWNWvqPoB77vHy5Zd5PPlkPvfeW8CTT+azalWONKWyhF19ZEwM4ZZCr16xz37buoJaa5LCZpsF\n29zV1L9/8vpnWlM2n8+Z+qO21sOFF1az224+lizxc+21ThPsyiuLWLy4hokTq6mtdfbPle62TGVJ\nwZgYWmophEUvlpOIYFSuibeSa4/KMJkVamsvrV24cB2Fhc5n3rlzHuXlNdx4YyEnn1xLVRVMm5bP\nAQf4+dOfSnj00UoeeSSfKVOqWj2GYZpn3UfGxBBuATS13GZYw5lT47XvvonPb5SrZ8gbbkijRYU+\n/ngdf/97FWPH1hAKefjTn0oAGDGimH//28f77+flzP0dmcaSgjFNePjhSv7wh6YHT5ctW8uWW7au\nZhowoO64e+wRX5eOx9M4MRx6aN3cSiNH1tCnT2KDvZlq001DFBRAv35Bfv55LZdeWs28eXV3R48c\nWcwttxRYYojTqlUe5szJq3fXeVMsKRjThCOP9EcGlGMpKWmfOA3Xh25KrJbCtGlVLe4TLS8L56Xz\neGDs2Br69g1x5pk19OzpJL677y6kf/9SLrigiBkzfFQnbww9q02YUMDgwaVMnFjInnuWMmlS89ON\nWFIwJktED2oPH964RRDPWfPWW9d/3WGH1TaxZ2YaP76a995bzxtvrOfnn9fywgsVDBgQ4NFH89ll\nl1KuvLKQzz+3ai1M1cu0afl8+OF63nijgnvvreKnn5o/c7CBZmOywDPPVLD77nUT85WVNd4nVlLo\n1i3I1lsH+fBD5089nnGJvn2DLF2auRXrBhvALrs4yW2rrUJstVUtI0bUsny5h2eeyef004spLQ0x\naFCAnXcOsssuAbbdNhjXKnnZKhiEl17ykZ/vdLntsovzfbjqqkLGjKmJzOq7++4Bdt89ADT9YVhS\nMCZNNt00GJluuikiAVTz2G8/Z9zB4wkBsWv2eFoK4UttmzN79nq23jpG1mlBuvv3t9gixGWX1XDJ\nJTV88omXBQvymD8/j4ceymfFCi8nnljLmWc63VC5JBCAM88s4scfvfTqFWTJkjxWrvTQt28JoRCc\ndVZirUFLCsakSTyV6IgRtVx+eXwDAdED0T/+uJbevcsaxYleGrQpnTvHFa7NuveIHag91lw42P2p\nZ6r7086xmhMs7UTFJVdQOWp00mK89JKPn3/2MnNmRaQ15PcXM316LfvuG2h2XCyWzG0jGpPj4k0K\nP/+8NvJ7rG6jsOgxh3grgosuqn85SllZcs+ig6VZuqxcK3nXr6Pk7+OTGmPSpAIuu6y6XvdYly5w\nwgn+Vt3LkXFJQUT2FJHH3J84zmuMyU7xJIWGl6HOmFHB/Pnrmn5BgnHC92G0dD9G2KhRcVzT2IyK\nS67okIkhWVS9rFnjYciQ9rtTPRO7j84GzgF2B06kUYPPmNzQmj74Xr2aflFrbm7bdNMQs2atZ/PN\nQ+ywQ6cWj3H99dVMmdL6EdvKUaOb7Urp3LmY8vLKVh8/EQ1jrVjh4Z138njrLR9vveWjW7cg++0X\n4NBD/ey2WyDhBZWa6h6Lxe93WnqJTpsyc6aPww/3t+v065mYFPJUtUZEVgD7pbswxiRLew3M/vTT\nWnr1KmO77QK88UbjP+lQqPmafrfdgvz2W9vL0SnLGwAbbxzipJP8nHSSn1AI5s/38uabPi64oIia\nGjj55Fp22ilIMAi9ewfZfvtgo7W7mxJOEO09hjEh/ODBGDGbe2EzX76UJgUR2R2YoKrDRMQLTAF2\nAqqBs1T1a6BCRAqA3sBPqSyfMal0zjm1/PBD6+auiP6b9nph5Upn3OHuuxufzob3/fLL9u3GWLBg\nHf37d+LYY2u56abqyGWPucDjgYEDgwwcWMNll9WwaJGX55/38eKLPvLyYPlyL4sXe9lssyDHHuvn\n9NNrG015EiztlNSuo2RJWVIQkUuBPwHhT+looEBV93STxe3utqnAA27ZRqaqfMak2ujRbeufj6VL\nlxBff11/WzgpNDdPU7jbaJtt4p8mI7zWxEUX1eRUQohl222DXH11/f+v2lpYuNDLQw8VMGBAKUVF\nITbfPMRGG4UIheC47tcyovIGSoPZlRhS2VL4CjgWeML9fW/gNQBVnSciu7mPPwZGpLBcxmSdvfcO\nMWNG4wr8yScrqKxsufWx884BPv208aWuzz9f0eJrJ02q5MILi1vcL9fl5zs30d1zTxWTJ8PKlR6+\n/97Db7958HohP/98ZnvPp3//ACUlbR8v+f13eOstH5MmFeDzOfEeeqiSgQMbfw9aitVc11LKkoKq\nTheRPlGbyoDyqN8DIuJV1dyY0cuYJDr11BBHHtm4At9oo/ivJIK6FkL43/B8ToccUsvy5c6CNwC7\n7RZg/nzn8XHH+dlmm/WtL3wO8nicadaTOZ33hhvCMcf4OeooP7Nn51Fd7YmZENoqnQPN5TiJISzh\nhODzeencOXVnLKmMZ7GyK1aq48Ufy0NenlPjR++fl+eNbPN667qYwvu89BJMnAjXXutsHzcODj+8\n7vmhQ+sidOpU2G43vGXmZ5h58YYPDz+KvfRfW2KlMynMBY4EnhORwcBniR7A7w+m7PI1SO/lchYr\ns2OlOl58scoIhUL4/SEgL2r/MoLBYGSb14s7w2hZvWNWVxcAeXTqFGLgwApWroTy8voRjjuuiC5d\nqhptT+77yr5YqY7XYvdR96bvgkzHzWvh9tWLQJWIzMUZZL44DWUxJufdcksV995bv4IIdxeFWwhF\nRXVXMIX95S8tD4Tfd19Vu00hbjJDSlsKqroM2NN9HALOS2V8YzqiXXcNsuuuifc9l5XBAw9U0qVL\nDk8vahrJxJvXjDEZ4phj/HTunN9u3UMm81lSMCZHHXywn403brqFsGjR2qxcic0klyUFY3LUE080\nP6jZpUuKCmKySsbNkmqMMSZ9LCkY0wG1ZkZV0zFYUjDGGBNhScGYDiiRqTBMx2IDzcZ0MAsWrEv6\nspsme1lSMKaDCU95bUws1n1kjDEmwpKCMcaYCEsKxhhjIiwpGGOMibCkYIwxJsKSgjHGmAhLCsYY\nYyIsKRhjjImwpGCMMSbCkoIxxpgISwrGGGMiLCkYY4yJsKRgjDEmwpKCMcaYCEsKxhhjIiwpGGOM\nibCkYIwxJiIjk4KI7CciD6a7HMYY09FkXFIQka2AXYCidJfFGGM6moxLCqr6tareke5yGGNMR+RL\nRRAR2R2YoKrDRMQLTAF2AqqBs1T1axG5EdgaOE9Vf09FuYwxxtSX9KQgIpcCfwLWuZuOBgpUdU83\nWdwOHK2q1yS7LMYYY5qXiu6jr4BjAY/7+97AawCqOg/YLdaLVPW0FJTNGGNMlKQnBVWdDvijNpUB\n5VG/B9wuJWOMMWmWkjGFBspxEkOYV1WDrTlQYaHP0717Wcs7tqNUxrNY2RUr1fEsVnbFSnW81sZK\nxxn6XOAwABEZDHyWhjIYY4yJIZUthZD774vAgSIy1/19RArLYIwxphmeUCjU8l7GGGM6BBvgNcYY\nE2FJwRhjTIQlBWOMMRHpuCQ1aURkP+BkVT071u/JiiUiewLnuE+NUdU1SYh3AnAQUANcpaq/tXeM\nqFiHA8cB+cDtqvpJsmK58cbgTILYD3hSVe9PYqztgDFAAXCbqi5MYqydgcnA18A0VX07WbHceD2B\nmao6MMlxBgAX4NyQeqmqrkxyvP2BE4ESYKKqJvWKxWTWG1Exkl5nxIgZ1/vKmZZCw9lVkznbatSx\nC91NZ+P8Bz+M8+VNhqOAkcBDbrxkWgX0BjYBvktyLFT1bpzPb2EyE4LrLOB7oApYluRYg4AVODdv\nJi35RLmE5L8ncL73FwGvAnukIF6xqp4D3IZzYpQ0KZylORV1RkQi7ytnkkLD2VWTOdtq1LHDU3fk\nqWoNTgWwcTJi4pxxPggcCXRLUoyws4ETgFuBw5McK+wU4IUUxNkK57N8HvhzkmO9h5OEJgJjkxlI\nRM4DnsRJdkmlqu8D2+G8p6S2It14M0WkFLgQeCzJsVI1S3Mq6oyIRN5XRncfpXJ21TbGqhCRApyz\n65/aOeYNON0qz+JUMEOAHZP0/sKxSoD1OC2G7RKNlUC88Oc5CthHVc9KYqzwe/sFqAB+oxUnRQnG\negXnj/53WvG3lmCsHu5zg0TkOFVNKMEmGOsOYD5wKHAdTndcMt5b+PsxBpgAXKuqq5Icq82zNMcT\nj1bWGW2IF7eMTQqpnF21HWJNBR7A+TxHtnPMa9399wUexekLPzdJ7y8cazBOszaE0yWRkNZ8niJS\nkmicVr63ATgtLg8JVmatiLUHTqukFhiXzFhRr3u8FQkh0fc1DHgEZ3zrgURiJRjvGnf/aTit4/Ei\n8lIi7y/VszTHG49W1BltjBe3jE0K1M2u+oT7e73ZVUUkrtlV45xttU2xVPVjEr8zO6GYqvoO8E6C\nMVob60Pgw1bGSjieu/2UVMRS1Y+A01MU6wPgg1TEiorZmi6xRN/XW8BbrYjT2nit/f9KOFZUzNbO\n0hxXvFbWGa2OFxbP+8rYMYVUzq6ajplcc/n95ep7s1jtw95bZsfL2KQQQ7vNrpphsdIRM9XvL1ff\nm8XKvni5/N7aJV42JYVUzq6ajplcc/n95ep7s1jZFy+X31u7xMvkMYWwVM6umo6ZXHP5/eXqe7NY\n2Rcvl99bu8azWVKNMcZEZFP3kTHGmCSzpGCMMSbCkoIxxpgISwrGGGMiLCkYY4yJsKRgjDEmwpKC\nMcaYCEsKxhhjIrLhjmZjEiYifYDFNF7xbKqq3pf6EjlE5Ayc9QhewVmLYKlbpnOj9tkF+BgYoarT\nmjjOmcBwVT2kwfZHgQU4d7JuB/RT1W+T8FZMjrKkYHLZD6raP92FaCAEvKSqf3ET12rgYBGJnrjs\nRJzFgJqbbuCfwO0i0l1Vf4HImhSHA39V1UkisjRp78LkLEsKpkMSkRXAczjzz/uBE1R1mYgMxDmT\nL8FZeW6ku/1tnAp8e5xKW3AWzqnAOav34cxpf6Oq7uXGOB3YXVVHRYX2ULeMKziLoyzAWVHvbXfb\ngcCc8H4icogbKx+nZXG2qv4qIi+6ZbnHfd3RwL9U9be2fj6m47IxBZPLeovIggY/27vP9QTmqOqu\nwLvABSKSDzwEnKyqA3CSw4Pu/iHgU1XdFvgRuBPYD9gN6AKEVPVNoJeI9HVf82ec1fJa8ixwPICb\nlAstXl8AAAIcSURBVD7DWdUMEekOjAcOcsv6Bs7a2bjHjl6c6M84K6IZ02rWUjC57McWuo9ec//9\nAudMfRtgS2CGiIT3iZ6bfp777z7A+6q6AiLLRR7jPjcNOE1EHgN6qup/4yjnTOBmEfHgnPn/EzjJ\nfW53YHPgbbdMeTgtFoB/A93cbqgqnPGD2XHEM6ZJlhRMh6WqNe7DEE5XTR7wTTiRuCtW9Yp6SaX7\nr9/dNyy6O+gxnGRThZMg4inHOhH5FCfZDAMuoy4peIH3VPUot0xFuIlKVUNuQjrFjfdEw2Mbkyjr\nPjKmrlJfBHQRkb3d3/8C/CPG/h8AA0Wkl3t2fxLuoLB7pc/3wHkkVkk/C0wA/quqAXdbCKd1soeI\n9HO3XQ1MjHrdNOA4nO6neLqqjGmWtRRMLustIgsabHtHVS+i/pU9IZwxgRoRGQ7c7Z6RrwEaLRqv\nqr+IyIXAbJwz9GXUtSLA6f45RlV/ilGmhlcUhX+fCTwMXNUg1s8i8hfgWRHJA74D/hT1/PcishJn\n2cXlMeIZkxBbZMeYBIlIF+BCYJzbhXM3sFhV7xWR8FVI/1TVl2K89gxgX1VN5mp+4VhL3Vh2n4KJ\nm3UfGZMgVf0V2BD4wh0LKAMedLuSfgD8sRKCKwT80R2ITgoRKRKRT4CNkxXD5C5rKRhjjImwloIx\nxpgISwrGGGMiLCkYY4yJsKRgjDEmwpKCMcaYCEsKxhhjIv4ffmxiGcOR8bsAAAAASUVORK5CYII=\n", - "text/plain": [ - "" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ "# Create a loglog plot of the U-235 continuous-energy fission cross section \n", "plt.loglog(u235.energy, fission.sigma, color='b', linewidth=1)\n", @@ -1877,22 +1010,11 @@ }, { "cell_type": "code", - "execution_count": 33, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "name": "stderr", - "output_type": "stream", - "text": [ - "/usr/local/lib/python2.7/dist-packages/numpy/lib/shape_base.py:872: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n", - " return c.reshape(shape_out)\n", - "/usr/local/lib/python2.7/dist-packages/numpy/lib/shape_base.py:872: DeprecationWarning: using a non-integer number instead of an integer will result in an error in the future\n", - " return c.reshape(shape_out)\n" - ] - } - ], + "outputs": [], "source": [ "# Construct a Pandas DataFrame for the microscopic nu-scattering matrix\n", "nuscatter = xs_library[moderator_cell.id]['nu-scatter']\n", @@ -1920,22 +1042,11 @@ }, { "cell_type": "code", - "execution_count": 34, + "execution_count": null, "metadata": { "collapsed": false }, - "outputs": [ - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAADSCAYAAABAbduaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGuRJREFUeJzt3Xm4HFWZx/HvDRN2ggqIIrvCqyACArIOCYJsigIyMgwC\nCaBBHh1mQLYgGLYBAWdcmIiCQEAGVGRHEYUQFlkDDMjyYxNG1gFEGFZJ0vPHOZ00l7t0bld1btX9\nfZ4nT7qqq0+d6vvW26dOVZ3qaTQamJlZ/Yya3xUwM7NyOMGbmdWUE7yZWU05wZuZ1ZQTvJlZTTnB\nm5nV1N/N7wqUKSJmA0tL+kvLvPHAFyXtMMDnVgBuAdZq/WyvZQ4BdgN6gAWAq4BJkt4eYl2PAu6W\ndFlEbADsLelr81jGROA9kr4zlDr0Kmtl4DHgBklje713FrAXvb7bPsrodzsiYj3gMEn/0GldR4qI\n2A/YDxgNNIA7gSMk/XmQz20NfEfSui3z1gJ+CIwBZgETJd3Zx2c3Av4NWIrUIPwz8E1J9w9xG94R\nExFxNfCPA8VRH2UsB/xS0qZDqUMf5V0HbA58WNKfWuaPBaYBB0v67iBl9LsdEXElcJCkB4uo77wY\niS34AS/8j4g9gRuADwywzD8AOwIbSVoHWB/4KDC5g3p9mrTjAqwJLD+vBUj6cRHJvcWbwGoRsWJz\nRkQsBmzGIN9j1u92SJrh5N6+iDgF2An4rKQ1gbWA3wE3R8SH+vnMIhFxHPBzUiOkOX9R4GrgREmf\nBI4Fzu/j8wsBVwAHSlpb0lrAecBvIqJniJvSOya2IjWS2ibp6aKSe9YAngC+3Gv+XsBztBfr/W6H\npM/Oj+QONW/BZ72/9H6DKbcMvgBsB9w3QJkfIO0wiwJvSXorIr4OLJPLWZzUOtoEmAlcIumIiFgd\n+E9gMWA54G5gV2BfYD3gpIhYBDgaWDIifippn4jYATgCWBB4ndSCuiUiJgMb5/rcAzxCalV/IyIe\nB84CtgRWBH4u6dBcv8OAvYH/I/2YfUHSKn1s5yxSctgdOCHP2xm4BDgolzUK+A9gQ2CJ/P3uC/wP\ncAwwJiJ+CpwD/AB4NW//IcB3gU8AvwfukHRoRGyV6/1JSc8P8DcYMSJieWAisLyklwEkNYBz85HQ\n4cDX+/jo1sAipL/1Mb3mPyzpqjx9OfAn3m1RYEnS35W83vMi4mVS7ng7IvYGDiTFygukpPg0g8fE\nmS3ruTYits+vf0iK19HABZJOyEeTNwD3Ayvndfxe0uJ5H1iZtA+sBDwP7CrpmYj4FDAll/Vofv9A\nSdP72NbzSHF+LMz5EdyUFJs9ed7nSN/1gsD7gamSjspHtM3t+CxwI6kH4BPAJOB7wBeBjwNH5fkA\ndwDHS/pZH/UpxEhowU+LiLua/0jJs89f5Nwy+KKkBwYpcyrwV+DZiPhDbl2tKOmO/P4xpCD4KLAO\nsGk+3NsXOEvSJsBHgFWA7SX9J+mPfXD+Yx9F6hrZJyJWA44HtsutrYnARTkAAVYA1pW0R55utPy/\nmKTNST8034iIlSJiG9IOsr6k9YDF+/s+snN5Z8tmT+DslukNgQ9I2ii3LM8hdb08CRzZ3A7STrIm\n6TB2HeBvMCdR7Q7sGRFfAM4EdnNyf4cNgQeayb2Xa0hHVO8i6VJJBwEv9XprdeC5iDgjIm4ntebf\n1diT9BLph/iqiHg0Is6JiAnANZLejoi1gROBbSStDVxGaoh8isFjYm9Je+dVbSHpKVKsnSlp/bzN\nn8lHywAfAo6RFMCzvDNmNwN2kfSxvK0TI2IB4FekLqy1SY2Ldeg/1u8C/pZ/FCA1ZC4jNdCanzkQ\n2FPSBqSG1eER8T5JE1q248m8/L2S1pB0SZ5uSJoK3AyclOszvczkDiOjBT+uVx/8XsAunRQo6RVg\nm4hYBdgCGAdcGRFTJB1GajX/a05eb+f3iYjrga0j4mAgSK34xftYRetRxmeAD5JaB815s0g/EA3g\nFkmz+6nqpbm+T0fE/5L6UbcHfpG3AdIRxZYDbOudETE7Ij5Jah0tIem+Zl0k3RwRL0bE14BV87Y2\ny+59tPTnvvqLJT0bEV8h7VBHSrqxv/qMYKP7mb8Q7XUh9C5re9K+cXtEfB74dUSs2PsckqT/iIif\nkP6umwOHAofmRLglcFVOzkj6fvNzEXFkmzHRXH4xYCzw3og4Ns9eDFgbuJ2UaG/uZ3umSXo1v74L\neB+pC6sh6be5btdFxB8H+V7OITVmbiM1ZP4V+GZLnXcAdoiI3YGP5fmLAX2dP7ihn3XsRzrafh34\n5CD16dhISPC9zQmwiLibuTvHPn2dZOpLRBwKXC/pZtKh7ZkRsSnpROthpGBsXf5DpP7sKaSunZ8D\nV5Ja34P1P44itZj+saW8FYEnSX2yrw3w2TdaXjfyut7mnUdu/f04tGq24p8n7QRz5EPS7wGnkLpu\nHuTdfZlNr/YzH9Lh67Oklpu90y2kcyHLSnqu13tbADflrpoz8rxGPtrrz1PAg5JuB8gn9s8gJWM1\nF8oxvYmkk0nxemVETALuJTU83vFjkPvsVwJWo/2YaGqeI9hY0pu5vKVJMbwMqSu0v1h9s+V1M85n\n8u59a9YA62+QumlmRMS/A2NaGjKNfMR8N+mo4AbSkeaOfayjqb9Y/wDpR3k06aikr66xwoyELpp+\nSVpH0rr5X1vJPVsYODEHYNMawIz8+vfAXhHRk4P+V6TWz9akw8xf5uU2ZG5gzyR160DacZottmmk\nVn8ARMS2pEBbmL7PLwz0g9Eg7ahfjIgxed4+DN4C/BnwJdL5gv/qtb6tgMsl/Zi0/Tv12qb+Wp5z\nRMSGwD+TzkO8JyL+ebDPjCS5hfwD4Px8ngiA3F2yE+kKmRktsTxYy/A3wMr5qIyI2Jz0Q9872TwP\nHJHfb/oQqdV6Dyk2t4qI5gUJXwNOZt5iYhawYD6ivIW553aWJCXSzw+yLb014/8B4K3cJUk+4liL\nAWJd0jN5u87knQ2ZHtKP1hKkI8wrSUclC7Vs1yzm7r99iojRpJPZR5K6cc+PiFIb2XVP8H39MRv9\nzG/ns03HkpL4jRHxQESIlMC/lN8/mtTH/N+kS9mulHQx6YTLxRHxB1I/+69IXS2QTnSdEhF7kA5F\nPxoRv5J0H/BV4IJ8xHEssIOk1/vYlkG3TdI04HTS1Re3ky6Te32g70DS06QTXA9J+muvdZ0GjM3n\nN35NurJj5bzMH5rb0U/dGvmE9HnA1/MONh44KvfvWiZpEumH9tKIuDciHiJdebVxX91efZjz3eej\ngB2BKRFxL+lk986S/tZrnQ/l5Y6NiD9FxH3ABcBXJD0s6Y/AwaQ++rtJDZiJtB8TABeR9qM1gH8C\nNoqIe4BbgfMlNa/ueVfstPz/rn1A0izSic3JEXEnqf/8WfqP9aZzSP3rrQ2ZBinxXwE8EBE3kI44\n72Du/nsRcENErNlHmc2jiuOBpyWdKel04EXguEHq05EeDxc8suRD+U0k/TBPHwhsIGm3+Vszs2JF\nxEnAKZL+N9K9LXcDq7Scf6q9kdgHP9I9RDpJ9lXmXv/71flbJbNSPAFcExFvk1rQ+4yk5A5uwZuZ\n1Vbd++DNzEYsJ3gzs5oaNn3wPTPm+WaNfq263kCjDMybx6b3dVLcqqgxdt7GPCnK3/O7wmL7xp4+\nh50Zol8UWJbNT43G5D5j2y14M7OacoI3M6spJ3gzs5pygjczq6lST7JGGit8Cmn847eAfSU9WuY6\nzcrmuLaqKLsFvyNpIKFNSKMsDvjYK7OKcFxbJZSd4JtD6CLpVtKj7cyqznFtlVB2gh/D3IH+AWbl\nw1uzKnNcWyWUHZSv0PI8R2DUAIP2m1WF49oqoewEfxPp0WBExEakMZXNqs5xbZVQ9lAFF5MenHtT\nnp4w0MJmFeG4tkooNcHnh05/rcx1mHWb49qqwieGzMxqygnezKymnODNzGrKCd7MrKac4M3MamrY\nPNGJV4sralFeL6ysjcdeW1hZN0//dGFlWXXc2HPT4Au16dscXVhZRxc6hM4rgy9iXecWvJlZTTnB\nm5nVlBO8mVlNOcGbmdWUE7yZWU11LcFHxIYRMa1b6zPrBse1DWdduUwyIg4BvkyhF0OazV+Oaxvu\nutWCfwTYGejp0vrMusFxbcNaVxK8pIuAmd1Yl1m3OK5tuPNJVjOzmnKCNzOrqW4n+EaX12fWDY5r\nG5a6NtiYpMeBTbq1PrNucFzbcOYuGjOzmnKCNzOrKSd4M7OacoI3M6spJ3gzs5oaPo/sK9Afp29Q\nWFnHjv1mYWXNGrtAYWXdduvYwsoCir0f0/d2luZovl1YWcdxUGFlfcuP/xuW3II3M6spJ3gzs5py\ngjczqykneDOzmnKCNzOrqdKuoomI0cCZwErAQsBxki4va31m3eLYtqooswW/O/C8pM2BbYFTS1yX\nWTc5tq0SyrwO/pfAhfn1KHx1tNWHY9sqobQEL+k1gIhYgrRDHFHWusy6ybFtVVHqSdaIWAG4FjhH\n0gVlrsusmxzbVgVlnmRdFrga2F/StLLWY9Ztjm2rijL74CcBSwJHRcRRed52kt4scZ1m3eDYtkoo\nsw/+AOCAsso3m18c21YVvtHJzKymnODNzGrKCd7MrKac4M3Maqqtk6wRMYZ01UBPc56k/ymrUmbd\n4ti2Ohs0wUfEJOAw4C9Ao+WtVcqq1HBy5PRTiivse8UVdcnF2xRXGPBvTCqsrEdmfaSwsv7y5PsL\nKwtGv2NqpMf2twp8/N8tBT7+byPuKKysZOSOA9dOC35f4MOSni+7MmZd5ti2WmunD/4J4KWyK2I2\nHzi2rdbaacE/AtwYEdcCb+V5DUnHlFcts65wbFuttZPgn8r/mnr6W9CsYhzbVmuDJnhJk7tQD7Ou\nc2xb3fWb4CNioFHyGpI+PVjhEbEAcDqwOukqhf0k3TfPtTQrUKex7bi2qhioBX/0AO81Bniv1eeA\n2ZI2i4ixwPHAju1Wzqwknca249oqod8EL+m6TguXdGlEXJEnV8ZXLNgw0GlsO66tKsocDx4ASbMi\n4mxgJ2CXstdn1g2Oa6uCroxFI2k8qb/y9IhYpBvrNCub49qGu3bHolkW2Ax4G7hBUluHpBGxB7C8\npBOAN4DZ+Z/ZsDCU2HZcW1UM2oKPiC8D/w38EzABuC8iPttm+RcC60TEdOAq4ABJbw3yGbOu6CC2\nHddWCe204I8E1pP0FEBErARcAVw52AclvQHs2lENzcozpNh2XFtVtNMH/wrwTHNC0hPMva3brMoc\n21Zr7bTg7wIui4jTgVnAbsBTEfElAEm/KLF+ZmVybFuttZPgFwReYO6NHG+Txs/eLk97J7Cqcmxb\nrbUzFs34LtTDrOsc21Z37TzR6U99zG5IWrWE+ph1jWPb6q6dLpotWl6PJh3OLlxOdcy6yrFttdbT\naLQ7bthcETFD0nqFVmR62wOYGcC4yYUWN/vFgcbfmjcnvu9fCivrtxT37Nnr2HbQ8d5Lie2eySMg\ntou7kbex9WGFlQXQ81iBX/8jk4srq0CNxuQ+Y7udLpqxzB1hrwf4OG7lWA04tq3u2umiOZq5O0GD\ndNXBXqXVyKx7HNtWa+1cRTMOICLGAAu0Ow6N2XDn2La6a6eL5sPA+cBHgJ6IeBzYVdJD5VbNrFyO\nbau7doYq+DFwkqT3SXovcALwk3ZXEBHvj4g/R8TqQ62kWUkc21Zr7ST4pSVd2JzIt28v1U7hETGa\ntBO9NrTqmZXKsW211k6Cfysi5lw2FhHr035Qnwz8iJYBncyGEce21Vo7V9EcAFwYEc0TUEvRxlCp\nETEeeF7S1RFxOOkyNLPhxLFttdZOgl8KCNKjyUYBavPhBhOARkRsBawDTI2IL0h6bsi1NSuWY9tq\nrZ0Ef7KkNYA/zkvBksY2X0fENGCidwAbZhzbVmvtJPhHI+JM4FbgzTyvIemc8qpl1hWObau1dhL8\ni6Q+xo3ydA/prr+2dwJJWwy+lFnXObat1jwevI1Yjm2ruwETfETsDzwj6eKIuA1YBpgJbCfpkW5U\n0KwMjm0bCfq9Dj5f/rUzcH+etTAwDvg+cHjpNTMriWPbRoqBbnTaC9hJkvL0rPzU+R8BG5deM7Py\nOLZtRBgowc+S9H8t08cDSJoFtHOtsNlw5di2EWGgBN+Th1EFoDlmR0QsCX76klWaY9tGhIFOsp5H\nukNvvKSXASJiCeCs/J7NT3dMLrS4UUuNGXyhNjVOO6iwsmKiBl+obds2Xzi2C/dGYSX1XN32gJ5t\naXynuJEkeh4u8Pf/jMnFldWPgRL8d4ApwNMRcT+pZbMGcC7w76XXzKw8jm0bEfpN8JJmAl+NiGOA\nT+XZM/LJKLPKcmzbSNHOjU5PAk92oS5mXeXYtrprZzx4MzOroHbGohmyiLgTeDlPPiZpnzLXZ9Yt\njm2rgtISfEQsDB6MyerHsW1VUWYLfm1g0Yj4bV7PJEm3lrg+s25xbFsllNkH/xrpgQrbAPsB50WE\n+/ytDhzbVgllBuVD5JtGJD1MGnv7gyWuz6xbHNtWCWUm+AnAdwEiYjlgDH4CvdWDY9sqocw++J8C\nZ0XE9Xl6gqTZJa7PrFsc21YJpSX4fLfgHmWVbza/OLatKnxiyMysppzgzcxqygnezKymnODNzGrK\nCd7MrKac4M3MaqrU0SStRK8WW9zCfx1fWFk93zuwsLIaHyvucWs8UFxRVqanCi2t59ALCyurManA\nx/9tW/7jf92CNzOrKSd4M7OacoI3M6spJ3gzs5oq+5F9hwM7AKOBUyVNLXN9Zt3i2LYqKK0FHxHj\ngI0lbQKMA1Yta11m3eTYtqooswW/NXBvRFxCGi/74BLXZdZNjm2rhDIT/DLACsDnSC2cy4CPlrg+\ns25xbFsllHmS9QXgakkzJT0EvBkRS5e4PrNucWxbJZSZ4G8EtoU5jzVbjPTsSrOqc2xbJZSW4CVd\nCdwVEbeRDmH3l1T+vblmJXNsW1WUepmkpEPLLN9sfnFsWxX4Riczs5pygjczqykneDOzmnKCNzOr\nKSd4M7OacoI3M6upnkZjeFy+2zOd4VGRkWrh4oq6ZsNNCivr+p6bCytrcqNR4PP/2tfTM9mxXRtH\nFFbSr1mwsLK26ye23YI3M6spJ3gzs5pygjczqykneDOzmir7kX17AePz5CLA2sCykl4pc71mZXJc\nW1WUPdjYVGAqQEScCpzhncCqznFtVdGVLpqIWB9YU9IZ3VifWTc4rm2461Yf/CRgcpfWZdYtjmsb\n1kpP8BHxHmB1SdPLXpdZtziurQq60YLfHLimC+sx6ybHtQ173UjwqwOPdmE9Zt3kuLZhr9SraAAk\nnVL2Osy6zXFtVeAbnczMasoJ3sysppzgzcxqygnezKymnODNzGrKCd7MrKaGzSP7zMysWG7Bm5nV\nlBO8mVlNOcGbmdWUE7yZWU05wZuZ1ZQTvJlZTZU+mmQRImIUMAX4BPAWsK+kjoZqjYgNgRMlbdFB\nGaOBM4GVgIWA4yRdPsSyFgBOJw1D2wD2k3TfUOuWy3w/MAPYUtJDHZRzJ/BynnxM0j4dlHU4sAMw\nGjg1P990qGVV/uHXju0h16+2sV1kXFelBb8jsKCkTYDDgO92UlhEHEIKuIU6rNfuwPOSNge2BU7t\noKzPAbMlbQZ8Czi+k4rlHfTHwGsdlrMwgKQt8r9OdoBxwMb57zgOWLWTukma2qwXcAfwjSol98yx\nPY/qHttFxnVVEvymwFUAkm4F1u+wvEeAnYGeDsv5JXBUfj0KmDnUgiRdCkzMkysDL3VUMzgZ+BHw\nTIflrA0sGhG/jYhrcutwqLYG7o2IS4DLgcs6rBtQ+YdfO7bn3YiI7SLiuioJfgzQ+gs2Kx/aDomk\ni+ggYFvKeU3SqxGxBGmHOKLD8mZFxNnAD4D/Gmo5ETGe1Pq6Os/qZGd/DThZ0jbAfsB5HXz3ywDr\nAbs0y+qgXq2q/PBrx/Y8GGGx3XFcVyXBvwIs0TI9StLs+VWZVhGxAnAtcI6kCzotT9J4Ul/l6RGx\nyBCLmQB8JiKmAesAUyNi2SGW9RA5WCU9DLwIfHCIZb0AXC1pZu43fTMilh5iWUAtHn7t2J43IyK2\ni4rrqiT4m4DtASJiI+Ce+VudJAfW1cAhks7usKw98kkagDeA2fnfPJM0VtK43Id3N7CnpOeGWLUJ\n5H7hiFiO1OIc6qHxjaT+3GZZi5F2qk5U/eHXju15MIJiu5C4rsRVNMDFpF/tm/L0hILK7XSktUnA\nksBREdHsr9xO0ptDKOtC4OyImE46C3+ApLc6rF8RfgqcFRHX5+kJQ21hSroyIjaPiNtIjYv9JXX6\nN6j6w68d2/PPcI7tQuLao0mamdVUVbpozMxsHjnBm5nVlBO8mVlNOcGbmdWUE7yZWU05wZuZ1VRV\nroOvnIj4O+BQ0qBNDWABYKqkE7pcj9WAU4A1gDcBAd+U9Pggnzsa+J2kG0uvpFWG47pa3IIvzxTS\nwFEbSVoT2ADYMiL271YF8t2I1wIXSFpN0lqkG2tuauM26s1JO69ZK8d1hfhGpxJExPKkFsVykl5u\nmR/AGpIuzgMvLQV8GDiYdFvz90nDvL4ATJT0aERcB3xb0vSIWBmYJmmV/Pm/AeuSbrE+VtLPetVj\nMrCipL17zf85cK+k4yJitqRRef54YCxp55lCum17p07H7rZ6cFxXj1vw5fgUcH/rTgCg5OI82SCN\nircG8DvgfNLtzesAp+Xp5nL9/QovB2wIfBo4pY9Bl9YHbuvjc9eTWl69NYCGpHNJ41DvOxJ2Amub\n47pinODLMyd4I2KXiLgrIu7JY1U03Zr/Xx14SdIMAEkXAh+JiDGDlH+6pNmSniINWrVZH8uM7uOz\nC9H/ztXTz2szcFxXihN8OWYAa+SxtJF0oaR1SY/zWqZluebATX39HXpIfYUN5gZk76Ce1fJ6FPB2\nr/dvBTbuo+yNgdv7mL9gr2n331krx3XFOMGXQNITwLmksaqXhDnPpdyBuQ9jaG1FCFgqP8GFiPgS\n8Likl0j9lh/Py+3Y8pkeYLe8/EqkQ9obelVlCrBZROzenBERe5J2hNPyrBciYs2I6AE+z9zgn0nf\nrSQboRzX1eMEX579SYeX0yLiLuBe0omj7fL7c/og89CpuwKnRsS9+bO75uVOAvaPiBnAwswN1Aaw\neETcAVwBfCXvOHNI+gvw98COEfFgRIgU7JtJao5VfVj+/B+AB1s+fhVwWh6j3KzJcV0hvoqmoiLi\nLOA3kn4xv+tiVhTHdbHcgjczqym34M3MasoteDOzmnKCNzOrKSd4M7OacoI3M6spJ3gzs5pygjcz\nq6n/BzrS9bjl9QXqAAAAAElFTkSuQmCC\n", - "text/plain": [ - "" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ "# Create plot of the H-1 scattering matrix\n", "fig = plt.subplot(121)\n", @@ -1956,6 +1067,15 @@ "# Show the plot on screen\n", "plt.show()" ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], + "source": [] } ], "metadata": { @@ -1974,7 +1094,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", - "version": "2.7.10" + "version": "2.7.6" } }, "nbformat": 4, diff --git a/docs/source/pythonapi/examples/mgxs-part-iii.ipynb b/docs/source/pythonapi/examples/mgxs-part-iii.ipynb index 119fd0b56..a541efbf0 100644 --- a/docs/source/pythonapi/examples/mgxs-part-iii.ipynb +++ b/docs/source/pythonapi/examples/mgxs-part-iii.ipynb @@ -1452,7 +1452,6 @@ ], "source": [ "# Generate tracks for OpenMOC\n", - "openmoc_geometry.initializeFlatSourceRegions()\n", "track_generator = openmoc.TrackGenerator(openmoc_geometry, num_azim=32, spacing=0.1)\n", "track_generator.generateTracks()\n", "\n", @@ -1638,7 +1637,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", - "version": "2.7.10" + "version": "2.7.6" } }, "nbformat": 4, diff --git a/openmc/tallies.py b/openmc/tallies.py index 980c5fa2b..9b75f35d4 100644 --- a/openmc/tallies.py +++ b/openmc/tallies.py @@ -1003,16 +1003,17 @@ class Tally(object): A list of filter type strings (e.g., ['mesh', 'energy']; default is []) filter_bins : list of Iterables - A list of the filter bins corresponding to the filter_types - parameter (e.g., [(1,), (0., 0.625e-6)]; default is []). Each bin - in the list is the integer ID for 'material', 'surface', 'cell', - 'cellborn', and 'universe' Filters. Each bin is an integer for the - cell instance ID for 'distribcell' Filters. Each bin is a 2-tuple of - floats for 'energy' and 'energyout' filters corresponding to the - energy boundaries of the bin of interest. The bin is a (x,y,z) - 3-tuple for 'mesh' filters corresponding to the mesh cell of - interest. The order of the bins in the list must correspond to the - filter_types parameter. + A list of tuples of filter bins corresponding to the filter_types + parameter (e.g., [(1,), ((0., 0.625e-6),)]; default is []). Each + tuple contains bins for the corresponding filter type in the filters + parameter. Each bins is the integer ID for 'material', 'surface', + 'cell', 'cellborn', and 'universe' Filters. Each bin is an integer + for the cell instance ID for 'distribcell' Filters. Each bin is a + 2-tuple of floats for 'energy' and 'energyout' filters corresponding + to the energy boundaries of the bin of interest. The bin is an + (x,y,z) 3-tuple for 'mesh' filters corresponding to the mesh cell + of interest. The order of the bins in the list must correspond to + the filter_types parameter. Returns ------- @@ -1174,16 +1175,17 @@ class Tally(object): A list of filter type strings (e.g., ['mesh', 'energy']; default is []) filter_bins : list of Iterables - A list of the filter bins corresponding to the filter_types - parameter (e.g., [(1,), (0., 0.625e-6)]; default is []). Each bin - in the list is the integer ID for 'material', 'surface', 'cell', - 'cellborn', and 'universe' Filters. Each bin is an integer for the - cell instance ID for 'distribcell' Filters. Each bin is a 2-tuple of - floats for 'energy' and 'energyout' filters corresponding to the - energy boundaries of the bin of interest. The bin is a (x,y,z) - 3-tuple for 'mesh' filters corresponding to the mesh cell of - interest. The order of the bins in the list must correspond to the - filter_types parameter. + A list of tuples of filter bins corresponding to the filter_types + parameter (e.g., [(1,), ((0., 0.625e-6),)]; default is []). Each + tuple contains bins for the corresponding filter type in the filters + parameter. Each bins is the integer ID for 'material', 'surface', + 'cell', 'cellborn', and 'universe' Filters. Each bin is an integer + for the cell instance ID for 'distribcell' Filters. Each bin is a + 2-tuple of floats for 'energy' and 'energyout' filters corresponding + to the energy boundaries of the bin of interest. The bin is an + (x,y,z) 3-tuple for 'mesh' filters corresponding to the mesh cell + of interest. The order of the bins in the list must correspond to + the filter_types parameter. nuclides : list of str A list of nuclide name strings (e.g., ['U-235', 'U-238']; default is []) @@ -2641,16 +2643,17 @@ class Tally(object): A list of filter type strings (e.g., ['mesh', 'energy']; default is []) filter_bins : list of Iterables - A list of the filter bins corresponding to the filter_types - parameter (e.g., [(1,), (0., 0.625e-6)]; default is []). Each bin - in the list is the integer ID for 'material', 'surface', 'cell', - 'cellborn', and 'universe' Filters. Each bin is an integer for the - cell instance ID for 'distribcell' Filters. Each bin is a 2-tuple of - floats for 'energy' and 'energyout' filters corresponding to the - energy boundaries of the bin of interest. The bin is a (x,y,z) - 3-tuple for 'mesh' filters corresponding to the mesh cell of - interest. The order of the bins in the list must correspond to the - filter_types parameter. + A list of tuples of filter bins corresponding to the filter_types + parameter (e.g., [(1,), ((0., 0.625e-6),)]; default is []). Each + tuple contains bins to slice for the corresponding filter type in + the filters parameter. Each bins is the integer ID for 'material', + 'surface', 'cell', 'cellborn', and 'universe' Filters. Each bin is + an integer for the cell instance ID for 'distribcell' Filters. Each + bin is a 2-tuple of floats for 'energy' and 'energyout' filters + corresponding to the energy boundaries of the bin of interest. The + bin is an (x,y,z) 3-tuple for 'mesh' filters corresponding to the + mesh cell of interest. The order of the bins in the list must + correspond to the filter_types parameter. nuclides : list of str A list of nuclide name strings (e.g., ['U-235', 'U-238']; default is []) @@ -2736,6 +2739,7 @@ class Tally(object): for filter_bin in filter_bins[i]: bin_index = find_filter.get_bin_index(filter_bin) if filter_type in ['energy', 'energyout']: + bin_indices.extend([bin_index]) bin_indices.extend([bin_index, bin_index+1]) num_bins += 1 elif filter_type == 'distribcell': @@ -2745,7 +2749,7 @@ class Tally(object): bin_indices.append(bin_index) num_bins += 1 - find_filter.bins = find_filter.bins[bin_indices] + find_filter.bins = set(find_filter.bins[bin_indices]) find_filter.num_bins = num_bins # Update the new tally's filter strides