diff --git a/docs/source/pythonapi/examples/multi-group-cross-sections.ipynb b/docs/source/pythonapi/examples/multi-group-cross-sections.ipynb index e7820bca8..5e1533781 100644 --- a/docs/source/pythonapi/examples/multi-group-cross-sections.ipynb +++ b/docs/source/pythonapi/examples/multi-group-cross-sections.ipynb @@ -150,9 +150,20 @@ "cell_type": "code", "execution_count": 6, "metadata": { - "collapsed": true + "collapsed": false }, - "outputs": [], + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/universe.py:223: DeprecationWarning: Cell.add_surface(...) has been deprecated and may be removed in a future version. The region for a Cell should be defined using the region property directly.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/universe.py:223: DeprecationWarning: Cell.add_surface(...) has been deprecated and may be removed in a future version. The region for a Cell should be defined using the region property directly.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/universe.py:223: DeprecationWarning: Cell.add_surface(...) has been deprecated and may be removed in a future version. The region for a Cell should be defined using the region property directly.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/universe.py:223: DeprecationWarning: Cell.add_surface(...) has been deprecated and may be removed in a future version. The region for a Cell should be defined using the region property directly.\n" + ] + } + ], "source": [ "# Instantiate a Cell\n", "cell = openmc.Cell(cell_id=1, name='cell')\n", @@ -451,8 +462,8 @@ " Copyright: 2011-2015 Massachusetts Institute of Technology\n", " License: http://mit-crpg.github.io/openmc/license.html\n", " Version: 0.7.0\n", - " Git SHA1: e0c2aace2e73367536fa03e153b67a2d038cd2b3\n", - " Date/Time: 2015-10-03 12:56:30\n", + " Git SHA1: 23535afa1c69644bb299bde18a094c3b99d53ae0\n", + " Date/Time: 2015-10-08 14:20:33\n", " MPI Processes: 1\n", "\n", " ===========================================================================\n", @@ -465,11 +476,11 @@ " Reading materials XML file...\n", " Reading tallies XML file...\n", " Building neighboring cells lists for each surface...\n", - " Loading ACE cross section table: 92238.71c\n", - " Loading ACE cross section table: 8016.71c\n", - " Loading ACE cross section table: 40090.71c\n", " Loading ACE cross section table: 1001.71c\n", + " Loading ACE cross section table: 8016.71c\n", " Loading ACE cross section table: 92235.71c\n", + " Loading ACE cross section table: 92238.71c\n", + " Loading ACE cross section table: 40090.71c\n", " Initializing source particles...\n", "\n", " ===========================================================================\n", @@ -478,56 +489,56 @@ "\n", " Bat./Gen. k Average k \n", " ========= ======== ==================== \n", - " 1/1 1.14249 \n", - " 2/1 1.18016 \n", - " 3/1 1.16083 \n", - " 4/1 1.09124 \n", - " 5/1 1.15214 \n", - " 6/1 1.13453 \n", - " 7/1 1.15552 \n", - " 8/1 1.18149 \n", - " 9/1 1.10404 \n", - " 10/1 1.15703 \n", - " 11/1 1.21224 \n", - " 12/1 1.14147 1.17686 +/- 0.03538\n", - " 13/1 1.12601 1.15991 +/- 0.02655\n", - " 14/1 1.11972 1.14986 +/- 0.02129\n", - " 15/1 1.15683 1.15125 +/- 0.01655\n", - " 16/1 1.15236 1.15144 +/- 0.01351\n", - " 17/1 1.17833 1.15528 +/- 0.01205\n", - " 18/1 1.13229 1.15241 +/- 0.01082\n", - " 19/1 1.22394 1.16035 +/- 0.01242\n", - " 20/1 1.15867 1.16019 +/- 0.01111\n", - " 21/1 1.13611 1.15800 +/- 0.01029\n", - " 22/1 1.14101 1.15658 +/- 0.00950\n", - " 23/1 1.20864 1.16059 +/- 0.00961\n", - " 24/1 1.13475 1.15874 +/- 0.00909\n", - " 25/1 1.10697 1.15529 +/- 0.00914\n", - " 26/1 1.20824 1.15860 +/- 0.00916\n", - " 27/1 1.16775 1.15914 +/- 0.00863\n", - " 28/1 1.15904 1.15913 +/- 0.00813\n", - " 29/1 1.16967 1.15969 +/- 0.00771\n", - " 30/1 1.12574 1.15799 +/- 0.00751\n", - " 31/1 1.16177 1.15817 +/- 0.00715\n", - " 32/1 1.18082 1.15920 +/- 0.00689\n", - " 33/1 1.19549 1.16078 +/- 0.00677\n", - " 34/1 1.18508 1.16179 +/- 0.00656\n", - " 35/1 1.17697 1.16240 +/- 0.00632\n", - " 36/1 1.16342 1.16244 +/- 0.00607\n", - " 37/1 1.17400 1.16286 +/- 0.00586\n", - " 38/1 1.19281 1.16393 +/- 0.00575\n", - " 39/1 1.15669 1.16368 +/- 0.00555\n", - " 40/1 1.17987 1.16422 +/- 0.00539\n", - " 41/1 1.14129 1.16348 +/- 0.00527\n", - " 42/1 1.18323 1.16410 +/- 0.00514\n", - " 43/1 1.13885 1.16334 +/- 0.00504\n", - " 44/1 1.17943 1.16381 +/- 0.00491\n", - " 45/1 1.20014 1.16485 +/- 0.00488\n", - " 46/1 1.16056 1.16473 +/- 0.00474\n", - " 47/1 1.20077 1.16570 +/- 0.00471\n", - " 48/1 1.15469 1.16541 +/- 0.00460\n", - " 49/1 1.18862 1.16601 +/- 0.00452\n", - " 50/1 1.18755 1.16655 +/- 0.00444\n", + " 1/1 1.19804 \n", + " 2/1 1.12945 \n", + " 3/1 1.15573 \n", + " 4/1 1.13929 \n", + " 5/1 1.16300 \n", + " 6/1 1.22117 \n", + " 7/1 1.19012 \n", + " 8/1 1.11299 \n", + " 9/1 1.16066 \n", + " 10/1 1.12566 \n", + " 11/1 1.20854 \n", + " 12/1 1.14691 1.17773 +/- 0.03082\n", + " 13/1 1.17204 1.17583 +/- 0.01789\n", + " 14/1 1.14148 1.16724 +/- 0.01529\n", + " 15/1 1.17272 1.16834 +/- 0.01189\n", + " 16/1 1.18575 1.17124 +/- 0.01014\n", + " 17/1 1.20498 1.17606 +/- 0.00983\n", + " 18/1 1.14754 1.17249 +/- 0.00923\n", + " 19/1 1.18141 1.17348 +/- 0.00820\n", + " 20/1 1.15074 1.17121 +/- 0.00768\n", + " 21/1 1.15914 1.17011 +/- 0.00703\n", + " 22/1 1.14586 1.16809 +/- 0.00673\n", + " 23/1 1.18999 1.16978 +/- 0.00642\n", + " 24/1 1.15101 1.16844 +/- 0.00609\n", + " 25/1 1.13791 1.16640 +/- 0.00602\n", + " 26/1 1.19791 1.16837 +/- 0.00597\n", + " 27/1 1.19818 1.17012 +/- 0.00587\n", + " 28/1 1.14160 1.16854 +/- 0.00576\n", + " 29/1 1.11487 1.16571 +/- 0.00614\n", + " 30/1 1.17538 1.16620 +/- 0.00584\n", + " 31/1 1.20210 1.16791 +/- 0.00581\n", + " 32/1 1.20078 1.16940 +/- 0.00574\n", + " 33/1 1.14624 1.16839 +/- 0.00558\n", + " 34/1 1.14618 1.16747 +/- 0.00542\n", + " 35/1 1.16866 1.16752 +/- 0.00520\n", + " 36/1 1.18565 1.16821 +/- 0.00504\n", + " 37/1 1.16824 1.16821 +/- 0.00485\n", + " 38/1 1.18299 1.16874 +/- 0.00471\n", + " 39/1 1.21418 1.17031 +/- 0.00480\n", + " 40/1 1.11167 1.16835 +/- 0.00504\n", + " 41/1 1.11545 1.16665 +/- 0.00516\n", + " 42/1 1.11114 1.16491 +/- 0.00529\n", + " 43/1 1.14227 1.16423 +/- 0.00517\n", + " 44/1 1.14104 1.16355 +/- 0.00506\n", + " 45/1 1.16756 1.16366 +/- 0.00492\n", + " 46/1 1.13065 1.16274 +/- 0.00487\n", + " 47/1 1.11251 1.16139 +/- 0.00492\n", + " 48/1 1.14731 1.16101 +/- 0.00481\n", + " 49/1 1.16691 1.16117 +/- 0.00469\n", + " 50/1 1.19679 1.16206 +/- 0.00465\n", " Creating state point statepoint.50.h5...\n", "\n", " ===========================================================================\n", @@ -537,27 +548,27 @@ "\n", " =======================> TIMING STATISTICS <=======================\n", "\n", - " Total time for initialization = 6.4800E-01 seconds\n", - " Reading cross sections = 1.5500E-01 seconds\n", - " Total time in simulation = 1.6951E+01 seconds\n", - " Time in transport only = 1.6927E+01 seconds\n", - " Time in inactive batches = 3.1560E+00 seconds\n", - " Time in active batches = 1.3795E+01 seconds\n", - " Time synchronizing fission bank = 7.0000E-03 seconds\n", - " Sampling source sites = 4.0000E-03 seconds\n", - " SEND/RECV source sites = 2.0000E-03 seconds\n", - " Time accumulating tallies = 0.0000E+00 seconds\n", - " Total time for finalization = 3.0000E-03 seconds\n", - " Total time elapsed = 1.7614E+01 seconds\n", - " Calculation Rate (inactive) = 7921.42 neutrons/second\n", - " Calculation Rate (active) = 7249.00 neutrons/second\n", + " Total time for initialization = 4.5400E-01 seconds\n", + " Reading cross sections = 1.0100E-01 seconds\n", + " Total time in simulation = 1.4106E+01 seconds\n", + " Time in transport only = 1.4092E+01 seconds\n", + " Time in inactive batches = 2.1000E+00 seconds\n", + " Time in active batches = 1.2006E+01 seconds\n", + " Time synchronizing fission bank = 4.0000E-03 seconds\n", + " Sampling source sites = 3.0000E-03 seconds\n", + " SEND/RECV source sites = 1.0000E-03 seconds\n", + " Time accumulating tallies = 3.0000E-03 seconds\n", + " Total time for finalization = 2.0000E-03 seconds\n", + " Total time elapsed = 1.4572E+01 seconds\n", + " Calculation Rate (inactive) = 11904.8 neutrons/second\n", + " Calculation Rate (active) = 8329.17 neutrons/second\n", "\n", " ============================> RESULTS <============================\n", "\n", - " k-effective (Collision) = 1.16600 +/- 0.00432\n", - " k-effective (Track-length) = 1.16655 +/- 0.00444\n", - " k-effective (Absorption) = 1.16281 +/- 0.00314\n", - " Combined k-effective = 1.16367 +/- 0.00307\n", + " k-effective (Collision) = 1.16131 +/- 0.00453\n", + " k-effective (Track-length) = 1.16206 +/- 0.00465\n", + " k-effective (Absorption) = 1.16096 +/- 0.00364\n", + " Combined k-effective = 1.16120 +/- 0.00325\n", " Leakage Fraction = 0.00000 +/- 0.00000\n", "\n" ] @@ -665,7 +676,7 @@ "name": "stderr", "output_type": "stream", "text": [ - "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/tallies.py:1485: RuntimeWarning: invalid value encountered in divide\n" + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/tallies.py:1486: RuntimeWarning: invalid value encountered in divide\n" ] } ], @@ -713,14 +724,14 @@ "\tDomain Type =\tcell\n", "\tDomain ID =\t1\n", "\tCross Sections [cm^-1]:\n", - " Group 1 [0.821 - 20.0 MeV]:\t1.11e-02 +/- 5.93e-01%\n", - " Group 2 [0.00553 - 0.821 MeV]:\t6.60e-04 +/- 3.04e-01%\n", - " Group 3 [4e-06 - 0.00553 MeV]:\t9.00e-03 +/- 4.10e-01%\n", - " Group 4 [6.25e-07 - 4e-06 MeV]:\t1.44e-02 +/- 6.58e-01%\n", - " Group 5 [2.8e-07 - 6.25e-07 MeV]:\t4.72e-02 +/- 9.80e-01%\n", - " Group 6 [1.4e-07 - 2.8e-07 MeV]:\t7.29e-02 +/- 8.59e-01%\n", - " Group 7 [5.8e-08 - 1.4e-07 MeV]:\t1.11e-01 +/- 7.92e-01%\n", - " Group 8 [0.0 - 5.8e-08 MeV]:\t2.39e-01 +/- 6.90e-01%\n", + " Group 1 [0.821 - 20.0 MeV]:\t1.11e-02 +/- 7.69e-01%\n", + " Group 2 [0.00553 - 0.821 MeV]:\t6.59e-04 +/- 2.97e-01%\n", + " Group 3 [4e-06 - 0.00553 MeV]:\t8.95e-03 +/- 5.12e-01%\n", + " Group 4 [6.25e-07 - 4e-06 MeV]:\t1.45e-02 +/- 7.10e-01%\n", + " Group 5 [2.8e-07 - 6.25e-07 MeV]:\t4.71e-02 +/- 1.02e+00%\n", + " Group 6 [1.4e-07 - 2.8e-07 MeV]:\t7.29e-02 +/- 8.86e-01%\n", + " Group 7 [5.8e-08 - 1.4e-07 MeV]:\t1.11e-01 +/- 6.67e-01%\n", + " Group 8 [0.0 - 5.8e-08 MeV]:\t2.38e-01 +/- 7.71e-01%\n", "\n", "\n", "\n" @@ -768,8 +779,8 @@ " 1\n", " 1\n", " total\n", - " 0.077460\n", - " 0.000890\n", + " 0.076970\n", + " 0.001012\n", " \n", " \n", " 62\n", @@ -777,8 +788,8 @@ " 1\n", " 2\n", " total\n", - " 0.087276\n", - " 0.000331\n", + " 0.087876\n", + " 0.000344\n", " \n", " \n", " 61\n", @@ -786,8 +797,8 @@ " 1\n", " 3\n", " total\n", - " 0.000450\n", - " 0.000026\n", + " 0.000418\n", + " 0.000023\n", " \n", " \n", " 60\n", @@ -849,8 +860,8 @@ " 2\n", " 2\n", " total\n", - " 0.266651\n", - " 0.001340\n", + " 0.266499\n", + " 0.001265\n", " \n", " \n", "\n", @@ -858,16 +869,16 @@ ], "text/plain": [ " cell group in group out nuclide mean std. dev.\n", - "63 1 1 1 total 0.077460 0.000890\n", - "62 1 1 2 total 0.087276 0.000331\n", - "61 1 1 3 total 0.000450 0.000026\n", + "63 1 1 1 total 0.076970 0.001012\n", + "62 1 1 2 total 0.087876 0.000344\n", + "61 1 1 3 total 0.000418 0.000023\n", "60 1 1 4 total 0.000000 0.000000\n", "59 1 1 5 total 0.000000 0.000000\n", "58 1 1 6 total 0.000000 0.000000\n", "57 1 1 7 total 0.000000 0.000000\n", "56 1 1 8 total 0.000000 0.000000\n", "55 1 2 1 total 0.000000 0.000000\n", - "54 1 2 2 total 0.266651 0.001340" + "54 1 2 2 total 0.266499 0.001265" ] }, "execution_count": 21, @@ -1020,133 +1031,133 @@ "[ NORMAL ] Ray tracing for track segmentation...\n", "[ NORMAL ] Dumping tracks to file...\n", "[ NORMAL ] Computing the eigenvalue...\n", - "[ NORMAL ] Iteration 0:\tk_eff = 0.685180\tres = 0.000E+00\n", - "[ NORMAL ] Iteration 1:\tk_eff = 0.785704\tres = 3.148E-01\n", - "[ NORMAL ] Iteration 2:\tk_eff = 0.750352\tres = 1.467E-01\n", - "[ NORMAL ] Iteration 3:\tk_eff = 0.729115\tres = 4.499E-02\n", - "[ NORMAL ] Iteration 4:\tk_eff = 0.696059\tres = 2.830E-02\n", - "[ NORMAL ] Iteration 5:\tk_eff = 0.663970\tres = 4.534E-02\n", - "[ NORMAL ] Iteration 6:\tk_eff = 0.633141\tres = 4.610E-02\n", - "[ NORMAL ] Iteration 7:\tk_eff = 0.605167\tres = 4.643E-02\n", - "[ NORMAL ] Iteration 8:\tk_eff = 0.580592\tres = 4.418E-02\n", - "[ NORMAL ] Iteration 9:\tk_eff = 0.559758\tres = 4.061E-02\n", - "[ NORMAL ] Iteration 10:\tk_eff = 0.542846\tres = 3.588E-02\n", - "[ NORMAL ] Iteration 11:\tk_eff = 0.529901\tres = 3.021E-02\n", - "[ NORMAL ] Iteration 12:\tk_eff = 0.520893\tres = 2.385E-02\n", - "[ NORMAL ] Iteration 13:\tk_eff = 0.515699\tres = 1.700E-02\n", - "[ NORMAL ] Iteration 14:\tk_eff = 0.514152\tres = 9.971E-03\n", - "[ NORMAL ] Iteration 15:\tk_eff = 0.516033\tres = 2.999E-03\n", - "[ NORMAL ] Iteration 16:\tk_eff = 0.521086\tres = 3.657E-03\n", - "[ NORMAL ] Iteration 17:\tk_eff = 0.529034\tres = 9.792E-03\n", - "[ NORMAL ] Iteration 18:\tk_eff = 0.539585\tres = 1.525E-02\n", - "[ NORMAL ] Iteration 19:\tk_eff = 0.552436\tres = 1.994E-02\n", - "[ NORMAL ] Iteration 20:\tk_eff = 0.567286\tres = 2.382E-02\n", - "[ NORMAL ] Iteration 21:\tk_eff = 0.583843\tres = 2.688E-02\n", - "[ NORMAL ] Iteration 22:\tk_eff = 0.601819\tres = 2.918E-02\n", - "[ NORMAL ] Iteration 23:\tk_eff = 0.620946\tres = 3.079E-02\n", - "[ NORMAL ] Iteration 24:\tk_eff = 0.640969\tres = 3.178E-02\n", - "[ NORMAL ] Iteration 25:\tk_eff = 0.661654\tres = 3.225E-02\n", - "[ NORMAL ] Iteration 26:\tk_eff = 0.682785\tres = 3.227E-02\n", - "[ NORMAL ] Iteration 27:\tk_eff = 0.704168\tres = 3.194E-02\n", - "[ NORMAL ] Iteration 28:\tk_eff = 0.725628\tres = 3.132E-02\n", - "[ NORMAL ] Iteration 29:\tk_eff = 0.747011\tres = 3.048E-02\n", - "[ NORMAL ] Iteration 30:\tk_eff = 0.768182\tres = 2.947E-02\n", - "[ NORMAL ] Iteration 31:\tk_eff = 0.789024\tres = 2.834E-02\n", - "[ NORMAL ] Iteration 32:\tk_eff = 0.809439\tres = 2.713E-02\n", - "[ NORMAL ] Iteration 33:\tk_eff = 0.829342\tres = 2.587E-02\n", - "[ NORMAL ] Iteration 34:\tk_eff = 0.848666\tres = 2.459E-02\n", - "[ NORMAL ] Iteration 35:\tk_eff = 0.867356\tres = 2.330E-02\n", - "[ NORMAL ] Iteration 36:\tk_eff = 0.885370\tres = 2.202E-02\n", - "[ NORMAL ] Iteration 37:\tk_eff = 0.902676\tres = 2.077E-02\n", - "[ NORMAL ] Iteration 38:\tk_eff = 0.919253\tres = 1.955E-02\n", - "[ NORMAL ] Iteration 39:\tk_eff = 0.935087\tres = 1.836E-02\n", - "[ NORMAL ] Iteration 40:\tk_eff = 0.950174\tres = 1.723E-02\n", - "[ NORMAL ] Iteration 41:\tk_eff = 0.964514\tres = 1.613E-02\n", - "[ NORMAL ] Iteration 42:\tk_eff = 0.978114\tres = 1.509E-02\n", - "[ NORMAL ] Iteration 43:\tk_eff = 0.990987\tres = 1.410E-02\n", - "[ NORMAL ] Iteration 44:\tk_eff = 1.003145\tres = 1.316E-02\n", - "[ NORMAL ] Iteration 45:\tk_eff = 1.014610\tres = 1.227E-02\n", - "[ NORMAL ] Iteration 46:\tk_eff = 1.025401\tres = 1.143E-02\n", - "[ NORMAL ] Iteration 47:\tk_eff = 1.035542\tres = 1.064E-02\n", - "[ NORMAL ] Iteration 48:\tk_eff = 1.045058\tres = 9.890E-03\n", - "[ NORMAL ] Iteration 49:\tk_eff = 1.053973\tres = 9.189E-03\n", - "[ NORMAL ] Iteration 50:\tk_eff = 1.062316\tres = 8.531E-03\n", - "[ NORMAL ] Iteration 51:\tk_eff = 1.070112\tres = 7.915E-03\n", - "[ NORMAL ] Iteration 52:\tk_eff = 1.077389\tres = 7.339E-03\n", - "[ NORMAL ] Iteration 53:\tk_eff = 1.084173\tres = 6.800E-03\n", - "[ NORMAL ] Iteration 54:\tk_eff = 1.090490\tres = 6.297E-03\n", - "[ NORMAL ] Iteration 55:\tk_eff = 1.096368\tres = 5.827E-03\n", - "[ NORMAL ] Iteration 56:\tk_eff = 1.101830\tres = 5.390E-03\n", - "[ NORMAL ] Iteration 57:\tk_eff = 1.106902\tres = 4.982E-03\n", - "[ NORMAL ] Iteration 58:\tk_eff = 1.111608\tres = 4.603E-03\n", - "[ NORMAL ] Iteration 59:\tk_eff = 1.115969\tres = 4.251E-03\n", - "[ NORMAL ] Iteration 60:\tk_eff = 1.120009\tres = 3.924E-03\n", - "[ NORMAL ] Iteration 61:\tk_eff = 1.123747\tres = 3.620E-03\n", - "[ NORMAL ] Iteration 62:\tk_eff = 1.127204\tres = 3.338E-03\n", - "[ NORMAL ] Iteration 63:\tk_eff = 1.130399\tres = 3.076E-03\n", - "[ NORMAL ] Iteration 64:\tk_eff = 1.133349\tres = 2.834E-03\n", - "[ NORMAL ] Iteration 65:\tk_eff = 1.136072\tres = 2.610E-03\n", - "[ NORMAL ] Iteration 66:\tk_eff = 1.138584\tres = 2.403E-03\n", - "[ NORMAL ] Iteration 67:\tk_eff = 1.140899\tres = 2.211E-03\n", - "[ NORMAL ] Iteration 68:\tk_eff = 1.143032\tres = 2.033E-03\n", - "[ NORMAL ] Iteration 69:\tk_eff = 1.144996\tres = 1.869E-03\n", - "[ NORMAL ] Iteration 70:\tk_eff = 1.146803\tres = 1.718E-03\n", - "[ NORMAL ] Iteration 71:\tk_eff = 1.148466\tres = 1.579E-03\n", - "[ NORMAL ] Iteration 72:\tk_eff = 1.149995\tres = 1.450E-03\n", - "[ NORMAL ] Iteration 73:\tk_eff = 1.151399\tres = 1.331E-03\n", - "[ NORMAL ] Iteration 74:\tk_eff = 1.152690\tres = 1.222E-03\n", - "[ NORMAL ] Iteration 75:\tk_eff = 1.153875\tres = 1.121E-03\n", - "[ NORMAL ] Iteration 76:\tk_eff = 1.154963\tres = 1.028E-03\n", - "[ NORMAL ] Iteration 77:\tk_eff = 1.155961\tres = 9.428E-04\n", - "[ NORMAL ] Iteration 78:\tk_eff = 1.156876\tres = 8.642E-04\n", - "[ NORMAL ] Iteration 79:\tk_eff = 1.157716\tres = 7.920E-04\n", - "[ NORMAL ] Iteration 80:\tk_eff = 1.158485\tres = 7.256E-04\n", - "[ NORMAL ] Iteration 81:\tk_eff = 1.159190\tres = 6.646E-04\n", - "[ NORMAL ] Iteration 82:\tk_eff = 1.159836\tres = 6.085E-04\n", - "[ NORMAL ] Iteration 83:\tk_eff = 1.160427\tres = 5.571E-04\n", - "[ NORMAL ] Iteration 84:\tk_eff = 1.160969\tres = 5.098E-04\n", - "[ NORMAL ] Iteration 85:\tk_eff = 1.161464\tres = 4.665E-04\n", - "[ NORMAL ] Iteration 86:\tk_eff = 1.161917\tres = 4.268E-04\n", - "[ NORMAL ] Iteration 87:\tk_eff = 1.162332\tres = 3.903E-04\n", - "[ NORMAL ] Iteration 88:\tk_eff = 1.162711\tres = 3.570E-04\n", - "[ NORMAL ] Iteration 89:\tk_eff = 1.163058\tres = 3.264E-04\n", - "[ NORMAL ] Iteration 90:\tk_eff = 1.163375\tres = 2.982E-04\n", - "[ NORMAL ] Iteration 91:\tk_eff = 1.163664\tres = 2.725E-04\n", - "[ NORMAL ] Iteration 92:\tk_eff = 1.163929\tres = 2.490E-04\n", - "[ NORMAL ] Iteration 93:\tk_eff = 1.164171\tres = 2.275E-04\n", - "[ NORMAL ] Iteration 94:\tk_eff = 1.164392\tres = 2.077E-04\n", - "[ NORMAL ] Iteration 95:\tk_eff = 1.164593\tres = 1.897E-04\n", - "[ NORMAL ] Iteration 96:\tk_eff = 1.164777\tres = 1.733E-04\n", - "[ NORMAL ] Iteration 97:\tk_eff = 1.164946\tres = 1.581E-04\n", - "[ NORMAL ] Iteration 98:\tk_eff = 1.165099\tres = 1.444E-04\n", - "[ NORMAL ] Iteration 99:\tk_eff = 1.165239\tres = 1.317E-04\n", - "[ NORMAL ] Iteration 100:\tk_eff = 1.165367\tres = 1.202E-04\n", - "[ NORMAL ] Iteration 101:\tk_eff = 1.165483\tres = 1.096E-04\n", - "[ NORMAL ] Iteration 102:\tk_eff = 1.165590\tres = 1.000E-04\n", - "[ NORMAL ] Iteration 103:\tk_eff = 1.165687\tres = 9.126E-05\n", - "[ NORMAL ] Iteration 104:\tk_eff = 1.165775\tres = 8.326E-05\n", - "[ NORMAL ] Iteration 105:\tk_eff = 1.165856\tres = 7.590E-05\n", - "[ NORMAL ] Iteration 106:\tk_eff = 1.165929\tres = 6.924E-05\n", - "[ NORMAL ] Iteration 107:\tk_eff = 1.165996\tres = 6.304E-05\n", - "[ NORMAL ] Iteration 108:\tk_eff = 1.166057\tres = 5.745E-05\n", - "[ NORMAL ] Iteration 109:\tk_eff = 1.166113\tres = 5.243E-05\n", - "[ NORMAL ] Iteration 110:\tk_eff = 1.166164\tres = 4.776E-05\n", - "[ NORMAL ] Iteration 111:\tk_eff = 1.166210\tres = 4.355E-05\n", - "[ NORMAL ] Iteration 112:\tk_eff = 1.166252\tres = 3.962E-05\n", - "[ NORMAL ] Iteration 113:\tk_eff = 1.166291\tres = 3.624E-05\n", - "[ NORMAL ] Iteration 114:\tk_eff = 1.166326\tres = 3.295E-05\n", - "[ NORMAL ] Iteration 115:\tk_eff = 1.166358\tres = 3.003E-05\n", - "[ NORMAL ] Iteration 116:\tk_eff = 1.166387\tres = 2.737E-05\n", - "[ NORMAL ] Iteration 117:\tk_eff = 1.166413\tres = 2.486E-05\n", - "[ NORMAL ] Iteration 118:\tk_eff = 1.166437\tres = 2.267E-05\n", - "[ NORMAL ] Iteration 119:\tk_eff = 1.166459\tres = 2.061E-05\n", - "[ NORMAL ] Iteration 120:\tk_eff = 1.166479\tres = 1.883E-05\n", - "[ NORMAL ] Iteration 121:\tk_eff = 1.166497\tres = 1.711E-05\n", - "[ NORMAL ] Iteration 122:\tk_eff = 1.166514\tres = 1.563E-05\n", - "[ NORMAL ] Iteration 123:\tk_eff = 1.166529\tres = 1.418E-05\n", - "[ NORMAL ] Iteration 124:\tk_eff = 1.166543\tres = 1.300E-05\n", - "[ NORMAL ] Iteration 125:\tk_eff = 1.166555\tres = 1.176E-05\n", - "[ NORMAL ] Iteration 126:\tk_eff = 1.166567\tres = 1.074E-05\n" + "[ NORMAL ] Iteration 0:\tk_eff = 0.685185\tres = 0.000E+00\n", + "[ NORMAL ] Iteration 1:\tk_eff = 0.785642\tres = 3.148E-01\n", + "[ NORMAL ] Iteration 2:\tk_eff = 0.750185\tres = 1.466E-01\n", + "[ NORMAL ] Iteration 3:\tk_eff = 0.728847\tres = 4.513E-02\n", + "[ NORMAL ] Iteration 4:\tk_eff = 0.695633\tres = 2.844E-02\n", + "[ NORMAL ] Iteration 5:\tk_eff = 0.663357\tres = 4.557E-02\n", + "[ NORMAL ] Iteration 6:\tk_eff = 0.632339\tres = 4.640E-02\n", + "[ NORMAL ] Iteration 7:\tk_eff = 0.604187\tres = 4.676E-02\n", + "[ NORMAL ] Iteration 8:\tk_eff = 0.579451\tres = 4.452E-02\n", + "[ NORMAL ] Iteration 9:\tk_eff = 0.558474\tres = 4.094E-02\n", + "[ NORMAL ] Iteration 10:\tk_eff = 0.541436\tres = 3.620E-02\n", + "[ NORMAL ] Iteration 11:\tk_eff = 0.528380\tres = 3.051E-02\n", + "[ NORMAL ] Iteration 12:\tk_eff = 0.519273\tres = 2.411E-02\n", + "[ NORMAL ] Iteration 13:\tk_eff = 0.513991\tres = 1.724E-02\n", + "[ NORMAL ] Iteration 14:\tk_eff = 0.512364\tres = 1.017E-02\n", + "[ NORMAL ] Iteration 15:\tk_eff = 0.514171\tres = 3.165E-03\n", + "[ NORMAL ] Iteration 16:\tk_eff = 0.519155\tres = 3.527E-03\n", + "[ NORMAL ] Iteration 17:\tk_eff = 0.527038\tres = 9.693E-03\n", + "[ NORMAL ] Iteration 18:\tk_eff = 0.537524\tres = 1.518E-02\n", + "[ NORMAL ] Iteration 19:\tk_eff = 0.550310\tres = 1.990E-02\n", + "[ NORMAL ] Iteration 20:\tk_eff = 0.565096\tres = 2.379E-02\n", + "[ NORMAL ] Iteration 21:\tk_eff = 0.581585\tres = 2.687E-02\n", + "[ NORMAL ] Iteration 22:\tk_eff = 0.599493\tres = 2.918E-02\n", + "[ NORMAL ] Iteration 23:\tk_eff = 0.618548\tres = 3.079E-02\n", + "[ NORMAL ] Iteration 24:\tk_eff = 0.638497\tres = 3.179E-02\n", + "[ NORMAL ] Iteration 25:\tk_eff = 0.659105\tres = 3.225E-02\n", + "[ NORMAL ] Iteration 26:\tk_eff = 0.680156\tres = 3.228E-02\n", + "[ NORMAL ] Iteration 27:\tk_eff = 0.701456\tres = 3.194E-02\n", + "[ NORMAL ] Iteration 28:\tk_eff = 0.722831\tres = 3.132E-02\n", + "[ NORMAL ] Iteration 29:\tk_eff = 0.744127\tres = 3.047E-02\n", + "[ NORMAL ] Iteration 30:\tk_eff = 0.765209\tres = 2.946E-02\n", + "[ NORMAL ] Iteration 31:\tk_eff = 0.785961\tres = 2.833E-02\n", + "[ NORMAL ] Iteration 32:\tk_eff = 0.806283\tres = 2.712E-02\n", + "[ NORMAL ] Iteration 33:\tk_eff = 0.826093\tres = 2.586E-02\n", + "[ NORMAL ] Iteration 34:\tk_eff = 0.845324\tres = 2.457E-02\n", + "[ NORMAL ] Iteration 35:\tk_eff = 0.863921\tres = 2.328E-02\n", + "[ NORMAL ] Iteration 36:\tk_eff = 0.881841\tres = 2.200E-02\n", + "[ NORMAL ] Iteration 37:\tk_eff = 0.899055\tres = 2.074E-02\n", + "[ NORMAL ] Iteration 38:\tk_eff = 0.915540\tres = 1.952E-02\n", + "[ NORMAL ] Iteration 39:\tk_eff = 0.931284\tres = 1.834E-02\n", + "[ NORMAL ] Iteration 40:\tk_eff = 0.946283\tres = 1.720E-02\n", + "[ NORMAL ] Iteration 41:\tk_eff = 0.960536\tres = 1.610E-02\n", + "[ NORMAL ] Iteration 42:\tk_eff = 0.974052\tres = 1.506E-02\n", + "[ NORMAL ] Iteration 43:\tk_eff = 0.986841\tres = 1.407E-02\n", + "[ NORMAL ] Iteration 44:\tk_eff = 0.998920\tres = 1.313E-02\n", + "[ NORMAL ] Iteration 45:\tk_eff = 1.010307\tres = 1.224E-02\n", + "[ NORMAL ] Iteration 46:\tk_eff = 1.021023\tres = 1.140E-02\n", + "[ NORMAL ] Iteration 47:\tk_eff = 1.031091\tres = 1.061E-02\n", + "[ NORMAL ] Iteration 48:\tk_eff = 1.040537\tres = 9.861E-03\n", + "[ NORMAL ] Iteration 49:\tk_eff = 1.049386\tres = 9.161E-03\n", + "[ NORMAL ] Iteration 50:\tk_eff = 1.057664\tres = 8.504E-03\n", + "[ NORMAL ] Iteration 51:\tk_eff = 1.065399\tres = 7.889E-03\n", + "[ NORMAL ] Iteration 52:\tk_eff = 1.072617\tres = 7.313E-03\n", + "[ NORMAL ] Iteration 53:\tk_eff = 1.079345\tres = 6.775E-03\n", + "[ NORMAL ] Iteration 54:\tk_eff = 1.085609\tres = 6.273E-03\n", + "[ NORMAL ] Iteration 55:\tk_eff = 1.091436\tres = 5.804E-03\n", + "[ NORMAL ] Iteration 56:\tk_eff = 1.096851\tres = 5.367E-03\n", + "[ NORMAL ] Iteration 57:\tk_eff = 1.101877\tres = 4.961E-03\n", + "[ NORMAL ] Iteration 58:\tk_eff = 1.106540\tres = 4.583E-03\n", + "[ NORMAL ] Iteration 59:\tk_eff = 1.110860\tres = 4.231E-03\n", + "[ NORMAL ] Iteration 60:\tk_eff = 1.114861\tres = 3.905E-03\n", + "[ NORMAL ] Iteration 61:\tk_eff = 1.118563\tres = 3.602E-03\n", + "[ NORMAL ] Iteration 62:\tk_eff = 1.121986\tres = 3.321E-03\n", + "[ NORMAL ] Iteration 63:\tk_eff = 1.125149\tres = 3.060E-03\n", + "[ NORMAL ] Iteration 64:\tk_eff = 1.128069\tres = 2.819E-03\n", + "[ NORMAL ] Iteration 65:\tk_eff = 1.130764\tres = 2.595E-03\n", + "[ NORMAL ] Iteration 66:\tk_eff = 1.133248\tres = 2.389E-03\n", + "[ NORMAL ] Iteration 67:\tk_eff = 1.135538\tres = 2.197E-03\n", + "[ NORMAL ] Iteration 68:\tk_eff = 1.137648\tres = 2.021E-03\n", + "[ NORMAL ] Iteration 69:\tk_eff = 1.139590\tres = 1.858E-03\n", + "[ NORMAL ] Iteration 70:\tk_eff = 1.141377\tres = 1.707E-03\n", + "[ NORMAL ] Iteration 71:\tk_eff = 1.143021\tres = 1.568E-03\n", + "[ NORMAL ] Iteration 72:\tk_eff = 1.144531\tres = 1.440E-03\n", + "[ NORMAL ] Iteration 73:\tk_eff = 1.145920\tres = 1.322E-03\n", + "[ NORMAL ] Iteration 74:\tk_eff = 1.147195\tres = 1.213E-03\n", + "[ NORMAL ] Iteration 75:\tk_eff = 1.148365\tres = 1.113E-03\n", + "[ NORMAL ] Iteration 76:\tk_eff = 1.149440\tres = 1.020E-03\n", + "[ NORMAL ] Iteration 77:\tk_eff = 1.150425\tres = 9.356E-04\n", + "[ NORMAL ] Iteration 78:\tk_eff = 1.151329\tres = 8.574E-04\n", + "[ NORMAL ] Iteration 79:\tk_eff = 1.152157\tres = 7.856E-04\n", + "[ NORMAL ] Iteration 80:\tk_eff = 1.152916\tres = 7.196E-04\n", + "[ NORMAL ] Iteration 81:\tk_eff = 1.153612\tres = 6.590E-04\n", + "[ NORMAL ] Iteration 82:\tk_eff = 1.154249\tres = 6.033E-04\n", + "[ NORMAL ] Iteration 83:\tk_eff = 1.154832\tres = 5.522E-04\n", + "[ NORMAL ] Iteration 84:\tk_eff = 1.155366\tres = 5.053E-04\n", + "[ NORMAL ] Iteration 85:\tk_eff = 1.155855\tres = 4.623E-04\n", + "[ NORMAL ] Iteration 86:\tk_eff = 1.156302\tres = 4.229E-04\n", + "[ NORMAL ] Iteration 87:\tk_eff = 1.156710\tres = 3.866E-04\n", + "[ NORMAL ] Iteration 88:\tk_eff = 1.157084\tres = 3.535E-04\n", + "[ NORMAL ] Iteration 89:\tk_eff = 1.157426\tres = 3.230E-04\n", + "[ NORMAL ] Iteration 90:\tk_eff = 1.157738\tres = 2.952E-04\n", + "[ NORMAL ] Iteration 91:\tk_eff = 1.158023\tres = 2.697E-04\n", + "[ NORMAL ] Iteration 92:\tk_eff = 1.158284\tres = 2.464E-04\n", + "[ NORMAL ] Iteration 93:\tk_eff = 1.158522\tres = 2.251E-04\n", + "[ NORMAL ] Iteration 94:\tk_eff = 1.158739\tres = 2.055E-04\n", + "[ NORMAL ] Iteration 95:\tk_eff = 1.158937\tres = 1.876E-04\n", + "[ NORMAL ] Iteration 96:\tk_eff = 1.159119\tres = 1.713E-04\n", + "[ NORMAL ] Iteration 97:\tk_eff = 1.159284\tres = 1.564E-04\n", + "[ NORMAL ] Iteration 98:\tk_eff = 1.159435\tres = 1.428E-04\n", + "[ NORMAL ] Iteration 99:\tk_eff = 1.159573\tres = 1.301E-04\n", + "[ NORMAL ] Iteration 100:\tk_eff = 1.159699\tres = 1.188E-04\n", + "[ NORMAL ] Iteration 101:\tk_eff = 1.159813\tres = 1.083E-04\n", + "[ NORMAL ] Iteration 102:\tk_eff = 1.159917\tres = 9.880E-05\n", + "[ NORMAL ] Iteration 103:\tk_eff = 1.160013\tres = 9.009E-05\n", + "[ NORMAL ] Iteration 104:\tk_eff = 1.160100\tres = 8.222E-05\n", + "[ NORMAL ] Iteration 105:\tk_eff = 1.160179\tres = 7.487E-05\n", + "[ NORMAL ] Iteration 106:\tk_eff = 1.160251\tres = 6.824E-05\n", + "[ NORMAL ] Iteration 107:\tk_eff = 1.160317\tres = 6.223E-05\n", + "[ NORMAL ] Iteration 108:\tk_eff = 1.160376\tres = 5.675E-05\n", + "[ NORMAL ] Iteration 109:\tk_eff = 1.160431\tres = 5.174E-05\n", + "[ NORMAL ] Iteration 110:\tk_eff = 1.160481\tres = 4.715E-05\n", + "[ NORMAL ] Iteration 111:\tk_eff = 1.160527\tres = 4.298E-05\n", + "[ NORMAL ] Iteration 112:\tk_eff = 1.160568\tres = 3.910E-05\n", + "[ NORMAL ] Iteration 113:\tk_eff = 1.160605\tres = 3.561E-05\n", + "[ NORMAL ] Iteration 114:\tk_eff = 1.160640\tres = 3.242E-05\n", + "[ NORMAL ] Iteration 115:\tk_eff = 1.160671\tres = 2.960E-05\n", + "[ NORMAL ] Iteration 116:\tk_eff = 1.160699\tres = 2.689E-05\n", + "[ NORMAL ] Iteration 117:\tk_eff = 1.160725\tres = 2.452E-05\n", + "[ NORMAL ] Iteration 118:\tk_eff = 1.160749\tres = 2.232E-05\n", + "[ NORMAL ] Iteration 119:\tk_eff = 1.160770\tres = 2.034E-05\n", + "[ NORMAL ] Iteration 120:\tk_eff = 1.160790\tres = 1.855E-05\n", + "[ NORMAL ] Iteration 121:\tk_eff = 1.160807\tres = 1.687E-05\n", + "[ NORMAL ] Iteration 122:\tk_eff = 1.160824\tres = 1.538E-05\n", + "[ NORMAL ] Iteration 123:\tk_eff = 1.160838\tres = 1.399E-05\n", + "[ NORMAL ] Iteration 124:\tk_eff = 1.160852\tres = 1.280E-05\n", + "[ NORMAL ] Iteration 125:\tk_eff = 1.160864\tres = 1.149E-05\n", + "[ NORMAL ] Iteration 126:\tk_eff = 1.160875\tres = 1.061E-05\n" ] } ], @@ -1179,9 +1190,9 @@ "name": "stdout", "output_type": "stream", "text": [ - "openmc keff = 1.163673\n", - "openmoc keff = 1.166567\n", - "bias [pcm]: 289.4\n" + "openmc keff = 1.161200\n", + "openmoc keff = 1.160875\n", + "bias [pcm]: -32.5\n" ] } ], @@ -1324,9 +1335,20 @@ "cell_type": "code", "execution_count": 31, "metadata": { - "collapsed": true + "collapsed": false }, - "outputs": [], + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/universe.py:223: DeprecationWarning: Cell.add_surface(...) has been deprecated and may be removed in a future version. The region for a Cell should be defined using the region property directly.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/universe.py:223: DeprecationWarning: Cell.add_surface(...) has been deprecated and may be removed in a future version. The region for a Cell should be defined using the region property directly.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/universe.py:223: DeprecationWarning: Cell.add_surface(...) has been deprecated and may be removed in a future version. The region for a Cell should be defined using the region property directly.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/universe.py:223: DeprecationWarning: Cell.add_surface(...) has been deprecated and may be removed in a future version. The region for a Cell should be defined using the region property directly.\n" + ] + } + ], "source": [ "# Create a Universe to encapsulate a fuel pin\n", "pin_cell_universe = openmc.Universe(name='1.6% Fuel Pin')\n", @@ -1364,7 +1386,18 @@ "metadata": { "collapsed": false }, - "outputs": [], + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/universe.py:223: DeprecationWarning: Cell.add_surface(...) has been deprecated and may be removed in a future version. The region for a Cell should be defined using the region property directly.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/universe.py:223: DeprecationWarning: Cell.add_surface(...) has been deprecated and may be removed in a future version. The region for a Cell should be defined using the region property directly.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/universe.py:223: DeprecationWarning: Cell.add_surface(...) has been deprecated and may be removed in a future version. The region for a Cell should be defined using the region property directly.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/universe.py:223: DeprecationWarning: Cell.add_surface(...) has been deprecated and may be removed in a future version. The region for a Cell should be defined using the region property directly.\n" + ] + } + ], "source": [ "# Create root Cell\n", "root_cell = openmc.Cell(name='root cell')\n", @@ -1600,25 +1633,25 @@ "\tDomain ID =\t10000\n", "\tNuclide =\tU-235\n", "\tCross Sections [barns]:\n", - " Group 1 [0.821 - 20.0 MeV]:\t3.31e+00 +/- 6.20e-01%\n", - " Group 2 [0.00553 - 0.821 MeV]:\t3.96e+00 +/- 3.40e-01%\n", - " Group 3 [4e-06 - 0.00553 MeV]:\t5.51e+01 +/- 5.07e-01%\n", - " Group 4 [6.25e-07 - 4e-06 MeV]:\t8.79e+01 +/- 7.27e-01%\n", - " Group 5 [2.8e-07 - 6.25e-07 MeV]:\t2.90e+02 +/- 1.13e+00%\n", - " Group 6 [1.4e-07 - 2.8e-07 MeV]:\t4.49e+02 +/- 1.11e+00%\n", - " Group 7 [5.8e-08 - 1.4e-07 MeV]:\t6.88e+02 +/- 9.03e-01%\n", - " Group 8 [0.0 - 5.8e-08 MeV]:\t1.44e+03 +/- 6.88e-01%\n", + " Group 1 [0.821 - 20.0 MeV]:\t3.30e+00 +/- 5.91e-01%\n", + " Group 2 [0.00553 - 0.821 MeV]:\t3.97e+00 +/- 4.03e-01%\n", + " Group 3 [4e-06 - 0.00553 MeV]:\t5.48e+01 +/- 5.56e-01%\n", + " Group 4 [6.25e-07 - 4e-06 MeV]:\t8.84e+01 +/- 8.48e-01%\n", + " Group 5 [2.8e-07 - 6.25e-07 MeV]:\t2.89e+02 +/- 1.25e+00%\n", + " Group 6 [1.4e-07 - 2.8e-07 MeV]:\t4.49e+02 +/- 1.09e+00%\n", + " Group 7 [5.8e-08 - 1.4e-07 MeV]:\t6.87e+02 +/- 7.98e-01%\n", + " Group 8 [0.0 - 5.8e-08 MeV]:\t1.44e+03 +/- 5.73e-01%\n", "\n", "\tNuclide =\tU-238\n", "\tCross Sections [barns]:\n", - " Group 1 [0.821 - 20.0 MeV]:\t1.07e+00 +/- 6.51e-01%\n", - " Group 2 [0.00553 - 0.821 MeV]:\t1.22e-03 +/- 6.61e-01%\n", - " Group 3 [4e-06 - 0.00553 MeV]:\t6.15e-04 +/- 9.95e+00%\n", - " Group 4 [6.25e-07 - 4e-06 MeV]:\t6.53e-06 +/- 6.29e-01%\n", - " Group 5 [2.8e-07 - 6.25e-07 MeV]:\t1.07e-05 +/- 1.08e+00%\n", - " Group 6 [1.4e-07 - 2.8e-07 MeV]:\t1.55e-05 +/- 1.11e+00%\n", - " Group 7 [5.8e-08 - 1.4e-07 MeV]:\t2.30e-05 +/- 9.03e-01%\n", - " Group 8 [0.0 - 5.8e-08 MeV]:\t4.25e-05 +/- 6.86e-01%\n", + " Group 1 [0.821 - 20.0 MeV]:\t1.06e+00 +/- 6.74e-01%\n", + " Group 2 [0.00553 - 0.821 MeV]:\t1.22e-03 +/- 8.28e-01%\n", + " Group 3 [4e-06 - 0.00553 MeV]:\t4.75e-04 +/- 7.97e+00%\n", + " Group 4 [6.25e-07 - 4e-06 MeV]:\t6.53e-06 +/- 7.56e-01%\n", + " Group 5 [2.8e-07 - 6.25e-07 MeV]:\t1.07e-05 +/- 1.22e+00%\n", + " Group 6 [1.4e-07 - 2.8e-07 MeV]:\t1.55e-05 +/- 1.09e+00%\n", + " Group 7 [5.8e-08 - 1.4e-07 MeV]:\t2.30e-05 +/- 7.97e-01%\n", + " Group 8 [0.0 - 5.8e-08 MeV]:\t4.25e-05 +/- 5.72e-01%\n", "\n", "\n", "\n" @@ -1653,14 +1686,14 @@ "\tDomain Type =\tcell\n", "\tDomain ID =\t10000\n", "\tCross Sections [cm^-1]:\n", - " Group 1 [0.821 - 20.0 MeV]:\t2.54e-02 +/- 6.20e-01%\n", - " Group 2 [0.00553 - 0.821 MeV]:\t1.51e-03 +/- 3.34e-01%\n", - " Group 3 [4e-06 - 0.00553 MeV]:\t2.07e-02 +/- 5.07e-01%\n", - " Group 4 [6.25e-07 - 4e-06 MeV]:\t3.30e-02 +/- 7.26e-01%\n", - " Group 5 [2.8e-07 - 6.25e-07 MeV]:\t1.09e-01 +/- 1.13e+00%\n", - " Group 6 [1.4e-07 - 2.8e-07 MeV]:\t1.69e-01 +/- 1.11e+00%\n", - " Group 7 [5.8e-08 - 1.4e-07 MeV]:\t2.58e-01 +/- 9.03e-01%\n", - " Group 8 [0.0 - 5.8e-08 MeV]:\t5.41e-01 +/- 6.88e-01%\n", + " Group 1 [0.821 - 20.0 MeV]:\t2.52e-02 +/- 6.42e-01%\n", + " Group 2 [0.00553 - 0.821 MeV]:\t1.52e-03 +/- 3.96e-01%\n", + " Group 3 [4e-06 - 0.00553 MeV]:\t2.06e-02 +/- 5.56e-01%\n", + " Group 4 [6.25e-07 - 4e-06 MeV]:\t3.32e-02 +/- 8.48e-01%\n", + " Group 5 [2.8e-07 - 6.25e-07 MeV]:\t1.09e-01 +/- 1.25e+00%\n", + " Group 6 [1.4e-07 - 2.8e-07 MeV]:\t1.69e-01 +/- 1.09e+00%\n", + " Group 7 [5.8e-08 - 1.4e-07 MeV]:\t2.58e-01 +/- 7.98e-01%\n", + " Group 8 [0.0 - 5.8e-08 MeV]:\t5.41e-01 +/- 5.73e-01%\n", "\n", "\n", "\n" @@ -1709,8 +1742,8 @@ " 1\n", " 1\n", " O-16\n", - " 1.570467\n", - " 0.018506\n", + " 1.560098\n", + " 0.017801\n", " \n", " \n", " 127\n", @@ -1718,8 +1751,8 @@ " 1\n", " 1\n", " H-1\n", - " 0.235674\n", - " 0.009063\n", + " 0.234877\n", + " 0.010096\n", " \n", " \n", " 124\n", @@ -1727,8 +1760,8 @@ " 1\n", " 2\n", " O-16\n", - " 0.288333\n", - " 0.003932\n", + " 0.288236\n", + " 0.004397\n", " \n", " \n", " 125\n", @@ -1736,8 +1769,8 @@ " 1\n", " 2\n", " H-1\n", - " 1.581295\n", - " 0.008248\n", + " 1.587815\n", + " 0.007847\n", " \n", " \n", " 122\n", @@ -1754,8 +1787,8 @@ " 1\n", " 3\n", " H-1\n", - " 0.010828\n", - " 0.000616\n", + " 0.010122\n", + " 0.000513\n", " \n", " \n", " 120\n", @@ -1799,12 +1832,12 @@ ], "text/plain": [ " cell group in group out nuclide mean std. dev.\n", - "126 10002 1 1 O-16 1.570467 0.018506\n", - "127 10002 1 1 H-1 0.235674 0.009063\n", - "124 10002 1 2 O-16 0.288333 0.003932\n", - "125 10002 1 2 H-1 1.581295 0.008248\n", + "126 10002 1 1 O-16 1.560098 0.017801\n", + "127 10002 1 1 H-1 0.234877 0.010096\n", + "124 10002 1 2 O-16 0.288236 0.004397\n", + "125 10002 1 2 H-1 1.587815 0.007847\n", "122 10002 1 3 O-16 0.000000 0.000000\n", - "123 10002 1 3 H-1 0.010828 0.000616\n", + "123 10002 1 3 H-1 0.010122 0.000513\n", "120 10002 1 4 O-16 0.000000 0.000000\n", "121 10002 1 4 H-1 0.000000 0.000000\n", "118 10002 1 5 O-16 0.000000 0.000000\n", @@ -1866,9 +1899,9 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWYAAADDCAYAAACxgLv/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFw5JREFUeJzt3XmYVNWZx/Fvs4mCgopxZWgSF3BcUBF3u6PG4EpmBo0a\nRyXqGJfRiXs0hiY6UTPGPY6ORsSVJIiOuzFqoyaCIosg6rjQCgooILiAitLzx3vKqq6u5VTVuVWn\n29/neerpru2tc6ve+95zz71VB0REREREREREREREREREREREREREquK/gV/WuhEBdJblkOTsBbxW\n60YE0FmWo40WYN+s244Dni3wnG2Ax4EPgdUerzECmAEsd895EqgvrZntNAF3ZN3WDBxfYdyk1GPv\n1bSs2/sBXwJzPeMcR+HP5tvuOGAW8BmwALgB6FPkOYcDf3fPeTrH/V2BS4D3gI+xzzBfzM2Ae7E8\nX+bacmwpC5BDPZY7XTJuO46486AZa/N2Wbff527f2zPOauC74ZpVui7FH5KIVncpxZfAePyK4ObA\nOODnWDIPBH4PfF3ia/oodTmyVeMzWBP4x4zrRwFvU3nbM9Uql2rtLOAy93cdYFdgAPAE0L3A85YA\nV7rn5jLGxdrVxT0a+DzPY+8A3gH+AVgP+FdgUSkLUUBdoDgA3QLGyqUVeB04JuO29YHdgA9KjFVo\nuZNejpqZC+yTddux+G2NN6d4j3kkML3A/V2AC4A3sd7IVGBTd981wLtYT3sqsKe7fTjwBbaB+ATr\njV8CfAWsdLdd6x47CFsxl2C7O4dlvPZt2C7/I8Cn2J7DbcDF7v5GYD5wJrZyvY/1VFLWBx507XvB\ntSHf+1aPvVcXAL/NuP1Fd1tmj/n8jPfjFeBH7vbBbvm+csu41HM5zgMmYz0/gJOB2UCPPG3tiNbB\n3pORWbf3wgrBKI8YJ9C+x7yuizvQsx2f0L6XmGlPrHf+EZbbqd70Qdh6stzdPjrjOe9iufMJlhO7\nYhuG7DxYA7gC2zAsxHKip7uvEcvlc7E9iXHutnkZr9OCbdRmYr398S5myrnYOjAfe68K9WafBi5y\n8VOF9TRsD2Ye6R7zMOB59368D1xHeiP6jHuNT91yHuaxHN/D1vUd3PVNsL0X3x56NOZS+lBGik9h\nHogVkyuxN7B31v3nAC8DW7jr22I9DYCfYCtGF6w4LiBdTEYDt2fFehr4acb1XtgHdqyLMQT7kAa7\n+2/DEnA3d30NYCzwa3e9EViFDZt0BQ7AdndTu7Hjgbux5B+MrUDP5HwX0oV5gHtcHbA18Cr2/mcW\n5pHARu7/w7HE3NBdz7XRLLYcdcAk7D3bAluRt8/Tzo5qOPZZ5dpbuA37nIrJVZj3xopGqhC8DpxS\nIMYTwHPAj7Fec6YBWGH9MZZP65H+HBpI70ltixXWERnPyx7KyJUHVwH3A32x9ewB4Dfuvkbs/bkU\nK3w9aV+Y52Ib8I2w9W4OcJK7bzi2/IOxvb47sb3eQoX5eGzIc7i7bQq2UckszDtixbmLW845wBkZ\ncbKLv89ynIB1aNZ0r5/ZESpZrXY/67AP86OMy+8Jt2s9F3vjNgX+hBXGsVjRBPvwLgTecNdnke4B\n3OXasxor7GsAW2W0O9cuTuZtB7vXH+dizAAm0rbXfD+2xQbrhWfHWIUVuK+BR7EiuRW2Yv0zVuw+\nxwrsuDxtyjQfW7l/gO3mZW9cACZgKybYe/YGsEuOtqW0FlmOVvdapwP/C1yO9Yo6k37AYnJ3FBa6\n+8uxGbYh3gLbuI7ENtT75Xn8YVjBvAgbopoODHX3HYUV7j9i+bSU9OcwCSsmYOvAeKxYQ/E8T10/\nEevALMPy9FLgiIzHrMbydRX5h2Kuxd6vj7C9wSHu9sOBW7E8X+ni+Ayt3I7l3iBsgzE56/5p2N7m\naqyn/z+klzufYstxC7bH+QLWobnQo5151XKMeQS2hUxdTiH9pv8E2434BHi4zNeYgvUSvoMdQd2b\n9JvVH3grz/POxragy7BE6UPxFSxzgzIAK2iZG52jSPc+W2m7pc1lCW1X9hVYb2QDbHwr8/nzi8RK\nvebt2K71EdiYZHaCH4Ot0Kk2b4MNmxRSbDnewQ7IDMA2vJ3NYiw3cq1HG2MdAoAbSefz+R5xV7q/\nv8Y2eKmieWCexy8DfoF9ZhtinYH73X39sWKdyy5YL/MDF+Mkin/mmTYA1gJeIp03j9J2ffkQG/4r\nZGHG/ytJd6A2prxcn4gNlZ5K7k7IlsBDWG98OfCfFF9un+W4BdsDuQ4r4GWL6YBNZqG4C1jbXQ4K\nEHsqdmQ2tds2DxsSybYXNsxxGLalXRf74DJ7gdmyb3sX64lkbnTWxpKkEJ+9hQ+xMb7+Gbf1z/PY\nbBOxFfst2if4AKzXcCq2q7suNh5caLl9HITtRj6JjUN2Ns9jhfNfsm7vje1KP+mu/4x0Pmcf7Mv1\n3r6c5/V8PoclwO+wcc71sHz8Xp7H3o0V8M2wfL+RdE3wyfXFWCHdmnSu98XG3ktpcz4LKC/XV2Ib\niJ/R/iwqsHHwOVgN6IN12IrVwmLL0Ru4GivOY7D3omwxFWYfPUmP965B24MEmfbAxnw2cNcHAYeQ\n3qW5BTtItTlWfLbDkrg3VvgWu9f5FW2TbCG2a5m5EVlE28R/CNsiH42NR3UHdnZtgPy7iD67aF9j\nBbYJG8sahB2B90n+z4DvY+9Ltl4uxmIsJ0Zhva+URdjKm3mWQbFd3X7Azdiw0XHY+3+ARzs7kuXY\nSngd8EPs/anHhoLmkbsopHTB8rm7+38N0u/vW9jQxIVYHg7G9v4eyhPrcqzT0Q0r/idjQ1FLseK7\nH9bZ6Ib1DFNjzL2xXu6X2JjrUaRzKXVaamZuL6RtHqzGPuOrSa9rmwL7F1huH6k8+hOWi4OwnvlF\nJcS4ABueeDfHfb2xvZcVLvbJWfdnr9M+rsGGMf4N28u/scTntxFTYS52Cl099kbOdo9biY095bIM\nOBTbBfwE23pOJD0gfyX2of8FW7luxlaSx4HHgP/DjhavpO0H+2f3dwnWCwf7QEZiK8HV2Djb/tiQ\nwXvYVv9S0huUXMuZfVuh9+E0bCu/EBtfvofCu1iZsabR9oBf6r45WC/reRd3G+xgUsqT2FjkQtKn\nHRVbjpuw3thj2HtzPLZBrKgnEaH/worAFVguTcaGcPal8O7sMVg+34Dtqa3E3rOUI7E9mSVYQf4l\nuc93BttI34cV2bewnuWh7r53sT2ls1ys6aTP4DgFGy75GCt6f8yIuQLbxf+bizsMeIr2eXAeNrY6\n2S3/E1jHJMWn5519X+r+x7Dx56exdTL7eEYhC7AzUXI5G9sIfYztKY7PalMTtm59hK3b+WpT6rYR\n2DqfKvBnYgcYj/Rop3RSl2MHNkU6u8HYHm1MHUoRwM7O2A7b3RuG7XIeWvAZIh3XP2HDPOtip+JN\nrG1zRHIbio0ffoYdbT+vts0RSdSj2NDkEuxr5xsWfriIiIiIiEhxOzakjljqoksyl+0bWqmBPg3b\n1H7ZdenElwGt5BHil6NamZo3ftpNTXBSU8GHDNjJ7+dNlzVdT9+m0wo+5p1Jgwre/42xTTCqcLu8\nKVYysRrrIOyvnPlq3af1waIPervpbr7bdFTBxzxVt1XB+9Ouxb7FXshdnrGasV8mCEGxwscaA3ny\nWqediIhERoVZRCQy1SvMOzUGC9WzcViwWAxpVKzOEKuG1m3cNmC0XYo/xFu9YnXQWD7jdsOxrxp3\nxb5Se3nW/X5jzB58x5h9eI8xS/ySG2Mumts+Y8w+/MeYffiOMUvcyh9j7gpcjyXw1th3vwcXfIZI\nx6DclmgVK8zDsB8oacF+kGU86RkORDoy5bZEq1hh3pT2P1S9aZ7HinQkym2JVrHZXv0Gj29qSv+/\nUyMMbSyzOSLA9GaY0Zz0q3jl9ttN6Wn71m3cNvCBPvl2aXGX4ooV5vdoP4NA++ldinxxRKQkOzTa\nJWXcmCRexSu3i31xRMRfPW3P1piU95HFhjKmkp4Qsgc2i8IDlTRNJBLKbYlWsR7zV9iMGY9jR7H/\nQP5ZQ0Q6EuW2RKtYYQb7PdRHk26ISA0otyVK+kq2iEhkVJhFRCKjwiwiEhkVZhGRyPgc/Cvu0yBR\nWJtPwgQCdmt4KlgsgOcn7RM0nnQMT9VNDRJnNIcEiQMwht8FiwUfB4wloajHLCISGRVmEZHIqDCL\niERGhVlEJDIqzCIikfEpzLcCi4BZCbdFpNqU2xIln8I8Fpt+R6SzUW5LlHwK87PAR0k3RKQGlNsS\nJY0xi4hEJsw3/8Y2pf8f0th29gmRUlVnailPzRn/19N2BgqRUrQQamopP6OagoQRAao1tZSnxhq+\ntnQu9YSaWkpERKrMpzDfA/wd2BKb7n1Uoi0SqR7ltkTJZyjjyMRbIVIbym2JkoYyREQio8IsIhIZ\nFWYRkcioMIuIRCbMecyBzJ60c7BYFzecHSwWwBcNPYLFmjZlz2Cx+CpcqKCxpI0xjA4WazRnBYul\naaripB6ziEhkVJhFRCKjwiwiEhkVZhGRyKgwi4hExqcw9weeBl4BZgOnJ9oikepRbkuUfE6XWwX8\nHJgB9AZeAp4AXk2wXSLVoNyWKPn0mBdiiQvwKZa0myTWIpHqUW5LlEodY64HdgCmhG+KSE3Vo9yW\nSJTyzb/ewATgDKx3kaappSSk6k8tlT+3NbWUBNNC6KmlugP3AncC97e7V1NLSUjVnVqqcG5raikJ\npp6QU0vVAX8A5gBXV9AqkdgotyVKPoV5D+Bo4PvAdHcZnmSjRKpEuS1R8hnKeA59EUU6J+W2RElJ\nKSISGRVmEZHIqDCLiERGhVlEJDJ1AWK00twaIEzkAp5Mdf99PwwW6zdcECzWm19vHiwWwNL53wkT\nqL4HhMnVUrUScEqoGE0l3DniQ5kaLJZ5MHC82IyBPHmtHrOISGRUmEVEIqPCLCISGRVmEZHIqDCL\niETGpzD3xH6jdgb2Yy+XJtoikepQXku0fH4r43PsR15WuMc/B+zp/op0VMpriZbvUMYK97cH0BVY\nmkxzRKpKeS1R8i3MXbBdvkXYrMJzEmuRSPUoryVKvjOYrAaGAH2Ax7FpHZq/uVdTS0lIz0+Cyfln\ndwiocF4DmlpKwmkh9NRSKcuBh4GhZGasppaSkHZrsEvKNZck/Yq58xrQ1FISTj0hp5bqB/R1/68J\n/ACb6UGkI1NeS7R8eswbA+OwIt4FuAN4MslGiVSB8lqi5VOYZwE7Jt0QkSpTXku09M0/EZHIqDCL\niERGhVlEJDIqzCIikVFhFhGJTKlfMPn2+o9woX5Ut1uwWKuXNAaLddl6ARcSeGTAgUHi6FeFkjOU\ny4LFam0YGiwWQN07AecSbWkKF6sK1GMWEYmMCrOISGRUmEVEIqPCLCISGd/C3BX7gZcHE2yLSC0o\ntyU6voX5DOxHxAMeJhWJgnJbouNTmDcDDgRuAeqSbY5IVSm3JUo+hfkq4BxstgeRzkS5LVEq9gWT\ng4EPsDG4xryP0tRSEtCy5pksb56Z9Mv45bamlpJgWgg1tdTuwKHY7l5PYB3gduCYNo/S1FISUN/G\n7enbuP031+eNuTOJl/HLbU0tJcHUE2pqqQuA/sBA4AjgKdolrkiHpNyWaJV6HrOOXEtnpdyWaJTy\nI0aTKNT3Fum4lNsSFX3zT0QkMirMIiKRUWEWEYmMCrOISGRUmEVEIqOppWphalOwUF3WXy9YrNab\nzwgWC6D+hLlB4mhqqSStDBapbtINwWIBtJ4V7udL6r4IeDbk9U3hYuWhHrOISGRUmEVEIqPCLCIS\nGRVmEZHI+B78awE+Br4GVgHDkmqQSBW1oLyWCPkW5lbs9w+XJtcUkapTXkuUShnK0NQ70hkpryU6\nvoW5FfgrMBU4MbnmiFSV8lqi5DuUsQewANgAeAJ4DXj2m3s1tZQENKd5MXOaF1fjpQrnNaCppSSc\nFkJNLZWywP39ELgPO0iSTmBNLSUBbd3Yj60b+31zfeKY15N6qcJ5DWhqKQmnnlBTSwGsBazt/u8F\n7A/MKq9hItFQXku0fHrMG2K9idTj7wL+kliLRKpDeS3R8inMc4EhSTdEpMqU1xItffNPRCQyKswi\nIpFRYRYRiYwKs4hIZFSYRUQiE+J3AlppDjhti5Sk55Bwv7/z+RXhpqkCaJ0Q5mco6l6zP0GClaYV\nRtfgZcXsGCxS67kjgsWqmxOo3j1UB3nyWj1mEZHIqDCLiERGhVlEJDIqzCIikfEpzH2BCcCrwBxg\n10RbJFI9ym2Jks9vZVwDPAKMdI/vlWiLRKpHuS1RKlaY+wB7Ace6618ByxNtkUh1KLclWsWGMgZi\nPyI+FpgG3Iz9jq1IR6fclmgV6zF3w87yPg14EbgaOB/4VZtHaWopCaj5M2hekfjL+OW2ppaSUBY3\nw5Jmr4cWK8zz3eVFd30ClrxtaWopCaixl11SxixJ5GX8cltTS0ko/RrtkvLGmLwPLTaUsRCYB2zp\nru8HvFJJ20QiodyWaPmclfHv2LQ7PYC3gFGJtkikepTbEiWfwjwT2DnphojUgHJboqRv/omIREaF\nWUQkMirMIiKRUWEWEYmMCrOISGR8zsqQiH3+WrjpoJ68ePdgsQCaLgkaTr51pgWLVPfbL4PFmhBo\nlrORBe5Tj1lEJDIqzCIikVFhFhGJjAqziEhkfArzVsD0jMty4PQkGyVSBcpriZbPWRmvAzu4/7sA\n7wH3JdYikepQXku0Sh3K2A/7Fa55CbRFpFaU1xKVUgvzEcDdSTREpIaU1xKVUr5g0gM4BDiv3T2a\nWkoCanGXKsmf14CmlpJQZuM/E0MphfkA4CVsAsu2NLWUBFRP2/I3KdmXy5/XgKaWklC2cZeUPxd4\nbClDGUcC95TVIpF4Ka8lOr6FuRd2gGRigm0RqTbltUTJdyjjM6Bfkg0RqQHltURJ3/wTEYlM9Qrz\n9GbFqlWsl8LFmtG8PFislmCRaq1FsTpFrHCHmWdX+PzqFeYZzYpVq1jTwsWa2fxxsFgtwSLVWoti\ndYpY4Qqz72lx+WgoQ0QkMirMIiKRCTFHSjPQECCOSD6TqM03PZpRbktyapXXIiIiIiIiIiId0XDg\nNeAN8v6Kl5dbgUXArABt6g88jZ3ZMpvKZq/oCUwBZgBzgEsrbFtXbFaNByuMA3Y+0csu3gsVxuoL\nTABexZZz1zLjdJbZQ5TXpQuV2y0oryvSFXgT+8Gw7tiHPLjMWHths06ESOCNgCHu/97YjBbltgtg\nLfe3GzAZ2LOCWGcCdwEPVBAjZS6wXoA4AOOAn7r/uwF9AsTsAizACkpHorwuT6jc7tR5XY3T5YZh\nCdwCrALGAyPKjPUs8FGYZrEQW5kAPsW2lptUEG+F+9sDW2mXlhlnM+BA4BbCnDVDoDh9sAJyq7v+\nFdYjqFRHnT1EeV260LndafO6GoV5U9o2br67LSb1WI9lSgUxumArxCJsV3JOmXGuAs4BVlfQlkyt\nwF+BqcCJFcQZiP1m8VhgGnAz6d5UJTrq7CHK69KFzO1OndfVKMytVXiNSvTGxpfOwHoY5VqN7UJu\nBuxNeecnHgx8gI1Pheot74GtnAcAp2K9g3J0A3YEbnB/PwPOr7BtqdlDCv1meKyU16UJndudOq+r\nUZjfo+04S3+sdxGD7sC9wJ3A/YFiLgceBoaW8dzdgUOx8bN7gH2A2ytszwL390NsFuhhZcaZ7y4v\nuusTsESuRJHZQ6KmvC5N6NxWXleoGzbWUo9tSSo5SIKLE+IgSR2WGFcFiNUPO7ILsCbwDLBvhTEb\nqPzI9VrA2u7/XsDfgP0riPcMsKX7vwm4vIJYYOOyx1YYo1aU1+WrNLeV14EcgB0dfhP4RQVx7gHe\nB77AxvdGVRBrT2w3bQbp01uGlxlrW2x8agZ2Cs85FbQrpYHKj1wPxNo0Azt1qpL3HmB7rGcxE5v1\no5Kj172AxaRXsI5IeV2eSnNbeS0iIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiKd2/8DY4XfFd+LX2AA\nAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWYAAADDCAYAAACxgLv/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFxZJREFUeJzt3XmYVNWZx/FvAyIKCijGBRiaxA0VRVTcpROX4BLNzKhR\n46jEOMZldOIel9gkJsRM4q6joxFxJRHRcVc0NGpcEVAEdVzoCAooi7iAcaHmj/cUVV3Ucqrq3KrT\nze/zPPV01/bWuVXvfe+5596qAyIiIiIiIiIiIiIiIiIiIiIiIiIiNfHfwIX1bkQAHWU5JDl7Am/U\nuxEBdJTlaKMV2DvntuOAp4s8ZxvgMeAjYIXHaxwCTAeWuuc8CTSW18xVNAO35dzWAhxfZdykNGLv\n1dSc2/sAXwKzPeMcR/HPZnV3HDAD+ByYB1wH9CzxnMOBZ91zJuW5vzNwCfA+8An2GRaK2Q+4B8vz\nj11bji1nAfJoxHKnU9ZtxxF3HrRgbd425/Z73e17ecZZAXw7XLPK16n0QxKRcpdyfAmMw68IbgqM\nBX6OJfNA4FrgmzJf00e5y5GrFp/BWsDWWdePAt6l+rZnq1cu1duZwO/c33WBXYABwERgjSLPWwRc\n5p6bzygXaxcX92jgiwKPvQ34O/BPwHrAvwELylmIIhoCxQHoEjBWPingTeCYrNvWB3YFPiwzVrHl\nTno56mY28L2c247Fb2u8KaV7zIcC04rc3wk4H3gb641MAfq6+64E3sN62lOAPdztI4B/YBuIT7He\n+CXA18Byd9tV7rFbYivmImx357Cs174F2+V/GPgM23O4Bfi1u78JmAucga1cH2A9lbT1gQdc+150\nbSj0vjVi79X5wO+zbn/J3ZbdYz4v6/2YCfzQ3T7ILd/XbhkXey7HucDzWM8P4CTgNaBrgba2R+ti\n78mhObd3xwrBSI8YP2XVHnNvF3egZzs+ZdVeYrY9sN75Eiy3073pA7H1ZKm7/eKs57yH5c6nWE7s\ngm0YcvNgTeAP2IZhPpYT3dx9TVgun4PtSYx1t83Jep1WbKP2CtbbH+dipp2DrQNzsfeqWG92EnCR\ni58urKdiezBzyPSYhwHPuffjA+BqMhvRp9xrfOaW8zCP5fgOtq5v765vgu29+PbQozGb8ocy0nwK\n80CsmFyGvYE9cu4/G3gV2MxdH4z1NAB+jK0YnbDiOI9MMbkYuDUn1iTgJ1nXu2Mf2LEuxhDsQxrk\n7r8FS8Bd3fU1gTHAr9z1JuArbNikM7A/trub3o0dB9yJJf8gbAV6Ku+7kCnMA9zjGoCtgNex9z+7\nMB8KbOT+PxxLzA3d9XwbzVLL0QBMxt6zzbAVebsC7WyvRmCfVb69hVuwz6mUfIV5L6xopAvBm8DJ\nRWJMBJ4BfoT1mrMNwArrj7B8Wo/M5zCczJ7UYKywHpL1vNyhjHx5cDlwH9ALW8/uB37r7mvC3p/R\nWOHrxqqFeTa2Ad8IW+9mASe6+0Zgyz8I2+u7HdvrLVaYj8eGPEe4217ANirZhXkoVpw7ueWcBZye\nFSe3+Pssx0+xDs1a7vWzO0Jlq9fuZwP2YS7JulxLuF3r2dgb1xf4C1YYx2BFE+zDuwB4y12fQaYH\ncIdrzwqssK8JbJHV7ny7ONm3HeRef6yLMR2YQNte833YFhusF54b4yuswH0DPIIVyS2wFetfsGL3\nBVZgxxZoU7a52Mq9L7abl7txARiPrZhg79lbwM552paWKrEcKfdapwH/C1yK9Yo6kj7AQvJ3FOa7\n+yvRD9sQb4ZtXA/FNtT7FHj8YVjBvAgbopoG7OjuOwor3H/G8mkxmc9hMlZMwNaBcVixhtJ5nr5+\nAtaB+RjL09HAEVmPWYHl61cUHoq5Cnu/lmB7g0Pc7YcDN2N5vtzF8RlauRXLvS2xDcbzOfdPxfY2\nV2A9/f8hs9yFlFqOm7A9zhexDs0FHu0sqJ5jzIdgW8j05WQyb/qPsd2IT4GHKnyNF7BewrewI6h7\nkXmz+gPvFHjeWdgW9GMsUXpSegXL3qAMwApa9kbnKDK9zxRtt7T5LKLtyr4M641sgI1vZT9/bolY\n6de8Fdu1PgIbk8xN8GOwFTrd5m2wYZNiSi3H37EDMgOwDW9HsxDLjXzr0cZYhwDgejL5fJ5H3OXu\n76+wDV66aB5Q4PEfA7/APrMNsc7Afe6+/lixzmdnrJf5oYtxIqU/82wbAGsDL5PJm0dou758hA3/\nFTM/6//lZDpQG1NZrk/AhkpPIX8nZHPgQaw3vhT4DaWX22c5bsL2QK7GCnjFYjpgk10o7gDWcZcD\nA8Segh2ZTe+2zcGGRHLtiQ1zHIZtaXtjH1x2LzBX7m3vYT2R7I3OOliSFOOzt/ARNsbXP+u2/gUe\nm2sCtmK/w6oJPgDrNZyC7er2xsaDiy23jwOx3cgnsXHIjuY5rHD+a87tPbBd6Sfd9Z+Ryefcg335\n3ttXC7yez+ewCPgjNs65HpaP3ynw2DuxAt4Py/frydQEn1xfiBXSrcjkei9s7L2cNhcyj8pyfTm2\ngfgZq55FBTYOPgurAT2xDlupWlhqOXoAV2DFeRT2XlQspsLsoxuZ8d41aXuQINvu2JjPBu76lsAP\nyOzS3IQdpNoUKz7bYkncAyt8C93r/JK2STYf27XM3ogsoG3iP4htkY/GxqPWAHZybYDCu4g+u2jf\nYAW2GRvL2hI7Au+T/J8D38Xel1zdXYyFWE6MxHpfaQuwlTf7LINSu7p9gBuxYaPjsPd/f492tidL\nsZXwauD72PvTiA0FzSF/UUjrhOXzGu7/Ncm8v+9gQxMXYHk4CNv7e7BArEuxTkcXrPifhA1FLcaK\n7z5YZ6ML1jNMjzH3wHq5X2JjrkeRyaX0aanZuT2ftnmwAvuMryCzrvUF9iuy3D7SefQXLBe3xHrm\nF5UR43xseOK9PPf1wPZelrnYJ+Xcn7tO+7gSG8b4d2wv//oyn99GTIW51Cl0jdgb+Zp73HJs7Cmf\nj4GDsV3AT7Gt5wQyA/KXYR/649jKdSO2kjwGPAr8H3a0eDltP9i73d9FWC8c7AM5FFsJrsDG2fbD\nhgzex7b6o8lsUPItZ+5txd6HU7Gt/HxsfPkuiu9iZceaStsDfun7ZmG9rOdc3G2wg0lpT2JjkfPJ\nnHZUajluwHpjj2LvzfHYBrGqnkSE/gsrAn/Acul5bAhnb4rvzh6D5fN12J7acuw9SzsS25NZhBXk\nC8l/vjPYRvperMi+g/UsD3b3vYftKZ3pYk0jcwbHydhwySdY0ftzVsxl2C7+31zcYcBfWTUPzsXG\nVp93yz8R65ik+fS8c+9L3/8oNv48CVsnc49nFDMPOxMln7OwjdAn2J7iuJw2NWPr1hJs3S5Um9K3\nHYKt8+kCfwZ2gPFIj3ZKB3UpdmBTpKMbhO3RxtShFAHs7Ixtsd29Ydgu58FFnyHSfv0zNszTGzsV\nb0J9myOS347Y+OHn2NH2c+vbHJFEPYINTS7Cvna+YfGHi4iIiIiIlDZ0ePqIpS66JHPZbniKOug5\nfHD9l12XDnwZkKKAEL8clWJKwfgZNzTDic1FH/LtHWYWvT9tcfO1rNdc/Psa707euuj9K41phpHF\n2+VNsZKJ1dQAYX/lzFdqeOqRkg9qbb6dxuajiz5mcoPvabFXYd9iL+YOz1gt2C8ThKBY4WONggJ5\nrdNOREQio8IsIhKZ2hXmHZqChVqraadgsRjSpFgdIVYd9Woq9lPI5dq59EO8NSpWO43lM243Avuq\ncWfsK7WX5tzvN8bswXeM2Yf3GLPEL7kx5pK57TPG7MN/jNmH7xizxK3yMebOwDVYAm+Fffd7UNFn\niLQPym2JVqnCPAz7gZJW7AdZxpGZ4UCkPVNuS7RKFea+rPpD1X0LPFakPVFuS7RKzfbqN3h8Q3Pm\n/x2aYMemCpsjAkxrgektSb+KV263Nt++8v9eTdsGPtAnq5dWdymtVGF+n1VnEFh1epcSXxwRKcv2\nTXZJGzsqiVfxyu1SXxwR8ddI27M1Jhd8ZKmhjClkJoTsis2icH81TROJhHJbolWqx/w1NmPGY9hR\n7D9ReNYQkfZEuS3RKlWYwX4PNczJnCJxUW5LlPSVbBGRyKgwi4hERoVZRCQyKswiIpHxOfhX2mdB\norAOn4YJBOwxfGKwWADPTN43aDxpHyY3PB8kzsXsHyQOwCj+GCwWfBIwloSiHrOISGRUmEVEIqPC\nLCISGRVmEZHIqDCLiETGpzDfDCwAZiTcFpFaU25LlHwK8xhs+h2Rjka5LVHyKcxPA0uSbohIHSi3\nJUoaYxYRiUyYb/6Nac78P6Sp7ewTIuWqzdRSnlqy/m+k7QwUIuVoJdTUUn5GNgcJIwLUamopT011\nfG3pWBoJNbWUiIjUmE9hvgt4Ftgcm+59ZKItEqkd5bZEyWco48jEWyFSH8ptiZKGMkREIqPCLCIS\nGRVmEZHIqDCLiEQmzHnMgbwyeZdgsX49/KxgsQC+HN41WKwXXx4eLBZfhAslyRnFxcFi/Y4zg8U6\njyuDxYLFAWOt3tRjFhGJjAqziEhkVJhFRCKjwiwiEhkVZhGRyPgU5v7AJGAm8BpwWqItEqkd5bZE\nyed0ua+AnwPTgR7Ay8BE4PUE2yVSC8ptiZJPj3k+lrgAn2FJu0liLRKpHeW2RKncMeZGYHvghfBN\nEamrRpTbEolyvvnXAxgPnI71LjI0tZSENLXFppeqncK5ramlJJhWQk8ttQZwD3A7cN8q92pqKQlp\naJNd0sYkOrVU8dzW1FISTCMhp5ZqAP4EzAKuqKJVIrFRbkuUfArz7sDRwHeBae4yIslGidSIclui\n5DOU8Qz6Iop0TMptiZKSUkQkMirMIiKRUWEWEYmMCrOISGQaAsRI0ZIKECZy14QLdd/d3w8W67ec\nHyzW299sGiwWwOL564cJ1G8tCJOr5UoRcEqoGD1NuHPE9+SVYLHMhMDxYjMKCuS1eswiIpFRYRYR\niYwKs4hIZFSYRUQio8IsIhIZn8LcDfuN2unYj72MTrRFIrWhvJZo+fxWxhfYj7wsc49/BtjD/RVp\nr5TXEi3foYxl7m9XoDOwOJnmiNSU8lqi5FuYO2G7fAuwWYVnJdYikdpRXkuUfGcwWQEMAXoCj2HT\nOrSsvFdTS0lIzz4Fzz1Vi1cqnteAppaScFoJPbVU2lLgIWBHsjNWU0tJSLvtZZe0y3+T9Cvmz2tA\nU0tJOI2EnFqqD9DL/b8WsC8204NIe6a8lmj59Jg3BsZiRbwTcBvwZJKNEqkB5bVEy6cwzwCGJt0Q\nkRpTXku09M0/EZHIqDCLiERGhVlEJDIqzCIikVFhFhGJTLlfMFl9/We4UD9s2DVYrNSCpmCxRn8r\n4EICT/TdO0icvwaJIvnsySXBYqUO2C5YLICGdwPOJfpGc7hYNaAes4hIZFSYRUQio8IsIhIZFWYR\nkcj4FubO2A+8PJBgW0TqQbkt0fEtzKdjPyIe8DCpSBSU2xIdn8LcDzgAuAloSLY5IjWl3JYo+RTm\ny4GzsdkeRDoS5bZEqdQXTA4CPsTG4JoKPkpTS0lAS1pmsKRlRtIv45fbmlpKgmkl1NRSuwEHY7t7\n3YB1gVuBY9o8SlNLSUC9mwbTu2nwyuuto+5K4mX8cltTS0kwjYSaWup8oD8wEDgC+3bsMUWfIdI+\nKLclWuWex6wj19JRKbclGuX8iNFkivW9Rdov5bZERd/8ExGJjAqziEhkVJhFRCKjwiwiEhkVZhGR\nyGhqKV9fB4w1pTlYqIYN1wsWK3XN6cFiAWx6yttB4mhqqSR9FSxSw8NjgsUCSF0T7udLGt4IeDbk\nNc3hYhWgHrOISGRUmEVEIqPCLCISGRVmEZHI+B78awU+Ab7BjhYMS6pBIjXUivJaIuRbmFPY7x8u\nTq4pIjWnvJYolTOUoal3pCNSXkt0fAtzCngCmAKckFxzRGpKeS1R8h3K2B2YB2wATATeAJ5eea+m\nlpKAZrYsZGbLolq8VPG8BjS1lITTSqippdLmub8fAfdiB0kyCayppSSgrZv6sHVTn5XXx496K6mX\nKp7XgKaWknAaCTW1FMDawDru/+7AfkDiM2WKJEx5LdHy6TFviPUm0o+/A3g8sRaJ1IbyWqLlU5hn\nA0OSbohIjSmvJVr65p+ISGRUmEVEIqPCLCISGRVmEZHIqDCLiEQmxO8EpGgJOG2LlKXbkHC/v/PF\nNeGmqQJI3R7mZyga3rA/QYKVJwUX1+FlxQwNFil1ziHBYjXMClTvHmyAAnmtHrOISGRUmEVEIqPC\nLCISGRVmEZHI+BTmXsB44HVgFrBLoi0SqR3ltkTJ57cyrgQeBg51j++eaItEake5LVEqVZh7AnsC\nx7rrXwNLE22RSG0otyVapYYyBmI/Ij4GmArciP2OrUh7p9yWaJXqMXfBzvI+FXgJuAI4D/hlm0dp\naikJqOVzaFmW+Mv45bamlpJQFrbAohavh5YqzHPd5SV3fTyWvG1paikJqKm7XdJGJTP9n19ua2op\nCaVPk13S3hpV8KGlhjLmA3OAzd31fYCZ1bRNJBLKbYmWz1kZ/4FNu9MVeAcYmWiLRGpHuS1R8inM\nrwA7Jd0QkTpQbkuU9M0/EZHIqDCLiERGhVlEJDIqzCIikVFhFhGJjM9ZGRKxL94ONx1UywU7B4sF\n0Hxh0HCy2pkaLFLD778MFuvxQLOc7VfkPvWYRUQio8IsIhIZFWYRkcioMIuIRManMG8BTMu6LAVO\nS7JRIjWgvJZo+ZyV8Sawvfu/E/A+cG9iLRKpDeW1RKvcoYx9sF/hmpNAW0TqRXktUSm3MB8B3JlE\nQ0TqSHktUSnnCyZdgR8A565yj6aWkoBa3aVGCuc1oKmlJJRX3MVHOYV5f+BlbALLtjS1lATUSNvy\nNznZlyuc14CmlpJQtnOXtNuLPLacoYwjgbsqapFIvJTXEh3fwtwdO0AyIcG2iNSa8lqi5DuU8TnQ\nJ8mGiNSB8lqipG/+iYhEpnaFeVqLYtUr1pRwsaa1fBIsVmuwSPXWqlgdIla4w8y+Z18UUrvCPL1F\nseoV6+VwsaarMOfRqlgdItbqWJhFRMSLCrOISGRCzJHSAgwPEEekkMnU55seLSi3JTn1ymsRERER\nERERkfZoBPAG8BYFf8XLy83AAmBGgDb1ByYBM4HXqG72im7AC8B0YBYwusq2dcZm1Xigyjhg5xO9\n6uK9WGWsXsB44HVsOXepME5HmT1EeV2+ULndivK6Kp2Bt7EfDFsD+5AHVRhrT2zWiRAJvBEwxP3f\nA5vRotJ2Aazt/nYBngf2qCLWGcAdwP1VxEibDawXIA7AWOAn7v8uQM8AMTsB87CC0p4orysTKrc7\ndF7X4nS5YVgCtwJfAeOAQyqM9TSwJEyzmI+tTACfYVvLTaqIt8z97YqttIsrjNMPOAC4iTBnzRAo\nTk+sgNzsrn+N9Qiq1V5nD1Fely90bnfYvK5FYe5L28bNdbfFpBHrsbxQRYxO2AqxANuVnFVhnMuB\ns4EVVbQlWwp4ApgCnFBFnIHYbxaPAaYCN5LpTVWjvc4eorwuX8jc7tB5XYvCnKrBa1SjBza+dDrW\nw6jUCmwXsh+wF5Wdn3gQ8CE2PhWqt7w7tnLuD5yC9Q4q0QUYClzn/n4OnFdl29Kzh9xdZZx6UF6X\nJ3Rud+i8rkVhfp+24yz9sd5FDNYA7sEmE7gvUMylwEPAjhU8dzfgYGz87C7ge8CtVbZnnvv7ETYL\n9LAK48x1l5fc9fFYIlejxOwhUVNelyd0biuvq9QFG2tpxLYk1RwkwcUJcZCkAUuMywPE6oMd2QVY\nC3gK2LvKmMOp/sj12sA67v/uwN+A/aqI9xSwufu/Gbi0ilhg47LHVhmjXpTXlas2t5XXgeyPHR1+\nG/hFFXHuAj4A/oGN742sItYe2G7adDKnt4yoMNZgbHxqOnYKz9lVtCttONUfuR6ItWk6dupUNe89\n2JRlL2E/njWB6o5edwcWklnB2iPldWWqzW3ltYiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiHRs/w9e\ned7lsJ91IAAAAABJRU5ErkJggg==\n", "text/plain": [ - "" + "" ] }, "metadata": {}, @@ -1936,18 +1969,18 @@ "\tDomain ID =\t10000\n", "\tNuclide =\tU-238\n", "\tCross Sections [cm^-1]:\n", - " Group 1 [6.25e-07 - 20.0 MeV]:\t2.16e-01 +/- 3.77e-01%\n", - " Group 2 [0.0 - 6.25e-07 MeV]:\t2.54e-01 +/- 6.46e-01%\n", + " Group 1 [6.25e-07 - 20.0 MeV]:\t2.17e-01 +/- 4.04e-01%\n", + " Group 2 [0.0 - 6.25e-07 MeV]:\t2.53e-01 +/- 5.85e-01%\n", "\n", "\tNuclide =\tO-16\n", "\tCross Sections [cm^-1]:\n", - " Group 1 [6.25e-07 - 20.0 MeV]:\t1.45e-01 +/- 4.03e-01%\n", - " Group 2 [0.0 - 6.25e-07 MeV]:\t1.75e-01 +/- 7.83e-01%\n", + " Group 1 [6.25e-07 - 20.0 MeV]:\t1.45e-01 +/- 4.10e-01%\n", + " Group 2 [0.0 - 6.25e-07 MeV]:\t1.75e-01 +/- 6.46e-01%\n", "\n", "\tNuclide =\tU-235\n", "\tCross Sections [cm^-1]:\n", - " Group 1 [6.25e-07 - 20.0 MeV]:\t7.72e-03 +/- 1.13e+00%\n", - " Group 2 [0.0 - 6.25e-07 MeV]:\t1.82e-01 +/- 5.18e-01%\n", + " Group 1 [6.25e-07 - 20.0 MeV]:\t7.91e-03 +/- 1.22e+00%\n", + " Group 2 [0.0 - 6.25e-07 MeV]:\t1.82e-01 +/- 4.98e-01%\n", "\n", "\n", "\n" @@ -1986,48 +2019,48 @@ " 10000\n", " 1\n", " U-238\n", - " 9.566947\n", - " 0.036112\n", + " 9.589323\n", + " 0.038756\n", " \n", " \n", " 4\n", " 10000\n", " 1\n", " O-16\n", - " 3.146780\n", - " 0.012666\n", + " 3.159101\n", + " 0.012939\n", " \n", " \n", " 5\n", " 10000\n", " 1\n", " U-235\n", - " 20.591253\n", - " 0.232675\n", + " 21.095256\n", + " 0.257787\n", " \n", " \n", " 0\n", " 10000\n", " 2\n", " U-238\n", - " 11.204912\n", - " 0.072348\n", + " 11.178844\n", + " 0.065428\n", " \n", " \n", " 1\n", " 10000\n", " 2\n", " O-16\n", - " 3.798407\n", - " 0.029742\n", + " 3.800027\n", + " 0.024538\n", " \n", " \n", " 2\n", " 10000\n", " 2\n", " U-235\n", - " 484.529684\n", - " 2.510940\n", + " 485.513530\n", + " 2.418761\n", " \n", " \n", "\n", @@ -2035,12 +2068,12 @@ ], "text/plain": [ " cell group in nuclide mean std. dev.\n", - "3 10000 1 U-238 9.566947 0.036112\n", - "4 10000 1 O-16 3.146780 0.012666\n", - "5 10000 1 U-235 20.591253 0.232675\n", - "0 10000 2 U-238 11.204912 0.072348\n", - "1 10000 2 O-16 3.798407 0.029742\n", - "2 10000 2 U-235 484.529684 2.510940" + "3 10000 1 U-238 9.589323 0.038756\n", + "4 10000 1 O-16 3.159101 0.012939\n", + "5 10000 1 U-235 21.095256 0.257787\n", + "0 10000 2 U-238 11.178844 0.065428\n", + "1 10000 2 O-16 3.800027 0.024538\n", + "2 10000 2 U-235 485.513530 2.418761" ] }, "execution_count": 46, @@ -2176,9 +2209,9 @@ "name": "stdout", "output_type": "stream", "text": [ - "openmc keff = 1.231090\n", - "openmoc keff = 1.229390\n", - "bias [pcm]: -170.1\n" + "openmc keff = 1.227616\n", + "openmoc keff = 1.225325\n", + "bias [pcm]: -229.1\n" ] } ], @@ -2206,7 +2239,23 @@ "metadata": { "collapsed": false }, - "outputs": [], + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/mgxs/mgxs.py:633: DeprecationWarning: elementwise comparison failed; this will raise the error in the future.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/mgxs/mgxs.py:633: DeprecationWarning: elementwise comparison failed; this will raise the error in the future.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/mgxs/mgxs.py:633: DeprecationWarning: elementwise comparison failed; this will raise the error in the future.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/mgxs/mgxs.py:633: DeprecationWarning: elementwise comparison failed; this will raise the error in the future.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/mgxs/mgxs.py:633: DeprecationWarning: elementwise comparison failed; this will raise the error in the future.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/mgxs/mgxs.py:633: DeprecationWarning: elementwise comparison failed; this will raise the error in the future.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/mgxs/mgxs.py:633: DeprecationWarning: elementwise comparison failed; this will raise the error in the future.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/mgxs/mgxs.py:633: DeprecationWarning: elementwise comparison failed; this will raise the error in the future.\n", + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.0-py2.7.egg/openmc/mgxs/mgxs.py:633: DeprecationWarning: elementwise comparison failed; this will raise the error in the future.\n" + ] + } + ], "source": [ "su.make_opencg_geometry()\n", "openmoc_geometry = get_openmoc_geometry(su.opencg_geometry)\n", @@ -2270,9 +2319,9 @@ "name": "stdout", "output_type": "stream", "text": [ - "openmc keff = 1.231090\n", - "openmoc keff = 1.232146\n", - "bias [pcm]: 105.6\n" + "openmc keff = 1.227616\n", + "openmoc keff = 1.227096\n", + "bias [pcm]: -52.0\n" ] } ], diff --git a/openmc/cross.py b/openmc/cross.py index 9b8a1d240..b91344671 100644 --- a/openmc/cross.py +++ b/openmc/cross.py @@ -9,7 +9,7 @@ if sys.version_info[0] >= 3: basestring = str # Acceptable tally arithmetic binary operations -TALLY_ARITHMETIC_OPS = ['+', '-', '*', '/', '^'] +_TALLY_ARITHMETIC_OPS = ['+', '-', '*', '/', '^'] class CrossScore(object): @@ -78,6 +78,11 @@ class CrossScore(object): else: return existing + def __repr__(self): + string = '({0} {1} {2})'.format(self.left_score, + self.binary_op, self.right_score) + return string + @property def left_score(self): return self._left_score @@ -103,14 +108,9 @@ class CrossScore(object): @binary_op.setter def binary_op(self, binary_op): cv.check_type('binary_op', binary_op, (basestring, CrossScore)) - cv.check_value('binary_op', binary_op, TALLY_ARITHMETIC_OPS) + cv.check_value('binary_op', binary_op, _TALLY_ARITHMETIC_OPS) self._binary_op = binary_op - def __repr__(self): - string = '({0} {1} {2})'.format(self.left_score, - self.binary_op, self.right_score) - return string - class CrossNuclide(object): """A special-purpose nuclide used to encapsulate all combinations of two @@ -178,6 +178,28 @@ class CrossNuclide(object): else: return existing + def __repr__(self): + + string = '' + + # If the Summary was linked, the left nuclide is a Nuclide object + if isinstance(self.left_nuclide, Nuclide): + string += '(' + self.left_nuclide.name + # If the Summary was not linked, the left nuclide is the ZAID + else: + string += '(' + str(self.left_nuclide) + + string += ' ' + self.binary_op + ' ' + + # If the Summary was linked, the right nuclide is a Nuclide object + if isinstance(self.right_nuclide, Nuclide): + string += self.right_nuclide.name + ')' + # If the Summary was not linked, the right nuclide is the ZAID + else: + string += str(self.right_nuclide) + ')' + + return string + @property def left_nuclide(self): return self._left_nuclide @@ -203,34 +225,9 @@ class CrossNuclide(object): @binary_op.setter def binary_op(self, binary_op): cv.check_type('binary_op', binary_op, basestring) - cv.check_value('binary_op', binary_op, TALLY_ARITHMETIC_OPS) + cv.check_value('binary_op', binary_op, _TALLY_ARITHMETIC_OPS) self._binary_op = binary_op - def __eq__(self, other): - return str(other) == str(self) - - def __repr__(self): - - string = '' - - # If the Summary was linked, the left nuclide is a Nuclide object - if isinstance(self.left_nuclide, Nuclide): - string += '(' + self.left_nuclide.name - # If the Summary was not linked, the left nuclide is the ZAID - else: - string += '(' + str(self.left_nuclide) - - string += ' ' + self.binary_op + ' ' - - # If the Summary was linked, the right nuclide is a Nuclide object - if isinstance(self.right_nuclide, Nuclide): - string += self.right_nuclide.name + ')' - # If the Summary was not linked, the right nuclide is the ZAID - else: - string += str(self.right_nuclide) + ')' - - return string - class CrossFilter(object): """A special-purpose filter used to encapsulate all combinations of two @@ -289,6 +286,18 @@ class CrossFilter(object): def __ne__(self, other): return not self == other + def __repr__(self): + + string = 'CrossFilter\n' + filter_type = '({0} {1} {2})'.format(self.left_filter.type, + self.binary_op, + self.right_filter.type) + filter_bins = '({0} {1} {2})'.format(self.left_filter.bins, + self.binary_op, + self.right_filter.bins) + string += '{0: <16}{1}{2}\n'.format('\tType', '=\t', filter_type) + string += '{0: <16}{1}{2}\n'.format('\tBins', '=\t', filter_bins) + return string def __deepcopy__(self, memo): existing = memo.get(id(self)) @@ -366,7 +375,7 @@ class CrossFilter(object): @binary_op.setter def binary_op(self, binary_op): cv.check_type('binary_op', binary_op, basestring) - cv.check_value('binary_op', binary_op, TALLY_ARITHMETIC_OPS) + cv.check_value('binary_op', binary_op, _TALLY_ARITHMETIC_OPS) self._binary_op = binary_op @stride.setter @@ -415,7 +424,7 @@ class CrossFilter(object): Parameters ---------- - data_size : Integral + datasize : Integral The total number of bins in the tally corresponding to this filter summary : None or Summary An optional Summary object to be used to construct columns for @@ -452,17 +461,4 @@ class CrossFilter(object): right_df = right_df.astype(str) df = '(' + left_df + ' ' + self.binary_op + ' ' + right_df + ')' - return df - - def __repr__(self): - - string = 'CrossFilter\n' - filter_type = '({0} {1} {2})'.format(self.left_filter.type, - self.binary_op, - self.right_filter.type) - filter_bins = '({0} {1} {2})'.format(self.left_filter.bins, - self.binary_op, - self.right_filter.bins) - string += '{0: <16}{1}{2}\n'.format('\tType', '=\t', filter_type) - string += '{0: <16}{1}{2}\n'.format('\tBins', '=\t', filter_bins) - return string \ No newline at end of file + return df \ No newline at end of file diff --git a/openmc/element.py b/openmc/element.py index a99d47127..56821b5d2 100644 --- a/openmc/element.py +++ b/openmc/element.py @@ -51,9 +51,17 @@ class Element(object): else: return False + def __ne__(self, other): + return not self == other + def __hash__(self): return hash((self._name, self._xs)) + def __repr__(self): + string = 'Element - {0}\n'.format(self._name) + string += '{0: <16}{1}{2}\n'.format('\tXS', '=\t', self._xs) + return string + @property def xs(self): return self._xs @@ -70,9 +78,4 @@ class Element(object): @name.setter def name(self, name): check_type('name', name, basestring) - self._name = name - - def __repr__(self): - string = 'Element - {0}\n'.format(self._name) - string += '{0: <16}{1}{2}\n'.format('\tXS', '=\t', self._xs) - return string + self._name = name \ No newline at end of file diff --git a/openmc/filter.py b/openmc/filter.py index 385ba20cf..63034fe57 100644 --- a/openmc/filter.py +++ b/openmc/filter.py @@ -104,6 +104,13 @@ class Filter(object): else: return existing + def __repr__(self): + string = 'Filter\n' + string += '{0: <16}{1}{2}\n'.format('\tType', '=\t', self.type) + string += '{0: <16}{1}{2}\n'.format('\tBins', '=\t', self.bins) + string += '{0: <16}{1}{2}\n'.format('\tOffset', '=\t', self.offset) + return string + @property def type(self): return self._type @@ -741,11 +748,4 @@ class Filter(object): filter_bins = filter_bins df = pd.concat([df, pd.DataFrame({self.type : filter_bins})]) - return df - - def __repr__(self): - string = 'Filter\n' - string += '{0: <16}{1}{2}\n'.format('\tType', '=\t', self.type) - string += '{0: <16}{1}{2}\n'.format('\tBins', '=\t', self.bins) - string += '{0: <16}{1}{2}\n'.format('\tOffset', '=\t', self.offset) - return string + return df \ No newline at end of file diff --git a/openmc/material.py b/openmc/material.py index 3ca8f2413..92c77858d 100644 --- a/openmc/material.py +++ b/openmc/material.py @@ -83,6 +83,38 @@ class Material(object): # If specified, this file will be used instead of composition values self._distrib_otf_file = None + def __repr__(self): + string = 'Material\n' + string += '{0: <16}{1}{2}\n'.format('\tID', '=\t', self._id) + string += '{0: <16}{1}{2}\n'.format('\tName', '=\t', self._name) + + string += '{0: <16}{1}{2}'.format('\tDensity', '=\t', self._density) + string += ' [{0}]\n'.format(self._density_units) + + string += '{0: <16}\n'.format('\tS(a,b) Tables') + + for sab in self._sab: + string += '{0: <16}{1}[{2}{3}]\n'.format('\tS(a,b)', '=\t', + sab[0], sab[1]) + + string += '{0: <16}\n'.format('\tNuclides') + + for nuclide in self._nuclides: + percent = self._nuclides[nuclide][1] + percent_type = self._nuclides[nuclide][2] + string += '{0: <16}'.format('\t{0}'.format(nuclide)) + string += '=\t{0: <12} [{1}]\n'.format(percent, percent_type) + + string += '{0: <16}\n'.format('\tElements') + + for element in self._elements: + percent = self._nuclides[element][1] + percent_type = self._nuclides[element][2] + string += '{0: >16}'.format('\t{0}'.format(element)) + string += '=\t{0: <12} [{1}]\n'.format(percent, percent_type) + + return string + @property def id(self): return self._id @@ -335,38 +367,6 @@ class Material(object): return nuclides - def __repr__(self): - string = 'Material\n' - string += '{0: <16}{1}{2}\n'.format('\tID', '=\t', self._id) - string += '{0: <16}{1}{2}\n'.format('\tName', '=\t', self._name) - - string += '{0: <16}{1}{2}'.format('\tDensity', '=\t', self._density) - string += ' [{0}]\n'.format(self._density_units) - - string += '{0: <16}\n'.format('\tS(a,b) Tables') - - for sab in self._sab: - string += '{0: <16}{1}[{2}{3}]\n'.format('\tS(a,b)', '=\t', - sab[0], sab[1]) - - string += '{0: <16}\n'.format('\tNuclides') - - for nuclide in self._nuclides: - percent = self._nuclides[nuclide][1] - percent_type = self._nuclides[nuclide][2] - string += '{0: <16}'.format('\t{0}'.format(nuclide)) - string += '=\t{0: <12} [{1}]\n'.format(percent, percent_type) - - string += '{0: <16}\n'.format('\tElements') - - for element in self._elements: - percent = self._nuclides[element][1] - percent_type = self._nuclides[element][2] - string += '{0: >16}'.format('\t{0}'.format(element)) - string += '=\t{0: <12} [{1}]\n'.format(percent, percent_type) - - return string - def _get_nuclide_xml(self, nuclide, distrib=False): xml_element = ET.Element("nuclide") xml_element.set("name", nuclide[0]._name) diff --git a/openmc/nuclide.py b/openmc/nuclide.py index a616edac9..f24e62a48 100644 --- a/openmc/nuclide.py +++ b/openmc/nuclide.py @@ -54,9 +54,19 @@ class Nuclide(object): else: return False + def __ne__(self, other): + return not self == other + def __hash__(self): return hash((self._name, self._xs)) + def __repr__(self): + string = 'Nuclide - {0}\n'.format(self._name) + string += '{0: <16}{1}{2}\n'.format('\tXS', '=\t', self._xs) + if self._zaid is not None: + string += '{0: <16}{1}{2}\n'.format('\tZAID', '=\t', self._zaid) + return string + @property def name(self): return self._name @@ -82,11 +92,4 @@ class Nuclide(object): @zaid.setter def zaid(self, zaid): check_type('zaid', zaid, Integral) - self._zaid = zaid - - def __repr__(self): - string = 'Nuclide - {0}\n'.format(self._name) - string += '{0: <16}{1}{2}\n'.format('\tXS', '=\t', self._xs) - if self._zaid is not None: - string += '{0: <16}{1}{2}\n'.format('\tZAID', '=\t', self._zaid) - return string + self._zaid = zaid \ No newline at end of file diff --git a/openmc/region.py b/openmc/region.py index 5baac22dd..f8f3a7378 100644 --- a/openmc/region.py +++ b/openmc/region.py @@ -213,6 +213,9 @@ class Intersection(Region): def __init__(self, *nodes): self.nodes = list(nodes) + def __str__(self): + return '(' + ' '.join(map(str, self.nodes)) + ')' + @property def nodes(self): return self._nodes @@ -222,9 +225,6 @@ class Intersection(Region): check_type('nodes', nodes, Iterable, Region) self._nodes = nodes - def __str__(self): - return '(' + ' '.join(map(str, self.nodes)) + ')' - class Union(Region): """Union of two or more regions. @@ -252,6 +252,9 @@ class Union(Region): def __init__(self, *nodes): self.nodes = list(nodes) + def __str__(self): + return '(' + ' | '.join(map(str, self.nodes)) + ')' + @property def nodes(self): return self._nodes @@ -261,9 +264,6 @@ class Union(Region): check_type('nodes', nodes, Iterable, Region) self._nodes = nodes - def __str__(self): - return '(' + ' | '.join(map(str, self.nodes)) + ')' - class Complement(Region): """Complement of a region. @@ -295,6 +295,9 @@ class Complement(Region): def __init__(self, node): self.node = node + def __str__(self): + return '~' + str(self.node) + @property def node(self): return self._node @@ -303,6 +306,3 @@ class Complement(Region): def node(self, node): check_type('node', node, Region) self._node = node - - def __str__(self): - return '~' + str(self.node) diff --git a/openmc/surface.py b/openmc/surface.py index afa426fa3..ae2e222ac 100644 --- a/openmc/surface.py +++ b/openmc/surface.py @@ -75,6 +75,22 @@ class Surface(object): def __pos__(self): return Halfspace(self, '+') + def __repr__(self): + string = 'Surface\n' + string += '{0: <16}{1}{2}\n'.format('\tID', '=\t', self._id) + string += '{0: <16}{1}{2}\n'.format('\tName', '=\t', self._name) + string += '{0: <16}{1}{2}\n'.format('\tType', '=\t', self._type) + string += '{0: <16}{1}{2}\n'.format('\tBoundary', '=\t', self._boundary_type) + + coeffs = '{0: <16}'.format('\tCoefficients') + '\n' + + for coeff in self._coeffs: + coeffs += '{0: <16}{1}{2}\n'.format(coeff, '=\t', self._coeffs[coeff]) + + string += coeffs + + return string + @property def id(self): return self._id @@ -120,22 +136,6 @@ class Surface(object): check_value('boundary type', boundary_type, _BC_TYPES) self._boundary_type = boundary_type - def __repr__(self): - string = 'Surface\n' - string += '{0: <16}{1}{2}\n'.format('\tID', '=\t', self._id) - string += '{0: <16}{1}{2}\n'.format('\tName', '=\t', self._name) - string += '{0: <16}{1}{2}\n'.format('\tType', '=\t', self._type) - string += '{0: <16}{1}{2}\n'.format('\tBoundary', '=\t', self._boundary_type) - - coeffs = '{0: <16}'.format('\tCoefficients') + '\n' - - for coeff in self._coeffs: - coeffs += '{0: <16}{1}{2}\n'.format(coeff, '=\t', self._coeffs[coeff]) - - string += coeffs - - return string - def create_xml_subelement(self): element = ET.Element("surface") element.set("id", str(self._id)) diff --git a/openmc/tallies.py b/openmc/tallies.py index a25004ff9..f990572e5 100644 --- a/openmc/tallies.py +++ b/openmc/tallies.py @@ -204,6 +204,32 @@ class Tally(object): return hash(tuple(hashable)) + def __repr__(self): + string = 'Tally\n' + string += '{0: <16}{1}{2}\n'.format('\tID', '=\t', self.id) + string += '{0: <16}{1}{2}\n'.format('\tName', '=\t', self.name) + + string += '{0: <16}{1}\n'.format('\tFilters', '=\t') + + for filter in self.filters: + string += '{0: <16}\t\t{1}\t{2}\n'.format('', filter.type, + filter.bins) + + string += '{0: <16}{1}'.format('\tNuclides', '=\t') + + for nuclide in self.nuclides: + if isinstance(nuclide, Nuclide): + string += '{0} '.format(nuclide.name) + else: + string += '{0} '.format(nuclide) + + string += '\n' + + string += '{0: <16}{1}{2}\n'.format('\tScores', '=\t', self.scores) + string += '{0: <16}{1}{2}\n'.format('\tEstimator', '=\t', self.estimator) + + return string + @property def id(self): return self._id @@ -518,32 +544,6 @@ class Tally(object): self._nuclides.remove(nuclide) - def __repr__(self): - string = 'Tally\n' - string += '{0: <16}{1}{2}\n'.format('\tID', '=\t', self.id) - string += '{0: <16}{1}{2}\n'.format('\tName', '=\t', self.name) - - string += '{0: <16}{1}\n'.format('\tFilters', '=\t') - - for filter in self.filters: - string += '{0: <16}\t\t{1}\t{2}\n'.format('', filter.type, - filter.bins) - - string += '{0: <16}{1}'.format('\tNuclides', '=\t') - - for nuclide in self.nuclides: - if isinstance(nuclide, Nuclide): - string += '{0} '.format(nuclide.name) - else: - string += '{0} '.format(nuclide) - - string += '\n' - - string += '{0: <16}{1}{2}\n'.format('\tScores', '=\t', self.scores) - string += '{0: <16}{1}{2}\n'.format('\tEstimator', '=\t', self.estimator) - - return string - def can_merge(self, tally): """Determine if another tally can be merged with this one diff --git a/openmc/trigger.py b/openmc/trigger.py index e695defde..5af477edd 100644 --- a/openmc/trigger.py +++ b/openmc/trigger.py @@ -58,6 +58,13 @@ class Trigger(object): else: return existing + def __repr__(self): + string = 'Trigger\n' + string += '{0: <16}{1}{2}\n'.format('\tType', '=\t', self._trigger_type) + string += '{0: <16}{1}{2}\n'.format('\tThreshold', '=\t', self._threshold) + string += '{0: <16}{1}{2}\n'.format('\tScores', '=\t', self._scores) + return string + @property def trigger_type(self): return self._trigger_type @@ -102,13 +109,6 @@ class Trigger(object): else: self._scores.append(score) - def __repr__(self): - string = 'Trigger\n' - string += '{0: <16}{1}{2}\n'.format('\tType', '=\t', self._trigger_type) - string += '{0: <16}{1}{2}\n'.format('\tThreshold', '=\t', self._threshold) - string += '{0: <16}{1}{2}\n'.format('\tScores', '=\t', self._scores) - return string - def get_trigger_xml(self, element): """Return XML representation of the trigger diff --git a/openmc/universe.py b/openmc/universe.py index e9f7f3284..a4ea38de9 100644 --- a/openmc/universe.py +++ b/openmc/universe.py @@ -73,6 +73,30 @@ class Cell(object): self._translation = None self._offsets = None + def __repr__(self): + string = 'Cell\n' + string += '{0: <16}{1}{2}\n'.format('\tID', '=\t', self._id) + string += '{0: <16}{1}{2}\n'.format('\tName', '=\t', self._name) + + if isinstance(self._fill, openmc.Material): + string += '{0: <16}{1}{2}\n'.format('\tMaterial', '=\t', + self._fill._id) + elif isinstance(self._fill, (Universe, Lattice)): + string += '{0: <16}{1}{2}\n'.format('\tFill', '=\t', + self._fill._id) + else: + string += '{0: <16}{1}{2}\n'.format('\tFill', '=\t', self._fill) + + string += '{0: <16}{1}{2}\n'.format('\tRegion', '=\t', self._region) + + string += '{0: <16}{1}{2}\n'.format('\tRotation', '=\t', + self._rotation) + string += '{0: <16}{1}{2}\n'.format('\tTranslation', '=\t', + self._translation) + string += '{0: <16}{1}{2}\n'.format('\tOffset', '=\t', self._offsets) + + return string + @property def id(self): return self._id @@ -298,30 +322,6 @@ class Cell(object): return universes - def __repr__(self): - string = 'Cell\n' - string += '{0: <16}{1}{2}\n'.format('\tID', '=\t', self._id) - string += '{0: <16}{1}{2}\n'.format('\tName', '=\t', self._name) - - if isinstance(self._fill, openmc.Material): - string += '{0: <16}{1}{2}\n'.format('\tMaterial', '=\t', - self._fill._id) - elif isinstance(self._fill, (Universe, Lattice)): - string += '{0: <16}{1}{2}\n'.format('\tFill', '=\t', - self._fill._id) - else: - string += '{0: <16}{1}{2}\n'.format('\tFill', '=\t', self._fill) - - string += '{0: <16}{1}{2}\n'.format('\tRegion', '=\t', self._region) - - string += '{0: <16}{1}{2}\n'.format('\tRotation', '=\t', - self._rotation) - string += '{0: <16}{1}{2}\n'.format('\tTranslation', '=\t', - self._translation) - string += '{0: <16}{1}{2}\n'.format('\tOffset', '=\t', self._offsets) - - return string - def create_xml_subelement(self, xml_element): element = ET.Element("cell") element.set("id", str(self._id)) @@ -836,6 +836,50 @@ class RectLattice(Lattice): self._lower_left = None self._offsets = None + def __repr__(self): + string = 'RectLattice\n' + string += '{0: <16}{1}{2}\n'.format('\tID', '=\t', self._id) + string += '{0: <16}{1}{2}\n'.format('\tName', '=\t', self._name) + string += '{0: <16}{1}{2}\n'.format('\tDimension', '=\t', + self._dimension) + string += '{0: <16}{1}{2}\n'.format('\tLower Left', '=\t', + self._lower_left) + string += '{0: <16}{1}{2}\n'.format('\tPitch', '=\t', self._pitch) + + if self._outer is not None: + string += '{0: <16}{1}{2}\n'.format('\tOuter', '=\t', + self._outer._id) + else: + string += '{0: <16}{1}{2}\n'.format('\tOuter', '=\t', + self._outer) + + string += '{0: <16}\n'.format('\tUniverses') + + # Lattice nested Universe IDs - column major for Fortran + for i, universe in enumerate(np.ravel(self._universes)): + string += '{0} '.format(universe._id) + + # Add a newline character every time we reach end of row of cells + if (i+1) % self._dimension[-1] == 0: + string += '\n' + + string = string.rstrip('\n') + + if self._offsets is not None: + string += '{0: <16}\n'.format('\tOffsets') + + # Lattice cell offsets + for i, offset in enumerate(np.ravel(self._offsets)): + string += '{0} '.format(offset) + + # Add a newline character when we reach end of row of cells + if (i+1) % self._dimension[-1] == 0: + string += '\n' + + string = string.rstrip('\n') + + return string + @property def dimension(self): return self._dimension @@ -893,50 +937,6 @@ class RectLattice(Lattice): return offset - def __repr__(self): - string = 'RectLattice\n' - string += '{0: <16}{1}{2}\n'.format('\tID', '=\t', self._id) - string += '{0: <16}{1}{2}\n'.format('\tName', '=\t', self._name) - string += '{0: <16}{1}{2}\n'.format('\tDimension', '=\t', - self._dimension) - string += '{0: <16}{1}{2}\n'.format('\tLower Left', '=\t', - self._lower_left) - string += '{0: <16}{1}{2}\n'.format('\tPitch', '=\t', self._pitch) - - if self._outer is not None: - string += '{0: <16}{1}{2}\n'.format('\tOuter', '=\t', - self._outer._id) - else: - string += '{0: <16}{1}{2}\n'.format('\tOuter', '=\t', - self._outer) - - string += '{0: <16}\n'.format('\tUniverses') - - # Lattice nested Universe IDs - column major for Fortran - for i, universe in enumerate(np.ravel(self._universes)): - string += '{0} '.format(universe._id) - - # Add a newline character every time we reach end of row of cells - if (i+1) % self._dimension[-1] == 0: - string += '\n' - - string = string.rstrip('\n') - - if self._offsets is not None: - string += '{0: <16}\n'.format('\tOffsets') - - # Lattice cell offsets - for i, offset in enumerate(np.ravel(self._offsets)): - string += '{0} '.format(offset) - - # Add a newline character when we reach end of row of cells - if (i+1) % self._dimension[-1] == 0: - string += '\n' - - string = string.rstrip('\n') - - return string - def create_xml_subelement(self, xml_element): # Determine if XML element already contains subelement for this Lattice path = './lattice[@id=\'{0}\']'.format(self._id) @@ -1052,6 +1052,34 @@ class HexLattice(Lattice): self._num_axial = None self._center = None + def __repr__(self): + string = 'HexLattice\n' + string += '{0: <16}{1}{2}\n'.format('\tID', '=\t', self._id) + string += '{0: <16}{1}{2}\n'.format('\tName', '=\t', self._name) + string += '{0: <16}{1}{2}\n'.format('\t# Rings', '=\t', self._num_rings) + string += '{0: <16}{1}{2}\n'.format('\t# Axial', '=\t', self._num_axial) + string += '{0: <16}{1}{2}\n'.format('\tCenter', '=\t', + self._center) + string += '{0: <16}{1}{2}\n'.format('\tPitch', '=\t', self._pitch) + + if self._outer is not None: + string += '{0: <16}{1}{2}\n'.format('\tOuter', '=\t', + self._outer._id) + else: + string += '{0: <16}{1}{2}\n'.format('\tOuter', '=\t', + self._outer) + + string += '{0: <16}\n'.format('\tUniverses') + + if self._num_axial is not None: + slices = [self._repr_axial_slice(x) for x in self._universes] + string += '\n'.join(slices) + + else: + string += self._repr_axial_slice(self._universes) + + return string + @property def num_rings(self): return self._num_rings @@ -1172,34 +1200,6 @@ class HexLattice(Lattice): 6*(self._num_rings - 1 - r)) raise ValueError(msg) - def __repr__(self): - string = 'HexLattice\n' - string += '{0: <16}{1}{2}\n'.format('\tID', '=\t', self._id) - string += '{0: <16}{1}{2}\n'.format('\tName', '=\t', self._name) - string += '{0: <16}{1}{2}\n'.format('\t# Rings', '=\t', self._num_rings) - string += '{0: <16}{1}{2}\n'.format('\t# Axial', '=\t', self._num_axial) - string += '{0: <16}{1}{2}\n'.format('\tCenter', '=\t', - self._center) - string += '{0: <16}{1}{2}\n'.format('\tPitch', '=\t', self._pitch) - - if self._outer is not None: - string += '{0: <16}{1}{2}\n'.format('\tOuter', '=\t', - self._outer._id) - else: - string += '{0: <16}{1}{2}\n'.format('\tOuter', '=\t', - self._outer) - - string += '{0: <16}\n'.format('\tUniverses') - - if self._num_axial is not None: - slices = [self._repr_axial_slice(x) for x in self._universes] - string += '\n'.join(slices) - - else: - string += self._repr_axial_slice(self._universes) - - return string - def create_xml_subelement(self, xml_element): # Determine if XML element already contains subelement for this Lattice path = './hex_lattice[@id=\'{0}\']'.format(self._id)