{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Multigroup Cross Section Generation Part II: Advanced Features\n",
    "This IPython Notebook illustrates the use of the `openmc.mgxs` module to calculate multi-group cross sections for a heterogeneous fuel pin cell geometry. In particular, this Notebook illustrates the following features:\n",
    "\n",
    "* Creation of multi-group cross sections on a **heterogeneous geometry**\n",
    "* Calculation of cross sections on a **nuclide-by-nuclide basis**\n",
    "* The use of **[tally precision triggers](../io_formats/settings.rst#trigger-element)** with multi-group cross sections\n",
    "* Built-in features for **energy condensation** in downstream data processing\n",
    "* The use of the **`openmc.data`** module to plot continuous-energy vs. multi-group cross sections\n",
    "* **Validation** of multi-group cross sections with **[OpenMOC](https://mit-crpg.github.io/OpenMOC/)**\n",
    "\n",
    "**Note:** This Notebook was created using [OpenMOC](https://mit-crpg.github.io/OpenMOC/) to verify the multi-group cross-sections generated by OpenMC. You must install [OpenMOC](https://mit-crpg.github.io/OpenMOC/) on your system in order to run this Notebook in its entirety. In addition, this Notebook illustrates the use of [Pandas](https://pandas.pydata.org/) `DataFrames` to containerize multi-group cross section data."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Generate Input Files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "# plt.style.use('seaborn-dark')\n",
    "\n",
    "import openmoc\n",
    "\n",
    "import openmc\n",
    "import openmc.mgxs as mgxs\n",
    "import openmc.data\n",
    "from openmc.openmoc_compatible import get_openmoc_geometry\n",
    "\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# create a model object to tie together geometry, materials, settings, and tallies\n",
    "model = openmc.Model()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "First we need to define materials that will be used in the problem. We'll create three distinct materials for water, clad and fuel."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# 1.6% enriched fuel\n",
    "fuel = openmc.Material(name='1.6% Fuel')\n",
    "fuel.set_density('g/cm3', 10.31341)\n",
    "fuel.add_nuclide('U235', 3.7503e-4)\n",
    "fuel.add_nuclide('U238', 2.2625e-2)\n",
    "fuel.add_nuclide('O16', 4.6007e-2)\n",
    "\n",
    "# borated water\n",
    "water = openmc.Material(name='Borated Water')\n",
    "water.set_density('g/cm3', 0.740582)\n",
    "water.add_nuclide('H1', 4.9457e-2)\n",
    "water.add_nuclide('O16', 2.4732e-2)\n",
    "\n",
    "# zircaloy\n",
    "zircaloy = openmc.Material(name='Zircaloy')\n",
    "zircaloy.set_density('g/cm3', 6.55)\n",
    "zircaloy.add_nuclide('Zr90', 7.2758e-3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "With our materials, we can now create a `Materials` object that can be exported to an actual XML file."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Instantiate a Materials collection\n",
    "model.materials = openmc.Materials([fuel, water, zircaloy])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now let's move on to the geometry. Our problem will have three regions for the fuel, the clad, and the surrounding coolant. The first step is to create the bounding surfaces -- in this case two cylinders and six reflective planes."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create cylinders for the fuel and clad\n",
    "fuel_outer_radius = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.39218)\n",
    "clad_outer_radius = openmc.ZCylinder(x0=0.0, y0=0.0, r=0.45720)\n",
    "\n",
    "# Create box to surround the geometry\n",
    "box = openmc.model.RectangularPrism(1.26, 1.26, boundary_type='reflective')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "With the surfaces defined, we can now create cells that are defined by intersections of half-spaces created by the surfaces."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create a Universe to encapsulate a fuel pin\n",
    "pin_cell_universe = openmc.Universe(name='1.6% Fuel Pin')\n",
    "\n",
    "# Create fuel Cell\n",
    "fuel_cell = openmc.Cell(name='1.6% Fuel')\n",
    "fuel_cell.fill = fuel\n",
    "fuel_cell.region = -fuel_outer_radius\n",
    "pin_cell_universe.add_cell(fuel_cell)\n",
    "\n",
    "# Create a clad Cell\n",
    "clad_cell = openmc.Cell(name='1.6% Clad')\n",
    "clad_cell.fill = zircaloy\n",
    "clad_cell.region = +fuel_outer_radius & -clad_outer_radius\n",
    "pin_cell_universe.add_cell(clad_cell)\n",
    "\n",
    "# Create a moderator Cell\n",
    "moderator_cell = openmc.Cell(name='1.6% Moderator')\n",
    "moderator_cell.fill = water\n",
    "moderator_cell.region = +clad_outer_radius & box\n",
    "pin_cell_universe.add_cell(moderator_cell)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We now must create a geometry with the pin cell universe and export it to XML."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create Geometry and set root Universe\n",
    "model.geometry = openmc.Geometry(pin_cell_universe)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, we must define simulation parameters. In this case, we will use 10 inactive batches and 40 active batches each with 10,000 particles."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# OpenMC simulation parameters\n",
    "batches = 50\n",
    "inactive = 10\n",
    "particles = 10000\n",
    "\n",
    "# Instantiate a Settings object\n",
    "settings = openmc.Settings()\n",
    "settings.batches = batches\n",
    "settings.inactive = inactive\n",
    "settings.particles = particles\n",
    "settings.output = {'tallies': True}\n",
    "\n",
    "# Create an initial uniform spatial source distribution over fissionable zones\n",
    "bounds = [-0.63, -0.63, -0.63, 0.63, 0.63, 0.63]\n",
    "uniform_dist = openmc.stats.Box(bounds[:3], bounds[3:], only_fissionable=True)\n",
    "settings.source = openmc.IndependentSource(space=uniform_dist)\n",
    "\n",
    "# Activate tally precision triggers\n",
    "settings.trigger_active = True\n",
    "settings.trigger_max_batches = settings.batches * 4\n",
    "\n",
    "model.settings = settings"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we are finally ready to make use of the `openmc.mgxs` module to generate multi-group cross sections! First, let's define \"coarse\" 2-group and \"fine\" 8-group structures using the built-in `EnergyGroups` class."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Instantiate a \"coarse\" 2-group EnergyGroups object\n",
    "coarse_groups = mgxs.EnergyGroups([0., 0.625, 20.0e6])\n",
    "\n",
    "# Instantiate a \"fine\" 8-group EnergyGroups object\n",
    "fine_groups = mgxs.EnergyGroups([0., 0.058, 0.14, 0.28,\n",
    "                                 0.625, 4.0, 5.53e3, 821.0e3, 20.0e6])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we will instantiate a variety of `MGXS` objects needed to run an OpenMOC simulation to verify the accuracy of our cross sections. In particular, we define transport, fission, nu-fission, nu-scatter and chi cross sections for each of the three cells in the fuel pin with the 8-group structure as our energy groups."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Extract all Cells filled by Materials\n",
    "openmc_cells = model.geometry.get_all_material_cells().values()\n",
    "\n",
    "# Create dictionary to store multi-group cross sections for all cells\n",
    "xs_library = {}\n",
    "\n",
    "# Instantiate 8-group cross sections for each cell\n",
    "for cell in openmc_cells:\n",
    "    xs_library[cell.id] = {}\n",
    "    xs_library[cell.id]['transport']  = mgxs.TransportXS(energy_groups=fine_groups)\n",
    "    xs_library[cell.id]['fission'] = mgxs.FissionXS(energy_groups=fine_groups)\n",
    "    xs_library[cell.id]['nu-fission'] = mgxs.FissionXS(energy_groups=fine_groups, nu=True)\n",
    "    xs_library[cell.id]['nu-scatter'] = mgxs.ScatterMatrixXS(energy_groups=fine_groups, nu=True)\n",
    "    xs_library[cell.id]['chi'] = mgxs.Chi(energy_groups=fine_groups)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, we showcase the use of OpenMC's [tally precision trigger](../io_formats/settings.rst#trigger-element) feature in conjunction with the `openmc.mgxs` module. In particular, we will assign a tally trigger of 1E-2 on the standard deviation for each of the tallies used to compute multi-group cross sections."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create a tally trigger for +/- 0.01 on each tally used to compute the multi-group cross sections\n",
    "tally_trigger = openmc.Trigger('std_dev', 1e-2)\n",
    "\n",
    "# Add the tally trigger to each of the multi-group cross section tallies\n",
    "for cell in openmc_cells:\n",
    "    for mgxs_type in xs_library[cell.id]:\n",
    "        xs_library[cell.id][mgxs_type].tally_trigger = tally_trigger"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now, we must loop over all cells to set the cross section domains to the various cells - fuel, clad and moderator - included in the geometry. In addition, we will set each cross section to tally cross sections on a per-nuclide basis through the use of the `MGXS` class' boolean `by_nuclide` instance attribute.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Instantiate an empty Tallies object\n",
    "tallies = openmc.Tallies()\n",
    "\n",
    "# Iterate over all cells and cross section types\n",
    "for cell in openmc_cells:\n",
    "    for rxn_type in xs_library[cell.id]:\n",
    "\n",
    "        # Set the cross sections domain to the cell\n",
    "        xs_library[cell.id][rxn_type].domain = cell\n",
    "\n",
    "        # Tally cross sections by nuclide\n",
    "        xs_library[cell.id][rxn_type].by_nuclide = True\n",
    "\n",
    "        # Add OpenMC tallies to the tallies file for XML generation\n",
    "        for tally in xs_library[cell.id][rxn_type].tallies.values():\n",
    "            tallies.append(tally, merge=True)\n",
    "\n",
    "model.tallies = tallies"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we a have a complete set of inputs, so we can go ahead and run our simulation."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/pshriwise/.pyenv/versions/3.9.1/lib/python3.9/site-packages/openmc/mixin.py:70: IDWarning: Another Filter instance already exists with id=53.\n",
      "  warn(msg, IDWarning)\n",
      "/home/pshriwise/.pyenv/versions/3.9.1/lib/python3.9/site-packages/openmc/mixin.py:70: IDWarning: Another Filter instance already exists with id=21.\n",
      "  warn(msg, IDWarning)\n",
      "/home/pshriwise/.pyenv/versions/3.9.1/lib/python3.9/site-packages/openmc/mixin.py:70: IDWarning: Another Filter instance already exists with id=2.\n",
      "  warn(msg, IDWarning)\n",
      "/home/pshriwise/.pyenv/versions/3.9.1/lib/python3.9/site-packages/openmc/mixin.py:70: IDWarning: Another Filter instance already exists with id=3.\n",
      "  warn(msg, IDWarning)\n",
      "/home/pshriwise/.pyenv/versions/3.9.1/lib/python3.9/site-packages/openmc/mixin.py:70: IDWarning: Another Filter instance already exists with id=4.\n",
      "  warn(msg, IDWarning)\n",
      "/home/pshriwise/.pyenv/versions/3.9.1/lib/python3.9/site-packages/openmc/mixin.py:70: IDWarning: Another Filter instance already exists with id=41.\n",
      "  warn(msg, IDWarning)\n",
      "/home/pshriwise/.pyenv/versions/3.9.1/lib/python3.9/site-packages/openmc/mixin.py:70: IDWarning: Another Filter instance already exists with id=15.\n",
      "  warn(msg, IDWarning)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                                %%%%%%%%%%%%%%%\n",
      "                           %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                      %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                                    %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                                     %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                 ###############      %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ##################     %%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ###################     %%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ####################     %%%%%%%%%%%%%%%%%%%%%%\n",
      "                #####################     %%%%%%%%%%%%%%%%%%%%%\n",
      "                ######################     %%%%%%%%%%%%%%%%%%%%\n",
      "                #######################     %%%%%%%%%%%%%%%%%%\n",
      "                 #######################     %%%%%%%%%%%%%%%%%\n",
      "                 ######################     %%%%%%%%%%%%%%%%%\n",
      "                  ####################     %%%%%%%%%%%%%%%%%\n",
      "                    #################     %%%%%%%%%%%%%%%%%\n",
      "                     ###############     %%%%%%%%%%%%%%%%\n",
      "                       ############     %%%%%%%%%%%%%%%\n",
      "                          ########     %%%%%%%%%%%%%%\n",
      "                                      %%%%%%%%%%%\n",
      "\n",
      "                 | The OpenMC Monte Carlo Code\n",
      "       Copyright | 2011-2023 MIT, UChicago Argonne LLC, and contributors\n",
      "         License | https://docs.openmc.org/en/latest/license.html\n",
      "         Version | 0.13.3\n",
      "        Git SHA1 | 50e39a4e20dc9e0f3d7ccf07333f6a5e6c797c8c\n",
      "       Date/Time | 2023-11-07 12:05:29\n",
      "  OpenMP Threads | 32\n",
      "\n",
      " Reading settings XML file...\n",
      " Reading cross sections XML file...\n",
      " Reading materials XML file...\n",
      " Reading geometry XML file...\n",
      " Reading U235 from /home/pshriwise/data/xs/openmc/nndc_hdf5/U235.h5\n",
      " Reading U238 from /home/pshriwise/data/xs/openmc/nndc_hdf5/U238.h5\n",
      " Reading O16 from /home/pshriwise/data/xs/openmc/nndc_hdf5/O16.h5\n",
      " Reading H1 from /home/pshriwise/data/xs/openmc/nndc_hdf5/H1.h5\n",
      " Reading Zr90 from /home/pshriwise/data/xs/openmc/nndc_hdf5/Zr90.h5\n",
      " Minimum neutron data temperature: 294 K\n",
      " Maximum neutron data temperature: 294 K\n",
      " Reading tallies XML file...\n",
      " Preparing distributed cell instances...\n",
      " Reading plot XML file...\n",
      " Writing summary.h5 file...\n",
      " Maximum neutron transport energy: 20000000 eV for U235\n",
      " Initializing source particles...\n",
      "\n",
      " ====================>     K EIGENVALUE SIMULATION     <====================\n",
      "\n",
      "  Bat./Gen.      k            Average k\n",
      "  =========   ========   ====================\n",
      "        1/1    1.22394\n",
      "        2/1    1.22784\n",
      "        3/1    1.25162\n",
      "        4/1    1.18933\n",
      "        5/1    1.24249\n",
      "        6/1    1.22774\n",
      "        7/1    1.22083\n",
      "        8/1    1.21236\n",
      "        9/1    1.24216\n",
      "       10/1    1.20525\n",
      "       11/1    1.22801\n",
      "       12/1    1.24750    1.23775 +/- 0.00975\n",
      "       13/1    1.23962    1.23838 +/- 0.00566\n",
      "       14/1    1.24706    1.24055 +/- 0.00455\n",
      "       15/1    1.22279    1.23700 +/- 0.00500\n",
      "       16/1    1.18058    1.22759 +/- 0.01025\n",
      "       17/1    1.21435    1.22570 +/- 0.00887\n",
      "       18/1    1.23701    1.22712 +/- 0.00781\n",
      "       19/1    1.23145    1.22760 +/- 0.00690\n",
      "       20/1    1.21727    1.22656 +/- 0.00626\n",
      "       21/1    1.22156    1.22611 +/- 0.00568\n",
      "       22/1    1.25529    1.22854 +/- 0.00573\n",
      "       23/1    1.23090    1.22872 +/- 0.00527\n",
      "       24/1    1.23231    1.22898 +/- 0.00489\n",
      "       25/1    1.23425    1.22933 +/- 0.00456\n",
      "       26/1    1.22500    1.22906 +/- 0.00428\n",
      "       27/1    1.21126    1.22801 +/- 0.00415\n",
      "       28/1    1.21880    1.22750 +/- 0.00395\n",
      "       29/1    1.25221    1.22880 +/- 0.00395\n",
      "       30/1    1.23143    1.22893 +/- 0.00375\n",
      "       31/1    1.24257    1.22958 +/- 0.00363\n",
      "       32/1    1.23452    1.22981 +/- 0.00347\n",
      "       33/1    1.20175    1.22859 +/- 0.00353\n",
      "       34/1    1.22931    1.22862 +/- 0.00338\n",
      "       35/1    1.23860    1.22902 +/- 0.00327\n",
      "       36/1    1.18381    1.22728 +/- 0.00359\n",
      "       37/1    1.20017    1.22627 +/- 0.00360\n",
      "       38/1    1.23447    1.22657 +/- 0.00348\n",
      "       39/1    1.20672    1.22588 +/- 0.00342\n",
      "       40/1    1.20439    1.22516 +/- 0.00339\n",
      "       41/1    1.21524    1.22484 +/- 0.00329\n",
      "       42/1    1.21467    1.22453 +/- 0.00320\n",
      "       43/1    1.21988    1.22439 +/- 0.00311\n",
      "       44/1    1.21610    1.22414 +/- 0.00302\n",
      "       45/1    1.20887    1.22371 +/- 0.00297\n",
      "       46/1    1.23855    1.22412 +/- 0.00291\n",
      "       47/1    1.21390    1.22384 +/- 0.00285\n",
      "       48/1    1.21694    1.22366 +/- 0.00278\n",
      "       49/1    1.24311    1.22416 +/- 0.00275\n",
      "       50/1    1.22343    1.22414 +/- 0.00268\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      " Creating state point statepoint.050.h5...\n",
      "       51/1    1.22495    1.22416 +/- 0.00261\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       52/1    1.21607    1.22397 +/- 0.00256\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       53/1    1.21531    1.22377 +/- 0.00251\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       54/1    1.22888    1.22388 +/- 0.00245\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       55/1    1.20765    1.22352 +/- 0.00242\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       56/1    1.22969    1.22366 +/- 0.00237\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       57/1    1.24581    1.22413 +/- 0.00237\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       58/1    1.20849    1.22380 +/- 0.00234\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       59/1    1.23903    1.22411 +/- 0.00232\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       60/1    1.20733    1.22378 +/- 0.00229\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       61/1    1.23517    1.22400 +/- 0.00226\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       62/1    1.22390    1.22400 +/- 0.00222\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       63/1    1.22695    1.22405 +/- 0.00217\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       64/1    1.20634    1.22373 +/- 0.00216\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       65/1    1.22490    1.22375 +/- 0.00212\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       66/1    1.21709    1.22363 +/- 0.00208\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       67/1    1.22494    1.22365 +/- 0.00205\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       68/1    1.21584    1.22352 +/- 0.00202\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       69/1    1.21189    1.22332 +/- 0.00199\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       70/1    1.21044    1.22310 +/- 0.00197\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       71/1    1.21279    1.22294 +/- 0.00194\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       72/1    1.24979    1.22337 +/- 0.00196\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       73/1    1.21196    1.22319 +/- 0.00194\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       74/1    1.23868    1.22343 +/- 0.00192\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       75/1    1.22632    1.22347 +/- 0.00189\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       76/1    1.20763    1.22323 +/- 0.00188\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       77/1    1.20654    1.22299 +/- 0.00187\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       78/1    1.21582    1.22288 +/- 0.00184\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       79/1    1.22156    1.22286 +/- 0.00182\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       80/1    1.26424    1.22345 +/- 0.00189\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       81/1    1.23856    1.22366 +/- 0.00187\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       82/1    1.22635    1.22370 +/- 0.00185\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       83/1    1.23682    1.22388 +/- 0.00183\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       84/1    1.21823    1.22381 +/- 0.00181\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       85/1    1.24048    1.22403 +/- 0.00180\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       86/1    1.22127    1.22399 +/- 0.00177\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       87/1    1.22862    1.22405 +/- 0.00175\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       88/1    1.24011    1.22426 +/- 0.00174\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       89/1    1.22234    1.22423 +/- 0.00172\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       90/1    1.23789    1.22440 +/- 0.00170\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       91/1    1.21690    1.22431 +/- 0.00169\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       92/1    1.23422    1.22443 +/- 0.00167\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       93/1    1.21595    1.22433 +/- 0.00165\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       94/1    1.23509    1.22446 +/- 0.00164\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       95/1    1.19459    1.22411 +/- 0.00166\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       96/1    1.22999    1.22417 +/- 0.00164\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       97/1    1.22043    1.22413 +/- 0.00162\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       98/1    1.22169    1.22410 +/- 0.00160\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "       99/1    1.22003    1.22406 +/- 0.00158\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      100/1    1.20257    1.22382 +/- 0.00158\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      101/1    1.21000    1.22367 +/- 0.00157\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      102/1    1.23888    1.22383 +/- 0.00157\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      103/1    1.21211    1.22371 +/- 0.00155\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      104/1    1.20659    1.22352 +/- 0.00155\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      105/1    1.21947    1.22348 +/- 0.00153\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      106/1    1.23843    1.22364 +/- 0.00152\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      107/1    1.22764    1.22368 +/- 0.00151\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      108/1    1.23291    1.22377 +/- 0.00150\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      109/1    1.22325    1.22377 +/- 0.00148\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      110/1    1.21052    1.22364 +/- 0.00147\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      111/1    1.23348    1.22373 +/- 0.00146\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      112/1    1.22456    1.22374 +/- 0.00145\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      113/1    1.20578    1.22357 +/- 0.00144\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      114/1    1.22042    1.22354 +/- 0.00143\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      115/1    1.22143    1.22352 +/- 0.00142\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      116/1    1.23161    1.22359 +/- 0.00140\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      117/1    1.22951    1.22365 +/- 0.00139\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      118/1    1.21103    1.22353 +/- 0.00138\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      119/1    1.24296    1.22371 +/- 0.00138\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      120/1    1.22464    1.22372 +/- 0.00137\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      121/1    1.24168    1.22388 +/- 0.00137\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      122/1    1.22327    1.22387 +/- 0.00136\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      123/1    1.22344    1.22387 +/- 0.00134\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      124/1    1.21236    1.22377 +/- 0.00134\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      125/1    1.25017    1.22400 +/- 0.00134\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      126/1    1.22034    1.22397 +/- 0.00133\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      127/1    1.20737    1.22383 +/- 0.00133\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      128/1    1.22039    1.22380 +/- 0.00132\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      129/1    1.19783    1.22358 +/- 0.00132\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      130/1    1.20356    1.22341 +/- 0.00132\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      131/1    1.24384    1.22358 +/- 0.00132\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      132/1    1.22925    1.22363 +/- 0.00131\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      133/1    1.22562    1.22364 +/- 0.00130\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      134/1    1.21676    1.22359 +/- 0.00129\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      135/1    1.20402    1.22343 +/- 0.00129\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      136/1    1.21821    1.22339 +/- 0.00128\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      137/1    1.25992    1.22368 +/- 0.00131\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      138/1    1.21551    1.22361 +/- 0.00130\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      139/1    1.21593    1.22355 +/- 0.00129\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      140/1    1.20545    1.22341 +/- 0.00129\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      141/1    1.23441    1.22350 +/- 0.00128\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      142/1    1.21756    1.22345 +/- 0.00127\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      143/1    1.24155    1.22359 +/- 0.00127\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      144/1    1.25159    1.22380 +/- 0.00127\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      145/1    1.21877    1.22376 +/- 0.00127\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      146/1    1.21616    1.22371 +/- 0.00126\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      147/1    1.23678    1.22380 +/- 0.00125\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      148/1    1.20358    1.22365 +/- 0.00125\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      149/1    1.21250    1.22357 +/- 0.00125\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      150/1    1.20860    1.22347 +/- 0.00124\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      151/1    1.22502    1.22348 +/- 0.00123\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      152/1    1.21802    1.22344 +/- 0.00122\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      153/1    1.24871    1.22362 +/- 0.00123\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      154/1    1.22950    1.22366 +/- 0.00122\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      155/1    1.19887    1.22349 +/- 0.00122\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      156/1    1.21147    1.22340 +/- 0.00122\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      157/1    1.20278    1.22326 +/- 0.00122\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      158/1    1.21427    1.22320 +/- 0.00121\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      159/1    1.21654    1.22316 +/- 0.00120\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      160/1    1.24186    1.22328 +/- 0.00120\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      161/1    1.21032    1.22320 +/- 0.00120\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      162/1    1.19739    1.22303 +/- 0.00120\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      163/1    1.22975    1.22307 +/- 0.00119\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      164/1    1.21494    1.22302 +/- 0.00119\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      165/1    1.20907    1.22293 +/- 0.00118\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      166/1    1.20143    1.22279 +/- 0.00118\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      167/1    1.22043    1.22278 +/- 0.00118\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      168/1    1.22880    1.22281 +/- 0.00117\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      169/1    1.24968    1.22298 +/- 0.00117\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      170/1    1.22392    1.22299 +/- 0.00117\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      171/1    1.23658    1.22307 +/- 0.00116\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      172/1    1.22377    1.22308 +/- 0.00116\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      173/1    1.21977    1.22306 +/- 0.00115\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      174/1    1.21514    1.22301 +/- 0.00114\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      175/1    1.20258    1.22288 +/- 0.00114\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      176/1    1.21552    1.22284 +/- 0.00114\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      177/1    1.20676    1.22274 +/- 0.00113\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      178/1    1.23262    1.22280 +/- 0.00113\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      179/1    1.20232    1.22268 +/- 0.00113\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      180/1    1.22200    1.22268 +/- 0.00112\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      181/1    1.21723    1.22265 +/- 0.00112\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      182/1    1.22808    1.22268 +/- 0.00111\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      183/1    1.22492    1.22269 +/- 0.00110\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      184/1    1.22204    1.22269 +/- 0.00110\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      185/1    1.22431    1.22270 +/- 0.00109\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      186/1    1.22800    1.22273 +/- 0.00109\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      187/1    1.22258    1.22273 +/- 0.00108\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      188/1    1.24411    1.22285 +/- 0.00108\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      189/1    1.23613    1.22292 +/- 0.00108\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      190/1    1.22197    1.22291 +/- 0.00107\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      191/1    1.21940    1.22289 +/- 0.00106\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      192/1    1.21190    1.22283 +/- 0.00106\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      193/1    1.21577    1.22280 +/- 0.00106\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      194/1    1.22491    1.22281 +/- 0.00105\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      195/1    1.21509    1.22277 +/- 0.00104\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      196/1    1.22257    1.22276 +/- 0.00104\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      197/1    1.21065    1.22270 +/- 0.00104\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      198/1    1.23695    1.22278 +/- 0.00103\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      199/1    1.21822    1.22275 +/- 0.00103\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      "      200/1    1.21757    1.22272 +/- 0.00102\n",
      " Triggers unsatisfied, no result tallied for score total in tally 23\n",
      " The estimated number of batches is -2147483637\n",
      " Creating state point statepoint.200.h5...\n",
      "\n",
      " =======================>     TIMING STATISTICS     <=======================\n",
      "\n",
      " Total time for initialization     = 5.0144e-01 seconds\n",
      "   Reading cross sections          = 4.5672e-01 seconds\n",
      " Total time in simulation          = 8.5119e+01 seconds\n",
      "   Time in transport only          = 3.2873e+01 seconds\n",
      "   Time in inactive batches        = 7.6939e-01 seconds\n",
      "   Time in active batches          = 8.4350e+01 seconds\n",
      "   Time synchronizing fission bank = 1.9959e-01 seconds\n",
      "     Sampling source sites         = 1.7881e-01 seconds\n",
      "     SEND/RECV source sites        = 2.0610e-02 seconds\n",
      "   Time accumulating tallies       = 5.1680e+01 seconds\n",
      "   Time writing statepoints        = 5.8805e-02 seconds\n",
      " Total time for finalization       = 4.1537e-03 seconds\n",
      " Total time elapsed                = 8.5672e+01 seconds\n",
      " Calculation Rate (inactive)       = 129973 particles/second\n",
      " Calculation Rate (active)         = 22525.3 particles/second\n",
      "\n",
      " ============================>     RESULTS     <============================\n",
      "\n",
      " k-effective (Collision)     = 1.22328 +/- 0.00091\n",
      " k-effective (Track-length)  = 1.22272 +/- 0.00102\n",
      " k-effective (Absorption)    = 1.22396 +/- 0.00079\n",
      " Combined k-effective        = 1.22359 +/- 0.00070\n",
      " Leakage Fraction            = 0.00000 +/- 0.00000\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# Run OpenMC\n",
    "sp_file = model.run()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Tally Data Processing"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Our simulation ran successfully and created statepoint and summary output files. We begin our analysis by instantiating a `StatePoint` object. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Load the last statepoint file\n",
    "sp = openmc.StatePoint(sp_file)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The statepoint is now ready to be analyzed by our multi-group cross sections. We simply have to load the tallies from the `StatePoint` into each object as follows and our `MGXS` objects will compute the cross sections for us under-the-hood."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Iterate over all cells and cross section types\n",
    "for cell in openmc_cells:\n",
    "    for rxn_type in xs_library[cell.id]:\n",
    "        xs_library[cell.id][rxn_type].load_from_statepoint(sp)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "That's it! Our multi-group cross sections are now ready for the big spotlight. This time we have cross sections in three distinct spatial zones - fuel, clad and moderator - on a per-nuclide basis."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Extracting and Storing MGXS Data"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's first inspect one of our cross sections by printing it to the screen as a microscopic cross section in units of barns."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Multi-Group XS\n",
      "\tReaction Type  =\tnu-fission\n",
      "\tDomain Type    =\tcell\n",
      "\tDomain ID      =\t1\n",
      "\tNuclide        =\tU235\n",
      "\tCross Sections [barns]:\n",
      "            Group 1 [821000.0   - 20000000.0eV]:\t3.30e+00 +/- 1.30e-01%\n",
      "            Group 2 [5530.0     - 821000.0  eV]:\t3.97e+00 +/- 8.86e-02%\n",
      "            Group 3 [4.0        - 5530.0    eV]:\t5.50e+01 +/- 1.26e-01%\n",
      "            Group 4 [0.625      - 4.0       eV]:\t8.84e+01 +/- 1.94e-01%\n",
      "            Group 5 [0.28       - 0.625     eV]:\t2.90e+02 +/- 2.66e-01%\n",
      "            Group 6 [0.14       - 0.28      eV]:\t4.49e+02 +/- 2.50e-01%\n",
      "            Group 7 [0.058      - 0.14      eV]:\t6.87e+02 +/- 1.76e-01%\n",
      "            Group 8 [0.0        - 0.058     eV]:\t1.44e+03 +/- 1.62e-01%\n",
      "\n",
      "\tNuclide        =\tU238\n",
      "\tCross Sections [barns]:\n",
      "            Group 1 [821000.0   - 20000000.0eV]:\t1.06e+00 +/- 1.57e-01%\n",
      "            Group 2 [5530.0     - 821000.0  eV]:\t1.21e-03 +/- 1.56e-01%\n",
      "            Group 3 [4.0        - 5530.0    eV]:\t5.81e-04 +/- 2.04e+00%\n",
      "            Group 4 [0.625      - 4.0       eV]:\t6.54e-06 +/- 1.72e-01%\n",
      "            Group 5 [0.28       - 0.625     eV]:\t1.07e-05 +/- 2.64e-01%\n",
      "            Group 6 [0.14       - 0.28      eV]:\t1.55e-05 +/- 2.51e-01%\n",
      "            Group 7 [0.058      - 0.14      eV]:\t2.30e-05 +/- 1.76e-01%\n",
      "            Group 8 [0.0        - 0.058     eV]:\t4.24e-05 +/- 1.61e-01%\n",
      "\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "nufission = xs_library[fuel_cell.id]['nu-fission']\n",
    "nufission.print_xs(xs_type='micro', nuclides=['U235', 'U238'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Our multi-group cross sections are capable of summing across all nuclides to provide us with macroscopic cross sections as well."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Multi-Group XS\n",
      "\tReaction Type  =\tnu-fission\n",
      "\tDomain Type    =\tcell\n",
      "\tDomain ID      =\t1\n",
      "\tCross Sections [cm^-1]:\n",
      "            Group 1 [821000.0   - 20000000.0eV]:\t2.52e-02 +/- 1.50e-01%\n",
      "            Group 2 [5530.0     - 821000.0  eV]:\t1.51e-03 +/- 8.70e-02%\n",
      "            Group 3 [4.0        - 5530.0    eV]:\t2.07e-02 +/- 1.26e-01%\n",
      "            Group 4 [0.625      - 4.0       eV]:\t3.32e-02 +/- 1.94e-01%\n",
      "            Group 5 [0.28       - 0.625     eV]:\t1.09e-01 +/- 2.66e-01%\n",
      "            Group 6 [0.14       - 0.28      eV]:\t1.69e-01 +/- 2.50e-01%\n",
      "            Group 7 [0.058      - 0.14      eV]:\t2.58e-01 +/- 1.76e-01%\n",
      "            Group 8 [0.0        - 0.058     eV]:\t5.40e-01 +/- 1.62e-01%\n",
      "\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "nufission = xs_library[fuel_cell.id]['nu-fission']\n",
    "nufission.print_xs(xs_type='macro', nuclides='sum')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Although a printed report is nice, it is not scalable or flexible. Let's extract the microscopic cross section data for the moderator as a [Pandas](https://pandas.pydata.org/) `DataFrame` ."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>cell</th>\n",
       "      <th>group in</th>\n",
       "      <th>group out</th>\n",
       "      <th>nuclide</th>\n",
       "      <th>mean</th>\n",
       "      <th>std. dev.</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>126</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>H1</td>\n",
       "      <td>0.233843</td>\n",
       "      <td>0.002374</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>127</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>O16</td>\n",
       "      <td>1.563140</td>\n",
       "      <td>0.003891</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>2</td>\n",
       "      <td>H1</td>\n",
       "      <td>1.590210</td>\n",
       "      <td>0.001807</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>125</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>2</td>\n",
       "      <td>O16</td>\n",
       "      <td>0.285848</td>\n",
       "      <td>0.001017</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>122</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>3</td>\n",
       "      <td>H1</td>\n",
       "      <td>0.010769</td>\n",
       "      <td>0.000138</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>3</td>\n",
       "      <td>O16</td>\n",
       "      <td>0.000003</td>\n",
       "      <td>0.000003</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>120</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>4</td>\n",
       "      <td>H1</td>\n",
       "      <td>0.000012</td>\n",
       "      <td>0.000004</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>121</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>4</td>\n",
       "      <td>O16</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>118</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>5</td>\n",
       "      <td>H1</td>\n",
       "      <td>0.000002</td>\n",
       "      <td>0.000002</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>119</th>\n",
       "      <td>3</td>\n",
       "      <td>1</td>\n",
       "      <td>5</td>\n",
       "      <td>O16</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     cell  group in  group out nuclide      mean  std. dev.\n",
       "126     3         1          1      H1  0.233843   0.002374\n",
       "127     3         1          1     O16  1.563140   0.003891\n",
       "124     3         1          2      H1  1.590210   0.001807\n",
       "125     3         1          2     O16  0.285848   0.001017\n",
       "122     3         1          3      H1  0.010769   0.000138\n",
       "123     3         1          3     O16  0.000003   0.000003\n",
       "120     3         1          4      H1  0.000012   0.000004\n",
       "121     3         1          4     O16  0.000000   0.000000\n",
       "118     3         1          5      H1  0.000002   0.000002\n",
       "119     3         1          5     O16  0.000000   0.000000"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "nuscatter = xs_library[moderator_cell.id]['nu-scatter']\n",
    "df = nuscatter.get_pandas_dataframe(xs_type='micro')\n",
    "df.head(10)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, we illustate how one can easily take multi-group cross sections and condense them down to a coarser energy group structure. The `MGXS` class includes a `get_condensed_xs(...)` method which takes an `EnergyGroups` parameter with a coarse(r) group structure and returns a new `MGXS` condensed to the coarse groups. We illustrate this process below using the 2-group structure created earlier."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Extract the 8-group transport cross section for the fuel\n",
    "fine_xs = xs_library[fuel_cell.id]['transport']\n",
    "\n",
    "# Condense to the 2-group structure\n",
    "condensed_xs = fine_xs.get_condensed_xs(coarse_groups)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Group condensation is as simple as that! We now have a new coarse 2-group `TransportXS` in addition to our original 8-group `TransportXS`. Let's inspect the 2-group `TransportXS` by printing it to the screen and extracting a Pandas `DataFrame` as we have already learned how to do."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Multi-Group XS\n",
      "\tReaction Type  =\ttransport\n",
      "\tDomain Type    =\tcell\n",
      "\tDomain ID      =\t1\n",
      "\tNuclide        =\tU235\n",
      "\tCross Sections [cm^-1]:\n",
      "            Group 1 [0.625      - 20000000.0eV]:\t7.79e-03 +/- 1.22e-01%\n",
      "            Group 2 [0.0        - 0.625     eV]:\t1.82e-01 +/- 1.12e-01%\n",
      "\n",
      "\tNuclide        =\tU238\n",
      "\tCross Sections [cm^-1]:\n",
      "            Group 1 [0.625      - 20000000.0eV]:\t2.17e-01 +/- 7.27e-02%\n",
      "            Group 2 [0.0        - 0.625     eV]:\t2.53e-01 +/- 1.11e-01%\n",
      "\n",
      "\tNuclide        =\tO16\n",
      "\tCross Sections [cm^-1]:\n",
      "            Group 1 [0.625      - 20000000.0eV]:\t1.45e-01 +/- 6.93e-02%\n",
      "            Group 2 [0.0        - 0.625     eV]:\t1.74e-01 +/- 1.22e-01%\n",
      "\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "condensed_xs.print_xs()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>cell</th>\n",
       "      <th>group in</th>\n",
       "      <th>nuclide</th>\n",
       "      <th>mean</th>\n",
       "      <th>std. dev.</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>U235</td>\n",
       "      <td>20.780474</td>\n",
       "      <td>0.025302</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>U238</td>\n",
       "      <td>9.584864</td>\n",
       "      <td>0.006964</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>O16</td>\n",
       "      <td>3.159245</td>\n",
       "      <td>0.002191</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1</td>\n",
       "      <td>2</td>\n",
       "      <td>U235</td>\n",
       "      <td>485.324537</td>\n",
       "      <td>0.544794</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>2</td>\n",
       "      <td>U238</td>\n",
       "      <td>11.195219</td>\n",
       "      <td>0.012455</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>1</td>\n",
       "      <td>2</td>\n",
       "      <td>O16</td>\n",
       "      <td>3.788412</td>\n",
       "      <td>0.004606</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   cell  group in nuclide        mean  std. dev.\n",
       "3     1         1    U235   20.780474   0.025302\n",
       "4     1         1    U238    9.584864   0.006964\n",
       "5     1         1     O16    3.159245   0.002191\n",
       "0     1         2    U235  485.324537   0.544794\n",
       "1     1         2    U238   11.195219   0.012455\n",
       "2     1         2     O16    3.788412   0.004606"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df = condensed_xs.get_pandas_dataframe(xs_type='micro')\n",
    "df"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Verification with OpenMOC"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now, let's verify our cross sections using OpenMOC. First, we construct an equivalent OpenMOC geometry."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Create an OpenMOC Geometry from the OpenMC Geometry\n",
    "openmoc_geometry = get_openmoc_geometry(sp.summary.geometry)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, we we can inject the multi-group cross sections into the equivalent fuel pin cell OpenMOC geometry."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Get all OpenMOC cells in the gometry\n",
    "openmoc_cells = openmoc_geometry.getRootUniverse().getAllCells()\n",
    "\n",
    "# Inject multi-group cross sections into OpenMOC Materials\n",
    "for cell_id, cell in openmoc_cells.items():\n",
    "\n",
    "    # Ignore the root cell\n",
    "    if cell.getName() == 'root cell':\n",
    "        continue\n",
    "\n",
    "    # Get a reference to the Material filling this Cell\n",
    "    openmoc_material = cell.getFillMaterial()\n",
    "\n",
    "    # Set the number of energy groups for the Material\n",
    "    openmoc_material.setNumEnergyGroups(fine_groups.num_groups)\n",
    "\n",
    "    # Extract the appropriate cross section objects for this cell\n",
    "    transport = xs_library[cell_id]['transport']\n",
    "    nufission = xs_library[cell_id]['nu-fission']\n",
    "    nuscatter = xs_library[cell_id]['nu-scatter']\n",
    "    chi = xs_library[cell_id]['chi']\n",
    "\n",
    "    # Inject NumPy arrays of cross section data into the Material\n",
    "    # NOTE: Sum across nuclides to get macro cross sections needed by OpenMOC\n",
    "    openmoc_material.setSigmaT(transport.get_xs(nuclides='sum').flatten())\n",
    "    openmoc_material.setNuSigmaF(nufission.get_xs(nuclides='sum').flatten())\n",
    "    openmoc_material.setSigmaS(nuscatter.get_xs(nuclides='sum').flatten())\n",
    "    openmoc_material.setChi(chi.get_xs(nuclides='sum').flatten())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We are now ready to run OpenMOC to verify our cross-sections from OpenMC."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[  NORMAL ]  Initializing a default angular quadrature...\n",
      "[  NORMAL ]  Initializing 2D tracks...\n",
      "[  NORMAL ]  Initializing 2D tracks reflections...\n",
      "[  NORMAL ]  Initializing 2D tracks array...\n",
      "[  NORMAL ]  Ray tracing for 2D track segmentation...\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 0.09 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 10.02 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 19.94 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 29.87 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 39.80 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 49.72 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 59.65 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 69.58 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 79.50 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 89.43 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 100.00 %\n",
      "[  NORMAL ]  Initializing FSR lookup vectors\n",
      "[  NORMAL ]  Total number of FSRs 3\n",
      "[  NORMAL ]  Initializing MOC eigenvalue solver...\n",
      "[  NORMAL ]  Initializing solver arrays...\n",
      "[  NORMAL ]  Centering segments around FSR centroid...\n",
      "[  NORMAL ]  Max boundary angular flux storage per domain =   0.42 MB\n",
      "[  NORMAL ]  Max scalar flux storage per domain =   0.00 MB\n",
      "[  NORMAL ]  Max source storage per domain =   0.00 MB\n",
      "[  NORMAL ]  Number of azimuthal angles = 128\n",
      "[  NORMAL ]  Azimuthal ray spacing = 0.100000\n",
      "[  NORMAL ]  Number of polar angles = 6\n",
      "[  NORMAL ]  Source type = Flat\n",
      "[  NORMAL ]  MOC transport undamped\n",
      "[  NORMAL ]  CMFD acceleration: OFF\n",
      "[  NORMAL ]  Using 1 threads\n",
      "[  NORMAL ]  Computing the eigenvalue...\n",
      "[  NORMAL ]  Iteration 0:  k_eff = 0.422947   res = 4.915E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -57705 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 1:  k_eff = 0.475805   res = 2.518E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  5285 D.R. = 5.1231\n",
      "[  NORMAL ]  Iteration 2:  k_eff = 0.491328   res = 4.930E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1552 D.R. = 1.9580\n",
      "[  NORMAL ]  Iteration 3:  k_eff = 0.487304   res = 8.137E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -402 D.R. = 1.6503\n",
      "[  NORMAL ]  Iteration 4:  k_eff = 0.483787   res = 2.904E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -351 D.R. = 0.3569\n",
      "[  NORMAL ]  Iteration 5:  k_eff = 0.477137   res = 2.435E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -664 D.R. = 0.8385\n",
      "[  NORMAL ]  Iteration 6:  k_eff = 0.468798   res = 3.811E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -833 D.R. = 1.5652\n",
      "[  NORMAL ]  Iteration 7:  k_eff = 0.460184   res = 1.231E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -861 D.R. = 3.2302\n",
      "[  NORMAL ]  Iteration 8:  k_eff = 0.450461   res = 4.870E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -972 D.R. = 0.3956\n",
      "[  NORMAL ]  Iteration 9:  k_eff = 0.441253   res = 2.027E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -920 D.R. = 0.4161\n",
      "[  NORMAL ]  Iteration 10:  k_eff = 0.431872   res = 4.507E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -938 D.R. = 2.2239\n",
      "[  NORMAL ]  Iteration 11:  k_eff = 0.422820   res = 3.781E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -905 D.R. = 0.8389\n",
      "[  NORMAL ]  Iteration 12:  k_eff = 0.414380   res = 2.450E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -843 D.R. = 0.6480\n",
      "[  NORMAL ]  Iteration 13:  k_eff = 0.406607   res = 3.085E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -777 D.R. = 1.2593\n",
      "[  NORMAL ]  Iteration 14:  k_eff = 0.399282   res = 4.537E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -732 D.R. = 1.4706\n",
      "[  NORMAL ]  Iteration 15:  k_eff = 0.392976   res = 1.815E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -630 D.R. = 0.0400\n",
      "[  NORMAL ]  Iteration 16:  k_eff = 0.387339   res = 0.000E+00  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -563 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 17:  k_eff = 0.382584   res = 4.235E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -475 D.R. = inf\n",
      "[  NORMAL ]  Iteration 18:  k_eff = 0.378660   res = 1.149E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -392 D.R. = 2.7143\n",
      "[  NORMAL ]  Iteration 19:  k_eff = 0.375564   res = 0.000E+00  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -309 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 20:  k_eff = 0.373412   res = 1.875E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -215 D.R. = inf\n",
      "[  NORMAL ]  Iteration 21:  k_eff = 0.372281   res = 1.875E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -113 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 22:  k_eff = 0.371898   res = 1.633E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -38 D.R. = 0.8710\n",
      "[  NORMAL ]  Iteration 23:  k_eff = 0.372505   res = 1.028E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  60 D.R. = 6.2963\n",
      "[  NORMAL ]  Iteration 24:  k_eff = 0.373979   res = 6.170E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  147 D.R. = 0.6000\n",
      "[  NORMAL ]  Iteration 25:  k_eff = 0.376305   res = 2.783E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  232 D.R. = 0.4510\n",
      "[  NORMAL ]  Iteration 26:  k_eff = 0.379483   res = 2.178E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  317 D.R. = 0.7826\n",
      "[  NORMAL ]  Iteration 27:  k_eff = 0.383500   res = 2.662E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  401 D.R. = 1.2222\n",
      "[  NORMAL ]  Iteration 28:  k_eff = 0.388295   res = 2.541E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  479 D.R. = 0.9545\n",
      "[  NORMAL ]  Iteration 29:  k_eff = 0.393849   res = 3.630E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  555 D.R. = 1.4286\n",
      "[  NORMAL ]  Iteration 30:  k_eff = 0.400141   res = 3.630E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  629 D.R. = 0.1000\n",
      "[  NORMAL ]  Iteration 31:  k_eff = 0.407137   res = 4.719E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  699 D.R. = 13.0000\n",
      "[  NORMAL ]  Iteration 32:  k_eff = 0.414782   res = 5.686E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  764 D.R. = 1.2051\n",
      "[  NORMAL ]  Iteration 33:  k_eff = 0.423064   res = 3.267E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  828 D.R. = 0.5745\n",
      "[  NORMAL ]  Iteration 34:  k_eff = 0.431936   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  887 D.R. = 0.5926\n",
      "[  NORMAL ]  Iteration 35:  k_eff = 0.441350   res = 1.452E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  941 D.R. = 0.7500\n",
      "[  NORMAL ]  Iteration 36:  k_eff = 0.451302   res = 1.210E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  995 D.R. = 0.0833\n",
      "[  NORMAL ]  Iteration 37:  k_eff = 0.461718   res = 1.694E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1041 D.R. = 14.0000\n",
      "[  NORMAL ]  Iteration 38:  k_eff = 0.472586   res = 1.089E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1086 D.R. = 0.6429\n",
      "[  NORMAL ]  Iteration 39:  k_eff = 0.483854   res = 3.751E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1126 D.R. = 3.4444\n",
      "[  NORMAL ]  Iteration 40:  k_eff = 0.495492   res = 7.259E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1163 D.R. = 0.1935\n",
      "[  NORMAL ]  Iteration 41:  k_eff = 0.507463   res = 3.025E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1197 D.R. = 4.1667\n",
      "[  NORMAL ]  Iteration 42:  k_eff = 0.519734   res = 3.751E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1227 D.R. = 1.2400\n",
      "[  NORMAL ]  Iteration 43:  k_eff = 0.532267   res = 3.025E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1253 D.R. = 0.8065\n",
      "[  NORMAL ]  Iteration 44:  k_eff = 0.545033   res = 5.445E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1276 D.R. = 1.8000\n",
      "[  NORMAL ]  Iteration 45:  k_eff = 0.557997   res = 2.904E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1296 D.R. = 0.5333\n",
      "[  NORMAL ]  Iteration 46:  k_eff = 0.571128   res = 2.178E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1313 D.R. = 0.7500\n",
      "[  NORMAL ]  Iteration 47:  k_eff = 0.584397   res = 7.259E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1326 D.R. = 0.3333\n",
      "[  NORMAL ]  Iteration 48:  k_eff = 0.597776   res = 9.074E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1337 D.R. = 1.2500\n",
      "[  NORMAL ]  Iteration 49:  k_eff = 0.611237   res = 4.235E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1346 D.R. = 0.4667\n",
      "[  NORMAL ]  Iteration 50:  k_eff = 0.624752   res = 5.263E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1351 D.R. = 12.4286\n",
      "[  NORMAL ]  Iteration 51:  k_eff = 0.638299   res = 4.779E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1354 D.R. = 0.9080\n",
      "[  NORMAL ]  Iteration 52:  k_eff = 0.651853   res = 3.993E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1355 D.R. = 0.8354\n",
      "[  NORMAL ]  Iteration 53:  k_eff = 0.665390   res = 4.053E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1353 D.R. = 1.0152\n",
      "[  NORMAL ]  Iteration 54:  k_eff = 0.678893   res = 3.025E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1350 D.R. = 0.0746\n",
      "[  NORMAL ]  Iteration 55:  k_eff = 0.692338   res = 2.904E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1344 D.R. = 9.6000\n",
      "[  NORMAL ]  Iteration 56:  k_eff = 0.705709   res = 3.448E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1337 D.R. = 1.1875\n",
      "[  NORMAL ]  Iteration 57:  k_eff = 0.718987   res = 4.840E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1327 D.R. = 0.1404\n",
      "[  NORMAL ]  Iteration 58:  k_eff = 0.732158   res = 3.327E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1317 D.R. = 6.8750\n",
      "[  NORMAL ]  Iteration 59:  k_eff = 0.745205   res = 1.754E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1304 D.R. = 0.5273\n",
      "[  NORMAL ]  Iteration 60:  k_eff = 0.758114   res = 1.512E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1290 D.R. = 0.8621\n",
      "[  NORMAL ]  Iteration 61:  k_eff = 0.770875   res = 3.267E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1276 D.R. = 2.1600\n",
      "[  NORMAL ]  Iteration 62:  k_eff = 0.783473   res = 2.541E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1259 D.R. = 0.7778\n",
      "[  NORMAL ]  Iteration 63:  k_eff = 0.795900   res = 1.815E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1242 D.R. = 0.7143\n",
      "[  NORMAL ]  Iteration 64:  k_eff = 0.808146   res = 1.089E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1224 D.R. = 0.6000\n",
      "[  NORMAL ]  Iteration 65:  k_eff = 0.820201   res = 1.391E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1205 D.R. = 1.2778\n",
      "[  NORMAL ]  Iteration 66:  k_eff = 0.832058   res = 2.057E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1185 D.R. = 1.4783\n",
      "[  NORMAL ]  Iteration 67:  k_eff = 0.843711   res = 2.601E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1165 D.R. = 1.2647\n",
      "[  NORMAL ]  Iteration 68:  k_eff = 0.855153   res = 2.783E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1144 D.R. = 1.0698\n",
      "[  NORMAL ]  Iteration 69:  k_eff = 0.866379   res = 3.327E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1122 D.R. = 1.1957\n",
      "[  NORMAL ]  Iteration 70:  k_eff = 0.877384   res = 6.049E-10  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1100 D.R. = 0.0182\n",
      "[  NORMAL ]  Iteration 71:  k_eff = 0.888166   res = 4.477E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1078 D.R. = 74.0000\n",
      "[  NORMAL ]  Iteration 72:  k_eff = 0.898721   res = 1.452E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1055 D.R. = 0.3243\n",
      "[  NORMAL ]  Iteration 73:  k_eff = 0.909046   res = 4.174E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1032 D.R. = 2.8750\n",
      "[  NORMAL ]  Iteration 74:  k_eff = 0.919140   res = 1.754E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1009 D.R. = 0.4203\n",
      "[  NORMAL ]  Iteration 75:  k_eff = 0.929002   res = 2.117E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  986 D.R. = 1.2069\n",
      "[  NORMAL ]  Iteration 76:  k_eff = 0.938631   res = 2.783E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  962 D.R. = 1.3143\n",
      "[  NORMAL ]  Iteration 77:  k_eff = 0.948026   res = 5.263E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  939 D.R. = 1.8913\n",
      "[  NORMAL ]  Iteration 78:  k_eff = 0.957187   res = 2.722E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  916 D.R. = 0.5172\n",
      "[  NORMAL ]  Iteration 79:  k_eff = 0.966117   res = 7.864E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  892 D.R. = 0.2889\n",
      "[  NORMAL ]  Iteration 80:  k_eff = 0.974815   res = 2.904E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  869 D.R. = 3.6923\n",
      "[  NORMAL ]  Iteration 81:  k_eff = 0.983284   res = 3.085E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  846 D.R. = 1.0625\n",
      "[  NORMAL ]  Iteration 82:  k_eff = 0.991523   res = 2.420E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  823 D.R. = 0.0784\n",
      "[  NORMAL ]  Iteration 83:  k_eff = 0.999538   res = 6.049E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  801 D.R. = 2.5000\n",
      "[  NORMAL ]  Iteration 84:  k_eff = 1.007328   res = 2.904E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  779 D.R. = 4.8000\n",
      "[  NORMAL ]  Iteration 85:  k_eff = 1.014897   res = 1.633E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  756 D.R. = 0.5625\n",
      "[  NORMAL ]  Iteration 86:  k_eff = 1.022248   res = 2.662E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  735 D.R. = 1.6296\n",
      "[  NORMAL ]  Iteration 87:  k_eff = 1.029384   res = 1.573E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  713 D.R. = 0.5909\n",
      "[  NORMAL ]  Iteration 88:  k_eff = 1.036308   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  692 D.R. = 1.2308\n",
      "[  NORMAL ]  Iteration 89:  k_eff = 1.043024   res = 1.512E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  671 D.R. = 0.7812\n",
      "[  NORMAL ]  Iteration 90:  k_eff = 1.049535   res = 2.964E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  651 D.R. = 1.9600\n",
      "[  NORMAL ]  Iteration 91:  k_eff = 1.055844   res = 5.384E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  630 D.R. = 1.8163\n",
      "[  NORMAL ]  Iteration 92:  k_eff = 1.061956   res = 3.267E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  611 D.R. = 0.6067\n",
      "[  NORMAL ]  Iteration 93:  k_eff = 1.067875   res = 2.299E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  591 D.R. = 0.7037\n",
      "[  NORMAL ]  Iteration 94:  k_eff = 1.073603   res = 3.630E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  572 D.R. = 0.1579\n",
      "[  NORMAL ]  Iteration 95:  k_eff = 1.079147   res = 4.114E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  554 D.R. = 11.3333\n",
      "[  NORMAL ]  Iteration 96:  k_eff = 1.084509   res = 7.804E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  536 D.R. = 1.8971\n",
      "[  NORMAL ]  Iteration 97:  k_eff = 1.089693   res = 3.146E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  518 D.R. = 0.4031\n",
      "[  NORMAL ]  Iteration 98:  k_eff = 1.094704   res = 6.836E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  501 D.R. = 2.1731\n",
      "[  NORMAL ]  Iteration 99:  k_eff = 1.099547   res = 4.598E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  484 D.R. = 0.6726\n",
      "[  NORMAL ]  Iteration 100:  k_eff = 1.104225   res = 2.480E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 467 D.R. = 0.5395\n",
      "[  NORMAL ]  Iteration 101:  k_eff = 1.108742   res = 5.505E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 451 D.R. = 2.2195\n",
      "[  NORMAL ]  Iteration 102:  k_eff = 1.113102   res = 5.082E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 436 D.R. = 0.9231\n",
      "[  NORMAL ]  Iteration 103:  k_eff = 1.117310   res = 4.779E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 420 D.R. = 0.9405\n",
      "[  NORMAL ]  Iteration 104:  k_eff = 1.121371   res = 1.089E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 406 D.R. = 0.2278\n",
      "[  NORMAL ]  Iteration 105:  k_eff = 1.125287   res = 3.569E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 391 D.R. = 3.2778\n",
      "[  NORMAL ]  Iteration 106:  k_eff = 1.129062   res = 3.630E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 377 D.R. = 0.1017\n",
      "[  NORMAL ]  Iteration 107:  k_eff = 1.132703   res = 6.654E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 364 D.R. = 1.8333\n",
      "[  NORMAL ]  Iteration 108:  k_eff = 1.136210   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 350 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 109:  k_eff = 1.139590   res = 1.149E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 337 D.R. = inf\n",
      "[  NORMAL ]  Iteration 110:  k_eff = 1.142845   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 325 D.R. = 2.5263\n",
      "[  NORMAL ]  Iteration 111:  k_eff = 1.145979   res = 2.178E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 313 D.R. = 0.7500\n",
      "[  NORMAL ]  Iteration 112:  k_eff = 1.148997   res = 1.210E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 301 D.R. = 0.5556\n",
      "[  NORMAL ]  Iteration 113:  k_eff = 1.151901   res = 4.840E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 290 D.R. = 0.4000\n",
      "[  NORMAL ]  Iteration 114:  k_eff = 1.154696   res = 2.117E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 279 D.R. = 4.3750\n",
      "[  NORMAL ]  Iteration 115:  k_eff = 1.157385   res = 1.210E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 268 D.R. = 0.5714\n",
      "[  NORMAL ]  Iteration 116:  k_eff = 1.159971   res = 7.864E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 258 D.R. = 0.6500\n",
      "[  NORMAL ]  Iteration 117:  k_eff = 1.162458   res = 4.658E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 248 D.R. = 5.9231\n",
      "[  NORMAL ]  Iteration 118:  k_eff = 1.164849   res = 4.235E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 239 D.R. = 0.9091\n",
      "[  NORMAL ]  Iteration 119:  k_eff = 1.167147   res = 3.509E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 229 D.R. = 0.8286\n",
      "[  NORMAL ]  Iteration 120:  k_eff = 1.169355   res = 9.074E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 220 D.R. = 0.2586\n",
      "[  NORMAL ]  Iteration 121:  k_eff = 1.171476   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 212 D.R. = 2.1333\n",
      "[  NORMAL ]  Iteration 122:  k_eff = 1.173514   res = 3.206E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 203 D.R. = 1.6563\n",
      "[  NORMAL ]  Iteration 123:  k_eff = 1.175471   res = 1.996E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 195 D.R. = 0.6226\n",
      "[  NORMAL ]  Iteration 124:  k_eff = 1.177351   res = 1.149E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 187 D.R. = 0.5758\n",
      "[  NORMAL ]  Iteration 125:  k_eff = 1.179155   res = 1.331E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 180 D.R. = 1.1579\n",
      "[  NORMAL ]  Iteration 126:  k_eff = 1.180887   res = 3.146E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 173 D.R. = 2.3636\n",
      "[  NORMAL ]  Iteration 127:  k_eff = 1.182549   res = 5.021E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 166 D.R. = 1.5962\n",
      "[  NORMAL ]  Iteration 128:  k_eff = 1.184143   res = 3.932E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 159 D.R. = 0.7831\n",
      "[  NORMAL ]  Iteration 129:  k_eff = 1.185673   res = 4.840E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 152 D.R. = 0.1231\n",
      "[  NORMAL ]  Iteration 130:  k_eff = 1.187140   res = 3.025E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 146 D.R. = 0.6250\n",
      "[  NORMAL ]  Iteration 131:  k_eff = 1.188547   res = 6.231E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 140 D.R. = 20.6000\n",
      "[  NORMAL ]  Iteration 132:  k_eff = 1.189896   res = 4.840E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 134 D.R. = 0.0777\n",
      "[  NORMAL ]  Iteration 133:  k_eff = 1.191189   res = 5.082E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 129 D.R. = 10.5000\n",
      "[  NORMAL ]  Iteration 134:  k_eff = 1.192429   res = 1.815E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 123 D.R. = 0.3571\n",
      "[  NORMAL ]  Iteration 135:  k_eff = 1.193617   res = 6.231E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 118 D.R. = 3.4333\n",
      "[  NORMAL ]  Iteration 136:  k_eff = 1.194755   res = 5.686E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 113 D.R. = 0.9126\n",
      "[  NORMAL ]  Iteration 137:  k_eff = 1.195846   res = 5.021E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 109 D.R. = 0.8830\n",
      "[  NORMAL ]  Iteration 138:  k_eff = 1.196890   res = 6.473E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 104 D.R. = 1.2892\n",
      "[  NORMAL ]  Iteration 139:  k_eff = 1.197891   res = 8.409E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 100 D.R. = 1.2991\n",
      "[  NORMAL ]  Iteration 140:  k_eff = 1.198849   res = 8.106E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 95 D.R. = 0.9640\n",
      "[  NORMAL ]  Iteration 141:  k_eff = 1.199767   res = 5.445E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 91 D.R. = 0.6716\n",
      "[  NORMAL ]  Iteration 142:  k_eff = 1.200645   res = 5.868E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 87 D.R. = 1.0778\n",
      "[  NORMAL ]  Iteration 143:  k_eff = 1.201486   res = 4.416E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 84 D.R. = 0.7526\n",
      "[  NORMAL ]  Iteration 144:  k_eff = 1.202291   res = 1.149E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 80 D.R. = 0.2603\n",
      "[  NORMAL ]  Iteration 145:  k_eff = 1.203060   res = 1.996E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 76 D.R. = 1.7368\n",
      "[  NORMAL ]  Iteration 146:  k_eff = 1.203797   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 73 D.R. = 3.8788\n",
      "[  NORMAL ]  Iteration 147:  k_eff = 1.204502   res = 4.598E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 70 D.R. = 0.5937\n",
      "[  NORMAL ]  Iteration 148:  k_eff = 1.205177   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 67 D.R. = 0.4211\n",
      "[  NORMAL ]  Iteration 149:  k_eff = 1.205822   res = 5.928E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 64 D.R. = 3.0625\n",
      "[  NORMAL ]  Iteration 150:  k_eff = 1.206439   res = 4.053E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 61 D.R. = 0.6837\n",
      "[  NORMAL ]  Iteration 151:  k_eff = 1.207028   res = 6.049E-10  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 58 D.R. = 0.0149\n",
      "[  NORMAL ]  Iteration 152:  k_eff = 1.207592   res = 2.843E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 56 D.R. = 47.0000\n",
      "[  NORMAL ]  Iteration 153:  k_eff = 1.208132   res = 3.085E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 53 D.R. = 1.0851\n",
      "[  NORMAL ]  Iteration 154:  k_eff = 1.208647   res = 4.235E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 51 D.R. = 0.1373\n",
      "[  NORMAL ]  Iteration 155:  k_eff = 1.209139   res = 1.512E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 49 D.R. = 3.5714\n",
      "[  NORMAL ]  Iteration 156:  k_eff = 1.209611   res = 2.117E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 47 D.R. = 1.4000\n",
      "[  NORMAL ]  Iteration 157:  k_eff = 1.210060   res = 5.989E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 44 D.R. = 2.8286\n",
      "[  NORMAL ]  Iteration 158:  k_eff = 1.210490   res = 4.053E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 42 D.R. = 0.6768\n",
      "[  NORMAL ]  Iteration 159:  k_eff = 1.210901   res = 1.089E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 41 D.R. = 0.2687\n",
      "[  NORMAL ]  Iteration 160:  k_eff = 1.211293   res = 6.049E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 39 D.R. = 0.5556\n",
      "[  NORMAL ]  Iteration 161:  k_eff = 1.211668   res = 2.541E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 37 D.R. = 4.2000\n",
      "[  NORMAL ]  Iteration 162:  k_eff = 1.212026   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 35 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 163:  k_eff = 1.212367   res = 7.864E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 34 D.R. = inf\n",
      "[  NORMAL ]  Iteration 164:  k_eff = 1.212693   res = 2.783E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 32 D.R. = 3.5385\n",
      "[  NORMAL ]  Iteration 165:  k_eff = 1.213005   res = 1.815E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 31 D.R. = 0.0652\n",
      "[  NORMAL ]  Iteration 166:  k_eff = 1.213302   res = 4.840E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 29 D.R. = 2.6667\n",
      "[  NORMAL ]  Iteration 167:  k_eff = 1.213587   res = 1.149E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 28 D.R. = 2.3750\n",
      "[  NORMAL ]  Iteration 168:  k_eff = 1.213857   res = 7.864E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 27 D.R. = 0.6842\n",
      "[  NORMAL ]  Iteration 169:  k_eff = 1.214116   res = 3.932E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 25 D.R. = 5.0000\n",
      "[  NORMAL ]  Iteration 170:  k_eff = 1.214363   res = 2.299E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 24 D.R. = 0.5846\n",
      "[  NORMAL ]  Iteration 171:  k_eff = 1.214598   res = 5.445E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 23 D.R. = 2.3684\n",
      "[  NORMAL ]  Iteration 172:  k_eff = 1.214823   res = 5.142E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 22 D.R. = 0.9444\n",
      "[  NORMAL ]  Iteration 173:  k_eff = 1.215038   res = 3.146E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 21 D.R. = 0.6118\n",
      "[  NORMAL ]  Iteration 174:  k_eff = 1.215242   res = 1.633E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 20 D.R. = 0.5192\n",
      "[  NORMAL ]  Iteration 175:  k_eff = 1.215438   res = 3.025E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 19 D.R. = 1.8519\n",
      "[  NORMAL ]  Iteration 176:  k_eff = 1.215624   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 18 D.R. = 0.3200\n",
      "[  NORMAL ]  Iteration 177:  k_eff = 1.215801   res = 8.469E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 17 D.R. = 8.7500\n",
      "[  NORMAL ]  Iteration 178:  k_eff = 1.215971   res = 6.352E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 16 D.R. = 0.7500\n",
      "[  NORMAL ]  Iteration 179:  k_eff = 1.216132   res = 5.021E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 16 D.R. = 0.7905\n",
      "[  NORMAL ]  Iteration 180:  k_eff = 1.216286   res = 6.231E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 15 D.R. = 1.2410\n",
      "[  NORMAL ]  Iteration 181:  k_eff = 1.216433   res = 6.594E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 14 D.R. = 1.0583\n",
      "[  NORMAL ]  Iteration 182:  k_eff = 1.216573   res = 3.630E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 14 D.R. = 0.5505\n",
      "[  NORMAL ]  Iteration 183:  k_eff = 1.216706   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 13 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 184:  k_eff = 1.216834   res = 2.541E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 12 D.R. = inf\n",
      "[  NORMAL ]  Iteration 185:  k_eff = 1.216955   res = 2.359E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 12 D.R. = 0.9286\n",
      "[  NORMAL ]  Iteration 186:  k_eff = 1.217071   res = 2.541E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 11 D.R. = 1.0769\n",
      "[  NORMAL ]  Iteration 187:  k_eff = 1.217181   res = 4.416E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 10 D.R. = 1.7381\n",
      "[  NORMAL ]  Iteration 188:  k_eff = 1.217285   res = 5.505E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 10 D.R. = 1.2466\n",
      "[  NORMAL ]  Iteration 189:  k_eff = 1.217385   res = 2.541E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 10 D.R. = 0.4615\n",
      "[  NORMAL ]  Iteration 190:  k_eff = 1.217481   res = 2.117E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 9 D.R. = 0.8333\n",
      "[  NORMAL ]  Iteration 191:  k_eff = 1.217572   res = 3.630E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 9 D.R. = 0.1714\n",
      "[  NORMAL ]  Iteration 192:  k_eff = 1.217658   res = 4.900E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 8 D.R. = 13.5000\n",
      "[  NORMAL ]  Iteration 193:  k_eff = 1.217741   res = 2.420E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 8 D.R. = 0.4938\n",
      "[  NORMAL ]  Iteration 194:  k_eff = 1.217819   res = 1.815E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 7 D.R. = 0.0750\n",
      "[  NORMAL ]  Iteration 195:  k_eff = 1.217894   res = 9.074E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 7 D.R. = 5.0000\n",
      "[  NORMAL ]  Iteration 196:  k_eff = 1.217965   res = 1.815E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 7 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 197:  k_eff = 1.218033   res = 8.469E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 6 D.R. = 0.4667\n",
      "[  NORMAL ]  Iteration 198:  k_eff = 1.218098   res = 1.210E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 6 D.R. = 0.1429\n",
      "[  NORMAL ]  Iteration 199:  k_eff = 1.218159   res = 3.025E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 6 D.R. = 2.5000\n",
      "[  NORMAL ]  Iteration 200:  k_eff = 1.218218   res = 6.654E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 5 D.R. = 2.2000\n",
      "[  NORMAL ]  Iteration 201:  k_eff = 1.218274   res = 2.662E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 5 D.R. = 4.0000\n",
      "[  NORMAL ]  Iteration 202:  k_eff = 1.218327   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 5 D.R. = 0.3636\n",
      "[  NORMAL ]  Iteration 203:  k_eff = 1.218377   res = 1.331E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 5 D.R. = 1.3750\n",
      "[  NORMAL ]  Iteration 204:  k_eff = 1.218425   res = 2.541E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = 1.9091\n",
      "[  NORMAL ]  Iteration 205:  k_eff = 1.218471   res = 1.512E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = 0.5952\n",
      "[  NORMAL ]  Iteration 206:  k_eff = 1.218515   res = 3.630E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = 0.2400\n",
      "[  NORMAL ]  Iteration 207:  k_eff = 1.218556   res = 3.085E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = 8.5000\n",
      "[  NORMAL ]  Iteration 208:  k_eff = 1.218596   res = 5.445E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 1.7647\n",
      "[  NORMAL ]  Iteration 209:  k_eff = 1.218634   res = 6.049E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 1.1111\n",
      "[  NORMAL ]  Iteration 210:  k_eff = 1.218669   res = 7.864E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 0.1300\n",
      "[  NORMAL ]  Iteration 211:  k_eff = 1.218704   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 4.9231\n",
      "[  NORMAL ]  Iteration 212:  k_eff = 1.218736   res = 6.654E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 0.1719\n",
      "[  NORMAL ]  Iteration 213:  k_eff = 1.218766   res = 2.722E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 4.0909\n",
      "[  NORMAL ]  Iteration 214:  k_eff = 1.218796   res = 1.210E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.0444\n",
      "[  NORMAL ]  Iteration 215:  k_eff = 1.218824   res = 8.046E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 66.5000\n",
      "[  NORMAL ]  Iteration 216:  k_eff = 1.218850   res = 9.982E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 1.2406\n",
      "[  NORMAL ]  Iteration 217:  k_eff = 1.218875   res = 3.630E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.3636\n",
      "[  NORMAL ]  Iteration 218:  k_eff = 1.218900   res = 1.633E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.4500\n",
      "[  NORMAL ]  Iteration 219:  k_eff = 1.218923   res = 2.420E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.1481\n",
      "[  NORMAL ]  Iteration 220:  k_eff = 1.218944   res = 3.085E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 12.7500\n",
      "[  NORMAL ]  Iteration 221:  k_eff = 1.218965   res = 1.452E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.4706\n",
      "[  NORMAL ]  Iteration 222:  k_eff = 1.218985   res = 2.057E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 1.4167\n",
      "[  NORMAL ]  Iteration 223:  k_eff = 1.219003   res = 1.089E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.5294\n",
      "[  NORMAL ]  Iteration 224:  k_eff = 1.219021   res = 5.686E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 5.2222\n",
      "[  NORMAL ]  Iteration 225:  k_eff = 1.219038   res = 1.089E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.1915\n",
      "[  NORMAL ]  Iteration 226:  k_eff = 1.219054   res = 2.117E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 1.9444\n",
      "[  NORMAL ]  Iteration 227:  k_eff = 1.219070   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 1.3714\n",
      "[  NORMAL ]  Iteration 228:  k_eff = 1.219084   res = 3.388E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 1.1667\n",
      "[  NORMAL ]  Iteration 229:  k_eff = 1.219098   res = 2.117E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.6250\n",
      "[  NORMAL ]  Iteration 230:  k_eff = 1.219111   res = 1.210E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.5714\n",
      "[  NORMAL ]  Iteration 231:  k_eff = 1.219123   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 2.4000\n",
      "[  NORMAL ]  Iteration 232:  k_eff = 1.219135   res = 1.633E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.5625\n",
      "[  NORMAL ]  Iteration 233:  k_eff = 1.219146   res = 6.654E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.4074\n",
      "[  NORMAL ]  Iteration 234:  k_eff = 1.219157   res = 1.028E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 1.5455\n",
      "[  NORMAL ]  Iteration 235:  k_eff = 1.219167   res = 6.654E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.6471\n",
      "[  NORMAL ]  Iteration 236:  k_eff = 1.219177   res = 5.445E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 0 D.R. = 8.1818\n"
     ]
    }
   ],
   "source": [
    "# Generate tracks for OpenMOC\n",
    "track_generator = openmoc.TrackGenerator(openmoc_geometry, num_azim=128, azim_spacing=0.1)\n",
    "track_generator.generateTracks()\n",
    "\n",
    "# Run OpenMOC\n",
    "solver = openmoc.CPUSolver(track_generator)\n",
    "solver.computeEigenvalue()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We report the eigenvalues computed by OpenMC and OpenMOC here together to summarize our results."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "openmc keff = 1.223594\n",
      "openmoc keff = 1.219177\n",
      "bias [pcm]: -441.7\n"
     ]
    }
   ],
   "source": [
    "# Print report of keff and bias with OpenMC\n",
    "openmoc_keff = solver.getKeff()\n",
    "openmc_keff = sp.keff.n\n",
    "bias = (openmoc_keff - openmc_keff) * 1e5\n",
    "\n",
    "print('openmc keff = {0:1.6f}'.format(openmc_keff))\n",
    "print('openmoc keff = {0:1.6f}'.format(openmoc_keff))\n",
    "print('bias [pcm]: {0:1.1f}'.format(bias))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As a sanity check, let's run a simulation with the coarse 2-group cross sections to ensure that they also produce a reasonable result."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "openmoc_geometry = get_openmoc_geometry(sp.summary.geometry)\n",
    "openmoc_cells = openmoc_geometry.getRootUniverse().getAllCells()\n",
    "\n",
    "# Inject multi-group cross sections into OpenMOC Materials\n",
    "for cell_id, cell in openmoc_cells.items():\n",
    "\n",
    "    # Ignore the root cell\n",
    "    if cell.getName() == 'root cell':\n",
    "        continue\n",
    "\n",
    "    openmoc_material = cell.getFillMaterial()\n",
    "    openmoc_material.setNumEnergyGroups(coarse_groups.num_groups)\n",
    "\n",
    "    # Extract the appropriate cross section objects for this cell\n",
    "    transport = xs_library[cell_id]['transport']\n",
    "    nufission = xs_library[cell_id]['nu-fission']\n",
    "    nuscatter = xs_library[cell_id]['nu-scatter']\n",
    "    chi = xs_library[cell_id]['chi']\n",
    "\n",
    "    # Perform group condensation\n",
    "    transport = transport.get_condensed_xs(coarse_groups)\n",
    "    nufission = nufission.get_condensed_xs(coarse_groups)\n",
    "    nuscatter = nuscatter.get_condensed_xs(coarse_groups)\n",
    "    chi = chi.get_condensed_xs(coarse_groups)\n",
    "\n",
    "    # Inject NumPy arrays of cross section data into the Material\n",
    "    openmoc_material.setSigmaT(transport.get_xs(nuclides='sum').flatten())\n",
    "    openmoc_material.setNuSigmaF(nufission.get_xs(nuclides='sum').flatten())\n",
    "    openmoc_material.setSigmaS(nuscatter.get_xs(nuclides='sum').flatten())\n",
    "    openmoc_material.setChi(chi.get_xs(nuclides='sum').flatten())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[  NORMAL ]  Initializing a default angular quadrature...\n",
      "[  NORMAL ]  Initializing 2D tracks...\n",
      "[  NORMAL ]  Initializing 2D tracks reflections...\n",
      "[  NORMAL ]  Initializing 2D tracks array...\n",
      "[  NORMAL ]  Ray tracing for 2D track segmentation...\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 0.09 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 10.02 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 19.94 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 29.87 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 39.80 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 49.72 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 59.65 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 69.58 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 79.50 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 89.43 %\n",
      "[  NORMAL ]  Progress Segmenting 2D tracks: 100.00 %\n",
      "[  NORMAL ]  Initializing FSR lookup vectors\n",
      "[  NORMAL ]  Total number of FSRs 3\n",
      "[  NORMAL ]  Initializing MOC eigenvalue solver...\n",
      "[  NORMAL ]  Initializing solver arrays...\n",
      "[  NORMAL ]  Centering segments around FSR centroid...\n",
      "[  NORMAL ]  Max boundary angular flux storage per domain =   0.10 MB\n",
      "[  NORMAL ]  Max scalar flux storage per domain =   0.00 MB\n",
      "[  NORMAL ]  Max source storage per domain =   0.00 MB\n",
      "[  NORMAL ]  Number of azimuthal angles = 128\n",
      "[  NORMAL ]  Azimuthal ray spacing = 0.100000\n",
      "[  NORMAL ]  Number of polar angles = 6\n",
      "[  NORMAL ]  Source type = Flat\n",
      "[  NORMAL ]  MOC transport undamped\n",
      "[  NORMAL ]  CMFD acceleration: OFF\n",
      "[  NORMAL ]  Using 1 threads\n",
      "[  NORMAL ]  Computing the eigenvalue...\n",
      "[  NORMAL ]  Iteration 0:  k_eff = 0.366745   res = 6.049E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -63325 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 1:  k_eff = 0.391062   res = 5.324E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  2431 D.R. = 0.8800\n",
      "[  NORMAL ]  Iteration 2:  k_eff = 0.392876   res = 4.840E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  181 D.R. = 0.0909\n",
      "[  NORMAL ]  Iteration 3:  k_eff = 0.380986   res = 5.324E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -1189 D.R. = 11.0000\n",
      "[  NORMAL ]  Iteration 4:  k_eff = 0.374905   res = 2.420E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -608 D.R. = 0.4545\n",
      "[  NORMAL ]  Iteration 5:  k_eff = 0.369482   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -542 D.R. = 0.8000\n",
      "[  NORMAL ]  Iteration 6:  k_eff = 0.365432   res = 4.840E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -404 D.R. = 2.5000\n",
      "[  NORMAL ]  Iteration 7:  k_eff = 0.362942   res = 6.775E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -248 D.R. = 1.4000\n",
      "[  NORMAL ]  Iteration 8:  k_eff = 0.361361   res = 5.807E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -158 D.R. = 0.8571\n",
      "[  NORMAL ]  Iteration 9:  k_eff = 0.361166   res = 6.775E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  -19 D.R. = 1.1667\n",
      "[  NORMAL ]  Iteration 10:  k_eff = 0.361886   res = 1.549E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  72 D.R. = 2.2857\n",
      "[  NORMAL ]  Iteration 11:  k_eff = 0.363597   res = 9.679E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  171 D.R. = 0.0625\n",
      "[  NORMAL ]  Iteration 12:  k_eff = 0.366212   res = 9.679E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  261 D.R. = 10.0000\n",
      "[  NORMAL ]  Iteration 13:  k_eff = 0.369672   res = 9.679E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  346 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 14:  k_eff = 0.373850   res = 1.355E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  417 D.R. = 1.4000\n",
      "[  NORMAL ]  Iteration 15:  k_eff = 0.378775   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  492 D.R. = 0.5714\n",
      "[  NORMAL ]  Iteration 16:  k_eff = 0.384321   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  554 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 17:  k_eff = 0.390468   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  614 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 18:  k_eff = 0.397157   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  668 D.R. = 4.0000\n",
      "[  NORMAL ]  Iteration 19:  k_eff = 0.404338   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  718 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 20:  k_eff = 0.411974   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  763 D.R. = 0.2500\n",
      "[  NORMAL ]  Iteration 21:  k_eff = 0.420027   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  805 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 22:  k_eff = 0.428443   res = 0.000E+00  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  841 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 23:  k_eff = 0.437200   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  875 D.R. = inf\n",
      "[  NORMAL ]  Iteration 24:  k_eff = 0.446255   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  905 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 25:  k_eff = 0.455580   res = 5.807E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  932 D.R. = 1.5000\n",
      "[  NORMAL ]  Iteration 26:  k_eff = 0.465142   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  956 D.R. = 0.6667\n",
      "[  NORMAL ]  Iteration 27:  k_eff = 0.474915   res = 5.807E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  977 D.R. = 1.5000\n",
      "[  NORMAL ]  Iteration 28:  k_eff = 0.484870   res = 1.742E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  995 D.R. = 3.0000\n",
      "[  NORMAL ]  Iteration 29:  k_eff = 0.494984   res = 1.549E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1011 D.R. = 0.8889\n",
      "[  NORMAL ]  Iteration 30:  k_eff = 0.505232   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1024 D.R. = 0.2500\n",
      "[  NORMAL ]  Iteration 31:  k_eff = 0.515593   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1036 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 32:  k_eff = 0.526047   res = 9.679E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1045 D.R. = 1.2500\n",
      "[  NORMAL ]  Iteration 33:  k_eff = 0.536574   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1052 D.R. = 0.2000\n",
      "[  NORMAL ]  Iteration 34:  k_eff = 0.547157   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1058 D.R. = 4.0000\n",
      "[  NORMAL ]  Iteration 35:  k_eff = 0.557779   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1062 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 36:  k_eff = 0.568425   res = 5.807E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1064 D.R. = 1.5000\n",
      "[  NORMAL ]  Iteration 37:  k_eff = 0.579080   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1065 D.R. = 0.6667\n",
      "[  NORMAL ]  Iteration 38:  k_eff = 0.589731   res = 1.161E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1065 D.R. = 3.0000\n",
      "[  NORMAL ]  Iteration 39:  k_eff = 0.600366   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1063 D.R. = 0.6667\n",
      "[  NORMAL ]  Iteration 40:  k_eff = 0.610974   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1060 D.R. = 0.2500\n",
      "[  NORMAL ]  Iteration 41:  k_eff = 0.621542   res = 9.679E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1056 D.R. = 5.0000\n",
      "[  NORMAL ]  Iteration 42:  k_eff = 0.632063   res = 5.807E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1052 D.R. = 0.6000\n",
      "[  NORMAL ]  Iteration 43:  k_eff = 0.642526   res = 0.000E+00  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1046 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 44:  k_eff = 0.652925   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1039 D.R. = inf\n",
      "[  NORMAL ]  Iteration 45:  k_eff = 0.663250   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1032 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 46:  k_eff = 0.673494   res = 0.000E+00  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1024 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 47:  k_eff = 0.683653   res = 0.000E+00  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1015 D.R. = -nan\n",
      "[  NORMAL ]  Iteration 48:  k_eff = 0.693719   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  1006 D.R. = inf\n",
      "[  NORMAL ]  Iteration 49:  k_eff = 0.703687   res = 1.355E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  996 D.R. = 1.7500\n",
      "[  NORMAL ]  Iteration 50:  k_eff = 0.713552   res = 9.679E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  986 D.R. = 0.7143\n",
      "[  NORMAL ]  Iteration 51:  k_eff = 0.723311   res = 5.807E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  975 D.R. = 0.6000\n",
      "[  NORMAL ]  Iteration 52:  k_eff = 0.732959   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  964 D.R. = 0.3333\n",
      "[  NORMAL ]  Iteration 53:  k_eff = 0.742493   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  953 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 54:  k_eff = 0.751909   res = 0.000E+00  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  941 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 55:  k_eff = 0.761205   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  929 D.R. = inf\n",
      "[  NORMAL ]  Iteration 56:  k_eff = 0.770378   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  917 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 57:  k_eff = 0.779427   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  904 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 58:  k_eff = 0.788350   res = 5.807E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  892 D.R. = 0.7500\n",
      "[  NORMAL ]  Iteration 59:  k_eff = 0.797144   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  879 D.R. = 0.6667\n",
      "[  NORMAL ]  Iteration 60:  k_eff = 0.805809   res = 1.355E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  866 D.R. = 3.5000\n",
      "[  NORMAL ]  Iteration 61:  k_eff = 0.814344   res = 1.161E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  853 D.R. = 0.8571\n",
      "[  NORMAL ]  Iteration 62:  k_eff = 0.822748   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  840 D.R. = 0.3333\n",
      "[  NORMAL ]  Iteration 63:  k_eff = 0.831020   res = 9.679E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  827 D.R. = 2.5000\n",
      "[  NORMAL ]  Iteration 64:  k_eff = 0.839161   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  814 D.R. = 0.2000\n",
      "[  NORMAL ]  Iteration 65:  k_eff = 0.847169   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  800 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 66:  k_eff = 0.855045   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  787 D.R. = 4.0000\n",
      "[  NORMAL ]  Iteration 67:  k_eff = 0.862789   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  774 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 68:  k_eff = 0.870401   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  761 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 69:  k_eff = 0.877882   res = 9.679E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  748 D.R. = 1.2500\n",
      "[  NORMAL ]  Iteration 70:  k_eff = 0.885232   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  735 D.R. = 0.2000\n",
      "[  NORMAL ]  Iteration 71:  k_eff = 0.892452   res = 5.807E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  722 D.R. = 3.0000\n",
      "[  NORMAL ]  Iteration 72:  k_eff = 0.899542   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  709 D.R. = 1.3333\n",
      "[  NORMAL ]  Iteration 73:  k_eff = 0.906504   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  696 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 74:  k_eff = 0.913339   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  683 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 75:  k_eff = 0.920047   res = 0.000E+00  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  670 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 76:  k_eff = 0.926629   res = 0.000E+00  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  658 D.R. = -nan\n",
      "[  NORMAL ]  Iteration 77:  k_eff = 0.933088   res = 9.679E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  645 D.R. = inf\n",
      "[  NORMAL ]  Iteration 78:  k_eff = 0.939423   res = 5.807E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  633 D.R. = 0.6000\n",
      "[  NORMAL ]  Iteration 79:  k_eff = 0.945637   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  621 D.R. = 1.3333\n",
      "[  NORMAL ]  Iteration 80:  k_eff = 0.951730   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  609 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 81:  k_eff = 0.957705   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  597 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 82:  k_eff = 0.963562   res = 2.904E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  585 D.R. = 0.7500\n",
      "[  NORMAL ]  Iteration 83:  k_eff = 0.969304   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  574 D.R. = 2.6667\n",
      "[  NORMAL ]  Iteration 84:  k_eff = 0.974931   res = 9.679E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  562 D.R. = 0.1250\n",
      "[  NORMAL ]  Iteration 85:  k_eff = 0.980444   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  551 D.R. = 4.0000\n",
      "[  NORMAL ]  Iteration 86:  k_eff = 0.985847   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  540 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 87:  k_eff = 0.991140   res = 9.679E-09  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  529 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 88:  k_eff = 0.996324   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  518 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 89:  k_eff = 1.001402   res = 1.161E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  507 D.R. = 6.0000\n",
      "[  NORMAL ]  Iteration 90:  k_eff = 1.006375   res = 1.258E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  497 D.R. = 1.0833\n",
      "[  NORMAL ]  Iteration 91:  k_eff = 1.011244   res = 0.000E+00  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  486 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 92:  k_eff = 1.016012   res = 0.000E+00  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  476 D.R. = -nan\n",
      "[  NORMAL ]  Iteration 93:  k_eff = 1.020679   res = 1.258E-07  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  466 D.R. = inf\n",
      "[  NORMAL ]  Iteration 94:  k_eff = 1.025249   res = 2.904E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  456 D.R. = 0.2308\n",
      "[  NORMAL ]  Iteration 95:  k_eff = 1.029722   res = 3.872E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  447 D.R. = 1.3333\n",
      "[  NORMAL ]  Iteration 96:  k_eff = 1.034099   res = 1.936E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  437 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 97:  k_eff = 1.038383   res = 5.807E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  428 D.R. = 3.0000\n",
      "[  NORMAL ]  Iteration 98:  k_eff = 1.042575   res = 7.743E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  419 D.R. = 1.3333\n",
      "[  NORMAL ]  Iteration 99:  k_eff = 1.046677   res = 4.840E-08  delta-k (pcm) =\n",
      "[  NORMAL ]  ...  410 D.R. = 0.6250\n",
      "[  NORMAL ]  Iteration 100:  k_eff = 1.050690   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 401 D.R. = 1.6000\n",
      "[  NORMAL ]  Iteration 101:  k_eff = 1.054616   res = 1.065E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 392 D.R. = 1.3750\n",
      "[  NORMAL ]  Iteration 102:  k_eff = 1.058457   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 384 D.R. = 0.4545\n",
      "[  NORMAL ]  Iteration 103:  k_eff = 1.062214   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 375 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 104:  k_eff = 1.065888   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 367 D.R. = inf\n",
      "[  NORMAL ]  Iteration 105:  k_eff = 1.069481   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 359 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 106:  k_eff = 1.072995   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 351 D.R. = 1.2500\n",
      "[  NORMAL ]  Iteration 107:  k_eff = 1.076432   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 343 D.R. = 1.2000\n",
      "[  NORMAL ]  Iteration 108:  k_eff = 1.079792   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 335 D.R. = 0.1667\n",
      "[  NORMAL ]  Iteration 109:  k_eff = 1.083078   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 328 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 110:  k_eff = 1.086291   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 321 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 111:  k_eff = 1.089431   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 314 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 112:  k_eff = 1.092501   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 307 D.R. = inf\n",
      "[  NORMAL ]  Iteration 113:  k_eff = 1.095502   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 300 D.R. = 0.2500\n",
      "[  NORMAL ]  Iteration 114:  k_eff = 1.098435   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 293 D.R. = 5.0000\n",
      "[  NORMAL ]  Iteration 115:  k_eff = 1.101302   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 286 D.R. = 1.6000\n",
      "[  NORMAL ]  Iteration 116:  k_eff = 1.104105   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 280 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 117:  k_eff = 1.106844   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 273 D.R. = 1.2500\n",
      "[  NORMAL ]  Iteration 118:  k_eff = 1.109520   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 267 D.R. = 1.2000\n",
      "[  NORMAL ]  Iteration 119:  k_eff = 1.112135   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 261 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 120:  k_eff = 1.114691   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 255 D.R. = 1.3333\n",
      "[  NORMAL ]  Iteration 121:  k_eff = 1.117188   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 249 D.R. = 0.3750\n",
      "[  NORMAL ]  Iteration 122:  k_eff = 1.119629   res = 9.679E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 244 D.R. = 3.3333\n",
      "[  NORMAL ]  Iteration 123:  k_eff = 1.122013   res = 9.679E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 238 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 124:  k_eff = 1.124342   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 232 D.R. = 0.4000\n",
      "[  NORMAL ]  Iteration 125:  k_eff = 1.126617   res = 1.065E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 227 D.R. = 2.7500\n",
      "[  NORMAL ]  Iteration 126:  k_eff = 1.128841   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 222 D.R. = 0.6364\n",
      "[  NORMAL ]  Iteration 127:  k_eff = 1.131013   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 217 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 128:  k_eff = 1.133134   res = 1.161E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 212 D.R. = 1.7143\n",
      "[  NORMAL ]  Iteration 129:  k_eff = 1.135206   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 207 D.R. = 0.7500\n",
      "[  NORMAL ]  Iteration 130:  k_eff = 1.137231   res = 1.258E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 202 D.R. = 1.4444\n",
      "[  NORMAL ]  Iteration 131:  k_eff = 1.139208   res = 1.065E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 197 D.R. = 0.8462\n",
      "[  NORMAL ]  Iteration 132:  k_eff = 1.141139   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 193 D.R. = 0.5455\n",
      "[  NORMAL ]  Iteration 133:  k_eff = 1.143025   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 188 D.R. = 0.6667\n",
      "[  NORMAL ]  Iteration 134:  k_eff = 1.144867   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 184 D.R. = 2.2500\n",
      "[  NORMAL ]  Iteration 135:  k_eff = 1.146666   res = 1.549E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 179 D.R. = 1.7778\n",
      "[  NORMAL ]  Iteration 136:  k_eff = 1.148423   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 175 D.R. = 0.3750\n",
      "[  NORMAL ]  Iteration 137:  k_eff = 1.150139   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 171 D.R. = 0.6667\n",
      "[  NORMAL ]  Iteration 138:  k_eff = 1.151815   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 167 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 139:  k_eff = 1.153451   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 163 D.R. = 1.5000\n",
      "[  NORMAL ]  Iteration 140:  k_eff = 1.155049   res = 1.355E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 159 D.R. = 4.6667\n",
      "[  NORMAL ]  Iteration 141:  k_eff = 1.156609   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 156 D.R. = 0.0714\n",
      "[  NORMAL ]  Iteration 142:  k_eff = 1.158132   res = 1.452E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 152 D.R. = 15.0000\n",
      "[  NORMAL ]  Iteration 143:  k_eff = 1.159620   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 148 D.R. = 0.2667\n",
      "[  NORMAL ]  Iteration 144:  k_eff = 1.161073   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 145 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 145:  k_eff = 1.162491   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 141 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 146:  k_eff = 1.163875   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 138 D.R. = 5.0000\n",
      "[  NORMAL ]  Iteration 147:  k_eff = 1.165227   res = 9.679E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 135 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 148:  k_eff = 1.166546   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 131 D.R. = 0.9000\n",
      "[  NORMAL ]  Iteration 149:  k_eff = 1.167835   res = 9.679E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 128 D.R. = 1.1111\n",
      "[  NORMAL ]  Iteration 150:  k_eff = 1.169094   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 125 D.R. = 0.6000\n",
      "[  NORMAL ]  Iteration 151:  k_eff = 1.170322   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 122 D.R. = 0.3333\n",
      "[  NORMAL ]  Iteration 152:  k_eff = 1.171521   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 119 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 153:  k_eff = 1.172691   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 117 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 154:  k_eff = 1.173833   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 114 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 155:  k_eff = 1.174948   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 111 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 156:  k_eff = 1.176037   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 108 D.R. = 5.0000\n",
      "[  NORMAL ]  Iteration 157:  k_eff = 1.177100   res = 1.065E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 106 D.R. = 2.2000\n",
      "[  NORMAL ]  Iteration 158:  k_eff = 1.178137   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 103 D.R. = 0.4545\n",
      "[  NORMAL ]  Iteration 159:  k_eff = 1.179149   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 101 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 160:  k_eff = 1.180137   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 98 D.R. = 0.4000\n",
      "[  NORMAL ]  Iteration 161:  k_eff = 1.181102   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 96 D.R. = 3.0000\n",
      "[  NORMAL ]  Iteration 162:  k_eff = 1.182044   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 94 D.R. = 0.3333\n",
      "[  NORMAL ]  Iteration 163:  k_eff = 1.182963   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 91 D.R. = 1.5000\n",
      "[  NORMAL ]  Iteration 164:  k_eff = 1.183859   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 89 D.R. = 2.6667\n",
      "[  NORMAL ]  Iteration 165:  k_eff = 1.184735   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 87 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 166:  k_eff = 1.185589   res = 1.355E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 85 D.R. = inf\n",
      "[  NORMAL ]  Iteration 167:  k_eff = 1.186424   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 83 D.R. = 0.0714\n",
      "[  NORMAL ]  Iteration 168:  k_eff = 1.187237   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 81 D.R. = 5.0000\n",
      "[  NORMAL ]  Iteration 169:  k_eff = 1.188031   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 79 D.R. = 0.2000\n",
      "[  NORMAL ]  Iteration 170:  k_eff = 1.188807   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 77 D.R. = 5.0000\n",
      "[  NORMAL ]  Iteration 171:  k_eff = 1.189563   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 75 D.R. = 0.4000\n",
      "[  NORMAL ]  Iteration 172:  k_eff = 1.190302   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 73 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 173:  k_eff = 1.191023   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 72 D.R. = inf\n",
      "[  NORMAL ]  Iteration 174:  k_eff = 1.191726   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 70 D.R. = 0.8000\n",
      "[  NORMAL ]  Iteration 175:  k_eff = 1.192413   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 68 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 176:  k_eff = 1.193083   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 67 D.R. = 0.7500\n",
      "[  NORMAL ]  Iteration 177:  k_eff = 1.193736   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 65 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 178:  k_eff = 1.194374   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 63 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 179:  k_eff = 1.194996   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 62 D.R. = inf\n",
      "[  NORMAL ]  Iteration 180:  k_eff = 1.195603   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 60 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 181:  k_eff = 1.196197   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 59 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 182:  k_eff = 1.196775   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 57 D.R. = 2.5000\n",
      "[  NORMAL ]  Iteration 183:  k_eff = 1.197340   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 56 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 184:  k_eff = 1.197890   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 55 D.R. = inf\n",
      "[  NORMAL ]  Iteration 185:  k_eff = 1.198428   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 53 D.R. = 2.3333\n",
      "[  NORMAL ]  Iteration 186:  k_eff = 1.198953   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 52 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 187:  k_eff = 1.199465   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 51 D.R. = 1.1429\n",
      "[  NORMAL ]  Iteration 188:  k_eff = 1.199964   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 49 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 189:  k_eff = 1.200452   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 48 D.R. = 0.3750\n",
      "[  NORMAL ]  Iteration 190:  k_eff = 1.200927   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 47 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 191:  k_eff = 1.201391   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 46 D.R. = 1.3333\n",
      "[  NORMAL ]  Iteration 192:  k_eff = 1.201843   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 45 D.R. = 2.2500\n",
      "[  NORMAL ]  Iteration 193:  k_eff = 1.202285   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 44 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 194:  k_eff = 1.202717   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 43 D.R. = 0.4444\n",
      "[  NORMAL ]  Iteration 195:  k_eff = 1.203138   res = 1.065E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 42 D.R. = 2.7500\n",
      "[  NORMAL ]  Iteration 196:  k_eff = 1.203549   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 41 D.R. = 0.7273\n",
      "[  NORMAL ]  Iteration 197:  k_eff = 1.203949   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 40 D.R. = 0.2500\n",
      "[  NORMAL ]  Iteration 198:  k_eff = 1.204339   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 39 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 199:  k_eff = 1.204721   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 38 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 200:  k_eff = 1.205093   res = 1.065E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 37 D.R. = inf\n",
      "[  NORMAL ]  Iteration 201:  k_eff = 1.205456   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 36 D.R. = 0.0909\n",
      "[  NORMAL ]  Iteration 202:  k_eff = 1.205810   res = 1.161E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 35 D.R. = 12.0000\n",
      "[  NORMAL ]  Iteration 203:  k_eff = 1.206156   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 34 D.R. = 0.3333\n",
      "[  NORMAL ]  Iteration 204:  k_eff = 1.206492   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 33 D.R. = 1.5000\n",
      "[  NORMAL ]  Iteration 205:  k_eff = 1.206822   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 32 D.R. = 1.1667\n",
      "[  NORMAL ]  Iteration 206:  k_eff = 1.207143   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 32 D.R. = 0.1429\n",
      "[  NORMAL ]  Iteration 207:  k_eff = 1.207456   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 31 D.R. = 9.0000\n",
      "[  NORMAL ]  Iteration 208:  k_eff = 1.207762   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 30 D.R. = 0.1111\n",
      "[  NORMAL ]  Iteration 209:  k_eff = 1.208060   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 29 D.R. = 7.0000\n",
      "[  NORMAL ]  Iteration 210:  k_eff = 1.208351   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 29 D.R. = 0.5714\n",
      "[  NORMAL ]  Iteration 211:  k_eff = 1.208634   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 28 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 212:  k_eff = 1.208911   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 27 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 213:  k_eff = 1.209182   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 27 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 214:  k_eff = 1.209445   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 26 D.R. = 1.5000\n",
      "[  NORMAL ]  Iteration 215:  k_eff = 1.209701   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 25 D.R. = 1.6667\n",
      "[  NORMAL ]  Iteration 216:  k_eff = 1.209953   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 25 D.R. = 0.6000\n",
      "[  NORMAL ]  Iteration 217:  k_eff = 1.210198   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 24 D.R. = 0.3333\n",
      "[  NORMAL ]  Iteration 218:  k_eff = 1.210436   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 23 D.R. = 6.0000\n",
      "[  NORMAL ]  Iteration 219:  k_eff = 1.210669   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 23 D.R. = 0.8333\n",
      "[  NORMAL ]  Iteration 220:  k_eff = 1.210897   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 22 D.R. = 0.2000\n",
      "[  NORMAL ]  Iteration 221:  k_eff = 1.211119   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 22 D.R. = 4.0000\n",
      "[  NORMAL ]  Iteration 222:  k_eff = 1.211335   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 21 D.R. = 0.7500\n",
      "[  NORMAL ]  Iteration 223:  k_eff = 1.211546   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 21 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 224:  k_eff = 1.211752   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 20 D.R. = 0.3333\n",
      "[  NORMAL ]  Iteration 225:  k_eff = 1.211953   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 20 D.R. = 3.0000\n",
      "[  NORMAL ]  Iteration 226:  k_eff = 1.212149   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 19 D.R. = 3.0000\n",
      "[  NORMAL ]  Iteration 227:  k_eff = 1.212341   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 19 D.R. = 0.4444\n",
      "[  NORMAL ]  Iteration 228:  k_eff = 1.212527   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 18 D.R. = 1.2500\n",
      "[  NORMAL ]  Iteration 229:  k_eff = 1.212709   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 18 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 230:  k_eff = 1.212886   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 17 D.R. = 1.2000\n",
      "[  NORMAL ]  Iteration 231:  k_eff = 1.213060   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 17 D.R. = 0.8333\n",
      "[  NORMAL ]  Iteration 232:  k_eff = 1.213229   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 16 D.R. = 1.4000\n",
      "[  NORMAL ]  Iteration 233:  k_eff = 1.213394   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 16 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 234:  k_eff = 1.213555   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 16 D.R. = 0.8571\n",
      "[  NORMAL ]  Iteration 235:  k_eff = 1.213711   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 15 D.R. = 1.5000\n",
      "[  NORMAL ]  Iteration 236:  k_eff = 1.213864   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 15 D.R. = 0.4444\n",
      "[  NORMAL ]  Iteration 237:  k_eff = 1.214014   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 14 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 238:  k_eff = 1.214160   res = 1.065E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 14 D.R. = 5.5000\n",
      "[  NORMAL ]  Iteration 239:  k_eff = 1.214302   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 14 D.R. = 0.5455\n",
      "[  NORMAL ]  Iteration 240:  k_eff = 1.214441   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 13 D.R. = 0.8333\n",
      "[  NORMAL ]  Iteration 241:  k_eff = 1.214576   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 13 D.R. = 0.4000\n",
      "[  NORMAL ]  Iteration 242:  k_eff = 1.214708   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 13 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 243:  k_eff = 1.214837   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 12 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 244:  k_eff = 1.214962   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 12 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 245:  k_eff = 1.215085   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 12 D.R. = 3.0000\n",
      "[  NORMAL ]  Iteration 246:  k_eff = 1.215204   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 11 D.R. = 1.6667\n",
      "[  NORMAL ]  Iteration 247:  k_eff = 1.215321   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 11 D.R. = 0.8000\n",
      "[  NORMAL ]  Iteration 248:  k_eff = 1.215434   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 11 D.R. = 1.2500\n",
      "[  NORMAL ]  Iteration 249:  k_eff = 1.215545   res = 1.161E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 11 D.R. = 2.4000\n",
      "[  NORMAL ]  Iteration 250:  k_eff = 1.215654   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 10 D.R. = 0.7500\n",
      "[  NORMAL ]  Iteration 251:  k_eff = 1.215759   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 10 D.R. = 0.8889\n",
      "[  NORMAL ]  Iteration 252:  k_eff = 1.215862   res = 1.355E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 10 D.R. = 1.7500\n",
      "[  NORMAL ]  Iteration 253:  k_eff = 1.215963   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 10 D.R. = 0.2143\n",
      "[  NORMAL ]  Iteration 254:  k_eff = 1.216061   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 9 D.R. = 3.0000\n",
      "[  NORMAL ]  Iteration 255:  k_eff = 1.216157   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 9 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 256:  k_eff = 1.216250   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 9 D.R. = 0.3333\n",
      "[  NORMAL ]  Iteration 257:  k_eff = 1.216341   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 9 D.R. = 0.6667\n",
      "[  NORMAL ]  Iteration 258:  k_eff = 1.216430   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 8 D.R. = 1.5000\n",
      "[  NORMAL ]  Iteration 259:  k_eff = 1.216516   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 8 D.R. = 2.6667\n",
      "[  NORMAL ]  Iteration 260:  k_eff = 1.216601   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 8 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 261:  k_eff = 1.216683   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 8 D.R. = inf\n",
      "[  NORMAL ]  Iteration 262:  k_eff = 1.216764   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 8 D.R. = 1.2500\n",
      "[  NORMAL ]  Iteration 263:  k_eff = 1.216842   res = 1.161E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 7 D.R. = 2.4000\n",
      "[  NORMAL ]  Iteration 264:  k_eff = 1.216919   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 7 D.R. = 0.5833\n",
      "[  NORMAL ]  Iteration 265:  k_eff = 1.216993   res = 1.065E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 7 D.R. = 1.5714\n",
      "[  NORMAL ]  Iteration 266:  k_eff = 1.217066   res = 1.355E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 7 D.R. = 1.2727\n",
      "[  NORMAL ]  Iteration 267:  k_eff = 1.217138   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 7 D.R. = 0.3571\n",
      "[  NORMAL ]  Iteration 268:  k_eff = 1.217207   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 6 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 269:  k_eff = 1.217275   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 6 D.R. = 1.2000\n",
      "[  NORMAL ]  Iteration 270:  k_eff = 1.217341   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 6 D.R. = 0.1667\n",
      "[  NORMAL ]  Iteration 271:  k_eff = 1.217405   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 6 D.R. = 5.0000\n",
      "[  NORMAL ]  Iteration 272:  k_eff = 1.217468   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 6 D.R. = 0.6000\n",
      "[  NORMAL ]  Iteration 273:  k_eff = 1.217529   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 6 D.R. = 2.6667\n",
      "[  NORMAL ]  Iteration 274:  k_eff = 1.217589   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 5 D.R. = 0.7500\n",
      "[  NORMAL ]  Iteration 275:  k_eff = 1.217647   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 5 D.R. = 0.1667\n",
      "[  NORMAL ]  Iteration 276:  k_eff = 1.217704   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 5 D.R. = 9.0000\n",
      "[  NORMAL ]  Iteration 277:  k_eff = 1.217760   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 5 D.R. = 0.8889\n",
      "[  NORMAL ]  Iteration 278:  k_eff = 1.217814   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 5 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 279:  k_eff = 1.217866   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 5 D.R. = 2.0000\n",
      "[  NORMAL ]  Iteration 280:  k_eff = 1.217918   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 5 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 281:  k_eff = 1.217968   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 5 D.R. = inf\n",
      "[  NORMAL ]  Iteration 282:  k_eff = 1.218017   res = 2.904E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = 3.0000\n",
      "[  NORMAL ]  Iteration 283:  k_eff = 1.218065   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = 1.3333\n",
      "[  NORMAL ]  Iteration 284:  k_eff = 1.218112   res = 1.452E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = 3.7500\n",
      "[  NORMAL ]  Iteration 285:  k_eff = 1.218158   res = 1.161E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = 0.8000\n",
      "[  NORMAL ]  Iteration 286:  k_eff = 1.218202   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = 0.3333\n",
      "[  NORMAL ]  Iteration 287:  k_eff = 1.218246   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = 1.7500\n",
      "[  NORMAL ]  Iteration 288:  k_eff = 1.218288   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 289:  k_eff = 1.218329   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = inf\n",
      "[  NORMAL ]  Iteration 290:  k_eff = 1.218370   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 4 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 291:  k_eff = 1.218409   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 9.0000\n",
      "[  NORMAL ]  Iteration 292:  k_eff = 1.218448   res = 0.000E+00  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 0.0000\n",
      "[  NORMAL ]  Iteration 293:  k_eff = 1.218485   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = inf\n",
      "[  NORMAL ]  Iteration 294:  k_eff = 1.218522   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 295:  k_eff = 1.218557   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 7.0000\n",
      "[  NORMAL ]  Iteration 296:  k_eff = 1.218592   res = 1.742E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 2.5714\n",
      "[  NORMAL ]  Iteration 297:  k_eff = 1.218626   res = 2.033E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 1.1667\n",
      "[  NORMAL ]  Iteration 298:  k_eff = 1.218658   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 0.3810\n",
      "[  NORMAL ]  Iteration 299:  k_eff = 1.218691   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 0.1250\n",
      "[  NORMAL ]  Iteration 300:  k_eff = 1.218722   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 301:  k_eff = 1.218753   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 3 D.R. = 8.0000\n",
      "[  NORMAL ]  Iteration 302:  k_eff = 1.218783   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.1250\n",
      "[  NORMAL ]  Iteration 303:  k_eff = 1.218812   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 7.0000\n",
      "[  NORMAL ]  Iteration 304:  k_eff = 1.218840   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.7143\n",
      "[  NORMAL ]  Iteration 305:  k_eff = 1.218868   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.8000\n",
      "[  NORMAL ]  Iteration 306:  k_eff = 1.218895   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 1.7500\n",
      "[  NORMAL ]  Iteration 307:  k_eff = 1.218922   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.5714\n",
      "[  NORMAL ]  Iteration 308:  k_eff = 1.218948   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.5000\n",
      "[  NORMAL ]  Iteration 309:  k_eff = 1.218973   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 4.0000\n",
      "[  NORMAL ]  Iteration 310:  k_eff = 1.218997   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.7500\n",
      "[  NORMAL ]  Iteration 311:  k_eff = 1.219021   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.6667\n",
      "[  NORMAL ]  Iteration 312:  k_eff = 1.219045   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.2500\n",
      "[  NORMAL ]  Iteration 313:  k_eff = 1.219067   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 314:  k_eff = 1.219090   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 5.0000\n",
      "[  NORMAL ]  Iteration 315:  k_eff = 1.219111   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 1.4000\n",
      "[  NORMAL ]  Iteration 316:  k_eff = 1.219133   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 0.1429\n",
      "[  NORMAL ]  Iteration 317:  k_eff = 1.219153   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 318:  k_eff = 1.219173   res = 7.743E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 2 D.R. = 8.0000\n",
      "[  NORMAL ]  Iteration 319:  k_eff = 1.219193   res = 9.679E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 1.2500\n",
      "[  NORMAL ]  Iteration 320:  k_eff = 1.219212   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.1000\n",
      "[  NORMAL ]  Iteration 321:  k_eff = 1.219231   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 7.0000\n",
      "[  NORMAL ]  Iteration 322:  k_eff = 1.219249   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.7143\n",
      "[  NORMAL ]  Iteration 323:  k_eff = 1.219267   res = 1.452E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 3.0000\n",
      "[  NORMAL ]  Iteration 324:  k_eff = 1.219284   res = 1.549E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 1.0667\n",
      "[  NORMAL ]  Iteration 325:  k_eff = 1.219301   res = 1.258E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.8125\n",
      "[  NORMAL ]  Iteration 326:  k_eff = 1.219318   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.5385\n",
      "[  NORMAL ]  Iteration 327:  k_eff = 1.219334   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.7143\n",
      "[  NORMAL ]  Iteration 328:  k_eff = 1.219349   res = 1.645E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 3.4000\n",
      "[  NORMAL ]  Iteration 329:  k_eff = 1.219365   res = 1.258E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.7647\n",
      "[  NORMAL ]  Iteration 330:  k_eff = 1.219380   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.0769\n",
      "[  NORMAL ]  Iteration 331:  k_eff = 1.219394   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 6.0000\n",
      "[  NORMAL ]  Iteration 332:  k_eff = 1.219409   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.6667\n",
      "[  NORMAL ]  Iteration 333:  k_eff = 1.219422   res = 1.161E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 3.0000\n",
      "[  NORMAL ]  Iteration 334:  k_eff = 1.219436   res = 1.065E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.9167\n",
      "[  NORMAL ]  Iteration 335:  k_eff = 1.219449   res = 1.065E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 1.0000\n",
      "[  NORMAL ]  Iteration 336:  k_eff = 1.219462   res = 1.452E-07  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 1.3636\n",
      "[  NORMAL ]  Iteration 337:  k_eff = 1.219475   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.4667\n",
      "[  NORMAL ]  Iteration 338:  k_eff = 1.219487   res = 5.807E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.8571\n",
      "[  NORMAL ]  Iteration 339:  k_eff = 1.219499   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.1667\n",
      "[  NORMAL ]  Iteration 340:  k_eff = 1.219510   res = 6.775E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 7.0000\n",
      "[  NORMAL ]  Iteration 341:  k_eff = 1.219522   res = 1.936E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.2857\n",
      "[  NORMAL ]  Iteration 342:  k_eff = 1.219533   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 2.5000\n",
      "[  NORMAL ]  Iteration 343:  k_eff = 1.219543   res = 3.872E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.8000\n",
      "[  NORMAL ]  Iteration 344:  k_eff = 1.219554   res = 8.711E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 2.2500\n",
      "[  NORMAL ]  Iteration 345:  k_eff = 1.219564   res = 4.840E-08  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 1 D.R. = 0.5556\n",
      "[  NORMAL ]  Iteration 346:  k_eff = 1.219574   res = 9.679E-09  delta-k (pcm)\n",
      "[  NORMAL ]  ...  = 0 D.R. = 0.2000\n"
     ]
    }
   ],
   "source": [
    "# Generate tracks for OpenMOC\n",
    "track_generator = openmoc.TrackGenerator(openmoc_geometry, num_azim=128, azim_spacing=0.1)\n",
    "track_generator.generateTracks()\n",
    "\n",
    "# Run OpenMOC\n",
    "solver = openmoc.CPUSolver(track_generator)\n",
    "solver.computeEigenvalue()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "openmc keff = 1.223594\n",
      "openmoc keff = 1.219574\n",
      "bias [pcm]: -402.0\n"
     ]
    }
   ],
   "source": [
    "# Print report of keff and bias with OpenMC\n",
    "openmoc_keff = solver.getKeff()\n",
    "openmc_keff = sp.keff.n\n",
    "bias = (openmoc_keff - openmc_keff) * 1e5\n",
    "\n",
    "print('openmc keff = {0:1.6f}'.format(openmc_keff))\n",
    "print('openmoc keff = {0:1.6f}'.format(openmoc_keff))\n",
    "print('bias [pcm]: {0:1.1f}'.format(bias))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "There is a non-trivial bias in both the 2-group and 8-group cases. In the case of a pin cell, one can show that these biases do not converge to <100 pcm with more particle histories. For heterogeneous geometries, additional measures must be taken to address the following three sources of bias:\n",
    "\n",
    "* Appropriate transport-corrected cross sections\n",
    "* Spatial discretization of OpenMOC's mesh\n",
    "* Constant-in-angle multi-group cross sections"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Visualizing MGXS Data"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "It is often insightful to generate visual depictions of multi-group cross sections. There are many different types of plots which may be useful for multi-group cross section visualization, only a few of which will be shown here for enrichment and inspiration.\n",
    "\n",
    "One particularly useful visualization is a comparison of the continuous-energy and multi-group cross sections for a particular nuclide and reaction type. We illustrate one option for generating such plots with the use of the `openmc.plotter` module to plot continuous-energy cross sections from the openly available cross section library distributed by NNDC.\n",
    "\n",
    "The MGXS data can also be plotted using the openmc.plot_xs command, however we will do this manually here to show how the openmc.Mgxs.get_xs method can be used to obtain data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(1e-05, 20000000.0)"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "[encoded data, 379880 characters]",
      "text/plain": [
       "<Figure size 3200x2400 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Create a figure of the U-235 continuous-energy fission cross section\n",
    "fig = openmc.plot_xs('U235', ['fission'])\n",
    "\n",
    "# Get the axis to use for plotting the MGXS\n",
    "ax = fig.gca()\n",
    "\n",
    "# Extract energy group bounds and MGXS values to plot\n",
    "fission = xs_library[fuel_cell.id]['fission']\n",
    "energy_groups = fission.energy_groups\n",
    "x = energy_groups.group_edges\n",
    "y = fission.get_xs(nuclides=['U235'], order_groups='decreasing', xs_type='micro')\n",
    "y = np.squeeze(y)\n",
    "\n",
    "# Fix low energy bound\n",
    "x[0] = 1.e-5\n",
    "\n",
    "# Extend the mgxs values array for matplotlib's step plot\n",
    "y = np.insert(y, 0, y[0])\n",
    "\n",
    "# Create a step plot for the MGXS\n",
    "ax.plot(x, y, drawstyle='steps', color='r', linewidth=3)\n",
    "\n",
    "ax.set_title('U-235 Fission Cross Section')\n",
    "ax.legend(['Continuous', 'Multi-Group'])\n",
    "ax.set_xlim((x.min(), x.max()))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Another useful type of illustration is scattering matrix sparsity structures. First, we extract Pandas `DataFrames` for the H-1 and O-16 scattering matrices."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Construct a Pandas DataFrame for the microscopic nu-scattering matrix\n",
    "nuscatter = xs_library[moderator_cell.id]['nu-scatter']\n",
    "df = nuscatter.get_pandas_dataframe(xs_type='micro')\n",
    "\n",
    "# Slice DataFrame in two for each nuclide's mean values\n",
    "h1 = df[df['nuclide'] == 'H1']['mean']\n",
    "o16 = df[df['nuclide'] == 'O16']['mean']\n",
    "\n",
    "# Cast DataFrames as NumPy arrays\n",
    "h1 = h1.values\n",
    "o16 = o16.values\n",
    "\n",
    "# Reshape arrays to 2D matrix for plotting\n",
    "h1.shape = (fine_groups.num_groups, fine_groups.num_groups)\n",
    "o16.shape = (fine_groups.num_groups, fine_groups.num_groups)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Matplotlib's `imshow` routine can be used to plot the matrices to illustrate their sparsity structures."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "[encoded data, 151468 characters]",
      "text/plain": [
       "<Figure size 3200x2400 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Create plot of the H-1 scattering matrix\n",
    "fig = plt.subplot(121)\n",
    "fig.imshow(h1, interpolation='nearest', cmap='jet')\n",
    "plt.title('H-1 Scattering Matrix')\n",
    "plt.xlabel('Group Out')\n",
    "plt.ylabel('Group In')\n",
    "\n",
    "# Create plot of the O-16 scattering matrix\n",
    "fig2 = plt.subplot(122)\n",
    "fig2.imshow(o16, interpolation='nearest', cmap='jet')\n",
    "plt.title('O-16 Scattering Matrix')\n",
    "plt.xlabel('Group Out')\n",
    "plt.ylabel('Group In')\n",
    "\n",
    "# Show the plot on screen\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [],
   "source": [
    "# close statepoint file to release HDF5 file handles\n",
    "sp.close()"
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.11.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
