GlobOptTest.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "from pycalphad.model import TestModel\n",
    "import pycalphad\n",
    "import pycalphad.variables as v\n",
    "from pycalphad.core.halton import halton\n",
    "import numpy as np\n",
    "import itertools\n",
    "rng = np.random.RandomState(1769)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "def sample_quasi(comp_count, pdof=10):\n",
    "    \"\"\"\n",
    "    Sample 'pdof * (sum(comp_count) - len(comp_count))' points in\n",
    "    composition space for the sublattice configuration specified\n",
    "    by 'comp_count'. Points are sampled quasi-randomly from a Halton sequence.\n",
    "    A Halton sequence is like a uniform random distribution, but the\n",
    "    result will always be the same for a given 'comp_count' and 'pdof'.\n",
    "    Note: For systems with only one component, only one point will be\n",
    "    returned, regardless of 'pdof'. This is because the degrees of freedom\n",
    "    are zero for that case.\n",
    "\n",
    "    Parameters\n",
    "    ----------\n",
    "    comp_count : list\n",
    "        Number of components in each sublattice.\n",
    "    pdof : int\n",
    "        Number of points to sample per degree of freedom.\n",
    "\n",
    "    Returns\n",
    "    -------\n",
    "    ndarray of generated points satisfying the mass balance.\n",
    "\n",
    "    Examples\n",
    "    --------\n",
    "    >>> comps = [8,1] # 8 components in sublattice 1; only 1 in sublattice 2\n",
    "    >>> pts = point_sample(comps, pdof=20) # 7 d.o.f, returns a 140x7 ndarray\n",
    "    \"\"\"\n",
    "    # Generate Halton sequence with appropriate dimensions and size\n",
    "    pts = halton(sum(comp_count),\n",
    "                 pdof * (sum(comp_count) - len(comp_count)), scramble=True)\n",
    "    # Convert low-discrepancy sequence to normalized exponential\n",
    "    # This will be uniformly distributed over the simplices\n",
    "    pts = -np.log(pts)\n",
    "    cur_idx = 0\n",
    "    for ctx in comp_count:\n",
    "        end_idx = cur_idx + ctx\n",
    "        pts[:, cur_idx:end_idx] /= pts[:, cur_idx:end_idx].sum(axis=1)[:, None]\n",
    "        cur_idx = end_idx\n",
    "\n",
    "    if len(pts) == 0:\n",
    "        pts = np.atleast_2d([1] * len(comp_count))\n",
    "    return pts\n",
    "def sample_pseudo(comp_count, pdof=10):\n",
    "    \"\"\"\n",
    "    Sample 'pdof * (sum(comp_count) - len(comp_count))' points in\n",
    "    composition space for the sublattice configuration specified\n",
    "    by 'comp_count'. Points are sampled pseudo-randomly from a symmetric\n",
    "    Dirichlet distribution.\n",
    "    Note: For systems with only one component, only one point will be\n",
    "    returned, regardless of 'pdof'. This is because the degrees of freedom\n",
    "    are zero for that case.\n",
    "\n",
    "    Parameters\n",
    "    ----------\n",
    "    comp_count : list\n",
    "        Number of components in each sublattice.\n",
    "    pdof : int\n",
    "        Number of points to sample per degree of freedom.\n",
    "\n",
    "    Returns\n",
    "    -------\n",
    "    ndarray of generated points satisfying the mass balance.\n",
    "\n",
    "    Examples\n",
    "    --------\n",
    "    >>> comps = [8,1] # 8 components in sublattice 1; only 1 in sublattice 2\n",
    "    >>> pts = point_sample(comps, pdof=20) # 7 d.o.f, returns a 140x7 ndarray\n",
    "    \"\"\"\n",
    "    # Generate Halton sequence with appropriate dimensions and size\n",
    "    pts = np.empty((pdof * (sum(comp_count) - len(comp_count)), sum(comp_count), ))\n",
    "    cur_idx = 0\n",
    "    for ctx in comp_count:\n",
    "        end_idx = cur_idx + ctx\n",
    "        pts[:, cur_idx:end_idx] = rng.dirichlet(np.ones(ctx, dtype=np.int),\n",
    "                                                size=(pdof * (sum(comp_count) - len(comp_count))))\n",
    "        cur_idx = end_idx\n",
    "    if len(pts) == 0:\n",
    "        pts = np.atleast_2d([1] * len(comp_count))\n",
    "    return pts\n",
    "def sample_baduniform(comp_count, pdof=10):\n",
    "    \"\"\"\n",
    "    Sample 'pdof * (sum(comp_count) - len(comp_count))' points in\n",
    "    composition space for the sublattice configuration specified\n",
    "    by 'comp_count'. Points are sampled pseudo-randomly from a uniform\n",
    "    distribution and then normalized to unity.\n",
    "    Note: For systems with only one component, only one point will be\n",
    "    returned, regardless of 'pdof'. This is because the degrees of freedom\n",
    "    are zero for that case.\n",
    "\n",
    "    Parameters\n",
    "    ----------\n",
    "    comp_count : list\n",
    "        Number of components in each sublattice.\n",
    "    pdof : int\n",
    "        Number of points to sample per degree of freedom.\n",
    "\n",
    "    Returns\n",
    "    -------\n",
    "    ndarray of generated points satisfying the mass balance.\n",
    "\n",
    "    Examples\n",
    "    --------\n",
    "    >>> comps = [8,1] # 8 components in sublattice 1; only 1 in sublattice 2\n",
    "    >>> pts = point_sample(comps, pdof=20) # 7 d.o.f, returns a 140x7 ndarray\n",
    "    \"\"\"\n",
    "    # Generate Halton sequence with appropriate dimensions and size\n",
    "    pts = np.empty((pdof * (sum(comp_count) - len(comp_count)), sum(comp_count), ))\n",
    "    cur_idx = 0\n",
    "    for ctx in comp_count:\n",
    "        end_idx = cur_idx + ctx\n",
    "        pts[:, cur_idx:end_idx] = rng.uniform(size=(pdof * (sum(comp_count) - len(comp_count)), ctx))\n",
    "        pts[:, cur_idx:end_idx] /= pts[:, cur_idx:end_idx].sum(axis=1)[:, None]\n",
    "        cur_idx = end_idx\n",
    "    if len(pts) == 0:\n",
    "        pts = np.atleast_2d([1] * len(comp_count))\n",
    "    return pts\n",
    "def sample_fixed_grid(comp_count, pdof=10):\n",
    "    \"\"\"\n",
    "    Sample _roughly_ 'pdof * (sum(comp_count) - len(comp_count))' points in\n",
    "    composition space for the sublattice configuration specified\n",
    "    by 'comp_count'. Points are generated on a uniform grid by rejection sampling.\n",
    "    Note: For systems with only one component, only one point will be\n",
    "    returned, regardless of 'pdof'. This is because the degrees of freedom\n",
    "    are zero for that case.\n",
    "\n",
    "    Parameters\n",
    "    ----------\n",
    "    comp_count : list\n",
    "        Number of components in each sublattice.\n",
    "    pdof : int\n",
    "        Number of points to sample per degree of freedom.\n",
    "\n",
    "    Returns\n",
    "    -------\n",
    "    ndarray of generated points satisfying the mass balance.\n",
    "    \"\"\"\n",
    "\n",
    "    print(result_simplex.shape)\n",
    "    for ctx in comp_count:\n",
    "        n = ctx # number of components\n",
    "        points_per_edge = 10\n",
    "        x = np.linspace(1e-12, 1-1e-12, points_per_edge)\n",
    "        result = np.array(np.meshgrid(*itertools.repeat(x, n-1)))\n",
    "        result_simplex = result[:, result.sum(axis=0) <= 1]\n",
    "        # Add dependent component\n",
    "        temp = np.empty((result_simplex.shape[0]+1, result_simplex.shape[1]))\n",
    "        temp[:-1, :] = result_simplex\n",
    "        temp[-1, :] = 1 - result_simplex.sum(axis=0)\n",
    "        result_simplex = temp\n",
    "def test_opt(num_components, num_samples, sampler, solutions):\n",
    "    result_array = np.empty((num_samples))\n",
    "    comps = ['X{}'.format(str(i+1).zfill(len(str(num_components)))) for i in range(num_components)]\n",
    "    print('{} COMPONENTS'.format(num_components))\n",
    "    for idx in range(num_samples):\n",
    "        print('TRIAL {}/{}'.format(idx+1, num_samples))\n",
    "        solution = solutions[idx]\n",
    "        tm = TestModel(None, comps, 'TEST', solution=solution, kmax=2)\n",
    "        dbf = pycalphad.Database()\n",
    "        dbf.add_phase('TEST', {}, [1])\n",
    "        dbf.add_phase_constituents('TEST', [comps])\n",
    "        conds = {v.P:101325, v.T:300}\n",
    "        conds.update({v.X(i):sol for i, sol in list(zip(comps, solution))[:-1]})\n",
    "        grid_opts = {'sampler': sampler, 'pdens': 100*num_components}\n",
    "        eq = pycalphad.equilibrium(dbf, comps, 'TEST', conds, model=tm, verbose=False, grid_opts=grid_opts)\n",
    "        if eq.GM < 1e-2:\n",
    "            result_array[idx] = eq.attrs['solve_iterations']\n",
    "        else:\n",
    "            print(eq.GM)\n",
    "            print(eq.X.values)\n",
    "            print(solution)\n",
    "            result_array[idx] = np.nan\n",
    "    return result_array"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false,
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "PSEUDO\n",
      "2 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "QUASI\n",
      "2 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "NORMUNITY\n",
      "2 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "PSEUDO\n",
      "3 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "QUASI\n",
      "3 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "NORMUNITY\n",
      "3 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "PSEUDO\n",
      "4 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "QUASI\n",
      "4 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "NORMUNITY\n",
      "4 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "PSEUDO\n",
      "5 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "QUASI\n",
      "5 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "NORMUNITY\n",
      "5 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1)>\n",
      "array([[[[[[ 7238.68588873]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.6361\n",
      "  * X_X2     (X_X2) float64 0.0401\n",
      "  * X_X3     (X_X3) float64 0.02001\n",
      "  * X_X4     (X_X4) float64 0.2639\n",
      "[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "          1.00000000e-12]\n",
      "       [  6.08307823e-01   1.34597549e-01   7.48176512e-06   2.46601915e-01\n",
      "          1.04852314e-02]\n",
      "       [  5.07932112e-01   3.97586318e-06   9.81424797e-02   2.84523475e-01\n",
      "          1.09397957e-01]\n",
      "       [  3.44334551e-01   5.76482064e-02   3.19692179e-05   5.52583078e-01\n",
      "          4.54021957e-02]\n",
      "       [  4.53327847e-01   2.71505425e-04   2.61954220e-02   3.65522236e-01\n",
      "          1.54682989e-01]]]]]]]]\n",
      "[ 0.63614401  0.04009614  0.02001309  0.26387137  0.0398754 ]\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "PSEUDO\n",
      "6 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>\n",
      "array([[[[[[[ 14325.51248918]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.005255\n",
      "  * X_X2     (X_X2) float64 0.01646\n",
      "  * X_X3     (X_X3) float64 0.2695\n",
      "  * X_X4     (X_X4) float64 0.4737\n",
      "  * X_X5     (X_X5) float64 0.0004482\n",
      "[[[[[[[[[  2.28994549e-02   2.63686386e-02   1.41746895e-01\n",
      "           5.49870809e-01   9.80551357e-07   2.59113222e-01]\n",
      "        [  2.52882907e-03   1.33567198e-01   4.90304457e-01\n",
      "           3.24080139e-01   2.07881999e-03   4.74405570e-02]\n",
      "        [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "           1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "        [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "           1.00000000e+00   1.00000000e-12   1.00000000e-12]\n",
      "        [  8.37532226e-03   4.41192924e-02   4.00108276e-01\n",
      "           4.77598536e-01   1.94258016e-02   5.03727722e-02]\n",
      "        [  2.11603941e-03   2.23942894e-02   3.30205564e-01\n",
      "           2.35238421e-01   3.94165622e-04   4.09651520e-01]]]]]]]]]\n",
      "[  5.25474112e-03   1.64646957e-02   2.69451682e-01   4.73661002e-01\n",
      "   4.48199795e-04   2.34719679e-01]\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>\n",
      "array([[[[[[[ 12555.60915059]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.336\n",
      "  * X_X2     (X_X2) float64 0.1077\n",
      "  * X_X3     (X_X3) float64 0.00534\n",
      "  * X_X4     (X_X4) float64 0.09086\n",
      "  * X_X5     (X_X5) float64 0.4575\n",
      "[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "           1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "        [  2.66080538e-01   4.09157822e-01   2.62032915e-05\n",
      "           4.10971526e-02   2.79313777e-01   4.32450698e-03]\n",
      "        [  4.82810025e-01   3.89000668e-02   4.55702050e-02\n",
      "           2.24143598e-02   4.09892565e-01   4.12778463e-04]\n",
      "        [  5.14111255e-01   5.31863628e-02   3.53283697e-04\n",
      "           2.13229225e-01   2.19082119e-01   3.77545277e-05]\n",
      "        [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "           1.00000000e-12   1.00000000e+00   1.00000000e-12]\n",
      "        [  4.36960260e-01   1.69089519e-02   1.78538181e-07\n",
      "           1.55286086e-01   3.76765658e-01   1.40788660e-02]]]]]]]]]\n",
      "[ 0.3360238   0.10770806  0.0053402   0.09086331  0.45751585  0.00254878]\n",
      "TRIAL 20/20\n",
      "QUASI\n",
      "6 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>\n",
      "array([[[[[[[ 23315.20975776]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.005255\n",
      "  * X_X2     (X_X2) float64 0.01646\n",
      "  * X_X3     (X_X3) float64 0.2695\n",
      "  * X_X4     (X_X4) float64 0.4737\n",
      "  * X_X5     (X_X5) float64 0.0004482\n",
      "[[[[[[[[[  4.04723993e-02   1.17201777e-03   1.02783905e-01\n",
      "           6.81601795e-01   1.56040473e-06   1.73968322e-01]\n",
      "        [  1.22040661e-02   4.83666910e-02   2.82503668e-01\n",
      "           2.93110740e-01   1.21864049e-04   3.63692972e-01]\n",
      "        [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "           1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "        [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "           1.00000000e+00   1.00000000e-12   1.00000000e-12]\n",
      "        [  1.49666711e-02   7.24832569e-03   2.85433043e-01\n",
      "           2.78381683e-01   6.88639087e-03   4.07083886e-01]\n",
      "        [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "           1.00000000e-12   1.00000000e-12   1.00000000e+00]]]]]]]]]\n",
      "[  5.25474112e-03   1.64646957e-02   2.69451682e-01   4.73661002e-01\n",
      "   4.48199795e-04   2.34719679e-01]\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>\n",
      "array([[[[[[[ 9395.60280102]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.336\n",
      "  * X_X2     (X_X2) float64 0.1077\n",
      "  * X_X3     (X_X3) float64 0.00534\n",
      "  * X_X4     (X_X4) float64 0.09086\n",
      "  * X_X5     (X_X5) float64 0.4575\n",
      "[[[[[[[[[  6.77486496e-01   3.45970907e-02   1.69615539e-02\n",
      "           2.36129985e-01   3.48244002e-02   4.74495647e-07]\n",
      "        [  5.48283355e-02   4.19898050e-01   5.63290984e-03\n",
      "           1.91002461e-01   3.28637368e-01   8.74960982e-07]\n",
      "        [  5.36472193e-01   2.81966942e-01   1.73428686e-02\n",
      "           7.60516103e-03   1.56611263e-01   1.57277379e-06]\n",
      "        [  2.46457641e-01   1.48773099e-02   2.12470196e-02\n",
      "           1.29319472e-01   5.88097310e-01   1.24760401e-06]\n",
      "        [  2.64258268e-01   3.45666704e-03   3.82536937e-03\n",
      "           4.30175287e-03   7.24156040e-01   1.90222818e-06]\n",
      "        [  6.01352316e-01   1.31506394e-01   5.58706897e-10\n",
      "           1.35043948e-01   1.23548989e-01   8.54835207e-03]]]]]]]]]\n",
      "[ 0.3360238   0.10770806  0.0053402   0.09086331  0.45751585  0.00254878]\n",
      "TRIAL 20/20\n",
      "NORMUNITY\n",
      "6 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>\n",
      "array([[[[[[[ 13818.1806527]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.005255\n",
      "  * X_X2     (X_X2) float64 0.01646\n",
      "  * X_X3     (X_X3) float64 0.2695\n",
      "  * X_X4     (X_X4) float64 0.4737\n",
      "  * X_X5     (X_X5) float64 0.0004482\n",
      "[[[[[[[[[  4.61554052e-02   8.27797832e-02   2.27583028e-01\n",
      "           3.70677232e-01   6.75748450e-04   2.72128803e-01]\n",
      "        [  1.52135241e-03   3.08921218e-01   2.90501438e-01\n",
      "           2.26290076e-01   2.94584845e-07   1.72765621e-01]\n",
      "        [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "           1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "        [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "           1.00000000e+00   1.00000000e-12   1.00000000e-12]\n",
      "        [             nan              nan              nan\n",
      "                      nan              nan              nan]\n",
      "        [             nan              nan              nan\n",
      "                      nan              nan              nan]]]]]]]]]\n",
      "[  5.25474112e-03   1.64646957e-02   2.69451682e-01   4.73661002e-01\n",
      "   4.48199795e-04   2.34719679e-01]\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>\n",
      "array([[[[[[[ 12488.49977685]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.5454\n",
      "  * X_X2     (X_X2) float64 0.005528\n",
      "  * X_X3     (X_X3) float64 0.135\n",
      "  * X_X4     (X_X4) float64 0.08867\n",
      "  * X_X5     (X_X5) float64 0.0082\n",
      "[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "           1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "        [  3.72980817e-01   4.65234656e-02   1.31283865e-01\n",
      "           9.55051871e-02   2.48863598e-03   3.51218029e-01]\n",
      "        [  4.56802832e-01   1.49080445e-03   3.94107230e-01\n",
      "           4.37233356e-02   3.29358153e-06   1.03872504e-01]\n",
      "        [  3.58038337e-01   1.11561921e-05   2.66068523e-02\n",
      "           2.91603860e-01   2.30479354e-02   3.00691859e-01]\n",
      "        [  3.24653212e-01   2.41257088e-05   2.07002562e-01\n",
      "           9.34844005e-02   4.24244391e-02   3.32411260e-01]\n",
      "        [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "           1.00000000e-12   1.00000000e-12   1.00000000e+00]]]]]]]]]\n",
      "[ 0.54544083  0.00552811  0.13496374  0.0886669   0.00820048  0.21719994]\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1)>\n",
      "array([[[[[[[ 21452.94040026]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.336\n",
      "  * X_X2     (X_X2) float64 0.1077\n",
      "  * X_X3     (X_X3) float64 0.00534\n",
      "  * X_X4     (X_X4) float64 0.09086\n",
      "  * X_X5     (X_X5) float64 0.4575\n",
      "[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "           1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "        [  2.25256905e-01   4.23827402e-01   2.56755900e-02\n",
      "           5.16846078e-02   2.73535480e-01   2.00150336e-05]\n",
      "        [  5.27543048e-01   2.21567822e-02   8.42526024e-02\n",
      "           5.54984906e-02   3.10532561e-01   1.65154555e-05]\n",
      "        [  1.08764846e-01   2.69558227e-01   2.12049496e-03\n",
      "           2.93867309e-01   3.25688992e-01   1.30132014e-07]\n",
      "        [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "           1.00000000e-12   1.00000000e+00   1.00000000e-12]\n",
      "        [  2.57116944e-01   1.70649815e-01   2.50484573e-08\n",
      "           1.67493091e-01   3.85881083e-01   1.88590413e-02]]]]]]]]]\n",
      "[ 0.3360238   0.10770806  0.0053402   0.09086331  0.45751585  0.00254878]\n",
      "TRIAL 20/20\n",
      "PSEUDO\n",
      "7 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1)>\n",
      "array([[[[[[[[ 12143.2016423]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.2864\n",
      "  * X_X2     (X_X2) float64 0.01361\n",
      "  * X_X3     (X_X3) float64 0.3705\n",
      "  * X_X4     (X_X4) float64 0.1148\n",
      "  * X_X5     (X_X5) float64 0.07925\n",
      "  * X_X6     (X_X6) float64 0.1344\n",
      "[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12]\n",
      "         [  4.06795681e-01   8.15709375e-02   3.53371408e-01\n",
      "            1.07040370e-02   2.33356673e-03   1.42483583e-01\n",
      "            2.74078647e-03]\n",
      "         [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12]\n",
      "         [  3.20983587e-01   2.11604841e-02   1.56648979e-01\n",
      "            7.38654361e-02   2.98633773e-01   1.28670610e-01\n",
      "            3.71309076e-05]\n",
      "         [  3.41347449e-01   4.21348088e-02   2.56336192e-01\n",
      "            5.33628045e-02   2.00197490e-01   1.02216075e-01\n",
      "            4.40518132e-03]\n",
      "         [  8.19245733e-02   6.94143385e-03   4.06373279e-01\n",
      "            4.20412062e-02   4.10129593e-02   4.21576857e-01\n",
      "            1.29691339e-04]\n",
      "         [  3.26122702e-01   1.27546656e-02   3.10233766e-01\n",
      "            2.25030638e-01   1.77543621e-02   1.05679822e-01\n",
      "            2.42404372e-03]]]]]]]]]]\n",
      "[ 0.2863812   0.0136059   0.37047733  0.114804    0.07925465  0.13435309\n",
      "  0.00112385]\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "QUASI\n",
      "7 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1)>\n",
      "array([[[[[[[[ 20768.12436207]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.2956\n",
      "  * X_X2     (X_X2) float64 0.000142\n",
      "  * X_X3     (X_X3) float64 0.06756\n",
      "  * X_X4     (X_X4) float64 0.03174\n",
      "  * X_X5     (X_X5) float64 0.1409\n",
      "  * X_X6     (X_X6) float64 0.2146\n",
      "[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12]\n",
      "         [  1.82319327e-01   1.65483781e-03   1.06801155e-01\n",
      "            1.03120283e-02   1.06629194e-01   2.97884072e-01\n",
      "            2.94399385e-01]\n",
      "         [  4.01622708e-01   4.33715605e-06   2.96283160e-01\n",
      "            1.31752327e-02   7.12994466e-02   2.01019883e-01\n",
      "            1.65952320e-02]\n",
      "         [  1.41232782e-01   5.85016252e-05   4.75439332e-02\n",
      "            7.55396312e-02   2.08764150e-01   2.62837737e-01\n",
      "            2.64023266e-01]\n",
      "         [  1.30011937e-01   2.14423163e-05   4.96361670e-02\n",
      "            3.78632117e-02   4.34443001e-01   8.79283090e-02\n",
      "            2.60095933e-01]\n",
      "         [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "            1.00000000e-12]\n",
      "         [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e+00]]]]]]]]]]\n",
      "[  2.95591933e-01   1.41958677e-04   6.75577666e-02   3.17416602e-02\n",
      "   1.40863083e-01   2.14621331e-01   2.49482267e-01]\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "NORMUNITY\n",
      "7 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1)>\n",
      "array([[[[[[[[ 30129.37048082]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.2956\n",
      "  * X_X2     (X_X2) float64 0.000142\n",
      "  * X_X3     (X_X3) float64 0.06756\n",
      "  * X_X4     (X_X4) float64 0.03174\n",
      "  * X_X5     (X_X5) float64 0.1409\n",
      "  * X_X6     (X_X6) float64 0.2146\n",
      "[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12]\n",
      "         [  2.53562387e-01   1.34072367e-03   4.46139593e-04\n",
      "            3.50878202e-02   3.31969453e-01   2.32949790e-01\n",
      "            1.44643687e-01]\n",
      "         [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12]\n",
      "         [  2.37993183e-01   8.00158467e-05   7.41135674e-02\n",
      "            7.85937575e-02   1.33041239e-01   2.50885301e-01\n",
      "            2.25292937e-01]\n",
      "         [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "            1.00000000e-12]\n",
      "         [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "            1.00000000e-12]\n",
      "         [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e+00]]]]]]]]]]\n",
      "[  2.95591933e-01   1.41958677e-04   6.75577666e-02   3.17416602e-02\n",
      "   1.40863083e-01   2.14621331e-01   2.49482267e-01]\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1)>\n",
      "array([[[[[[[[ 15543.65030879]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.0001097\n",
      "  * X_X2     (X_X2) float64 0.2581\n",
      "  * X_X3     (X_X3) float64 0.2186\n",
      "  * X_X4     (X_X4) float64 0.04642\n",
      "  * X_X5     (X_X5) float64 0.2661\n",
      "  * X_X6     (X_X6) float64 0.1149\n",
      "[[[[[[[[[[  3.11675739e-03   3.05939450e-01   1.60983103e-01\n",
      "            7.21524643e-02   2.77349720e-01   1.90515173e-02\n",
      "            1.61406988e-01]\n",
      "         [  1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12]\n",
      "         [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12]\n",
      "         [  7.78576033e-07   2.58993420e-01   3.08435205e-01\n",
      "            1.40567349e-01   1.35289094e-01   5.48925683e-02\n",
      "            1.01821586e-01]\n",
      "         [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "            1.00000000e-12]\n",
      "         [  6.15393783e-07   2.57676965e-01   1.97439094e-01\n",
      "            5.88892870e-02   2.12131437e-01   2.12841783e-01\n",
      "            6.10208189e-02]\n",
      "         [  9.59530652e-06   2.02220894e-01   1.20179140e-01\n",
      "            2.85381378e-02   2.53040728e-01   1.78094773e-01\n",
      "            2.17916732e-01]]]]]]]]]]\n",
      "[  1.09720560e-04   2.58123592e-01   2.18574574e-01   4.64156446e-02\n",
      "   2.66086284e-01   1.14907390e-01   9.57827936e-02]\n",
      "TRIAL 20/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1)>\n",
      "array([[[[[[[[ 9276.12132551]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.5347\n",
      "  * X_X2     (X_X2) float64 0.2542\n",
      "  * X_X3     (X_X3) float64 0.08581\n",
      "  * X_X4     (X_X4) float64 0.05228\n",
      "  * X_X5     (X_X5) float64 0.01141\n",
      "  * X_X6     (X_X6) float64 0.03852\n",
      "[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12]\n",
      "         [  1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "            1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "            1.00000000e-12]\n",
      "         [  3.68784363e-01   2.99669911e-01   2.09561728e-01\n",
      "            6.16606033e-02   5.35216108e-07   3.26033780e-02\n",
      "            2.77194820e-02]\n",
      "         [  4.01529562e-01   3.44737716e-01   3.45576770e-02\n",
      "            1.18296457e-01   3.42961326e-05   9.38142043e-02\n",
      "            7.03008768e-03]\n",
      "         [  4.64047990e-01   2.68104176e-01   2.48725422e-02\n",
      "            1.01905726e-01   1.67509930e-02   2.40809625e-02\n",
      "            1.00237611e-01]\n",
      "         [  4.17028951e-01   2.89627588e-01   1.46138096e-01\n",
      "            1.06610025e-04   8.62343173e-02   4.03585411e-02\n",
      "            2.05058962e-02]\n",
      "         [  3.15252340e-01   3.25542784e-01   1.60642141e-01\n",
      "            3.87306275e-02   5.72606478e-08   9.67421065e-02\n",
      "            6.30899434e-02]]]]]]]]]]\n",
      "[ 0.53468409  0.25423942  0.08580726  0.05227689  0.01140856  0.03851943\n",
      "  0.02306436]\n",
      "PSEUDO\n",
      "8 COMPONENTS\n",
      "TRIAL 1/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>\n",
      "array([[[[[[[[[ 24925.26079685]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.132\n",
      "  * X_X2     (X_X2) float64 0.3108\n",
      "  * X_X3     (X_X3) float64 0.001656\n",
      "  * X_X4     (X_X4) float64 0.007117\n",
      "  * X_X5     (X_X5) float64 0.05303\n",
      "  * X_X6     (X_X6) float64 0.1749\n",
      "  * X_X7     (X_X7) float64 0.01697\n",
      "[[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  2.45350891e-01   2.53024277e-03   2.19646487e-02\n",
      "             6.64647951e-06   1.84315590e-01   2.60247963e-01\n",
      "             5.57802357e-04   2.85026216e-01]\n",
      "          [  2.08126209e-01   1.61910448e-01   1.31610565e-07\n",
      "             1.52374550e-02   2.59663396e-02   1.70396324e-01\n",
      "             2.38899459e-02   3.94473147e-01]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  4.56930588e-02   2.99722910e-01   5.12090565e-04\n",
      "             1.41072319e-05   2.07422880e-01   1.57307182e-01\n",
      "             9.25090928e-02   1.96818679e-01]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e+00]]]]]]]]]]]\n",
      "[ 0.13198901  0.31075186  0.00165636  0.00711656  0.05302571  0.17494918\n",
      "  0.01697014  0.30354119]\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>\n",
      "array([[[[[[[[[ 29020.59743032]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.01202\n",
      "  * X_X2     (X_X2) float64 0.002629\n",
      "  * X_X3     (X_X3) float64 0.06178\n",
      "  * X_X4     (X_X4) float64 0.1926\n",
      "  * X_X5     (X_X5) float64 0.001547\n",
      "  * X_X6     (X_X6) float64 0.2561\n",
      "  * X_X7     (X_X7) float64 0.2868\n",
      "[[[[[[[[[[[  3.06804350e-02   9.66445553e-06   2.54806365e-01\n",
      "             2.71455875e-02   4.90583095e-03   1.84318101e-02\n",
      "             2.22672374e-01   4.41347933e-01]\n",
      "          [  2.63840442e-02   8.12518019e-03   1.12363698e-01\n",
      "             2.21027081e-01   1.56394501e-06   1.77915799e-01\n",
      "             2.34147510e-01   2.20035124e-01]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  4.66985013e-03   1.90365211e-05   1.93950589e-02\n",
      "             6.56178124e-01   5.26194990e-03   1.46536685e-01\n",
      "             4.30737026e-02   1.24865593e-01]\n",
      "          [  8.44995295e-03   6.80495372e-05   1.03343849e-02\n",
      "             2.44006807e-01   1.95509978e-02   2.26380340e-01\n",
      "             3.28080674e-01   1.63128793e-01]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e+00   1.00000000e-12]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e+00]]]]]]]]]]]\n",
      "[ 0.01202487  0.00262912  0.06178035  0.19263393  0.00154735  0.25611985\n",
      "  0.28678588  0.18647865]\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "QUASI\n",
      "8 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>\n",
      "array([[[[[[[[[ 26711.88091663]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.01202\n",
      "  * X_X2     (X_X2) float64 0.002629\n",
      "  * X_X3     (X_X3) float64 0.06178\n",
      "  * X_X4     (X_X4) float64 0.1926\n",
      "  * X_X5     (X_X5) float64 0.001547\n",
      "  * X_X6     (X_X6) float64 0.2561\n",
      "  * X_X7     (X_X7) float64 0.2868\n",
      "[[[[[[[[[[[  2.68793576e-02   2.61567032e-04   9.13722020e-02\n",
      "             1.01270066e-01   1.87975902e-03   2.42202336e-01\n",
      "             1.08268326e-01   4.27866386e-01]\n",
      "          [  1.30322758e-03   2.04288838e-02   4.76188312e-02\n",
      "             3.78046249e-02   1.40714227e-07   1.86153776e-01\n",
      "             6.09296668e-01   9.73938478e-02]\n",
      "          [  1.17118219e-02   3.86411293e-04   4.18077159e-01\n",
      "             4.95621115e-02   1.32493090e-02   1.76350675e-01\n",
      "             2.82255743e-01   4.84067699e-02]\n",
      "          [  2.41304820e-02   6.55383823e-04   6.14534678e-02\n",
      "             4.07421770e-01   8.11797282e-07   1.04007343e-01\n",
      "             1.04706854e-01   2.97623887e-01]\n",
      "          [  1.10913058e-05   6.86044693e-03   2.31760827e-01\n",
      "             3.73398133e-02   2.20013913e-02   2.25373639e-01\n",
      "             3.11472468e-01   1.65180323e-01]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e+00   1.00000000e-12]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e+00]]]]]]]]]]]\n",
      "[ 0.01202487  0.00262912  0.06178035  0.19263393  0.00154735  0.25611985\n",
      "  0.28678588  0.18647865]\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "NORMUNITY\n",
      "8 COMPONENTS\n",
      "TRIAL 1/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>\n",
      "array([[[[[[[[[ 36582.43384545]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.132\n",
      "  * X_X2     (X_X2) float64 0.3108\n",
      "  * X_X3     (X_X3) float64 0.001656\n",
      "  * X_X4     (X_X4) float64 0.007117\n",
      "  * X_X5     (X_X5) float64 0.05303\n",
      "  * X_X6     (X_X6) float64 0.1749\n",
      "  * X_X7     (X_X7) float64 0.01697\n",
      "[[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  3.10938021e-01   1.53445960e-01   2.44556899e-02\n",
      "             1.01461956e-02   1.80289338e-03   2.62122842e-01\n",
      "             3.28575053e-02   2.04230892e-01]\n",
      "          [  2.18353347e-01   2.13403942e-01   4.95389178e-04\n",
      "             7.20272596e-02   1.92592930e-02   2.54290171e-01\n",
      "             6.83001625e-02   1.53870435e-01]\n",
      "          [  1.14242449e-01   1.99545072e-01   3.14120203e-04\n",
      "             5.41035957e-03   2.26866191e-01   2.47855708e-01\n",
      "             2.27734152e-02   1.82992684e-01]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  2.21314783e-01   1.93722590e-01   1.27473530e-03\n",
      "             3.61336483e-03   1.35766167e-02   2.42844156e-01\n",
      "             1.86729123e-01   1.36924631e-01]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e+00]]]]]]]]]]]\n",
      "[ 0.13198901  0.31075186  0.00165636  0.00711656  0.05302571  0.17494918\n",
      "  0.01697014  0.30354119]\n",
      "TRIAL 2/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>\n",
      "array([[[[[[[[[ 11926.98627586]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.00885\n",
      "  * X_X2     (X_X2) float64 0.04943\n",
      "  * X_X3     (X_X3) float64 0.1276\n",
      "  * X_X4     (X_X4) float64 0.01819\n",
      "  * X_X5     (X_X5) float64 0.22\n",
      "  * X_X6     (X_X6) float64 0.2321\n",
      "  * X_X7     (X_X7) float64 0.01136\n",
      "[[[[[[[[[[[  1.16752330e-02   1.04250710e-02   1.20763895e-01\n",
      "             9.25891505e-02   1.44203206e-01   3.34224302e-01\n",
      "             4.55614035e-05   2.86073581e-01]\n",
      "          [  7.11727074e-03   1.07663566e-01   3.76451511e-01\n",
      "             3.95537060e-03   6.86824701e-02   1.95410549e-01\n",
      "             1.69328577e-07   2.40719093e-01]\n",
      "          [  2.07324560e-02   1.04925040e-02   1.80631466e-01\n",
      "             1.52932000e-03   8.64077363e-02   3.87146177e-01\n",
      "             5.24730335e-02   2.60587307e-01]\n",
      "          [  1.11282438e-01   1.33940705e-02   1.84401208e-01\n",
      "             1.18986031e-02   2.23817579e-01   2.73018199e-01\n",
      "             1.30848454e-05   1.82174818e-01]\n",
      "          [  7.38646546e-02   1.24286333e-01   1.09249905e-01\n",
      "             1.45298742e-07   2.30374058e-01   2.40698386e-01\n",
      "             2.80464999e-04   2.21246052e-01]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  1.37591238e-06   5.44590777e-02   8.71072338e-02\n",
      "             2.55573236e-02   3.63978451e-01   9.67690412e-02\n",
      "             1.06368041e-02   3.61490693e-01]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e+00]]]]]]]]]]]\n",
      "[ 0.00885012  0.04943079  0.12760323  0.01819249  0.22004347  0.23211777\n",
      "  0.01136232  0.3323998 ]\n",
      "TRIAL 3/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>\n",
      "array([[[[[[[[[ 24547.35806647]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.04442\n",
      "  * X_X2     (X_X2) float64 0.3364\n",
      "  * X_X3     (X_X3) float64 0.004386\n",
      "  * X_X4     (X_X4) float64 0.03402\n",
      "  * X_X5     (X_X5) float64 0.491\n",
      "  * X_X6     (X_X6) float64 0.0195\n",
      "  * X_X7     (X_X7) float64 0.008923\n",
      "[[[[[[[[[[[  2.02814262e-01   2.67469671e-01   7.72620840e-03\n",
      "             2.85150903e-02   2.10278425e-01   4.87925142e-03\n",
      "             7.52373872e-02   2.03079704e-01]\n",
      "          [  1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  1.76930236e-02   2.95818610e-01   2.28377433e-01\n",
      "             1.77384750e-05   2.94607456e-01   4.51802180e-02\n",
      "             6.04574040e-03   1.12259781e-01]\n",
      "          [  1.03108584e-01   3.50225910e-01   2.84071408e-03\n",
      "             1.21792349e-01   3.16964852e-01   4.45465079e-02\n",
      "             1.06145239e-05   6.05104687e-02]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  7.58496988e-02   1.66031071e-01   1.76688893e-03\n",
      "             2.49963341e-04   4.67322178e-01   6.49287863e-02\n",
      "             6.45353627e-02   1.59316051e-01]\n",
      "          [  9.81651754e-02   2.40820809e-01   3.90078223e-02\n",
      "             3.19462279e-02   3.24115032e-01   8.38429704e-05\n",
      "             1.73096624e-01   9.27644660e-02]\n",
      "          [  1.02327975e-01   1.46922057e-01   1.57648926e-03\n",
      "             7.53443205e-04   3.16702839e-01   1.08870542e-05\n",
      "             9.33689101e-03   4.22369419e-01]]]]]]]]]]]\n",
      "[ 0.04442347  0.33642228  0.00438618  0.03401742  0.49100443  0.01949908\n",
      "  0.00892316  0.06132397]\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1)>\n",
      "array([[[[[[[[[ 44246.02970379]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.01202\n",
      "  * X_X2     (X_X2) float64 0.002629\n",
      "  * X_X3     (X_X3) float64 0.06178\n",
      "  * X_X4     (X_X4) float64 0.1926\n",
      "  * X_X5     (X_X5) float64 0.001547\n",
      "  * X_X6     (X_X6) float64 0.2561\n",
      "  * X_X7     (X_X7) float64 0.2868\n",
      "[[[[[[[[[[[  5.56878341e-02   3.54251350e-03   1.46870307e-02\n",
      "             2.94141587e-01   4.03335091e-04   9.17200204e-02\n",
      "             2.84088988e-01   2.55728691e-01]\n",
      "          [  6.89594046e-03   3.17844302e-02   3.41487518e-01\n",
      "             2.61129143e-01   2.82067120e-04   1.68831556e-01\n",
      "             4.73453276e-02   1.42244018e-01]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  1.06276494e-04   4.45265368e-03   2.45366581e-01\n",
      "             2.30643063e-01   3.28399307e-02   1.66489557e-01\n",
      "             1.00770205e-01   2.19331733e-01]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "             1.00000000e-12   1.00000000e-12]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e+00   1.00000000e-12]\n",
      "          [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "             1.00000000e-12   1.00000000e+00]]]]]]]]]]]\n",
      "[ 0.01202487  0.00262912  0.06178035  0.19263393  0.00154735  0.25611985\n",
      "  0.28678588  0.18647865]\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "PSEUDO\n",
      "9 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>\n",
      "array([[[[[[[[[[ 22351.08264421]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.00283\n",
      "  * X_X2     (X_X2) float64 0.158\n",
      "  * X_X3     (X_X3) float64 0.006395\n",
      "  * X_X4     (X_X4) float64 0.009869\n",
      "  * X_X5     (X_X5) float64 0.2794\n",
      "  * X_X6     (X_X6) float64 0.2289\n",
      "  * X_X7     (X_X7) float64 0.06501\n",
      "  * X_X8     (X_X8) float64 0.1835\n",
      "[[[[[[[[[[[[  1.88635980e-05   4.14667245e-01   1.01248930e-03\n",
      "              1.15237908e-02   1.65610570e-01   2.28220611e-01\n",
      "              1.15068280e-01   3.20651772e-02   3.18129734e-02]\n",
      "           [  7.47166899e-05   2.02856552e-01   1.89631513e-02\n",
      "              2.88676023e-03   2.96893520e-01   1.58372861e-01\n",
      "              1.66552429e-02   1.23812493e-01   1.79484703e-01]\n",
      "           [  7.48148142e-03   7.06868878e-02   2.54706622e-03\n",
      "              9.20269885e-02   1.65583645e-01   4.49637066e-01\n",
      "              1.47156951e-02   1.56861469e-01   4.04597013e-02]\n",
      "           [  7.85949870e-03   8.28920592e-03   6.85691087e-06\n",
      "              4.38913540e-02   2.68776444e-01   2.74402939e-01\n",
      "              9.43394301e-02   1.10110475e-01   1.92323797e-01]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  2.05863462e-04   2.48613957e-02   1.57773917e-02\n",
      "              2.66265783e-03   3.23541866e-01   2.17030867e-01\n",
      "              1.61565756e-01   1.52883514e-01   1.01470688e-01]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e+00   1.00000000e-12]\n",
      "           [  1.27827640e-02   2.63674075e-01   9.71426803e-08\n",
      "              9.84036312e-03   2.00336310e-02   4.57602593e-01\n",
      "              8.81016931e-02   1.07389786e-01   4.05749970e-02]]]]]]]]]]]]\n",
      "[ 0.00282968  0.15800187  0.00639466  0.00986923  0.27939811  0.22891086\n",
      "  0.06500936  0.18352608  0.06606016]\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>\n",
      "array([[[[[[[[[[ 27853.15998024]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.2817\n",
      "  * X_X2     (X_X2) float64 0.0007679\n",
      "  * X_X3     (X_X3) float64 0.3302\n",
      "  * X_X4     (X_X4) float64 0.2655\n",
      "  * X_X5     (X_X5) float64 0.01434\n",
      "  * X_X6     (X_X6) float64 0.009814\n",
      "  * X_X7     (X_X7) float64 0.008701\n",
      "  * X_X8     (X_X8) float64 0.02087\n",
      "[[[[[[[[[[[[  4.53072223e-01   1.00000000e-12   3.10862158e-01\n",
      "              3.90503127e-02   2.16292438e-02   2.73578118e-02\n",
      "              2.85883642e-02   1.97978570e-02   9.96420292e-02]\n",
      "           [  7.61935055e-01   2.34086371e-03   1.19631493e-01\n",
      "              4.97718074e-02   1.32086261e-02   1.10246063e-05\n",
      "              8.66218111e-03   1.52039705e-02   2.92349777e-02]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  4.63165792e-02   1.41015257e-04   3.26997956e-01\n",
      "              4.61739831e-01   1.02728291e-01   1.90861723e-02\n",
      "              1.73606934e-02   3.14107671e-06   2.56263207e-02]\n",
      "           [  2.49153940e-01   3.16829138e-02   3.66476705e-01\n",
      "              1.89254098e-01   3.44386537e-02   4.00141249e-04\n",
      "              4.45125538e-02   2.84115954e-02   5.56693994e-02]\n",
      "           [  2.34915758e-02   6.82461048e-08   5.24923059e-01\n",
      "              2.02424382e-01   7.41610300e-03   1.02710726e-02\n",
      "              2.19432630e-01   5.58269889e-03   6.45841078e-03]\n",
      "           [  2.64126330e-01   1.95286539e-06   5.56801498e-01\n",
      "              5.82959408e-02   6.76872063e-03   2.47907623e-02\n",
      "              2.02202073e-03   6.10893330e-02   2.61034417e-02]\n",
      "           [  1.99322244e-01   3.28286823e-06   4.46886686e-01\n",
      "              8.44958981e-02   3.33925162e-02   8.00680450e-03\n",
      "              1.01966388e-02   2.11664072e-02   1.96529522e-01]]]]]]]]]]]]\n",
      "[ 0.28171799  0.00076793  0.33020705  0.26546957  0.01433725  0.00981422\n",
      "  0.00870094  0.02086778  0.06811727]\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "QUASI\n",
      "9 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>\n",
      "array([[[[[[[[[[ 27130.52965317]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.2817\n",
      "  * X_X2     (X_X2) float64 0.0007679\n",
      "  * X_X3     (X_X3) float64 0.3302\n",
      "  * X_X4     (X_X4) float64 0.2655\n",
      "  * X_X5     (X_X5) float64 0.01434\n",
      "  * X_X6     (X_X6) float64 0.009814\n",
      "  * X_X7     (X_X7) float64 0.008701\n",
      "  * X_X8     (X_X8) float64 0.02087\n",
      "[[[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  5.44747482e-01   5.11107711e-02   7.75277005e-02\n",
      "              1.59300939e-01   2.54185534e-02   1.49869532e-03\n",
      "              5.35822023e-03   5.19889092e-02   8.30487293e-02]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  2.33487631e-01   1.98865335e-06   4.44565558e-02\n",
      "              3.46186651e-01   6.23758072e-02   6.73565508e-03\n",
      "              1.62227144e-03   1.01811750e-02   2.94952264e-01]\n",
      "           [  1.00010211e-01   1.27030860e-03   2.47213449e-01\n",
      "              4.76712724e-01   3.99958761e-02   6.40957662e-05\n",
      "              5.92502536e-02   9.50530523e-03   6.59777766e-02]\n",
      "           [  1.06305319e-01   4.78529776e-03   8.71312188e-02\n",
      "              4.65190409e-01   1.44192214e-02   1.20718239e-04\n",
      "              1.13341776e-02   3.14251950e-02   2.79288443e-01]\n",
      "           [  3.05843814e-01   1.32689315e-05   4.00680731e-01\n",
      "              1.03029987e-01   7.90990919e-03   7.74854915e-03\n",
      "              8.61569512e-03   1.00345461e-01   6.58125860e-02]\n",
      "           [  5.80904825e-01   9.43147913e-07   1.09710141e-01\n",
      "              2.20955384e-01   3.39797375e-03   2.36307198e-02\n",
      "              3.45458874e-03   3.49580559e-02   2.29873677e-02]]]]]]]]]]]]\n",
      "[ 0.28171799  0.00076793  0.33020705  0.26546957  0.01433725  0.00981422\n",
      "  0.00870094  0.02086778  0.06811727]\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "NORMUNITY\n",
      "9 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>\n",
      "array([[[[[[[[[[ 18845.35210356]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.02335\n",
      "  * X_X2     (X_X2) float64 0.06201\n",
      "  * X_X3     (X_X3) float64 0.0399\n",
      "  * X_X4     (X_X4) float64 0.03426\n",
      "  * X_X5     (X_X5) float64 0.3044\n",
      "  * X_X6     (X_X6) float64 0.0172\n",
      "  * X_X7     (X_X7) float64 0.38\n",
      "  * X_X8     (X_X8) float64 0.1372\n",
      "[[[[[[[[[[[[  1.38470880e-01   3.75910811e-02   8.29710176e-02\n",
      "              7.85944496e-04   3.05692290e-01   1.05376609e-01\n",
      "              2.83576918e-01   4.54897921e-02   4.54674600e-05]\n",
      "           [  6.38642361e-03   7.49551392e-02   3.04119255e-01\n",
      "              3.69701112e-07   2.05764674e-01   1.56615833e-02\n",
      "              1.99463695e-01   1.82134667e-01   1.15141929e-02]\n",
      "           [  5.75022704e-02   4.62050137e-03   9.79584711e-02\n",
      "              5.08338298e-02   2.46979785e-01   2.49614713e-02\n",
      "              2.54960957e-01   2.62182413e-01   3.01692083e-07]\n",
      "           [  2.31628775e-03   3.64822978e-02   4.94493563e-03\n",
      "              2.13758097e-01   3.37032976e-01   4.51031308e-02\n",
      "              3.13586045e-01   3.99255322e-02   6.85069769e-03]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  2.45838153e-03   1.07952168e-02   4.25684339e-02\n",
      "              6.33233694e-02   2.10801668e-01   2.46123724e-01\n",
      "              2.79600566e-01   1.44328527e-01   1.13394214e-07]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e+00   1.00000000e-12   1.00000000e-12]\n",
      "           [  5.51169099e-09   2.24522430e-01   2.38209963e-03\n",
      "              3.09627411e-02   3.36091529e-01   1.24883447e-02\n",
      "              2.53679182e-01   1.38679511e-01   1.19415693e-03]\n",
      "           [  7.32878038e-06   8.58043598e-02   3.22480630e-02\n",
      "              6.46635289e-03   3.12293445e-01   1.29966531e-02\n",
      "              3.05240926e-01   5.56875890e-02   1.89255283e-01]]]]]]]]]]]]\n",
      "[ 0.02334898  0.06201314  0.03989883  0.03426188  0.3043589   0.01719865\n",
      "  0.37995435  0.13717176  0.00179351]\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>\n",
      "array([[[[[[[[[[ 17229.16218313]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.005729\n",
      "  * X_X2     (X_X2) float64 0.2758\n",
      "  * X_X3     (X_X3) float64 0.04265\n",
      "  * X_X4     (X_X4) float64 0.01342\n",
      "  * X_X5     (X_X5) float64 0.04879\n",
      "  * X_X6     (X_X6) float64 0.06658\n",
      "  * X_X7     (X_X7) float64 0.3424\n",
      "  * X_X8     (X_X8) float64 0.02674\n",
      "[[[[[[[[[[[[  8.69841068e-02   1.71778179e-01   1.33808650e-06\n",
      "              4.20061211e-03   1.75310765e-01   5.58322147e-02\n",
      "              2.87204061e-01   2.20819570e-03   2.16480527e-01]\n",
      "           [  1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  4.41070407e-03   2.59354668e-01   2.45190260e-08\n",
      "              6.93565693e-05   3.53542177e-02   1.80890893e-01\n",
      "              2.24746879e-01   1.21827671e-01   1.73345586e-01]\n",
      "           [  1.89479075e-02   2.05710293e-01   2.99295827e-05\n",
      "              6.36906645e-02   1.78370853e-01   1.96414509e-02\n",
      "              1.96446024e-01   6.09502329e-02   2.56212643e-01]\n",
      "           [  1.65397598e-05   2.46260170e-01   6.15121840e-02\n",
      "              4.40088784e-03   1.50386842e-01   1.39016888e-01\n",
      "              2.38657416e-01   1.59323769e-02   1.43816696e-01]\n",
      "           [  9.55388533e-08   1.41062047e-01   9.87435994e-02\n",
      "              1.04178342e-02   2.51380054e-02   5.16631863e-02\n",
      "              3.85285871e-01   8.26732641e-03   2.79422035e-01]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e+00   1.00000000e-12   1.00000000e-12]\n",
      "           [  4.80530670e-03   2.01155838e-01   2.15616077e-04\n",
      "              1.38203502e-01   3.89232302e-02   1.57996886e-01\n",
      "              2.22134688e-01   6.47996457e-02   1.71765287e-01]\n",
      "           [  8.43277284e-06   9.12344747e-02   3.77618278e-02\n",
      "              9.25555555e-03   1.53410813e-01   1.18325447e-01\n",
      "              2.52802517e-01   1.82089349e-02   3.18991997e-01]]]]]]]]]]]]\n",
      "[ 0.00572916  0.27580921  0.04265009  0.01342425  0.04878618  0.06658208\n",
      "  0.34240152  0.02674375  0.17787376]\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>\n",
      "array([[[[[[[[[[ 21237.55386294]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.2367\n",
      "  * X_X2     (X_X2) float64 0.02297\n",
      "  * X_X3     (X_X3) float64 0.006571\n",
      "  * X_X4     (X_X4) float64 0.3339\n",
      "  * X_X5     (X_X5) float64 0.08391\n",
      "  * X_X6     (X_X6) float64 0.00423\n",
      "  * X_X7     (X_X7) float64 0.05765\n",
      "  * X_X8     (X_X8) float64 0.1325\n",
      "[[[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  2.61280345e-01   2.50884548e-03   2.50971286e-05\n",
      "              1.22001364e-01   1.61454770e-02   5.21866265e-03\n",
      "              1.36395347e-01   2.41875877e-01   2.14548984e-01]\n",
      "           [  1.99043674e-01   8.46513687e-02   1.78488420e-02\n",
      "              2.10263033e-01   7.67601052e-02   1.02006476e-06\n",
      "              3.45780183e-02   1.75158191e-01   2.01695748e-01]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  1.10435544e-01   5.41510293e-02   3.42579749e-03\n",
      "              3.41523389e-03   2.97605171e-01   3.46558805e-05\n",
      "              8.09757987e-03   1.66758243e-01   3.56076746e-01]\n",
      "           [  2.00706994e-01   1.16943625e-02   1.79307717e-07\n",
      "              1.95186127e-01   1.92135919e-01   2.08540779e-02\n",
      "              8.84845103e-02   1.52427984e-01   1.38509845e-01]\n",
      "           [  2.58389865e-01   5.83934714e-02   7.50950020e-08\n",
      "              1.98174023e-01   1.17143369e-01   8.16756596e-03\n",
      "              1.34095877e-01   1.19114427e-01   1.06521327e-01]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e+00   1.00000000e-12]\n",
      "           [  3.18023778e-01   2.18592598e-02   3.02135104e-02\n",
      "              1.97055195e-01   1.46568197e-01   9.60574129e-07\n",
      "              5.23903405e-02   1.07938613e-01   1.25950146e-01]]]]]]]]]]]]\n",
      "[ 0.23665472  0.0229676   0.00657122  0.33394364  0.08390696  0.00422985\n",
      "  0.05765071  0.13248851  0.12158678]\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>\n",
      "array([[[[[[[[[[ 22714.1544021]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.04242\n",
      "  * X_X2     (X_X2) float64 0.007688\n",
      "  * X_X3     (X_X3) float64 0.002395\n",
      "  * X_X4     (X_X4) float64 0.03897\n",
      "  * X_X5     (X_X5) float64 0.05962\n",
      "  * X_X6     (X_X6) float64 0.2164\n",
      "  * X_X7     (X_X7) float64 0.3567\n",
      "  * X_X8     (X_X8) float64 0.08278\n",
      "[[[[[[[[[[[[  1.55860083e-01   2.24433685e-02   1.93698328e-04\n",
      "              5.17020701e-03   5.56081772e-02   2.14183690e-01\n",
      "              1.65223056e-01   8.52643689e-02   2.96053351e-01]\n",
      "           [  3.78208281e-02   1.61055479e-02   2.97884298e-04\n",
      "              8.11677031e-03   1.27880673e-01   2.98259518e-01\n",
      "              2.02556959e-01   1.37018104e-01   1.71943716e-01]\n",
      "           [  2.37544097e-02   2.32386898e-06   4.02441712e-02\n",
      "              7.89375554e-02   3.36242735e-02   3.55310292e-01\n",
      "              3.76335808e-01   5.51274875e-02   3.66636784e-02]\n",
      "           [  5.84899745e-03   7.85509743e-04   4.43242957e-03\n",
      "              2.17887089e-01   6.78627470e-02   2.52146182e-01\n",
      "              2.59422149e-01   5.64763332e-03   1.85967263e-01]\n",
      "           [  1.10867089e-02   7.09962104e-07   4.77729651e-03\n",
      "              9.38480165e-02   1.05196278e-01   3.71429744e-01\n",
      "              1.35005463e-01   2.01465267e-01   7.71905164e-02]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e+00   1.00000000e-12   1.00000000e-12]\n",
      "           [  5.16645352e-02   4.67254542e-02   6.44319867e-03\n",
      "              9.19397437e-02   1.26993232e-05   1.16450644e-01\n",
      "              2.39934723e-01   2.15500589e-01   2.31328413e-01]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e+00]]]]]]]]]]]]\n",
      "[ 0.04242133  0.00768829  0.00239485  0.03896605  0.05961641  0.2164378\n",
      "  0.35667422  0.08278424  0.19301679]\n",
      "TRIAL 14/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>\n",
      "array([[[[[[[[[[ 30379.06864119]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.00283\n",
      "  * X_X2     (X_X2) float64 0.158\n",
      "  * X_X3     (X_X3) float64 0.006395\n",
      "  * X_X4     (X_X4) float64 0.009869\n",
      "  * X_X5     (X_X5) float64 0.2794\n",
      "  * X_X6     (X_X6) float64 0.2289\n",
      "  * X_X7     (X_X7) float64 0.06501\n",
      "  * X_X8     (X_X8) float64 0.1835\n",
      "[[[[[[[[[[[[  3.61163653e-02   1.45635676e-01   8.76563458e-04\n",
      "              1.19819003e-03   2.76769221e-01   9.04968428e-02\n",
      "              1.05361720e-01   3.39366111e-01   4.17931040e-03]\n",
      "           [  1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  2.20017521e-07   2.36206741e-01   5.40237248e-02\n",
      "              9.61529998e-03   1.28720063e-01   1.73024791e-01\n",
      "              9.41448647e-02   2.50784268e-01   5.34800277e-02]\n",
      "           [  4.18321526e-03   1.78722848e-01   3.95401471e-06\n",
      "              7.10895770e-02   1.00289014e-01   2.56274323e-01\n",
      "              7.15545066e-02   2.17225868e-01   1.00656695e-01]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  9.86002829e-07   1.65909542e-01   9.66264355e-04\n",
      "              7.92454894e-04   1.34336649e-01   1.68737104e-01\n",
      "              2.09153510e-01   1.78762304e-01   1.41341186e-01]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e+00   1.00000000e-12]\n",
      "           [  2.69187535e-04   2.21965048e-01   1.22727952e-02\n",
      "              1.19947527e-02   1.84964679e-01   1.83535729e-01\n",
      "              5.10774946e-02   1.64619637e-01   1.69300678e-01]]]]]]]]]]]]\n",
      "[ 0.00282968  0.15800187  0.00639466  0.00986923  0.27939811  0.22891086\n",
      "  0.06500936  0.18352608  0.06606016]\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X1: 1, X_X2: 1, X_X3: 1, X_X4: 1, X_X5: 1, X_X6: 1, X_X7: 1, X_X8: 1)>\n",
      "array([[[[[[[[[[ 48780.99231019]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X1     (X_X1) float64 0.2817\n",
      "  * X_X2     (X_X2) float64 0.0007679\n",
      "  * X_X3     (X_X3) float64 0.3302\n",
      "  * X_X4     (X_X4) float64 0.2655\n",
      "  * X_X5     (X_X5) float64 0.01434\n",
      "  * X_X6     (X_X6) float64 0.009814\n",
      "  * X_X7     (X_X7) float64 0.008701\n",
      "  * X_X8     (X_X8) float64 0.02087\n",
      "[[[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  2.42921399e-01   6.83067510e-02   2.41270201e-01\n",
      "              1.50118177e-01   1.15174153e-02   3.33431800e-02\n",
      "              1.30536663e-02   1.77151406e-02   2.21754069e-01]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12]\n",
      "           [  3.51726388e-01   1.61249531e-05   6.10922172e-02\n",
      "              8.60933635e-02   9.07572905e-02   3.31300194e-03\n",
      "              4.17503580e-02   1.40716148e-02   3.51179641e-01]\n",
      "           [  1.28061724e-01   3.44103099e-04   1.43617125e-01\n",
      "              2.67312792e-01   4.79024662e-02   5.88413942e-02\n",
      "              9.94614714e-04   1.20403769e-01   2.32522011e-01]\n",
      "           [  2.33027583e-01   5.22866946e-05   1.19297223e-01\n",
      "              2.06664912e-01   9.61852251e-03   9.91220718e-02\n",
      "              7.12543029e-02   5.29279094e-02   2.08035189e-01]\n",
      "           [  2.34747950e-01   6.81532920e-04   2.52842441e-01\n",
      "              1.13472497e-01   8.64418639e-03   4.62828716e-02\n",
      "              1.19118499e-01   1.57041647e-01   6.71683747e-02]\n",
      "           [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "              1.00000000e-12   1.00000000e-12   1.00000000e+00]]]]]]]]]]]]\n",
      "[ 0.28171799  0.00076793  0.33020705  0.26546957  0.01433725  0.00981422\n",
      "  0.00870094  0.02086778  0.06811727]\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "PSEUDO\n",
      "10 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "QUASI\n",
      "10 COMPONENTS\n",
      "TRIAL 1/20\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "TRIAL 8/20\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "NORMUNITY\n",
      "10 COMPONENTS\n",
      "TRIAL 1/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X01: 1, X_X02: 1, X_X03: 1, X_X04: 1, X_X05: 1, X_X06: 1, X_X07: 1, X_X08: 1, X_X09: 1)>\n",
      "array([[[[[[[[[[[ 32992.27139106]]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X01    (X_X01) float64 0.06489\n",
      "  * X_X02    (X_X02) float64 0.02209\n",
      "  * X_X03    (X_X03) float64 0.08005\n",
      "  * X_X04    (X_X04) float64 0.02539\n",
      "  * X_X05    (X_X05) float64 0.002807\n",
      "  * X_X06    (X_X06) float64 0.2206\n",
      "  * X_X07    (X_X07) float64 0.3263\n",
      "  * X_X08    (X_X08) float64 0.1224\n",
      "  * X_X09    (X_X09) float64 0.005118\n",
      "[[[[[[[[[[[[[  9.28319273e-02   1.12761786e-01   2.15646733e-01\n",
      "               1.32299314e-02   7.36246031e-03   7.59389544e-02\n",
      "               1.91644900e-01   2.33460356e-01   1.86335353e-05\n",
      "               5.71043187e-02]\n",
      "            [  8.25166714e-02   2.76327247e-03   1.58166329e-01\n",
      "               2.21489913e-02   1.51720532e-04   2.36911395e-01\n",
      "               2.09283602e-01   2.25419923e-01   2.54921061e-02\n",
      "               3.71459893e-02]\n",
      "            [  1.94918877e-01   2.22781758e-03   2.38002395e-01\n",
      "               4.95664475e-03   7.31740448e-04   8.57676219e-03\n",
      "               1.84525940e-01   1.41064936e-01   2.58138390e-03\n",
      "               2.22413503e-01]\n",
      "            [  6.59858807e-02   1.40062627e-01   1.24390384e-01\n",
      "               9.23242504e-02   3.21058098e-05   8.41272960e-02\n",
      "               1.84474408e-01   1.58008662e-01   9.28513967e-04\n",
      "               1.49665873e-01]\n",
      "            [  2.48505716e-02   7.28459165e-02   2.08625197e-01\n",
      "               8.84015698e-03   3.67285303e-02   8.49977447e-02\n",
      "               2.54332356e-01   1.73203619e-01   8.09777923e-06\n",
      "               1.35567810e-01]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  8.94638922e-02   2.56407669e-02   4.34686382e-02\n",
      "               7.63292427e-02   9.67864954e-04   1.03662969e-01\n",
      "               2.43247000e-01   2.86652030e-01   1.34947787e-02\n",
      "               1.17072817e-01]\n",
      "            [  2.14050967e-01   9.25548789e-04   1.52001262e-02\n",
      "               6.15170509e-02   2.24791634e-05   2.21769992e-01\n",
      "               2.09324243e-01   6.76362252e-02   1.12105510e-02\n",
      "               1.98342816e-01]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e+00]]]]]]]]]]]]]\n",
      "[ 0.06488928  0.02209467  0.08005051  0.02538987  0.00280749  0.22061855\n",
      "  0.32629228  0.12243823  0.00511819  0.13030094]\n",
      "TRIAL 2/20\n",
      "TRIAL 3/20\n",
      "TRIAL 4/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X01: 1, X_X02: 1, X_X03: 1, X_X04: 1, X_X05: 1, X_X06: 1, X_X07: 1, X_X08: 1, X_X09: 1)>\n",
      "array([[[[[[[[[[[ 35906.97331219]]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X01    (X_X01) float64 0.2537\n",
      "  * X_X02    (X_X02) float64 0.2332\n",
      "  * X_X03    (X_X03) float64 0.1363\n",
      "  * X_X04    (X_X04) float64 0.04557\n",
      "  * X_X05    (X_X05) float64 0.01003\n",
      "  * X_X06    (X_X06) float64 0.1446\n",
      "  * X_X07    (X_X07) float64 0.009234\n",
      "  * X_X08    (X_X08) float64 0.1441\n",
      "  * X_X09    (X_X09) float64 0.01678\n",
      "[[[[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.66598660e-01   1.48727186e-01   1.98077174e-01\n",
      "               1.94907352e-01   2.08502521e-02   2.43630074e-01\n",
      "               3.22203084e-07   2.34995689e-02   9.68890280e-04\n",
      "               2.74052044e-03]\n",
      "            [  2.47661566e-01   1.00401392e-01   6.97649547e-02\n",
      "               7.15243160e-02   2.28845287e-06   2.48195992e-01\n",
      "               3.80687841e-02   2.11423076e-01   1.15366080e-02\n",
      "               1.42102276e-03]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.90025154e-01   7.34768819e-02   2.01056539e-01\n",
      "               9.85015056e-03   4.59820764e-02   1.79292945e-01\n",
      "               1.17500302e-06   2.35286548e-01   1.32448158e-02\n",
      "               5.17837139e-02]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e+00   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  2.05680919e-01   2.05086444e-01   1.44232066e-02\n",
      "               1.06670381e-01   2.17649849e-02   2.19287502e-01\n",
      "               1.34422111e-03   1.47097996e-01   7.86266993e-02\n",
      "               1.76457982e-05]\n",
      "            [  1.20587086e-01   7.79626369e-02   2.56469127e-01\n",
      "               6.42765056e-02   5.67070499e-02   2.31781352e-01\n",
      "               1.75491315e-03   1.38989660e-01   5.14469198e-02\n",
      "               2.47493634e-05]]]]]]]]]]]]]\n",
      "[ 0.253678    0.23320191  0.13628781  0.04556506  0.0100317   0.14460173\n",
      "  0.00923435  0.14410382  0.01677506  0.00652055]\n",
      "TRIAL 5/20\n",
      "TRIAL 6/20\n",
      "TRIAL 7/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X01: 1, X_X02: 1, X_X03: 1, X_X04: 1, X_X05: 1, X_X06: 1, X_X07: 1, X_X08: 1, X_X09: 1)>\n",
      "array([[[[[[[[[[[ 44334.2276809]]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X01    (X_X01) float64 0.5539\n",
      "  * X_X02    (X_X02) float64 0.001489\n",
      "  * X_X03    (X_X03) float64 0.006258\n",
      "  * X_X04    (X_X04) float64 0.01409\n",
      "  * X_X05    (X_X05) float64 0.02235\n",
      "  * X_X06    (X_X06) float64 0.1907\n",
      "  * X_X07    (X_X07) float64 0.1355\n",
      "  * X_X08    (X_X08) float64 0.03399\n",
      "  * X_X09    (X_X09) float64 0.004718\n",
      "[[[[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.29381233e-01   6.17042257e-04   1.44450715e-02\n",
      "               1.65719829e-02   9.31488358e-03   3.03914733e-01\n",
      "               2.15848090e-01   1.78035377e-01   4.04707504e-02\n",
      "               9.14008363e-02]\n",
      "            [  2.83992648e-01   5.50126245e-04   1.78171230e-02\n",
      "               2.29944091e-02   1.35887742e-01   1.87106553e-01\n",
      "               6.23213480e-02   2.15302755e-01   3.18445918e-02\n",
      "               4.21827038e-02]\n",
      "            [  1.26216164e-01   2.04161309e-02   2.50789842e-03\n",
      "               9.59769815e-02   8.78495894e-02   1.95847912e-01\n",
      "               1.08664660e-01   1.64018417e-01   1.35283133e-06\n",
      "               1.98500894e-01]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e+00]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  8.19885197e-02   1.99736376e-05   2.28176529e-03\n",
      "               9.07380841e-02   2.12067943e-01   1.55349900e-01\n",
      "               2.31126501e-01   2.65467904e-03   4.05685712e-02\n",
      "               1.83204063e-01]\n",
      "            [  1.04825050e-01   9.41733983e-03   6.49080302e-02\n",
      "               4.84212721e-02   6.99447259e-03   2.20980034e-01\n",
      "               1.74831367e-01   1.67767689e-01   1.17702312e-06\n",
      "               2.01853568e-01]\n",
      "            [  1.30633062e-01   7.99857971e-03   8.70160487e-03\n",
      "               1.44581187e-01   1.28871506e-01   1.17551674e-01\n",
      "               1.40014939e-01   1.58358571e-01   1.15513617e-07\n",
      "               1.63288763e-01]]]]]]]]]]]]]\n",
      "[ 0.55389931  0.00148919  0.00625847  0.01409256  0.02235222  0.19074092\n",
      "  0.13547748  0.03398572  0.00471823  0.03698591]\n",
      "TRIAL 8/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X01: 1, X_X02: 1, X_X03: 1, X_X04: 1, X_X05: 1, X_X06: 1, X_X07: 1, X_X08: 1, X_X09: 1)>\n",
      "array([[[[[[[[[[[ 25161.85918291]]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X01    (X_X01) float64 0.09572\n",
      "  * X_X02    (X_X02) float64 0.02795\n",
      "  * X_X03    (X_X03) float64 0.1509\n",
      "  * X_X04    (X_X04) float64 0.1646\n",
      "  * X_X05    (X_X05) float64 0.03749\n",
      "  * X_X06    (X_X06) float64 0.3075\n",
      "  * X_X07    (X_X07) float64 0.006549\n",
      "  * X_X08    (X_X08) float64 0.04803\n",
      "  * X_X09    (X_X09) float64 0.1564\n",
      "[[[[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.74857248e-01   2.96173558e-02   1.74856487e-01\n",
      "               1.77947578e-01   4.61011568e-02   1.96937302e-01\n",
      "               4.81331395e-09   1.42543234e-01   4.82827745e-02\n",
      "               8.85685981e-03]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.11826974e-01   2.28554216e-02   1.67839497e-01\n",
      "               1.84592309e-01   1.29616678e-01   1.38526357e-01\n",
      "               9.52754526e-03   5.08795955e-02   1.84335343e-01\n",
      "               2.80282346e-07]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.80131744e-02   6.27719523e-02   2.72634297e-01\n",
      "               1.31059257e-01   7.07206622e-02   2.81895541e-01\n",
      "               1.96623695e-02   8.80558396e-03   1.34411590e-01\n",
      "               2.55731459e-05]\n",
      "            [  1.06760933e-01   4.33499537e-03   1.68161042e-01\n",
      "               2.08222986e-01   5.79010217e-03   1.92143339e-01\n",
      "               1.79480472e-02   2.03455406e-01   9.27956882e-02\n",
      "               3.87461275e-04]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "               1.00000000e-12]\n",
      "            [  1.89372742e-01   2.07566903e-02   1.92335927e-01\n",
      "               2.14661048e-01   2.96712689e-02   1.80872445e-01\n",
      "               2.60699766e-06   9.22777721e-03   1.36420161e-01\n",
      "               2.66793332e-02]]]]]]]]]]]]]\n",
      "[ 0.09571528  0.0279492   0.15093519  0.16460131  0.03749186  0.30750794\n",
      "  0.0065485   0.04802951  0.15641141  0.00480979]\n",
      "TRIAL 9/20\n",
      "TRIAL 10/20\n",
      "TRIAL 11/20\n",
      "TRIAL 12/20\n",
      "TRIAL 13/20\n",
      "TRIAL 14/20\n",
      "TRIAL 15/20\n",
      "TRIAL 16/20\n",
      "<xarray.DataArray 'GM' (P: 1, T: 1, X_X01: 1, X_X02: 1, X_X03: 1, X_X04: 1, X_X05: 1, X_X06: 1, X_X07: 1, X_X08: 1, X_X09: 1)>\n",
      "array([[[[[[[[[[[ 15143.63158122]]]]]]]]]]])\n",
      "Coordinates:\n",
      "  * P        (P) float64 1.013e+05\n",
      "  * T        (T) float64 300.0\n",
      "  * X_X01    (X_X01) float64 0.1467\n",
      "  * X_X02    (X_X02) float64 0.08028\n",
      "  * X_X03    (X_X03) float64 0.05512\n",
      "  * X_X04    (X_X04) float64 0.06812\n",
      "  * X_X05    (X_X05) float64 0.004017\n",
      "  * X_X06    (X_X06) float64 0.07684\n",
      "  * X_X07    (X_X07) float64 0.1361\n",
      "  * X_X08    (X_X08) float64 0.01094\n",
      "  * X_X09    (X_X09) float64 0.3856\n",
      "[[[[[[[[[[[[[  1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.61789072e-01   1.45742581e-02   4.92647815e-03\n",
      "               2.08922862e-01   5.79231218e-05   9.30901223e-02\n",
      "               1.01954225e-01   1.03521159e-02   2.37300591e-01\n",
      "               1.67032353e-01]\n",
      "            [  1.82555591e-01   1.73741629e-01   1.39450608e-01\n",
      "               5.72234304e-02   8.75693696e-06   2.10634692e-02\n",
      "               1.24848504e-01   1.81913744e-03   2.77958720e-01\n",
      "               2.13301534e-02]\n",
      "            [  2.08639885e-01   1.33930760e-01   7.03464244e-03\n",
      "               1.46408215e-01   1.22733032e-04   9.48127269e-02\n",
      "               1.30640584e-01   1.36160838e-02   2.02274441e-01\n",
      "               6.25199282e-02]\n",
      "            [  9.19589496e-02   6.65839887e-02   1.51409610e-02\n",
      "               2.40449827e-02   5.20503261e-08   1.33658365e-01\n",
      "               2.71054421e-01   3.58197781e-02   2.10785919e-01\n",
      "               1.50952583e-01]\n",
      "            [  2.02930234e-01   8.94445084e-03   5.50923689e-02\n",
      "               1.23589063e-01   2.72600754e-05   2.04430013e-01\n",
      "               1.45437736e-01   4.47213236e-02   2.04248070e-01\n",
      "               1.05794809e-02]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e+00   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12]\n",
      "            [  1.85130462e-01   1.02335599e-01   3.25421169e-02\n",
      "               1.68911543e-02   4.35370755e-04   8.15001958e-02\n",
      "               2.03382202e-01   6.64794763e-02   2.67806634e-01\n",
      "               4.34967891e-02]\n",
      "            [  1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e-12\n",
      "               1.00000000e-12   1.00000000e-12   1.00000000e+00\n",
      "               1.00000000e-12]\n",
      "            [  2.03913693e-01   7.87279928e-02   1.82021907e-02\n",
      "               5.39568482e-02   2.70794604e-02   2.07898770e-01\n",
      "               2.12264418e-01   2.81805390e-04   1.91656541e-01\n",
      "               6.01828052e-03]]]]]]]]]]]]]\n",
      "[ 0.14668039  0.08028181  0.0551166   0.068121    0.00401697  0.07684048\n",
      "  0.13608608  0.01094241  0.38563009  0.03628417]\n",
      "TRIAL 17/20\n",
      "TRIAL 18/20\n",
      "TRIAL 19/20\n",
      "TRIAL 20/20\n",
      "CPU times: user 1h 11min 7s, sys: 1h 6min 39s, total: 2h 17min 47s\n",
      "Wall time: 1h 13min 4s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "comp_iters_pseudo = []\n",
    "comp_iters_quasi = []\n",
    "comp_iters_normunity = []\n",
    "max_comps = 10\n",
    "num_trials = 20\n",
    "for num_comps in range(2, max_comps+1):\n",
    "    solutions = [rng.dirichlet(np.ones(num_comps, dtype=np.int)) for i in range(num_trials)]\n",
    "    from pycalphad.core.lower_convex_hull import lower_convex_hull, _initialize_array\n",
    "    print('PSEUDO')\n",
    "    results_pseudo = test_opt(num_comps, num_trials, sample_pseudo, solutions)\n",
    "    print('QUASI')\n",
    "    results_quasi = test_opt(num_comps, num_trials, sample_quasi, solutions)\n",
    "    print('NORMUNITY')\n",
    "    results_normunity = test_opt(num_comps, num_trials, sample_baduniform, solutions)\n",
    "    comp_iters_pseudo.append(results_pseudo)\n",
    "    comp_iters_quasi.append(results_quasi)\n",
    "    comp_iters_normunity.append(results_normunity)\n",
    "comp_iters_pseudo = np.array(comp_iters_pseudo, dtype=np.float)\n",
    "comp_iters_quasi = np.array(comp_iters_quasi, dtype=np.float)\n",
    "comp_iters_normunity = np.array(comp_iters_normunity, dtype=np.float)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAESCAYAAAAVLtXjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXd4HNXZt++zq96bZcmyilXce8eSjbGNMc2hmQ5JwMmX\nN4SEBBwIeSEhhSRUEwJ5A4EAgZhmgjEBBxcwWDbuRbZsbFm9WL13ac/3x65WxZa9K205ss99XXtp\nZnZmzk+zu8+c+c0zzxFSSjQajUZzYWFwtwCNRqPRuB4d/DUajeYCRAd/jUajuQDRwV+j0WguQHTw\n12g0mgsQHfw1Go3mAsRlwV8IMVoIsV8Isc/yt1YI8WNXta/RaDSaboQ78vyFEAagEJgjpSxwuQCN\nRqO5wHGX7bMEOKkDv0aj0bgHdwX/m4A1bmpbo9FoLnhcbvsIITyBYmC8lLLcpY1rNBqNBgAPN7R5\nObC3v8C/fPly2dLSQlRUFAD+/v4kJyczdepUAA4cOADg1PmsrCxuuOEGl7XX33zXtLvaB3j//fdd\nfvz7zuvPQ63Po+cxuNA/D3d9PwEOHjzIqVOnAEhKSuKvf/2rwA7c0fNfA2yQUr5+pvfvvPNO+dxz\nz7lUU1/++Mc/8tBDD7lVgyo6tAa1dKigQRUdKmhQRcdPfvIT3njjDbuCv0s9fyGEH+abvR/0t07X\nmcyd5Ofnu1sCoIYOraEbFXSooAHU0KGCBlBHh7241PaRUjYBw1zZpkaj0WhOx/jrX//a3Rp6UVZW\n9utp06a5VUNwcDBxcXFu1aCKDq1BLR0qaFBFhwoaVNFRUlLCvHnzHrNnG7c85HU2Nm/eLKdPn+5u\nGRqNRjNk2LdvH4sXL7bL83dHts9ZOXDgAGcK/lJKysrK6OzsdLqG2tpagoODnd7OUNBhrwaj0Uhk\nZCRC2PU9PCvbtm0jLS3NYfsbyjpU0KCKDhU0qKTDXpQL/v1RVlZGYGAgfn5+Tm9rxIgRTm/DFlTQ\nYa+GpqYmysrKGD58uJMUaTQaRzBkbJ/i4mIlgqHm3OjPSqNxLeeF7aPRaDT9kZ9dSUF2FTu2nATg\nokVJAMQmhhGXGO5OaUMO5er593yC7VwcLK7njb0lLP37fpb+fT9v7C3hjb0lHCyud+k+NM5l27Zt\n7pYAqKFDBQ3gPh1xieGkLkkBIK8ok9QlKaQuSXFr4FflM7GXId3znzIikCkjAnlzv/nBsDtnRLtl\nH+6koKCAqVOnUl5ejsGg3Llco9EoinLRoquGxVBiypQpxMTEEBcXx7hx47jnnntoampyWfuOzKxR\nEVUyKVTQoYIGUENHfMx4d0sA1DgWA0G54D8UEULw9ttvk5+fz+eff86BAwd4+umn3S1Lo9Fo+kW5\n4G+P568SXVlTUVFRLFmyhKNHj7JmzRqmT59OXFwc06dPZ+3atdb133zzTebOnUtSUhIrVqygsLAQ\nMNs44eHhmEwm67rLly/nzTffBMBkMvHII4+QkpLCjBkz+Oyzz3rpOHXqFLfddhtJSUnMmjWLN954\nw9n/utNRxVNVQYcKGkANHXlFme6WAKhxLAbCkPb8+7L07/vdLYHCwkI2btzI0qVLeeihh/j8889J\nTEykrKyM6upqAD755BOee+451qxZQ2JiIqtXr2blypVs2LABOLuN8/rrr7Nx40a+/PJL/Pz8uPPO\nO3u9f/fddzNx4kSOHTvGN998w3XXXUdiYuKQvTTVaDTOQbme/1D0/AFuv/12EhMTueqqq5g/fz4/\n//nPMRqNZGZm0tLSQmRkJGPGjAHgtdde47777iM5ORmDwcB9993H4cOHrb3/s7Fu3Tp+8IMfEB0d\nTXBwMPfdd5/1vcLCQnbv3s2vfvUrPD09mThxInfccQdvv/220/5vV6DKiUsFHSpoADV0aM9/cJxX\nPf/PVg68INxgrxreeust5s+f32vZK6+8wvPPP8+9997L3Llz+e1vf0tycjIFBQX84he/4JFHHgHM\nlpEQgpKSEusgNv1RUlJCTEyMdT42NtY6XVpaSmhoaK+noGNjY4eslabRaJyHcj3/oRqozvSk9CWX\nXMIHH3zAsWPHSE5OtvbSY2JiePbZZ8nOziY7O5ucnBwKCgqYNWuWNXD3zBYqLS21TkdFRVFUVGSd\nLygo6PVedXU1jY2N1mWFhYVERw+t9NW+qOKpqqBDBQ2ghg7t+Q8O5YL/+UJ5eTmffvopTU1NeHp6\n4u/vb83D/+53v8szzzzDsWPHAKirq2PdunUAhIeHEx0dzXvvvYfJZOLNN98kNzfXut9rrrmGl156\nieLiYmpqavjzn/9sfS8mJobZs2fz29/+ltbWVo4cOcKbb77JTTfd5Lp/XKPRDAmUC/5D0fM/0w1a\nk8nEiy++yIQJE0hOTmbHjh089dRTAFx55ZXcd999rFy5koSEBNLS0ti8ebN129WrV/PnP/+Z5ORk\njh8/zpw5c6zv3XnnnSxatIgFCxawaNEirr766l7tvvzyy+Tl5TF+/Hi+/e1v84tf/OI0O2qooYqn\nqoIOFTSAGjq05z84zovCbl1+vSM8/8HsQ2NGF3bTOJunHjZnxj3w+DI3K1GDgRR2U67nP5DaPrdP\ni+L2aVGDqu0zmH1onIsqnqoKOlTQAGro0J7/4BjS2T5ddXncvQ+NRqMZaijX8x+Knr/Guajiqaqg\nQwUNoIYO7fkPDuWCv0aj0Wicj3LBf6jm+Wuchyqeqgo6VNAAaujQnv/gcGnwF0IECyHeE0IcFUIc\nEULMOfdWGo1Go3E0rr7h+xzwiZRyhRDCAzhtNHbt+Wv6ooqnqoIOFTSAGjq05z84XBb8hRBBwHwp\n5XcApJQdQN1g9lmZvo+q7fs4+fSrACTdfxcAYfOmE556+rMCztqHRqPRDDVcafuMAiqEEP8QQuwT\nQrwkhPDtu5I9nn946nRSVq20zqesWknKqpV2BW1H7MOdxMXFkZ+f724ZTkUVT1UFHSpoADV0aM9/\ncLjS9vEApgP3SCn3CCFWAw8Bv+q50tatW9mzZw9xcXEABAcHM2nSJBITE10odWD861//4sUXXyQ3\nN5egoCCuvPJKHnnkEYKCgpzWpoqBv7a21vqEb9cPo+vSeCDzGRkZg9r+fJrPyMhQQk8Xrm7/9Q8/\nI6uqmfLQMMrqfHjkVXNNrBuWLWLKiEC3HA93fD+7prt+/zNnzmTx4sXYg8vKOwghhgM7pJSJlvk0\n4EEpZa/iNAMp77Ahah4Ay05tH7C+we7jL3/5Cy+88AIvvvgiCxYsoKSkhPvvv5+qqio2bNiA0Wgc\nsLahhi7voHE2urxDb5Qu7yClLAUKhBCjLYsWA2pctw2S+vp6nnjiCf70pz9xySWXYDQaGTlyJK++\n+ir5+fm899573HPPPTz++OPWbdLT05k4caJ1/rnnnmPGjBnExcUxb948/vOf/1jfy8nJ4eqrryYh\nIYHRo0ezcmW3TRUeHt6r6qdGo9HYgquzfX4MvCWE8ASyge/2XeHAgQOcqedvC129d1eza9cuWltb\nueqqq3ot9/f3Z8mSJXz++ed4enqetl3PaqCjRo3i008/JTIykg8//JAf/OAH7N27l8jISB5//HEW\nLVrE+vXraWtrY//+/Wfcx/nKtm3blMioUEGHChpU0WH2/N3f81fhWAwElwZ/KeVBYJYr23QFVVVV\nhIeHW+v192T48OEcOnTonCN0LV++3Dp9zTXX8Oyzz7Jv3z6WLVuGp6cnBQUFVjulZ4ln1aqyajQX\nAvnZlRRkV7Fjy0nyijIRLcMBiE0MIy4x3M3qbEO5wm6DyfN3hOc/EMLCwqisrMRkMp12AigtLSU8\n/Nxfhrfffpu//vWv1hs4TU1NVFZWAvDYY4/x+9//nksvvZSQkBB++MMfcttttw1Y71BDlV6VCjpU\n0ABq6HBnnn9cYjhxieHs2HKS+JjxpC5JcZuWgaJceYehyKxZs/D29mb9+vW9ljc0NLBp0yZSU1Px\n8/OjubnZ+t6pU6es04WFhfz0pz/lySefJCcnh5ycHMaOHWvt1Q8bNozVq1dz5MgRnn76aVatWqV9\nfo1GMyiUC/5DsbZPUFAQDzzwAA899BCbN2+mo6OD/Px87r77biIiIlixYgUTJ05k48aN1NTUUFpa\nyt/+9jfr9o2NjRgMBsLDwzGZTLz11lscPXrU+v66desoLi4GzKmvBoPhjBbT+YoqedQq6FBBA6ih\nQ5U8f1V02Ityts9Q5cc//jHh4eE8+uij5OTk0NraSlpaGh988AG+vr7cdNNNbN26lSlTphAfH8+t\nt97KCy+8AMCYMWP44Q9/yNKlSzEajdx0003MnTvXuu/9+/fz8MMPU19fT2RkJH/4wx+sz0FcCDd8\nNRqN4zkvhnFUIc+/L2vWrOEPf/gDn376KTExMQ7Z51BB5/lrnI278/wPFtdzsKSB7VuyAJi3KBmA\nKdEBbhkcaiB5/kM6+HfV5enLQGr7DGYf/fHee+/h4eHBtddeO6j9DDV08Nc4i67f647N5qB70WJz\n0HVXLS53n4S6GEjwV872sSfPPzx18B+4I/bRHytWrHDKfi80VMmjVkGHChrcqaPr97quegN5RZnc\n2aMul7tQ5XkDe7lw7hpqNBqNxopywV/X89f0RYWeLqihQwUNoIYOVer5q6LDXpQL/hqNRqNxPsoF\n/6GY569xLirklIMaOlTQAGroUCW/XhUd9qJc8NdoNBqN81Eu28cez79ncSWAixYlAfYVV3LEPjTO\nRQV/GdTQoYIGUENHfMx4OjtNGI3u7cNqz98NxCWG9yqolLokhdQlKXYFbUfsY6hz44038s4777hb\nhkZjEx3tndbpV57+kv078mjvsUxjGwMK/kKIS4QQFztaDAxNz3/KlCmMGTOmV+G2f/7zn73KNDuL\ngoICa02gnvQdPOZsvPvuu9x0002A+cnkK664wuE6B4MK/jKooUMFDeA+HZ2dJta/fRAwe+11NS1s\nXn+Ul57YytdfnKSlud3lms5rz18IsVUIkWqZfhB4G/iXEOJhZ4obKgghMJlM/N///d9py+1lIE9c\nO7K+j5RS1wvSKIk0STaszeDk0TLrsuW3TmV4TBDNjW1s++wELz2xlS83fENjfasblQ4NbO35TwS+\ntkx/D7gEmAv8wNGChmqe/7333ssLL7xAXV3dae/t3LmTJUuWMGrUKJYsWcKuXbus7y1fvpzf//73\nXH755YwcOZK8vDzrsmXLlhEXF8dtt91GdXU1/+///T/i4+NZsmQJhYWFNmvr6s0/+uijJCYmMn36\ndDZt2tRLw5tvvsnx48d54IEH2L17N3FxcSQmJrJ///5e5aUB1q9fz4IFCwZ4pOxHBX8Z1NChggZw\nvQ4pJZs+yuTogRI8vczjYcfHjGf0xChu/+FF3PDdmcQlhtHW2sGuL3N46cmtbPzwCDVVTU7XNlQ9\nf1tv+BoAKYRIwlwPKBNACBHqNGUDoKvOhjuYOnUqqampPP/88/zyl7+0Lq+pqeGWW27hiSee4Lrr\nruPf//43N998M/v27SMkJAQw2y7vvfceycnJVvvmww8/ZO3atYSGhrJ06VKWLVvGU089xV//+ld+\n9KMf8ac//Ynnn3/eZn379u3j1ltv5eTJk7z22mv85Cc/4ciRI73WGT16NE8//TRvvvlmrzGEw8LC\n2LJlC4sXLwbMNYtuueWWAR8rjcZevvrvcQ7uKsDoYeDaO6fz7t93W98TQpCQEkFCSgQlBTXs3JpN\nVmYZB3cVcGh3AWMnRzP74kSGRTmu4FpXjaHIfeYaQyeeNHfG3FVjaCDY2vPfBvwFeAr4N4DlRFDh\naEFD0fPv4qGHHuLvf/87VVVV1mWfffYZSUlJ3HDDDRgMBq6//npSUlLYsKH7RHXLLbcwevRoDAYD\nHh7m8/Gtt95KXFwcgYGBLFmyhISEBObPn4/BYOBb3/oWGRkZdmmLjY3l9ttvRwjBzTffzKlTpygv\nL7dp25tvvpl3330XgOrqarZs2cL1119vV/uD4UL3uVXTAK7VsfOLk+z6MgeDQbD81qnWZIwzee3R\nsSFcc/t0vvOTNMZPG4EQgqMHS3j9z+l88PpeivKqHaIpPHU6KatWUjZ9IbuHR5KyaiUpq1YOmcAP\ntvf8vwPcD5QDT1qWjQWec4KmATOYynqOuGoYN24cl112Gc8++yyjR49GSsmpU6eIjY3ttV5sbCwl\nJSXW+TOVfB42bJh12sfH57T5xsZGAIxG8yVwe3s73t7e1nU6Ojp6DRofGRlpnfb19QXMg8j03G9/\nrFixgmeeeYbm5mY+/PBDLrrool770ziX82G82IGyf0ceX312AgRcsWIySWNt+95FDA/gihWTSV2S\nwp5tOWTsKST7m3KyvylnZEIosy9OZNToiAv6/pZNPX8pZaWU8mEp5a+klA2WZf+RUq52tKCh6vl3\n8eCDD/LGG29QUlKCEILo6GjruLxdFBYWEh0dbZ0fzBcwKioKT0/P09rIy8s77aRjC2fSEh0dzaxZ\ns1i/fn2vzCBXcaH63F30TEfuGi/W3enIrjgWR/YXsXm9eUS7S781gbFTonu9b4vXHhzqy+Krx/P9\nVQuZe0kS3j4eFOZW88Hre3njL9s5drAEU6fpnPs5G+e15y+E8MLc+58KBPR8T0p5p+NlDV1GjRrF\ntddey0svvcSECRNYsmQJDz74IGvXruWaa65h3bp1HD9+nGXLHFMC1mAwcPXVV/O73/2O1atXExgY\naG1jyZIldu9v2LBhFBcX097e3uvK4aabbmL16tUUFhZy1VVXOUS7ZmjhygciT2SWsmHtYQAWLBvD\nlNn2d2R64hfgRdqlKcyaP4qDuwrYm55LeUk9H79zkJCNfsyan8CE6TF4eBodIX9IYKvn/zpwH1AP\nnOzzcihD0fPv21tetWoVzc3NCCEIDQ1lzZo1vPDCCyQnJ/PCCy/w9ttvW2/2nqmnbe+VwFNPPUVo\naChpaWmMGTOGV155hXfeeYeIiAibNPecXrBgAWPHjmXs2LGMHj3auvzKK6+koKCAq666Ch8fH7v0\nDZYL0efuD3fmlPe8AskrynTaFUheVgUfrzmANEnmLkxk9oJRZ15vAMfC28eD2QtG8b0HFnDpNRMI\nCfOjpqqJjesyefmpL9n1ZQ5trR326R2ief62ev7LgFFSyprBNCaEyAVqARPQLqWcPZj9qcL+/ft7\nzcfExFBUVGSdnzNnDlu2bDnjtuvWrTvnsp7ZQwAXX3wxu3d3ZzsEBQWxenX/Dtwtt9xyWnZORUX3\nvfqe7Xl6erJmzZrT9uHr60tERAQ33nhjv+1oNIOlOL+af/9zP52dkmkXxZF6acq5NxoAHp5GpsyO\nZdKMGI4fLmXnl9mUl9Tz5YZv2PnFSaZdFM/0i+LxC/BySvsqYGvwzwe8z7nWuTEBC6WU/d5yH0ht\nn67Lz/RNJ4CB1fYZzD4uBD766CMMBoNL8/u7uNA9/56401/umd4YCZx48u+A49Iby0rqWPvaXjra\nO5kwfQSLrhx31qtgRxwLg9HA2CnRjJkcRc7xCnZtzaYwt5qvPz/Jnm05TJo5klnzRxEU4utUHe7A\n1uD/BrBOCPEcUNrzDSnlmbu0Z0bgwHpCcYnhgw7QjtjH+c7y5cs5fvz4aU8way4seg6hCJCyynFD\nF1ZVNPL+q3tobekgZcJwLrt2IsLgukwcIQSJY4aROGYYRXnV7NyaTfaxcvbvyOfgzgLGTolm9oJE\nIoYHnHtnQwRbg/+PLH/7FouRQKId7UlgoxCiE3hJSvly3xXsGcNX4xo++ugjt7Z/oY9b2xNVxot1\npI66mmbee2U3TY1tJKSEc+VNUzDYUKnTWcciJj6U6+6cQXlJPbu+zObYoRIy9xeTub+Y5PGRzLk4\nkejYEKfrcDY2BX8p5ZnvuNhPqpSyRAgxDPNJ4KiUstddtK1bt7Jnzx7i4uIACA4OZtKkSSQm2nOO\n0biT2tpaRowYAXTfJO0KmgOZz8jIGNT2Q33+ZEUTInYSWaH+fHOimEdeXUfKlNlMiQ6gPvugy/Qc\nLK7n/Q1bOFJXBIE+vLG3hBMHd5Ec5su3r1k6oP1v+uxztnycSVhAEiPiQghPaOLrr7efc/sunPn/\nDosOJDimngmBHhjbojm8r4jNG79g88YvmD8/jdkXJ5JXlElpRa5L9PT9/7dt22ZN8Z45c6b1CXxb\nEbYWEhNCeADzgBigENghpbTvtnjv/f0KqJdSPtNz+ebNm+WZev7FxcXWgKJRG/1ZOYeuBxEH8zCj\nSjpamtt55+VdlJ+qJzI6kBtXzsbH1/PcGzpQgz001reyNz2XAzvzaWvtXUL6Z7+7DIMLbaq+7Nu3\nj8WLF9slwNY8/7HAesAXKABigRYhxNVSyqM27sMPMEgpG4QQ/sBS4DFbhRqNRpqamvDz87N1E40b\naGpqsj51rNH0R1trB2tf20P5qXrCIvy5/rszbQ787sI/0JsFy8Yw++JEDu7MZ096Hs2NbQC0t3Xg\n7aO2/r7Y6vm/CLwEPCUtlwpCiAcsyy+xcR/DgX8LIaSl3beklJ/1Xak/zz8yMpKysjJqagaVbWoT\ntbW1BAcHO72doaDDXg1Go9HhpR9U8NpV0aGKvzwYHR3tnXz45n5KCmoJDPHhhrtm4h9gfzKhu46F\nj68ncxYmMT01ged+tZG8oky8fdz/mdiLrcF/KnCp7O0RrQZ+2c/6pyGlzLHsZ0AIIRg+fPhAN7eL\n7Oxsxo0b55K2VNehggbN+UNnp4mP3z5I/slK/AK8uPGuWWdNo1QZzyH+NLCtaZfFQN+Ru+ZbljsU\nFWr7uLt314UKOrSGblTQoUpO+UB0SJPkv2sPk3W0DB9fT1bcNYvQCH+XanAGquiwF1t7/g8DHwkh\nPgbygHjgSuB2ZwnTaDTnD1JKNq3PJPNAMZ5eRq7/zgyH1tfX2I+tVT0/AqYDh4FAy98ZUsrTaxMM\nEhVq+6hQwwXU0KE1dKOCDlXqyNirY9tnJzi40zIYyx3Te+XJu0qDs1BFh73Y2vNHSnkc+J0TtWg0\nmvOQnVuz2bk1G2EQXH3LVOKS9BP1KtBv8BdCvCSl/L5l+p+Yn849DUeXdNaefzcq6NAaulFBhyr+\nsq06Dnydz1f/PW4ejOWGSSSPc1wm2FA7Fqpxtp5/To/pLGcL0Wg05xeZ+4vZtN5siVy6fDzjpuoH\n/1SiX89fSvmHHrN/k1I+1vcF/M3RgrTn340KOrSGblTQoYq/fC4dWZmlfLo2AyQsWDaaKXPiXK7B\nVaiiw15sTfU83s/yoflfazQap5GXVcl6y2AscxYmMnuBrsulIrbe8D2tZoQQIghzfX6Hoj3/blTQ\ncaFr6Dt0oWhx/XgPzq6jPxD687mL86v58M19dHZKps6NI81Jg7GcTYOrUUWHvZw1+AshCjDf6PUV\nQuT3eTscOH3IJ43mPKJrvIeu4N81jKErCU+dTti8aayr/i8A9aljiRoZTFB0kMu1nI3yknrWvraX\n9rZOxk8bweKrzj4Yi8a9nKvnfzvmXv8nwB09lkugVEr5jaMFqVDPX4UaLqro0Bq6cUctGSklBdlV\npG/KsmrY8rH5PWEQRAwPICommKiR5lfE8ACMNtTCHyx9j0V1RSPv/WM3rS0dJI+PZNl1zhuMpaSu\nldxgP7JLjuK1rYAwPw9C/TwJ8/U0T/t6EurrgacLjgOoU2/JXs4a/KWUWwGEEBFSyibXSNJoNGC2\nnLZvyqIwt/eopxNnxHCqqJbK0gbKS+opL6knY08hAB4eBoZFB5pPBpaTQliEv1NHxaqraebdV3fT\n1NBGfHI4V9081abBWGxFSklOVQvpeTWk59aQXdUC4YHU1Xjx8bGKfrcL8jaedlII8/MkzNfD8teT\nUD8PAryMF+QViq2DuTQJIaZirucTQY97AFLKRx0pSHv+3aigQ2voxlXebkFOFds3ZVGQUwWYq0jO\nTEtg28YTxMeMZ9n1kwBoa+ugrLieU4W1lBbVcqqwlurKJkoKaikpqLXuz8vbyPAR3VcHUSODCArx\nHVTA6zoWjfWtvPfqbuprWhgRF8K3bp+Gh8fgA79JSo6WNZKeW0t6bg0l9W3W9/w8DYSU1zKrM5DZ\nKb5UeXhT3WmkqrmdqqZ2qps7qG5up661k7rWTvKqW87alqdRmE8EPU4K9lxNnJeefxdCiO8DzwKf\nAZcDn2Kux+/w8g4azYVKYW412zedID/bHPS9fTyYmZbA9HnxePt4sm3jiV7re3l5MDIhlJEJodZl\nLc3t1hPBqcI6ThXVUl/bQkFOlfVkAuDr58lwy9VBtOWk4B9oX1nlluZ23n9tD9UVTQyLDuS6b8/A\ny8vmogGn0d5p4mBJA+m5NezIq6Wq2TxWlLGjnaSqUmY3l5FQVoRn1knqM08ipAnegGggxscLz+Ag\nPIMD8QwNwhgUAIEBdPj70+bnT7OPLw0+vtR7+lLj4UOlhw8VBi9K8aKhE0ob2ihtaDu7QHpeTZhP\nFFlhAfh0dJ5zOxWx9ZP6ObBMSvmVEKJaSnmtEOJy4GZHC9Kev1o6tIZunOXtFuVVs31zFnlZlYA5\n6M9INQf9vgOcnEuDj68n8ckRxCdHWJc11reaTwbWk0ItzU3t5B6vIPd4t20SGOxjsYqCGG6xjPob\nYCWvKJMPXg+hvKSe0Ag/bhjgYCzN7Z3sLaxnW24NOwvqaGlqIaK0mJiifNLKiogtK8Q7vwA6zCeC\nVssL4KC3YFZgCO219Zha2mhtqaC1tH8byN/yiuqz3OjviyEoEAID6AwIoN3Pj1ZfP5p8/Gjw9qXO\nw4cqT2+qjD40+/hR4+tHqa8frT6+yBB/Wo7v4/Six+pja/CPlFJ+ZZk2CSEMUspPhRBvOUuYRnO+\nU1JQQ/qmLHJPmAOWl7eR6fMSmJmW4NBRrfwDvUkaF0mSpbSClJK6mmbzlYHlZFBabL5CqK9t4URm\nqXXbkHC/XjeUI0d0V+Iszq8hMMSHFXfNsmswlrqWDr7Or2V7VgX5e48RWpDH8OICrinKZ1hZMYbO\nPj1pIfBPSSB4yhiCJo8leMpY/vF+HoVlWTz4+s+QUtLZ1EJ7TR0dtfW019TTXlNn/ltbR3ttPe3V\n5r8dPabg98lsAAAgAElEQVTbLet2NjbT2dgMJWUAGAE/yyuir/g+utq8vDnoDTxxt83/vyrYGvwL\nhRAJUspczA98fUsIUQGc+zrJTrTn340KOrSGbhzl7ZYU1LB9cxY5ll63p5eRGfPimZGWgK+fl9M1\nCCEIDvUjONSPMZPM/WBpklRVNPa6QigrqaemsomayiaOHSqxbNutw8/fixU2DsZSVtXAzi8zOLk9\ng/bME0QW5zOrtJi5Zwz08QRNHkPwlHEETR5D0MQUPAJ61/03rSuxHgshBB7+vnj4+0KMfQM+SZOJ\njoYm6wnDfPKo63Hy6D6RdNTW91reUdeAV2sLUzwHX6HUHdga/J8AxgG5wG+A9wEv4MfOkaXRnH+c\nKqxl++Yssr8pB8xBf/pF8cycf+6g72yEQRAeGUB4ZAATpscA5lG3KkobrFcHp4pqqShtAMuAfivu\nmkXYGQZjMbW1U3/0JHk7D5Pz9WFaDn9DQFEhnp2djO2xnhQCn8RYQqeNI3jyWHOgnzT6tEDvTITB\ngGdQAJ5BARAXbde2po4Onvv5Rxg62p2kzrmcM/gLc0rAl0A+gMXuCQW8pJQNjhakPX+1dGgN3QzU\n8y8tMgf9k8e6g/60uXHMnD8KP3/7gr4rc8qNRgPDRwQxfEQQU2bHAtDe1slzvzaPWzssepk50B/L\npu7QMWoPHKN831FajmcjLB69t+UlhaA5Ohq/iaOJnzORyBnjBx3o3Z1fb/DwoNPHj+whWtvnnMFf\nSimFEBmYB3HpWtaGEywfjeZ8oqy4ju2bs8g6avaSPTzNQX/W/FH4Bbi3pz9QZF0dIScO0pC9l+2X\n3UX90ZPItt49XwFURURSOTIev0mjSb5oMtMXTiEwTK0nki90bLV99gOjgWNO1AJoz78nKujQGrqx\n1W8vL6ln++Ys641TD08DU+fEMWvBKLtujA5Gg6OQUtJwLJuyjemUf7aNmr1HGCklI4G6IvM6VRGR\nlI2Io3RELE2JiaTMnchF46OZMiIQDyc+XKZKfr0qOuzF1uD/BbBBCPEa0FXvBwAp5auOl6XRDD3K\nT1mC/hFL0PcwMGVOLLMXJNqdQ+8umts7Ka9povirvdRu3k5H+m6Mp7qzfzqNHhQkppCXNJbSEXGU\njYglIjKY1PgQ7kgIYWykH4YL8GnZoYitwT8V8+AufZNZJeDQ4K89f7V0aA3d9OcxV5TWs33zSY4f\nPgWA0cPAlFmxzL54FAFBPi7RcC6klDS2dVLe2E5FYzsVTe1UNLZR0dhOeWMbdaVVBOw7wIjDB4nP\nOop3q/mpWCPQ5B9A9uiJZI+dRF7yWDq8vGn/Zi933bCA1IQQEkJ9XFYeoavK6rD9W8mqKyZ9UxLg\n2iqrfXH3vYeBYmt5h0sc1aAQwgDsAQqllMsdtV+NxtVUljWwY0sWxzJOgQSjUTB5ViyzL04kMNix\nQf9smKSktqWDysZ2S3C3BPWmdiob26wBv6WjRwV2KQkrP0XSsQySjmUQXZCDQXaP1FofM5LGmdMR\n82YTNHUcc4N8ucrfk2H+Xrz6+83kl9dz+3T7smMcQVeV1frbt1JpanRLldXzBZufxRZChANXAFFS\nyieFECMAg5Sy0M42f4J5EJgz3v3Rnn83KujQGrrp8naryhvYseUkRw+VWIP+pJmxzFno3KDfZhCI\ncTN4aWcRlU3mHntFYzuVje20m844xHYv/IRkXEkOiccyiMw4gHdpWfebHh4Ez51K9LI0Ii9Nwy++\n/yEXBWr43OMNrksJPRsqHIuBYGttn4uBtZh77KnAk0AK8ABwta2NCSFGYj6B/B74mb1iNRp388m7\nhzh6sBgpwWAUTJoxkjkLE2160GmgSCnZeKKK9NgI2o0GjmWUnbZOoLeRYf6eRPh7Ee7naZ0Oa2/G\nZ+8+2r7aRc3WnXTUN1q38QwLZtjieUQuTSVi4Rw8AtUIphrXYGvPfzVwk5RysxCiq77sTmC2ne09\nC6wCgvtbQXv+aum40DVUpu+jLH0/PhWCb1rNX32BZHScLxffNIfgUOcFfYCi2laeS8/nQHEDGA3I\nzD3c/e2rifDzsgR4T8L9vfCxVNKUUtKYlUf5Z1so27iNql0ZYOq2ewJGj2LY0lQil6YRMmMCwmgc\nkC4VfO5MU6MSTrsKx2Ig2Br8E6SUmy3TXdeXbXZsjxDiSswDwBwQQizkDENDAmzdupU9e/YQF2ce\n8Dk4OJhJkyZZf/xdg2g7cz4jI8Ol7ak8n5GR4XY97vw8dlRWsL2wlbCIRCiqwRhQzvhp0Vx2+UKn\ntj93XirvHSrlxfc20G6SxE6YSWx+JQXfZDCyfhxpU7rXP9nRyQSPAMo2bmPLh+tpLSm3WiJHRQuB\nE5NZctP1RC5NY19BDuXAmNmTB6yv54Dl7vp+urv9rvm8okxKK3JdrqdrOj/fPMDizJkzWbx4MfYg\npDy3VyiESAd+I6X8rxCiSkoZJoRYCjwspVxoU0NCPI55ZLAOwBfzQ2MfSCnv7Lne5s2bpbt7/hqN\nlJLDe4vY/FEmHT1ulD7wuPN7eEfLGnn2q3xyLXXoL00J4/tzYnj5N5usGtpr6ijf8jVln22j4vOd\ndNTWW7f3DA1i2OKLGHZpGhGXzDGXLnAgTz28warDXWyImgfAslPb3dJ+V9ZR1tPmZMfk++8C3Jd1\ntG/fPhYvXmxXypWtPff7gY+FEP/BPJ7v3zB7/d+ytSEp5cPAw2C9h3B/38Cv0ahAW2sHmz7KJHN/\nMWAeOevw3iKnt9vY1sk/9hSzPrMCCYwI8uInqXFMizE/XO/R1EBwdgY7r11Pza5DyB5F0fxT4om8\nNI1hS1MJmTkRg8fA6+przk3PrCOA1CV/cLMi+7E11fNrIcRkzD33VzE/6DV7AJk+50R7/mrpuNA0\nlJ+qZ/2aA1SVN+LhaeTSb41nwnRz8Hemt5ueW8ML2wupaGrHKGDF5OHcNi0Kbw8DtfszyXvlPUZ/\nsJFjHfVUG/wRHkbCUqcTuTSNYUvT8B810im6+kMFn1sVz18VHfZia7bPA1LKpzBX9+y5/GdSymfs\nbdQyNvBWe7fTXDhUpu+javs+Tj79KpmmRoavuheAsHnTCU91TucgY2+h2eZpNxEeGcDVt0ylqFPy\nxt4SskL9Kavz4Y295tLGU6IDmNKjtv1AqWhs44XthaTnmYddHDvMj5/OjyM+wINTH20i7+/vUbvv\nCGAuXdwYGceU3/2YiEvm4hk8+PY1Fy62Xhs+Cjx1huX/C9gd/M+GzvPvRgUd7tIQnmoO8ieffpXx\nBn9SVq10WlttbR1s/iiTI/vMNs+E6SNYvHw8Xl4eRABTRgRStvYgyUGJ3DnDMQ82dZok/zlWwau7\ni2lqN+HnaeCuWSNYEm6g+B9r2Pr6v2ktM4/s5RkSyMhbl7OhPBzfwBCir7nUIRoGgwq57ark+aui\nw17OGvyFEIssk0YhxCX0ztBJBOpP30qjGTxdN9QO3/UoAIGbzOPXOvqGWkVpA+vXHKCyrAEPTwNL\nlo9n4gznWig5Vc2s3pbP0bImAObFB/OdgAbqXnyRr9ZtRrabyyEHjE0kfuUKOqfPpqi4gRDLzUUV\nShpohj7n6vm/YvnrQ+8aPhI4BdzraEHa81dLh7s0dN1Q27HlJHlFmTywxPGu6uF9RWxal0lHeydh\nw/xZfutUIob3b6UM1udu7TDx1v5TvHeolE4Jw7wEK5tz8Hn2E45ZrB0MBiIvX0D83SsIS51urZmT\nMB7qb99KpqmRbytwc1F7/t2oosNezhr8pZSjAIQQb+jMHM35QntbJ5vXZ1ozeMZPG8GS5ePx8nZe\nhsz+onqeS8+nuK4N/4Y6bsjZR9SWzbSWV9EKeAQHEnvbcmK/cx1+do4opdEMBFuzfVwW+LXn340K\nOlTQ4Eh/ubKsgY/+ZbF5PAwsXj6eiTNibKpKORAdtS0d/G1nEZtOVDG8MI/r935J/P490NFBOxZr\n5+4biL7uMvMYtOdAFX9Ze/7dqKLDXs7l+X9Fj9r9Z0JKucChijQaJ3FkfxEbP+y2ea6+ZSrDopyT\nMSOlZHNWNS9ty2X4vr3c8vUXRBfkmN/sx9rRaFzJuXr+f3eJih5oz18tHSpoGKy/3NfmGTc1mku/\nNcFum8dWHcV1rbz4n8MY1m/gxl1fEdBQB5itnZG3Xk3cd647a9XMs6GKv6w9/25U0WEv5/L8X3eV\nEI3GGVSWmbN5KkrNNs+iq8cxaeZIp/S2O0ySD99Pp/Af7zP30D48Oi1ZO2NGEb9yhc3WjkbjCpR7\nBlx7/t2ooEMFDQP1lzMPFLPxwyO0t3USGuHH8lumMSx64DZPfzpMbe3s+9cGjv/tHcJyshkLSCEI\nXZpG8vdvJCx1hsNONqr4y9rz70YVHfaiXPDXaAZLe3snn398lEO7zdVHxk6OZum19ts856K1vIrs\n1/7NyVc/wKO6mjCgzdeP4BsuZ9aPbhmwtaPRuAKDuwX05cCBA+6WcFrZWHehgg4VNPQsIXwuqsob\neOuvOzi0uxCjh4FLr5nAlTdNdkjg79JRe+Aoh370Gz6ffg15T7+CR3U1FZHRVP/w+yw+uI6Ln7zf\naYE/09R47pVcgD2fibNQ5VioosNe+v1FCCG+llLOtUz/Skr5mOtkaTT2c/RgMZ/922LzhPtx9a1T\niYw+42ihdiM6O/EvzuHrq75PzZ7DgNnaOTluMhXLLuP2u5eSHDE0L/81FyZn6w6NFkL4SClbMJd0\ndknw155/NyroUEHDufzlvjbPmElRLL12It4+A+vtdzQ2UX/0JPVHsqg/kkXdkeOMPXicCR3t1ACt\nvr5kzJjH0XkLuWHpJP5n/DCMBteka6riL2vPvxtVdNjL2X4d64DjQohczDX8vzzTSjrPX+NOqioa\nWb/mAOUl9Rg9DFxy5VimzI616QarlJLWknLqDp+gPvOEJdCfoCmnEPoMcmQEqoZFs/eiizk6dTYz\nkobxbGoskQFeTvrPNBrn0m/wl1J+VwiRBiQAs+iu8+NUdJ6/WjpU0NBfTvmxQyX894PDtLd1EhJm\ntnmGjzizzWNqa6fhRC71h09QZwn09UdO0F5dd9q6wtODgNEJkDSKkxGjyPUIo9zbn+Lib4hPnsLt\nkf6kJYe4JfC7M6e8q9jesP1byaordnuBOVXy61XRYS/nyvPfBmwTQnjpnH+NKnS0d/L5f45xcFcB\nAKMnRnHZdd02T1tVbXdP/vAJ6jOzaDieY62W2RPP0CACJ6QQOCGZgPHJVEbHsssQxLaiRk7Vt1nX\n825qJLjgG37zzPcJcGINIJXpOXpVpamR1CUp7pbkNrrGm9hxyRUUVhVw4knz87DOHG/C0dha2+dV\ny6DrdwIxQBHwTynl544WpD3/blTQoYKGnv5ydWUj6/91gLKSeoxGwbzZw4kzVZK3+hXqLYG+pbjs\n9J0IgV9iLIHjkwmamELg+BSCJqbgMTyCw6WNbM6tIT23lsrj7UA1ACE+HqQmBJOWEMKp2VdgNJnc\nHvhV8ZdV0OFODV3jTez4+34AUlZOc5uWgWLrSF4rgccxl3vYCcQBa4QQj0gpX3aiPo0GMN+EPbgx\ng6/2VNFhEvi0NhD7xVqqX86zhOpujL4+BIxPImhCirVXHzguCQ9/PwDaOk0cKK5nW1YtO7Ycobal\n+4ogMsCT1IQQ0hJCGB/pb72Ru8Fk4kJn10fbKPhiD9x6IwBrf7YagNiFM5m93P2dBI192NqN+Tlw\nqZTyYNcCIcQ7wFrAocFfe/5q6XCXBikltfuOELXzvxTk7ueNo7upGjcLEATlZBKzbT3G9lZ8RkQS\nOD6ZwAnJBE0YTeCEZPwSYhBGY6/9tXSY+Dqnhm25NXydX0tTe3cwjwnyJm1UCPMTQkiJ8O33ZrEK\n3q47NcxenmYN8ip8N1X4PADqTh4AztOePxAO9H2q4xsgzLFyNBcyXQG/5KPNFKzbQkdZJd4jEqmY\nNA/PcbMQpk5iG/OYuSCa4HueJHB8Cl5hwf3ur7Gtk535tWzLrWV3YR2tHd0Bf1SoD2mjzD38hFAf\nXVlTc8Fha/DfBjwjhHhQStkkhPAH/gBsd7Qg7fl3o4IOZ2voGfBLP/6C5uIyGqPiqU2aTd2SCXR6\n+RANBIf4cPWt04gaeeVZ91fX0sH2vFrSc2vYV1RPu6k7ZXPMMD/SEkJISwgmJtjHbq0Xus/dExW+\nm6oci6Ak98esgWBr8P8B8A5QK4Sowtzj3w7c4ixhFypdWQQnLeO1Jt1/F+DaLAJnazgt4BeV0hQZ\nS23iNOoWTqTD2++0be64NxUfX88z621qJz23hvTcGg6WNNAV7wUwKSqAtIRgUhPck5qp0aiKrdk+\nJcACIcRIYARQLKUsdIagC93zLxyVwkGfaLL2n6KwsoCFi8w93SnRAbgqk7oxOp6yaYEcvmskeUWZ\n3DztYgC8o8MGrOFMAb8lPJraxEnULbiVNt/uapsh4X6MnRzNmElRvP7ndPKKMvHx7e3unqpvZVuu\nuYefWdpoHXHIKGBGTCBpo0KYFxdMqN+ZTxgDwV0ec9+0wiQF0grd9RtR8Vic754/AJaAP6CgL4Tw\nBr4EvCztvq/rBZ3OlBGBTBkRyIYtn5BpauTOGa4fz7Xn4OnAgPO5+wb8lqJSWkKGUZs4gbq0m2n1\n7/brA4N9GDM5irGToxk+IuiMHnx+TQvpuTV8lVNDVmWzdbmnUTAzJoi0UcHMiQ0maIBlHVSlZ1ph\n3ckDpKz6rrsluQ19LByHy34lUspWIcQllnsGRiBdCPGplHJXz/W059+NCp6mvTVczhTwWwNDqU2c\nQP28G2gO7L528PP3InnicBLGDScoKpDmDhPVHSaKi+tpbjfR3G6iINCXtomz+d77R8mrabFu6+Nh\nYE5sEGmjQpg1Mgg/L+OZ5DgUd30eB4vrOVjSwO3TomDaMt7YWwKYrwanjHDOMJTnQoXfiCpeuyo6\n7MWlXSQpZZNl0tvS9lnHB9aoRVuniRZLUG7u6LQE6E6a2jppPniU9k1fwdZ0DOUVtPkHUTdqAtWz\nr6E1dLh1H51GA9VBvpwK8KHU00hneRuUF/Tf6DBLuYaaFgK8jMyND2Z+QgjTYwLx9lCuIrlT6Loa\n1KhDrxMyKHFCtpdzBn8hhAFYCGyTUradY3Vb9rUXSAJekFLu7rvOhe7592SgHrOUknaTPGOg7upR\nt3T0nO+kucOyvN1EU3snLR0mCkeGU5l7iB3/PERzu4mOHpkzSEl0YS6jM/aRcmQ/QbXVmHz8qR01\nnqqLrqE1oruefYcQlPl7cyrAh0pfL2QPS8cowM/LiK+nAV8PIz6eBuu0r6eB4/uLqMo5xP8+eBtT\nRwTi4aLqmWdChbxyVb6bKuhwp9fe84SswrEYCOcM/lJKkxBinZRy0KczKaUJmCaECAI+FEKMl1L2\nen5g69at7Nmzh7i4OACCg4OZNGmS9eB2DS7izPmMjAyXttdzfu27n1BWXEfF935NflEmj/72dTol\njL0oFRnmz5G9X9PaaWL42Bm0tJs4eWgXbZ2SoKSpNLd3UnJ0L60dJvwTzZei5h9I96WpXfNeHtSU\nZtNmNBCUNBUDEv/dnzEy+zgXF5bhV11FhqGdzKh4IudfQ9OwkeQVH4XWGhKMMQSMDKGyJZeQ6EAu\nS52Bj6eBI3u/xtvDwMXz5+PraWDnju1nPR73vvYBDRW5zBwZ5JbPo+9gNu5qv2s+IyPDre27+3i8\n/uFnZFU1c/u02ZwwhPHIq+sAuGHZIqaMCHTL8XBHvOiazs/PB2DmzJksXrwYexBSntt5EUL8B/it\nlPJru/Z+9n0+AjRKKZ/puXzz5s3S3T1/d1LW0MYTX+Rx6FTDoPZzrh61+WXEx6N72rfHOj6eBt55\nYTtGk4m7Lo+i5tMvKP/E7OF3enpRHzeG+nHTqYuMQ2LujRuMglEpEYydHE3SuEiHjJ711MMbAHjg\ncff2uTdEzQNg2SmHP9qi0Qyaffv2sXjxYrsui239deYBnwoh1gEF9PDqpZSP2rIDIUQE0C6lrBVC\n+AKXAn+0R+z5zpfZ1azeVkDAySwu//pL/FpbCAr2wSgEHgaBQQiMBjAaBEYhevztu4xBP7HaAow7\nfArf8iIyXq7DZPSgPjaFhiuWUBs1CpNlBFBhECQkhTFmcjQp44f3m4uv0WjUwtbg7wt8aJkeOcC2\nooHXLb6/AXhHSvlJ35UuRM+/ub2TF3cU8sXBQtI2rmPy7nSElGSaGonvJ8NEAh2Wl7MINBrZGxpM\nzPxbqB2RRGfXkM8CRsaHMnZyNCkTh+Mf4O1EFf3X83cFXTf2shZdQWFlAWVuvrGnir+sgg4VNKik\nw15sfchr0Mm0UsoM4ML1c/rhm/JG/rAll6Bt6Xz307X4NdQjPIyUjZtNqaGdW++5xqntd5okDa1Q\n3yJ7vepaJI0tkrLio/iOMOf5R40Mtj58FTiA8ghDERWeu9BonIFNnj+AEGIssAIYLqX8kRBiDOAt\npTzkSEEXiuffaZK8e6iUj/57kIUfvUP8yWMAhMyezIQ/reJv/zwBOMbr7ugwUVvVRE1lE9WVTVRX\nNlqn62ua+45YeBrzLxvNmElRhISdXnbBWXT1uLdvyQJg3qJkwPU97q4nSvsylAbt0Jz/OM3zF0Ks\nAF7EXML5VuBHQCBmz36JnToveMoa2nhyYxa+733AbV9+hkdHBx4hQYx99B5ibr4SYTAAJ+za50AD\nvBAQHOZLaLgfIeH+lr9+hEb48+ozXwEw5+LEQf7H9tPV4y5ba64i7q4ed9cTpRrN+Yatnv9vgCVS\nyoNCiJssyw4CUxwt6Hz3/L/Mrubd1zYy74N/EVZhHnFqxI1XMPbRe/CKCO21bl+v2xkBPjjEF+NZ\nHpZyp9+ukgZQw9tVQYMqOlTQoJIOe7E1+EcCXfaO7PFXP6FrI83tnfxtwxFMf3mFqw6an23zSYpj\n0hM/P61n2dnZXXd+47ojgwvw4X4Eh/qdNcBrNJoLD1uD/17gDuCNHstuBnadefWBcz7W9jlWWs97\nj7/JpHXv49PSjPTyIuWn3yHxntswePVOjdy8I59DX5gLqsXHjOfgzu7SB0JAcKgvoRGuC/D21vZx\nBipoADXq2aigAdTQoYIGUEeHvdga/H8MfCaEuBvwF0L8FxgNLHWaMhfjjBr2nSbJ2g930vDH55mV\nnwOAf9osZjy1Cr+E3hmzlWUNfPHJMXKOV/RafsmVY3UPXqPROBybIomU8hgwFngB+F/gH8AkKaV9\ndyVt4MCBA47epU2Ep04nZdVKwFzDJWXVSlJWrRxw4C8pq+Xl7/6OgB89wIj8HDpCQ5jw18dIe291\nr8Df3NTG5vWZvPbndHKOV+Dl7cHFl48BzF73jNQEksZGEjYswC2B3+y3uxcVNMDppQ0uVA2ghg4V\nNIA6OuzF5ufvLaWY04EczIO5DK7+wHnMltc/pfIPL5BQU4UUAp8bl5P6ux/hGRRgXaez08TBnQVs\n35xFS3M7QsCU2bHMW5KMf4A3Wz/9xo3/gfvpuhKL3JdFbV0xJxQYtEOjOZ+wtbZPHPAWMBeoBkKB\nr4HbpZR5jhTkrjx/65OcT78CQPL9dwP25ZVX5xWz8Ud/JGj3HgAa4uOZ89wviJ07udd62d+U88Un\nx6gqbwQgLjGMS64cx7DowO6Rijab89svWmzOb3fHMI7u1JCfXUlBdtVpy2MTw4hLdNWYZhrN0MCZ\ntX1ex3zTd5mUslEIEQD81rJ8oV0qFaXnk5wAy/71vzZva+roYM/qtyj982sEtbXS5uWN4Xt3cO1D\nd2D07L6h29fXDwnzY+EVY0gaF2mtxdOVV76u2lzQLGWV61Mcu4ZxTLIM31hmWT6YYRztpWs0MY1G\n4xxsDf4zgKVSynYAKWWDEOJBoNLRgtyV59/V0zx816PkFWUSuMl8O+NcPc2qPYfZ/uPH8cjOxRMo\nmjqDBatXkTI2zrpOc1Mb2zdncWBnAdIk8fL24KJFSUy7KB4PBXPsewZeFXKYVdCgig4VNKiiQwUN\nKumwF1uD/9fAbCC9x7KZwA6HK3IT9o5b215bz8HHXqB8zXo8pKQ2JJzme77HHf9zFV6WgH4uX1+j\n0WjcRb/BXwjxmx6zJ4FPLHX9C4BY4ArgX44WpEKe/9nyyqWUlHy4kUO/XA1VNZgMBo5cchkLH/sB\ns1Mireudzdd3hA5XoUKPRgUNoIYOFTSAGjpU0ADq6LCXs/X8Y/vMf2D5Gwm0Av8GLozSjhYasws4\n/OBTVH9lfkK3MD6Jsu+v5H9unUeopY69Lb6+RqPRuJt+g78jyjgPBBVq+/T12k2tbWT/5U1OPvc6\nsq2dZl9/tl9xLfP/53runmAO6gP19fvSde/hokVJZBzZS7qN9x6chQp+pgoaVNGhggZVdKigQSUd\n9mJznr8Qwg9IBgJ6LpdSntfj2lVu28ORnz9JU7a5zMLhaXPJu+lm7v/WJBJCfS2+fr7DfP2eN1ul\nTympaWe/96DRaDQDwdY8/zuBvwBtQHOPt6SUMu7MWw0Md9fz7xoz9t6fzuabx56n+P3/AlA5bDib\nl9/CrKvmcfesEXh5GBzi62s0Gs1gcWae/xPA9VLKjfbLGmJISejx/Wyb/yztNfV0eHjy9cJlnLx0\nGT9blMSs2CDt62s0miGPrWZ0G/CFE3VYcVdtH4CW0gpG/ecf1H71Du019eSmjOP1e3+JuONGXrxx\nIhPDfc5Yh+c796WRPH64wwO/CjVDtIZuVNChggZQQ4cKGkAdHfZia8//EeAZIcRjUsqKc649RPEK\nC4H2Nlp8/fn4mrvInTKD788dyRWjwzi0q1Dn62s0mvMGWz3/i4C3gZ51iAVmz9/oSEHu9PxP1bdy\n3/ObaQgKISY6lF8sSsBksXi0r6/RaFTFmZ7/PzEP5PIOvW/4nldEBXoT4B1MQKuJx+bFkL7uiPb1\nNRrNeYmtnn848KiU8rCU8mTPl6MFudPzBxhTWY/voV289cIOl/j6Z0MFL1Fr6EYFHSpoADV0qKAB\n1H+S8KAAABEvSURBVNFhL7b2/P/B6cM42oUQYqRl++GACXhZSvnnge7PGdTVNGMN7VJqX1+j0Zy3\n2Or5b8Nc2C0HKO35npRygU0NCREFREkpD1hKQu8FvmUZJcyKOz1/KSVP/9Kc1//te1O1r6/RaIYE\nzvT8X7a8BoyU8hRwyjLdIIQ4CsQAx866oQvpaenowK/RaM5nbB3D9/X+XgNpVAiRAEwFdvZ9z92e\nP+gxY7WGM6OCDhU0gBo6VNAA6uiwF5t6/kKIu/p7T0r5qj0NWiyf94GfnGkc4K1bt7Jnzx7i4sxV\nI4KDg5k0aZK1cFLXgXbWfF5RJqUVuVY9zm5P9fmMjAy368nIyFDmeLh7XoXPoyfuPh4qzLvj+9k1\nnZ+fD8DMmTNZvHgx9mCr5/95n0VRQBKQLqW8xObGhPAAPgY+lVI+d6Z1VKnt88Djrh9BS6PRaAaC\n0zz/MwV4y9XAOHsaA14FMvsL/BqNRqNxDbYXmj+d14C7bV1ZCJEK3AYsEkLsF0LsE0Kc1r3Wnn83\nKniJWkM3KuhQQQOooUMFDaCODnux1fPve5LwA24HamxtSEqZDji0FIRGo9FoBoatnr8J6LtiEfA9\nKeV/HSlIe/4ajUZjH87M8x/VZ77xfKvuueujbRR8sQePnCoA1v7M/PhB7MKZzF4+9IZo02g0mrNh\na55/Xp+X0wK/uzz/2cvTuP6Z++gYN4+TQSFc/8x9XP/MfW4N/Cp4iVpDNyroUEEDqKFDBQ2gjg57\nOWvP35LieTZfSEop7Usu1Wg0Go3bOavnL4ToL5snBvgx4Cel9HOkIO35azQajX043POXUr7Sc14I\nEQ78Avge5tr+v7FXpEaj0Wjcj02evxAiSAjxWyALc0nm6VLK70spCx0tyF2ef352JembTnDRoiQC\nhteQvukE6ZtOkJ9d6RY9oIaXqDV0o4IOFTSAGjpU0ADq6LCXc3n+vsB9wP2YB3BPk1IecYEulxOX\nGE5cYjgA0qeU1LQUNyvSaDQa53Euz78U89XBk8CeM60jpdziSEHu9vw1Go1mqOGMPP9mzNk+/9PP\n+xJItKdBjUaj0bifs3r+UsoEKeWos7wcHvhVqO2jioengg6toRsVdKigAdTQoYIGUEeHvQymsJtG\no9Fohig21fZxJdrz12g0GvsYiOeve/4ajUZzAaJc8Neefzcq6NAaulFBhwoaQA0dKmgAdXTYi3LB\nX6PRaDTOR3v+Go1GM8TRnr9Go9FobEK54K89/25U0KE1dKOCDhU0gBo6VNAA6uiwF+WCv0aj0Wic\nj/b8NRqNZoijPX+NRqPR2IRywV97/t2ooENr6EYFHSpoADV0qKAB1NFhLy4L/kKIV4QQpUKIQ65q\nU6PRaDRnxmWevxAiDWgA3pBSTu5vPe35azQajX0o7flLKbcB1a5qT6PRaDT9oz3/M6CKh6eCDq2h\nGxV0qKAB1NChggZQR4e9nGskL5ezdetW9uzZQ1xcHADBwcFMmjSJtLQ0oPtAO3M+IyPDpe2pPJ+R\nkeF2PfrzUOvz6Im7j4cK8+74fnZN5+fnAzBz5kwWL16MPbg0z18IEQ+s156/RqPROA6lPX8LwvLS\naDQajRtxZarnv4DtwGghRL4Q4rtnWk97/t2ooENr6EYFHSpoADV0qKAB1NFhLy7z/KWUt7qqLY1G\no9GcHV3bR6PRaIY4Q8Hz12g0Go0CKBf8teffjQo6tIZuVNChggZQQ4cKGkAdHfaiXPDXaDQajfPR\nnr9Go9EMcbTnr9FoNBqbUC74a8+/GxV0aA3dqKBDBQ2ghg4VNIA6OuxFueCv0Wg0GuejPX+NRqMZ\n4mjPX6PRaDQ2oVzw155/Nyro0Bq6UUGHChpADR0qaAB1dNiLcsFfo9FoNM5He/4ajUYzxNGev0aj\n0WhsQrngrz3/blTQoTV0o4IOFTSAGjpU0ADq6LAX5YK/RqPRaJyP9vw1Go1miKM9f41Go9HYhHLB\nX3v+3aigQ2voRgUdKmgANXSooAHU0WEvygV/jUaj0Tgf7flrNBrNEEd7/v+/vTMPumu84/jni6AJ\nItQWkcRSS5WgGYRqKbUX7TDEKBo1nU6R2mprJ8q01XZoTbfR0hQt6VAGxVRUWrXFktUasSQhC1kk\nxBbJt388z5uc9/a+b+5Lc87R9/eZuXPPec55nud7n3Pu7z7Lub9fEARB0BKlGn9JB0t6VtIUSec1\nOyfm/FdQBx2hYQV10FEHDVAPHXXQAPXR0VVKM/6SVgN+BRwE7AgMlbR943lTp04tS1KHTJ48uWoJ\nQD10hIYV1EFHHTRAPXTUQQPUQ8eH6TSX2fPfHXje9jTbS4BRwJGNJy1evLhESc1ZuHBh1RKAeugI\nDSuog446aIB66KiDBqiHjokTJ3Y5T5nGf3NgRmH/lZwWBEEQlEztFnxnz55dtQSmT59etQSgHjpC\nwwrqoKMOGqAeOuqgAeqjo6usUWJdrwL9C/v9clo7tt56a4YPH758f9CgQeyyyy6rXl2BwYMHM27c\nuFLrrKuO0FAvHXXQUBcdddBQlY4JEya0m+rp1atXl8so7Tl/SasDzwH7A7OAR4Ghtp8pRUAQBEGw\nnNJ6/raXSjoNuIc03XRNGP4gCIJqqN0/fIMgCIJVT20WfCVdI2mOpEkVaugn6T5JT0maLOmMCjSs\nJWmspPFZw4iyNRS0rCZpnKTbK9TwsqSJuT0erVBHb0k3SXom3x97lFz/trkNxuX3hRXdn2dKelLS\nJEl/lrRm2RqyjuH5+1Hq97SZnZLUR9I9kp6T9HdJvSvQcHS+LkslteQfpzbGHxhJ+gNYlXwAnGV7\nR2AI8O1mf0Rbldh+D9jP9q7ALsAhknYvU0OB4cDTFdXdxjJgX9u72q6qHQCuBO6yvQMwCCh1ytL2\nlNwGuwGfBRYDt5apQVJf4HRgN9s7k6aNjytTQ9axI3AKMJj0HTlc0lYlVd/MTp0P3Gt7O+A+4IIK\nNEwGvgL8q9VCamP8bT8ALKhYw2zbE/L2W6QveOn/RbD9dt5ci/QFK31uTlI/4FDg6rLrbpRCxfep\npPWAfWyPBLD9ge1FFUo6AHjB9oyVnvm/Z3Wgl6Q1gJ7AzAo07ACMtf2e7aXA/cBXy6i4Azt1JHBt\n3r4WOKpsDbafs/086fvSErUx/nVD0kBSr2JsBXWvJmk8MBsYbfuxsjUAPwfOpYIfngYMjJb0mKRT\nK9KwJTBX0sg87fI7SZ+oSAvAscCNZVdqeyZwOTCd9Jj2G7bvLVsH8CSwT55u6UnqpGxRgY42NrY9\nB1IHEti4Qi0tE8a/CZLWAW4GhucRQKnYXpanffoBe0j6dJn1SzoMmJNHQaILvYlVwN55quNQ0jTc\n5yrQsAawG/DrrOVt0lC/dCT1AI4Abqqg7vVJvdwBQF9gHUnHl63D9rPAT4DRwF3AeGBp2To6oeoO\nU0uE8W8gD2dvBq63fVuVWvLUwhjg4JKr3hs4QtKLpB7mfpKuK1kDALZn5ffXSXPcVcz7vwLMsP14\n3r+Z9GNQBYcAT+T2KJsDgBdtz8/TLbcAe1WgA9sjbQ+2vS/wBjClCh2ZOZI2AZC0KfBahVpapm7G\nv+peJsAfgKdtX1lF5ZI+2fa0QJ5a+BLwbJkabF9ou7/trUgLevfZPrFMDQCSeuZRGJJ6AQeShvyl\nkof0MyRtm5P2p7qF8KFUMOWTmQ7sKWltSSK1QyX/1ZG0UX7vT1rovKHM6mlvp24HTs7bJwFldBo7\ns5Wt2VDbtXiRLt5M4D3STfb1CjTsTRo+TiANJccBB5esYadc7wRgEnBRxdflC8DtFdW9ZeFaTAbO\nr7AdBgGPZT23AL0r0NATeB1Yt8J2GEEy+JNIi5s9KtJxP6kjMJ70NFhZ9f6XnQL6APeSPBjcA6xf\ngYajSI4z3yF5ULh7ZeXEn7yCIAi6IXWb9gmCIAhKIIx/EARBNySMfxAEQTckjH8QBEE3JIx/EARB\nNySMfxAEQTckjH9QKdlfziUV1z9f0iNVaQiCKgjjH7Qj+8+fU3RcJukUSWOq1LUqyH6C9gf62t6z\ng3M2lXS1pJnZh/7TkkZU7NitNPJnrcS1R7BqCeMfNGLSffGdJum1RlJX7+eBwMu23+2gvD7AwyTX\n2nvY7k1yt9Eb2PojSA2CygnjHzTjZ8DZ2Y99OyQNkLSsaGgljZE0LG+fJOkBSVdIWiBpqqQhOX26\npNmSGv0EbZQjIS3KZfUvlL19PjYvR9E6pnBspKTfSLpT0pvAvk30bibptpx/iqRv5PRhwO+BIbne\nZhHTzgYW2f6as+9826/aPsv2k7mcvSQ9mj/rWElDGtrlUkkPSnoz69hA0p/yKGJsw2ddJul0SS9I\nek3STwvHJOl7eWQ2W9If265P4ZqcKGlaznthQ97z87V4XdKo7KGz07ySDgIuBI7N+sfn9JOzxkX5\nfWiTtgvqTlU+QuJVzxfwEvBFkufKS3PaKSTnbpDc+S4FVivkGQMMy9snAe8DJ5IcTF0KTAN+CfQg\n9ZwXAT3z+SOBhSS/Sj2AXwD/zsd6knyXtJU1iOTbZvtC3gXAnnl/zSaf5/5C3YNIHhf3LWi9v5O2\neBgY0cnxPsB84HhSR+q4vN+n0C5TSCOMdYGnSE769svnXwtcUyhvGfAP0siiH8lXTFu7DstlDcjt\n8lfgusI1WQZcBawJ7Ay8C2yXjw8HHgI2y+3wW+CGFvOOaKuncE0WAtvk/U2AHaq+b+PV9Vf0/IOO\nGAGcJmnDD5H3JdvXOVmHv5AM2Q9sL7E9mvTjsE3h/DttP2h7CXARyXPk5sDhxbJsTyQZvWMKeW+z\n/QiA7feLIpSikQ0Bzst1TyRFJmvVQ+mGJCdZHXEYMMX2DU4xGEaRjPuXC+eMtP2y7TeBu0kRuMbY\nXkbyyb9rQ5mX2V5o+xXSD2Fbr/p44Arb05wivV0AHFcYgRm42Pb7ticBE0k/dgDfJDkInJXb+BLg\n6BbzNmMpsJOktW3PsV2JZ8/goxHGP2iK7aeAv/Hh4pHOKWy/k8ub25C2TmF/eThC24tJvfm+pF7p\nnvlpnPmSFpCM4CbN8jahLzDfK8JiQhqFtBqacx6pt9xZ+dMa0hrLb2yLxv1iO0CKHVAsq28HdU0j\nBZkptkWx7LcLZQ8Abm1rR5I76iUt5m1HbstjgW8BsyTdIWm7ZucG9SaMf9AZFwOn0t6YLc7vPQtp\nm37EepaH4FPy39+H5LJ2BvBP2xvkVx/b69k+rZC3s4XomcAGSrEA2uhPCkHYCveSfMV3Vv7AhrSu\nlN+MYjjCAayIkTsz7xePLaG90e6I6cAhDe3YyzlQzkr4r/a1Pdr2gaTr/hxp7ST4mBHGP+gQ2y+Q\npm3OKKTNJRm3E5RiDQ9j5U++rCy4xKF54XRN0hrBI7ZfJY08tpV0gqQ1JPWQNLjVnmaeOnkI+LGk\ntSTtTFq/uL6V/MAVwHqSrm1bmJW0uaTLJX2GFELwU5KOk7S6pGNJwcXvaLH8ZpwraX1JW5DafVRO\nvxE4U9LA/AP5Q2BUnj6Cztv4KuBHhc+wkaQjCsc7yzsHGChJOe/Gko5Qip27BHiLeoVQDFokjH/Q\nSGNP7xJSL7+YfirwXWAuydg92MUy3bB9A2mUMY80B34CgFP85ANJC6kz8+sy0qOXrTKUFBRmJmm9\n4Pu2W/rPgu0FpDCFS4CxkhaS4sa+AUy1PZ+0LnEOqS3OAQ7L+Ro/Z6vcBjxBCuhzBymyHPn9etIC\n9gukqZkzCvk6a+Mrc7n35M/wEO3DYXaW9ybSj8M8SY/n7bNIHYC5wOdJU0DBx4wI5hIENUHSMtJT\nNC9WrSX4/yd6/kEQBN2QMP5BUB9iGB6URkz7BEEQdEOi5x8EQdANCeMfBEHQDQnjHwRB0A0J4x8E\nQdANCeMfBEHQDQnjHwRB0A35D6JlxxFkC0buAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7faa1c05d710>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib\n",
    "import matplotlib.pyplot as plt\n",
    "matplotlib.style.use('bmh')\n",
    "ax = plt.figure()\n",
    "plt.errorbar(np.arange(2, max_comps+1), np.nanmean(comp_iters_pseudo, axis=1),\n",
    "             yerr=np.nanstd(comp_iters_pseudo, axis=1), label='Pseudo')\n",
    "plt.errorbar(np.arange(2, max_comps+1), np.nanmean(comp_iters_quasi, axis=1),\n",
    "             yerr=np.nanstd(comp_iters_quasi, axis=1), label='Quasi')\n",
    "plt.errorbar(np.arange(2, max_comps+1), np.nanmean(comp_iters_normunity, axis=1),\n",
    "             yerr=np.nanstd(comp_iters_quasi, axis=1), label='NormUnity')\n",
    "plt.ylim((0,None))\n",
    "plt.xticks(list(range(1, max_comps+2)))\n",
    "plt.xlabel('Number of Components')\n",
    "plt.ylabel('Number of Iterations')\n",
    "plt.legend(loc='best')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAESCAYAAAD5d3KwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XucVWW9+PHPd5BMSBCREXREQEPoOKF2SMlBH4XMUtFU\nVLyUmb9fZmbeKtSfly7HsjK11GOe1OMNSNATaJY3fFK6qAXUqChHAWFE8ALeL4nz/P5Ya4Y9e561\n9wL2nu+G+b5fr3mx15q91/OZvYd5Zq89ey0JIWCMMca0qdMOMMYYU1tsYjDGGNOBTQzGGGM6sInB\nGGNMBzYxGGOM6cAmBmOMMR102cQgIteLyEoR+WfBun4icp+IPCMi94pI34LPnSsi/ysiC0TkgK7q\nNMaY7q4rnzHcCHyuaN1k4IEQwi7AbOBcABH5BHAUMBL4PHCNiEgXthpjTLfVZRNDCGEOsLpo9aHA\nTenlm4DD0ssTgGkhhDUhhCXA/wKf7opOY4zp7rRfY6gPIawECCGsAOrT9dsDywqu90K6zhhjTJVp\nTwzF7PgcxhijbDPl8VeKyLYhhJUiMhB4KV3/ArBDwfUa0nWdfP3rXw/PPfccAwcOBKB3797svPPO\n7LbbbgDMnz8foNNy27qsz2ssF7dp97QtP/vssxx55JE101N4H9VKT60+fjNmzMj1/6G7fz91h8dv\n/vz53HvvvQAMHDiQUaNGcfbZZ0dfu5WuPIieiAwB7gohNKbLlwKrQgiXish3gX4hhMnpi8+3AXuS\n7EK6H/h4iMQ++OCDYdz48WXHvvuWZRxyQkPBmovTj7VueeRuTnjwkJLbmXX4LJoamsqO12POHPpM\nmLAJNHXuytOUt6tSTQCzZr1BU9OHZbuq5cc//jGTJ09WGz/GmvKrZlctfp/PnTuXcePGRSeGLnvG\nICJTAAf0F5GlwEXAj4HpInIS8DzJXyIRQnhKRG4HngI+AE6NTQobZkllN1cRS7QDMizRDohYoh3Q\nydKlS7UTOrGm/Gqza4nKqF02MYQQjs34VPTX/RDCj4AfVa/IGGNMTK29+NyFTtQOiDhROyDDidoB\nESdqB3Ry7LFZv/vosab8arPrRJVRu/HE4LQDIpx2QAanHRDhtAM6aWoq/zpPV7Om/Gqzy6mM2o0n\nBq8dEOG1AzJ47YAIrx3QyZw5c7QTOrGm/Gqzy6uM2o0nBmOMMTHdeGJw2gERTjsgg9MOiHDaAZ3U\n4q4Ia8qvNrucyqjdeGIwxhgT040nBq8dEOG1AzJ47YAIrx3QSS3uo7am/Gqzy6uMqn1IDGMqZrMB\ni5jTsqzkdRq2bGBI3yFdE1RBdUuWIC0tZa/3XN/daXm9T/tyc3Md0KPDdRoaAkOGtHZZV2hooHXI\nkIqMl0ctNm1suvHE4LQDIpx2QAanHRDhOq15dU0LJ9xZ+rADsw6fVbWJoZr7qKWlpewhFQBeuGUZ\nE07oU7DmoE7XmTXrDSr1MzFP1xuzZlE4YLX35a9PE9hrDIW68a4kY4wxMd14YvDaARFeOyCD1w6I\n8NoBndg+6nxq836q1S6vMmo3nhiMMcbEdOOJwWkHRDjtgAxOOyDCaQd0Yvuo86nN+6lWu5zKqN14\nYjDGGBPTjScGrx0Q4bUDMnjtgAivHdCJ7aPOpzbvp1rt8iqjduOJwRhjTEw3nhicdkCE0w7I4LQD\nIpx2QCe2jzqf2ryfarXLqYzajScGY4wxMd34nc+e2vttylN7TVCbXZ5qNeU5pELxoScAmpsfobFx\nbId1lTz8xPrxaD92izYbzrI5aw/LUQv3U3FTrEv/sQOtx68bTwzGxOU5pELnQ08A9AY6rqvk4Sc2\nVi2v9i66r/Tvp85NUNzVnR+7brwryWkHRDjtgAxOOyDCaQdEOO2ACKcdEOG0AzI47YAIpzJqN54Y\njDHGxHTjicFrB0R47YAMXjsgwmsHRHjtgAivHRDhtQMyeO2ACK8yajeeGIwxxsR044nBaQdEOO2A\nDE47IMJpB0Q47YAIpx0Q4bQDMjjtgAinMmo3nhiMMcbEdOOJwWsHRHjtgAxeOyDCawdEeO2ACK8d\nEOG1AzJ47YAIrzJqN54YjDHGxHTjicFpB0Q47YAMTjsgwmkHRDjtgAinHRDhtAMyOO2ACKcyqr3z\n2Zgq2mzAIua0LCt5nYYtGxjSd0jXBFGbTaa2dOOJwVN7vyF4aq8JarPLszE0vbqmhRPuLH14jVmH\nz6riD+GNo6k2eGqvy6PR1I13JRljjInpxhOD0w6IcNoBGZx2QITTDohw2gERTjsgwmkHZHDaARFO\nZdRuPDEYY4yJ6cYTg9cOiPDaARm8dkCE1w6I8NoBEV47IMJrB2Tw2gERXmXUbjwxGGOMiamJiUFE\nzhSRJ0TknyJym4h8RET6ich9IvKMiNwrIn0rO6qr7OYqwmkHZHDaARFOOyDCaQdEOO2ACKcdkMFp\nB0Q4lVHVJwYR2Q74JrBHCOGTJH9COwmYDDwQQtgFmA2cq1dpjDHdR+bEICJ1eT4q1NED6C0imwFb\nAC8AhwI3pZ+/CTisQmOlfGU3VxFeOyCD1w6I8NoBEV47IMJrB0R47YAMXjsgwquMWuoH+xrggxwf\nGySEsBy4DFhKMiG8HkJ4ANg2hLAyvc4KoH5DxzLGGFNeqXc+Dy24fBBwJPAj4HlgR+C7wB0bGiAi\nW5E8O9gReB2YLiLHAaHoqsXLAMyYMYNbgCHp8lbAbqzdM+c73aJtTfE10uXF6WLbV1+03Px4MyyB\npqYmAObMmQN0Xt4332hFa1ynayyY15w0ZPSwOGlqaijd07Zcrqd5wV+AAZFrrO1dMK957WKZ+6tc\nzyPNzfQu0ePbmyaWuEbBcgUev7rmZg4qM9pahddwnW5Rqccv7/dT58ev7Vprb5Hr8UtV7/Gjw/Jm\nA7ZiTsuy5PEBGkc3Jl9PwXLDlg20NLeU7Kn041fu68+7vH4/DwrXrb1Fc/PbNDV9Zp175syZw5Qp\nUwAYPHgw9fX1jBs3jpjMiSGE8HzbZRE5C/j3EMJr6aqFIvI34G/Af2ZtI6fxwKIQwqp0rP8BPgOs\nFJFtQwgrRWQg8FLsxkceeSTjrr02c+Mu/ffuTmsylodScrlxdGP7f2JY+wB0Wk4fkDKjlb3GyN0b\nYVXpvrb/RCV7co0GjSPHAA0lbzFy97fgweyeUuMXL49tbKRPmb63Ro4pc43s8dfn8euRY7S8308V\ne/xyfj9trI/fq2vu7niYjuJDOS1LDtNRrqfSj1+e8XItV+jnATgaG98APlznnqampg7Lc+fO7TR6\nm7yvEfQFehWt65Wu31BLgb1E5KMiIsA44ClgFnBiep0vAzMrMFYBX9nNVYTXDsjgtQMivHZAhNcO\niPDaARFeOyCD1w6I8Cqj5j2I3k3AAyJyBclcvgNwOmtfHF5vIYTHRGQGMI/kNYt5wHXAlsDtInIS\nye6rozZ0LGOMMeXlnRi+AzwLHA1sB7wIXAX8VyUiQgjfA75XtHoVyW6mKnHV2/R6c9oBGZx2QITT\nDohw2gERTjsgwmkHZHDaARFOZdRcE0MIoRW4Nv0wxhizCcv9PgQROUBEviMi3y/8qGZcdXntgAiv\nHZDBawdEeO2ACK8dEOG1AyK8dkAGrx0Q4VVGzfWMQUSuItnH/xDwTsGnon9CaowxZuOV9zWGY4FR\nIYTS5wPcqDjtgAinHZDBaQdEOO2ACKcdEOG0AyKcdkAGpx0Q4VRGzbsr6RXgtbLXMsYYs9HLOzFc\nBtwmImNEZFjhRzXjqstrB0R47YAMXjsgwmsHRHjtgAivHRDhtQMyeO2ACK8yat5dSW3vbj64aH2A\nDm80NMYYU0WbDVjEnJbSe/UbtmxgSN8h6z9GniuFENQPz115TjsgwmkHZHDaARFOOyDCaQdEOO2A\nCKcdkMFpB0S4TmteXdPS8fAhEbMOn7VBE8Mm+APfGGPMhsg1MYjIIyLycOyj2oHV47UDIrx2QAav\nHRDhtQMivHZAhNcOiPDaARm8dkCEVxk172sMvy5aHgh8Fbi1sjnGGGO05X2NodPB8kTkDuBGYCN9\n97PTDohw2gEZnHZAhNMOiHDaARFOOyDCaQdkcNoBEU5l1A15jeEF4JOVCjHGGFMb8r7GcFLRx2nA\n74C/Vjevmrx2QITXDsjgtQMivHZAhNcOiPDaARFeOyCD1w6I8Cqj5n2N4YSi5beBPwOXVzbHGGOM\ntryvMexX7ZCu57QDIpx2QAanHRDhtAMinHZAhNMOiHDaARmcdkCEUxk17zMGROTjwCRge5LXF6aG\nEP63WmHGGGN05H2N4RDg78AIkjOr7QL8TURKv/2upnntgAivHZDBawdEeO2ACK8dEOG1AyK8dkAG\nrx0Q4VVGzfuM4RLg0BDCQ20rRMSRnN5zVhW6jDFGVZ5jEsGGH5eoFuWdGBqAR4rWzUnXb6ScdkCE\n0w7I4LQDIpx2QITTDohw2gERTjsgg+uwlOeYRLDhxyUqzVVpu6XlfR/DfODsonVnpeuNMcZsQvJO\nDF8HThaR5SLyqIgsB/5vun4j5bUDIrx2QAavHRDhtQMivHZAhNcOiPDaARm8dkCEVxk175+rPi0i\nI4ExwCBgOfBoCOGDasYZY4zpeiUnBhGZGEKYDhBCWEPn1xk2Yk47IMJpB2Rw2gERTjsgwmkHRDjt\ngAinHZDBaQdEOJVRy+1Kur5wQUReqmKLMcaYGlBuYpCi5Z7VCul6XjsgwmsHZPDaARFeOyDCawdE\neO2ACK8dkMFrB0R4lVHLTQyhzLIxxphNTLkXnzcXkcLzLWxRtEwI4cLKZ3UFpx0Q4bQDMjjtgAin\nHRDhtAMinHZAhNMOyOC0AyKcyqjlJoYpwA4Fy9OKlu0ZhDHGbGJKTgwhhK90VUjX89Tebwie2muC\n2uzyWFMeHmvKy1N7XR6Npg05g5sxxphNUDeeGJx2QITTDsjgtAMinHZAhNMOiHDaARFOOyCD0w6I\ncCqjduOJwRhjTEw3nhi8dkCE1w7I4LUDIrx2QITXDojw2gERXjsgg9cOiPAqo67LGdw+QnKCnm0o\neONbCGF2FbqMMcYoyTUxiEgTMB3YHOgDvAFsCSwDhm1ohIj0BX4N7Aq0AicBC4HfADsCS4CjQgiv\nb+hYa7nKbapinHZABqcdEOG0AyKcdkCE0w6IcNoBGZx2QIRTGTXvrqTLgZ+EELYG3kz//QFwTYU6\nrgTuCSGMBEYBTwOTgQdCCLsAs4FzKzSWMcaYEvJODMNJfngX+jFw5oYGiEgfYGwI4UZIjuKaPjM4\nFLgpvdpNwGEbOlZHvrKbqwivHZDBawdEeO2ACK8dEOG1AyK8dkAGrx0Q4VVGzTsxvE6yCwngRRH5\nBNAP+FgFGoYCr4jIjSIyV0SuE5FewLYhhJUAIYQVQH0FxjLGGFNG3onhTuAL6eUbgIeAvwMzKtCw\nGbAHcHUIYQ/gbZLdSFU+gJ+r7OYqwmkHZHDaARFOOyDCaQdEOO2ACKcdkMFpB0Q4lVHznsHtjILL\nPxORR0meLdxbgYYWYFkI4W/p8h0kE8NKEdk2hLBSRAYC0XNBzJgxg1uAIenyVsBurL07fadbtK0p\nvka6vDhdHEp0ufnxZlgCTU1NAMyZMwfovLxvvtHKXmPBvOakIaOHxUlTU0Ppnrblcj3NC/4CDChx\nDZ80tSlzf5XreaS5md4lR2trmljiGgXLFXj86pqbOajMaGuV7qnU45f3+6lij1+q6x6/Mj05v58q\n/fiV66n041f+Guv/+M2ZM4cpU6YAMHjwYOrr6xk3bhwxuf9ctVAIoWJnckt/8C8TkeEhhIXAOODJ\n9ONE4FLgy8DM2O2PPPJIxl17beb2Xfrv3Z3W+PSyo4OhlFxuHN3Y/k0Aa78BOi2nD0jR1iPzf+Ea\n3+kaI3dvhFWl+xpHN5bvyRi/eLlx5BigoegavsPyyN3fggeze0qNX7w8trGxfR9lVt9bI8dEruGp\n1uPXo0wPxL6fCpvWqtjjl/P7qfPjV3ytWnj8fOQWd1ekp9KPX7meSj9+8WcIrsPl9X38mpqaOizP\nnTs3MlYi75+rDgYuAnan6HWFEMLwPNso43TgNhHpCSwCvgL0AG4XkZOA54GjKjCOMcaYMvI+Y5hO\n8iekFwLvVjoihPAPYHTkU+MrPdZarnqbXm9OOyCD0w6IcNoBEU47IMJpB0Q47YAMTjsgwqmMmndi\nGAGMCSG0VjPGGGOMvrx/lXQXa18/2UR47YAIrx2QwWsHRHjtgAivHRDhtQMivHZABq8dEOFVRs37\njOF04M8i8hywsvATIYSTKl5ljDFGTd6J4UbgQ2ABVXiNQYfTDohw2gEZnHZAhNMOiHDaARFOOyDC\naQdkcNoBEU5l1LwTw/7AdiGEN6sZY4wxRl/e1xj+CfSvZkjX89oBEV47IIPXDojw2gERXjsgwmsH\nRHjtgAxeOyDCq4ya9xnDbOA+EbmRzq8x3FDxKmOMMWryTgxNwAvAAUXrA8mxkzZCTjsgwmkHZHDa\nARFOOyDCaQdEOO2ACKcdkMFpB0Q4lVHzHitpv2qHGGOMqQ2ZrzGIiBRcrsv66JrMavDaARFeOyCD\n1w6I8NoBEV47IMJrB0R47YAMXjsgwquMWuoZQ+E5GNbQ+bDXkq7rgTHGmE1GqYnh3wouFx/DbxPg\ntAMinHZABqcdEOG0AyKcdkCE0w6IcNoBGZx2QIRTGTVzYgghLCu4/HzX5BhjjNGW6zUCEdlaRC4R\nkXtE5OHCj2oHVo/XDojw2gEZvHZAhNcOiPDaARFeOyDCawdk8NoBEV5l1Lx/rjoF2By4HXinejnG\nGGO05Z0YPgMMCCG8X82YruW0AyKcdkAGpx0Q4bQDIpx2QITTDohw2gEZnHZAhFMZdV0OiVF8vkBj\njDGboLwTw2zgDyJynoicVPhRzbjq8toBEV47IIPXDojw2gERXjsgwmsHRHjtgAxeOyDCq4yad1fS\nWKAF+GzR+o34kBjGGGNiuvEhMZx2QITTDsjgtAMinHZAhNMOiHDaARFOOyCD0w6IcCqjZk4MIiIh\nhJBeztzlZOeBNsaYTUup1xheL7i8Bvig6KNt3UbKawdEeO2ADF47IMJrB0R47YAIrx0Q4bUDMnjt\ngAivMmo3PiSGMcaYmG58SAynHRDhtAMyOO2ACKcdEOG0AyKcdkCE0w7I4LQDIpzKqCVffBaRsn9x\nFELYiP9k1RhjTLFy72M4ERgDvExyBrfYx0bKawdEeO2ADF47IMJrB0R47YAIrx0Q4bUDMnjtgAiv\nMmq5P1f9IvCl9OO3wM0hhL9UvcoYY4yaks8YQggzQwhHAJ8gOSzGz0XkGRH5fyKyVZcUVo3TDohw\n2gEZnHZAhNMOiHDaARFOOyDCaQdkcNoBEU5l1FyHxAghrA4h/CfwOWAmcDGwexW7jDHGKCk7MaTn\ndv68iEwDngG2BsaFEB6qel1Vee2ACK8dkMFrB0R47YAIrx0Q4bUDIrx2QAavHRDhVUYt91dJlwFH\nAc3AzcCJIYT3uiLMGGOMjnIvPp8JPAdsCZwKnCoiHa4QQtinOmnV5rQDIpx2QAanHRDhtAMinHZA\nhNMOiHDaARmcdkCEUxm13MTwlS6pMMYYUzNKTgwhhJu6KqTreWrvNwRP7TVBbXZ5rCkPjzXl5am9\nLo9GU94T9RhjjOkmuvHE4LQDIpx2QAanHRDhtAMinHZAhNMOiHDaARmcdkCEUxm1ZiaG9M9i54rI\nrHS5n4jcl76h7l4R6avdaIwx3UHNTAzAt4CnCpYnAw+EEHYhOef0uZUdzld2cxXhtQMyeO2ACK8d\nEOG1AyK8dkCE1w7I4LUDIrzKqKXO4Pb9PBsIIVy4oREi0gB8AfgP4Kx09aHAvunlm0juockbOpYx\nxpjSSv1V0g5dVgGXA98GCncXbRtCWAkQQlghIvWVHdJVdnMV4bQDMjjtgAinHRDhtAMinHZAhNMO\nyOC0AyKcyqilTtTTJe9hEJGDgJUhhPki4kpcNXRFjzHGdHfl3uDWgYhsCWwDtL/9OYSwaAMb9gYm\niMgXgC2ALUXkFmCFiGwbQlgpIgOBl2I3njFjBrcAQ9LlrYDdWDvP+k63KFzjCpbTWyxOF9tOZlq0\n3Px4MyyBpqYmAObMmQN0Xm7bB1a09U7LHdcUt8GCec1JQ0YPi5OmpobSPW3L5XqaF/wFGFB0jfnA\nGe3LC+Y1r80sc3+V63mkuZneJXp8e9PEomu0XavoFhV4/OqamzmoRE9Hhdcobqvc45f3+6nz43cF\nxf8jcj1+qeo8fh2/nzrYwO+nSj9+5Xoq/fh1XlO5x2/OnDlMmTIFgMGDB1NfX8+4ceOIyTUxiMgn\ngNuAUSS/uQtrf4PvkWcbWUII5wHnpePsC5wdQjhBRH5CcqKgS4EvkxzVtZMjjzyScddem7l9l/57\nd6c1vmg5VXx266LlxtGN7d8EsPYboNNy+oAUbT3yxLD0NUbu3girSvc1jm4s35NrNGgcOQZoKHmL\nkbu/BQ9m95Qav3h5bGMjfcr0vTVyTOQaPn6LCjx+hd/QsR6IfT/Flyv2+OX8fur8+O1WdK1aefyK\nl++uSE+lH79yPZV+/Dqvqdzj19TU1GF57ty5nUZvk/evkq4BHiI5suobQD/gVyQ/sKvlx8BnReQZ\nYFy6XEGuspurCKcdkMFpB0Q47YAIpx0Q4bQDIpx2QAanHRDhVEbNuytpFPDZEMIHIiIhhNdF5NvA\nE8CtlYoJIfwR+GN6eRUwvlLbNsYYk0/eZwzvAT3Ty6+IyOD0tv2rUtUlvHZAhNcOyOC1AyK8dkCE\n1w6I8NoBEV47IIPXDojwKqPmnRgeITkvA8AM4Pckv9nPrkaUMcYYPbl2JYUQjipYPA94EvgYycl7\nNlJOOyDCaQdkcNoBEU47IMJpB0Q47YAIpx2QwWkHRDiVUdf1z1WFZPfRrSEEe1+BMcZsgnLtShKR\nrdL3FrwLrATeFZFbRGTrqtZVldcOiPDaARm8dkCE1w6I8NoBEV47IMJrB2Tw2gERXmXUvK8x3Ejy\n5rPdSXYh7Q5sDtxQpS5jjDFK8u5K2h8YGEJ4N11eICInAsurUtUlnHZAhNMOyOC0AyKcdkCE0w6I\ncNoBEU47IIPTDohwKqPmfcbwNGuPOtFmMPBMRWuMMcaoy5wYROSktg+SN2DfJyKXiMjXReQS4D7g\nga4KrTyvHRDhtQMyeO2ACK8dEOG1AyK8dkCE1w7I4LUDIrzKqKV2JZ1QtPwsMCb9AHiu4LIxxphN\nRKnDbu/XlSFdz2kHRDjtgAxOOyDCaQdEOO2ACKcdEOG0AzI47YAIpzJq7vcxiEg/4BBge+AF4O70\neEbGGGM2IXnfxzCGZNfRKcAnga8Bz6brN1JeOyDCawdk8NoBEV47IMJrB0R47YAIrx2QwWsHRHiV\nUfM+Y7gCODWEMK1thYgcDfwCGF2NMGOMMTry/rnqcOD2onUzgJ0rm9OVnHZAhNMOyOC0AyKcdkCE\n0w6IcNoBEU47IIPTDohwKqPmnRj+FzimaN1Ekt1LxhhjNiF5J4YzgKtE5K8i8hsReZTkrG6nVy+t\n2rx2QITXDsjgtQMivHZAhNcOiPDaARFeOyCD1w6I8Cqj5j3s9p9FZCfgIGA74C7gHvurJGOM2fTk\n/nPVEMJqKngaz67wbn09y848k7DVVmw9eBW/+tXbBZ8dRPERPQZ/ZGt+NeZXJbfZ/4P+LF9e/hBR\n0r8/PX5VelsbR1PnrjxNubv69GFYfT1bvPRS2e115Nbx+l3BaQdEOO2ACKcdkMFpB0Q4lVEzJwYR\neQQoe86FEMI+FS2qkHfr61l01VUM3H9/6urqqAdGNpa7VT2NjKxMwHbbwcjS27ImaB04kEVXXcWw\n005bj8nBGFMNpZ4x/LrLKqpg2Zlntk8KpnbV1dUxcP/9WXbWWQyfPHkdbumpvd/wPNaUh6f2mqA2\nuzwaTaUOiXFTV4ZUWthqK5sUNhJ1dXWEvn21M4wxqZKvMYhIX2BUCOHhdPm8otv8Mn3tofbYpLBx\nWefHy1WjYgM57YAIpx0Q4bQDMjjtgAinMmq5F5/PAdYAD6fL5wH/k14eCvQELqhOmjHGGA3lfk37\nInB9wfIHIYQTQggnkLzh7dCqlZmqGDx4MEuXLtXOqACvHRDhtQMivHZAhNcOyOC1AyK8yqjlnjFs\nF0JoKVi+ru1CCKFFRBqqk1UddUuWIC0t5a+4nkJDA61DhuS+/pQpU7jmmmtYsmQJffr04aCDDuKC\nCy6gT58+VWvcNCYFY0w1lX0fg4gMCCG8DBBC+G7h+mqGVYO0tNBnwoSqbf+NWbMg58Rw1VVXcfXV\nV3PNNdewzz778OKLL3L22WdzxBFH8Ic//IEePXpUrXPT4LQDIpx2QITTDohw2gEZnHZAhFMZtdyu\npD8BX8n43FeAv1Q2p3t48803+clPfsKll17KfvvtR48ePWhoaOCGG25g6dKlTJ8+nW984xtccskl\n7bf505/+xK677tq+fOWVV/KpT32KwYMH85nPfIbf/e537Z9bvHgxhxxyCEOGDGH48OGcfPLJ7Z/r\n378/S5Ys6ZKv0xizcSr3jOF7wEMish1wJ7CC5G2whwMnAftXN2/T9Nhjj/H+++9z8MEHd1jfu3dv\nxo8fz0MPPUTPnj073U5E2i8PHTqU3//+99TX1/Pb3/6WU045hb///e/U19dzySWXsP/++3PXXXfx\nr3/9i3nz5kW3sXHz1N5veB5rysNTe01Qm10ejaaSzxhCCH8DPgd8CpgNPA08CPw7cGAI4fGqF26C\nVq1aRf/+/aPvs9h222159dVXy25jwoQJ1NfXA3DYYYcxbNgw5s6dC0DPnj1ZtmwZy5cv5yMf+Qh7\n7rln++1CKPtmdmNMN1f2j8dDCH8OIYwF+gA7AH1CCHuHEP5U9bpN1NZbb82rr75Ka2trp8+tXLmS\n/v37l93GtGnT2HfffRk6dChDhw7l6aefbp9Qvve979Ha2spnP/tZ9t57b2677baKfw36nHZAhNMO\niHDaARFOOyCD0w6IcCqj5n5XUQjhnRDCCyGEd6oZ1B2MHj2azTffnLvuuqvD+rfeeosHHniAvffe\nm169evHFIb6AAAAbeUlEQVTuu++2f27FihXtl1taWjjzzDP56U9/yuLFi1m8eDEjRoxofzYwYMAA\nrrjiCp588kkuu+wyvv3tb9vrCsaY3OztwQr69OnDOeecw+TJk3nwwQdZs2YNS5cu5atf/SrbbLMN\nEydOZNddd+X+++/ntddeY+XKlfyq4Kiob7/9NnV1dfTv35/W1lZuu+02FixY0P75mTNnth/VtG/f\nvtTV1W2Chwfx2gERXjsgwmsHRHjtgAxeOyDCq4ya+7Dbm4LQ0JD8SWkVt5/X6aefTv/+/bnwwgtZ\nvHgx77//Pk1NTdx5551sscUWHH300fzxj39k1KhR7Ljjjhx77LFcffXVAOyyyy6ceuqpHHDAAfTo\n0YOjjz6avfbaq33b8+bN47zzzuPNN9+kvr6eH/3oRwwePBjYlF58NsZUS7eaGFqHDMn9PoOucNxx\nx3HccccBMHXqVH70ox+xZs0aADbffHOuv/76Dtc/5ZRT2i+ff/75nH/++dHtXnzxxVx88cXRz73y\nyisVKK8FTjsgwmkHRDjtgAinHZDBaQdEOJVR12liEJE+wLlAI7AI+HEIofwZYkxZkyZNYrPNNuOx\nxx7ji1/8onaOMaYbW9cdz1cDbwG/AN4GZmxogIg0iMhsEXlSRJpF5PR0fT8RuU9EnhGRe9MjvW7S\nJk6caJNCbl47IMJrB0R47YAIrx2QwWsHRHiVUUtODCJyhYh8rGDVYJJnCfcBPwRGVKBhDXBWCOHf\ngDHAN0RkBDAZeCCEsAvJeyjOrcBYxhhjyij3jOEx4I8icnS6fAcwT0RuBeYCG3wynxDCihDC/PTy\nW8ACoIHkyK1t278JOGxDxzKbEqcdEOG0AyKcdkCE0w7I4LQDIpzKqOXe+TyF5LAXTSJyL3AvyeG2\nZwHHhxDOrGSMiAwBdgP+CmwbQliZdqwgOR2xMcaYKiv74nMI4XXgmyLyKeAG4I/A90MI71UyJN1l\nNQP4VgjhLREpPnZD9FgOM2bM4BZgSLq8FcnMMqiScabL+PRfV7DcvOAvwMSia7Rdq+gWi9PFoUSX\nmx9vhiXQ1NQEwJw5c4COy3XNzRxUoie7uLgNFsxrThoyelicNDU1ZPcA7Fump205ua8GFKy5guR/\nxNpbLJjXvDYz6/5KZfW0LT/S3EzvEj2e2OM3Hzij6BZleobm66n041eup9KPX+c1lXv85syZw5Qp\nU4DkvCz19fWMGzeOmHKn9tyOZF//MOBJkt07RwN/EZGLQggVeVOAiGxGMincEkKYma5eKSLbhhBW\nishA4KXYbY888kjGXXttp/XPVCLMdDkXWX5r5JjINXz8FkMpudw4urH9PzGs/Q9cuFx4wPNYD8Dd\nZa+RGLl7I6wq3dc4urFkDwDpf/DSo0HjyDEke2Lb7FZ0LcfI3d9KjniW0RMdP2N5bGMjhWcPifXF\nH7/i5bsr0lPpx69cT6Ufv85rKvf4NTU1dVhuO7ZaTLnXGKYD7wG/BAT4RQjhauBA4CgRuavUjdfB\nDcBTIYQrC9bNAk5ML38ZmFl8I9OdOe2ACKcdEOG0AyKcdkAGpx0Q4VRGLbcraSTgQggfiMgfSfb9\nk+77P15E9tvQABHZGzgOaBaReSS7jM4DLgVuF5GTgOeBozZ0LJNYtmwZu+22Gy+//PImeKgMY8yG\nKjcx3Aw8ICJzgLHAfxd+MoTw0IYGpEdpzTpd2fgN3X6hJUvqaGmp3iEhGhoCQ4Z0PmJqzKhRo3jl\nlVfo0aMHvXv3Zv/99+enP/0pvXr1qlpfoY3/0Bie2vsNz2NNeXhqrwlqs8uj0VRyYgghnCEio0n2\nXE0JITzZNVnV0dIiTJhQvfMpz5r1Ru4jbogI06ZNY+zYsaxYsYIjjjiCyy67jAsuuKBqfcYYk0ee\n8zE8HkK4fWOfFGpR22GyBw4cyPjx41mwYAFTp05ljz32YPDgweyxxx7ccccd7de/9dZb2Wuvvdhp\np52YOHEiLS0tQLJrqO1Iq20mTJjArbfeCkBraysXXHABH//4x/nUpz7Ffffd16FjxYoVHHfccey0\n006MHj2am2++udpfegU47YAIpx0Q4bQDIpx2QAanHRDhVEbtVgfRq1UtLS3cf//9HHDAAUyePJmH\nHnqIYcOG8dJLL7F69WoA7rnnHq688kqmTp3KsGHDuOKKKzj55JP5wx/+AJTeNXTTTTdx//338/DD\nD9OrVy++9KUvdfj8V7/6VXbddVeefvppnnnmGQ4//HCGDRvW6S8sjDHdg73yqOj4449n2LBhHHzw\nwYwdO5bvfOc79OjRg6eeeor33nuP+vp6dtllFwD++7//mzPOOIOdd96Zuro6zjjjDJ544on2Zw2l\nzJw5k1NOOYVBgwbRt29fzjjjjPbPtbS08Pjjj3PRRRfRs2dPdt11V0444QSmTZtWta+7Mrx2QITX\nDojw2gERXjsgg9cOiPAqo9rEoOi2225j0aJFzJ8/n0svvZRevXpx/fXXc8MNNzBy5EgmTZrEs88+\nCyS7i84991yGDRvGsGHD2GmnnRARXnzxxbLjvPjii2y//fbtyzvssEP75ZUrV9KvX78OL3rvsMMO\nubZrjNk02cSgqO01hkL77bcfd955J08//TQ777xz+2/322+/PZdffjmLFi1i0aJFLF68mGXLljF6\n9Oj2H+rvvLP2rKsrV65svzxw4EBeeOGF9uVly5Z1+Nzq1at5++2329e1tLQwaFCtv3fcaQdEOO2A\nCKcdEOG0AzI47YAIpzKqTQw15OWXX+b3v/8977zzDj179qR3797t7zP4yle+ws9//nOefvppAN54\n4w1mzkze89e/f38GDRrE9OnTaW1t5dZbb+1wjufDDjuM6667juXLl/Paa6/xi1/8ov1z22+/PZ/+\n9Kf5wQ9+wPvvv8+TTz7JrbfeytFHH40xpnvqVi8+NzQEZs16o6rbzyv2YnFrayvXXHMNp556KiJC\nY2MjP/vZzwA46KCDeOeddzj55JNpaWmhT58+OOc49NBDAbjiiis455xz+MEPfsDxxx/Pnnvu2b7d\nL33pSzz33HPss88+9OnTh9NOO41HHnmk/fP/9V//xVlnncUnPvEJ+vXrx7nnnsvYsWPX927oIp7a\n+w3PY015eGqvCWqzy1Nz72PY1AwZ0lozZ/acN29ep3Xbbrstd92VfZSRiRMnMnHixOjnxo0bF90m\nQI8ePfjhD3/ID3/4w/Z1J510UvvlQYMGMXXq1LzpxphNnO1KMhsppx0Q4bQDIpx2QITTDsjgtAMi\nnMqoNjEYY4zpwCYGs5Hy2gERXjsgwmsHRHjtgAxeOyDCq4xqE4MxxpgObGIwGymnHRDhtAMinHZA\nhNMOyOC0AyKcyqg2MRhjjOnAJgazkfLaARFeOyDCawdEeO2ADF47IMKrjGoTgzHGmA5sYjDr5aij\njuI3v/mNYoFTHDuL0w6IcNoBEU47IIPTDohwKqN2q3c+L3l9CS1vlj9M9fpq2LKBIX2H5LruqFGj\neO+995g/fz5bbLEFALfccgvTp09n1qxZVWuE7HM+f+Mb32D77bfnvPPOK7uN22+/vf3y1KlTueWW\nW7jnnnuq0muM6VrdamJoebOFCXdOqNr2Zx0+K/fEICK0trZy7bXXcuaZZ3ZYv65CCOt8u0qe83l9\nxt9wntr7Dc9jTXl4aq8JarPLo9Fku5IUffOb3+Tqq6/mjTc6H9jv0UcfZfz48QwdOpTx48fz2GOP\ntX9uwoQJ/Md//Aef//znaWho4Pnnn29fd+CBBzJ48GCOO+44Vq9ezde+9jV23HFHxo8fn+ukPm2m\nTp3KF77wBS688EKGDRvGHnvswQMPPNCh4dZbb2XhwoWcc845PP744wwePJhhw4Yxb948RowY0eGw\n4nfddRf77LPPet5TxpiuZBODot122429996bX/7ylx3Wv/baa0yaNIlTTjmF5557jq9//escc8wx\nvPbaa+3Xuf3227nyyitZunQpDQ0NAPz2t7/luuuu48knn2TRokUceOCBHH/88SxevJjhw4dz6aWX\nrlPf3LlzGT58OM899xzf/OY3+da3vtXpOsOHD+eyyy5j9OjRLF26lEWLFrH77ruz9dZbM3v27Pbr\nTZ8+nUmTJq3T+KW5Cm6rUpx2QITTDohw2gEZnHZAhFMZ1SYGZZMnT+bXv/41q1atal933333sdNO\nO3HkkUdSV1fHEUccwcc//vH28zsDTJo0ieHDh1NXV8dmmyV7BI899lgGDx7Mlltuyfjx4xkyZAhj\nx46lrq6OQw89lObm5nVq22GHHTj++OMREY455hhWrFjByy+/nOu2xxxzTPvrEKtXr2b27NkcccQR\n6zS+MUaHTQzKRo4cyec+9zkuv/xyINlfv2LFig6n34TOp9ssPFVnmwEDBrRf/uhHP9ppue0sbT16\n9ADggw8+6HD7NWvW0LNnz/bl+vr69sttL5AXnumtlIkTJ3Lvvffy7rvv8tvf/pYxY8Z02N6G8xXc\nVqV47YAIrx0Q4bUDMnjtgAivMqpNDDXgu9/9LjfffDMvvvgiIsKgQYNYunRph+sUn25zQ17sHThw\nID179uw0xvPPP99pQsoj1jJo0CBGjx7NXXfdxe23325nhDNmI2ITQw0YOnQoX/ziF7nuuusAGD9+\nPIsWLeKOO+7gww8/5M4772ThwoUceOCBFRmvrq6OQw45hB/+8IesXr2aNWvWcMcdd7Bw4ULGjx+/\nztsbMGAAy5cv7/QM5Oijj+YXv/gFCxYs4OCDD65I+1quwturBKcdEOG0AyKcdkAGpx0Q4VRG7VZ/\nrtqwZQOzDq/eewQatmzIfd3i37K//e1vM336dESEfv36MXXqVM4991zOOecchg0bxrRp09hqq62i\nt81aV8rPfvYzLrzwQpqamnjvvffYZZdd+M1vfsM222yTq7nw8j777MOIESMYMWIEPXr0YOHChUBy\nOtKzzz6bQw45hI9+9KPr1GeM0dOtJoYhfYfkfp9BtRWfhnP77bfnhRdeaF/ec889O/xVT6GZM2eW\nXXf++ed3WN533315/PHH25f79OnDFVdckdk3adKkTn9F9Morr0TH69mzZ/TUoFtssQXbbLMNRx11\nVOY4689Te7/heawpD0/tNUFtdnnsfQxmkzJr1izq6urs/QvGbGS61TMG03UmTJjAwoULufbaa6s0\ngqvSdjeE0w6IcNoBEU47IIPTDohwKqPaxGCqotrHezLGVM+muyuptVW7wKyLdX68fDUqNpDXDojw\n2gERXjsgg9cOiPAqo26yE4O89hqtNjlsFFpbW5HXX9fOMMakNtmJYYfLL2fF7Nk2OdS41tZWVsye\nzQ4///k63tJVI2cDOe2ACKcdEOG0AzI47YAIpzLqJvsawxYvvcSw005j2VlnEfr25dVPfIYnnt68\n5G1G7PYqT7/2RMnr7LrNrvTdvG/Z8eX11+nxROltvTpir022KXfXK68w7LTT2OKll8puzxjTNWp+\nYhCRA4ErSJ7dXB9CyH2I0C1eeonhkycDcPcty/ja1wrfgOYpno1veeRuvvaXr5Xc5qzDZzFyu5Fl\nx+6xaBF9vlZ6WxtHU+euPE15u3osWrSek0LHptrgsaY8PLXXBLXZ5bH3MRQRkTrgKuBzwL8Bk0Rk\nRGW2Pr8ym6moWmyC2uyypnysKb9a7NJpqumJAfg08L8hhOdDCB8A04BDK7Pp18pfpcvVYhPUZpc1\n5WNN+dVil05TrU8M2wPLCpZb0nXGGGOqpNYnhipaoh0QsUQ7IMMS7YCIJdoBEUu0AyKWaAdELNEO\nyLBEOyBiicqoUnhe3lojInsBF4cQDkyXJwOh8AXoyy67LPzjH/9ov82oUaPYbbfdym57/vz5ua7X\nlWqxCWqzy5rysab8arGrkk3z58+n+Gfl2WefHT0sc61PDD2AZ4BxwIvAY8CkEMIC1TBjjNmE1fSf\nq4YQPhSR04D7WPvnqjYpGGNMFdX0MwZjjDFdr9u9+CwiDSIyW0SeFJFmETm9Bpo2F5FHRWRe2nSR\ndlMbEakTkbkiUjOHSxWRJSLyj/T+eky7B0BE+orIdBFZkH5v7ancMzy9f+am/75eI9/rZ4rIEyLy\nTxG5TUQ+UgNN30r/36n+PBCR60VkpYj8s2BdPxG5T0SeEZF7RaT8IQ4qoNtNDMAa4KwQwr8BY4Bv\nVO5Nc+snhPA+sF8IYXdgN+DzIvJpzaYC3wKe0o4o0gq4EMLuIYRauZ+uBO4JIYwERgGquzxDCAvT\n+2cP4FPA28D/aDaJyHbAN4E9QgifJNmVfYxy078BXwX+neT/3sEiMkwp50aSN/MWmgw8EELYBZgN\nnNsVId1uYgghrAghzE8vv0XyH1j9vREhhHfSi5uT/IdR38cnIg3AF4Bfa7cUEWroe1dE+gBjQwg3\nAoQQ1oQQ3lDOKjQeeC6EsKzsNauvB9BbRDYDegHLlXtGAo+GEN4PIXwIPAwcrhESQpgDrC5afShw\nU3r5JuCwrmipmf9cGkRkCMlvCY/qlrTvspkHrADuDyE8Xu42XeBy4NvUwCRVJAD3i8jjIvJ/tGOA\nocArInJjuuvmOhHZQjuqwNFA55Nyd7EQwnLgMmAp8ALwWgjhAd0qngDGprtsepH8IrSDclOh+hDC\nSkh+qQXqu2LQbjsxiMjHgBnAt9JnDqpCCK3prqQGYE8R+YRmj4gcBKxMn11J+lEr9k53kXyBZFdg\nk3LPZsAewNVp1zskuwDUiUhPYAIwvQZatiL5DXhHYDvgYyJyrGZTCOFp4FLgfuAeYB7woWZTGV3y\nS1q3nBjSp7EzgFtCCDO1ewqluyAeAg5UTtkbmCAii0h+29xPRG5WbgIghPBi+u/LJPvNtV9naAGW\nhRD+li7PIJkoasHngb+n95W28cCiEMKqdLfNncBnlJsIIdwYQvj3EIIjOTjRQuWkQitFZFsAERkI\ndMnx6bvlxADcADwVQrhSOwRARLZp+2uDdBfEZ4GnNZtCCOeFEAaHEIaRvEA4O4TwJc0mABHplT7b\nQ0R6AweQ7A5Qkz7VXyYiw9NV46idF+wnUQO7kVJLgb1E5KMiIiT3k/r7kkRkQPrvYOCLwBTNHDo+\nO58FnJhe/jLQJb/I1vQb3KpBRPYGjgOa0336ATgvhPAHxaxBwE3pYcbrgN+EEO5R7Kll2wL/IyKB\n5Pv3thDCfcpNAKcDt6W7bhYBX1HuId1nPh74v9otACGEx0RkBsnumg/Sf6/TrQLgDhHZmqTpVK0/\nHBCRKSQnX+gvIkuBi4AfA9NF5CTgeeCoLmmxN7gZY4wp1F13JRljjMlgE4MxxpgObGIwxhjTgU0M\nxhhjOrCJwRhjTAc2MRhjjOnAJgaz0UmPSfR95fFXichftRqMqSabGMwGS8+PsLLwwHEi8lUReUiz\nqxrS4zKNA7YLIeyVcZ2BIvJrEVmengfhKRG5qMYOrFc16ddaE4dPMevHJgZTCYHke+mMyPqalr7b\nfF0MAZaEEN7L2F4/4C8kh0/fM4TQl+QQJ32BnTYg1ZguYxODqZSfAmen5yboQER2FJHWwh/CIvJQ\n+jZ/ROTLIjJHRH4uIqtF5FkRGZOuXyoiK0Sk+DhNA9IzW72RbmtwwbZHpJ97NT2j2sSCz90oIteI\nyO9E5E2SQxAU9w4SkZnp7ReKyMnp+pOA/wLGpOPGzrR3NvBGCOGEtvMfhBBeCCGcFUJ4It3OZ0Tk\nsfRrfVRExhTdLz8QkT+JyJtpx9Yicmv67OPRoq+1VUS+KSLPichLIvKTgs+JiPy/9BndChH577bH\np+Ax+ZKIPJ/e9ryi205OH4uXRWSaJEdHLXlbEfkccB5wdNo/L11/Ytr4RvrvpMh9Z2pFCME+7GOD\nPoDFwP4kRxX9QbruqyQH3oPkMMsfAnUFt3kIOCm9/GXgX8CXSA4g9gOS48L8EuhJ8hv3G0Cv9Po3\nAq+THAG2J3AF8Ej6uV4kB2tr29Yo4GVgRMFtVwN7pcsfiXw9DxeMPYrkiJauoPXhEvfFX4CLSny+\nH7AKOJbkF7Nj0uV+BffLQpJnJlsCT5IcUHG/9Po3AdcXbK8VeJDkGUkD8EzB/XpSuq0d0/vlDuDm\ngsekFfgV8BHgk8B7wC7p578F/JnkOF49gf8EpuS87UVt4xQ8Jq8DO6fL2wIjtb9v7SP7w54xmEq6\nCDhNRPqvx20XhxBuDslPjt+Q/JD7XgjhgxDC/SQTx84F1/9dCOFPIYQPgPNJjtq5PXBw4bZCCP8g\n+YE4seC2M0MIfwUIIfyrMEKSs9aNAb6bjv0PkjPY5T2ybH/gxRKfPwhYGEKYEpJzcEwj+cF/SMF1\nbgwhLAkhvAn8nuTsaw+FEFpJzquwe9E2fxxCeD2E0EIySbb9Nn4s8PMQwvMhOUPgucAxBc/cAnBx\nCOFfIYR/Av8gmQgBvgacH0J4Mb2Pvw8cmfO2MR8CjSLy0RDCyhCC+lFVTTabGEzFhBCeBO5m/c5L\nu7Lg8rvp9l4pWvexguX201SGEN4meRawHclvs3ulfzW0SkRWk/yA3DZ224jtgFVh7alWIXn2kvf0\nr6+S/JZdavvPF60r3n7xfVG8XHg/QHI+iMJtbZcx1vMkR6QtvC8Kt/1OwbZ3JDmK7SoRWUVyGPEP\nct62g/S+PBr4OvCiiNwlIrvErmtqg00MptIuBv4PHX/QvZ3+26tg3cANHKf99IuSnJ+hH8n5g5cB\nPoSwdfrRL4TQJ4RwWsFtS70ovhzYWpJzPbQZTHIqyjweIDmmf6ntDylaty7bjyk8FeWOrD2P8vJ0\nufBzH9DxB3qWpcDni+7H3iE9SVIZne7fEML9IYQDSB73Z0heqzE1yiYGU1EhhOdIdgWdXrDuFZIf\nfMdLcm7rkyj/FzrlTiX6hfRF3I+QvCbx1xDCCyTPWIaLyPEispmI9BSRf8/7G2q6O+bPwI9EZHMR\n+STJ6yW35Lk98HOgj4jc1PYisYhsLyKXiciuJKeP/LiIHCMiPUTkaJIT0t+Vc/sx3xaRrURkB5L7\nfVq6fipwpogMSSfP/wCmpbukoPR9/CvgkoKvYYCITCj4fKnbrgSGiIikt60XkQmSnB/iA+Atavv0\nmd2eTQymEop/Q/w+ybODwvX/B/gO8ArJD8I/reM2Q9HlKSTPTl4l2ed+PEBIzt99AMmLusvTjx+T\n/PloXpOAoelt7wAuCCHkek9GCGE1yekqPwAeFZHXSc4n/BrwbAhhFcnrIOeQ3BfnAAeltyv+OvOa\nCfwdmEsywdyQrr+BZEJ7GHiOZHfP6QW3K3UfX5lu9770a/gzHU+hWuq200kmjldF5G/p5bNIfjl4\nBdiHZLeSqVF2oh5jNmIi0kry1z6LtFvMpsOeMRhjjOnAJgZjNm72lN9UnO1KMsYY04E9YzDGGNOB\nTQzGGGM6sInBGGNMBzYxGGOM6cAmBmOMMR3YxGCMMaaD/w9zdN1rsLtswQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7faa146ea3c8>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "ind = np.arange(2, num_comps+1)  # the x locations for the groups\n",
    "width = 0.5       # the width of the bars\n",
    "fig, ax = plt.subplots()\n",
    "rects1 = ax.bar(2*ind, 100*np.invert(np.isnan(comp_iters_quasi)).sum(1)/comp_iters_quasi.shape[-1],\n",
    "                width, color='r', label='Quasi')\n",
    "rects2 = ax.bar(2*ind + width, 100*np.invert(np.isnan(comp_iters_pseudo)).sum(1)/comp_iters_pseudo.shape[-1],\n",
    "                width, color='b', label='Pseudo')\n",
    "rects3 = ax.bar(2*ind + 2*width, 100*np.invert(np.isnan(comp_iters_normunity)).sum(1)/comp_iters_normunity.shape[-1],\n",
    "                width, color='g', label='NormUnity')\n",
    "# add some text for labels, title and axes ticks\n",
    "ax.set_xticks(2*ind + 1.5*width)\n",
    "ax.set_xticklabels(tuple(ind))\n",
    "ax.set_xlabel('Number of Components')\n",
    "ax.set_ylabel('% Global Minima Found')\n",
    "plt.legend(loc=3)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.5.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}