diff --git a/docs/source/pythonapi/examples/mgxs-part-i.ipynb b/docs/source/pythonapi/examples/mgxs-part-i.ipynb index e4c976718..2f2a80177 100644 --- a/docs/source/pythonapi/examples/mgxs-part-i.ipynb +++ b/docs/source/pythonapi/examples/mgxs-part-i.ipynb @@ -372,6 +372,7 @@ "\n", "* `TotalXS`\n", "* `TransportXS`\n", + "* `NuTransportXS`\n", "* `AbsorptionXS`\n", "* `CaptureXS`\n", "* `FissionXS`\n", @@ -409,7 +410,7 @@ }, { "cell_type": "code", - "execution_count": 13, + "execution_count": 26, "metadata": { "collapsed": false }, @@ -418,25 +419,27 @@ "data": { "text/plain": [ "OrderedDict([('flux', Tally\n", - " \tID =\t10000\n", - " \tName =\t\n", - " \tFilters =\t\n", - " \t\tcell\t[1]\n", - " \t\tenergy\t[ 0.00000000e+00 6.25000000e-07 2.00000000e+01]\n", - " \tNuclides =\ttotal \n", - " \tScores =\t['flux']\n", - " \tEstimator =\ttracklength), ('absorption', Tally\n", - " \tID =\t10001\n", - " \tName =\t\n", - " \tFilters =\t\n", - " \t\tcell\t[1]\n", - " \t\tenergy\t[ 0.00000000e+00 6.25000000e-07 2.00000000e+01]\n", - " \tNuclides =\ttotal \n", - " \tScores =\t['absorption']\n", - " \tEstimator =\ttracklength)])" + "\tID =\t10012\n", + "\tName =\t\n", + "\tFilters =\t\n", + " \t\tcell\t[1]\n", + " \t\tenergy\t[ 0.00000000e+00 6.25000000e-07 2.00000000e+01]\n", + "\tNuclides =\ttotal \n", + "\tScores =\t[u'flux']\n", + "\tEstimator =\ttracklength\n", + "), ('absorption', Tally\n", + "\tID =\t10013\n", + "\tName =\t\n", + "\tFilters =\t\n", + " \t\tcell\t[1]\n", + " \t\tenergy\t[ 0.00000000e+00 6.25000000e-07 2.00000000e+01]\n", + "\tNuclides =\ttotal \n", + "\tScores =\t[u'absorption']\n", + "\tEstimator =\ttracklength\n", + ")])" ] }, - "execution_count": 13, + "execution_count": 26, "metadata": {}, "output_type": "execute_result" } @@ -510,8 +513,8 @@ " Copyright: 2011-2016 Massachusetts Institute of Technology\n", " License: http://openmc.readthedocs.io/en/latest/license.html\n", " Version: 0.7.1\n", - " Git SHA1: 502482dcf630ee6e290c15b8535e6e850a351c88\n", - " Date/Time: 2016-05-10 20:52:19\n", + " Git SHA1: ae588276014a905ecc6e0967bf08288ecec5b550\n", + " Date/Time: 2016-05-12 20:41:27\n", " MPI Processes: 1\n", "\n", " ===========================================================================\n", @@ -597,20 +600,20 @@ "\n", " =======================> TIMING STATISTICS <=======================\n", "\n", - " Total time for initialization = 5.3200E-01 seconds\n", - " Reading cross sections = 1.3300E-01 seconds\n", - " Total time in simulation = 2.3438E+01 seconds\n", - " Time in transport only = 2.3419E+01 seconds\n", - " Time in inactive batches = 2.9490E+00 seconds\n", - " Time in active batches = 2.0489E+01 seconds\n", - " Time synchronizing fission bank = 6.0000E-03 seconds\n", - " Sampling source sites = 3.0000E-03 seconds\n", - " SEND/RECV source sites = 1.0000E-03 seconds\n", - " Time accumulating tallies = 0.0000E+00 seconds\n", + " Total time for initialization = 4.7500E-01 seconds\n", + " Reading cross sections = 9.7000E-02 seconds\n", + " Total time in simulation = 1.8074E+01 seconds\n", + " Time in transport only = 1.8055E+01 seconds\n", + " Time in inactive batches = 2.1180E+00 seconds\n", + " Time in active batches = 1.5956E+01 seconds\n", + " Time synchronizing fission bank = 4.0000E-03 seconds\n", + " Sampling source sites = 4.0000E-03 seconds\n", + " SEND/RECV source sites = 0.0000E+00 seconds\n", + " Time accumulating tallies = 2.0000E-03 seconds\n", " Total time for finalization = 0.0000E+00 seconds\n", - " Total time elapsed = 2.3985E+01 seconds\n", - " Calculation Rate (inactive) = 8477.45 neutrons/second\n", - " Calculation Rate (active) = 4880.67 neutrons/second\n", + " Total time elapsed = 1.8559E+01 seconds\n", + " Calculation Rate (inactive) = 11803.6 neutrons/second\n", + " Calculation Rate (active) = 6267.23 neutrons/second\n", "\n", " ============================> RESULTS <============================\n", "\n", @@ -1121,7 +1124,7 @@ " 6.250000e-07\n", " total\n", " (((absorption / flux) / (total / flux)) + ((sc...\n", - " 1\n", + " 1.0\n", " 0.007763\n", " \n", " \n", @@ -1131,7 +1134,7 @@ " 2.000000e+01\n", " total\n", " (((absorption / flux) / (total / flux)) + ((sc...\n", - " 1\n", + " 1.0\n", " 0.003739\n", " \n", " \n", @@ -1178,7 +1181,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", - "version": "2.7.11" + "version": "2.7.6" } }, "nbformat": 4, diff --git a/docs/source/pythonapi/examples/mgxs-part-ii.ipynb b/docs/source/pythonapi/examples/mgxs-part-ii.ipynb index 7b313da74..ca07519e5 100644 --- a/docs/source/pythonapi/examples/mgxs-part-ii.ipynb +++ b/docs/source/pythonapi/examples/mgxs-part-ii.ipynb @@ -445,8 +445,8 @@ " Copyright: 2011-2016 Massachusetts Institute of Technology\n", " License: http://openmc.readthedocs.io/en/latest/license.html\n", " Version: 0.7.1\n", - " Git SHA1: 7b20f8ad4aa9e6f02f8b1d51e002f9f56ba7aa15\n", - " Date/Time: 2016-05-09 13:34:05\n", + " Git SHA1: ae588276014a905ecc6e0967bf08288ecec5b550\n", + " Date/Time: 2016-05-12 21:00:03\n", " MPI Processes: 1\n", "\n", " ===========================================================================\n", @@ -523,7 +523,7 @@ " 48/1 1.21610 1.22612 +/- 0.00251\n", " 49/1 1.22199 1.22602 +/- 0.00245\n", " 50/1 1.20860 1.22558 +/- 0.00243\n", - " Triggers unsatisfied, max unc./thresh. is 1.25496 for flux in tally 10050\n", + " Triggers unsatisfied, max unc./thresh. is 1.25496 for flux in tally 10051\n", " The estimated number of batches is 73\n", " Creating state point statepoint.050.h5...\n", " 51/1 1.21850 1.22541 +/- 0.00237\n", @@ -549,7 +549,7 @@ " 71/1 1.19720 1.22444 +/- 0.00195\n", " 72/1 1.23770 1.22465 +/- 0.00193\n", " 73/1 1.23894 1.22488 +/- 0.00191\n", - " Triggers unsatisfied, max unc./thresh. is 1.00243 for flux in tally 10050\n", + " Triggers unsatisfied, max unc./thresh. is 1.00243 for flux in tally 10051\n", " The estimated number of batches is 74\n", " 74/1 1.22437 1.22487 +/- 0.00188\n", " Triggers satisfied for batch 74\n", @@ -562,20 +562,20 @@ "\n", " =======================> TIMING STATISTICS <=======================\n", "\n", - " Total time for initialization = 3.8900E-01 seconds\n", - " Reading cross sections = 8.3000E-02 seconds\n", - " Total time in simulation = 2.2066E+02 seconds\n", - " Time in transport only = 2.2061E+02 seconds\n", - " Time in inactive batches = 1.5872E+01 seconds\n", - " Time in active batches = 2.0478E+02 seconds\n", - " Time synchronizing fission bank = 1.9000E-02 seconds\n", - " Sampling source sites = 1.2000E-02 seconds\n", - " SEND/RECV source sites = 6.0000E-03 seconds\n", - " Time accumulating tallies = 3.0000E-03 seconds\n", - " Total time for finalization = 1.1000E-02 seconds\n", - " Total time elapsed = 2.2111E+02 seconds\n", - " Calculation Rate (inactive) = 6300.40 neutrons/second\n", - " Calculation Rate (active) = 1953.28 neutrons/second\n", + " Total time for initialization = 4.1000E-01 seconds\n", + " Reading cross sections = 8.6000E-02 seconds\n", + " Total time in simulation = 2.2903E+02 seconds\n", + " Time in transport only = 2.2897E+02 seconds\n", + " Time in inactive batches = 1.4619E+01 seconds\n", + " Time in active batches = 2.1441E+02 seconds\n", + " Time synchronizing fission bank = 2.5000E-02 seconds\n", + " Sampling source sites = 1.6000E-02 seconds\n", + " SEND/RECV source sites = 8.0000E-03 seconds\n", + " Time accumulating tallies = 1.0000E-03 seconds\n", + " Total time for finalization = 1.2000E-02 seconds\n", + " Total time elapsed = 2.2951E+02 seconds\n", + " Calculation Rate (inactive) = 6840.41 neutrons/second\n", + " Calculation Rate (active) = 1865.57 neutrons/second\n", "\n", " ============================> RESULTS <============================\n", "\n", @@ -786,9 +786,10 @@ " group in\n", " group out\n", " nuclide\n", - " score\n", + " moment\n", " mean\n", " std. dev.\n", + " moment\n", " \n", " \n", " \n", @@ -798,9 +799,10 @@ " 1\n", " 1\n", " H-1\n", - " ((nu-scatter-0 - scatter-1) / flux)\n", + " P0\n", " 0.234115\n", " 0.003568\n", + " P0\n", " \n", " \n", " 127\n", @@ -808,9 +810,10 @@ " 1\n", " 1\n", " O-16\n", - " ((nu-scatter-0 - scatter-1) / flux)\n", + " P0\n", " 1.563707\n", " 0.005953\n", + " P0\n", " \n", " \n", " 124\n", @@ -818,9 +821,10 @@ " 1\n", " 2\n", " H-1\n", - " ((nu-scatter-0 - scatter-1) / flux)\n", + " P0\n", " 1.594129\n", " 0.002369\n", + " P0\n", " \n", " \n", " 125\n", @@ -828,9 +832,10 @@ " 1\n", " 2\n", " O-16\n", - " ((nu-scatter-0 - scatter-1) / flux)\n", + " P0\n", " 0.285761\n", " 0.001676\n", + " P0\n", " \n", " \n", " 122\n", @@ -838,9 +843,10 @@ " 1\n", " 3\n", " H-1\n", - " ((nu-scatter-0 - scatter-1) / flux)\n", + " P0\n", " 0.011089\n", " 0.000248\n", + " P0\n", " \n", " \n", " 123\n", @@ -848,9 +854,10 @@ " 1\n", " 3\n", " O-16\n", - " ((nu-scatter-0 - scatter-1) / flux)\n", + " P0\n", " 0.000000\n", " 0.000000\n", + " P0\n", " \n", " \n", " 120\n", @@ -858,9 +865,10 @@ " 1\n", " 4\n", " H-1\n", - " ((nu-scatter-0 - scatter-1) / flux)\n", + " P0\n", " 0.000000\n", " 0.000000\n", + " P0\n", " \n", " \n", " 121\n", @@ -868,9 +876,10 @@ " 1\n", " 4\n", " O-16\n", - " ((nu-scatter-0 - scatter-1) / flux)\n", + " P0\n", " 0.000000\n", " 0.000000\n", + " P0\n", " \n", " \n", " 118\n", @@ -878,9 +887,10 @@ " 1\n", " 5\n", " H-1\n", - " ((nu-scatter-0 - scatter-1) / flux)\n", + " P0\n", " 0.000000\n", " 0.000000\n", + " P0\n", " \n", " \n", " 119\n", @@ -888,38 +898,27 @@ " 1\n", " 5\n", " O-16\n", - " ((nu-scatter-0 - scatter-1) / flux)\n", + " P0\n", " 0.000000\n", " 0.000000\n", + " P0\n", " \n", " \n", "\n", "" ], "text/plain": [ - " cell group in group out nuclide score \\\n", - "126 10002 1 1 H-1 ((nu-scatter-0 - scatter-1) / flux) \n", - "127 10002 1 1 O-16 ((nu-scatter-0 - scatter-1) / flux) \n", - "124 10002 1 2 H-1 ((nu-scatter-0 - scatter-1) / flux) \n", - "125 10002 1 2 O-16 ((nu-scatter-0 - scatter-1) / flux) \n", - "122 10002 1 3 H-1 ((nu-scatter-0 - scatter-1) / flux) \n", - "123 10002 1 3 O-16 ((nu-scatter-0 - scatter-1) / flux) \n", - "120 10002 1 4 H-1 ((nu-scatter-0 - scatter-1) / flux) \n", - "121 10002 1 4 O-16 ((nu-scatter-0 - scatter-1) / flux) \n", - "118 10002 1 5 H-1 ((nu-scatter-0 - scatter-1) / flux) \n", - "119 10002 1 5 O-16 ((nu-scatter-0 - scatter-1) / flux) \n", - "\n", - " mean std. dev. \n", - "126 0.234115 0.003568 \n", - "127 1.563707 0.005953 \n", - "124 1.594129 0.002369 \n", - "125 0.285761 0.001676 \n", - "122 0.011089 0.000248 \n", - "123 0.000000 0.000000 \n", - "120 0.000000 0.000000 \n", - "121 0.000000 0.000000 \n", - "118 0.000000 0.000000 \n", - "119 0.000000 0.000000 " + " cell group in group out nuclide moment mean std. dev. moment\n", + "126 10002 1 1 H-1 P0 0.234115 0.003568 P0\n", + "127 10002 1 1 O-16 P0 1.563707 0.005953 P0\n", + "124 10002 1 2 H-1 P0 1.594129 0.002369 P0\n", + "125 10002 1 2 O-16 P0 0.285761 0.001676 P0\n", + "122 10002 1 3 H-1 P0 0.011089 0.000248 P0\n", + "123 10002 1 3 O-16 P0 0.000000 0.000000 P0\n", + "120 10002 1 4 H-1 P0 0.000000 0.000000 P0\n", + "121 10002 1 4 O-16 P0 0.000000 0.000000 P0\n", + "118 10002 1 5 H-1 P0 0.000000 0.000000 P0\n", + "119 10002 1 5 O-16 P0 0.000000 0.000000 P0" ] }, "execution_count": 19, @@ -1019,7 +1018,6 @@ " cell\n", " group in\n", " nuclide\n", - " score\n", " mean\n", " std. dev.\n", " \n", @@ -1030,7 +1028,6 @@ " 10000\n", " 1\n", " U-235\n", - " ((total - scatter-1) / flux)\n", " 20.611692\n", " 0.104237\n", " \n", @@ -1039,7 +1036,6 @@ " 10000\n", " 1\n", " U-238\n", - " ((total - scatter-1) / flux)\n", " 9.585358\n", " 0.013808\n", " \n", @@ -1048,7 +1044,6 @@ " 10000\n", " 1\n", " O-16\n", - " ((total - scatter-1) / flux)\n", " 3.164190\n", " 0.005049\n", " \n", @@ -1057,7 +1052,6 @@ " 10000\n", " 2\n", " U-235\n", - " ((total - scatter-1) / flux)\n", " 485.413426\n", " 0.996410\n", " \n", @@ -1066,7 +1060,6 @@ " 10000\n", " 2\n", " U-238\n", - " ((total - scatter-1) / flux)\n", " 11.190386\n", " 0.028731\n", " \n", @@ -1075,7 +1068,6 @@ " 10000\n", " 2\n", " O-16\n", - " ((total - scatter-1) / flux)\n", " 3.794859\n", " 0.011139\n", " \n", @@ -1084,13 +1076,13 @@ "" ], "text/plain": [ - " cell group in nuclide score mean std. dev.\n", - "3 10000 1 U-235 ((total - scatter-1) / flux) 2.06e+01 1.04e-01\n", - "4 10000 1 U-238 ((total - scatter-1) / flux) 9.59e+00 1.38e-02\n", - "5 10000 1 O-16 ((total - scatter-1) / flux) 3.16e+00 5.05e-03\n", - "0 10000 2 U-235 ((total - scatter-1) / flux) 4.85e+02 9.96e-01\n", - "1 10000 2 U-238 ((total - scatter-1) / flux) 1.12e+01 2.87e-02\n", - "2 10000 2 O-16 ((total - scatter-1) / flux) 3.79e+00 1.11e-02" + " cell group in nuclide mean std. dev.\n", + "3 10000 1 U-235 20.611692 0.104237\n", + "4 10000 1 U-238 9.585358 0.013808\n", + "5 10000 1 O-16 3.164190 0.005049\n", + "0 10000 2 U-235 485.413426 0.996410\n", + "1 10000 2 U-238 11.190386 0.028731\n", + "2 10000 2 O-16 3.794859 0.011139" ] }, "execution_count": 22, @@ -1202,161 +1194,161 @@ "[ NORMAL ] Iteration 5:\tk_eff = 0.625810\tres = 2.417E-02\n", "[ NORMAL ] Iteration 6:\tk_eff = 0.606678\tres = 2.675E-02\n", "[ NORMAL ] Iteration 7:\tk_eff = 0.587485\tres = 3.057E-02\n", - "[ NORMAL ] Iteration 8:\tk_eff = 0.569029\tres = 3.164E-02\n", + "[ NORMAL ] Iteration 8:\tk_eff = 0.569028\tres = 3.164E-02\n", "[ NORMAL ] Iteration 9:\tk_eff = 0.551707\tres = 3.142E-02\n", - "[ NORMAL ] Iteration 10:\tk_eff = 0.536035\tres = 3.044E-02\n", - "[ NORMAL ] Iteration 11:\tk_eff = 0.522275\tres = 2.841E-02\n", - "[ NORMAL ] Iteration 12:\tk_eff = 0.510610\tres = 2.567E-02\n", - "[ NORMAL ] Iteration 13:\tk_eff = 0.501106\tres = 2.234E-02\n", - "[ NORMAL ] Iteration 14:\tk_eff = 0.493832\tres = 1.861E-02\n", - "[ NORMAL ] Iteration 15:\tk_eff = 0.488781\tres = 1.452E-02\n", - "[ NORMAL ] Iteration 16:\tk_eff = 0.485924\tres = 1.023E-02\n", - "[ NORMAL ] Iteration 17:\tk_eff = 0.485211\tres = 5.846E-03\n", - "[ NORMAL ] Iteration 18:\tk_eff = 0.486571\tres = 1.467E-03\n", - "[ NORMAL ] Iteration 19:\tk_eff = 0.489905\tres = 2.802E-03\n", - "[ NORMAL ] Iteration 20:\tk_eff = 0.495105\tres = 6.853E-03\n", - "[ NORMAL ] Iteration 21:\tk_eff = 0.502056\tres = 1.061E-02\n", - "[ NORMAL ] Iteration 22:\tk_eff = 0.510630\tres = 1.404E-02\n", - "[ NORMAL ] Iteration 23:\tk_eff = 0.520696\tres = 1.708E-02\n", - "[ NORMAL ] Iteration 24:\tk_eff = 0.532120\tres = 1.971E-02\n", - "[ NORMAL ] Iteration 25:\tk_eff = 0.544768\tres = 2.194E-02\n", - "[ NORMAL ] Iteration 26:\tk_eff = 0.558505\tres = 2.377E-02\n", - "[ NORMAL ] Iteration 27:\tk_eff = 0.573200\tres = 2.522E-02\n", - "[ NORMAL ] Iteration 28:\tk_eff = 0.588723\tres = 2.631E-02\n", - "[ NORMAL ] Iteration 29:\tk_eff = 0.604951\tres = 2.708E-02\n", - "[ NORMAL ] Iteration 30:\tk_eff = 0.621765\tres = 2.756E-02\n", - "[ NORMAL ] Iteration 31:\tk_eff = 0.639053\tres = 2.779E-02\n", - "[ NORMAL ] Iteration 32:\tk_eff = 0.656709\tres = 2.780E-02\n", - "[ NORMAL ] Iteration 33:\tk_eff = 0.674632\tres = 2.763E-02\n", - "[ NORMAL ] Iteration 34:\tk_eff = 0.692730\tres = 2.729E-02\n", - "[ NORMAL ] Iteration 35:\tk_eff = 0.710919\tres = 2.683E-02\n", - "[ NORMAL ] Iteration 36:\tk_eff = 0.729118\tres = 2.626E-02\n", - "[ NORMAL ] Iteration 37:\tk_eff = 0.747258\tres = 2.560E-02\n", - "[ NORMAL ] Iteration 38:\tk_eff = 0.765273\tres = 2.488E-02\n", - "[ NORMAL ] Iteration 39:\tk_eff = 0.783104\tres = 2.411E-02\n", - "[ NORMAL ] Iteration 40:\tk_eff = 0.800701\tres = 2.330E-02\n", - "[ NORMAL ] Iteration 41:\tk_eff = 0.818017\tres = 2.247E-02\n", - "[ NORMAL ] Iteration 42:\tk_eff = 0.835012\tres = 2.163E-02\n", - "[ NORMAL ] Iteration 43:\tk_eff = 0.851651\tres = 2.078E-02\n", - "[ NORMAL ] Iteration 44:\tk_eff = 0.867906\tres = 1.993E-02\n", - "[ NORMAL ] Iteration 45:\tk_eff = 0.883750\tres = 1.909E-02\n", - "[ NORMAL ] Iteration 46:\tk_eff = 0.899164\tres = 1.826E-02\n", - "[ NORMAL ] Iteration 47:\tk_eff = 0.914131\tres = 1.744E-02\n", - "[ NORMAL ] Iteration 48:\tk_eff = 0.928638\tres = 1.665E-02\n", - "[ NORMAL ] Iteration 49:\tk_eff = 0.942676\tres = 1.587E-02\n", - "[ NORMAL ] Iteration 50:\tk_eff = 0.956239\tres = 1.512E-02\n", - "[ NORMAL ] Iteration 51:\tk_eff = 0.969323\tres = 1.439E-02\n", - "[ NORMAL ] Iteration 52:\tk_eff = 0.981927\tres = 1.368E-02\n", - "[ NORMAL ] Iteration 53:\tk_eff = 0.994054\tres = 1.300E-02\n", - "[ NORMAL ] Iteration 54:\tk_eff = 1.005705\tres = 1.235E-02\n", - "[ NORMAL ] Iteration 55:\tk_eff = 1.016886\tres = 1.172E-02\n", - "[ NORMAL ] Iteration 56:\tk_eff = 1.027604\tres = 1.112E-02\n", - "[ NORMAL ] Iteration 57:\tk_eff = 1.037866\tres = 1.054E-02\n", - "[ NORMAL ] Iteration 58:\tk_eff = 1.047682\tres = 9.987E-03\n", - "[ NORMAL ] Iteration 59:\tk_eff = 1.057062\tres = 9.458E-03\n", - "[ NORMAL ] Iteration 60:\tk_eff = 1.066016\tres = 8.953E-03\n", - "[ NORMAL ] Iteration 61:\tk_eff = 1.074556\tres = 8.471E-03\n", - "[ NORMAL ] Iteration 62:\tk_eff = 1.082693\tres = 8.011E-03\n", - "[ NORMAL ] Iteration 63:\tk_eff = 1.090442\tres = 7.573E-03\n", - "[ NORMAL ] Iteration 64:\tk_eff = 1.097813\tres = 7.156E-03\n", - "[ NORMAL ] Iteration 65:\tk_eff = 1.104820\tres = 6.760E-03\n", - "[ NORMAL ] Iteration 66:\tk_eff = 1.111477\tres = 6.383E-03\n", - "[ NORMAL ] Iteration 67:\tk_eff = 1.117795\tres = 6.025E-03\n", - "[ NORMAL ] Iteration 68:\tk_eff = 1.123789\tres = 5.685E-03\n", - "[ NORMAL ] Iteration 69:\tk_eff = 1.129470\tres = 5.362E-03\n", - "[ NORMAL ] Iteration 70:\tk_eff = 1.134853\tres = 5.056E-03\n", - "[ NORMAL ] Iteration 71:\tk_eff = 1.139950\tres = 4.766E-03\n", - "[ NORMAL ] Iteration 72:\tk_eff = 1.144772\tres = 4.491E-03\n", - "[ NORMAL ] Iteration 73:\tk_eff = 1.149333\tres = 4.230E-03\n", - "[ NORMAL ] Iteration 74:\tk_eff = 1.153643\tres = 3.984E-03\n", - "[ NORMAL ] Iteration 75:\tk_eff = 1.157716\tres = 3.751E-03\n", - "[ NORMAL ] Iteration 76:\tk_eff = 1.161561\tres = 3.530E-03\n", - "[ NORMAL ] Iteration 77:\tk_eff = 1.165190\tres = 3.321E-03\n", - "[ NORMAL ] Iteration 78:\tk_eff = 1.168613\tres = 3.124E-03\n", - "[ NORMAL ] Iteration 79:\tk_eff = 1.171841\tres = 2.938E-03\n", - "[ NORMAL ] Iteration 80:\tk_eff = 1.174883\tres = 2.762E-03\n", - "[ NORMAL ] Iteration 81:\tk_eff = 1.177749\tres = 2.596E-03\n", - "[ NORMAL ] Iteration 82:\tk_eff = 1.180447\tres = 2.439E-03\n", - "[ NORMAL ] Iteration 83:\tk_eff = 1.182988\tres = 2.291E-03\n", - "[ NORMAL ] Iteration 84:\tk_eff = 1.185378\tres = 2.152E-03\n", - "[ NORMAL ] Iteration 85:\tk_eff = 1.187627\tres = 2.021E-03\n", - "[ NORMAL ] Iteration 86:\tk_eff = 1.189741\tres = 1.897E-03\n", - "[ NORMAL ] Iteration 87:\tk_eff = 1.191728\tres = 1.780E-03\n", - "[ NORMAL ] Iteration 88:\tk_eff = 1.193595\tres = 1.670E-03\n", - "[ NORMAL ] Iteration 89:\tk_eff = 1.195349\tres = 1.567E-03\n", - "[ NORMAL ] Iteration 90:\tk_eff = 1.196996\tres = 1.469E-03\n", - "[ NORMAL ] Iteration 91:\tk_eff = 1.198543\tres = 1.378E-03\n", - "[ NORMAL ] Iteration 92:\tk_eff = 1.199994\tres = 1.292E-03\n", - "[ NORMAL ] Iteration 93:\tk_eff = 1.201355\tres = 1.211E-03\n", - "[ NORMAL ] Iteration 94:\tk_eff = 1.202632\tres = 1.134E-03\n", - "[ NORMAL ] Iteration 95:\tk_eff = 1.203829\tres = 1.063E-03\n", - "[ NORMAL ] Iteration 96:\tk_eff = 1.204952\tres = 9.956E-04\n", - "[ NORMAL ] Iteration 97:\tk_eff = 1.206004\tres = 9.324E-04\n", - "[ NORMAL ] Iteration 98:\tk_eff = 1.206989\tres = 8.730E-04\n", - "[ NORMAL ] Iteration 99:\tk_eff = 1.207913\tres = 8.173E-04\n", - "[ NORMAL ] Iteration 100:\tk_eff = 1.208777\tres = 7.650E-04\n", - "[ NORMAL ] Iteration 101:\tk_eff = 1.209587\tres = 7.159E-04\n", - "[ NORMAL ] Iteration 102:\tk_eff = 1.210345\tres = 6.698E-04\n", - "[ NORMAL ] Iteration 103:\tk_eff = 1.211054\tres = 6.266E-04\n", - "[ NORMAL ] Iteration 104:\tk_eff = 1.211718\tres = 5.861E-04\n", - "[ NORMAL ] Iteration 105:\tk_eff = 1.212339\tres = 5.481E-04\n", - "[ NORMAL ] Iteration 106:\tk_eff = 1.212920\tres = 5.125E-04\n", - "[ NORMAL ] Iteration 107:\tk_eff = 1.213463\tres = 4.792E-04\n", - "[ NORMAL ] Iteration 108:\tk_eff = 1.213971\tres = 4.479E-04\n", - "[ NORMAL ] Iteration 109:\tk_eff = 1.214446\tres = 4.186E-04\n", - "[ NORMAL ] Iteration 110:\tk_eff = 1.214890\tres = 3.912E-04\n", - "[ NORMAL ] Iteration 111:\tk_eff = 1.215305\tres = 3.655E-04\n", - "[ NORMAL ] Iteration 112:\tk_eff = 1.215692\tres = 3.414E-04\n", - "[ NORMAL ] Iteration 113:\tk_eff = 1.216055\tres = 3.189E-04\n", - "[ NORMAL ] Iteration 114:\tk_eff = 1.216393\tres = 2.979E-04\n", - "[ NORMAL ] Iteration 115:\tk_eff = 1.216709\tres = 2.782E-04\n", - "[ NORMAL ] Iteration 116:\tk_eff = 1.217004\tres = 2.597E-04\n", - "[ NORMAL ] Iteration 117:\tk_eff = 1.217279\tres = 2.425E-04\n", - "[ NORMAL ] Iteration 118:\tk_eff = 1.217536\tres = 2.263E-04\n", - "[ NORMAL ] Iteration 119:\tk_eff = 1.217776\tres = 2.113E-04\n", - "[ NORMAL ] Iteration 120:\tk_eff = 1.218000\tres = 1.971E-04\n", - "[ NORMAL ] Iteration 121:\tk_eff = 1.218209\tres = 1.840E-04\n", - "[ NORMAL ] Iteration 122:\tk_eff = 1.218404\tres = 1.716E-04\n", - "[ NORMAL ] Iteration 123:\tk_eff = 1.218587\tres = 1.601E-04\n", - "[ NORMAL ] Iteration 124:\tk_eff = 1.218756\tres = 1.494E-04\n", - "[ NORMAL ] Iteration 125:\tk_eff = 1.218915\tres = 1.393E-04\n", - "[ NORMAL ] Iteration 126:\tk_eff = 1.219062\tres = 1.299E-04\n", - "[ NORMAL ] Iteration 127:\tk_eff = 1.219200\tres = 1.212E-04\n", - "[ NORMAL ] Iteration 128:\tk_eff = 1.219329\tres = 1.130E-04\n", - "[ NORMAL ] Iteration 129:\tk_eff = 1.219448\tres = 1.054E-04\n", - "[ NORMAL ] Iteration 130:\tk_eff = 1.219560\tres = 9.822E-05\n", - "[ NORMAL ] Iteration 131:\tk_eff = 1.219664\tres = 9.156E-05\n", - "[ NORMAL ] Iteration 132:\tk_eff = 1.219761\tres = 8.534E-05\n", - "[ NORMAL ] Iteration 133:\tk_eff = 1.219851\tres = 7.954E-05\n", - "[ NORMAL ] Iteration 134:\tk_eff = 1.219936\tres = 7.413E-05\n", - "[ NORMAL ] Iteration 135:\tk_eff = 1.220014\tres = 6.908E-05\n", - "[ NORMAL ] Iteration 136:\tk_eff = 1.220087\tres = 6.437E-05\n", - "[ NORMAL ] Iteration 137:\tk_eff = 1.220156\tres = 5.997E-05\n", - "[ NORMAL ] Iteration 138:\tk_eff = 1.220219\tres = 5.587E-05\n", - "[ NORMAL ] Iteration 139:\tk_eff = 1.220278\tres = 5.205E-05\n", - "[ NORMAL ] Iteration 140:\tk_eff = 1.220333\tres = 4.848E-05\n", - "[ NORMAL ] Iteration 141:\tk_eff = 1.220385\tres = 4.516E-05\n", - "[ NORMAL ] Iteration 142:\tk_eff = 1.220433\tres = 4.206E-05\n", - "[ NORMAL ] Iteration 143:\tk_eff = 1.220477\tres = 3.917E-05\n", - "[ NORMAL ] Iteration 144:\tk_eff = 1.220518\tres = 3.648E-05\n", - "[ NORMAL ] Iteration 145:\tk_eff = 1.220557\tres = 3.397E-05\n", - "[ NORMAL ] Iteration 146:\tk_eff = 1.220593\tres = 3.163E-05\n", - "[ NORMAL ] Iteration 147:\tk_eff = 1.220627\tres = 2.945E-05\n", - "[ NORMAL ] Iteration 148:\tk_eff = 1.220658\tres = 2.742E-05\n", - "[ NORMAL ] Iteration 149:\tk_eff = 1.220687\tres = 2.552E-05\n", - "[ NORMAL ] Iteration 150:\tk_eff = 1.220714\tres = 2.376E-05\n", - "[ NORMAL ] Iteration 151:\tk_eff = 1.220739\tres = 2.212E-05\n", - "[ NORMAL ] Iteration 152:\tk_eff = 1.220762\tres = 2.059E-05\n", - "[ NORMAL ] Iteration 153:\tk_eff = 1.220784\tres = 1.916E-05\n", - "[ NORMAL ] Iteration 154:\tk_eff = 1.220804\tres = 1.783E-05\n", - "[ NORMAL ] Iteration 155:\tk_eff = 1.220823\tres = 1.660E-05\n", - "[ NORMAL ] Iteration 156:\tk_eff = 1.220841\tres = 1.545E-05\n", - "[ NORMAL ] Iteration 157:\tk_eff = 1.220857\tres = 1.437E-05\n", - "[ NORMAL ] Iteration 158:\tk_eff = 1.220872\tres = 1.338E-05\n", - "[ NORMAL ] Iteration 159:\tk_eff = 1.220886\tres = 1.245E-05\n", - "[ NORMAL ] Iteration 160:\tk_eff = 1.220899\tres = 1.158E-05\n", - "[ NORMAL ] Iteration 161:\tk_eff = 1.220912\tres = 1.077E-05\n", - "[ NORMAL ] Iteration 162:\tk_eff = 1.220923\tres = 1.002E-05\n" + "[ NORMAL ] Iteration 10:\tk_eff = 0.536034\tres = 3.044E-02\n", + "[ NORMAL ] Iteration 11:\tk_eff = 0.522274\tres = 2.841E-02\n", + "[ NORMAL ] Iteration 12:\tk_eff = 0.510609\tres = 2.567E-02\n", + "[ NORMAL ] Iteration 13:\tk_eff = 0.501105\tres = 2.234E-02\n", + "[ NORMAL ] Iteration 14:\tk_eff = 0.493831\tres = 1.861E-02\n", + "[ NORMAL ] Iteration 15:\tk_eff = 0.488780\tres = 1.452E-02\n", + "[ NORMAL ] Iteration 16:\tk_eff = 0.485922\tres = 1.023E-02\n", + "[ NORMAL ] Iteration 17:\tk_eff = 0.485209\tres = 5.846E-03\n", + "[ NORMAL ] Iteration 18:\tk_eff = 0.486569\tres = 1.467E-03\n", + "[ NORMAL ] Iteration 19:\tk_eff = 0.489903\tres = 2.801E-03\n", + "[ NORMAL ] Iteration 20:\tk_eff = 0.495102\tres = 6.853E-03\n", + "[ NORMAL ] Iteration 21:\tk_eff = 0.502053\tres = 1.061E-02\n", + "[ NORMAL ] Iteration 22:\tk_eff = 0.510627\tres = 1.404E-02\n", + "[ NORMAL ] Iteration 23:\tk_eff = 0.520692\tres = 1.708E-02\n", + "[ NORMAL ] Iteration 24:\tk_eff = 0.532116\tres = 1.971E-02\n", + "[ NORMAL ] Iteration 25:\tk_eff = 0.544763\tres = 2.194E-02\n", + "[ NORMAL ] Iteration 26:\tk_eff = 0.558500\tres = 2.377E-02\n", + "[ NORMAL ] Iteration 27:\tk_eff = 0.573195\tres = 2.522E-02\n", + "[ NORMAL ] Iteration 28:\tk_eff = 0.588717\tres = 2.631E-02\n", + "[ NORMAL ] Iteration 29:\tk_eff = 0.604945\tres = 2.708E-02\n", + "[ NORMAL ] Iteration 30:\tk_eff = 0.621759\tres = 2.756E-02\n", + "[ NORMAL ] Iteration 31:\tk_eff = 0.639046\tres = 2.779E-02\n", + "[ NORMAL ] Iteration 32:\tk_eff = 0.656701\tres = 2.780E-02\n", + "[ NORMAL ] Iteration 33:\tk_eff = 0.674624\tres = 2.763E-02\n", + "[ NORMAL ] Iteration 34:\tk_eff = 0.692722\tres = 2.729E-02\n", + "[ NORMAL ] Iteration 35:\tk_eff = 0.710909\tres = 2.683E-02\n", + "[ NORMAL ] Iteration 36:\tk_eff = 0.729109\tres = 2.626E-02\n", + "[ NORMAL ] Iteration 37:\tk_eff = 0.747248\tres = 2.560E-02\n", + "[ NORMAL ] Iteration 38:\tk_eff = 0.765262\tres = 2.488E-02\n", + "[ NORMAL ] Iteration 39:\tk_eff = 0.783093\tres = 2.411E-02\n", + "[ NORMAL ] Iteration 40:\tk_eff = 0.800689\tres = 2.330E-02\n", + "[ NORMAL ] Iteration 41:\tk_eff = 0.818004\tres = 2.247E-02\n", + "[ NORMAL ] Iteration 42:\tk_eff = 0.834999\tres = 2.163E-02\n", + "[ NORMAL ] Iteration 43:\tk_eff = 0.851638\tres = 2.078E-02\n", + "[ NORMAL ] Iteration 44:\tk_eff = 0.867891\tres = 1.993E-02\n", + "[ NORMAL ] Iteration 45:\tk_eff = 0.883735\tres = 1.909E-02\n", + "[ NORMAL ] Iteration 46:\tk_eff = 0.899148\tres = 1.826E-02\n", + "[ NORMAL ] Iteration 47:\tk_eff = 0.914114\tres = 1.744E-02\n", + "[ NORMAL ] Iteration 48:\tk_eff = 0.928621\tres = 1.664E-02\n", + "[ NORMAL ] Iteration 49:\tk_eff = 0.942659\tres = 1.587E-02\n", + "[ NORMAL ] Iteration 50:\tk_eff = 0.956221\tres = 1.512E-02\n", + "[ NORMAL ] Iteration 51:\tk_eff = 0.969305\tres = 1.439E-02\n", + "[ NORMAL ] Iteration 52:\tk_eff = 0.981909\tres = 1.368E-02\n", + "[ NORMAL ] Iteration 53:\tk_eff = 0.994034\tres = 1.300E-02\n", + "[ NORMAL ] Iteration 54:\tk_eff = 1.005685\tres = 1.235E-02\n", + "[ NORMAL ] Iteration 55:\tk_eff = 1.016866\tres = 1.172E-02\n", + "[ NORMAL ] Iteration 56:\tk_eff = 1.027583\tres = 1.112E-02\n", + "[ NORMAL ] Iteration 57:\tk_eff = 1.037845\tres = 1.054E-02\n", + "[ NORMAL ] Iteration 58:\tk_eff = 1.047661\tres = 9.986E-03\n", + "[ NORMAL ] Iteration 59:\tk_eff = 1.057040\tres = 9.458E-03\n", + "[ NORMAL ] Iteration 60:\tk_eff = 1.065993\tres = 8.952E-03\n", + "[ NORMAL ] Iteration 61:\tk_eff = 1.074533\tres = 8.470E-03\n", + "[ NORMAL ] Iteration 62:\tk_eff = 1.082670\tres = 8.011E-03\n", + "[ NORMAL ] Iteration 63:\tk_eff = 1.090418\tres = 7.573E-03\n", + "[ NORMAL ] Iteration 64:\tk_eff = 1.097789\tres = 7.156E-03\n", + "[ NORMAL ] Iteration 65:\tk_eff = 1.104796\tres = 6.760E-03\n", + "[ NORMAL ] Iteration 66:\tk_eff = 1.111452\tres = 6.383E-03\n", + "[ NORMAL ] Iteration 67:\tk_eff = 1.117770\tres = 6.025E-03\n", + "[ NORMAL ] Iteration 68:\tk_eff = 1.123764\tres = 5.685E-03\n", + "[ NORMAL ] Iteration 69:\tk_eff = 1.129445\tres = 5.362E-03\n", + "[ NORMAL ] Iteration 70:\tk_eff = 1.134828\tres = 5.056E-03\n", + "[ NORMAL ] Iteration 71:\tk_eff = 1.139924\tres = 4.766E-03\n", + "[ NORMAL ] Iteration 72:\tk_eff = 1.144746\tres = 4.491E-03\n", + "[ NORMAL ] Iteration 73:\tk_eff = 1.149306\tres = 4.230E-03\n", + "[ NORMAL ] Iteration 74:\tk_eff = 1.153617\tres = 3.984E-03\n", + "[ NORMAL ] Iteration 75:\tk_eff = 1.157689\tres = 3.750E-03\n", + "[ NORMAL ] Iteration 76:\tk_eff = 1.161534\tres = 3.530E-03\n", + "[ NORMAL ] Iteration 77:\tk_eff = 1.165163\tres = 3.321E-03\n", + "[ NORMAL ] Iteration 78:\tk_eff = 1.168586\tres = 3.124E-03\n", + "[ NORMAL ] Iteration 79:\tk_eff = 1.171813\tres = 2.938E-03\n", + "[ NORMAL ] Iteration 80:\tk_eff = 1.174855\tres = 2.762E-03\n", + "[ NORMAL ] Iteration 81:\tk_eff = 1.177721\tres = 2.596E-03\n", + "[ NORMAL ] Iteration 82:\tk_eff = 1.180419\tres = 2.439E-03\n", + "[ NORMAL ] Iteration 83:\tk_eff = 1.182960\tres = 2.291E-03\n", + "[ NORMAL ] Iteration 84:\tk_eff = 1.185350\tres = 2.152E-03\n", + "[ NORMAL ] Iteration 85:\tk_eff = 1.187598\tres = 2.021E-03\n", + "[ NORMAL ] Iteration 86:\tk_eff = 1.189712\tres = 1.897E-03\n", + "[ NORMAL ] Iteration 87:\tk_eff = 1.191699\tres = 1.780E-03\n", + "[ NORMAL ] Iteration 88:\tk_eff = 1.193567\tres = 1.670E-03\n", + "[ NORMAL ] Iteration 89:\tk_eff = 1.195320\tres = 1.567E-03\n", + "[ NORMAL ] Iteration 90:\tk_eff = 1.196967\tres = 1.469E-03\n", + "[ NORMAL ] Iteration 91:\tk_eff = 1.198513\tres = 1.378E-03\n", + "[ NORMAL ] Iteration 92:\tk_eff = 1.199964\tres = 1.292E-03\n", + "[ NORMAL ] Iteration 93:\tk_eff = 1.201326\tres = 1.211E-03\n", + "[ NORMAL ] Iteration 94:\tk_eff = 1.202602\tres = 1.134E-03\n", + "[ NORMAL ] Iteration 95:\tk_eff = 1.203800\tres = 1.063E-03\n", + "[ NORMAL ] Iteration 96:\tk_eff = 1.204922\tres = 9.955E-04\n", + "[ NORMAL ] Iteration 97:\tk_eff = 1.205974\tres = 9.323E-04\n", + "[ NORMAL ] Iteration 98:\tk_eff = 1.206959\tres = 8.730E-04\n", + "[ NORMAL ] Iteration 99:\tk_eff = 1.207883\tres = 8.173E-04\n", + "[ NORMAL ] Iteration 100:\tk_eff = 1.208747\tres = 7.649E-04\n", + "[ NORMAL ] Iteration 101:\tk_eff = 1.209557\tres = 7.159E-04\n", + "[ NORMAL ] Iteration 102:\tk_eff = 1.210315\tres = 6.698E-04\n", + "[ NORMAL ] Iteration 103:\tk_eff = 1.211024\tres = 6.266E-04\n", + "[ NORMAL ] Iteration 104:\tk_eff = 1.211688\tres = 5.861E-04\n", + "[ NORMAL ] Iteration 105:\tk_eff = 1.212309\tres = 5.481E-04\n", + "[ NORMAL ] Iteration 106:\tk_eff = 1.212890\tres = 5.125E-04\n", + "[ NORMAL ] Iteration 107:\tk_eff = 1.213433\tres = 4.791E-04\n", + "[ NORMAL ] Iteration 108:\tk_eff = 1.213941\tres = 4.479E-04\n", + "[ NORMAL ] Iteration 109:\tk_eff = 1.214416\tres = 4.186E-04\n", + "[ NORMAL ] Iteration 110:\tk_eff = 1.214860\tres = 3.912E-04\n", + "[ NORMAL ] Iteration 111:\tk_eff = 1.215275\tres = 3.655E-04\n", + "[ NORMAL ] Iteration 112:\tk_eff = 1.215662\tres = 3.414E-04\n", + "[ NORMAL ] Iteration 113:\tk_eff = 1.216024\tres = 3.189E-04\n", + "[ NORMAL ] Iteration 114:\tk_eff = 1.216362\tres = 2.979E-04\n", + "[ NORMAL ] Iteration 115:\tk_eff = 1.216678\tres = 2.781E-04\n", + "[ NORMAL ] Iteration 116:\tk_eff = 1.216973\tres = 2.597E-04\n", + "[ NORMAL ] Iteration 117:\tk_eff = 1.217249\tres = 2.425E-04\n", + "[ NORMAL ] Iteration 118:\tk_eff = 1.217506\tres = 2.263E-04\n", + "[ NORMAL ] Iteration 119:\tk_eff = 1.217746\tres = 2.112E-04\n", + "[ NORMAL ] Iteration 120:\tk_eff = 1.217970\tres = 1.971E-04\n", + "[ NORMAL ] Iteration 121:\tk_eff = 1.218179\tres = 1.840E-04\n", + "[ NORMAL ] Iteration 122:\tk_eff = 1.218374\tres = 1.716E-04\n", + "[ NORMAL ] Iteration 123:\tk_eff = 1.218556\tres = 1.601E-04\n", + "[ NORMAL ] Iteration 124:\tk_eff = 1.218726\tres = 1.494E-04\n", + "[ NORMAL ] Iteration 125:\tk_eff = 1.218884\tres = 1.393E-04\n", + "[ NORMAL ] Iteration 126:\tk_eff = 1.219032\tres = 1.299E-04\n", + "[ NORMAL ] Iteration 127:\tk_eff = 1.219170\tres = 1.212E-04\n", + "[ NORMAL ] Iteration 128:\tk_eff = 1.219298\tres = 1.130E-04\n", + "[ NORMAL ] Iteration 129:\tk_eff = 1.219418\tres = 1.053E-04\n", + "[ NORMAL ] Iteration 130:\tk_eff = 1.219529\tres = 9.821E-05\n", + "[ NORMAL ] Iteration 131:\tk_eff = 1.219633\tres = 9.155E-05\n", + "[ NORMAL ] Iteration 132:\tk_eff = 1.219730\tres = 8.534E-05\n", + "[ NORMAL ] Iteration 133:\tk_eff = 1.219821\tres = 7.954E-05\n", + "[ NORMAL ] Iteration 134:\tk_eff = 1.219905\tres = 7.412E-05\n", + "[ NORMAL ] Iteration 135:\tk_eff = 1.219984\tres = 6.907E-05\n", + "[ NORMAL ] Iteration 136:\tk_eff = 1.220057\tres = 6.436E-05\n", + "[ NORMAL ] Iteration 137:\tk_eff = 1.220125\tres = 5.997E-05\n", + "[ NORMAL ] Iteration 138:\tk_eff = 1.220188\tres = 5.587E-05\n", + "[ NORMAL ] Iteration 139:\tk_eff = 1.220248\tres = 5.205E-05\n", + "[ NORMAL ] Iteration 140:\tk_eff = 1.220303\tres = 4.848E-05\n", + "[ NORMAL ] Iteration 141:\tk_eff = 1.220354\tres = 4.516E-05\n", + "[ NORMAL ] Iteration 142:\tk_eff = 1.220402\tres = 4.206E-05\n", + "[ NORMAL ] Iteration 143:\tk_eff = 1.220446\tres = 3.917E-05\n", + "[ NORMAL ] Iteration 144:\tk_eff = 1.220488\tres = 3.648E-05\n", + "[ NORMAL ] Iteration 145:\tk_eff = 1.220526\tres = 3.397E-05\n", + "[ NORMAL ] Iteration 146:\tk_eff = 1.220562\tres = 3.163E-05\n", + "[ NORMAL ] Iteration 147:\tk_eff = 1.220596\tres = 2.945E-05\n", + "[ NORMAL ] Iteration 148:\tk_eff = 1.220627\tres = 2.742E-05\n", + "[ NORMAL ] Iteration 149:\tk_eff = 1.220656\tres = 2.552E-05\n", + "[ NORMAL ] Iteration 150:\tk_eff = 1.220683\tres = 2.376E-05\n", + "[ NORMAL ] Iteration 151:\tk_eff = 1.220708\tres = 2.212E-05\n", + "[ NORMAL ] Iteration 152:\tk_eff = 1.220732\tres = 2.059E-05\n", + "[ NORMAL ] Iteration 153:\tk_eff = 1.220753\tres = 1.916E-05\n", + "[ NORMAL ] Iteration 154:\tk_eff = 1.220774\tres = 1.783E-05\n", + "[ NORMAL ] Iteration 155:\tk_eff = 1.220792\tres = 1.660E-05\n", + "[ NORMAL ] Iteration 156:\tk_eff = 1.220810\tres = 1.545E-05\n", + "[ NORMAL ] Iteration 157:\tk_eff = 1.220826\tres = 1.437E-05\n", + "[ NORMAL ] Iteration 158:\tk_eff = 1.220841\tres = 1.337E-05\n", + "[ NORMAL ] Iteration 159:\tk_eff = 1.220856\tres = 1.244E-05\n", + "[ NORMAL ] Iteration 160:\tk_eff = 1.220869\tres = 1.158E-05\n", + "[ NORMAL ] Iteration 161:\tk_eff = 1.220881\tres = 1.077E-05\n", + "[ NORMAL ] Iteration 162:\tk_eff = 1.220892\tres = 1.002E-05\n" ] } ], @@ -1389,8 +1381,8 @@ "output_type": "stream", "text": [ "openmc keff = 1.223474\n", - "openmoc keff = 1.220923\n", - "bias [pcm]: -255.0\n" + "openmoc keff = 1.220892\n", + "bias [pcm]: -258.1\n" ] } ], @@ -1467,235 +1459,235 @@ "[ NORMAL ] Computing the eigenvalue...\n", "[ NORMAL ] Iteration 0:\tk_eff = 0.495816\tres = 0.000E+00\n", "[ NORMAL ] Iteration 1:\tk_eff = 0.557477\tres = 5.042E-01\n", - "[ NORMAL ] Iteration 2:\tk_eff = 0.518301\tres = 1.244E-01\n", - "[ NORMAL ] Iteration 3:\tk_eff = 0.509212\tres = 7.027E-02\n", + "[ NORMAL ] Iteration 2:\tk_eff = 0.518300\tres = 1.244E-01\n", + "[ NORMAL ] Iteration 3:\tk_eff = 0.509211\tres = 7.027E-02\n", "[ NORMAL ] Iteration 4:\tk_eff = 0.496489\tres = 1.754E-02\n", "[ NORMAL ] Iteration 5:\tk_eff = 0.488581\tres = 2.498E-02\n", "[ NORMAL ] Iteration 6:\tk_eff = 0.482897\tres = 1.593E-02\n", "[ NORMAL ] Iteration 7:\tk_eff = 0.479775\tres = 1.163E-02\n", - "[ NORMAL ] Iteration 8:\tk_eff = 0.478835\tres = 6.464E-03\n", - "[ NORMAL ] Iteration 9:\tk_eff = 0.479872\tres = 1.960E-03\n", - "[ NORMAL ] Iteration 10:\tk_eff = 0.482685\tres = 2.166E-03\n", - "[ NORMAL ] Iteration 11:\tk_eff = 0.487085\tres = 5.861E-03\n", - "[ NORMAL ] Iteration 12:\tk_eff = 0.492901\tres = 9.116E-03\n", - "[ NORMAL ] Iteration 13:\tk_eff = 0.499973\tres = 1.194E-02\n", - "[ NORMAL ] Iteration 14:\tk_eff = 0.508155\tres = 1.435E-02\n", - "[ NORMAL ] Iteration 15:\tk_eff = 0.517314\tres = 1.637E-02\n", - "[ NORMAL ] Iteration 16:\tk_eff = 0.527327\tres = 1.802E-02\n", - "[ NORMAL ] Iteration 17:\tk_eff = 0.538081\tres = 1.936E-02\n", - "[ NORMAL ] Iteration 18:\tk_eff = 0.549475\tres = 2.039E-02\n", - "[ NORMAL ] Iteration 19:\tk_eff = 0.561413\tres = 2.117E-02\n", - "[ NORMAL ] Iteration 20:\tk_eff = 0.573810\tres = 2.173E-02\n", - "[ NORMAL ] Iteration 21:\tk_eff = 0.586589\tres = 2.208E-02\n", - "[ NORMAL ] Iteration 22:\tk_eff = 0.599678\tres = 2.227E-02\n", - "[ NORMAL ] Iteration 23:\tk_eff = 0.613012\tres = 2.231E-02\n", - "[ NORMAL ] Iteration 24:\tk_eff = 0.626532\tres = 2.224E-02\n", - "[ NORMAL ] Iteration 25:\tk_eff = 0.640186\tres = 2.206E-02\n", - "[ NORMAL ] Iteration 26:\tk_eff = 0.653925\tres = 2.179E-02\n", - "[ NORMAL ] Iteration 27:\tk_eff = 0.667706\tres = 2.146E-02\n", - "[ NORMAL ] Iteration 28:\tk_eff = 0.681489\tres = 2.107E-02\n", - "[ NORMAL ] Iteration 29:\tk_eff = 0.695240\tres = 2.064E-02\n", - "[ NORMAL ] Iteration 30:\tk_eff = 0.708927\tres = 2.018E-02\n", - "[ NORMAL ] Iteration 31:\tk_eff = 0.722522\tres = 1.969E-02\n", - "[ NORMAL ] Iteration 32:\tk_eff = 0.736000\tres = 1.918E-02\n", - "[ NORMAL ] Iteration 33:\tk_eff = 0.749339\tres = 1.865E-02\n", - "[ NORMAL ] Iteration 34:\tk_eff = 0.762519\tres = 1.812E-02\n", - "[ NORMAL ] Iteration 35:\tk_eff = 0.775522\tres = 1.759E-02\n", - "[ NORMAL ] Iteration 36:\tk_eff = 0.788335\tres = 1.705E-02\n", - "[ NORMAL ] Iteration 37:\tk_eff = 0.800944\tres = 1.652E-02\n", - "[ NORMAL ] Iteration 38:\tk_eff = 0.813338\tres = 1.599E-02\n", - "[ NORMAL ] Iteration 39:\tk_eff = 0.825507\tres = 1.547E-02\n", - "[ NORMAL ] Iteration 40:\tk_eff = 0.837444\tres = 1.496E-02\n", - "[ NORMAL ] Iteration 41:\tk_eff = 0.849142\tres = 1.446E-02\n", - "[ NORMAL ] Iteration 42:\tk_eff = 0.860595\tres = 1.397E-02\n", - "[ NORMAL ] Iteration 43:\tk_eff = 0.871801\tres = 1.349E-02\n", - "[ NORMAL ] Iteration 44:\tk_eff = 0.882755\tres = 1.302E-02\n", - "[ NORMAL ] Iteration 45:\tk_eff = 0.893455\tres = 1.256E-02\n", - "[ NORMAL ] Iteration 46:\tk_eff = 0.903901\tres = 1.212E-02\n", - "[ NORMAL ] Iteration 47:\tk_eff = 0.914091\tres = 1.169E-02\n", - "[ NORMAL ] Iteration 48:\tk_eff = 0.924027\tres = 1.127E-02\n", - "[ NORMAL ] Iteration 49:\tk_eff = 0.933708\tres = 1.087E-02\n", - "[ NORMAL ] Iteration 50:\tk_eff = 0.943137\tres = 1.048E-02\n", - "[ NORMAL ] Iteration 51:\tk_eff = 0.952314\tres = 1.010E-02\n", - "[ NORMAL ] Iteration 52:\tk_eff = 0.961243\tres = 9.731E-03\n", - "[ NORMAL ] Iteration 53:\tk_eff = 0.969927\tres = 9.376E-03\n", - "[ NORMAL ] Iteration 54:\tk_eff = 0.978367\tres = 9.033E-03\n", - "[ NORMAL ] Iteration 55:\tk_eff = 0.986568\tres = 8.702E-03\n", - "[ NORMAL ] Iteration 56:\tk_eff = 0.994534\tres = 8.382E-03\n", - "[ NORMAL ] Iteration 57:\tk_eff = 1.002267\tres = 8.074E-03\n", - "[ NORMAL ] Iteration 58:\tk_eff = 1.009773\tres = 7.776E-03\n", - "[ NORMAL ] Iteration 59:\tk_eff = 1.017055\tres = 7.489E-03\n", - "[ NORMAL ] Iteration 60:\tk_eff = 1.024118\tres = 7.212E-03\n", - "[ NORMAL ] Iteration 61:\tk_eff = 1.030966\tres = 6.944E-03\n", - "[ NORMAL ] Iteration 62:\tk_eff = 1.037603\tres = 6.687E-03\n", - "[ NORMAL ] Iteration 63:\tk_eff = 1.044036\tres = 6.438E-03\n", - "[ NORMAL ] Iteration 64:\tk_eff = 1.050267\tres = 6.199E-03\n", - "[ NORMAL ] Iteration 65:\tk_eff = 1.056302\tres = 5.968E-03\n", - "[ NORMAL ] Iteration 66:\tk_eff = 1.062145\tres = 5.746E-03\n", - "[ NORMAL ] Iteration 67:\tk_eff = 1.067802\tres = 5.532E-03\n", - "[ NORMAL ] Iteration 68:\tk_eff = 1.073276\tres = 5.326E-03\n", - "[ NORMAL ] Iteration 69:\tk_eff = 1.078573\tres = 5.127E-03\n", - "[ NORMAL ] Iteration 70:\tk_eff = 1.083698\tres = 4.935E-03\n", - "[ NORMAL ] Iteration 71:\tk_eff = 1.088654\tres = 4.751E-03\n", - "[ NORMAL ] Iteration 72:\tk_eff = 1.093447\tres = 4.573E-03\n", - "[ NORMAL ] Iteration 73:\tk_eff = 1.098080\tres = 4.402E-03\n", - "[ NORMAL ] Iteration 74:\tk_eff = 1.102560\tres = 4.238E-03\n", - "[ NORMAL ] Iteration 75:\tk_eff = 1.106889\tres = 4.079E-03\n", - "[ NORMAL ] Iteration 76:\tk_eff = 1.111072\tres = 3.926E-03\n", - "[ NORMAL ] Iteration 77:\tk_eff = 1.115114\tres = 3.779E-03\n", - "[ NORMAL ] Iteration 78:\tk_eff = 1.119018\tres = 3.638E-03\n", - "[ NORMAL ] Iteration 79:\tk_eff = 1.122790\tres = 3.501E-03\n", - "[ NORMAL ] Iteration 80:\tk_eff = 1.126432\tres = 3.370E-03\n", - "[ NORMAL ] Iteration 81:\tk_eff = 1.129949\tres = 3.244E-03\n", - "[ NORMAL ] Iteration 82:\tk_eff = 1.133344\tres = 3.122E-03\n", - "[ NORMAL ] Iteration 83:\tk_eff = 1.136622\tres = 3.005E-03\n", - "[ NORMAL ] Iteration 84:\tk_eff = 1.139786\tres = 2.892E-03\n", - "[ NORMAL ] Iteration 85:\tk_eff = 1.142840\tres = 2.784E-03\n", - "[ NORMAL ] Iteration 86:\tk_eff = 1.145786\tres = 2.679E-03\n", - "[ NORMAL ] Iteration 87:\tk_eff = 1.148630\tres = 2.579E-03\n", - "[ NORMAL ] Iteration 88:\tk_eff = 1.151373\tres = 2.482E-03\n", - "[ NORMAL ] Iteration 89:\tk_eff = 1.154020\tres = 2.388E-03\n", - "[ NORMAL ] Iteration 90:\tk_eff = 1.156573\tres = 2.299E-03\n", - "[ NORMAL ] Iteration 91:\tk_eff = 1.159035\tres = 2.212E-03\n", - "[ NORMAL ] Iteration 92:\tk_eff = 1.161410\tres = 2.129E-03\n", - "[ NORMAL ] Iteration 93:\tk_eff = 1.163700\tres = 2.049E-03\n", - "[ NORMAL ] Iteration 94:\tk_eff = 1.165909\tres = 1.972E-03\n", - "[ NORMAL ] Iteration 95:\tk_eff = 1.168038\tres = 1.898E-03\n", - "[ NORMAL ] Iteration 96:\tk_eff = 1.170091\tres = 1.826E-03\n", - "[ NORMAL ] Iteration 97:\tk_eff = 1.172070\tres = 1.758E-03\n", - "[ NORMAL ] Iteration 98:\tk_eff = 1.173977\tres = 1.691E-03\n", - "[ NORMAL ] Iteration 99:\tk_eff = 1.175816\tres = 1.628E-03\n", - "[ NORMAL ] Iteration 100:\tk_eff = 1.177589\tres = 1.566E-03\n", - "[ NORMAL ] Iteration 101:\tk_eff = 1.179297\tres = 1.507E-03\n", - "[ NORMAL ] Iteration 102:\tk_eff = 1.180943\tres = 1.451E-03\n", - "[ NORMAL ] Iteration 103:\tk_eff = 1.182529\tres = 1.396E-03\n", - "[ NORMAL ] Iteration 104:\tk_eff = 1.184058\tres = 1.343E-03\n", - "[ NORMAL ] Iteration 105:\tk_eff = 1.185530\tres = 1.293E-03\n", - "[ NORMAL ] Iteration 106:\tk_eff = 1.186949\tres = 1.244E-03\n", - "[ NORMAL ] Iteration 107:\tk_eff = 1.188316\tres = 1.197E-03\n", - "[ NORMAL ] Iteration 108:\tk_eff = 1.189633\tres = 1.152E-03\n", - "[ NORMAL ] Iteration 109:\tk_eff = 1.190902\tres = 1.108E-03\n", - "[ NORMAL ] Iteration 110:\tk_eff = 1.192124\tres = 1.066E-03\n", - "[ NORMAL ] Iteration 111:\tk_eff = 1.193301\tres = 1.026E-03\n", - "[ NORMAL ] Iteration 112:\tk_eff = 1.194435\tres = 9.874E-04\n", - "[ NORMAL ] Iteration 113:\tk_eff = 1.195527\tres = 9.501E-04\n", - "[ NORMAL ] Iteration 114:\tk_eff = 1.196578\tres = 9.142E-04\n", - "[ NORMAL ] Iteration 115:\tk_eff = 1.197591\tres = 8.797E-04\n", - "[ NORMAL ] Iteration 116:\tk_eff = 1.198566\tres = 8.464E-04\n", - "[ NORMAL ] Iteration 117:\tk_eff = 1.199506\tres = 8.144E-04\n", - "[ NORMAL ] Iteration 118:\tk_eff = 1.200410\tres = 7.836E-04\n", - "[ NORMAL ] Iteration 119:\tk_eff = 1.201281\tres = 7.540E-04\n", - "[ NORMAL ] Iteration 120:\tk_eff = 1.202119\tres = 7.255E-04\n", - "[ NORMAL ] Iteration 121:\tk_eff = 1.202927\tres = 6.980E-04\n", - "[ NORMAL ] Iteration 122:\tk_eff = 1.203704\tres = 6.716E-04\n", - "[ NORMAL ] Iteration 123:\tk_eff = 1.204452\tres = 6.462E-04\n", - "[ NORMAL ] Iteration 124:\tk_eff = 1.205173\tres = 6.217E-04\n", - "[ NORMAL ] Iteration 125:\tk_eff = 1.205867\tres = 5.982E-04\n", - "[ NORMAL ] Iteration 126:\tk_eff = 1.206534\tres = 5.755E-04\n", - "[ NORMAL ] Iteration 127:\tk_eff = 1.207177\tres = 5.537E-04\n", - "[ NORMAL ] Iteration 128:\tk_eff = 1.207796\tres = 5.327E-04\n", - "[ NORMAL ] Iteration 129:\tk_eff = 1.208391\tres = 5.126E-04\n", - "[ NORMAL ] Iteration 130:\tk_eff = 1.208965\tres = 4.931E-04\n", - "[ NORMAL ] Iteration 131:\tk_eff = 1.209517\tres = 4.744E-04\n", - "[ NORMAL ] Iteration 132:\tk_eff = 1.210048\tres = 4.565E-04\n", - "[ NORMAL ] Iteration 133:\tk_eff = 1.210559\tres = 4.391E-04\n", - "[ NORMAL ] Iteration 134:\tk_eff = 1.211051\tres = 4.225E-04\n", - "[ NORMAL ] Iteration 135:\tk_eff = 1.211525\tres = 4.065E-04\n", - "[ NORMAL ] Iteration 136:\tk_eff = 1.211980\tres = 3.910E-04\n", - "[ NORMAL ] Iteration 137:\tk_eff = 1.212419\tres = 3.762E-04\n", - "[ NORMAL ] Iteration 138:\tk_eff = 1.212841\tres = 3.619E-04\n", - "[ NORMAL ] Iteration 139:\tk_eff = 1.213247\tres = 3.482E-04\n", - "[ NORMAL ] Iteration 140:\tk_eff = 1.213638\tres = 3.350E-04\n", - "[ NORMAL ] Iteration 141:\tk_eff = 1.214015\tres = 3.222E-04\n", - "[ NORMAL ] Iteration 142:\tk_eff = 1.214377\tres = 3.100E-04\n", - "[ NORMAL ] Iteration 143:\tk_eff = 1.214725\tres = 2.982E-04\n", - "[ NORMAL ] Iteration 144:\tk_eff = 1.215060\tres = 2.869E-04\n", - "[ NORMAL ] Iteration 145:\tk_eff = 1.215383\tres = 2.760E-04\n", - "[ NORMAL ] Iteration 146:\tk_eff = 1.215693\tres = 2.655E-04\n", - "[ NORMAL ] Iteration 147:\tk_eff = 1.215992\tres = 2.554E-04\n", - "[ NORMAL ] Iteration 148:\tk_eff = 1.216280\tres = 2.457E-04\n", - "[ NORMAL ] Iteration 149:\tk_eff = 1.216556\tres = 2.364E-04\n", - "[ NORMAL ] Iteration 150:\tk_eff = 1.216822\tres = 2.274E-04\n", - "[ NORMAL ] Iteration 151:\tk_eff = 1.217078\tres = 2.187E-04\n", - "[ NORMAL ] Iteration 152:\tk_eff = 1.217325\tres = 2.104E-04\n", - "[ NORMAL ] Iteration 153:\tk_eff = 1.217562\tres = 2.024E-04\n", - "[ NORMAL ] Iteration 154:\tk_eff = 1.217790\tres = 1.947E-04\n", - "[ NORMAL ] Iteration 155:\tk_eff = 1.218009\tres = 1.873E-04\n", - "[ NORMAL ] Iteration 156:\tk_eff = 1.218220\tres = 1.802E-04\n", - "[ NORMAL ] Iteration 157:\tk_eff = 1.218423\tres = 1.733E-04\n", - "[ NORMAL ] Iteration 158:\tk_eff = 1.218619\tres = 1.667E-04\n", - "[ NORMAL ] Iteration 159:\tk_eff = 1.218807\tres = 1.604E-04\n", - "[ NORMAL ] Iteration 160:\tk_eff = 1.218988\tres = 1.543E-04\n", - "[ NORMAL ] Iteration 161:\tk_eff = 1.219162\tres = 1.484E-04\n", - "[ NORMAL ] Iteration 162:\tk_eff = 1.219329\tres = 1.428E-04\n", - "[ NORMAL ] Iteration 163:\tk_eff = 1.219490\tres = 1.373E-04\n", - "[ NORMAL ] Iteration 164:\tk_eff = 1.219645\tres = 1.321E-04\n", - "[ NORMAL ] Iteration 165:\tk_eff = 1.219794\tres = 1.271E-04\n", - "[ NORMAL ] Iteration 166:\tk_eff = 1.219938\tres = 1.222E-04\n", - "[ NORMAL ] Iteration 167:\tk_eff = 1.220076\tres = 1.176E-04\n", - "[ NORMAL ] Iteration 168:\tk_eff = 1.220208\tres = 1.131E-04\n", - "[ NORMAL ] Iteration 169:\tk_eff = 1.220336\tres = 1.088E-04\n", - "[ NORMAL ] Iteration 170:\tk_eff = 1.220459\tres = 1.046E-04\n", - "[ NORMAL ] Iteration 171:\tk_eff = 1.220577\tres = 1.006E-04\n", - "[ NORMAL ] Iteration 172:\tk_eff = 1.220691\tres = 9.681E-05\n", - "[ NORMAL ] Iteration 173:\tk_eff = 1.220800\tres = 9.312E-05\n", - "[ NORMAL ] Iteration 174:\tk_eff = 1.220905\tres = 8.957E-05\n", - "[ NORMAL ] Iteration 175:\tk_eff = 1.221006\tres = 8.615E-05\n", - "[ NORMAL ] Iteration 176:\tk_eff = 1.221104\tres = 8.287E-05\n", - "[ NORMAL ] Iteration 177:\tk_eff = 1.221197\tres = 7.971E-05\n", - "[ NORMAL ] Iteration 178:\tk_eff = 1.221287\tres = 7.667E-05\n", - "[ NORMAL ] Iteration 179:\tk_eff = 1.221374\tres = 7.374E-05\n", - "[ NORMAL ] Iteration 180:\tk_eff = 1.221457\tres = 7.093E-05\n", - "[ NORMAL ] Iteration 181:\tk_eff = 1.221537\tres = 6.823E-05\n", - "[ NORMAL ] Iteration 182:\tk_eff = 1.221615\tres = 6.562E-05\n", - "[ NORMAL ] Iteration 183:\tk_eff = 1.221689\tres = 6.312E-05\n", - "[ NORMAL ] Iteration 184:\tk_eff = 1.221760\tres = 6.071E-05\n", - "[ NORMAL ] Iteration 185:\tk_eff = 1.221829\tres = 5.840E-05\n", - "[ NORMAL ] Iteration 186:\tk_eff = 1.221895\tres = 5.617E-05\n", - "[ NORMAL ] Iteration 187:\tk_eff = 1.221958\tres = 5.402E-05\n", - "[ NORMAL ] Iteration 188:\tk_eff = 1.222019\tres = 5.196E-05\n", - "[ NORMAL ] Iteration 189:\tk_eff = 1.222078\tres = 4.998E-05\n", - "[ NORMAL ] Iteration 190:\tk_eff = 1.222134\tres = 4.807E-05\n", - "[ NORMAL ] Iteration 191:\tk_eff = 1.222189\tres = 4.624E-05\n", - "[ NORMAL ] Iteration 192:\tk_eff = 1.222241\tres = 4.447E-05\n", - "[ NORMAL ] Iteration 193:\tk_eff = 1.222291\tres = 4.277E-05\n", - "[ NORMAL ] Iteration 194:\tk_eff = 1.222340\tres = 4.114E-05\n", - "[ NORMAL ] Iteration 195:\tk_eff = 1.222386\tres = 3.957E-05\n", - "[ NORMAL ] Iteration 196:\tk_eff = 1.222431\tres = 3.806E-05\n", - "[ NORMAL ] Iteration 197:\tk_eff = 1.222474\tres = 3.661E-05\n", - "[ NORMAL ] Iteration 198:\tk_eff = 1.222515\tres = 3.521E-05\n", - "[ NORMAL ] Iteration 199:\tk_eff = 1.222555\tres = 3.386E-05\n", - "[ NORMAL ] Iteration 200:\tk_eff = 1.222594\tres = 3.257E-05\n", - "[ NORMAL ] Iteration 201:\tk_eff = 1.222630\tres = 3.133E-05\n", - "[ NORMAL ] Iteration 202:\tk_eff = 1.222666\tres = 3.013E-05\n", - "[ NORMAL ] Iteration 203:\tk_eff = 1.222700\tres = 2.898E-05\n", - "[ NORMAL ] Iteration 204:\tk_eff = 1.222733\tres = 2.787E-05\n", - "[ NORMAL ] Iteration 205:\tk_eff = 1.222764\tres = 2.681E-05\n", - "[ NORMAL ] Iteration 206:\tk_eff = 1.222795\tres = 2.578E-05\n", - "[ NORMAL ] Iteration 207:\tk_eff = 1.222824\tres = 2.480E-05\n", - "[ NORMAL ] Iteration 208:\tk_eff = 1.222852\tres = 2.385E-05\n", - "[ NORMAL ] Iteration 209:\tk_eff = 1.222879\tres = 2.294E-05\n", - "[ NORMAL ] Iteration 210:\tk_eff = 1.222905\tres = 2.206E-05\n", - "[ NORMAL ] Iteration 211:\tk_eff = 1.222930\tres = 2.122E-05\n", - "[ NORMAL ] Iteration 212:\tk_eff = 1.222954\tres = 2.041E-05\n", - "[ NORMAL ] Iteration 213:\tk_eff = 1.222977\tres = 1.963E-05\n", - "[ NORMAL ] Iteration 214:\tk_eff = 1.222999\tres = 1.888E-05\n", - "[ NORMAL ] Iteration 215:\tk_eff = 1.223020\tres = 1.816E-05\n", - "[ NORMAL ] Iteration 216:\tk_eff = 1.223041\tres = 1.747E-05\n", - "[ NORMAL ] Iteration 217:\tk_eff = 1.223061\tres = 1.680E-05\n", - "[ NORMAL ] Iteration 218:\tk_eff = 1.223080\tres = 1.616E-05\n", - "[ NORMAL ] Iteration 219:\tk_eff = 1.223098\tres = 1.554E-05\n", - "[ NORMAL ] Iteration 220:\tk_eff = 1.223116\tres = 1.495E-05\n", - "[ NORMAL ] Iteration 221:\tk_eff = 1.223132\tres = 1.437E-05\n", - "[ NORMAL ] Iteration 222:\tk_eff = 1.223149\tres = 1.382E-05\n", - "[ NORMAL ] Iteration 223:\tk_eff = 1.223164\tres = 1.330E-05\n", - "[ NORMAL ] Iteration 224:\tk_eff = 1.223179\tres = 1.279E-05\n", - "[ NORMAL ] Iteration 225:\tk_eff = 1.223194\tres = 1.230E-05\n", - "[ NORMAL ] Iteration 226:\tk_eff = 1.223208\tres = 1.183E-05\n", - "[ NORMAL ] Iteration 227:\tk_eff = 1.223221\tres = 1.138E-05\n", - "[ NORMAL ] Iteration 228:\tk_eff = 1.223234\tres = 1.094E-05\n", - "[ NORMAL ] Iteration 229:\tk_eff = 1.223246\tres = 1.052E-05\n", - "[ NORMAL ] Iteration 230:\tk_eff = 1.223258\tres = 1.012E-05\n" + "[ NORMAL ] Iteration 8:\tk_eff = 0.478834\tres = 6.465E-03\n", + "[ NORMAL ] Iteration 9:\tk_eff = 0.479871\tres = 1.960E-03\n", + "[ NORMAL ] Iteration 10:\tk_eff = 0.482684\tres = 2.166E-03\n", + "[ NORMAL ] Iteration 11:\tk_eff = 0.487084\tres = 5.861E-03\n", + "[ NORMAL ] Iteration 12:\tk_eff = 0.492900\tres = 9.116E-03\n", + "[ NORMAL ] Iteration 13:\tk_eff = 0.499971\tres = 1.194E-02\n", + "[ NORMAL ] Iteration 14:\tk_eff = 0.508153\tres = 1.435E-02\n", + "[ NORMAL ] Iteration 15:\tk_eff = 0.517312\tres = 1.637E-02\n", + "[ NORMAL ] Iteration 16:\tk_eff = 0.527324\tres = 1.802E-02\n", + "[ NORMAL ] Iteration 17:\tk_eff = 0.538079\tres = 1.935E-02\n", + "[ NORMAL ] Iteration 18:\tk_eff = 0.549472\tres = 2.039E-02\n", + "[ NORMAL ] Iteration 19:\tk_eff = 0.561410\tres = 2.117E-02\n", + "[ NORMAL ] Iteration 20:\tk_eff = 0.573807\tres = 2.173E-02\n", + "[ NORMAL ] Iteration 21:\tk_eff = 0.586585\tres = 2.208E-02\n", + "[ NORMAL ] Iteration 22:\tk_eff = 0.599674\tres = 2.227E-02\n", + "[ NORMAL ] Iteration 23:\tk_eff = 0.613007\tres = 2.231E-02\n", + "[ NORMAL ] Iteration 24:\tk_eff = 0.626528\tres = 2.223E-02\n", + "[ NORMAL ] Iteration 25:\tk_eff = 0.640181\tres = 2.206E-02\n", + "[ NORMAL ] Iteration 26:\tk_eff = 0.653920\tres = 2.179E-02\n", + "[ NORMAL ] Iteration 27:\tk_eff = 0.667700\tres = 2.146E-02\n", + "[ NORMAL ] Iteration 28:\tk_eff = 0.681483\tres = 2.107E-02\n", + "[ NORMAL ] Iteration 29:\tk_eff = 0.695234\tres = 2.064E-02\n", + "[ NORMAL ] Iteration 30:\tk_eff = 0.708921\tres = 2.018E-02\n", + "[ NORMAL ] Iteration 31:\tk_eff = 0.722515\tres = 1.969E-02\n", + "[ NORMAL ] Iteration 32:\tk_eff = 0.735993\tres = 1.918E-02\n", + "[ NORMAL ] Iteration 33:\tk_eff = 0.749331\tres = 1.865E-02\n", + "[ NORMAL ] Iteration 34:\tk_eff = 0.762510\tres = 1.812E-02\n", + "[ NORMAL ] Iteration 35:\tk_eff = 0.775514\tres = 1.759E-02\n", + "[ NORMAL ] Iteration 36:\tk_eff = 0.788326\tres = 1.705E-02\n", + "[ NORMAL ] Iteration 37:\tk_eff = 0.800934\tres = 1.652E-02\n", + "[ NORMAL ] Iteration 38:\tk_eff = 0.813328\tres = 1.599E-02\n", + "[ NORMAL ] Iteration 39:\tk_eff = 0.825497\tres = 1.547E-02\n", + "[ NORMAL ] Iteration 40:\tk_eff = 0.837433\tres = 1.496E-02\n", + "[ NORMAL ] Iteration 41:\tk_eff = 0.849131\tres = 1.446E-02\n", + "[ NORMAL ] Iteration 42:\tk_eff = 0.860584\tres = 1.397E-02\n", + "[ NORMAL ] Iteration 43:\tk_eff = 0.871789\tres = 1.349E-02\n", + "[ NORMAL ] Iteration 44:\tk_eff = 0.882742\tres = 1.302E-02\n", + "[ NORMAL ] Iteration 45:\tk_eff = 0.893442\tres = 1.256E-02\n", + "[ NORMAL ] Iteration 46:\tk_eff = 0.903888\tres = 1.212E-02\n", + "[ NORMAL ] Iteration 47:\tk_eff = 0.914078\tres = 1.169E-02\n", + "[ NORMAL ] Iteration 48:\tk_eff = 0.924013\tres = 1.127E-02\n", + "[ NORMAL ] Iteration 49:\tk_eff = 0.933694\tres = 1.087E-02\n", + "[ NORMAL ] Iteration 50:\tk_eff = 0.943122\tres = 1.048E-02\n", + "[ NORMAL ] Iteration 51:\tk_eff = 0.952299\tres = 1.010E-02\n", + "[ NORMAL ] Iteration 52:\tk_eff = 0.961228\tres = 9.731E-03\n", + "[ NORMAL ] Iteration 53:\tk_eff = 0.969911\tres = 9.376E-03\n", + "[ NORMAL ] Iteration 54:\tk_eff = 0.978351\tres = 9.033E-03\n", + "[ NORMAL ] Iteration 55:\tk_eff = 0.986552\tres = 8.702E-03\n", + "[ NORMAL ] Iteration 56:\tk_eff = 0.994517\tres = 8.382E-03\n", + "[ NORMAL ] Iteration 57:\tk_eff = 1.002250\tres = 8.074E-03\n", + "[ NORMAL ] Iteration 58:\tk_eff = 1.009756\tres = 7.776E-03\n", + "[ NORMAL ] Iteration 59:\tk_eff = 1.017037\tres = 7.488E-03\n", + "[ NORMAL ] Iteration 60:\tk_eff = 1.024100\tres = 7.211E-03\n", + "[ NORMAL ] Iteration 61:\tk_eff = 1.030948\tres = 6.944E-03\n", + "[ NORMAL ] Iteration 62:\tk_eff = 1.037585\tres = 6.687E-03\n", + "[ NORMAL ] Iteration 63:\tk_eff = 1.044017\tres = 6.438E-03\n", + "[ NORMAL ] Iteration 64:\tk_eff = 1.050248\tres = 6.199E-03\n", + "[ NORMAL ] Iteration 65:\tk_eff = 1.056282\tres = 5.968E-03\n", + "[ NORMAL ] Iteration 66:\tk_eff = 1.062125\tres = 5.746E-03\n", + "[ NORMAL ] Iteration 67:\tk_eff = 1.067782\tres = 5.532E-03\n", + "[ NORMAL ] Iteration 68:\tk_eff = 1.073256\tres = 5.325E-03\n", + "[ NORMAL ] Iteration 69:\tk_eff = 1.078553\tres = 5.127E-03\n", + "[ NORMAL ] Iteration 70:\tk_eff = 1.083677\tres = 4.935E-03\n", + "[ NORMAL ] Iteration 71:\tk_eff = 1.088633\tres = 4.751E-03\n", + "[ NORMAL ] Iteration 72:\tk_eff = 1.093425\tres = 4.573E-03\n", + "[ NORMAL ] Iteration 73:\tk_eff = 1.098059\tres = 4.402E-03\n", + "[ NORMAL ] Iteration 74:\tk_eff = 1.102538\tres = 4.238E-03\n", + "[ NORMAL ] Iteration 75:\tk_eff = 1.106867\tres = 4.079E-03\n", + "[ NORMAL ] Iteration 76:\tk_eff = 1.111050\tres = 3.926E-03\n", + "[ NORMAL ] Iteration 77:\tk_eff = 1.115091\tres = 3.779E-03\n", + "[ NORMAL ] Iteration 78:\tk_eff = 1.118996\tres = 3.638E-03\n", + "[ NORMAL ] Iteration 79:\tk_eff = 1.122767\tres = 3.501E-03\n", + "[ NORMAL ] Iteration 80:\tk_eff = 1.126408\tres = 3.370E-03\n", + "[ NORMAL ] Iteration 81:\tk_eff = 1.129925\tres = 3.244E-03\n", + "[ NORMAL ] Iteration 82:\tk_eff = 1.133320\tres = 3.122E-03\n", + "[ NORMAL ] Iteration 83:\tk_eff = 1.136598\tres = 3.005E-03\n", + "[ NORMAL ] Iteration 84:\tk_eff = 1.139762\tres = 2.892E-03\n", + "[ NORMAL ] Iteration 85:\tk_eff = 1.142815\tres = 2.784E-03\n", + "[ NORMAL ] Iteration 86:\tk_eff = 1.145762\tres = 2.679E-03\n", + "[ NORMAL ] Iteration 87:\tk_eff = 1.148605\tres = 2.578E-03\n", + "[ NORMAL ] Iteration 88:\tk_eff = 1.151348\tres = 2.482E-03\n", + "[ NORMAL ] Iteration 89:\tk_eff = 1.153995\tres = 2.388E-03\n", + "[ NORMAL ] Iteration 90:\tk_eff = 1.156548\tres = 2.299E-03\n", + "[ NORMAL ] Iteration 91:\tk_eff = 1.159010\tres = 2.212E-03\n", + "[ NORMAL ] Iteration 92:\tk_eff = 1.161384\tres = 2.129E-03\n", + "[ NORMAL ] Iteration 93:\tk_eff = 1.163674\tres = 2.049E-03\n", + "[ NORMAL ] Iteration 94:\tk_eff = 1.165883\tres = 1.972E-03\n", + "[ NORMAL ] Iteration 95:\tk_eff = 1.168012\tres = 1.898E-03\n", + "[ NORMAL ] Iteration 96:\tk_eff = 1.170064\tres = 1.826E-03\n", + "[ NORMAL ] Iteration 97:\tk_eff = 1.172043\tres = 1.757E-03\n", + "[ NORMAL ] Iteration 98:\tk_eff = 1.173951\tres = 1.691E-03\n", + "[ NORMAL ] Iteration 99:\tk_eff = 1.175790\tres = 1.628E-03\n", + "[ NORMAL ] Iteration 100:\tk_eff = 1.177562\tres = 1.566E-03\n", + "[ NORMAL ] Iteration 101:\tk_eff = 1.179270\tres = 1.507E-03\n", + "[ NORMAL ] Iteration 102:\tk_eff = 1.180916\tres = 1.450E-03\n", + "[ NORMAL ] Iteration 103:\tk_eff = 1.182502\tres = 1.396E-03\n", + "[ NORMAL ] Iteration 104:\tk_eff = 1.184030\tres = 1.343E-03\n", + "[ NORMAL ] Iteration 105:\tk_eff = 1.185503\tres = 1.292E-03\n", + "[ NORMAL ] Iteration 106:\tk_eff = 1.186922\tres = 1.244E-03\n", + "[ NORMAL ] Iteration 107:\tk_eff = 1.188289\tres = 1.197E-03\n", + "[ NORMAL ] Iteration 108:\tk_eff = 1.189605\tres = 1.152E-03\n", + "[ NORMAL ] Iteration 109:\tk_eff = 1.190874\tres = 1.108E-03\n", + "[ NORMAL ] Iteration 110:\tk_eff = 1.192096\tres = 1.066E-03\n", + "[ NORMAL ] Iteration 111:\tk_eff = 1.193273\tres = 1.026E-03\n", + "[ NORMAL ] Iteration 112:\tk_eff = 1.194407\tres = 9.873E-04\n", + "[ NORMAL ] Iteration 113:\tk_eff = 1.195498\tres = 9.500E-04\n", + "[ NORMAL ] Iteration 114:\tk_eff = 1.196550\tres = 9.141E-04\n", + "[ NORMAL ] Iteration 115:\tk_eff = 1.197563\tres = 8.796E-04\n", + "[ NORMAL ] Iteration 116:\tk_eff = 1.198538\tres = 8.464E-04\n", + "[ NORMAL ] Iteration 117:\tk_eff = 1.199477\tres = 8.144E-04\n", + "[ NORMAL ] Iteration 118:\tk_eff = 1.200381\tres = 7.836E-04\n", + "[ NORMAL ] Iteration 119:\tk_eff = 1.201252\tres = 7.539E-04\n", + "[ NORMAL ] Iteration 120:\tk_eff = 1.202091\tres = 7.254E-04\n", + "[ NORMAL ] Iteration 121:\tk_eff = 1.202898\tres = 6.980E-04\n", + "[ NORMAL ] Iteration 122:\tk_eff = 1.203675\tres = 6.715E-04\n", + "[ NORMAL ] Iteration 123:\tk_eff = 1.204423\tres = 6.461E-04\n", + "[ NORMAL ] Iteration 124:\tk_eff = 1.205144\tres = 6.217E-04\n", + "[ NORMAL ] Iteration 125:\tk_eff = 1.205837\tres = 5.981E-04\n", + "[ NORMAL ] Iteration 126:\tk_eff = 1.206505\tres = 5.755E-04\n", + "[ NORMAL ] Iteration 127:\tk_eff = 1.207148\tres = 5.537E-04\n", + "[ NORMAL ] Iteration 128:\tk_eff = 1.207766\tres = 5.327E-04\n", + "[ NORMAL ] Iteration 129:\tk_eff = 1.208362\tres = 5.125E-04\n", + "[ NORMAL ] Iteration 130:\tk_eff = 1.208935\tres = 4.931E-04\n", + "[ NORMAL ] Iteration 131:\tk_eff = 1.209487\tres = 4.744E-04\n", + "[ NORMAL ] Iteration 132:\tk_eff = 1.210018\tres = 4.564E-04\n", + "[ NORMAL ] Iteration 133:\tk_eff = 1.210529\tres = 4.391E-04\n", + "[ NORMAL ] Iteration 134:\tk_eff = 1.211021\tres = 4.225E-04\n", + "[ NORMAL ] Iteration 135:\tk_eff = 1.211495\tres = 4.064E-04\n", + "[ NORMAL ] Iteration 136:\tk_eff = 1.211950\tres = 3.910E-04\n", + "[ NORMAL ] Iteration 137:\tk_eff = 1.212389\tres = 3.762E-04\n", + "[ NORMAL ] Iteration 138:\tk_eff = 1.212811\tres = 3.619E-04\n", + "[ NORMAL ] Iteration 139:\tk_eff = 1.213217\tres = 3.482E-04\n", + "[ NORMAL ] Iteration 140:\tk_eff = 1.213608\tres = 3.349E-04\n", + "[ NORMAL ] Iteration 141:\tk_eff = 1.213984\tres = 3.222E-04\n", + "[ NORMAL ] Iteration 142:\tk_eff = 1.214346\tres = 3.100E-04\n", + "[ NORMAL ] Iteration 143:\tk_eff = 1.214695\tres = 2.982E-04\n", + "[ NORMAL ] Iteration 144:\tk_eff = 1.215030\tres = 2.869E-04\n", + "[ NORMAL ] Iteration 145:\tk_eff = 1.215353\tres = 2.760E-04\n", + "[ NORMAL ] Iteration 146:\tk_eff = 1.215663\tres = 2.655E-04\n", + "[ NORMAL ] Iteration 147:\tk_eff = 1.215962\tres = 2.554E-04\n", + "[ NORMAL ] Iteration 148:\tk_eff = 1.216249\tres = 2.457E-04\n", + "[ NORMAL ] Iteration 149:\tk_eff = 1.216526\tres = 2.364E-04\n", + "[ NORMAL ] Iteration 150:\tk_eff = 1.216792\tres = 2.274E-04\n", + "[ NORMAL ] Iteration 151:\tk_eff = 1.217048\tres = 2.187E-04\n", + "[ NORMAL ] Iteration 152:\tk_eff = 1.217294\tres = 2.104E-04\n", + "[ NORMAL ] Iteration 153:\tk_eff = 1.217531\tres = 2.024E-04\n", + "[ NORMAL ] Iteration 154:\tk_eff = 1.217759\tres = 1.947E-04\n", + "[ NORMAL ] Iteration 155:\tk_eff = 1.217979\tres = 1.873E-04\n", + "[ NORMAL ] Iteration 156:\tk_eff = 1.218190\tres = 1.802E-04\n", + "[ NORMAL ] Iteration 157:\tk_eff = 1.218393\tres = 1.733E-04\n", + "[ NORMAL ] Iteration 158:\tk_eff = 1.218588\tres = 1.667E-04\n", + "[ NORMAL ] Iteration 159:\tk_eff = 1.218776\tres = 1.604E-04\n", + "[ NORMAL ] Iteration 160:\tk_eff = 1.218957\tres = 1.543E-04\n", + "[ NORMAL ] Iteration 161:\tk_eff = 1.219131\tres = 1.484E-04\n", + "[ NORMAL ] Iteration 162:\tk_eff = 1.219298\tres = 1.427E-04\n", + "[ NORMAL ] Iteration 163:\tk_eff = 1.219459\tres = 1.373E-04\n", + "[ NORMAL ] Iteration 164:\tk_eff = 1.219614\tres = 1.321E-04\n", + "[ NORMAL ] Iteration 165:\tk_eff = 1.219763\tres = 1.270E-04\n", + "[ NORMAL ] Iteration 166:\tk_eff = 1.219907\tres = 1.222E-04\n", + "[ NORMAL ] Iteration 167:\tk_eff = 1.220045\tres = 1.176E-04\n", + "[ NORMAL ] Iteration 168:\tk_eff = 1.220177\tres = 1.131E-04\n", + "[ NORMAL ] Iteration 169:\tk_eff = 1.220305\tres = 1.088E-04\n", + "[ NORMAL ] Iteration 170:\tk_eff = 1.220428\tres = 1.046E-04\n", + "[ NORMAL ] Iteration 171:\tk_eff = 1.220546\tres = 1.006E-04\n", + "[ NORMAL ] Iteration 172:\tk_eff = 1.220660\tres = 9.680E-05\n", + "[ NORMAL ] Iteration 173:\tk_eff = 1.220769\tres = 9.311E-05\n", + "[ NORMAL ] Iteration 174:\tk_eff = 1.220874\tres = 8.956E-05\n", + "[ NORMAL ] Iteration 175:\tk_eff = 1.220975\tres = 8.614E-05\n", + "[ NORMAL ] Iteration 176:\tk_eff = 1.221073\tres = 8.286E-05\n", + "[ NORMAL ] Iteration 177:\tk_eff = 1.221166\tres = 7.970E-05\n", + "[ NORMAL ] Iteration 178:\tk_eff = 1.221256\tres = 7.666E-05\n", + "[ NORMAL ] Iteration 179:\tk_eff = 1.221343\tres = 7.373E-05\n", + "[ NORMAL ] Iteration 180:\tk_eff = 1.221426\tres = 7.092E-05\n", + "[ NORMAL ] Iteration 181:\tk_eff = 1.221506\tres = 6.822E-05\n", + "[ NORMAL ] Iteration 182:\tk_eff = 1.221583\tres = 6.562E-05\n", + "[ NORMAL ] Iteration 183:\tk_eff = 1.221658\tres = 6.311E-05\n", + "[ NORMAL ] Iteration 184:\tk_eff = 1.221729\tres = 6.070E-05\n", + "[ NORMAL ] Iteration 185:\tk_eff = 1.221797\tres = 5.839E-05\n", + "[ NORMAL ] Iteration 186:\tk_eff = 1.221863\tres = 5.616E-05\n", + "[ NORMAL ] Iteration 187:\tk_eff = 1.221927\tres = 5.402E-05\n", + "[ NORMAL ] Iteration 188:\tk_eff = 1.221988\tres = 5.196E-05\n", + "[ NORMAL ] Iteration 189:\tk_eff = 1.222047\tres = 4.997E-05\n", + "[ NORMAL ] Iteration 190:\tk_eff = 1.222103\tres = 4.807E-05\n", + "[ NORMAL ] Iteration 191:\tk_eff = 1.222158\tres = 4.623E-05\n", + "[ NORMAL ] Iteration 192:\tk_eff = 1.222210\tres = 4.447E-05\n", + "[ NORMAL ] Iteration 193:\tk_eff = 1.222260\tres = 4.277E-05\n", + "[ NORMAL ] Iteration 194:\tk_eff = 1.222308\tres = 4.114E-05\n", + "[ NORMAL ] Iteration 195:\tk_eff = 1.222355\tres = 3.957E-05\n", + "[ NORMAL ] Iteration 196:\tk_eff = 1.222400\tres = 3.805E-05\n", + "[ NORMAL ] Iteration 197:\tk_eff = 1.222443\tres = 3.660E-05\n", + "[ NORMAL ] Iteration 198:\tk_eff = 1.222484\tres = 3.520E-05\n", + "[ NORMAL ] Iteration 199:\tk_eff = 1.222524\tres = 3.386E-05\n", + "[ NORMAL ] Iteration 200:\tk_eff = 1.222562\tres = 3.257E-05\n", + "[ NORMAL ] Iteration 201:\tk_eff = 1.222599\tres = 3.132E-05\n", + "[ NORMAL ] Iteration 202:\tk_eff = 1.222635\tres = 3.013E-05\n", + "[ NORMAL ] Iteration 203:\tk_eff = 1.222669\tres = 2.898E-05\n", + "[ NORMAL ] Iteration 204:\tk_eff = 1.222701\tres = 2.787E-05\n", + "[ NORMAL ] Iteration 205:\tk_eff = 1.222733\tres = 2.680E-05\n", + "[ NORMAL ] Iteration 206:\tk_eff = 1.222763\tres = 2.578E-05\n", + "[ NORMAL ] Iteration 207:\tk_eff = 1.222792\tres = 2.480E-05\n", + "[ NORMAL ] Iteration 208:\tk_eff = 1.222820\tres = 2.385E-05\n", + "[ NORMAL ] Iteration 209:\tk_eff = 1.222847\tres = 2.294E-05\n", + "[ NORMAL ] Iteration 210:\tk_eff = 1.222873\tres = 2.206E-05\n", + "[ NORMAL ] Iteration 211:\tk_eff = 1.222898\tres = 2.122E-05\n", + "[ NORMAL ] Iteration 212:\tk_eff = 1.222922\tres = 2.041E-05\n", + "[ NORMAL ] Iteration 213:\tk_eff = 1.222945\tres = 1.963E-05\n", + "[ NORMAL ] Iteration 214:\tk_eff = 1.222968\tres = 1.888E-05\n", + "[ NORMAL ] Iteration 215:\tk_eff = 1.222989\tres = 1.816E-05\n", + "[ NORMAL ] Iteration 216:\tk_eff = 1.223009\tres = 1.746E-05\n", + "[ NORMAL ] Iteration 217:\tk_eff = 1.223029\tres = 1.680E-05\n", + "[ NORMAL ] Iteration 218:\tk_eff = 1.223048\tres = 1.615E-05\n", + "[ NORMAL ] Iteration 219:\tk_eff = 1.223067\tres = 1.554E-05\n", + "[ NORMAL ] Iteration 220:\tk_eff = 1.223084\tres = 1.494E-05\n", + "[ NORMAL ] Iteration 221:\tk_eff = 1.223101\tres = 1.437E-05\n", + "[ NORMAL ] Iteration 222:\tk_eff = 1.223117\tres = 1.382E-05\n", + "[ NORMAL ] Iteration 223:\tk_eff = 1.223133\tres = 1.329E-05\n", + "[ NORMAL ] Iteration 224:\tk_eff = 1.223148\tres = 1.279E-05\n", + "[ NORMAL ] Iteration 225:\tk_eff = 1.223162\tres = 1.230E-05\n", + "[ NORMAL ] Iteration 226:\tk_eff = 1.223176\tres = 1.183E-05\n", + "[ NORMAL ] Iteration 227:\tk_eff = 1.223190\tres = 1.137E-05\n", + "[ NORMAL ] Iteration 228:\tk_eff = 1.223203\tres = 1.094E-05\n", + "[ NORMAL ] Iteration 229:\tk_eff = 1.223215\tres = 1.052E-05\n", + "[ NORMAL ] Iteration 230:\tk_eff = 1.223227\tres = 1.012E-05\n" ] } ], @@ -1721,8 +1713,8 @@ "output_type": "stream", "text": [ "openmc keff = 1.223474\n", - "openmoc keff = 1.223258\n", - "bias [pcm]: -21.5\n" + "openmoc keff = 1.223227\n", + "bias [pcm]: -24.7\n" ] } ], @@ -1813,7 +1805,7 @@ "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEfCAYAAACqKwpQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XeYVOX1wPHvzGwvgOKCoqigeBQEC4ian6ighhXEEgli\nAwNEQMUSNYq9RkNiogbEBipWsEYEwYqQWIKiKMVXEFERFaTJ9p3y++POsNN3Zpm+5/M8POzcvXPP\nO7O798zbbR6PB6WUUsqfPd0FUEoplXk0OSillAqhyUEppVQITQ5KKaVCaHJQSikVQpODUkqpEHnp\nLoDKLiKyL7DaGJMXdPwC4DxjzIlhntMGeADog/WB5DljzE3e7x0HTALaAjXA5caYhd7r3Qf86Hep\nycaYyUHXPh54A1gTFPZ54EFgvjHm4Ba8zkuAjsaYG+N9bpRrng/8CSgGCoAPgKuNMesTFSPGcpwE\n3ArsinUPWAtcaoxZ0cLrHQnUGmM+T8b7ptJDk4NKhb8ADUB3oBT4TEQWAf8BXgQGGmM+EZHTgFki\nsof3eS8bYy6I4frfGWMOjPC9uBMDQHAS2lkiMh4rMZxqjFkpIvnADcBCETnYGFPnd67dGONOZHy/\na7fDSpz9jTGfiojNW64XRaS7MaYlE5/+gPWz/DzR75tKH00OKhVeAlZ5b3jbRWQp0AP4HzDaGPOJ\n97y3gY5Au0QE9a/liMiewAxgD6xP7TONMddHOX4LsJcxZoyI7A08AuwLNAKTjDEzvNf/ALgL+CPW\nJ/E/GWNmBpXDDtwMjDDGrAQwxjQCN4vIp4DHW1MaglWD+gy4SkQuBcZh1bYMMMYYs9Fb2/onUATY\ngJuMMc9HOh70tnQDPMDn3nJ4ROQ+4Env1zbgRuBc73Ve8b4ml4h0BR4HOgFbgLFAX2AEcKqIdADa\nJOp9U+mlfQ4q6Ywx7xhjvocdTUy/AT4yxmwzxvzbe9wGjAYWGWO2eJ96qIgsEJGvRGSaiLTdiWJc\nDiw0xnTHqk3s7a2hRDru72FggTFGgMHA/d4bHMBugNsY09N7rTvCxD4Q2AV4M/gbxphXjDH13oe/\nBcYbY64SkaOAq4HjvbWi77BupgB/B67wlvlk4IxmjvtbDvwKLBCRc0RkD2OM0xizwfv984BhWDf9\n/bz/xvu9D88aY/YH7sRKKA9iJfk/G2P+keD3TaWRJgeVMiJSADwDvGqM+cDv+FCsvoVLgIu8h78C\n/o31afpQrE+k/4xw6b1F5Mugf38MOmcDMFBEjgEajTHnG2N+jHLcV7Z84CSsPhOMMd8C7wIDvKfk\nAY95v14C7B2mfLsCG2NosllljFnl/Xow8ILfTftRrOThey0jRORAY8xaY8w5zRzfwRhTAxyNdUO/\nFVgvIh95ax1gvd/TvYnb6Y37OxEpAvoDz3rP+zdwZKQXkqD3TaWRNiupeLkBm4jYgm52DsAFICJv\nA3sC+PoCRKQMq3lpHVZTyQ7GmBeAF0RkAPC2iBxqjHkfeN93jojcBcyLUKawfQ5+n1LBSiwOrJtV\nJxGZAtwS5bhPe8BmjNnmd2wL0MH7tcsYU+372nutYL8AHUUkz3vDjWSz39cVgH9HtX/MUVj9FW+J\nSC0w0fseRjoewNsBfiVwpfc9uhiYKyKdsZr0rhKRC72n5wEbsRKcHdjmvYYHqIryWhLxvqk00uSg\n4vULVpt1Z6ymDp8DfI+NMSf4P0FE8oCXgWXGmCv8jncGDjXGzPY+7x0RWQcc6W2LrzHGbPSenofV\nbt0i3pvy3cDdInIA8DrwH2PMm+GOB71et4js4tfc1R74OY7wX2F9qj8VK0HuICI3AVPDPOdnbxyf\nHTGNMT8DE4AJIvJb4CURmRfl+I6buIh0A8qMMZ96r7UWuFpERgFdsRLSq2FGhRVi/dzbA794mwH3\nA76O8JoT8b6pNNJmJRUXb7PEE8Bt3mYiROQwYCTwrwhPuxTY7p8YvAqAJ0Wkh/c6AuyP1S4+FnhQ\nRPJExIF105vT0nKLyEPeIZxg3dB+wuoIDnvc7/U6gfne8iAi+wHHAm/FGtvbEX8DVpv7Ed7r5IvI\nHVj9Ar+GedocrOYcX4IYC8zxPm+BX7/IJ1hJ0xHhePCop8Oxksb+fu/NYMAJrMRqLjpfREq83xsr\nIiO9/SJvABd4nzYQmOutQTQSNIggEe+bSi9NDqolLsVqAvlMRFYCk4FzjDGfRzh/LNA3qE/gdmPM\n11ijVZ4VkS+xRsZc4m13vwPYjnXDWoF187p6J8r8IHCnN84KrNEyb0c57m8ccLz3nJexRg19H09w\nY8xj3vI/IiJfAV9gdcoO8OuQ9j//f1g1mkXeuO2A672jnB7Fan5bAbwHTPA234Q7XhN03ZlYHdsv\niogRka+xfp6V3maeV4DZwBJv3FOxbvIAY4AhIrIG6+fj69N4GfiriAR3SO/0+6bSx6b7OSillAqm\nNQellFIhNDkopZQKoclBKaVUCE0OSimlQmT9PAen0+XZsqWm+RMTYJddSsjFWKmOp7GyL57Gyq5Y\nscSrqCi3RXt+1tcc8vJSN7EyV2OlOp7Gyr54Giu7YiUiXtYnB6WUUomnyUEppVQITQ5KKaVCaHJQ\nSikVQpODUkqpEJoclFJKhdDkoJRSKkTWT4LTRWWVUonw/fffcf/997B16xZcLjc9e/bi4osvp6Cg\nIOZrvPvuW/TvfyKrVhkWLlzA6NFjk1ji5Mr6mkOfPrBuXdSJfkopFZXL5eKGG/7MOeeM4JFHZjBt\n2pMAPPbYI3Fd56mnngCgWzfJ6sQAOVBzOPdcOPnkEh55pI6jjnKluzhKqSy0ePFH7L33vhx2WG8A\nbDYbF110KTabnVmznuXtt98AoF+/4zjvvAu4885b2G23CoxZyc8//8RNN93BJ5/8j9Wrv+K6665m\n6NCzeOmlWdxxxyTOOut0+vU7ni+/XEZhYQl/+9u9PPbYI7Rr144zzzyLNWtW849/TGLy5Id5++03\nmTnzaRwOByIHcfnlVzFt2kNhz7333r/x5ZcrcblcnHHGUAYNGpLQ9yTraw5/+hPcd18do0YV8eST\n+ekujlIqC3333Vq6dTsg4FhhYRG//LKR11+fzZQpjzBlyiO8886b/PDDOgAaGhr4xz8m8/vfD2fe\nvDmcc84IysrK+Mtf/hZwnfXrf6CycjAzZ85k+/Zf+frrVWHLUFNTw8MPT+Heex9g6tRprF//A0uW\nfBz23F9/3cb77/+HBx+cztSp03A6nQl4FwJlfc0BYMAAF6+9VsOIEcUsX27n9tvrydc8oVTWOvbY\nEr78MnFrER14oIuFC6MtemfD7Q7ebhtWrTL06NGTvDzrVtmz5yGsXv0VAIccchgAFRUdWbFiecQr\nl5aWsv/+3QDo0KEDVVVVYc/7/vvv2GuvvSkpKQHgsMN689VXX4Y9t02btnTuvA/XXvsn+vc/kcrK\nwVFeW8vkRHIA6NrVw+uv1zB+fDHDhhXz6KN1tG+vvdVKZaPoN/LE22effXnxxVkBxxoaGvjmmzX4\nb6Xc2NiIzWY1uDgcTckr2nbL/uf5zrXZmvpJfZ/6bbbA6zidjRQWFoY9F+Cee+7HmC958815zJs3\nh3/+c0pMrzVWWd+s5K+8HJ54opY+fVwMHFjC8uU59fKUUklyxBFH8vPPP/Kf/ywEwO12M3Xqv/j+\n+29ZtuwLnE4nTqeTFSuWc8ABEvE6bndsH0hLS0v55ZdfAPj8888A6Nx5H9at+46ammoAPv10CSLd\nw57744/ref755xA5kEsuuZxt27a17IVHkTM1Bx+HA66/voHu3d0MHVrMpEn1DBmS+PY4pVTusNvt\n3HPPZCZNupPHHnuE/Px8jjjiSCZMuIKXX36BCRMuxO32MGTIaey++x4Rr3PAAcIf/ziC8eMvjRrv\nuOMGcPXVl7Fy5XIOPfRwAIqLi7n44su48soJ2Gx2evU6lEMOOZSOHTuGnLvbbhUsW7aUt99+g/z8\nfAYPPjVxb4aXLVp1KEt4Nm7cHvYbn39u54ILihk+vJGrrmrAvpMViYqKciLFSrRUxkp1PI2VffE0\nVnbFiiVezm/2E02vXm7mzath4UIHo0YVEaEfSCmlVJCcTg4AHTp4eOmlWnbd1cPgwSWsXasT5pRS\nqjk5nxwACgrgnnvqGTGikcGDS1i0KLXb9SmlVLZpFckBrGFio0c38tBDdYwbV8S0afm6LpNSSkXQ\napKDzzHHuJg7t4YZM/K58spCGhrSXSKllMo8rS45AOyzj4c5c2rYvNnG735XzIYN2g+hlFL+WmVy\nACgrg+nT6zjuOBeVlSUsXdpq3wqlWr0ff1zPMcf0YcWKZQHH//jHEdx55y1hnzN37mwmT74XsJbq\nBmu5jWnTHgp7/qJFixg/fhTjx49i1KhzeeihKbhcmbtYaKu+I9rtcPXVDdx6az3Dhxfz0ks5NydQ\nKRWjTp325J133trx+KeffuTXX3+N6bnNLdX944/rufvuu7n99r8ydep0Hn74Cb755mtee+3fAedl\n0rwzvRsCQ4Y46drVzciRxaxYYWfixAYcOqBJqValR4+efPLJ/3Y8fvfdtzniiKOor69j6NAhzJgx\nk5KSEiZPvpeuXffbcd4zz8wIu1S3v1deeZGRI0ey224VAOTl5XHHHZN2LOg3fPgZHH30MbRt25ZB\ng4Zw11230djYiN1u59prb8Rms3HDDdfs2Gdi9OjzueOOvzJ9+sOUlJSwdu1atm3bynXX3cQBBxyY\nkPejVdcc/PXo4Wb+/Bo++cTBiBHFxPiBQSmVaPfcQ/sunajo0CZh/9p36UTxA/+KGjYvL49u3YRl\ny74A4P33F3H00f/XbHEjLdXt77vv1nLAAYFLgvsSA1gL6h155NFccMEYHn30QU455TQmT36YM84Y\nyvTpD0eN73Q6ue++BxgzZhyPPfZos+WNVcYlBxEpFZGPReSUVMdu397DrFm1dO7s5uSTS1izRjuq\nlUq5e+7BXp3Y5Qzs1VUUT42eHAD69z+Bd999k59//ony8jYUFxcnJL7NZt+xour69T9wySUXMn78\naK655ood53Tv3gMAY1bu2HTo8MP7sGqViXrtPn36AnDwwb34/vtvE1JeSEFyEJHpIrJBRJYFHa8U\nESMiq0XkWr9vXQMErp2bQvn5cPfd9Ywb18gpp5TwzjvavqRUSl15Je7SsoRe0l1aRu34Cc2e16fP\nkXz88WLee+9djjuu/47jkZbNjsSXAC655EK+/HIlXbp0Zdky6xbYqdOeTJ78MDfddPuO1VYB8vJ8\nm9DYdvQ9NDY6sdnsAfGDy+BbCdZ6TuI+0Kaiz+FxYDIww3dARBzAFOAkYB2wWEReBfYEVgBFKShX\nVOef30i3bm7++Mcixo9vYPz4xnQXSanW4cor2TTiwrSEzs/P54ADhDlz/s2UKY/u2GynpKSUTZt+\nobBwT5Yv/yJk2e7gpbp9CcCnffv2XHbZOHr1OoLOnfcG4OOP/0dBQUFIGQ46qDtLlnzMSSdV8tln\nn3DggQdRUlLKli2b8Xg8bN68ifXr1+04//PPP+WEE05i+fLP2XffLgl7L5KeHIwxC0Vk36DDfYHV\nxpg1ACLyHHAaUAaUAt2BWhGZa4wJ3Z4pRY46ysXrr9cwcmQxy5c7mDGj+ecopbJb//4nsnXrFsrK\nmmovZ545jGuuuYK9996HLl26hjynuaW6Kyo68M9//pPbbrsdl8uF0+lkn3325ZZb7gw5d8yYcdx1\n1+3Mnv0KeXn5TJx4I23atKFPn76MGTOC/ffvRrduTcmpoaGBP//5cn7++Wduuun2BLwDlpQs2e1N\nDq8ZYw72Ph4KVBpjxngfnw8caYy5xPv4AuAXY8xrMVw+6S+gpgZGjYJvvoGXX4ZOnZIdUSmlmnft\ntdcycOBA+vfv3/zJoaK2QWXkUFZjzOPxnJ+KNdLvvx+mTSunTx8306fX0rt3cis0mbb2u8bKrFip\njqexMjNWXV0j27bVhr1uDPs5RL12upLDD0Bnv8d7eY9lLJsNJk6Ezp3rOP/8Ym66qZ7hw3WHOaVU\n+lx//S1Ju3a6ksNioJuIdMFKCsOBc9JUlrgMHOji5ZdrGTHC6oe4+eZ68jKy/qWUUi2XiqGszwIf\nWF/KOhEZbYxxApcA84GVwCxjzPJklyVRRNzMn1+NMXbOPruYrVvTXSKllEqsVIxWOjvC8bnA3GTH\nT5Z27eCZZ2q57bZCBg4s5cknaznggLQNrFJKqYTKuBnS2SQvD267rZ4rrqjn9NOLeeMNnTCnlMoN\n2lqeAMOHO+nWzc2oUcWsWNHIZZc1YNOVN5RSWUxrDgnSu7e1cN+8eXmMHVtETU26S6SUUi2nySGB\ndt/dwyuv1JCfD0OGlLBunVYflFLZSZNDghUVweTJdQwd2sjJJ5fw4YfaD6GUyj6aHJLAZoPx4xu5\n7746Ro0q4skn85t/klJKZRBNDkk0YICL2bNrePDBfK69tpBGXdhVKZUlNDkk2X77eXj99Rq+/97O\nsGHFbNqk/RBKqcynySEF2rSBGTNq6d3bxcCBJSxfrm+7Uiqz6V0qRRwOuOGGBq67rp6hQ4t57TWd\nYqKUylx6h0qx3/3OyX77ubnggmJWrLBz1VUN2DVFK6UyjN6W0uCQQ6wJcwsXOhg1qoiqxO6lrpRS\nO02TQ5p06ODhxRdr2XVXD4MHl7B2rXZUK6UyhyaHNCoshHvuqWfEiEYGDy5h0SKdMKeUygyaHNLM\nZoPRoxt58ME6xo0rYtq0fFKwrbdSSkWlySFD9OvnYs6cGmbMyOfKKwtpaEh3iZRSrZkmhwyy774e\n5sypYfNmGwMGwIYN2g+hlEoPTQ4ZpqwMpk+v48QTobKyhKVL9UeklEo9vfNkILsdbrkFbr21nuHD\ni3npJZ2OopRKLb3rZLAhQ5x07epm5EhrwtzEiQ04dECTUioFtOaQ4Xr0sCbMffKJgxEjivn113SX\nSCnVGmhyyALt23uYNauWzp3dnHxyCV9/rR3VSqnk0uSQJfLz4e676xk7tpEhQ0pYuFDbl5RSyaPJ\nIcuMGNHII4/UMX58EY8/rjvMKaWSQ5NDFvq//7N2mHvkkXwmTizE6Ux3iZRSuUaTQ5bq2tXaYW7N\nGjvnnFPMtm3pLpFSKpdocshibdrA00/XcsABVkf1mjXaUa2USgxNDlkuLw/uuKOpo/o//9GOaqXU\nztPkkCNGjmzkoYfqGDu2iCee0I5qpdTO0eSQQ445xuqofuihfG64QTuqlVItp8khx/g6qr/6ys65\n5+qMaqVUy2hyyEFt28Izz9TStaubQYNK+OYb7ahWSsVHk0OOysuDu+6qZ8yYRk45pYT339eOaqVU\n7OJKDiLSTkT0Y2gWueCCRqZOrWPMmCKeeko7qpVSsYmYHESkl4i86Pf4aWA9sF5E+iajMCJykIg8\nKCLPi8iYZMRojY491uqonjKlgBtvLMTlSneJlFKZLlrN4X7gCQARORY4GugIDAD+EmsAEZkuIhtE\nZFnQ8UoRMSKyWkSuBTDGrDTGjAPOAgbG91JUNPvt5+H116tZscLOyJHFVFWlu0RKqUwWLTnYjTGv\ner8eAjxnjNlujFkJxNO09DhQ6X9ARBzAFOBkoDtwtoh0937vVGAu8FwcMVQM2rWD556rpUMHN6ee\nWsL69dpCqJQKL1pyaPT7uj+wIMbnBTDGLAQ2Bx3uC6w2xqwxxjRgJYLTvOe/aoypBEbGGkPFLj8f\n7rmnnjPOcDJoUAlffKFjEpRSoaJtE1orIqcBbYC9gXfB6hcAdnboy57A936P1wFHisjxwO+AIgKT\nkUogmw0mTGhg333dDBtWzL331nHeeekulVIqk0RLDpcBU4FdgHOMMY0iUgwsBIYlozDGmAW0IClU\nVJQnvCytIdaoUdCjB5xxRgmbN8Oll+bOa2sNsVIdT2NlV6ydjRcxORhjvgZ+G3SsVkS6GWO2tjii\n5Qegs9/jvbzHWmTjxu07WZzYVFSU51ysrl1h9mwbI0aU8cUXDdx+ez2OJE+JyMX3MdWxUh1PY2VX\nrFjiNZc4og1lvSjK956KpXBRLAa6iUgXESkAhgOvNvMclSR77+3hv/+Fr76yM2KEjmRSSkXvWK4U\nkTdEpJPvgHck0afA8lgDiMizwAfWl7JOREYbY5zAJcB8YCUwyxgT8zVV4rVrB88+W0vHjm6GDNGR\nTEq1dtGalU4VkXOABSIyCTgW6AJUGmNMrAGMMWdHOD4Xa8iqyhC+kUyTJxcwaFAJM2bU0quXO93F\nUkqlQbQOaYwxz4jIj8AbgAGONMZUp6RkKi38RzKddZY1kmngQJ1SrVRrE63PwS4i1wEPACdhTWb7\nSET6pahsKo2GDHHy1FO1XHVVEQ8/nI/Hk+4SZa5333WwYUP0ZriqKvjxR22qU9kjWp/DR8B+QF9j\nzAJjzN+xOo7vFZF/paR0Kq1693YzZ04NTz6Zz403FuLWFqawzjqrhL/+tSDqOZdeWsQhh5SlqERK\n7bxoyeEOY8xoY8yOsVDGmGVYayylbjyWSqu99/Ywe3YNn39uZ+zYIurr012izNRc4vzlF601qOwS\nrUP63xGONwDXJa1E8SovpyKFYy8rUhYptbGixavAGm4GQNjfilDu0jJqrp5I7UUTdr5gWcDt1pu/\nyi3Zv7CODsrPSPbqKkr+dle6i5EyzdUcbJo7VJbJ/uRQpu24mcpe3XoSt/bHqFwTdSirj4i0BXbF\nb6luY8yaZBUqLtu35+T090ybah/sqafyufvuAp58spbDDgu8M1Z0aJPo4mW85kZzac1BZZtmaw4i\ncj/Wqqlv+/17K8nlUhnuvPMa+fvf6zj33GI+/FD3p96ZmsNHHzm44IKixBVGqQSIpebQH6gwxtQl\nuzAqu1RWuigurmPUqCKmTq3juONa72S5nak5vPZaHnPn5gP6J6YyRyx9Dqs0MahIjjvOxfTpdYwf\nX8Sbb7beGoROElS5JpaawzoRWQj8B3D6DhpjbkpaqVRWOeooF08+Wcv55xczaVI9f0h3gdLAPzm4\nXGC3B9YWtM9BZZtYag6bsPoZ6gGX3z+ldujd283MmbVce21huouSdiJl3Hhj7O+D1jpUJmq25mCM\nuVVESgEBPNYhU5P0kqms07Onm+eeq4UB6S5J6vnf4H/91cann7beJjaVG2IZrXQ6sBp4EHgE+EpE\nTk52wVR2Ovjg1jngX4eyqlwTS7PS1UAvY0xfY0wfoC9wY3KLpXLFokWt4xN0cHLQpiKV7WJJDg3G\nmI2+B8aY9Vj9D0o1a+zYolYxD0KTg8o1sYxWqhKRK4E3vY8Hoquyqhg98IA1D+KZZ2o59NDcbXLS\n5KByTSzJYTRwG3AeVof0h95jSjXr98NK+T3AbwOP+68A29pWcFUqG8QyWmkDMC4FZVE5wl1aFtei\ne74VXLM5OSSj5tClSxm33FLPyJGNO38xpeIUbZvQmd7/vxeR7/z+fS8i36WuiCrb1Fw9EXdpfKvl\nZvsKrsloRqqutvHJJ7nfX6MyU7Saw6Xe/49JRUFU7qi9aELEWsDUqfk8/XQRL71URYcOnlazgmu0\noazaP6EyUcSagzHmZ++XNqCzMeZbrJbjm4CSFJRN5aDx4xs5+2wYNqyYLVvSXZrUW7w48gDBoUOL\nqQtaxUwTh0qXWIayPgY0iMhhwBjgReD+pJZK5bSbb4Zjj3Vxzjm58xkj1m1CBw8uDXi8dKmdF17I\nB2Dhwjw2bdLZciozxJIcPMaY/wFnAJONMXPx2/RHqXjZbHDrrfV07567S3TFOiP6hhsK2bIl9OS9\n97b6bLTmoNIlluRQJiJHAEOBeSJSCOyS3GKpXGezwaRJzc+lXLvWxpIlmb+bbaKXz6irs7XoeUol\nSizzHO7BWlPpIWPMRhG5C3gmucVSrYEjaCBOuM7pCqDKVsbHg6+nx/SLU1OwFmjpJ3xdk0llqmY/\nkhljZgKHGWPu89YaHjDG3JP8oqnWIJYhr2WeKnq/difV1SkoUIbR5KDSJZZVWScCl4lICfAp8IKI\n3Jb0kqlWIdY5EeVU8eKL+SkoUcskqwZgs2mng0qPWBpzhwD3Ab8HZhtjjkTnPqgEqb1oApu+Wc/G\nDb+G/efv5ZdjaQVND3eClo3SDmiVKWJJDo3GGA9wMvCK95hO21Qp98UXDn75JTvbWVpac2hoyM7X\nq7JfLMlhq4jMAQ4yxnwgIqcAubu8pspY/fs7ef31zKw9JOoTf3AS0T4HlS6xJIdzsEYrneh9XA+M\nTFqJlIpg8GAnr72W/uTw2mt5PPNMYDni6XNwxTG9w2aDqiqorIxtwqDHA198kflDf1Xmi7bwnm8r\n0LOAXYEhIjIK6ExTolAqZU480cnixQ62bk1vOa68sojLLy+Oek60T/x77FFOQ0Nssex2WL/ezpIl\nsbXkLljg4IQTSkOOacJQ8Yr2MawX8DrQL8z3PMD0ZBTIu2f1YKANMM0Y80Yy4qjsU1YG/fo5mTcv\nj+HDnWkrhzWCKPDuH2+zUmMjFBTEEiv69086qYTjjnNyww1WtqkPM69w2LAS9tnHzeLFrXAssGqx\naMnhdQBjzB8ARKS9MWZTS4KIyHTgFGCDMeZgv+OVWCOhHMCjxpi7jTGvAK+IyC7A3wFNDgqwJsnN\nBes389LQ76dq0yB7mA/hX3+dnk/mS5c6cLnYkRyUSpRov9H3Bj1+fifiPA5U+h8QEQcwBWsUVHfg\nbBHp7nfKDd7vq1Ysnn0hfJsGBYu1CWdn/PBD9OQQXAOIVNMIPu4/z+G995pvWvJ4tAdbJUa03+jg\n37IW/9YZYxYCm4MO9wVWG2PWGGMagOeA00TEJiJ/BV43xixpaUyVG+LdOCh406Bt22CvvcpZty75\nN8145zps29b8zfx//2tKCL//fQnnnVeMM4YWtTVrbGGbmJSKVbTkEPzZJtHTc/YEvvd7vM57bAJW\nh/dQEdHtSVu5cJPkJv+rhhNPaIw4Wc6fMdb/q1YlrtknUhKI5abt89NPNrp1Kw/7veOOaxqZtHq1\nI6DW8cYbedTUBJ6/bJmDgw4K7IQ+6qgy7r8/hk4NpSJI/7jAIMaY+4lzv4iKivB/ZMmQq7FSHW9n\nYo0aZe1YgIvuAAAgAElEQVQJUVtbzt57R7/2G94eq9raEioq4ovj8YTvEPY1/ZSUlFPqd092Opti\nOxyOgHIUFgZeo6DAqg3l5wc2FbVvX8bKlYHn7rJLadA55bRrF3jOpk12KirKaeO3dqHLVUhFhRU4\nL8++0z/fbPn90FiJiRctOfwmaK/oDt7HNqw9HsL8WcblB6xhsT57eY/FbePG7TtZlNhUVJTnZKxU\nx0tErNNPL+Rf//Lw5z9bHQr+933/a69aZf1xrF9fx8aNjTFff9EiB2eeWcKGDaHldLvLABtlZXi/\nb8VwOn2xy3G5XGzcaH3Eb2iAefOs5/hs3lwNlNLY6MJ/wYHNm6uAwGY037n+r6+xMfQPf+PG7Wzd\nmgdYw2xrahrYuLEeKOfrr2Hlyip2261lDQDZ9vvR2mPFEq+5xBEtOUgLyxSrxUA3EemClRSGY024\nU6pZ553XyLnnFnP55Q1Rh4SuWgW77+5m+/b4+hyi9VFE6kz2b1ZatcrOuecW8/TTtcydmxeyU5xv\nv4ZYxNqZ3ZwffrC1ODmo1idicvDuGZ0QIvIscDywm4isA242xkwTkUuA+VgfnaYbY5YnKqbKbQcf\n7Gb//d289FL0OQ9ffw29ern59df4kkO0+QXR+hx8z/v1VxtvvpnHmjU2LrwwdMLcqaeGn/H8/POJ\nW3n2wQcL6Ns3d3fbU8mVkj4HY8zZEY7PBWvoulLxuuyyBiZOLGTYsOjJYeRIFxs3xpccoo088v+e\n/6d4pxNqawPPPeqo2EdaAfzlL4XNnvPppw4GDAi96d92WwGHHx5Y8IsuKgp4/O67DqqrbZxySvom\nEarsoHPqVdbq189FmzbwwgvhP+P8+qt1s+7a1U1VVbzJIbZmpeDkMHFiUegTEmz48BJeeSWPjz4K\nPD55cmhiCW6+uvDCYkaNir70h1IQY81BRPoBR2ANZ/3QGPNBUkulVAxsNrjllnrGjSsi3Aaiq1fb\n2X9/a9mN6urk1Bzcbmuimt0OTqctZJhpslx4YTGHHBJ6fNs2nQSnEiOWneBuA/4G7IE1D+F+7+5w\nSqXdkUe6OOyw8O3qS5c66N0bSks9cd+0oyUH/9qC2w15ebDPPp645jmEu1a8li4NPXbFFdFrLrqZ\nkIpVLDWH/sBvjDFuABHJAxYCoesUKJUGd9xRD6+FHl+yxMHxx1vJIZE1h+DkYLdDXp6HxthHyiqV\n8WLpc7D7EgOAMcaJbvajMkinTqEfh51OeOstB5WVUFIC1XEuSBrtE7b/fgy+5OBwxDdDOh7Juq7P\no4/m88kn2v2oAsVSc1giIq8Cb3kfn4Q1R0GpjLR+vY3XX8/jwAPd7LuvnU2b4q85REsO/p3VvlnU\n+fnJu4mvX5+YfoTJk0MnhDz2WD7XXVfEgAFOnnuuNsyzVGsVS3K4DBgGHInVIf0kO7dCq1JJ9X//\nV0qbNh5mzaoF8igtja9DeuLEQkpKYmuctzqkrX6HZCWHRPUT/Pvf+bRpE3ixa64pSmgMlTtiSQ4T\njTF3Yq2aqlTG++qrKhyOpn0XrD6H2J8/bVoBBxwQ2+Qxl8tqUmpps9LHH8eyDHf8143lWlu2NH3t\ncllzIu67r478xM3DU1kslobGg0Rk/6SXRKkEyc8P3JCnsNC6+cXTYRzr8tsulw2Hw+qQTlbNId6l\nwGMl0rS2zsKFebzwQj4bNuhQWGWJpebQC1gpIpuABhK38J5SKWGzQWkp1NRA27aJvbZVc/BkRbMS\nsGONqe+/j5wE/vtfB8uW2Rk7VodftWaxJIchSS+FUknmG87atm18d9pIy3b7BI9WSkbbfTJ2d+vd\nO/yyHh6PtYTH4sUOTQ6tXCzNSqXAOGPMt97F+G4heE1hpTJcrHMdfDd334gkVzNdD03zHKCyMr7V\nVmOVrs7iL77Q4a2tWSw//SkELo43HXggOcVRKjlinevg21rTt4Bec8nB1yGdl2fdwX/+OfuTgy/e\nlCm6k1xrFktyyDPGLPI98P9aqWwRa83BlxR85zbXGew/WimW81silclh6VLHjhFU9fVN78eoUbBp\nk3ZWtyax9DlsE5HxwAKsZFIJpG47I6XiVNGhTeBj4H2AM2J4Lt7N0n3bUu8D7tIyaq6eSO1FE0LO\n929WgmT1OST+mpH84Q9NK7bOmZPPCSfYef/9Gh57DAYMsDNwoO4P0VrEUnP4A9AbmAU8C3TzHlMq\nY7hLk9cNZq+uouRv4ZcS8w1ldTQ/XaHFPvkkiRdvxurV6Yut0qvZmoMxZiMwJgVlUarFaq6eSMnf\n7sJeXZWU6/uuG/wp3jeU9aWXrJljyWhWMkY7hlXqRUwOIjLTGHOWiHyPt6btT+c5qExSe9GEsM0+\nvk3Wr7++kH32cXPhhdGHZ378sZ1Bg0p3PPYQ2M4ePJHO1+dQWdnIvHn5uFzZ3yEdyS+/2AFtVmot\notUcLvX+f0wqCqJUMsXaId3cjnG+0Uw+vj6Ho45yMW9eflImwr3zTkp2823WFVdYC/TtsUeGZCuV\nVNF+60REJMr3v010YZRKltJS2B5mGMWXX9o58MCmtqAtW5pLDoHf99UcCgqaHueaN95o6neoq0tj\nQVRKRUsOC4Avgf9h7d/g/1fhwdrwR6msUFLi4aefQtvujz22lFWrtu9YVmPLFht2uyfiHtINDYGP\ng4eyJnvvhXQ477ySdBdBpUG05HAMcB5wLPAG8JQxZklKSqVUgoVrVvL1H2zb1rSsxpYtNjp29PDj\nj7EnB/+hrLmYHPydfXYJH35YzV13FbBgQR7z56do02yVchGTgzHmfeB977agg4CJIrIf8ALwtHcp\nDaWywi67wC+/BN7wAye8Wclh61Ybu+/u4ccfw1+noSHwGk6nDYfDg8NhPT8ZHdKZZM0aO3/4QxFz\n5ui63rkulqGsTuBV4FURGQj8E/gTsFuSy6ZUwvTo4WLZssKAY7W11o3cf1mNzZttdOzoBkLH91d0\naNM0Sc7ndDgN4H/WrlhsCXla7pnj93WHll0i2sRClRmaHUAtIvuKyE0ishwYB9wIdEp6yZRKoM6d\nPdTV2QLWPvLVHGpqmo5t3Wo1KwHY7R5cJbrGZDJEm1ioMkO0eQ5jgPO95zwF9DPGbE5VwZRKJJsN\nevZ0sWyZnY4drSFFvhVUa/yazTdvtnH44VZyaNfOw3dnX8c+j/8laZPrWjN9TzNbtJrDw8DuWBv8\nDANeEJF3fP9SUjqlEqhnTzdffNHUXBSu5rBli40997SGtrZrB+vOupRN36zHhofzz6vnxReqcdjd\n2PBgw8P0aTUM6N/I9Gk12PBgtzV9r2I3146vc/XfHrtbr3HshfVs3PBrTP9UdojW59AlZaVQKgV6\n9XLx2mtNv/JNNYem5LBtG/Tr5+L++2uZOrUgYN6Cx2ON8y8qaqptNDYGDmX135gnLzPmriVVuOHB\nKjdEG62ko5FUTunZ081f/hJac/DvkK6utrHLLh6GD3fy0EMFAWslbd9u47zzSthjD/eOhOJLDr79\nHPzl+w3oOfxwF0uW5O4idmvWaJLINfoTVa1G165u6upg5Urr1953g/f973Ra8xiKvatW2+2BC+n9\n+KP1PP8hsU6nNWku3KqsyVypNdO89VYrqCa1MpocVKvhcMAf/9jI3XcX4PFYk9+gqYmopsbaMc63\nZ7TDEbgcRpW3/7RLl6aMEdys5K+oqKk24b943tixDaEn54AanQ+XUzQ5qFZl3LgGvv3WzpQp+WzZ\nYg1X9a9BlJY23cUdjsCaQ1WVjT59XPzjH00LDDmd1nnHHuuiZ8/AWP59Dv7JwZ6jf3UPP6zbiuaS\nHP01VSq8wkJ4+ulaHn+8gClTCujZ072jz6G62lqgz8dmC5zxXFVl9Uf41xJ8NQe7HQ49tOn4M8/U\nREwOthydRL1+fY6+sFZKk4Nqdfbc08OcOTWMH9/I+ec37qg5VFcH1xw8ATf17duhrCw4Odh21ARe\nfLHp+IknugKSQGDNITeXvH78ca055BJNDqpV6tjRw6WXNrDvvu4di+xVV9soKWm6cdvtgX0OTqct\nTHJo6m945hnr/2nTrGFQkZJDrtYcVG7JqCEGItIVuB5oa4wZmu7yqNzXu7eL5csdbNliDW0tKmr6\nXnCHNAR2WIOVHHzDWE87DQ46yIWI1VHh31+hyUFlm6QnBxGZDpwCbDDGHOx3vBK4D2uFs0eNMXcb\nY9YAo0XkhWSXSymw+hhOPtnJU08V0LWrm+LiwJpD8J7QwTWH+npbQN/Ce+81DdmJlBxytUNa5ZZU\n/Jo+DlT6HxARBzAFOBnoDpwtIt1TUBalQowZ08CMGfnU1jbNcYDYkkNDQ+BkN3/+CcHjgbfftnq+\nteagskHSk4MxZiEQvGBfX2C1MWaNMaYBeA7vysdKpdqhh7opLfXw7rt5Ic1KockhsEO5vj7yMhnB\nNYf27a3naXJQ2SBdfQ57At/7PV4HHCki7YE7gcNEZKIxJqY1fSsqypNQxNYVK9XxMi3W6afDAw84\nOOccqKiwqgJFRVBWlkdFRdN5nToVBTy22wto2xYqKgpCYvk3H+XlOdhtN2v579LSQpzO3Fx7Kd6f\na6TzM+33Ixtj7Wy8jPr1NMZswtozIi4bN4bZOT4JKirKczJWquNlYqxevRxs3VqCx9PAxo31ALhc\nRWze7GTjRidg/ZG53TVs3eoGrBv9r782suuubjZubAiJVVdXgm/TIKfTxebNtUAZtbX1bN7csOOa\nuSSW99ovt4Y9PxN/P7ItVizxmksc6eoa+wHo7Pd4L+8xpdKiZ0+rDaj5DunAGkFdXeQ+h+3bA9uP\nfM1J2qykskG6ksNioJuIdBGRAmA41lakSqVFRYWVFKz9pC2xdUjbIiaHQYOcYY9rclDZIOnJQUSe\nBT6wvpR1IjLauy/1JcB8YCUwyxizPNllUSqa3r1d9OnTNLEhfId08FDW8Mt1A9x5Z/2Or12upqQQ\nz1DWu+6qa/4kpZIg6X0OxpizIxyfC8xNdnylYvX664HLigbPkAarWcn/k399feRmpeBlNpprVurW\nzcWqVdaTzjyzkdmz86isdDJxYjyvQqnE0Ok4SkUQLjmUloY2K8Uy6qixMXqfw/LlVbz7bs2OpTem\nTq1j3TrdY1mljyYHpSII16xUWmotyOdTXw8FMaw3558cwjUrtW3roaDAatryp/0TKl00OSgVQbgO\naZst8OZuLZ/R/CqrTifYbJHP0ySgMo0mB6UicDg8uFy2gGUwIDA5RFs+w5/LZfOrOcS+ZHe2JY3g\nZjiVvTQ5KBWBr+YQfMPz72OI1iHt784763bc6MOd7/tetESUDXr3LuWbb7Iso6mwsuxXT6nU8XVI\nN3i3fPZ1FhcWNp0Ta4f073/v3JEA4lk2I9tqDj/9ZGPhwoxaeEG1kCYHpSJwOKxP8o2NUF7uYcgQ\na1Kb74ZdUODxNivF1kwUS80h1uOZ6k9/auCNN/Koq4P//MfB9u3wm9+U7NiKVWUPTQ5KReCrOWzf\nbqO8PDABfPZZFUcc4Yq6fEawpuQQOZnssYeHyZNrA8qQTUaPbqSuDv7v/0q58MIirruuiNWrHQwY\nUMpDD8X4RqmMkGW/ekqljm8nuC1bbLRtG3hD79TJQ3ExeDyRl8/wue++wG1DozUr2e0wbJjT73F2\n7Tfdvr2HF1+spbLSSX29jZkzrTfnm2/sPP+8JodsoslBqQh8HdI//WTbsfaSv112sY5FG8r65pvV\n/P73gWssxdOslG01B58bbqjnqadqGTCg6bVv25ZlbWStnPYcKRWB1axk45NPHBx+eOgYTV9y8O+g\nDnbIIU0TJXy7zMUyac6/DNmouBiOPtrF/vvX0aOHtbz5t9/a+fFHW8CS3SpzZemvnlLJ53B48Hhg\n8WJHwIJ8Pnvuad34Y73Z+5bdCFfTyJUO6WAVFR42bLD2FPjNb5y88YZ+Hs0WmhyUisBut4axfvaZ\ng9693SHfL/fulVJYGF+/QLSaRrgy5IING7YzcmQjV19d1PzJKiPkyK+eUomXlwcrVtipqPDs2P/Z\nny8pxNNMNH9+Nf36hdZCcq3PIZzKSifXXlsfcGzSpAK+/DKHXmQO0Z+KUhEUFVlNSr16hV8Tosj7\nITiemsBhh7njuuHnUnIoLrbmQfj7+98LefJJHcWUiXLoV0+pxCoq8rBhg53u3UOblKBpvkK8zUrx\nyKXkEMkjjxTgcsHEiYXMnKl9EplCfxJKReCrGXTqFD45+JqT4mlWildrSA4A115byBNPWG/kxRdb\n80v8981QqddKfvWUil9RkVUj2G238DUD32S2ZCaH4L6I/fYLn6jiETzbO53WrdvORRc17EgMYCXE\nPfYoZ+lSOx06lKexdK2bJgelIvDNSwg3AQ6aVlBN5qf74OTwwQc7v0jR6ac37vQ1EqWgAC6+uCHs\n915+2eqLePTRfOp0K+2U0+SgVATN1RyCNwJSLeNLvnvu6ebpp2tYvRp69HDxwANWbeK664pYulTb\nmFJNk4NSEfhGIYUbxgqhey+ollu6tIolS6o56SQX++0Hzz1nrUc1eXItnTq5eestBx9/bDUz/fyz\njRtvjGOImGoRTQ5KReBryojUp6A1h8TZYw9PQBNax44e/vnPOs4800llpZP77ivk/POtdr4jjihl\n+vR8nM4IF1MJoclBqQj239/NYYdF3veyXz8X06fXRvx+c+65J/EN6Z9/XpXwa6bLuec24nBYi/hN\nn16L223joINc1NXZaGy0cdllRZogkkiTg1IRdO3qYf78mojfLyyEU05p+d3Jt8bSiSe2/BqFhR5m\nzmwqo6+fJJeUlVnv86efVjFrVlMyfv75fK67rpDNm9NYuCyxZQtMmZLPjTfG/n5pclAqTXx9Fu3a\ntfyG3ratJ2CviZYs1Ne9e+TaUSYpKYEOHTzcfXcdJSXWa1661MERR5Rx+unFzJqVR2PmDMTKGLW1\n0L9/KcuXO2hshGOOKY2pz0aTg1Jp4na3fMnV3r1bdkMP3rQIrOazbGGzwahRjSxdWsU332xn/vwa\nli2r4sILG5k5M5++fUuZOjWfqtxpXdtpzz2XT8+ebh54oI67767n1VdrIg7P9qfJQak0+e1vA/ek\nbs6zzzY1Hx1ySFNy8B81FW4E1V//2tS3sTO1lEzSti2UllpfFxfDoEFOXnyxlscfr+XTTx306VPK\nlVcW8vzzeXz3na1VjCzzeMAYe8B+3fX1MGVKARMmNC14uP/+Hi69NPzcEn+6fIZSadKhg3XHijU5\nnHBC+NpCPDc+m81qnmnOvHnVVFaWxn7hDHHIIW4efriOb7+1MX9+HvPm5XHLLYXk58OQIU4uuaSB\njh1zL1O43XDNNYXMmZNHdbWNTp08HHUU5OUV0q2bm759468danJQKs1iSQ677Rb4xx0pIfiuVVnZ\nyGWXNXDyyYE3eI8Hrr22noceir7mx+GHp6apqaJDm/DHd/a6QJ/ggw95/yU4VnPcpWXUXD0Rbr4u\naTGefz6Pzz5z8NFH1RQXw6pVdpYvL+XZZ+38/e8tGxWnyUGpDDdrVg1lZfF92i0uJuwGRdDUHOMT\nvMDdzTcnd60Kd2kZ9urW0ylgr66i5G93JS05eDzw8MMFTJxYv2MDqoMOcnPssTB0aMuHWmufg1Jp\n1lzN4fjjXfTpE3ij32cfd9jn+h43V7Pwlx+0nUKyV0OtuXoi7tKy5AbJMMlMhl99ZWfzZhsDBiR2\n1JnWHJRKs3iHn37zzXaKi+GWW6w1xf0TQUuGsoYbwRRNv35OFi1q+a2j9qIJ1F40IeL3KyrK2bhx\ne4uvH4+KinK++WY7M2fmM3t2HsuWOejf38mppzoZNMi504kyUrNZIi1Y4GDAAGfCF4DUmoNSaWaz\nxXdzLi0NXAm2XTv/azUXK/TY9dfX89FHVYwZ0xDTNa68svmRLtmkrAxGj27klVdq+eADaxvXBx8s\nYL/9yhgzpojnn8/L6KVSFizI4/jjEz9XRWsOSmW5bt3cfP31dvbbrzxss1Jzo5lKSqBLl9gTVEtq\nJ9miosLDiBGNjBjRyC+/2JgzJ49HHy3g8suLGDGikQkTGth9dw9ud9N+Hi3ldltNQh995ODzz+18\n9ZWdXXbx0KWLhy5d3JxwgpPOnaP/XOrr4cMPHTzwQMv7FiLJqOQgIqXAA0ADsMAY83Sai6RU0kVa\nEjwe5eXw3nvVIZ3N/oqLPTz55M7fRPyTw5//XM+kSYX07OkK2R862+22m4eRI61EsWKFneefz+eY\nY0ppbLT6ZXr2dHHEES6OPtrFCSe4mm/WsdlCRkZ1BPrtZDmrACT896KOxGrmU0PSk4OITAdOATYY\nYw72O14J3Ac4gEeNMXcDvwNeMMbMFpGZgCYHldM++6wqIckBrBEq0RxxhIsDDoh8ju+mH0/NoFs3\n63pvvx15DapsZ7NBjx5uevSo56ab6qmttT71L1niYPFiB5MmFTJpEowY0Ui3bm66dHGz664eqqpg\n15IyHDXZOTIrFX0OjwOV/gdExAFMAU4GugNni0h3YC/ge+9p2bHgi1I7oVMnz05tMxruRl5S4uHw\nw0P/fILPjbTi7BFHRP/T8/9+LjcxhWO3W30+5eVw3HEurrqqgTfeqOHiixtYvNjBrbcW0r9/CZ07\nl9G3bxm3cjNVtuwcmZX0moMxZqGI7Bt0uC+w2hizBkBEngNOA9ZhJYjP0M5ypaK6/npo0ya0KWft\n2pZ9UvXd6A87LHoNxOGAww93sWSJ7s4G1vt2+ulOTj+9aXVdj8f3fo6llrHUkrhRWMuW2bnmmiK2\nb4ft220sXlwdtv+juXjNTf5LV5/DnjTVEMBKCkcC9wOTRWQwMDsdBVMqW9xxB2zcGPsypPF+yi8r\n81BVFf5JM2fWUF9v48MPNUGEk8wa1cEHu5k9u4YFCxx07uzZ6Y7xSDKqQ9oYUw38Id7nVVSUJ6E0\nrStWquNprNTFKyuz5kMUFuYFnJ+X5wi4hm/NJd/j7dutmdZ1ddZEucMOg/Xrre9XeD92rlgRezni\nkas/s0TGOuus5MZLV3L4Aejs93gv77EWSeWEmVyMlep4GiuV8crZvr0OKKKhwcnGjbU7jrtcLsCx\n4xoNDYVAQdA1ywAbd90FQ4dux+2GjRubvrvrrnagNKGvO1d/Zpn2+9Fc4khXu/5ioJuIdBGRAmA4\n8GqayqJUTvM1cYi4wx73ufLKeubOrQ44NmhQUzt6SYk1Ycxfr15uNmxI3Q1PpU7Sk4OIPAt8YH0p\n60RktDHGCVwCzAdWArOMMcuTXRalWiOPB777bjs331wf9bzyckLWcHrwQWsRvsLmNw5TOSYVo5XO\njnB8LjA32fGVUlBU1PLnvv12NcccU8q2bYkrj8p8OlxUqVZq991jWzCoZ0/3Ts3FUNlJk4NSOW7P\nPcMngXPOacQY7S9Q4WXUUFalVGKtXbs94ragNhvssktqy6Oyh9YclMphsewXrVQ4mhyUUkqF0OSg\nVCvV2hbNU/HR5KCUUiqEJgelWqnmdohTrZvNo78hSimlgmjNQSmlVAhNDkoppUJoclBKKRVCk4NS\nSqkQmhyUUkqF0OSglFIqhCYHpZRSITQ5KKWUCpGTS3aLSFfgeqCtMWZopGNJjFUKPAA0AAuMMU8n\nKp73+t2BW4BNwNvGmBcSef2gWHsB/wK2AF8ZY+5OVixvvH7AuVi/m92NMb9JYiw7cDvQBvjYGPNE\nEmMd7421HHjOGLMgWbG88UqB94BbjDGvJTHOQcBlQHtgvjHm0WTF8sY7HRiM9TObZox5I4mxknLP\n8Lt+Uu8TQbHifi0ZlxxEZDpwCrDBGHOw3/FK4D7AATwa7SZljFkDjBaRF6IdS1Ys4HfAC8aY2SIy\nE9jxQ09ETOBk4F/GmEUi8ioQNjkkKFYv4EVjzFPe1xJRgt7PRcAi701gcTJjAacBe2El2XVJjuUB\nqoCiFMQCuAaYFe2EBP28VgLjvIl2JhAxOSQo3ivAKyKyC/B3IGxySOLfdlRxxo14n0h0rJa8loxL\nDsDjwGRghu+AiDiAKcBJWH9Yi703RQdwV9DzRxljNqQ51l7AF96vXYmOCTwJ3Cwip2J9Ykva6wP+\nC8wWEV/caHY6nt/7eQ4wOsmvTYD3jTEPef9o3k5irEXGmPdEpCPwD6zaUbJiHQKswEpE0ex0LGPM\nBu/v4UXAI6mI5/36Bu/zUhErHvHEjXafSGgsY8yKeC+eccnBGLNQRPYNOtwXWO3NfojIc8Bpxpi7\nsDJnpsVah/WD/4ygfp0ExrzY+4vwUqRCJCKWiFwB3OC91gvAY8mM5z1nb2CbibKHZYJe2zqsKj1A\nxA2VE/x7sgUoTPLrOh4oBboDtSIy1xgT8voS9bqMMa8Cr3pveC8m+bXZgLuB140xS5IZqyXiiUuU\n+0QSYsWdHLKlQ3pP4Hu/x+u8x8ISkfYi8iBwmIhMjHQsWbGwbthnishUYHaUWC2Nua+IPIz1ieFv\nMVy/xbGAd4DLvK9xbZyxWhIPrBpDxCSUwFgvAQNF5F9Y7fNJiyUivxORh7BqX5OTGcsYc70x5nLg\nGeCRcIkhUbFE5HgRud/7+7ggjjgtigdMAE4EhorIuGTGiuOe0dK48d4nWhyrJa8l42oOiWCM2QSM\na+5YEmNVA39IdCy/668FLkzW9YNiLQXOTEUsv5g3pyhODdGbrhIZ6yWi1PKSFPPxFMRYQMuSQkvj\n3Q/cn6JYSbln+F0/qfeJoFhxv5ZsqTn8AHT2e7yX91i2x0pHzFS/vlx9bRor++Kl42871XETFitb\nag6LgW4i0gXrhQ7H6rDM9ljpiJnq15err01jZV+8dPxtpzpuwmJlXM1BRJ4FPrC+lHUiMtoY4wQu\nAeYDK4FZxpjl2RQrHTFT/fpy9bVpLP39yMS4yY6lO8EppZQKkXE1B6WUUumnyUEppVQITQ5KKaVC\naHJQSikVQpODUkqpEJoclFJKhdDkoJRSKkS2zJBWKi7e1SoN1iQhf3OMMfEuVpgwInIB1kZNr3j/\nvSwX0sAAAAMlSURBVAsMNH6b1ojIOVhr+3fxrqMV7jozgE+MMfcFHf8KaynnU4E6Y8zxiX4NqnXQ\n5KBy2cZE3xxFxGaM2dmZo48bY27xLq39FTCCwE1rzvUej2Ya8E+sTV18ZfsN4DLG/EVEnsFKEkq1\niCYH1SqJyDbgTqAS2AMYZoz5QkR6AfcA+d5/lxhjPhWRBVjr7vf23tQvxNrg5kfgQ2BvrI2RjjHG\njPTGGA78zhgzLEpRPgKOEpEyY0yViHQAdvFe11fWCcAwrL/XL71xFwLlItLTGOPbMGYEVtJQaqdp\nn4NqrdoAXxhjBgDPAWO8x58GxnlrHBcRuO1llTGmH1AG/AXoDwwCjvN+/1ngtyJS7n18NlG2zfRy\nA/+maVn0s/Hb3lNE+gJnAMcaY44GtgJjvLWX6YAvERV6z5uBUgmgNQeVyyq8n/j9/dkY8z/v1+96\n//8W2N/7qV2AaSLiO7+NWPsjA7zv/b8b8I0x5hcAEZkNHOz95P8KMFxEZgEHAm/FUM4nsZqInsBK\nDqcBp3u/dzywP/Cut0ylQKP3e08AH4nINVh9DP9t4daWSoXQ5KByWXN9Dk6/r21APVAf7jneG7Nv\nS1E7kbcVfQhrD18X8Ewsu7AZYz4XkV1FZACw1Rjzs19yqgdeNcZcEuZ560XkM+C3wPne2EolhDYr\nKeVljNkGrBWRQQAicoCI3BTm1K+BriJSLtY+3qf4XeMzrA3rryC+rU6fxkoqTwcd/y9wsoiUect0\nkYgc7ff9aVi72R0MzIsjnlJRac1B5bJwzUrfGGOibc04ArhfRK7F6pD+U/AJxphNIvI3rGGya4HP\ngRK/U2YApxpjvoujrM8ANwEvB8X6WESmAAtEpA5YT+AopNeAB4FpxhhXHPGUikr3c1CqBURkBFZz\nz1YReQBYa4yZJCI2rM3i7/efu+D3vAuAfY0xtyS5fPtiDZk9PplxVO7SZiWlWqYd8J6ILAL2BB4U\nkcOBT7BGQYUkBj8XiMi9ySqYiFRijcBSqsW05qCUUiqE1hyUUkqF0OSglFIqhCYHpZRSITQ5KKWU\nCqHJQSmlVAhNDkoppUL8Pzlt5uQccjZkAAAAAElFTkSuQmCC\n", "text/plain": [ - "" + "" ] }, "metadata": {}, @@ -1896,7 +1888,7 @@ "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAADUCAYAAACWNDiHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGXJJREFUeJzt3Xm0HFW59/HvAUEmgTDPQgAfkNcLIoogMYREBkUQ1Jfx\nhsj0iiBXERmuKOEiAgoBFZTBdUH0vQyiUbyARqKEIKAMIiL4Y7ohhkQkRAhh0iTn/rF3Q6c9Q59Q\n1UOd32etrNVdXbXr6ZOnntq1q7qqp7e3FzMzq55l2h2AmZmVwwXezKyiXODNzCrKBd7MrKJc4M3M\nKsoF3sysot7U7gDKFBG9wMaSZtVNmwAcKmlcP8tsCVwLzOtvnjzf8cBRwHLA8sB04DhJLyxlrHsA\nD0uaGRHrAjtKumGIbRwHrCvpi0sTQx/tzQDmStqhYfppwJnAZpJmDNLGUZIu7+ezqcDnJd1XRLxV\nFhE9wPHAkaScWwb4FXCapGcGWObzwFnAGEm31322C3AJsCLwJGmbmN1HGwcBJwIr5/X+AfhUX/M2\n+T12BF6W9EBEvBk4QNJVQ2xjP+DDkg5fmhj6aO9WYCtgQ0mL6qYfCnyP9Le7dZA2Bsrzq4AfSPpp\nEfEOhXvwdSJiK2AycOcg8+0JHAPsKmkr4O2kDeBrb2D1nwU2ya/HAPsMZeGIWEbSRUUV9zrr5J1e\nvY8AfRaVhpjWA07q57MeSWNd3Jt2FnAIsFddzj0H3BoRK/azzKWknHq6fmJErApcBxwpaXPg58BB\njQtHxNuBC4GP5nW+DZgB/Ocb+B6fAP4lv34nMH4oC+c8n1xUca/zd2C3hmkHAn9uIqZl6Wfbz/GO\nb0dxh4r34JfCAmBXYG8gBpjvHcBjkp4FkPRKRBwOLAaIiLWAK4BtcpsnSpqSe+bfBTYF3gx8U9Kk\niDgTGAtsHRHfIvW63hQRq0g6MCL2Bb5M2ok8BhwsaW5ETAQ2ALYDrskb7kaSjsy9khuA/YHNSEcY\nB0nqzUcx55A2/AuAKyT19PNdbyZt/P+Rv9v/IRWWtWozRMQ+wFdIRzIvAEdIuh+4A9goIv5E2qgf\nAb5D2qh3zzEeCuwIjJa0T25vCvATSRcP8H8wbETEGsBngO1qR6OSFgInR8RY4F+By/pY9NuSfhcR\nezdM3xe4T9Jdua1z+1n1NsDTtaM0SYsi4lRghRzXiqSdyCjgFeAsSd+PiJVI+b8dKSd+KOnEiPgk\n6f9+n4jYMH+nVSNiuqRREfE+0g5lBDCXlOdP5Hz9MLAacH9EPEg+Co+IK0lHIDuTdkCPAPtKeikf\nFX+HtA1eQCrC2/Zz1FnL81/k7zYCGAk8UZshInYCLiJth4uB4yXdkpdZLef5Xvm7Twc+BhwZEWfl\nOF4GTgPeJWlxRFwGPC/p8/38/d8w9+DrSJolaW4Ts95CKlDfjYi9IuItkuZLWpA/Pwd4SNJI4DDg\n6nw4+kVgVu4NjQXOjoiNc6/7KeCQvLFdBFyfi/tI0mHiQbm9X5EOrWs+BHxQ0qQ+4vww8AFS4o8B\nds7F4lvAOFIPao9BvusPSD2ZmgPzNAAi4k2kndYnJb0N+AlwXv74cGCmpK0k/T1P21jS2xo2sguB\nDSNi97wzewvw7UHiGk7eS/o7PtLHZz8FRve1kKTf9dPetsDciJgcEY9ExDW5U9Lo18AmEXFDROwX\nEWtIelnS3/LnnwOWl7QZKc8uiogNgE+RivTWwPbAhIjYRdIlwG+Bk3KenwrcmYv7W/J3+XdJWwBf\nJx1l1OwOHCPpxD7i/DhwALA5sDawX+5Vfxc4WtLWwJbAKv38PQD+G9gzb6cAHyXlcr3LgQvy9nsO\nr2+HhwOLcp7/T562A7CNpF/XFpb0Q2Amqei/k3TEcPoAMb1hw6HA3xoRf6r9A85+ow3mDed9pL/f\nd4Fn88ZSG2L5IHB13bybSnqVNIZ6bJ7+BPAXUu96IHsCt0p6ML+/hNQDWja//80AO6Xr8wb5Iqln\nswmpt/yIpAclLWbwQvoY8GJEbJfffxT4Yd3fYiGwQd347nRSz6c/NzZOyOOeRwHnkzaco3JslqxB\n/0NiT+fPh2J1UsH8PKmX/ippJ7uEPM7+HmAO8A3gmYi4JSJqQywfBK7J884iHT3OJv0/7itpcd4Z\n/JGBcwLSUcAsSb/I7V0NbFG3TT0q6dF+lr1R0ryci38g5fnbgDdLujnP800GrncvkHZoe+X3B5LO\nxdXbnrxdM3ie39xPDh8LnEza7o6V9NIAbbxhw2GIZte+TrLm11eREhhgrKSnmm1U0j3Av+YTWduT\nTjpeC+xEGr54rm7e2onXHUi99k2ARcD6DL6TXR14f9451TwPrJlfzxtg2efrXi8CliX1rOqXaeY7\nXw0cnHvrM/PwUP3nn4qIw0jDTisAA93gqM94Jd0XEfNJPaEH+5pnGJtLGorry7rAXyPiPUDtZOVk\nSacO0N7zwFRJjwFExNeBn/U1Yz5q+H95vq2BU4CbI2Jj/jnPa0ewWwCT8jmtRcDGpGGLgawObN6Q\n56+SeuSwdHn+t7rpzZwUruX5ncD6ku5vyPMDgePz0cayQH/Dmv3GK2lWRNxF6iD+oomY3pDhUOD7\nJWlIJ3hq8ljhk3lIpxe4NyJO5vWTs3NJyT8jz78pqZB+nzQWeEkeC2+muM4GbpH0sT7iWJrw57Pk\noer6TSxzLXAradzxmoYYdib1SN4jaUZEfIB0KDskEfEhYCGwQkR8UNJNQ22jwu4E1oiIbSX9vuGz\nvUnncn5LuhKkGU+ShixqFuV/S8hHbS9LEoCkhyNdqTWfdNRQy/Pa/BuRCtvFwL3AR/K4/a8b2+7D\nbNJVZDs0fhAR72jye9VrzPP1mljmJlLsBwHXN8SwISmvd8yFf0vSUfGQRMS2pA7h/aQLNUo9zzQc\nhmjKcChwSUSsBq+NQx8ETMuf3wBMyJ+9HbiPtDNdB7g3F/fDSCdrakn4D1IvpvH1z4FReSyeiHhP\n7nEtrXuBf4mILSJiGdJldwPKRzazgf9Lusqo3jrAX4GZ+eTaYcDK+cjmH8Aq+e/Tr4hYmTTmehzw\naeDiPM0ASc+TrqL5XkRsBinnIuJsUk/ymoGW78OPgdF1hfNo0nmlRnsA3490NVTtsstDSeeX5pLy\nfHxE9OR5fkcq+OsAv8vF/QMsOf7dmOer5nZ/A6wf6TJKImJkRHwvf7Y0HgWWi4hd8/tPMvCRJZJe\nAaaQhq4ah2fWBl4E/pTz+egc5yr5eyyTe/b9ytvbZcAJpOHa0/KOozQu8HUi4pSIeIW0px4TEa9E\nuqKj0WeAPwF3R4RIe/J1SZeAQerRbhTpOvJrSVcDvEw6yTo5Ih4gJfylwOW5eF9PuhLmBFKS7RYR\nd0uaQxqfnhwRD5NOwDYmX9Nye/9OOln7G9JYYjOuBv4o6bmG6T8jFf/Hc9wXkg6ZfwA8QOrR/aVu\nLLUvZwD/LekPuSc6lXTVkGWSziMVh5/mYYyHSL3ocXUnsJcQEQtyPr8VmJrz+f2SZpJydXJEPEoa\n/jmhjya+Stqh/zLn+eOkiwM+nD+/gLRzf5J0hHdibvvLwPmRrnYZTfr/PSMf7U0Gzo2IScDted2z\nSZcpfgz4Zs7zyaRrx5fqfub5nNcxwJURcT9pG13MIEWelOfzJD3UMP33pB7+I6Qjqp8Cd+XvPSd/\nl5n5O/bnU8AcSTfnv9PFpO25ND2+H/zwE+ka9N78ehvgdkkj2hyWWWnyEeECYPV8RDQsuAc/zOTD\ny6dqh8Kky8sG/GGXWTeKiLsj4oD89gDSGP+wKe7gHvywFOmn3meTdvBzSD9Meqy9UZkVK9LtGC4m\n3Y5hPuk6+rvbG1VrucCbmVWUh2jMzCqqo66D77l30DPcTRv5rj8W1RRPTNumsLasvXpHD/jjlFLs\nys8Ky+tpPQNdjDRU1w0+i3WN3t6J/5Tb7sGbmVWUC7yZWUW5wJuZVZQLvJlZRZV+kjUiLiDdz7oX\n+Lfhdh2qVZPz2rpBqT34iBgNbClpJ+AI0j2lzbqa89q6RdlDNGNJd65D0sPAiEiPlTPrZs5r6wpl\nF/j1WPJJNM/Q3H2ZzTqZ89q6QqtPsrb8RyZmLeC8to5UdoGfzZI9mw1IN7cy62bOa+sKZRf4KaSb\n+BMR2wOz655PatatnNfWFUot8JLuID2v9A7SlQbHlrk+s1ZwXlu3KP06eEmnlL0Os1ZzXls38C9Z\nzcwqygXezKyiXODNzCrKBd7MrKI66olOLCiuqZV4qbC2dhr9y8LaunPaboW1Zd1hWs9dhbV1OnsV\n1tYZnF9YW8n8gtuzN8o9eDOzinKBNzOrKBd4M7OKcoE3M6soF3gzs4pygTczqygXeDOzinKBNzOr\nKBd4M7OKcoE3M6soF3gzs4pygTczqygXeDOzinKBNzOrKBd4M7OKcoE3M6soF3gzs4pygTczq6jO\nemRfgR6c9u7C2jpz9ImFtfXq6OULa+u+3+xSWFsALOzQtuw1Z3B6YW2dw+cKawvglEIfAejH/xXB\nPXgzs4pygTczqygXeDOzinKBNzOrKBd4M7OKKv0qmoj4KjAqr+tsST8qe51mZXNeWzcotQcfEWOA\nd0jaCdgTuLDM9Zm1gvPaukXZQzTTgY/n188BK0fEsiWv06xszmvrCqUO0UhaCCzIb48AbpK0qMx1\nmpXNeW3doiW/ZI2IfUkbwu6tWJ9ZKzivrdO14iTrHsAXgD0lPV/2+sxawXlt3aDUAh8RqwFfA8ZJ\nmlfmusxaxXlt3aLsHvwBwFrAdRFRmzZe0syS12tWJue1dYWyT7JeBlxW5jrMWs15bd3Cv2Q1M6so\nF3gzs4pygTczqygXeDOziurp7e0ddKZ8WdgaQE9tmqQnCg9mGoMH0+0KvGvJjyfvUVxjwH/wxcLa\nmrFos8LamjdrncLa6n3rcj3171uR2z09E6uf18B0ziisrVHcVVhbcHOBbXWu3t6JPY3TBr2KJiK+\nAXwCeIbXN4JeYGSh0Zm1mHPbqq6ZyyTHAGtLeqXsYMxazLltldbMGPyj3gCsopzbVmnN9OBnRcRt\nwO3AwtpESV8qLSqz1nBuW6U1U+CfBaaWHYhZGzi3rdIGLfCSijs1btZBnNtWdf0W+IiYDv1ftijp\n/aVEZFYy57YNFwP14E9rWRRmreXctmGh3wIvaVorAzFrFee2DRe+VYGZWUW5wJuZVVRTD/yIiBHA\nlqQTU5I0v9SozFrEuW1VNmgPPiI+CzxGuk3WN4HHI+KYsgMzK5tz26qumR78YcDI2pPjc4/nV8C3\nywzMrAWc21ZpzYzB/6W2AQBI+htQ+K2CzdrAuW2V1kwP/vGI+DEwhbRDGAM8GxGHA0j6zxLjMyuT\nc9sqrZkCvxLwN+Dd+f38vNwo0okpbwTWrZzbVmnN3IvmE60IxKzVnNtWdc080enP9HHfDkmblBKR\nWYs4t63qmhmi2aXu9fLAWNKhrS2NzxTX1Ed6diquMWDxs6MKa+usNT5XWFs/f2uRz579QP0b53aB\nRnFOYW317v7ewtrqeaTAR+LOmFhcWy3QzBDNkw2THo2InwOTygnJrDWc21Z1zQzR7NYwaWNg83LC\nMWsd57ZVXTNDNF+se91LutLgk+WEY9ZSzm2rtGaGaMa0IhCzVnNuW9U1M0SzFfAtYAdSL+cu4FhJ\njzWzgohYEXgQOFPSlUsfqlmxnNtWdc3cquAi4HxgfWBD4BKGdq+O04B5Qw/NrHTObau0ZsbgeyTd\nWPd+ckR8upnGcw9pa+DGweY1awPntlVaMz345SNi+9qbiHg3Td5HHjgPOGFpAjNrAee2VVozyXwi\n8F8RsU5+PwcYP9hCETEeuE3SjIh4AyGalca5bZXWTIH/s6StImI1oHcIT7z5EDAyIvYHNgJejYhZ\nkm5Z2mDNCubctkprpsD/f2BM/X2zmyHpgNrriJgIzPAGYB3GuW2V1kyBV0RcBdwB/P21ib5XtnU/\n57ZVWjMF/s3AImDHumlDule2pIlDC8usJZzbVmm+H7wNW85tq7oBC3xE7Cdpcn59LekHIS8DB0t6\ntgXxmZXCuW3DQb/XwUfE8cAZEVHbCWxCujnTPcAXWhCbWSmc2zZcDPRDpwnAOEkL8/tXJE0DJpKe\nWWnWrSbg3LZhYKACv0DSX+ve/xeApH8AL5YalVm5nNs2LAw0Br9K/RtJl9e9Xa2ccGxI7plYaHPL\nrLlqYW31Xl7cI/vefuRDhbWVH9nn3C7Fy4W11DPlssLa6j23p7C2ev6nwMf/AVwysdj2GgzUg38g\nIo5qnBgRJwO/Ki8ks9I5t21YGKgHfzLwk3zfjXvyvDsDc4F9WhCbWVmc2zYs9FvgJT0NvDcixgLb\nkH4Qcp2k6a0KzqwMzm0bLpr5odNUYGoLYjFrKee2VV0z94M3M7Mu5AJvZlZRLvBmZhXlAm9mVlEu\n8GZmFeUCb2ZWUS7wZmYV5QJvZlZRLvBmZhXlAm9mVlEu8GZmFeUCb2ZWUS7wZmYV5QJvZlZRg94u\n2DrYgmKbW+G5CYW11XPhCYW11bt1cY9c4+HimrIyPVVYSz0nX19YW70nFZiLQM+4gh8B2MA9eDOz\ninKBNzOrKBd4M7OKcoE3M6uo0k+yRsQhwEnAQuBLkm4se51mZXNeWzcotQcfEWsCpwO7AHsD+5a5\nPrNWcF5btyi7Bz8OuEXSC8ALwNElr8+sFZzX1hXKLvCbAitFxA3ACGCipKklr9OsbJvivLYuUPZJ\n1h5gTWB/YAJwRUQU+0sBs9ZzXltXKLvAPw3cIWmhpMdJh7Nrl7xOs7I5r60rlF3gpwC7RcQy+cTU\nKsDcktdpVjbntXWFUgu8pKeA64G7gJuBT0taXOY6zcrmvLZuUfp18JIuBS4tez1mreS8tm7gX7Ka\nmVWUC7yZWUW5wJuZVZQLvJlZRbnAm5lVVE9vb7mPjBqKnml0TjDD0QrFNTV1x50La+u2njsLa2ti\nb2/Lf3Ha0zPReV0ZXyi0tZtYvrC29uojt92DNzOrKBd4M7OKcoE3M6soF3gzs4pygTczqygXeDOz\ninKBNzOrKBd4M7OKcoE3M6soF3gzs4pygTczqygXeDOzinKBNzOrKBd4M7OKcoE3M6soF3gzs4py\ngTczqygXeDOziuqoR/aZmVlx3IM3M6soF3gzs4pygTczqygXeDOzinKBNzOrKBd4M7OKelO7AxiK\niLgAeC/QC/ybpLvbHBIAEfFVYBTp73m2pB+1OaTXRMSKwIPAmZKubHM4AETEIcBJwELgS5JubHNI\nbefcHppOzGvovNzumh58RIwGtpS0E3AE8I02hwRARIwB3pHj2hO4sM0hNToNmNfuIGoiYk3gdGAX\nYG9g3/ZG1H7O7aXSUXkNnZnbXVPggbHAjwEkPQyMiIhV2xsSANOBj+fXzwErR8SybYznNRGxFbA1\n0Ek95HHALZJekDRH0tHtDqgDOLeHoEPzGjowt7tpiGY94N6698/kafPbE04iaSGwIL89ArhJ0qI2\nhlTvPOA4YEKb46i3KbBSRNwAjAAmSpra3pDazrk9NJ2Y19CBud1NPfhGPe0OoF5E7EvaCI5rdywA\nETEeuE3SjHbH0qAHWBPYn7SBXhERHfV/2QE66u/RSbndwXkNHZjb3dSDn03q1dRsAMxpUyxLiIg9\ngC8Ae0p6vt3xZB8CRkbE/sBGwKsRMUvSLW2O62ngjtw7fDwiXgDWBv7a3rDayrndvE7Na+jA3O6m\nAj8FOAO4NCK2B2ZLeqHNMRERqwFfA8ZJ6piTPpIOqL2OiInAjA7ZCKYAV0bEuaTD2FWAue0Nqe2c\n203q4LyGDsztrinwku6IiHsj4g5gMXBsu2PKDgDWAq6LiNq08ZJmti+kziXpqYi4HrgrT/q0pMXt\njKndnNvV0Im57dsFm5lVVDefZDUzswG4wJuZVZQLvJlZRbnAm5lVlAu8mVlFucCbmVVU11wH320i\nYj3gXGBb4AXgLcAVkr7e4jjeBXyF9Is6SPc5OVXSfYMstzPwF0lPlByidRnndvdwD74E+f4TPwHu\nlLSdpFHAHsBREfHRfuYvI451chxflrS9pO1JG8QNEbHWIIt/AhhZRlzWvZzb3cU/dCpBRIwDzpD0\nvobpy0v6e359JfAqsBVwCLAhcD7wD9JDH46T9FBE3EpK4lsiYlPgdkkb5eVfAjYH1geulDSpYX1f\nAZaVdHLD9EnAS5JOi4heYDlJCyNiAumWpz8ErgCeBD4r6ZfF/GWs2zm3u4t78OXYBrincWJtA6iz\nsqTRkmYBV5ESbgwwCbi4ifVsJGkP4P3AafmBA/XeCfy2j+XuBLbvr1FJk4H7gc9VfQOwIXNudxGP\nwZdjEXV/24g4GjgYWAH4s6TaQxTuyJ+vDqxb95i2W4FrmljPFABJz0XEI8CWwLN1n79I/zvxYX3/\nF1tqzu0u4h58OR4Adqq9kXSZpF2BU0iHnDW1Xk/jOFlP3bT6z5ZvmK/+/6+Hf25niTjqvJu+ez+N\n7Zs1cm53ERf4Eki6DXg2Ik6tTYuI5YDdgZf7mP95YE5E7JgnjeP1O9LNBzbOr3drWHRMbnsEsAWg\nhs8vBj6en61Zi2Nn0gMJalc81Lc/pm7ZxcByA35RG3ac293FQzTl2Qf4SkTcT0q0lUnPuDy4n/nH\nA5MiYhHpMPiYPP0i4JKIOBj4WcMy8yJiMulk1OmSnqv/UNKzEbEr8I2IOI/UC3oa2K/u4Q3nAFMi\n4lHg97y+QfyCdH/yz0j60dC/vlWYc7tL+CqaLpWvNLhd0nfaHYtZkZzbxfEQjZlZRbkHb2ZWUe7B\nm5lVlAu8mVlFucCbmVWUC7yZWUW5wJuZVZQLvJlZRf0vhqRibo7514YAAAAASUVORK5CYII=\n", "text/plain": [ - "" + "" ] }, "metadata": {}, @@ -1923,6 +1915,15 @@ "# Show the plot on screen\n", "plt.show()" ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], + "source": [] } ], "metadata": { diff --git a/docs/source/pythonapi/examples/mgxs-part-iii.ipynb b/docs/source/pythonapi/examples/mgxs-part-iii.ipynb index b807ab4a9..842c334a2 100644 --- a/docs/source/pythonapi/examples/mgxs-part-iii.ipynb +++ b/docs/source/pythonapi/examples/mgxs-part-iii.ipynb @@ -32,7 +32,7 @@ "name": "stderr", "output_type": "stream", "text": [ - "/home/wboyd/anaconda2/lib/python2.7/site-packages/matplotlib/__init__.py:1350: UserWarning: This call to matplotlib.use() has no effect\n", + "/usr/local/lib/python2.7/dist-packages/matplotlib-1.5.1+1178.ga40c9ec-py2.7-linux-x86_64.egg/matplotlib/__init__.py:1362: 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", @@ -459,7 +459,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAD6AgMAAAD1grKuAAAABGdBTUEAALGPC/xhBQAAACBjSFJN\nAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEX///9yEhLpgJFNv8Tq\nQYT7AAAAAWJLR0QAiAUdSAAAAAd0SU1FB+AFCwA1IrIgOgEAAAWFSURBVGje7Zs7cttADIZ9CSvX\ncrP0iCxUqbBc8Ag6xR6BhV2EvYvwFD4CCx1ABT1jMdgndpegRQnOrCbjpPlGESISC4A/gd27e8H5\n83CX3b4+iKJrRHkS4vkghMPBonRYWGwtfgD2YN+dRDUOoh6lACw0Noi9w2fESuEoAR/uVuMolX03\n9oXGT7F3eFL2iEfhUX1f4cPdL/ishs+68ai+udE4xPhexbjX2FfjGNoPj/DPNX4Tsd+EODr8FvsV\ndf1Hd9P2VvCi4+s/aXvrf+upAD+1/9GV1mkOH5X9vV6THtfvACslcaUCbESL61drBPtdI8SrFMWr\nELsXCkuFDYW75gbiP7d9Cf7bAYI/aCwUShrBvh30+lWQkzVgZ/HD4OixNCgcQpJ3BxU/Ln91elKo\nM5VEE38QtJ+Yv6cQ9xjKNYayyl8TypP8DfJnQ2H/b/N3ye9P83cT33SQv/sQh9gV7zZ/0dNj5HQa\nC5vVzv9+/WFN2w8KVaZ2BwL1+pv4g0x1QRfjq0dB4Q3kT277oP6VNL6gKxNU9a8zK+WLbi/Wwpdi\nhbboKqyxFOulHMj6v4W/AXbmUeAxrv9J/CqEBXaRKsXaodD4nsYvkT/G6H1D4SR/iPy1Roj9JsQ5\ne18/7EUHv1+Fvx/Xj5V9Ugb5K8TW4TZEEdcvoz/up0VTe9qsVIppKVX6a7D6y9ZvwEKjrtQxPtv6\nfXII9vCxKOGaIeAIfEF8IvAG8ie3vRK9rRQl+PPpSctbhfpTUCpviH+kxsZgpT91+snoX1l49KK3\niUQvICRy5aUw6l8leoVwoo3Uv1rKreF/UFLY6d9QP4L9Wf2r7EP9GOSfcsjZ56f60kz+XmVPXv+R\nuP49ff0T/53Rv6n/7m2lvXT9Wqd/VUz8hvh5M/ED6ILmt4mfHYZSaePnTWpsf/SvqV9O6dLYYClL\nEetnoH/LBLFoBvrX189uTv8++kot5vTvQD4/9jP690g9P/4z/bvo/XVG/xYoZZx+8fr3MxAtsf7t\nUOkG2JqsTtCIpgCt/qX1226KqZS7gfzJbe+c9jLrtIZ8lXD+s4umlW6AKIVrlML2/cXjgPFjlJqI\nRC+Fj0bVJe+vSh56pSdR6YkQ1ygF10Wqf0FeLta/iKn9Mv1L24ti2e+7W4n1b3T/W+L+t9H9T/Sv\nVboUmqJJon1/hZq8LnzRDlDrX1u0xRT1+6vEpomMmyYkqi95vIH8yW1PN+122KkLcNLKi/WTF01z\n/cNASrWE/l3ev6T17zX909z9X27/euK/Rf3zWP+Waf9eEv37KkWJ+rfDl6ZglNDa+cEBhwYDvkoN\nP/rX69814NaI3imq0l7OYDy/qSdDGwr7r+Y3VbzoKZr6XX2lfxfOb87qXzr+b1j/Xlp/nP6dn98M\ncdH7cn7zjPObKsYWS3Eb9w8n85smHtqQuPuZ30T2dlIT6F9xFl+n8xslegL9a4c2KRr9W4rp/GYq\numiM9Nec/j2v/yj9u1h//hv9e93vc++f63/u+rPjL3f+5Lbn1j9m/eXWf+7zh/v8+2b9e/Hzn6s/\nuPqHrb8g71n6L3f+5Lbnvn8w33+4718/+5d47//c/gO7/5E7/nPbc/tv3P4fs//I7X9y+6/fqH+v\n6j9z+9/c/ju3/8+eP+TOn9z23PkXc/7Gnf9x5483q38Xzn+582fu/Js9fy8kb/6fO39y23P3n3S8\n/S/c/Tfc/T83uX/pgv1XE/9duP+Lu/+Mvf8td/znti8kb/8ld/9nx9t/Sjw/Ltr/yt1/+337f6/b\nf0zoB3nJ/ucVc/81d/83e/957vzJbc89/8A8f8E9/5HE78XnT/4H/cs5f8Q9/8Q9f8U+/5U7f3Lb\nc88fdrzzjyvm+cuf/Uu887/c88fs88954/8vO4SjPC+2QRIAAAAldEVYdGRhdGU6Y3JlYXRlADIw\nMTYtMDUtMTBUMjA6NTM6MzQtMDQ6MDDpUcfgAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDE2LTA1LTEw\nVDIwOjUzOjM0LTA0OjAwmAx/XAAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAD6AgMAAAD1grKuAAAABGdBTUEAALGPC/xhBQAAAAFzUkdC\nAK7OHOkAAAAgY0hSTQAAeiYAAICEAAD6AAAAgOgAAHUwAADqYAAAOpgAABdwnLpRPAAAAAxQTFRF\n////chIS6YCRTb/E6kGE+wAAAAFiS0dEAIgFHUgAAAAJcEhZcwAAAEgAAABIAEbJaz4AAAWFSURB\nVGje7Zs7cttADIZ9CSvXcrP0iCxUqbBc8Ag6xR6BhV2EvYvwFD4CCx1ABT1jMdgndpegRQnOrCbj\npPlGESISC4A/gd27e8H583CX3b4+iKJrRHkS4vkghMPBonRYWGwtfgD2YN+dRDUOoh6lACw0Noi9\nw2fESuEoAR/uVuMolX039oXGT7F3eFL2iEfhUX1f4cPdL/ishs+68ai+udE4xPhexbjX2FfjGNoP\nj/DPNX4Tsd+EODr8FvsVdf1Hd9P2VvCi4+s/aXvrf+upAD+1/9GV1mkOH5X9vV6THtfvACslcaUC\nbESL61drBPtdI8SrFMWrELsXCkuFDYW75gbiP7d9Cf7bAYI/aCwUShrBvh30+lWQkzVgZ/HD4Oix\nNCgcQpJ3BxU/Ln91elKoM5VEE38QtJ+Yv6cQ9xjKNYayyl8TypP8DfJnQ2H/b/N3ye9P83cT33SQ\nv/sQh9gV7zZ/0dNj5HQaC5vVzv9+/WFN2w8KVaZ2BwL1+pv4g0x1QRfjq0dB4Q3kT277oP6VNL6g\nKxNU9a8zK+WLbi/WwpdihbboKqyxFOulHMj6v4W/AXbmUeAxrv9J/CqEBXaRKsXaodD4nsYvkT/G\n6H1D4SR/iPy1Roj9JsQ5e18/7EUHv1+Fvx/Xj5V9Ugb5K8TW4TZEEdcvoz/up0VTe9qsVIppKVX6\na7D6y9ZvwEKjrtQxPtv6fXII9vCxKOGaIeAIfEF8IvAG8ie3vRK9rRQl+PPpSctbhfpTUCpviH+k\nxsZgpT91+snoX1l49KK3iUQvICRy5aUw6l8leoVwoo3Uv1rKreF/UFLY6d9QP4L9Wf2r7EP9GOSf\ncsjZ56f60kz+XmVPXv+RuP49ff0T/53Rv6n/7m2lvXT9Wqd/VUz8hvh5M/ED6ILmt4mfHYZSaePn\nTWpsf/SvqV9O6dLYYClLEetnoH/LBLFoBvrX189uTv8++kot5vTvQD4/9jP690g9P/4z/bvo/XVG\n/xYoZZx+8fr3MxAtsf7tUOkG2JqsTtCIpgCt/qX1226KqZS7gfzJbe+c9jLrtIZ8lXD+s4umlW6A\nKIVrlML2/cXjgPFjlJqIRC+Fj0bVJe+vSh56pSdR6YkQ1ygF10Wqf0FeLta/iKn9Mv1L24ti2e+7\nW4n1b3T/W+L+t9H9T/SvVboUmqJJon1/hZq8LnzRDlDrX1u0xRT1+6vEpomMmyYkqi95vIH8yW1P\nN+122KkLcNLKi/WTF01z/cNASrWE/l3ev6T17zX909z9X27/euK/Rf3zWP+Waf9eEv37KkWJ+rfD\nl6ZglNDa+cEBhwYDvkoNP/rX69814NaI3imq0l7OYDy/qSdDGwr7r+Y3VbzoKZr6XX2lfxfOb87q\nXzr+b1j/Xlp/nP6dn98McdH7cn7zjPObKsYWS3Eb9w8n85smHtqQuPuZ30T2dlIT6F9xFl+n8xsl\negL9a4c2KRr9W4rp/GYqumiM9Nec/j2v/yj9u1h//hv9e93vc++f63/u+rPjL3f+5Lbn1j9m/eXW\nf+7zh/v8+2b9e/Hzn6s/uPqHrb8g71n6L3f+5Lbnvn8w33+4718/+5d47//c/gO7/5E7/nPbc/tv\n3P4fs//I7X9y+6/fqH+v6j9z+9/c/ju3/8+eP+TOn9z23PkXc/7Gnf9x5483q38Xzn+582fu/Js9\nfy8kb/6fO39y23P3n3S8/S/c/Tfc/T83uX/pgv1XE/9duP+Lu/+Mvf8td/znti8kb/8ld/9nx9t/\nSjw/Ltr/yt1/+337f6/bf0zoB3nJ/ucVc/81d/83e/957vzJbc89/8A8f8E9/5HE78XnT/4H/cs5\nf8Q9/8Q9f8U+/5U7f3Lbc88fdrzzjyvm+cuf/Uu887/c88fs88954/8vO4SjPC+2QRIAAAAldEVY\ndGRhdGU6Y3JlYXRlADIwMTYtMDUtMTJUMjE6MDQ6MzMtMDQ6MDDlS7KzAAAAJXRFWHRkYXRlOm1v\nZGlmeQAyMDE2LTA1LTEyVDIxOjA0OjMzLTA0OjAwlBYKDwAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] @@ -539,6 +539,7 @@ "\n", "* `TotalXS` (`\"total\"`)\n", "* `TransportXS` (`\"transport\"`)\n", + "* `NuTransportXS` (`\"nu-transport\"`)\n", "* `AbsorptionXS` (`\"absorption\"`)\n", "* `CaptureXS` (`\"capture\"`)\n", "* `FissionXS` (`\"fission\"`)\n", @@ -725,8 +726,8 @@ " Copyright: 2011-2016 Massachusetts Institute of Technology\n", " License: http://openmc.readthedocs.io/en/latest/license.html\n", " Version: 0.7.1\n", - " Git SHA1: 502482dcf630ee6e290c15b8535e6e850a351c88\n", - " Date/Time: 2016-05-10 20:53:35\n", + " Git SHA1: ae588276014a905ecc6e0967bf08288ecec5b550\n", + " Date/Time: 2016-05-12 21:04:33\n", " MPI Processes: 1\n", "\n", " ===========================================================================\n", @@ -813,20 +814,20 @@ "\n", " =======================> TIMING STATISTICS <=======================\n", "\n", - " Total time for initialization = 5.8400E-01 seconds\n", - " Reading cross sections = 1.4100E-01 seconds\n", - " Total time in simulation = 7.7003E+01 seconds\n", - " Time in transport only = 7.6958E+01 seconds\n", - " Time in inactive batches = 6.4820E+00 seconds\n", - " Time in active batches = 7.0521E+01 seconds\n", - " Time synchronizing fission bank = 8.0000E-03 seconds\n", - " Sampling source sites = 6.0000E-03 seconds\n", + " Total time for initialization = 4.5500E-01 seconds\n", + " Reading cross sections = 1.1200E-01 seconds\n", + " Total time in simulation = 5.6386E+01 seconds\n", + " Time in transport only = 5.6351E+01 seconds\n", + " Time in inactive batches = 4.3700E+00 seconds\n", + " Time in active batches = 5.2016E+01 seconds\n", + " Time synchronizing fission bank = 5.0000E-03 seconds\n", + " Sampling source sites = 2.0000E-03 seconds\n", " SEND/RECV source sites = 1.0000E-03 seconds\n", - " Time accumulating tallies = 6.0000E-03 seconds\n", + " Time accumulating tallies = 0.0000E+00 seconds\n", " Total time for finalization = 0.0000E+00 seconds\n", - " Total time elapsed = 7.7616E+01 seconds\n", - " Calculation Rate (inactive) = 3856.83 neutrons/second\n", - " Calculation Rate (active) = 1418.02 neutrons/second\n", + " Total time elapsed = 5.6857E+01 seconds\n", + " Calculation Rate (inactive) = 5720.82 neutrons/second\n", + " Calculation Rate (active) = 1922.49 neutrons/second\n", "\n", " ============================> RESULTS <============================\n", "\n", @@ -952,8 +953,7 @@ "name": "stderr", "output_type": "stream", "text": [ - "/home/wboyd/Documents/NSE-CRPG-Codes/openmc/openmc/tallies.py:1996: RuntimeWarning: invalid value encountered in true_divide\n", - " self_rel_err = data['self']['std. dev.'] / data['self']['mean']\n" + "/usr/local/lib/python2.7/dist-packages/openmc-0.7.1-py2.7.egg/openmc/tallies.py:1988: RuntimeWarning: invalid value encountered in true_divide\n" ] }, { @@ -1301,122 +1301,122 @@ "[ NORMAL ] Computing the eigenvalue...\n", "[ NORMAL ] Iteration 0:\tk_eff = 0.854370\tres = 0.000E+00\n", "[ NORMAL ] Iteration 1:\tk_eff = 0.801922\tres = 1.521E-01\n", - "[ NORMAL ] Iteration 2:\tk_eff = 0.761746\tres = 6.349E-02\n", - "[ NORMAL ] Iteration 3:\tk_eff = 0.732367\tres = 5.029E-02\n", - "[ NORMAL ] Iteration 4:\tk_eff = 0.711075\tres = 3.869E-02\n", - "[ NORMAL ] Iteration 5:\tk_eff = 0.696557\tres = 2.912E-02\n", - "[ NORMAL ] Iteration 6:\tk_eff = 0.687673\tres = 2.044E-02\n", - "[ NORMAL ] Iteration 7:\tk_eff = 0.683470\tres = 1.277E-02\n", - "[ NORMAL ] Iteration 8:\tk_eff = 0.683129\tres = 6.141E-03\n", - "[ NORMAL ] Iteration 9:\tk_eff = 0.685949\tres = 7.889E-04\n", - "[ NORMAL ] Iteration 10:\tk_eff = 0.691329\tres = 4.181E-03\n", - "[ NORMAL ] Iteration 11:\tk_eff = 0.698755\tres = 7.875E-03\n", - "[ NORMAL ] Iteration 12:\tk_eff = 0.707786\tres = 1.077E-02\n", - "[ NORMAL ] Iteration 13:\tk_eff = 0.718050\tres = 1.295E-02\n", - "[ NORMAL ] Iteration 14:\tk_eff = 0.729230\tres = 1.452E-02\n", - "[ NORMAL ] Iteration 15:\tk_eff = 0.741058\tres = 1.559E-02\n", - "[ NORMAL ] Iteration 16:\tk_eff = 0.753310\tres = 1.624E-02\n", - "[ NORMAL ] Iteration 17:\tk_eff = 0.765800\tres = 1.655E-02\n", - "[ NORMAL ] Iteration 18:\tk_eff = 0.778371\tres = 1.660E-02\n", - "[ NORMAL ] Iteration 19:\tk_eff = 0.790897\tres = 1.643E-02\n", - "[ NORMAL ] Iteration 20:\tk_eff = 0.803273\tres = 1.611E-02\n", - "[ NORMAL ] Iteration 21:\tk_eff = 0.815415\tres = 1.566E-02\n", - "[ NORMAL ] Iteration 22:\tk_eff = 0.827256\tres = 1.513E-02\n", - "[ NORMAL ] Iteration 23:\tk_eff = 0.838747\tres = 1.453E-02\n", - "[ NORMAL ] Iteration 24:\tk_eff = 0.849847\tres = 1.390E-02\n", - "[ NORMAL ] Iteration 25:\tk_eff = 0.860527\tres = 1.324E-02\n", - "[ NORMAL ] Iteration 26:\tk_eff = 0.870770\tres = 1.258E-02\n", - "[ NORMAL ] Iteration 27:\tk_eff = 0.880562\tres = 1.191E-02\n", - "[ NORMAL ] Iteration 28:\tk_eff = 0.889897\tres = 1.125E-02\n", - "[ NORMAL ] Iteration 29:\tk_eff = 0.898776\tres = 1.061E-02\n", - "[ NORMAL ] Iteration 30:\tk_eff = 0.907202\tres = 9.986E-03\n", - "[ NORMAL ] Iteration 31:\tk_eff = 0.915181\tres = 9.382E-03\n", - "[ NORMAL ] Iteration 32:\tk_eff = 0.922724\tres = 8.803E-03\n", - "[ NORMAL ] Iteration 33:\tk_eff = 0.929843\tres = 8.249E-03\n", - "[ NORMAL ] Iteration 34:\tk_eff = 0.936550\tres = 7.721E-03\n", - "[ NORMAL ] Iteration 35:\tk_eff = 0.942861\tres = 7.220E-03\n", - "[ NORMAL ] Iteration 36:\tk_eff = 0.948791\tres = 6.744E-03\n", - "[ NORMAL ] Iteration 37:\tk_eff = 0.954357\tres = 6.295E-03\n", - "[ NORMAL ] Iteration 38:\tk_eff = 0.959575\tres = 5.871E-03\n", - "[ NORMAL ] Iteration 39:\tk_eff = 0.964461\tres = 5.472E-03\n", - "[ NORMAL ] Iteration 40:\tk_eff = 0.969033\tres = 5.097E-03\n", - "[ NORMAL ] Iteration 41:\tk_eff = 0.973306\tres = 4.744E-03\n", - "[ NORMAL ] Iteration 42:\tk_eff = 0.977297\tres = 4.414E-03\n", - "[ NORMAL ] Iteration 43:\tk_eff = 0.981021\tres = 4.104E-03\n", - "[ NORMAL ] Iteration 44:\tk_eff = 0.984493\tres = 3.814E-03\n", - "[ NORMAL ] Iteration 45:\tk_eff = 0.987729\tres = 3.543E-03\n", - "[ NORMAL ] Iteration 46:\tk_eff = 0.990742\tres = 3.290E-03\n", - "[ NORMAL ] Iteration 47:\tk_eff = 0.993546\tres = 3.053E-03\n", - "[ NORMAL ] Iteration 48:\tk_eff = 0.996153\tres = 2.833E-03\n", - "[ NORMAL ] Iteration 49:\tk_eff = 0.998577\tres = 2.627E-03\n", - "[ NORMAL ] Iteration 50:\tk_eff = 1.000829\tres = 2.436E-03\n", - "[ NORMAL ] Iteration 51:\tk_eff = 1.002920\tres = 2.257E-03\n", - "[ NORMAL ] Iteration 52:\tk_eff = 1.004860\tres = 2.091E-03\n", - "[ NORMAL ] Iteration 53:\tk_eff = 1.006661\tres = 1.937E-03\n", - "[ NORMAL ] Iteration 54:\tk_eff = 1.008330\tres = 1.793E-03\n", - "[ NORMAL ] Iteration 55:\tk_eff = 1.009877\tres = 1.660E-03\n", - "[ NORMAL ] Iteration 56:\tk_eff = 1.011311\tres = 1.536E-03\n", - "[ NORMAL ] Iteration 57:\tk_eff = 1.012639\tres = 1.421E-03\n", - "[ NORMAL ] Iteration 58:\tk_eff = 1.013868\tres = 1.314E-03\n", - "[ NORMAL ] Iteration 59:\tk_eff = 1.015006\tres = 1.215E-03\n", - "[ NORMAL ] Iteration 60:\tk_eff = 1.016059\tres = 1.124E-03\n", - "[ NORMAL ] Iteration 61:\tk_eff = 1.017033\tres = 1.039E-03\n", - "[ NORMAL ] Iteration 62:\tk_eff = 1.017933\tres = 9.596E-04\n", - "[ NORMAL ] Iteration 63:\tk_eff = 1.018766\tres = 8.865E-04\n", - "[ NORMAL ] Iteration 64:\tk_eff = 1.019535\tres = 8.188E-04\n", - "[ NORMAL ] Iteration 65:\tk_eff = 1.020246\tres = 7.562E-04\n", - "[ NORMAL ] Iteration 66:\tk_eff = 1.020903\tres = 6.981E-04\n", - "[ NORMAL ] Iteration 67:\tk_eff = 1.021509\tres = 6.445E-04\n", - "[ NORMAL ] Iteration 68:\tk_eff = 1.022069\tres = 5.948E-04\n", - "[ NORMAL ] Iteration 69:\tk_eff = 1.022586\tres = 5.489E-04\n", - "[ NORMAL ] Iteration 70:\tk_eff = 1.023063\tres = 5.064E-04\n", - "[ NORMAL ] Iteration 71:\tk_eff = 1.023503\tres = 4.671E-04\n", - "[ NORMAL ] Iteration 72:\tk_eff = 1.023909\tres = 4.308E-04\n", - "[ NORMAL ] Iteration 73:\tk_eff = 1.024284\tres = 3.973E-04\n", - "[ NORMAL ] Iteration 74:\tk_eff = 1.024629\tres = 3.663E-04\n", - "[ NORMAL ] Iteration 75:\tk_eff = 1.024948\tres = 3.377E-04\n", - "[ NORMAL ] Iteration 76:\tk_eff = 1.025241\tres = 3.113E-04\n", - "[ NORMAL ] Iteration 77:\tk_eff = 1.025512\tres = 2.869E-04\n", - "[ NORMAL ] Iteration 78:\tk_eff = 1.025761\tres = 2.644E-04\n", - "[ NORMAL ] Iteration 79:\tk_eff = 1.025991\tres = 2.436E-04\n", - "[ NORMAL ] Iteration 80:\tk_eff = 1.026203\tres = 2.244E-04\n", - "[ NORMAL ] Iteration 81:\tk_eff = 1.026398\tres = 2.067E-04\n", - "[ NORMAL ] Iteration 82:\tk_eff = 1.026578\tres = 1.904E-04\n", - "[ NORMAL ] Iteration 83:\tk_eff = 1.026743\tres = 1.754E-04\n", - "[ NORMAL ] Iteration 84:\tk_eff = 1.026895\tres = 1.615E-04\n", - "[ NORMAL ] Iteration 85:\tk_eff = 1.027036\tres = 1.487E-04\n", - "[ NORMAL ] Iteration 86:\tk_eff = 1.027165\tres = 1.369E-04\n", - "[ NORMAL ] Iteration 87:\tk_eff = 1.027284\tres = 1.260E-04\n", - "[ NORMAL ] Iteration 88:\tk_eff = 1.027393\tres = 1.160E-04\n", - "[ NORMAL ] Iteration 89:\tk_eff = 1.027494\tres = 1.068E-04\n", - "[ NORMAL ] Iteration 90:\tk_eff = 1.027587\tres = 9.825E-05\n", - "[ NORMAL ] Iteration 91:\tk_eff = 1.027672\tres = 9.041E-05\n", - "[ NORMAL ] Iteration 92:\tk_eff = 1.027751\tres = 8.319E-05\n", - "[ NORMAL ] Iteration 93:\tk_eff = 1.027823\tres = 7.654E-05\n", - "[ NORMAL ] Iteration 94:\tk_eff = 1.027889\tres = 7.042E-05\n", - "[ NORMAL ] Iteration 95:\tk_eff = 1.027950\tres = 6.478E-05\n", - "[ NORMAL ] Iteration 96:\tk_eff = 1.028007\tres = 5.959E-05\n", - "[ NORMAL ] Iteration 97:\tk_eff = 1.028058\tres = 5.481E-05\n", - "[ NORMAL ] Iteration 98:\tk_eff = 1.028106\tres = 5.041E-05\n", - "[ NORMAL ] Iteration 99:\tk_eff = 1.028150\tres = 4.636E-05\n", - "[ NORMAL ] Iteration 100:\tk_eff = 1.028190\tres = 4.263E-05\n", - "[ NORMAL ] Iteration 101:\tk_eff = 1.028227\tres = 3.920E-05\n", - "[ NORMAL ] Iteration 102:\tk_eff = 1.028261\tres = 3.604E-05\n", - "[ NORMAL ] Iteration 103:\tk_eff = 1.028292\tres = 3.314E-05\n", - "[ NORMAL ] Iteration 104:\tk_eff = 1.028321\tres = 3.047E-05\n", - "[ NORMAL ] Iteration 105:\tk_eff = 1.028347\tres = 2.801E-05\n", - "[ NORMAL ] Iteration 106:\tk_eff = 1.028371\tres = 2.575E-05\n", - "[ NORMAL ] Iteration 107:\tk_eff = 1.028394\tres = 2.367E-05\n", - "[ NORMAL ] Iteration 108:\tk_eff = 1.028414\tres = 2.175E-05\n", - "[ NORMAL ] Iteration 109:\tk_eff = 1.028433\tres = 1.999E-05\n", - "[ NORMAL ] Iteration 110:\tk_eff = 1.028450\tres = 1.838E-05\n", - "[ NORMAL ] Iteration 111:\tk_eff = 1.028466\tres = 1.689E-05\n", - "[ NORMAL ] Iteration 112:\tk_eff = 1.028481\tres = 1.552E-05\n", - "[ NORMAL ] Iteration 113:\tk_eff = 1.028494\tres = 1.426E-05\n", - "[ NORMAL ] Iteration 114:\tk_eff = 1.028507\tres = 1.310E-05\n", - "[ NORMAL ] Iteration 115:\tk_eff = 1.028518\tres = 1.204E-05\n", - "[ NORMAL ] Iteration 116:\tk_eff = 1.028528\tres = 1.106E-05\n", - "[ NORMAL ] Iteration 117:\tk_eff = 1.028538\tres = 1.017E-05\n" + "[ NORMAL ] Iteration 2:\tk_eff = 0.761745\tres = 6.349E-02\n", + "[ NORMAL ] Iteration 3:\tk_eff = 0.732366\tres = 5.029E-02\n", + "[ NORMAL ] Iteration 4:\tk_eff = 0.711073\tres = 3.869E-02\n", + "[ NORMAL ] Iteration 5:\tk_eff = 0.696554\tres = 2.912E-02\n", + "[ NORMAL ] Iteration 6:\tk_eff = 0.687670\tres = 2.044E-02\n", + "[ NORMAL ] Iteration 7:\tk_eff = 0.683465\tres = 1.277E-02\n", + "[ NORMAL ] Iteration 8:\tk_eff = 0.683124\tres = 6.142E-03\n", + "[ NORMAL ] Iteration 9:\tk_eff = 0.685943\tres = 7.897E-04\n", + "[ NORMAL ] Iteration 10:\tk_eff = 0.691322\tres = 4.180E-03\n", + "[ NORMAL ] Iteration 11:\tk_eff = 0.698747\tres = 7.873E-03\n", + "[ NORMAL ] Iteration 12:\tk_eff = 0.707777\tres = 1.076E-02\n", + "[ NORMAL ] Iteration 13:\tk_eff = 0.718040\tres = 1.295E-02\n", + "[ NORMAL ] Iteration 14:\tk_eff = 0.729218\tres = 1.452E-02\n", + "[ NORMAL ] Iteration 15:\tk_eff = 0.741045\tres = 1.559E-02\n", + "[ NORMAL ] Iteration 16:\tk_eff = 0.753296\tres = 1.624E-02\n", + "[ NORMAL ] Iteration 17:\tk_eff = 0.765785\tres = 1.655E-02\n", + "[ NORMAL ] Iteration 18:\tk_eff = 0.778355\tres = 1.659E-02\n", + "[ NORMAL ] Iteration 19:\tk_eff = 0.790879\tres = 1.643E-02\n", + "[ NORMAL ] Iteration 20:\tk_eff = 0.803254\tres = 1.610E-02\n", + "[ NORMAL ] Iteration 21:\tk_eff = 0.815394\tres = 1.566E-02\n", + "[ NORMAL ] Iteration 22:\tk_eff = 0.827235\tres = 1.513E-02\n", + "[ NORMAL ] Iteration 23:\tk_eff = 0.838724\tres = 1.453E-02\n", + "[ NORMAL ] Iteration 24:\tk_eff = 0.849823\tres = 1.390E-02\n", + "[ NORMAL ] Iteration 25:\tk_eff = 0.860503\tres = 1.324E-02\n", + "[ NORMAL ] Iteration 26:\tk_eff = 0.870744\tres = 1.258E-02\n", + "[ NORMAL ] Iteration 27:\tk_eff = 0.880535\tres = 1.191E-02\n", + "[ NORMAL ] Iteration 28:\tk_eff = 0.889870\tres = 1.125E-02\n", + "[ NORMAL ] Iteration 29:\tk_eff = 0.898748\tres = 1.061E-02\n", + "[ NORMAL ] Iteration 30:\tk_eff = 0.907172\tres = 9.985E-03\n", + "[ NORMAL ] Iteration 31:\tk_eff = 0.915151\tres = 9.382E-03\n", + "[ NORMAL ] Iteration 32:\tk_eff = 0.922693\tres = 8.802E-03\n", + "[ NORMAL ] Iteration 33:\tk_eff = 0.929811\tres = 8.248E-03\n", + "[ NORMAL ] Iteration 34:\tk_eff = 0.936517\tres = 7.720E-03\n", + "[ NORMAL ] Iteration 35:\tk_eff = 0.942827\tres = 7.219E-03\n", + "[ NORMAL ] Iteration 36:\tk_eff = 0.948757\tres = 6.744E-03\n", + "[ NORMAL ] Iteration 37:\tk_eff = 0.954322\tres = 6.295E-03\n", + "[ NORMAL ] Iteration 38:\tk_eff = 0.959539\tres = 5.871E-03\n", + "[ NORMAL ] Iteration 39:\tk_eff = 0.964425\tres = 5.472E-03\n", + "[ NORMAL ] Iteration 40:\tk_eff = 0.968996\tres = 5.096E-03\n", + "[ NORMAL ] Iteration 41:\tk_eff = 0.973268\tres = 4.744E-03\n", + "[ NORMAL ] Iteration 42:\tk_eff = 0.977259\tres = 4.413E-03\n", + "[ NORMAL ] Iteration 43:\tk_eff = 0.980982\tres = 4.104E-03\n", + "[ NORMAL ] Iteration 44:\tk_eff = 0.984454\tres = 3.814E-03\n", + "[ NORMAL ] Iteration 45:\tk_eff = 0.987689\tres = 3.543E-03\n", + "[ NORMAL ] Iteration 46:\tk_eff = 0.990702\tres = 3.289E-03\n", + "[ NORMAL ] Iteration 47:\tk_eff = 0.993505\tres = 3.053E-03\n", + "[ NORMAL ] Iteration 48:\tk_eff = 0.996112\tres = 2.832E-03\n", + "[ NORMAL ] Iteration 49:\tk_eff = 0.998536\tres = 2.627E-03\n", + "[ NORMAL ] Iteration 50:\tk_eff = 1.000787\tres = 2.435E-03\n", + "[ NORMAL ] Iteration 51:\tk_eff = 1.002878\tres = 2.257E-03\n", + "[ NORMAL ] Iteration 52:\tk_eff = 1.004818\tres = 2.091E-03\n", + "[ NORMAL ] Iteration 53:\tk_eff = 1.006618\tres = 1.937E-03\n", + "[ NORMAL ] Iteration 54:\tk_eff = 1.008287\tres = 1.793E-03\n", + "[ NORMAL ] Iteration 55:\tk_eff = 1.009834\tres = 1.660E-03\n", + "[ NORMAL ] Iteration 56:\tk_eff = 1.011268\tres = 1.536E-03\n", + "[ NORMAL ] Iteration 57:\tk_eff = 1.012595\tres = 1.421E-03\n", + "[ NORMAL ] Iteration 58:\tk_eff = 1.013824\tres = 1.314E-03\n", + "[ NORMAL ] Iteration 59:\tk_eff = 1.014962\tres = 1.215E-03\n", + "[ NORMAL ] Iteration 60:\tk_eff = 1.016015\tres = 1.123E-03\n", + "[ NORMAL ] Iteration 61:\tk_eff = 1.016988\tres = 1.038E-03\n", + "[ NORMAL ] Iteration 62:\tk_eff = 1.017889\tres = 9.595E-04\n", + "[ NORMAL ] Iteration 63:\tk_eff = 1.018721\tres = 8.864E-04\n", + "[ NORMAL ] Iteration 64:\tk_eff = 1.019490\tres = 8.187E-04\n", + "[ NORMAL ] Iteration 65:\tk_eff = 1.020201\tres = 7.560E-04\n", + "[ NORMAL ] Iteration 66:\tk_eff = 1.020858\tres = 6.980E-04\n", + "[ NORMAL ] Iteration 67:\tk_eff = 1.021464\tres = 6.444E-04\n", + "[ NORMAL ] Iteration 68:\tk_eff = 1.022024\tres = 5.947E-04\n", + "[ NORMAL ] Iteration 69:\tk_eff = 1.022541\tres = 5.488E-04\n", + "[ NORMAL ] Iteration 70:\tk_eff = 1.023017\tres = 5.063E-04\n", + "[ NORMAL ] Iteration 71:\tk_eff = 1.023458\tres = 4.670E-04\n", + "[ NORMAL ] Iteration 72:\tk_eff = 1.023863\tres = 4.308E-04\n", + "[ NORMAL ] Iteration 73:\tk_eff = 1.024238\tres = 3.972E-04\n", + "[ NORMAL ] Iteration 74:\tk_eff = 1.024583\tres = 3.663E-04\n", + "[ NORMAL ] Iteration 75:\tk_eff = 1.024902\tres = 3.376E-04\n", + "[ NORMAL ] Iteration 76:\tk_eff = 1.025195\tres = 3.112E-04\n", + "[ NORMAL ] Iteration 77:\tk_eff = 1.025466\tres = 2.868E-04\n", + "[ NORMAL ] Iteration 78:\tk_eff = 1.025715\tres = 2.643E-04\n", + "[ NORMAL ] Iteration 79:\tk_eff = 1.025945\tres = 2.435E-04\n", + "[ NORMAL ] Iteration 80:\tk_eff = 1.026157\tres = 2.244E-04\n", + "[ NORMAL ] Iteration 81:\tk_eff = 1.026352\tres = 2.067E-04\n", + "[ NORMAL ] Iteration 82:\tk_eff = 1.026531\tres = 1.904E-04\n", + "[ NORMAL ] Iteration 83:\tk_eff = 1.026697\tres = 1.753E-04\n", + "[ NORMAL ] Iteration 84:\tk_eff = 1.026849\tres = 1.614E-04\n", + "[ NORMAL ] Iteration 85:\tk_eff = 1.026989\tres = 1.487E-04\n", + "[ NORMAL ] Iteration 86:\tk_eff = 1.027118\tres = 1.369E-04\n", + "[ NORMAL ] Iteration 87:\tk_eff = 1.027237\tres = 1.260E-04\n", + "[ NORMAL ] Iteration 88:\tk_eff = 1.027347\tres = 1.160E-04\n", + "[ NORMAL ] Iteration 89:\tk_eff = 1.027447\tres = 1.067E-04\n", + "[ NORMAL ] Iteration 90:\tk_eff = 1.027540\tres = 9.823E-05\n", + "[ NORMAL ] Iteration 91:\tk_eff = 1.027625\tres = 9.039E-05\n", + "[ NORMAL ] Iteration 92:\tk_eff = 1.027704\tres = 8.317E-05\n", + "[ NORMAL ] Iteration 93:\tk_eff = 1.027776\tres = 7.652E-05\n", + "[ NORMAL ] Iteration 94:\tk_eff = 1.027843\tres = 7.040E-05\n", + "[ NORMAL ] Iteration 95:\tk_eff = 1.027904\tres = 6.476E-05\n", + "[ NORMAL ] Iteration 96:\tk_eff = 1.027960\tres = 5.957E-05\n", + "[ NORMAL ] Iteration 97:\tk_eff = 1.028012\tres = 5.479E-05\n", + "[ NORMAL ] Iteration 98:\tk_eff = 1.028059\tres = 5.039E-05\n", + "[ NORMAL ] Iteration 99:\tk_eff = 1.028103\tres = 4.635E-05\n", + "[ NORMAL ] Iteration 100:\tk_eff = 1.028143\tres = 4.262E-05\n", + "[ NORMAL ] Iteration 101:\tk_eff = 1.028180\tres = 3.919E-05\n", + "[ NORMAL ] Iteration 102:\tk_eff = 1.028214\tres = 3.603E-05\n", + "[ NORMAL ] Iteration 103:\tk_eff = 1.028245\tres = 3.313E-05\n", + "[ NORMAL ] Iteration 104:\tk_eff = 1.028274\tres = 3.046E-05\n", + "[ NORMAL ] Iteration 105:\tk_eff = 1.028300\tres = 2.800E-05\n", + "[ NORMAL ] Iteration 106:\tk_eff = 1.028324\tres = 2.574E-05\n", + "[ NORMAL ] Iteration 107:\tk_eff = 1.028347\tres = 2.366E-05\n", + "[ NORMAL ] Iteration 108:\tk_eff = 1.028367\tres = 2.175E-05\n", + "[ NORMAL ] Iteration 109:\tk_eff = 1.028386\tres = 1.999E-05\n", + "[ NORMAL ] Iteration 110:\tk_eff = 1.028403\tres = 1.837E-05\n", + "[ NORMAL ] Iteration 111:\tk_eff = 1.028419\tres = 1.688E-05\n", + "[ NORMAL ] Iteration 112:\tk_eff = 1.028434\tres = 1.551E-05\n", + "[ NORMAL ] Iteration 113:\tk_eff = 1.028447\tres = 1.426E-05\n", + "[ NORMAL ] Iteration 114:\tk_eff = 1.028460\tres = 1.310E-05\n", + "[ NORMAL ] Iteration 115:\tk_eff = 1.028471\tres = 1.204E-05\n", + "[ NORMAL ] Iteration 116:\tk_eff = 1.028481\tres = 1.106E-05\n", + "[ NORMAL ] Iteration 117:\tk_eff = 1.028491\tres = 1.016E-05\n" ] } ], @@ -1449,8 +1449,8 @@ "output_type": "stream", "text": [ "openmc keff = 1.028263\n", - "openmoc keff = 1.028538\n", - "bias [pcm]: 27.5\n" + "openmoc keff = 1.028491\n", + "bias [pcm]: 22.8\n" ] } ], @@ -1558,7 +1558,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 43, @@ -1567,9 +1567,9 @@ }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW0AAADDCAYAAABJYEAIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcFNW1B/DfmRl2ZpAdBoGRVQXZlyjgjHHlqSBuwSWB\nGJen5sU8TRSjzwHME6LGuBtfHgFjUCIqQV8SxY1BFlkFBMEFGEA2QZBF1pk574+qgZ6hu071LN19\n8ff9fPjQ03W67u2q06e7q+vWFVUFERG5IS3ZHSAiovBYtImIHMKiTUTkEBZtIiKHsGgTETmERZuI\nyCEs2kkkIs+JyH2VePy9IvI/VdknouogIgNFZFUlHt9aRPaIiFRlv1yUUkVbREaKyHIR+U5ENovI\nsyLSIEFtF4rIQRFpVO7+j0WkRETaRNzXT0T+ISK7RGSHiHwkIiNjrHeEiBT5CbfX//9JAFDVW1X1\nvyvaZ1Udp6o3V/TxsZTr87f+Nrg4jsdPFJGxVd0vVziUx2eJyHv+ft4lItNF5LRyj8sUkcdFZL0f\n94WIPFZ+/RHxJRF5vldEdgKAqs5W1dOiPSYMVd2oqllaDQNLyvV5o4j8Puybg4jkisjGqu5TkJQp\n2iJyF4BxAO4CkAXgBwDaAnhHRDIS0AUFsA7ANRF96gqgjr+s9L4zAbwH4AMA7VW1CYBbAVwYsO65\nfsJl+v//ojqeQBUr7fNJAJ4DMEVEspLdqVTnWB6/DWAagJYATgGwHMAcEcnxY2oAeB/AaQAuUNUs\nAGcC2AGgX0D73SLyPWpxTzFH+wwgF8CPANwQ8rGCiO2aEKqa9H8AMgHsBXBFufvrAfgawEj/73wA\nUwFMAbAHwCJ4G7s0viWAV/3HrAHwHxHL8gH8DcAL/mM/AdArYvk6AL8BsCDivkcA3AugGEAb/74P\nATwZx3MbAWBWjGUTAYz1bzcG8CaAXQC+AVAQEXcPgK/8fq8CcE7Ec3oxIm4IgBUAdsJ7sZ1a7vnd\nBWCZ38bLAGqG6TO8F3wJgN4R970CYIu/rpkATvPvvwnAYQAH/f5OD7Fv+gJYCGC3v85Hk52T34M8\nngXgqSjP4Z8AJvm3b/T3R504tkEJgHZR7s8FsDFETkfNBXhvfCUA0iK20XT/tfI5gBvDbiOrz/5j\nn4r4eySAT/11fQngZv/+ugD2Ayjy9/seAC3gFfJRfux2fz+f5D+mFoAX4b3x7QIwH0DTuPIs2Ynu\nP5EL4b3Q06IsmwRgcsTOOARgGIB0eEVorX9b/OS/z/87x99o50c8dr/flgB4CMC8csn+Qz+BOsP7\nFrIBQGt/p7aBV7yKAOTG8dzCFu2HADzrt5sOYIB/fye/H839v9sAOCXiOf0lIm6f/xzSAfwawBcA\nMiKe30cAmgM4yU/Cm60+++u6HV4RblIukesCqAHgMQAfR3te/t/WvpkL4LqIF0K/ZOfk9zWP/f26\nyb/9MoCJcW6DoKK9IUROR80FeEW7GMeK9iwAT/n51x3eG1xemG0U1GcApwLYDOAXEcsHA8jxbw8C\n8B2AHuWfV0T8Hf7zaOn37zkAL/nLbob3ZlPL71tPAPXj2capcnikCYAdqloSZdkWf3mpxao6TVWL\n4RWLWvC+gvaFV1T+W1WLVbUQwP8CGB7x2Nmq+rZ6W+9FAN2itPcivKJ1PrzE3xyxrCG8F8GWOJ/f\nmSKy0z9uuFNEon21PAL/a6rf/zn+/cUAagLoKiIZqrpBVddFefzVAP5PVd/3t82j8F6cZ0XEPKGq\n21T1W3if6ntYfQZwAMDDAK5X1R2lC1V1kqruV9UjAMYC6C4imTHWZe2bIwA6iEhjf50LAvqVylzJ\n40aInceR/WwcI8ayJCLXH4+yPCinD8PIBRFpDe8wzT2qekRVl8HbRj+JCAuzjcr3eR+8DzMfwCu0\nAABV/Ze/H6CqHwKYAa94x3ILgPtUdUvE6+NKEUmDl+uNAXRSz8equs/oWxmpUrR3AGjiP6nyWvrL\nSx096O/vkE0AsuG9E7fyE2WniOyC95WwWcRjt0bc3g+gdpQ2/wrgWnifOP5SbtkueO/KLUM+r1Lz\nVLWRqjb0/49WlB6B91V4hoh8KSL3+M9xDYBfAhgNYJuIvCQiLaI8PhvA+tI//G2zEUCriJhtEbf3\nA6hv9Rnep/I3AJxdukBE0kRkvN/Pb+F9ulOULUqRrH1zA7xPhatFZH48P3qmmBMhjyP7+U2MGEvP\niFz/ZfmFMXK6tJ2fwc6FlgB2qur+iPvWo2yuh9lG5ftcH96Hn/7wDmkBAERksIjME5Fv/P0xGLFz\nHfD24bTSfQjvjeAIvG+5L8L7LWGKiHzlv47SA9Z1nFQp2vPgfV28PPJOEakPbwO9G3F364jlAuBk\neJ8iNgJY6ydKaYFsoKqXxtMRVd0ArwgNBvB6uWUH/L5eEc86Q7a7T1V/part4R2bvlNEzvGXTVHV\nQfCSAQB+F2UVmyOWl2oN77hhZfq1H8BtAH4sIt39u68FcCmAH6r3Q2UOvK96pb+4a7nVBO4bVV2j\nqteqalN4n+pfFZE6lel3kriSx/v9vl4V5aFXR/TzXQAXVmBfmGdeRMnp8f79YXJhM4BGIlIv4r42\n8N74Kkr89l+FdxgxHwBEpCa83xcehnfsuSGAfyF2rgPeoZ/B5fZhPf+Td5GqPqiqXeB9C74UZb8h\nmFKiaKvqHnhfIZ4SkQtFJMP/Bftv8DbAXyPCe4vIZf6703/CO9b6EYAFAPaKyN0iUltE0kWki4j0\nCWg6VnLdAK8gHYiy7G4AI0XkrtLTnkSku4i8HP4ZR+mIyMUi0t7/cy+8Y44lItJJRM7xk+cwvMMV\n0b5+vwLgYj82Q0R+BW/bzKtMvwBAVXfB+/qZ79+VCa847fJfOONQNnm3AWgX8XfgvhGR60Sk9JPL\nbn9d0Z5jSnMsj0cBGCEiPxeR+iLSUER+C+8QTenpmi/CexN5TUQ6i6exeOMDLgqxSaJ3NiCnjVwo\nLaxfwTtmPE5EaolIN3if0F8MajaOLo4HcJOININ3GKcm/MNeIjIYwAURsdsANJayZ1Y9D+Ah8U+v\nFJGmIjLEv50nIl39T/374H0CjyvXU6JoA4CqPgLvV+9H4e2sefC+8pznHxcqNR3eKTm7AFwHYJh/\n7K8EwCXwjtOug/fDxJ/gnXYVs9lot1V1naouibFsHrwfes4FsEZEdgD4I4B/xPWEj9cRwLsishfA\nHADPqGoBvGOd4+H9Cr0ZQFN4X5fLPhHVzwFcD+BpP/ZiAJeqalH551BBjwMYLN7pY3+BV4Q2wTtb\nZW652AkAuvhfD18PsW8uArBSRPYA+AOAH6nqoUr2NykcyuM58H6ouwLecet18H7QG+AfvoCqHgZw\nHoDVAN7xn89H8I7Jzg/Rl1iCcjooFyLXfQ280xQ3A3gNwH+p6gcBbQb1q8wyVV0BoADAr/3jzXcA\nmOof6hgOb9+Vxn4G7wfbtX6+twDwhB8zQ0R2w3t9lP6O1QLeJ/fdAFbCO34e9GZzHPEOp7lBRPLh\nnRsd19cJolTCPKbKSJlP2kREZGPRJiJyiFOHR4iIvu/4SZuIyCGVuoCNf9rP4/CK/wRVPe78YRHh\nR3mqVqpa5ZfrZG5TKoiW2xU+POKfZ/g5vFPfNsO7yMtwVV1dLk6xNOI0xOdGA7eOLrOuId2mmO29\nMXe4GXPTgCfNmDdxiRnzRLlBXFNHf4arRncuc9+PnnzTXA9y7G37s6FPBy5vpDvNdTzyRP7xd741\nGrhodMTfIfbz22PMkIztd5oxRS1DXAzwR1H6s3w00G30sb8nv2KvB8OrvGjHldvYEHHPYwAit88E\nu7EfjLZjRtn7rsuQRWbMl7vbH3df0fjfIWPUPUf/PrS+od1WN7utlZ8GnVbur+f0EOtZHmU95WpI\n7Rz7NdI+a63d1vS+ZgzGh3gdzS//OpoJIK/cfTcGriI3tyYKCppHze3KHB7pB+ALVV3vn386BcDQ\nSqyPKFUwtyllVaZot0LE9RPgDZduFSOWyCXMbUpZibgou/d1plTmSQlpsiqdntc42V2IX4e8ZPcg\nfs3zQgSthHf9nVTxWMRt9+aISBs4INldiF+fvGT3IE45IePmofSqE4WFsa8hVZmivQneRVpKnYxY\nF2wpdwzbNV3ygi7olaJO2KLdxf9X6rXq6En43IZ9jD+VpQ0cmOwuxK9vXrJ7EKeckHFn+v+AnJya\nWL/+0ahRlTk8shDedW/b+hd+GQ7vEp5ErmNuU8qq8CdtVS0WkZ/DuyB46WlRFZ5tmShVMLcplVXq\nmLaqvgXvguVEJxTmNqWqah/GLiKKYcGXix34+jvmesKcijt37LlmTEG+fS7md4ETunhOiTrjV1mn\nrik0Y4qzgt83uzZbaK7jLv29GbNVo012U9b9/7TXE+oCtEUhTpu+LUTe9foiRGOdq2VwTRjeedpR\nzpE/ys5H9LCPKbdYYp9j3FOWmjET9admzGtl52+Iam/MWeWOuUynmzF/F/ssykzda8ZcUXaOh6hG\nyiQz5uOSoNn3PNt6nWLGYNlsOwbvBS7NzW2LgoIbqvw8bSIiSjAWbSIih7BoExE5hEWbiMghLNpE\nRA5h0SYicgiLNhGRQ1i0iYgckpDBNU1L1gfGZEuMa/FE6ITPzZgStd+DXvlihBkjHYvNmMm40oy5\n6ttpZkzNk4y2ZtnPSTvYY0sk235OujL2lcVKZXf50oypqwfMmHUnn27GjNoUNHDFMz7tweQOrukf\nMHDsUIiVhLhKd8ngEDmwKUQODLNzYFmIQaDdFtk5IH3strDYfl5Le3cyY3pgtRmDv4fYhi3tbZj+\nVvBAQa+tEDW1dvDi3J5AwR/TOLiGiMh1LNpERA5h0SYicgiLNhGRQ1i0iYgcwqJNROQQFm0iIock\nZDb2x/GLwOXDO9vT701ffYEZcxn+ZcZs6GBPxNqmg32+8nV3hzgXM8Tk3IcHB7c162x7tuzzNs8x\nYxan28+pt31aNLZ0sS8Cr+PtzwJfbGptxtwmz9gdSrZRsRe1GGJPXrBpfke7jX72ucH76tjbvH43\nOweOLLDPn/+6jz0JQrMVdlvbejcwY4pQw4zRfnZbe1fY52Bn7bfPLS8O8Tk3+wH7PPZtb7QLDmgM\n4I/RF/GTNhGRQ1i0iYgcwqJNROQQFm0iIoewaBMROYRFm4jIISzaREQOYdEmInJIQgbX9Jf5gct7\nr55trmNJoT3IBOvtk+zbNA8xKOY2O+TUm5eYMasf6WXG7GhwUuDy8560B87IYDMEvZ+3n7fssdej\n79vbWELsqiky3IxJR4gLzidZl6GLYi5rgw3m42WLvV/21LG3+cEQEy7UX2u31anInmwk65UjZoxM\ntvvT4vrdZkzdq/bbKwrxvA4etFeDuvZ2znrJbqtn2lIzZuPQbwKX5yATBTGW8ZM2EZFDWLSJiBzC\nok1E5BAWbSIih7BoExE5hEWbiMghLNpERA5h0SYickilBteISCGA3QBKABxR1X7R4jpO3xS4nt8O\nuctuLMceaJH+Z/vE9wfG3GvG5J863oxZ/WWI97ve9mwZ2Qg+yV562+3oIbsd3GjPyqEP2m0djLqH\ny6pd395XW/EHMybEs6o2YXN77Z7YM/m8l3Wu3dAwe1tldQ0xK806e2vJdjsHsn4dIq8X223p+3Zb\nco7dVuaSIns9O+xt2LRJiOdlTCYDALjMbmuS2jPy5OxZF7g8Oz12aa7siMgSAHmququS6yFKNcxt\nSkmVPTwiVbAOolTE3KaUVNmkVADviMhCEbmpKjpElCKY25SSKnt4ZICqbhGRpvASfJWq2ld/Ikp9\nzG1KSZUq2qq6xf9/u4hMA9APwHGJrS+POfZH11zIGXmVaZa+x3bO/AQ7Z66o9nbC5vaRcQ8fvZ02\ncADSB4W4xCFRFMUfzkbJbO+qnmvSYh8EqXDRFpG6ANJUdZ+I1ANwAYAxUWOvya9oM0RlNMo7A43y\nzjj699oxf6vyNuLJ7Rr33l3l7dP3U/qggUgfNBAA0D49A2vH/S5qXGU+aTcHME1E1F/PZFWdUYn1\nEaUK5jalrAoXbVVdB6BHFfaFKCUwtymViWqImVwq04CIYlXwCenFjewZIwY3fd2M2YMsM+Ygapsx\ni58ZaMbcfPsTZsyfl99uxtRuuzNw+XfFTc11yCQzBCVz7QERC149w4zpv2y5GfPv3e2BM22x3oy5\nb0P0r4dl5NSCqiZlHI6IKJbGzu1nu48013GmzjNjjqCmGdO56DMzpv5v7AEvcx+236sy1F5Pvx6f\nmDELltn5VqR2bTjrHnummL3j7M+nn6d3MmMyYM/aMw9nmjG3L58YuDy3HlDQMS1qbvM8VCIih7Bo\nExE5hEWbiMghLNpERA5h0SYicgiLNhGRQ1i0iYgcwqJNROSQyl7lL5zJwYufH/tjcxX/3Hm5GZNR\naJ/0n9vrLTNGDtkDjhZJHzNmTrdeZkz/WcGDVf5x9jnmOi7u+YEZ85s7HzBj7j/8oBkzq3tfM+Y/\n8JQZc/r9wTN3AMB9pz9mxiRb124LYy7bq5nm47sv/sKM2drbngklc6o9wwtid/WoDNivof4j7AFW\nY+2xNXggxHrmv9DNjElbaL9es6bag2JOHv6VGdNi8W4zZkafC8yYLt0WBS7PQSYKYizjJ20iIoew\naBMROYRFm4jIISzaREQOYdEmInIIizYRkUNYtImIHMKiTUTkkMQMrhkRfPJ7zRCzQaQ1sk/6Lx5n\nz3KxpNdpZgzuDJ5pBwD6qz1zTf8r7cED8mrw87q4d4j31SvsiVvGnWNPrqzfRZ27towPatqDffIx\n3oy5+EF7JiKMSMqENHFZ+WnsQVZDT7/SXkFvO9darAiRAy+F2FYfhBg4081ua8xKu638ErutsQEz\njpf6r6X2KB1dZm9DucRuq3nXPWZMmP11mbY1Y0Z9+mTg8iZ1Yy/jJ20iIoewaBMROYRFm4jIISza\nREQOYdEmInIIizYRkUNYtImIHMKiTUTkkIQMrrmv/f2Byw9JTXMdt6s9i8k1j/Q0Y/5HbjFj7tV2\nZkxjud6M+WZqwBnypeuZFDwg6L3FA+x14BszpoFmmzGnmBHAKWLPOLOqxN5+/5q1NkRrqa/L6bFn\nIJmOIebjf7XIHhD2dZ8sM6bFtfbAkJIf2m0tWGbPFJP/E3vQ2Jh0u618+yWE+S+cYcb0D/G89Ca7\nra+72jMNNQuxv/7e93YzJihvAM5cQ0R0wmDRJiJyCIs2EZFDWLSJiBzCok1E5BAWbSIih7BoExE5\nhEWbiMgh5uAaEZkA4BIA21S1m39fQwB/A9AWQCGAq1V1d8x1IHjmmlsLXjA7+sfcn5gx9bDfjHlc\n7zBjGj550IyZdMcIM2aYTDNjGg9ZGbj8NKwy19Hq/Z1mDILHNwEApsy1B4P8eOqrZkzfK2MNC4jQ\n0A7J+IM9YKTor/Z6YqmK3F6xvG/M9Wd2f8bswyd92psxh1HLjKl79WdmTOaSIjOmKM0ePDL/LyEG\n4CyzB+CEWU8R7P6gd3B9AYA9V9cwYzZKazNma5/DZkym7jVjVi6PPeMRADSpF3tZmE/aEwFcWO6+\nUQDeVdXOAN4HcG+I9RClGuY2Occs2qo6G8CucncPBVD68fgFAJdVcb+Iqh1zm1xU0WPazVR1GwCo\n6lYAzaquS0RJxdymlFZVP0TaB5WI3MTcppRS0av8bROR5qq6TURaAPg6KLhg9IdHb7fNa4OcPHuK\neaJoSuZ8CJ0zuzqbiCu38dzoY7f75AF986qxa3RCWzgTWDQTAFAYcOHTsEVb/H+l3gAwEsDvAIwA\nMD3owbmjB4VshihY2oBBwIBj+VT0yPjKrrJSuY1bR1e2fSJP37yjb/o59YD1T42NGmYeHhGRlwDM\nBdBJRDaIyE8BjAdwvoh8BuBc/28ipzC3yUXmJ21VvTbGovOquC9ECcXcJheJavX+ziIieoc+FBiz\nsCT2AIVSs+Vcu7HP7N9VNU3MGOlYbLe1OkRb/wjR1l3Bbb0L+9DSeVPnmjG4yn5Om9DEjGn1fPkz\n5KK4xW5rBnLNmGuKXzZjdtU4Gapqb+hqICJae9eOmMsLG9hzATVDzHE7x/Szc61knb0J0raHyOtf\nh8jrxSHy+v0QbZ0Toq2+Idp62G5Lm4Y456JdiLbm2219jQZmTNvdhYHLB6Vn4N2sBlFzm8PYiYgc\nwqJNROQQFm0iIoewaBMROYRFm4jIISzaREQOYdEmInIIizYRkUMqesGouEw6NDJw+e9r3WWuY2XJ\nLWZM16ftvkhHezBRyUv2bBkT8q8zYw6das84cn1x7cDlOzIuMtex66qAq8v4Gj5hP6dW7ext88Yt\n55sxlz5kt/X0fVPMmEHps+z+mBHVq32DtTGX3YA/m49/83V7W+37xO7HgUP2vmvS1G7ru612Sch6\nxZ4BRy8NMePMzXbIvqvs9dRrYsfsCDG5U5199jasP81u66eXTzVjOjRYE7i8FTJjLuMnbSIih7Bo\nExE5hEWbiMghLNpERA5h0SYicgiLNhGRQ1i0iYgcwqJNROSQhMxcg44lgTE15u0x13NDo4lmzHO4\nw+7P6BDvUyEmE9EQE3NMe8oeGHOargpcfrc8bK5jDPLNmMNqD8DZo1lmzPliD3iZg95mzIAFH5sx\nm/o3MmNay86kzlyDv8fO7ZZDggdQAMCm+Z3shvrbyba3jp3XmWfYTS1a0MWMaY2NZkzzFfZremtX\ne4aXr3CyGdOn70ozZs8KO0WyDoR4US+wt3N2vy/NmK1vtAtcntsYKBiUxplriIhcx6JNROQQFm0i\nIoewaBMROYRFm4jIISzaREQOYdEmInIIizYRkUMSMnONFOwPXH7kr/agjoO/sAeH6EZ7VomX84ea\nMdfINDNmTJr9ftfwmbfNmGHFwSf0v7HJbkf/aQ8ckBvtgQPvpJ1txrxYcrUZM3LVIjPm8n6TzZgw\ngziA+0LEVKPxsbf9ljHtzYenD7VngSmGndf1J4cYX3S5nQM1YQ/2abZor91Wn+ABdQDQYrGd29t6\nN7PbWmi3lTUtxOvoI3s7p79tt4V/t0NQ29hfPWMv4idtIiKHsGgTETmERZuIyCEs2kREDmHRJiJy\nCIs2EZFDWLSJiBzCok1E5BBzcI2ITABwCYBtqtrNvy8fwE0AvvbDfqOqb8VaR2bD4JPx95xa1+zo\nZrQyYzKm2gMVpt1pzyajve2T7NeUPG/GNNHtZszh3cFt3Zhtt3PgxjpmTC5uMmOWqT1w5tXDV5ox\n2tP+LPDaL683Y9o9bM9IUpnBNVWR2/hodEAL55l9UAwwY1o98LkZ0xNLzZg/w54pZq4MM2Pe7nOh\nGTMUbc2Y6b1vN2MyxR7I01Lt5/XTYa+aMR9rDzMGt9khWDo7RNB7wYtrxd5+YT5pTwQQbS89pqq9\n/H+xk5oodTG3yTlm0VbV2QB2RVmUlHn5iKoKc5tcVJlj2j8XkaUi8r8iYn8/IXIHc5tSVkUvGPUs\ngLGqqiLyWwCPAfhZrOCDv330WINnn4WMs8+qYLP0fXdg5kIcmLmwOpuIK7eBmRG3c/x/RBVR6P8D\nCgtjf1aoUNFWLfML258AvBkUX/v+X1WkGaLj1Mnrizp5fY/+/e2Y56p0/fHmNpBXpe3T91kOSt/0\nc3LaYv36N6JGhT08Iog4ziciLSKWXQ5gRQV6SJQKmNvklDCn/L0E7+NEYxHZACAfwDki0gNACbzP\n87dUYx+JqgVzm1xkFm1VvTbK3ROroS9ECcXcJheJqlZvAyKKT4zZHs4KcYbVO3Y/L+przzjz9o2X\nmTHjJtxhxty74VEz5vW2Q8yYdA0eEHTt/pfNdZxZb54Z84l2NWO2Nm5nxtTbsMOM+e7SJmZMrWnR\nzrQr69B39qAhnFwXqpqUU/RERIENARET7JX8YLQdM8rO/S5D7NmC1uy29+/B9Y3MmK7d7B+CV37a\nx4zpcrrd5xXL+5oxdXK+MWPaZa0zY1ZOt/uM8XYI5o8JERTw2zaA3NxaKChoHjW3OYydiMghLNpE\nRA5h0SYicgiLNhGRQxJftBfOTHiTlbVm5lfJ7kLcds78JNldiFvJ7DBXR0tl9g/CqabYxW3uXA0p\nrNK1sWiHsHbmpmR3IW67Zro3JqRk9pxkd6GS3CvaTm7zRTOT3YM4FVbp2nh4hIjIIRW9YFRcetU+\ndntzBpBdu1xAiGuPw54nAR1wkhmz3b42O5qjdZm/62P1cff1qmmfGtwAHcyYdBQHLu+RZu+iDlEu\nbr8LtcvdX9NcT3Z3MwR1QvTnQEd7PTXTj5/8YaMIWkfcf7iGvY2X2E1Vq169ahy9vXlzOrKza0Qs\nbWmvoHOIRkJcZ7B9iBdIVtRtnlZmmx8KcWp8mLZqlX+NR9EuxHpqRunP5hpAdsT9tdLsSUtODtPn\nMNdzDLO/jpTd75s310d2dvlcqIEgnTploKAg+rLEDK4hqkbJHVxDVH2i5Xa1F20iIqo6PKZNROQQ\nFm0iIocktGiLyEUislpEPheRexLZdkWJSKGILBORj0VkQbL7E42ITBCRbSKyPOK+hiIyQ0Q+E5G3\nU2narBj9zReRr0Rkif/vomT2MR7M6+rhWl4DicnthBVtEUkD8DS82a+7ALhGRE5NVPuVUAIgT1V7\nqmq/ZHcmhmizio8C8K6qdgbwPoB7E96r2E6YWdCZ19XKtbwGEpDbifyk3Q/AF6q6XlWPAJgCYGgC\n268oQYofRooxq/hQAC/4t18AYF+TNkFOsFnQmdfVxLW8BhKT24ncaa0AbIz4+yv/vlSnAN4RkYUi\nclOyOxOHZqq6DQBUdSuAZknuTxguzoLOvE4sF/MaqMLcTul32hQxQFV7Afg3ALeLyMBkd6iCUv3c\nzmcBtFPVHgC2wpsFnaoP8zpxqjS3E1m0NwFoE/H3yf59KU1Vt/j/bwcwDd7XYRdsE5HmwNHJar9O\ncn8Cqep2PTZo4E8A7ClLUgPzOrGcymug6nM7kUV7IYAOItJWRGoCGA4g+hzxKUJE6opIff92PQAX\nIHVn5y4zqzi8bTvSvz0CwPREd8hwosyCzryuXq7lNVDNuZ2Qa48AgKoWi8jPAcyA92YxQVVXJar9\nCmoOYJoneEpvAAAAZElEQVQ/XDkDwGRVnZHkPh0nxqzi4wFMFZEbAKwHcHXyeljWiTQLOvO6+riW\n10BicpvD2ImIHMIfIomIHMKiTUTkEBZtIiKHsGgTETmERZuIyCEs2kREDmHRJiJyCIs2EZFD/h/n\n2ajR7Q6wgAAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAADHCAYAAADmiMMCAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHxdJREFUeJzt3Xu8VHW9//HXZ29uKmwUUARviKAGBAgImhqaVlomamrS\nsbQ07WRpDz0VpUc21S/JrLxk2akIy9vxRuDJMi/9IFMRMVTwkoqoIBcR5abcP+ePtbZn2O49n7Wv\nM7N6Px+P/dh7z/rM+n7XWp/5zJo1852vuTsiIpIPVaXugIiItB4VdRGRHFFRFxHJERV1EZEcUVEX\nEckRFXURkRxRUS8jZna9mf1nC+7/HTP7dWv2SaQtmNkRZvZ8C+6/t5mtM7Pq1uxXHpR1UTezs8zs\naTN7x8yWmdkvzGzndmp7kZltMrNe9W7/h5m5mfUruG20md1jZm+b2Soze8zMvtDIes8ys61pQtb9\n/AzA3b/s7t9rbp/d/Qfufk5z79+Yen1ekx6Tk5tw/6lm9v3W7lelqKA8/pCZPWhma81stZndbWaD\n6t2vxsyuMrNX03x4Kf1/u/UXxLuZrS/I9bcB3P1v7n5Ac7fL3V91967uvrW562hMvT6/bmbXmlnH\njPc90swWt3afmqJsi7qZXQz8EPgG0B04BNgHuM/MOrVTN14Gxhf06YPAjvX6eSjwIDATGAD0BP4d\nOLbIeh9JE7Lu56ut3vPW94i7dwV2Bn4G3GxmPUvcp7JXYXn8F2A60BfYF3gS+LuZ9U9jOgEPAINJ\n8rsGOBRYCYwu0v6wglxvlyezVjAszfcPAycD55a4P9m5e9n9kCTLOuC0erd3Bd4Avpj+XwvcAfw3\nsBZ4guRg1MX3Be5M7/MycEHBslrgNuB36X0XAKMKli8CLgXmFNx2JXAJ4EC/9LaHgOuasG1nAQ81\nsmwq8P30717A/wBvA6uAvwFV6bJvAUvSfj8PHF2wTTcWrO+EdLveBv4/8IF62/cfwFPA6nQfdsnS\nZ5KC4MDogttuB5al65oFDE5vPxfYDGxKj+ndGY7NaOBxYA2wHPhJqXPyXyCP/wb8vIFt+BPwu/Tv\nc9Lj0bUJ+8CBAQ3cfiSwuOD/xnK6wVwA+qXr7lCwj2aQPFZeBL6UdR9FfU7v+/OC/78APJuuayFw\nXnr7TsC7wLb0uK9L+1UFTABeAt5M19cjvU8X4Mb09reBOUDvFuVdqRO/kZ16LLCl7oDVW3YDcEvB\nwdoMnAJ0JClSL6d/VwFzgcuATkD/9AB8vOC+G4BPANXA5cCj9R4Mx6QJ9oE0ZjHJWZanSbUjsBU4\nqgnbdhbZivrlwPXptnQEjgAMOAB4DehbkNz7FWzTjenf+wPrgY+m9/9mmuydCrbvsTTpeqRJ+uWo\nz+l+OD9NwO4FMV8EugGdgauAeQ1tV/p/dGweAT6X/t0VOKTUOfmvmsckBWxp+vetwA1N3AdhUQ9y\nusFc4P1FfRbwc5IiOZzkCfAjWfZRsT4DBwJLgbMKln8S2I/k8TgWeAcYUX+7CuIvBB4F9iR5fPyy\n4NifB9ydHoNqYCRQ05K8K9fLL72Ale6+pYFlS9Pldea6+x3uvhn4CclBPQQ4GNjV3b/r7pvcfSHw\nK+D0gvs+5O73eHJd7vfAsAba+z3weZLi+CzJ2USdXUgedEubuH2HpNff634OaSBmM9AH2MfdN3ty\nDdJJHnydgUFm1tHdF7n7Sw3c/zPAH939vnTfXAnsAHyoIOYad3/d3VeRJNbwqM8kD44rgU+5++q6\nhe4+xd3XuvtGkgfRMDPr3si6omOzGRhgZr3cfZ27P1qkX+WsUvK4B43ncWE/ezYSE3miINevaWB5\nsZwOc8HM9gIOA77l7hvcfR7wa5LtrZNlH9Xv83qSfXWnu0+tW+Duf3T3lzwxk+Sy1RFF1vVl4BJ3\nX1zw+DjFzDqk29eT5Elkq7vPdfc1Qd+KKteivhLolW50fX3S5XVeq/vD3beRnIX0JTkT6VtYPIHv\nAL0L7rus4O93gC4NtPl74LMkZ6u/q7fsLZKXWn0ybledR91954KfhorWj0jOrP9iZgvNbEK6jS8C\nXydJjBVmdquZ9W3g/n2BV+r+SffNa8AeBTH1t79r1GeSJ7IZJGf+AJhZtZlNTt80W0NydgjbF61C\n0bE5m+SVxnNmNsfMji/Sr3KWhzwu7OebjcRERhTk+gX1FwY5nSUX+gKr3H1twW2vUDzXG9pH2/WZ\n5PHwGeBz9d5QPs7MHk0/FPE2ySuAxnIdkmM4reD4PUvyRNab5LjcC9yavil7RdY3ZRtTrkX9EWAj\nyRsU7zGzrsBxJG/W1NmrYHkVyUuc10keJC/XK57d3P0TTemIu79C8lL4E8Bd9Za9k/b1001ZZ8Z2\n17r7xe7en+Ta+EVmdnS67GZ3P5z/ewn9wwZW8Xq6HAAzM5J9taSB2Kb0ax3JG8FjzezI9ObPAuNI\nXuZ3J3lpDMnLU9I+Fip6bNz9BXcfD+yWbtsdZrZTS/pdIpWSx+vTvp7awF1PK+jn/cDH2+JYNJbT\nGXPhdaCHmXUruG1vWp7r7u63kby3VQtgZp1J3t+4kuTa987APTSe65Acw+PqHcMu7r4kfRU+yd0H\nkbyKPp7tX2E0WVkW9fRl/STgWjM71sw6ps+Ut5Gcwfy+IHykmZ2cPut+neRB9CjJ9eK1ZvYtM9sh\nPZscYmYHN6NLZ5Ncn1vfwLJvAmeZ2TfqPg1iZsPM7NZmtPMeMzvezAakxXg1yTP7NjM7wMw+kibX\nBv7vjZn6bgM+aWZHp8/8F5Psm4db0i+A9HLNf5G8+QPJtfSNJGdyOwI/qHeX5STXgusUPTZmdoaZ\n7Zqesb6d3qehbSxrFZbHE4AzzewCM+tmZrtY8jHUQ9NtIO3va8CdZnagmVWZWU9Lxkc06UmmULGc\nzpIL7v4aSV5fbmZdzGxouq03NrdP9UwGxqeXeTqRXCp6A9hiZscBHyuIXQ70rHfp8Xrg/5nZPuk2\n7Wpm49K/jzKzD1ryefs1JJdjWpTrZVnUAdz9CpKXmVeSbOxskoQ6Or0uVWc6yUukt4DPASenz35b\nSZ71hpOcoawkuc7W2HXeYn15yd0fb2TZw8BH0p+FZlZX8O5pajv1DCQ5M1pHchb1c3f/K0lCTSbZ\nnmUkZzDfbqBfzwNnANemsZ8iuQ6+qYX9qnMVcJSZDSd5Of8KyZnRMyTFqNBvSK6Xvm1mf8hwbI4F\nFpjZOuBq4HR3f7eV+t2uKiiPHwI+TvKqYinJ8TwIONzdX0hjNpK8GnsOuC/dnsdILj3Mbmp/ChTL\n6ay5MJ7kFeLrwDRgorvf34I+vcfdnyb52PLF6SWeC0iemN8ieZU6oyD2OeAWklrwdnoZ6eo05i9m\ntpbk8TEmvcvuJJ98WkNyWWYm2z/ZN5m5V+4kGWZWS/IGwxml7otIcymPpTWV/Ew9TeiKoj63j0rs\nc6FK7L/63Pbaur8lP1M3M3d3iyMbvG8tJTjDaUmfS0V9bn9Z+19OZ+qVuM8rrc9t3d+KLuqloj63\nj0rsc6FK7L/63Pbaur8lv/wiIiKtp0Vn6mZ2LMk7u9XAr919cob7VO47s1IRWuMsSLkt5SjT5bzm\nFvX0c5X/JBl2vJjki2jGu/szwf2cecU/hnnC0Pgj3jMePj2M+dJhDY1I3t7dxIMVr/avhzGfuebu\nMIZ+8b4+e9zPii7v4avCdfzo6olxX/6c4bjfOykM6fDGRWHMlj41cVufydCfm26LYzi9xUW9RbnN\nq0UifhM3fkhtHDMh3leDT2jwk4vbeXH1fmHMxld2idsaGre14JlR8XoGZVjPU/F6uvSLHyP71SyM\n25qeYSjA5Ax5Ozt+HCXflda4sWM7MXNm70y53ZLLL6OBF919YfrZ51tJRhWKVDrltlSslhT1PSj4\nvgqSM5o9GokVqSTKbalYbf5GqZnVWjKTiOuao7SHwnxry88EK7elvWXJ7WLfUhZZQsGXEJF8AdH7\nvkDH3WtJvwynrlMtaFMk1ApvlCq3pSy19TX1OcBAM9vXkmmuTqfgOxBEKphyWypWs8/U3X2LmX2V\n5LuAq4Ep7r6g1XomUiLKbalkLbn8grvfQ8u/jVCk7Ci3pVK1+9cEmJlzUvHPqR9+133herJcNn34\nu0eHMTMnxp9FXV90QqDEvv5yGHPgS4vCmK01xZ9nh+w2J1zHxf7jMGaZ7x7GXHpPvB7+GIewJcMl\n7q9kyMMRL2Ro7IBWGXzUHMk19WJjBOJ8ZPjhYcjuT8SfsT7I5oUxv/UvhDF3bj+/R4PWbjc3RcNO\n9OlhzB8s/tRot+0mN2rYp7efA6RBZ9nUMOYf24rN7phYPmLfMIYnH4pjtpsv5f3Gjt2HmTO/2ObX\n1EVEpMyoqIuI5IiKuohIjqioi4jkiIq6iEiOqKiLiOSIirqISI6oqIuI5EhJBh/tuu2VojF97X3f\nnfQ++/PPMGabx89Zt71wZhhjA7eGMTdxShhz6tvTwphOOwdtzYq3yQfEY2+sb7xNvqA6jOk7+MUw\nZkd/N4x5ec9BYcyEJfHkH5OrvlfawUdjigys25hhJRm+tX3bcRlyYEmGHDgpzoEnOSCMGfp4nAM2\nKm6LufF2zRu5fxgznOfitv6QYR/2ifdh9Z+LD6RM2spQY7sUXzz2IJh5fZUGH4mI/KtRURcRyREV\ndRGRHFFRFxHJERV1EZEcUVEXEckRFXURkRxp0cxHzXUVFxRdfvoB8XSQ05/7WBhzIn8KY14dcFEY\ns/eA+PPa//bNDJ9FrYlDNh1XvK1ZHz4sXMcxr/89jJlbHW/TyPhj4SwdHE8S4JPjc4cXluwVxnzF\nros7VGoTGl+0+wnx5BZLZg+M2xgdfzZ63Q7xPu86NM6BzY/F4wdWjIonydhtftzW8pHdw5gtdAxj\nfHTc1tr58WfQa96JP1u/NcN5cd/L4s/xL5/Rv3hAT+D6cDWAztRFRHJFRV1EJEdU1EVEckRFXUQk\nR1TURURyREVdRCRHVNRFRHJERV1EJEdKMvhojM0uunzkcw+F63hiUTwIh1fiQQh7984waOgrcciB\n5z4Rxjz3oxFhzMruOxddfsw18cAiOy4MYeQv4+22NfF6/MF4H1uGQ3WrnR7GVJNhQoISGzzu8UaX\n7c2r4f1taXxc1uwQ7/MNGSbk6Lowbmv/LfFkNDW3bQ5j7Ka4P7ufsTqM2fHUd+IVZdiuDRvi1bBj\nvJ9rbo7bOqhqXhjz2rg3iy7vRzdmhmtJ6ExdRCRHVNRFRHJERV1EJEdU1EVEckRFXUQkR1TURURy\nREVdRCRHVNRFRHKkRYOPzGwRsBbYCmxx91FZ7jdw+pKiy79/wsXxSvrFA1Gqp8QDAy6b9O0wZuKB\nk8OY517M8Pw4Mp5tpS/FByHYyLgd3xi3wznxrC7+vbitDaPjprp0jY/VMn4axmTYqlbT3NxeuKbx\nmaAeqDk6XsFJ8b6qGZJhVqOX471lb8Q5UPONDHk9N27LH4zbsqPitro9sSVez8p4H+7aK8N2BZMR\nAXBi3NZUj2d06rfm5aLL+1ZnL9WtMaL0KHdf2QrrESk3ym2pOLr8IiKSIy0t6g7cb2Zzzezc1uiQ\nSJlQbktFaunll8PdfYmZ7QbcZ2bPufus1uiYSIkpt6UitehM3d2XpL9XANOA971tZma1ZuZ1Py1p\nTySLwnwzs9rmrEO5LeUoS243u6ib2U5m1q3ub+BjwPz6ce5e6+5W99Pc9kSyKsw3d69t6v2V21Ku\nsuR2Sy6/9AammVndem529z+3YH0i5UK5LRWr2UXd3RcCw1qxLyJlQbktlczc2/dSoJk5zxb/wP7W\nHvGMI8ftelcYs4aaMGYDXcKYudcdHsace/7VYcyUp84PY7rss6ro8vVbdw3XYVPDELY9HF8teOyO\nD4YxY558Koz58rB4YNE+vBLGXPLqD8MY+nWmVJdCzMyZ13hu/3zYWeE6DvVHwpjNdApjDtjyfBjT\n9TvxgKCHrxgexnTweD2jhz8dxjz2ZJxvWzyuDR/6VjzT0NrL4/PZf1bvH8Z0IJ716REODWPOf+q3\nRZeP3QlmDqzKlNv6nLqISI6oqIuI5IiKuohIjqioi4jkiIq6iEiOqKiLiOSIirqISI6oqIuI5Ehr\nTJLRdDcVX/zL734uXMU9q04OYzosigdFjB0Rj/62jfEArcctnhjn70NHhDFjZhUfzPPHDx8VruOT\nB/01jPnORZeFMZdu+l4YM2vYwWHM17g2jBl0afGZXwAuGfSTMKbUhgyd0+iytd4tvP+wuS+EMctG\nxjPpdLs9niGIxrv6ng7Ej6ExZ8YD0L4bjz3isgzrmX3D0DCmak78eK25PR40tOfpi8OY3eeuDmP+\nMupjYczgoY8XXd6PbswM15LQmbqISI6oqIuI5IiKuohIjqioi4jkiIq6iEiOqKiLiOSIirqISI6o\nqIuI5EhpBh+dWXxwQKcMs4lU9YgHRWy9PJ4l5YkRHwhjuKj4TE0AYzye+WjMKfHgCruj+HZ9cmSG\n5+FPxxP/XH7UxDDG108KY/7aKR4MNZHJYcwnvxfPZMWZ5T+384JnGh+ENm7QKfEKRsa5tvv8DDlw\nc4Z99dcMA4uGxm1NWhC3NXFb3NZ3q+K2/nNePIrJn4z3oR0ft9V7yJowJsvxOtH3CWMmPHNN0eW9\ndoy7Ukdn6iIiOaKiLiKSIyrqIiI5oqIuIpIjKuoiIjmioi4ikiMq6iIiOaKiLiKSIyUZfHTJfpcW\nXb7ROoXrON/jWXDG/+igMOa/7Lww5tveP4zpaWeEMW/eHo8g6Dm1+ICpB+YeFq+DN8OY7t43jNk3\njIB9LZ6x6Nlt8f7706yFGVorf4MHNT6DzXROCO//H4/HA+ZWjKoJY3b/bDxwZttH4rYeezKeaWji\n5+NBdZOq47Ymxg8hZt/wwTBmTIbt8i/Fba0YEs9UtVuG4/WHg88PY4rlDWjmIxGRf1kq6iIiOaKi\nLiKSIyrqIiI5oqIuIpIjKuoiIjmioi4ikiMq6iIiORIOPjKzKcDxwAp3H5Le1gP4b6AfsAg4zd3f\nytqoUXzmo3+feUO4juvHfj6M2Yl3wpir/MIwZpdrNoQxUy88M4w5yaaFMT1PWFB0+Qd4NlzHHg+u\nCmMoPv4LgFsfjgfLfO72O8KYg0/JMGxilzikw0/jATVbbozXU6ctcnv+Uwc3uqzbsOvC+z89ar8w\nZhOdw5gdT3s+jOn2xJYwZktVPLhm9u8yDFB6Mh6glGU9W4j7w8ji9QVgzWkdw5jXbK8wZtmoTWFM\nN18bxix4qvEZswB67RSu4j1ZztSnAsfWu20C8IC7DwQeSP8XqTRTUW5LzoRF3d1nAfVP/cYBdafT\nNwAntnK/RNqcclvyqLnX1Hu7+9L072VA71bqj0ipKbelorX4jVJ3dwgukotUIOW2VKLmFvXlZtYH\nIP29orFAM6s1M6/7aWZ7IpkV5puZ1Tbx7sptKVtZcru5RX0GUPdxjzOB6Y0Funutu1vdTzPbE8ms\nMN/cvbaJd1duS9nKktthUTezW4BHgAPMbLGZnQ1MBj5qZi8Ax6T/i1QU5bbkUfg5dXcf38iio1u5\nLyLtSrkteWTJe0Ht2KCZX+g/KBozZ1vjAzjqPGQZHnfPx1eXvCp+1WwDt8ZtPZehrT9maOvi4m3d\nzxHhOo65/eEwhlPjbVpCrzBmj19mGJdzXtzWXxgbxozfeksY81bHPSnVpRAz8y5vrWx0+aLu8VxS\nu7E6bmh0nGvbXo53QdUbGfL6Gxnyem6GvH4wQ1tHZWjr4AxtXRG35btmuPLcP0Nbs+O2VtA9jNln\n9aKiy4+o7sD9Nd0z5ba+JkBEJEdU1EVEckRFXUQkR1TURURyREVdRCRHVNRFRHJERV1EJEdU1EVE\nciQcUdoWpm48q+jyH3e+OFzHgm3nhTFDfhb3xQbGg6+23RzPtvKbif8Wxmw8MJ6x5oytXYouX9mh\n/pwO7/fWqZ3CmF2ujrdpj/7xvplx3kfDmE/9IG7rZ5fcGsYcUT0r7k8Y0bb2676w0WVfZEp4/7vv\nivfVuqfjfry7MT52vXaN21q/LC4RNbfFMyj5pzLMWHRuHLLu1Hg9O/WKY1ZmmBxsh3XxPuw6LW7r\nCyffHsYM6P5S0eV70C1cRx2dqYuI5IiKuohIjqioi4jkiIq6iEiOqKiLiOSIirqISI6oqIuI5IiK\nuohIjpRk5iMGbisa0/GRNeF6vtjjt2HML7gw7k9thue1DJPReIaJXaZdGw8c+oA/W3T5N+2KcB2T\nmBjGbPJ4gNIarwljPmrxgKC/MzKMOeyxf4QxS8b0CGP2slUlnfmIPzSe231OKD7ABGDJ7P3jhsbE\nybZ2hzivu30wburxxwaHMXvxWhjTe378mF42JJ4haDF7hjGjDl4QxqyZH6dIzbsZHtSPxfu57+gX\nw5hlM/oXXT62J8w8okozH4mI/KtRURcRyREVdRGRHFFRFxHJERV1EZEcUVEXEckRFXURkRxRURcR\nyZGSzHxkM98punzzjfGglw0XxINn/LV4VpJbJo4LY8bbtDBmUlX8/LjLdfeGMSdtLT7gYcaSuB2/\nJx5YYefEAyvuq/pwGPP7baeFMWc9+3gYc/Lom8KYLINc4JIMMW1ocuP7fumk/cK7V4+LZxHaSpzX\nXW/KMP7q5DgHOhEPhtrt8bVxW6OKDzgE2H1unNvLR+4WtzUnbqtmWobH0aPxfq6+N26LL8chdAmO\n10EZ1pHSmbqISI6oqIuI5IiKuohIjqioi4jkiIq6iEiOqKiLiOSIirqISI6oqIuI5Eg4+MjMpgDH\nAyvcfUh6Wy3wJeCNNOw77n5P1ka77VJ8sMKaA3cM1/E6e4QxHW6PB3JMuyiejchHxoMQXtr2yzCm\nl78RxmxaXbytc/rG7bx7zg5hzFi+FMY86fHAojs2nRLG+EHxucOdXz8jjOl/RTyjTVMGH7VFbvNo\nbZGFx4R3dw4LY/a47J9hzEHMC2OmEM809LCdFMbcO+rjYcw49gljpo88P4zpZvFApz4eb9cXTroj\njPmHDw9j+EocwryHMgQ9UHxx53j/1clypj4VaKjy/dTdh6c/2ZNepHxMRbktORMWdXefBaxqh76I\ntCvltuRRS66pf83MnjKzKWa2S6v1SKT0lNtSsZpb1H8B9AeGA0uBHzcWaGa1ZuZ1P81sTySzwnxL\nr5E3hXJbylaW3G7WtzS6+/KCRn4F/E+R2FqgtiBeyS9tyt0zfEVho/dVbkvZypLbzTpTN7M+Bf+e\nBMxvznpEyo1yWypdlo803gIcCfQys8XAROBIMxsOOLAIOK8N+yjSJpTbkkdhUXf38Q3c/Js26ItI\nu1JuSx6Ze/teBjQz5+lgtpAPZbgkel/c72MPjmcsuvecE8OYy39zYRjz7VevDGPu2ueEMKbaiw+Y\n+uw7t4TrOHSnR8KYp5OxNkUt69k/jNnp1ZVhzPpP9QpjOk97K4zZuD4eVMWeO7bomnpLJNfUXy0S\nkeH54pDaOGZCnPuDT4hnm3ppdXx8N7zSI4wZMnROGLPgmVFhzOBBcZ/nP3VwGLNDvzfDmP41L4cx\nC6bHfWZyHMLsSRmCzi66dOzYzsyc2bvtrqmLiEh5UlEXEckRFXURkRxRURcRyREVdRGRHFFRFxHJ\nERV1EZEcadZ3v7TUiC5BQIbvpieeR4MB7BzGvJHhu+d7s1cYM6JT/NHo7gwIY6rZWnT58Kr4kA3I\nMPkBdAoj+g6L17JDhv68OzBDb6rjiUg2dYz38RNxU21qxIiORZb2KbIsdUCGRjIc3v0yPEBqMuzz\njRmGBmRpq3P0mAf6Z1hPpwz96VwVb9eeWfqc5WGU5XhtznDcKZY3sP/+HZg5M8NqKNXgI5E2VNrB\nRyJtJ0tut3tRf18HzLxUD8LmUp/bRyX2uVAl9l99bntt3V9dUxcRyREVdRGRHCmHop7l227Kjfrc\nPiqxz4Uqsf/qc9tr0/6W/Jq6iIi0nnI4UxcRkVaioi4ikiMlLepmdqyZPW9mL5rZhFL2JSszW2Rm\nT5vZPDOLv9W/BMxsipmtMLP5Bbf1MLP7zOyF9PcupexjoUb6W2tmS9L9PM/MPlHKPjaF8rptVFpe\nQ2lyu2RF3cyqgeuA44BBwHgzG1Sq/jTRUe4+3N0zTI1SElOBY+vdNgF4wN0HAg+k/5eLqby/vwA/\nTffzcHe/p5371CzK6zY1lcrKayhBbpfyTH008KK7L3T3TcCtwLgS9ic33H0WsKrezeOAG9K/bwDi\nefzaSSP9rVTK6zZSaXkNpcntUhb1PYDXCv5fnN5W7hy438zmmtm5pe5ME/R296Xp38uA3qXsTEZf\nM7On0pewZfWyugjldfuqxLyGNsxtvVHadIe7+3CSl9fnm9mHS92hpvLkc6zl/lnWXwD9Sb7ebSnw\n49J2J/eU1+2nTXO7lEV9CWz39Yd7preVNXdfkv5eAUwjebldCZabWR+A9PeKEvenKHdf7u5b3X0b\n8CsqZz8rr9tXReU1tH1ul7KozwEGmtm+ZtYJOB2YUcL+hMxsJzPrVvc38DFgfvF7lY0ZwJnp32cC\n00vYl1DdAzV1EpWzn5XX7aui8hraPrdL8n3qAO6+xcy+CtwLVANT3H1BqfqTUW9gmplBsu9udvc/\nl7ZL72dmtwBHAr3MbDEwEZgM3GZmZwOvAKeVrofba6S/R5rZcJKX04uA80rWwSZQXredSstrKE1u\n62sCRERyRG+UiojkiIq6iEiOqKiLiOSIirqISI6oqIuI5IiKuohIjqioi4jkiIq6iEiO/C/zDkyC\nZclMUgAAAABJRU5ErkJggg==\n", "text/plain": [ - "" + "" ] }, "metadata": {}, @@ -1609,7 +1609,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", - "version": "2.7.11" + "version": "2.7.6" } }, "nbformat": 4, diff --git a/openmc/mgxs/mgxs.py b/openmc/mgxs/mgxs.py index 89f12d5ed..5b49005ec 100644 --- a/openmc/mgxs/mgxs.py +++ b/openmc/mgxs/mgxs.py @@ -227,29 +227,30 @@ class MGXS(object): domain_filter = openmc.Filter(self.domain_type, self.domain.id) # Create each Tally needed to compute the multi group cross section - for score, key, filters in zip(self.scores, self.keys, self.filters): - self.tallies[key] = openmc.Tally(name=self.name) - self.tallies[key].scores = [score] - self.tallies[key].estimator = self.estimator - self.tallies[key].filters = [domain_filter] + tally_metadata = zip(self.scores, self.tally_keys, self.filters) + for score, key, filters in tally_metadata: + self._tallies[key] = openmc.Tally(name=self.name) + self._tallies[key].scores = [score] + self._tallies[key].estimator = self.estimator + self._tallies[key].filters = [domain_filter] # If a tally trigger was specified, add it to each tally if self.tally_trigger: trigger_clone = copy.deepcopy(self.tally_trigger) trigger_clone.scores = [score] - self.tallies[key].triggers.append(trigger_clone) + self._tallies[key].triggers.append(trigger_clone) # Add non-domain specific Filters (e.g., 'energy') to the Tally for add_filter in filters: - self.tallies[key].filters.append(add_filter) + self._tallies[key].filters.append(add_filter) # If this is a by-nuclide cross-section, add nuclides to Tally if self.by_nuclide and score != 'flux': all_nuclides = self.get_all_nuclides() for nuclide in all_nuclides: - self.tallies[key].nuclides.append(nuclide) + self._tallies[key].nuclides.append(nuclide) else: - self.tallies[key].nuclides.append('total') + self._tallies[key].nuclides.append('total') return self._tallies @@ -311,7 +312,7 @@ class MGXS(object): def filters(self): group_edges = self.energy_groups.group_edges energy_filter = openmc.Filter('energy', group_edges) - return [[energy_filter] * len(self.scores)] + return [[energy_filter]] * len(self.scores) @property def tally_keys(self): @@ -1544,7 +1545,7 @@ class NuTransportXS(TransportXS): return ['flux', 'total', 'nu-scatter-1'] @property - def keys(self): + def tally_keys(self): return ['flux', 'total', 'scatter-1'] @@ -1674,50 +1675,38 @@ class ScatterMatrixXS(MGXS): @property def scores(self): - return ['flux', 'total', 'nu-scatter-1'] + scores = ['flux'] + + for moment in range(self.legendre_order+1): + scores.append('scatter-{}'.format(moment)) + + if self.correction == 'P0' and self.legendre_order == 0: + scores.append('scatter-1') + + return scores @property - def keys(self): - return ['flux', 'total', 'scatter-1'] + def filters(self): + group_edges = self.energy_groups.group_edges + energy = openmc.Filter('energy', group_edges) + energyout = openmc.Filter('energyout', group_edges) + filters = [[energy]] + + for moment in range(self.legendre_order+1): + filters.append([energy, energyout]) + + if self.correction == 'P0' and self.legendre_order == 0: + filters.append([energyout]) + + return filters @property - def tallies(self): - """Construct the OpenMC tallies needed to compute this cross section. + def tally_keys(self): + return ['flux', 'scatter-0', 'scatter-1'] - This method constructs three analog tallies to compute the 'flux' - and Legendre scattering moment reaction rates in the spatial domain and - energy groups of interest. - - """ - - # Instantiate tallies if they do not exist - if self._tallies is None: - - group_edges = self.energy_groups.group_edges - energy = openmc.Filter('energy', group_edges) - energyout = openmc.Filter('energyout', group_edges) - - # Create lists of scores, filters for each Tally to be created - scores = ['flux'] - filters = [[energy]] - - # Create separate tallies for each moment - for moment in range(self.legendre_order+1): - scores.append('scatter-{}'.format(moment)) - filters.append([energy, energyout]) - - # Append to the lists for the P0 approximation if needed - if self.correction == 'P0' and self.legendre_order == 0: - scores.append('scatter-1') - filters.append([energyout]) - - estimator = 'analog' - keys = scores - - # Initialize the Tallies - self._create_tallies(scores, filters, keys, estimator) - - return self._tallies + @property + def estimator(self): + return 'analog' @property def rxn_rate_tally(self): @@ -1727,7 +1716,7 @@ class ScatterMatrixXS(MGXS): # If using P0 correction subtract scatter-1 from the diagonal if self.correction == 'P0' and self.legendre_order == 0: scatter_p1 = self.tallies['scatter-1'] - scatter_p1 = scatter_p1.get_slice(scores=['scatter-1']) + scatter_p1 = scatter_p1.get_slice(scores=[self.scores[-1]]) energy_filter = self.tallies['scatter-0'].find_filter('energy') energy_filter = copy.deepcopy(energy_filter) scatter_p1 = scatter_p1.diagonalize_filter(energy_filter) @@ -1737,8 +1726,7 @@ class ScatterMatrixXS(MGXS): else: rxn_rate_tally = self.tallies['scatter-0'] for moment in range(1, self.legendre_order+1): - scatter_key = 'scatter-{}'.format(moment) - scatter_pn = self.tallies[scatter_key] + scatter_pn = self.tallies['scatter-{}'.format(moment)] rxn_rate_tally = rxn_rate_tally.merge(scatter_pn) self._rxn_rate_tally = rxn_rate_tally @@ -2188,45 +2176,16 @@ class NuScatterMatrixXS(ScatterMatrixXS): self._rxn_type = 'nu-scatter matrix' @property - def tallies(self): - """Construct the OpenMC tallies needed to compute this cross section. + def scores(self): + scores = ['flux'] - This method constructs three analog tallies to compute the 'flux', - 'nu-scatter' and 'scatter-P1' reaction rates in the spatial domain and - energy groups of interest. + for moment in range(self.legendre_order+1): + scores.append('nu-scatter-{}'.format(moment)) - """ + if self.correction == 'P0' and self.legendre_order == 0: + scores.append('nu-scatter-1') - # Instantiate tallies if they do not exist - if self._tallies is None: - - group_edges = self.energy_groups.group_edges - energy = openmc.Filter('energy', group_edges) - energyout = openmc.Filter('energyout', group_edges) - - # Create lists of scores, filters for each Tally to be created - scores = ['flux'] - filters = [[energy]] - keys = ['flux'] - - # Create separate tallies for each moment - for moment in range(self.legendre_order+1): - scores.append('nu-scatter-{}'.format(moment)) - filters.append([energy, energyout]) - keys.append('scatter-{}'.format(moment)) - - # Append to the lists for the P0 approximation if needed - if self.correction == 'P0' and self.legendre_order == 0: - scores.append('scatter-1') - filters.append([energyout]) - keys.append('scatter-1') - - estimator = 'analog' - - # Intialize the Tallies - self._create_tallies(scores, filters, keys, estimator) - - return self._tallies + return scores class Chi(MGXS): @@ -2238,33 +2197,24 @@ class Chi(MGXS): self._rxn_type = 'chi' @property - def tallies(self): - """Construct the OpenMC tallies needed to compute this cross section. + def scores(self): + return ['nu-fission', 'nu-fission'] - This method constructs two analog tallies to compute 'nu-fission' - reaction rates with 'energy' and 'energyout' filters in the spatial - domain and energy groups of interest. + @property + def filters(self): + # Create the non-domain specific Filters for the Tallies + group_edges = self.energy_groups.group_edges + energyout = openmc.Filter('energyout', group_edges) + energyin = openmc.Filter('energy', [group_edges[0], group_edges[-1]]) + return [[energyin], [energyout]] - """ + @property + def tally_keys(self): + return ['nu-fission-in', 'nu-fission-out'] - # Instantiate tallies if they do not exist - if self._tallies is None: - - # Create a list of scores for each Tally to be created - scores = ['nu-fission', 'nu-fission'] - estimator = 'analog' - keys = ['nu-fission-in', 'nu-fission-out'] - - # Create the non-domain specific Filters for the Tallies - group_edges = self.energy_groups.group_edges - energyout = openmc.Filter('energyout', group_edges) - energyin = openmc.Filter('energy', [group_edges[0], group_edges[-1]]) - filters = [[energyin], [energyout]] - - # Intialize the Tallies - self._create_tallies(scores, filters, keys, estimator) - - return self._tallies + @property + def estimator(self): + return 'analog' @property def rxn_rate_tally(self): diff --git a/tests/test_mgxs_library_condense/inputs_true.dat b/tests/test_mgxs_library_condense/inputs_true.dat index 708ec114e..064981fa9 100644 --- a/tests/test_mgxs_library_condense/inputs_true.dat +++ b/tests/test_mgxs_library_condense/inputs_true.dat @@ -1 +1 @@ -e3834da92fc6ae57ce109621e3f692a186a03820b61332fa9ed898bc07fb8a63484ace095713d5b88196b1d2f1430d2e7b27a505944c7c3027f6365801f58146 \ No newline at end of file +ee40a2b826dea8323249c7261502f8339c78a5dc236e019842cc5244c048d5978fe66e036b86d46b262260556fbd62b19cbb2f0d70325b92b8c40275e75afe4f \ No newline at end of file diff --git a/tests/test_mgxs_library_condense/results_true.dat b/tests/test_mgxs_library_condense/results_true.dat index 5f1d886b0..89e4dbb3e 100644 --- a/tests/test_mgxs_library_condense/results_true.dat +++ b/tests/test_mgxs_library_condense/results_true.dat @@ -1,7 +1,7 @@ material group in nuclide mean std. dev. 0 1 1 total 0.412084 0.02359 material group in nuclide mean std. dev. 0 1 1 total 0.076425 0.003691 material group in group out nuclide moment mean std. dev. moment -0 1 1 1 total P0 0.345643 0.021487 P0 material group out nuclide mean std. dev. +0 1 1 1 total P0 0.345503 0.021465 P0 material group out nuclide mean std. dev. 0 1 1 total 1.0 0.055333 material group in nuclide mean std. dev. 0 2 1 total 0.241262 0.00841 material group in nuclide mean std. dev. 0 2 1 total 0.0 0.0 material group in group out nuclide moment mean std. dev. moment diff --git a/tests/test_mgxs_library_distribcell/inputs_true.dat b/tests/test_mgxs_library_distribcell/inputs_true.dat index 4ab730274..5ffea7f8f 100644 --- a/tests/test_mgxs_library_distribcell/inputs_true.dat +++ b/tests/test_mgxs_library_distribcell/inputs_true.dat @@ -1 +1 @@ -aadb1e94492741c091bff4b5e17634ee327c718fc9fd1f27aa22fe8406fb70f732750dfdceb30eb40b5e4f406061bea6bd5235ba613c3c81009f5857a9051728 \ No newline at end of file +c46381a2d86bd849ca20dc64022ffcf836ba0f236f392bba6335c42559df61d14d09a616bc3a9590d954a5bf099610eb071982296b75b10d1c168cc3e343d383 \ No newline at end of file diff --git a/tests/test_mgxs_library_hdf5/inputs_true.dat b/tests/test_mgxs_library_hdf5/inputs_true.dat index 708ec114e..064981fa9 100644 --- a/tests/test_mgxs_library_hdf5/inputs_true.dat +++ b/tests/test_mgxs_library_hdf5/inputs_true.dat @@ -1 +1 @@ -e3834da92fc6ae57ce109621e3f692a186a03820b61332fa9ed898bc07fb8a63484ace095713d5b88196b1d2f1430d2e7b27a505944c7c3027f6365801f58146 \ No newline at end of file +ee40a2b826dea8323249c7261502f8339c78a5dc236e019842cc5244c048d5978fe66e036b86d46b262260556fbd62b19cbb2f0d70325b92b8c40275e75afe4f \ No newline at end of file diff --git a/tests/test_mgxs_library_hdf5/results_true.dat b/tests/test_mgxs_library_hdf5/results_true.dat index e19b9ffa5..93aceba7d 100644 --- a/tests/test_mgxs_library_hdf5/results_true.dat +++ b/tests/test_mgxs_library_hdf5/results_true.dat @@ -5,9 +5,9 @@ domain=1 type=nu-fission [ 0.02178897 0.71407658] [ 0.00118187 0.04055185] domain=1 type=nu-scatter matrix -[[ 0.3373971 0.00155945] +[[ 0.33724504 0.00155945] [ 0. 0.42205129]] -[[ 0.02303884 0.00051015] +[[ 0.02301463 0.00051015] [ 0. 0.02161702]] domain=1 type=chi [ 1. 0.] diff --git a/tests/test_mgxs_library_no_nuclides/inputs_true.dat b/tests/test_mgxs_library_no_nuclides/inputs_true.dat index 708ec114e..064981fa9 100644 --- a/tests/test_mgxs_library_no_nuclides/inputs_true.dat +++ b/tests/test_mgxs_library_no_nuclides/inputs_true.dat @@ -1 +1 @@ -e3834da92fc6ae57ce109621e3f692a186a03820b61332fa9ed898bc07fb8a63484ace095713d5b88196b1d2f1430d2e7b27a505944c7c3027f6365801f58146 \ No newline at end of file +ee40a2b826dea8323249c7261502f8339c78a5dc236e019842cc5244c048d5978fe66e036b86d46b262260556fbd62b19cbb2f0d70325b92b8c40275e75afe4f \ No newline at end of file diff --git a/tests/test_mgxs_library_no_nuclides/results_true.dat b/tests/test_mgxs_library_no_nuclides/results_true.dat index 4e6d88208..0ad8e04aa 100644 --- a/tests/test_mgxs_library_no_nuclides/results_true.dat +++ b/tests/test_mgxs_library_no_nuclides/results_true.dat @@ -3,7 +3,7 @@ 0 1 2 total 0.861607 0.032349 material group in nuclide mean std. dev. 1 1 1 total 0.021789 0.001182 0 1 2 total 0.714077 0.040552 material group in group out nuclide moment mean std. dev. moment -3 1 1 1 total P0 0.337397 0.023039 P0 +3 1 1 1 total P0 0.337245 0.023015 P0 2 1 1 2 total P0 0.001559 0.000510 P0 1 1 2 1 total P0 0.000000 0.000000 P0 0 1 2 2 total P0 0.422051 0.021617 P0 material group out nuclide mean std. dev. diff --git a/tests/test_mgxs_library_nuclides/inputs_true.dat b/tests/test_mgxs_library_nuclides/inputs_true.dat index 5b7a83037..d2c11978a 100644 --- a/tests/test_mgxs_library_nuclides/inputs_true.dat +++ b/tests/test_mgxs_library_nuclides/inputs_true.dat @@ -1 +1 @@ -f1c203fb7f0b141ee608d7bb9223aa5f7ab84966b6a80525879df730d19179f6c6a1a4bc7038d84e6b366b360f43a1ca17a0af8d02f96eb23e53b93a2445380e \ No newline at end of file +f4abbd7867b0f0d2d9d93ed089c95904541f970522e1ef3a843373b60094ef4571a64a7b5f68efe9e51fd49754bc9e20a3bcc85c0bde3a8224608f8b97c01b85 \ No newline at end of file diff --git a/tests/test_mgxs_library_nuclides/results_true.dat b/tests/test_mgxs_library_nuclides/results_true.dat index ff34c9fff..9cef6fdd8 100644 --- a/tests/test_mgxs_library_nuclides/results_true.dat +++ b/tests/test_mgxs_library_nuclides/results_true.dat @@ -138,7 +138,7 @@ 102 1 1 1 U-234 P0 0.000000 0.000000 P0 103 1 1 1 U-235 P0 0.003226 0.001139 P0 104 1 1 1 U-236 P0 0.001697 0.000923 P0 -105 1 1 1 U-238 P0 0.194620 0.013297 P0 +105 1 1 1 U-238 P0 0.194468 0.013279 P0 106 1 1 1 Np-237 P0 0.000000 0.000000 P0 107 1 1 1 Pu-238 P0 0.000000 0.000000 P0 108 1 1 1 Pu-239 P0 0.001005 0.000477 P0