{ "cells": [ { "cell_type": "markdown", "id": "aa7f6d5e", "metadata": {}, "source": [ "3D TV-Regularized Sparse-View CT Reconstruction (Proximal ADMM Solver)\n", "======================================================================\n", "\n", "This example demonstrates solution of a sparse-view, 3D CT\n", "reconstruction problem with isotropic total variation (TV)\n", "regularization\n", "\n", " $$\\mathrm{argmin}_{\\mathbf{x}} \\; (1/2) \\| \\mathbf{y} - C \\mathbf{x}\n", " \\|_2^2 + \\lambda \\| D \\mathbf{x} \\|_{2,1} \\;,$$\n", "\n", "where $C$ is the X-ray transform (the CT forward projection operator),\n", "$\\mathbf{y}$ is the sinogram, $D$ is a 3D finite difference operator,\n", "and $\\mathbf{x}$ is the reconstructed image.\n", "\n", "In this example the problem is solved via proximal ADMM, while standard\n", "ADMM is used in a [companion example](ct_astra_3d_tv_admm.rst)." ] }, { "cell_type": "code", "execution_count": 1, "id": "cae513b3", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2025-08-19T21:04:19.447332Z", "iopub.status.busy": "2025-08-19T21:04:19.446299Z", "iopub.status.idle": "2025-08-19T21:04:26.183953Z", "shell.execute_reply": "2025-08-19T21:04:26.182002Z" } }, "outputs": [], "source": [ "import numpy as np\n", "\n", "from mpl_toolkits.axes_grid1 import make_axes_locatable\n", "\n", "import scico.numpy as snp\n", "from scico import functional, linop, loss, metric, plot\n", "from scico.examples import create_tangle_phantom\n", "from scico.linop.xray.astra import XRayTransform3D, angle_to_vector\n", "from scico.optimize import ProximalADMM\n", "from scico.util import device_info\n", "plot.config_notebook_plotting()" ] }, { "cell_type": "markdown", "id": "543b45b5", "metadata": {}, "source": [ "Create a ground truth image and projector." ] }, { "cell_type": "code", "execution_count": 2, "id": "7bbbeaa1", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2025-08-19T21:04:26.195968Z", "iopub.status.busy": "2025-08-19T21:04:26.194859Z", "iopub.status.idle": "2025-08-19T21:04:26.771776Z", "shell.execute_reply": "2025-08-19T21:04:26.769986Z" } }, "outputs": [], "source": [ "Nx = 128\n", "Ny = 256\n", "Nz = 64\n", "\n", "tangle = snp.array(create_tangle_phantom(Nx, Ny, Nz))\n", "\n", "n_projection = 10 # number of projections\n", "angles = np.linspace(0, np.pi, n_projection, endpoint=False) # evenly spaced projection angles\n", "det_spacing = [1.0, 1.0]\n", "det_count = [Nz, max(Nx, Ny)]\n", "vectors = angle_to_vector(det_spacing, angles)\n", "\n", "# It would have been more straightforward to use the det_spacing and angles keywords\n", "# in this case (since vectors is just computed directly from these two quantities), but\n", "# the more general form is used here as a demonstration.\n", "C = XRayTransform3D(tangle.shape, det_count=det_count, vectors=vectors) # CT projection operator\n", "y = C @ tangle # sinogram" ] }, { "cell_type": "markdown", "id": "45faa6e0", "metadata": {}, "source": [ "Set up problem and solver. We want to minimize the functional\n", "\n", " $$\\mathrm{argmin}_{\\mathbf{x}} \\; (1/2) \\| \\mathbf{y} - C \\mathbf{x}\n", " \\|_2^2 + \\lambda \\| D \\mathbf{x} \\|_{2,1} \\;,$$\n", "\n", "where $C$ is the X-ray transform and $D$ is a finite difference\n", "operator. This problem can be expressed as\n", "\n", " $$\\mathrm{argmin}_{\\mathbf{x}, \\mathbf{z}} \\; (1/2) \\| \\mathbf{y} -\n", " \\mathbf{z}_0 \\|_2^2 + \\lambda \\| \\mathbf{z}_1 \\|_{2,1} \\;\\;\n", " \\text{such that} \\;\\; \\mathbf{z}_0 = C \\mathbf{x} \\;\\; \\text{and} \\;\\;\n", " \\mathbf{z}_1 = D \\mathbf{x} \\;,$$\n", "\n", "which can be written in the form of a standard ADMM problem\n", "\n", " $$\\mathrm{argmin}_{\\mathbf{x}, \\mathbf{z}} \\; f(\\mathbf{x}) + g(\\mathbf{z})\n", " \\;\\; \\text{such that} \\;\\; A \\mathbf{x} + B \\mathbf{z} = \\mathbf{c}$$\n", "\n", "with\n", "\n", " $$f = 0 \\qquad g = g_0 + g_1$$\n", " $$g_0(\\mathbf{z}_0) = (1/2) \\| \\mathbf{y} - \\mathbf{z}_0 \\|_2^2 \\qquad\n", " g_1(\\mathbf{z}_1) = \\lambda \\| \\mathbf{z}_1 \\|_{2,1}$$\n", " $$A = \\left( \\begin{array}{c} C \\\\ D \\end{array} \\right) \\qquad\n", " B = \\left( \\begin{array}{cc} -I & 0 \\\\ 0 & -I \\end{array} \\right) \\qquad\n", " \\mathbf{c} = \\left( \\begin{array}{c} 0 \\\\ 0 \\end{array} \\right) \\;.$$\n", "\n", "This is a more complex splitting than that used in the\n", "[companion example](ct_astra_3d_tv_admm.rst), but it allows the use of a\n", "proximal ADMM solver in a way that avoids the need for the conjugate\n", "gradient sub-iterations used by the ADMM solver in the\n", "[companion example](ct_astra_3d_tv_admm.rst)." ] }, { "cell_type": "code", "execution_count": 3, "id": "8ed237a6", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2025-08-19T21:04:26.779737Z", "iopub.status.busy": "2025-08-19T21:04:26.779094Z", "iopub.status.idle": "2025-08-19T21:04:34.531034Z", "shell.execute_reply": "2025-08-19T21:04:34.529979Z" } }, "outputs": [], "source": [ "𝛼 = 1e2 # improve problem conditioning by balancing C and D components of A\n", "λ = 2e0 / 𝛼 # ℓ2,1 norm regularization parameter\n", "ρ = 5e-3 # ADMM penalty parameter\n", "maxiter = 1000 # number of ADMM iterations\n", "\n", "f = functional.ZeroFunctional()\n", "g0 = loss.SquaredL2Loss(y=y)\n", "g1 = λ * functional.L21Norm()\n", "g = functional.SeparableFunctional((g0, g1))\n", "D = linop.FiniteDifference(input_shape=tangle.shape, append=0)\n", "\n", "A = linop.VerticalStack((C, 𝛼 * D))\n", "mu, nu = ProximalADMM.estimate_parameters(A)\n", "\n", "solver = ProximalADMM(\n", " f=f,\n", " g=g,\n", " A=A,\n", " B=None,\n", " rho=ρ,\n", " mu=mu,\n", " nu=nu,\n", " maxiter=maxiter,\n", " itstat_options={\"display\": True, \"period\": 50},\n", ")" ] }, { "cell_type": "markdown", "id": "a82bec9a", "metadata": {}, "source": [ "Run the solver." ] }, { "cell_type": "code", "execution_count": 4, "id": "6b205530", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2025-08-19T21:04:34.536763Z", "iopub.status.busy": "2025-08-19T21:04:34.536579Z", "iopub.status.idle": "2025-08-19T21:06:08.551083Z", "shell.execute_reply": "2025-08-19T21:06:08.550169Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Solving on GPU (NVIDIA GeForce RTX 2080 Ti)\n", "\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Iter Time Objective Prml Rsdl Dual Rsdl\n", "-----------------------------------------------\n", " 0 1.19e+00 1.361e+04 3.267e+04 3.267e+04\n", " 1 1.67e+00 5.171e+04 3.165e+04 1.566e+02\r", " 2 1.72e+00 1.117e+05 3.023e+04 1.511e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 3 1.79e+00 1.899e+05 2.845e+04 2.691e+02\r", " 4 1.90e+00 2.864e+05 2.634e+04 9.520e+02\r", " 5 1.95e+00 3.921e+05 2.394e+04 9.383e+02\r", " 6 2.00e+00 4.989e+05 2.131e+04 8.218e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 7 2.06e+00 6.031e+05 1.850e+04 8.527e+02\r", " 8 2.11e+00 7.014e+05 1.558e+04 9.520e+02\r", " 9 2.16e+00 7.885e+05 1.266e+04 9.844e+02\r", " 10 2.21e+00 8.595e+05 9.914e+03 9.809e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 11 2.26e+00 9.114e+05 7.657e+03 9.649e+02\r", " 12 2.31e+00 9.437e+05 6.488e+03 9.623e+02\r", " 13 2.36e+00 9.562e+05 6.936e+03 9.502e+02\r", " 14 2.41e+00 9.499e+05 8.655e+03 9.310e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 15 2.46e+00 9.270e+05 1.094e+04 9.056e+02\r", " 16 2.52e+00 8.903e+05 1.336e+04 8.740e+02\r", " 17 2.57e+00 8.433e+05 1.572e+04 8.371e+02\r", " 18 2.62e+00 7.900e+05 1.791e+04 7.962e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 19 2.67e+00 7.343e+05 1.987e+04 7.526e+02\r", " 20 2.72e+00 6.801e+05 2.155e+04 7.082e+02\r", " 21 2.77e+00 6.311e+05 2.292e+04 6.654e+02\r", " 22 2.82e+00 5.901e+05 2.396e+04 6.268e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 23 2.87e+00 5.596e+05 2.466e+04 5.952e+02\r", " 24 2.92e+00 5.413e+05 2.501e+04 5.736e+02\r", " 25 2.99e+00 5.359e+05 2.501e+04 5.644e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 26 3.09e+00 5.434e+05 2.467e+04 5.665e+02\r", " 27 3.15e+00 5.628e+05 2.399e+04 5.792e+02\r", " 28 3.20e+00 5.927e+05 2.300e+04 6.003e+02\r", " 29 3.25e+00 6.308e+05 2.172e+04 6.344e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 30 3.30e+00 6.746e+05 2.016e+04 6.685e+02\r", " 31 3.35e+00 7.212e+05 1.836e+04 7.049e+02\r", " 32 3.40e+00 7.675e+05 1.636e+04 7.363e+02\r", " 33 3.45e+00 8.107e+05 1.420e+04 7.753e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 34 3.50e+00 8.480e+05 1.192e+04 8.096e+02\r", " 35 3.55e+00 8.772e+05 9.604e+03 8.302e+02\r", " 36 3.60e+00 8.966e+05 7.346e+03 8.501e+02\r", " 37 3.65e+00 9.050e+05 5.374e+03 8.633e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 38 3.70e+00 9.019e+05 4.205e+03 8.643e+02\r", " 39 3.76e+00 8.875e+05 4.501e+03 8.730e+02\r", " 40 3.81e+00 8.624e+05 5.981e+03 8.574e+02\r", " 41 3.86e+00 8.279e+05 7.919e+03 8.383e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 42 3.91e+00 7.855e+05 9.947e+03 8.208e+02\r", " 43 3.96e+00 7.371e+05 1.192e+04 7.755e+02\r", " 44 4.01e+00 6.852e+05 1.375e+04 7.243e+02\r", " 45 4.06e+00 6.320e+05 1.540e+04 6.735e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 46 4.11e+00 5.802e+05 1.684e+04 6.217e+02\r", " 47 4.16e+00 5.319e+05 1.804e+04 5.690e+02\r", " 48 4.22e+00 4.893e+05 1.899e+04 5.185e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 49 4.34e+00 4.540e+05 1.968e+04 4.731e+02\r", " 50 4.39e+00 4.273e+05 2.009e+04 4.357e+02\n", " 51 4.45e+00 4.099e+05 2.024e+04 4.068e+02\r", " 52 4.50e+00 4.020e+05 2.013e+04 3.865e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 53 4.55e+00 4.033e+05 1.976e+04 3.788e+02\r", " 54 4.60e+00 4.130e+05 1.914e+04 3.831e+02\r", " 55 4.65e+00 4.299e+05 1.829e+04 3.977e+02\r", " 56 4.70e+00 4.525e+05 1.724e+04 4.195e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 57 4.75e+00 4.791e+05 1.600e+04 4.462e+02\r", " 58 4.80e+00 5.078e+05 1.462e+04 4.752e+02\r", " 59 4.85e+00 5.366e+05 1.311e+04 5.027e+02\r", " 60 4.90e+00 5.637e+05 1.154e+04 5.285e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 61 4.95e+00 5.876e+05 9.964e+03 5.503e+02\r", " 62 5.00e+00 6.068e+05 8.461e+03 5.684e+02\r", " 63 5.06e+00 6.204e+05 7.159e+03 5.817e+02\r", " 64 5.11e+00 6.277e+05 6.240e+03 5.899e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 65 5.16e+00 6.284e+05 5.899e+03 5.921e+02\r", " 66 5.21e+00 6.228e+05 6.209e+03 5.896e+02\r", " 67 5.26e+00 6.113e+05 7.034e+03 5.807e+02\r", " 68 5.31e+00 5.948e+05 8.154e+03 5.660e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 69 5.36e+00 5.744e+05 9.391e+03 5.483e+02\r", " 70 5.41e+00 5.514e+05 1.063e+04 5.268e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 71 5.49e+00 5.271e+05 1.180e+04 5.028e+02\r", " 72 5.69e+00 5.031e+05 1.284e+04 4.749e+02\r", " 73 5.76e+00 4.806e+05 1.374e+04 4.466e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 74 5.82e+00 4.608e+05 1.447e+04 4.188e+02\r", " 75 5.87e+00 4.448e+05 1.501e+04 3.927e+02\r", " 76 5.92e+00 4.333e+05 1.535e+04 3.687e+02\r", " 77 5.97e+00 4.267e+05 1.549e+04 3.513e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 78 6.02e+00 4.252e+05 1.542e+04 3.409e+02\r", " 79 6.07e+00 4.287e+05 1.516e+04 3.395e+02\r", " 80 6.12e+00 4.367e+05 1.471e+04 3.438e+02\r", " 81 6.17e+00 4.485e+05 1.407e+04 3.561e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 82 6.22e+00 4.633e+05 1.326e+04 3.740e+02\r", " 83 6.28e+00 4.801e+05 1.230e+04 3.946e+02\r", " 84 6.33e+00 4.977e+05 1.121e+04 4.169e+02\r", " 85 6.38e+00 5.149e+05 1.000e+04 4.393e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 86 6.43e+00 5.308e+05 8.718e+03 4.621e+02\r", " 87 6.48e+00 5.443e+05 7.386e+03 4.823e+02\r", " 88 6.53e+00 5.545e+05 6.059e+03 4.995e+02\r", " 89 6.59e+00 5.608e+05 4.819e+03 5.130e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 90 6.64e+00 5.628e+05 3.817e+03 5.224e+02\r", " 91 6.70e+00 5.603e+05 3.310e+03 5.280e+02\r", " 92 6.77e+00 5.532e+05 3.506e+03 5.294e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 93 7.00e+00 5.420e+05 4.271e+03 5.276e+02\r", " 94 7.10e+00 5.270e+05 5.311e+03 5.218e+02\r", " 95 7.18e+00 5.089e+05 6.436e+03 5.115e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 96 7.23e+00 4.885e+05 7.546e+03 4.971e+02\r", " 97 7.29e+00 4.665e+05 8.589e+03 4.799e+02\r", " 98 7.34e+00 4.441e+05 9.531e+03 4.607e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 99 7.42e+00 4.219e+05 1.035e+04 4.407e+02\r", " 100 7.54e+00 4.010e+05 1.103e+04 4.197e+02\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 101 7.65e+00 3.819e+05 1.157e+04 3.989e+02\r", " 102 7.77e+00 3.655e+05 1.195e+04 3.801e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 103 7.88e+00 3.520e+05 1.218e+04 3.637e+02\r", " 104 8.00e+00 3.417e+05 1.225e+04 3.511e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 105 8.13e+00 3.348e+05 1.216e+04 3.433e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 106 8.39e+00 3.312e+05 1.193e+04 3.402e+02\r", " 107 8.58e+00 3.306e+05 1.155e+04 3.415e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 108 8.75e+00 3.325e+05 1.105e+04 3.456e+02\r", " 109 8.94e+00 3.365e+05 1.043e+04 3.517e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 110 9.11e+00 3.420e+05 9.709e+03 3.596e+02\r", " 111 9.28e+00 3.483e+05 8.914e+03 3.692e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 112 9.44e+00 3.547e+05 8.069e+03 3.797e+02\r", " 113 9.58e+00 3.608e+05 7.206e+03 3.903e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 114 9.85e+00 3.659e+05 6.366e+03 4.000e+02\r", " 115 1.00e+01 3.695e+05 5.607e+03 4.078e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 116 1.02e+01 3.714e+05 5.000e+03 4.123e+02\r", " 117 1.04e+01 3.714e+05 4.628e+03 4.129e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 118 1.05e+01 3.694e+05 4.552e+03 4.113e+02\r", " 119 1.07e+01 3.655e+05 4.770e+03 4.054e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 120 1.08e+01 3.599e+05 5.220e+03 3.972e+02\r", " 121 1.10e+01 3.528e+05 5.812e+03 3.860e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 122 1.13e+01 3.447e+05 6.467e+03 3.761e+02\r", " 123 1.14e+01 3.359e+05 7.129e+03 3.598e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 124 1.16e+01 3.268e+05 7.756e+03 3.405e+02\r", " 125 1.18e+01 3.180e+05 8.321e+03 3.225e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 126 1.20e+01 3.099e+05 8.801e+03 3.046e+02\r", " 127 1.22e+01 3.027e+05 9.184e+03 2.839e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 128 1.23e+01 2.970e+05 9.460e+03 2.669e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 129 1.26e+01 2.929e+05 9.622e+03 2.530e+02\r", " 130 1.28e+01 2.905e+05 9.668e+03 2.428e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 131 1.30e+01 2.900e+05 9.597e+03 2.363e+02\r", " 132 1.32e+01 2.911e+05 9.411e+03 2.338e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 133 1.33e+01 2.937e+05 9.114e+03 2.358e+02\r", " 134 1.35e+01 2.977e+05 8.713e+03 2.417e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 135 1.37e+01 3.026e+05 8.217e+03 2.504e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 136 1.40e+01 3.081e+05 7.636e+03 2.607e+02\r", " 137 1.42e+01 3.138e+05 6.984e+03 2.716e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 138 1.44e+01 3.192e+05 6.277e+03 2.832e+02\r", " 139 1.45e+01 3.240e+05 5.534e+03 2.955e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 140 1.47e+01 3.279e+05 4.782e+03 3.072e+02\r", " 141 1.49e+01 3.306e+05 4.058e+03 3.159e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 142 1.51e+01 3.319e+05 3.417e+03 3.246e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 143 1.54e+01 3.316e+05 2.942e+03 3.328e+02\r", " 144 1.56e+01 3.297e+05 2.731e+03 3.383e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 145 1.58e+01 3.262e+05 2.835e+03 3.390e+02\r", " 146 1.60e+01 3.213e+05 3.204e+03 3.375e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 147 1.62e+01 3.152e+05 3.730e+03 3.355e+02\r", " 148 1.64e+01 3.082e+05 4.323e+03 3.299e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 149 1.65e+01 3.005e+05 4.925e+03 3.237e+02\r", " 150 1.68e+01 2.925e+05 5.500e+03 3.164e+02\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 151 1.69e+01 2.846e+05 6.025e+03 3.043e+02\r", " 152 1.70e+01 2.769e+05 6.484e+03 2.924e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 153 1.71e+01 2.699e+05 6.867e+03 2.815e+02\r", " 154 1.73e+01 2.637e+05 7.166e+03 2.722e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 155 1.74e+01 2.586e+05 7.378e+03 2.640e+02\r", " 156 1.76e+01 2.546e+05 7.500e+03 2.530e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 157 1.77e+01 2.519e+05 7.533e+03 2.410e+02\r", " 158 1.78e+01 2.504e+05 7.478e+03 2.323e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 159 1.79e+01 2.500e+05 7.340e+03 2.273e+02\r", " 160 1.82e+01 2.508e+05 7.124e+03 2.255e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 161 1.84e+01 2.524e+05 6.838e+03 2.260e+02\r", " 162 1.85e+01 2.548e+05 6.491e+03 2.290e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 163 1.87e+01 2.577e+05 6.096e+03 2.320e+02\r", " 164 1.89e+01 2.609e+05 5.666e+03 2.364e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 165 1.90e+01 2.642e+05 5.220e+03 2.415e+02\r", " 166 1.92e+01 2.672e+05 4.779e+03 2.467e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 167 1.94e+01 2.700e+05 4.368e+03 2.516e+02\r", " 168 1.96e+01 2.723e+05 4.019e+03 2.561e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 169 1.98e+01 2.739e+05 3.763e+03 2.593e+02\r", " 170 1.99e+01 2.750e+05 3.628e+03 2.611e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 171 2.00e+01 2.753e+05 3.627e+03 2.607e+02\r", " 172 2.02e+01 2.749e+05 3.751e+03 2.592e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 173 2.03e+01 2.739e+05 3.976e+03 2.571e+02\r", " 174 2.04e+01 2.724e+05 4.267e+03 2.535e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 175 2.05e+01 2.706e+05 4.591e+03 2.476e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 176 2.06e+01 2.685e+05 4.921e+03 2.413e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 177 2.09e+01 2.664e+05 5.235e+03 2.347e+02\r", " 178 2.11e+01 2.644e+05 5.516e+03 2.264e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 179 2.13e+01 2.626e+05 5.753e+03 2.178e+02\r", " 180 2.14e+01 2.611e+05 5.937e+03 2.088e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 181 2.16e+01 2.602e+05 6.060e+03 2.007e+02\r", " 182 2.18e+01 2.597e+05 6.120e+03 1.935e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 183 2.19e+01 2.599e+05 6.114e+03 1.880e+02\r", " 184 2.20e+01 2.607e+05 6.041e+03 1.846e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 185 2.22e+01 2.620e+05 5.903e+03 1.826e+02\r", " 186 2.23e+01 2.637e+05 5.703e+03 1.825e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 187 2.25e+01 2.659e+05 5.445e+03 1.846e+02\r", " 188 2.25e+01 2.683e+05 5.134e+03 1.878e+02\r", " 189 2.26e+01 2.709e+05 4.776e+03 1.923e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 190 2.27e+01 2.734e+05 4.382e+03 1.975e+02\r", " 191 2.27e+01 2.757e+05 3.960e+03 2.036e+02\r", " 192 2.28e+01 2.778e+05 3.525e+03 2.104e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 193 2.29e+01 2.795e+05 3.094e+03 2.169e+02\r", " 194 2.30e+01 2.807e+05 2.691e+03 2.215e+02\r", " 195 2.30e+01 2.813e+05 2.347e+03 2.269e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 196 2.31e+01 2.813e+05 2.105e+03 2.318e+02\r", " 197 2.32e+01 2.807e+05 2.005e+03 2.360e+02\r", " 198 2.32e+01 2.794e+05 2.063e+03 2.391e+02\r", " 199 2.33e+01 2.776e+05 2.254e+03 2.398e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 200 2.33e+01 2.753e+05 2.533e+03 2.388e+02\n", " 201 2.34e+01 2.725e+05 2.856e+03 2.357e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 202 2.36e+01 2.694e+05 3.190e+03 2.330e+02\r", " 203 2.37e+01 2.661e+05 3.513e+03 2.306e+02\r", " 204 2.38e+01 2.626e+05 3.812e+03 2.263e+02\r", " 205 2.39e+01 2.592e+05 4.074e+03 2.208e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 206 2.39e+01 2.558e+05 4.296e+03 2.120e+02\r", " 207 2.40e+01 2.526e+05 4.471e+03 2.043e+02\r", " 208 2.40e+01 2.497e+05 4.597e+03 1.977e+02\r", " 209 2.41e+01 2.470e+05 4.672e+03 1.915e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 210 2.41e+01 2.447e+05 4.698e+03 1.872e+02\r", " 211 2.42e+01 2.428e+05 4.675e+03 1.854e+02\r", " 212 2.42e+01 2.412e+05 4.607e+03 1.812e+02\r", " 213 2.43e+01 2.400e+05 4.496e+03 1.784e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 214 2.43e+01 2.390e+05 4.350e+03 1.770e+02\r", " 215 2.44e+01 2.382e+05 4.173e+03 1.751e+02\r", " 216 2.44e+01 2.376e+05 3.975e+03 1.726e+02\r", " 217 2.45e+01 2.371e+05 3.764e+03 1.723e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 218 2.45e+01 2.366e+05 3.551e+03 1.721e+02\r", " 219 2.46e+01 2.361e+05 3.349e+03 1.723e+02\r", " 220 2.47e+01 2.355e+05 3.171e+03 1.734e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 221 2.48e+01 2.348e+05 3.030e+03 1.741e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 222 2.51e+01 2.341e+05 2.937e+03 1.748e+02\r", " 223 2.53e+01 2.331e+05 2.899e+03 1.752e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 224 2.54e+01 2.321e+05 2.918e+03 1.762e+02\r", " 225 2.56e+01 2.310e+05 2.988e+03 1.792e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 226 2.57e+01 2.299e+05 3.098e+03 1.825e+02\r", " 227 2.58e+01 2.287e+05 3.235e+03 1.826e+02\r", " 228 2.59e+01 2.275e+05 3.387e+03 1.799e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 229 2.60e+01 2.263e+05 3.541e+03 1.768e+02\r", " 230 2.62e+01 2.251e+05 3.687e+03 1.728e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 231 2.64e+01 2.241e+05 3.816e+03 1.688e+02\r", " 232 2.66e+01 2.232e+05 3.921e+03 1.617e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 233 2.67e+01 2.225e+05 3.998e+03 1.542e+02\r", " 234 2.68e+01 2.219e+05 4.043e+03 1.488e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 235 2.69e+01 2.216e+05 4.053e+03 1.446e+02\r", " 236 2.71e+01 2.215e+05 4.026e+03 1.412e+02\r", " 237 2.72e+01 2.217e+05 3.964e+03 1.386e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 238 2.72e+01 2.220e+05 3.865e+03 1.373e+02\r", " 239 2.73e+01 2.224e+05 3.733e+03 1.369e+02\r", " 240 2.74e+01 2.230e+05 3.569e+03 1.378e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 241 2.75e+01 2.237e+05 3.377e+03 1.392e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 242 2.78e+01 2.244e+05 3.161e+03 1.414e+02\r", " 243 2.79e+01 2.251e+05 2.927e+03 1.440e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 244 2.80e+01 2.256e+05 2.683e+03 1.470e+02\r", " 245 2.82e+01 2.260e+05 2.435e+03 1.500e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 246 2.83e+01 2.263e+05 2.194e+03 1.530e+02\r", " 247 2.84e+01 2.263e+05 1.974e+03 1.560e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 248 2.85e+01 2.261e+05 1.790e+03 1.585e+02\r", " 249 2.86e+01 2.257e+05 1.659e+03 1.605e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 250 2.87e+01 2.249e+05 1.595e+03 1.623e+02\n", " 251 2.88e+01 2.240e+05 1.604e+03 1.635e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 252 2.90e+01 2.228e+05 1.679e+03 1.641e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 253 2.93e+01 2.215e+05 1.801e+03 1.643e+02\r", " 254 2.94e+01 2.200e+05 1.954e+03 1.643e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 255 2.96e+01 2.184e+05 2.120e+03 1.639e+02\r", " 256 2.97e+01 2.167e+05 2.286e+03 1.631e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 257 2.99e+01 2.151e+05 2.443e+03 1.617e+02\r", " 258 3.00e+01 2.135e+05 2.583e+03 1.591e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 259 3.02e+01 2.119e+05 2.702e+03 1.555e+02\r", " 260 3.03e+01 2.105e+05 2.797e+03 1.526e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 261 3.06e+01 2.092e+05 2.867e+03 1.499e+02\r", " 262 3.07e+01 2.081e+05 2.909e+03 1.468e+02\r", " 263 3.07e+01 2.071e+05 2.924e+03 1.438e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 264 3.08e+01 2.063e+05 2.914e+03 1.412e+02\r", " 265 3.09e+01 2.057e+05 2.880e+03 1.398e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 266 3.11e+01 2.052e+05 2.824e+03 1.395e+02\r", " 267 3.12e+01 2.049e+05 2.751e+03 1.395e+02\r", " 268 3.13e+01 2.047e+05 2.663e+03 1.382e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 269 3.14e+01 2.046e+05 2.567e+03 1.357e+02\r", " 270 3.15e+01 2.045e+05 2.468e+03 1.342e+02\r", " 271 3.15e+01 2.044e+05 2.372e+03 1.341e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 272 3.16e+01 2.043e+05 2.286e+03 1.340e+02\r", " 273 3.17e+01 2.042e+05 2.215e+03 1.334e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 274 3.20e+01 2.041e+05 2.166e+03 1.324e+02\r", " 275 3.22e+01 2.040e+05 2.143e+03 1.313e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 276 3.23e+01 2.038e+05 2.146e+03 1.300e+02\r", " 277 3.25e+01 2.035e+05 2.174e+03 1.279e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 278 3.27e+01 2.032e+05 2.226e+03 1.259e+02\r", " 279 3.29e+01 2.028e+05 2.294e+03 1.237e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 280 3.30e+01 2.024e+05 2.373e+03 1.215e+02\r", " 281 3.31e+01 2.021e+05 2.458e+03 1.194e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 282 3.34e+01 2.017e+05 2.542e+03 1.178e+02\r", " 283 3.35e+01 2.014e+05 2.619e+03 1.156e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 284 3.37e+01 2.012e+05 2.687e+03 1.128e+02\r", " 285 3.38e+01 2.010e+05 2.741e+03 1.104e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 286 3.40e+01 2.008e+05 2.779e+03 1.083e+02\r", " 287 3.41e+01 2.008e+05 2.799e+03 1.065e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 288 3.42e+01 2.009e+05 2.799e+03 1.054e+02\r", " 289 3.43e+01 2.010e+05 2.779e+03 1.045e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 290 3.45e+01 2.013e+05 2.739e+03 1.037e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 291 3.47e+01 2.016e+05 2.680e+03 1.032e+02\r", " 292 3.49e+01 2.020e+05 2.602e+03 1.034e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 293 3.50e+01 2.024e+05 2.508e+03 1.031e+02\r", " 294 3.51e+01 2.029e+05 2.400e+03 1.033e+02\r", " 295 3.52e+01 2.034e+05 2.280e+03 1.041e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 296 3.53e+01 2.039e+05 2.152e+03 1.054e+02\r", " 297 3.53e+01 2.043e+05 2.020e+03 1.072e+02\r", " 298 3.54e+01 2.047e+05 1.889e+03 1.088e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 299 3.55e+01 2.050e+05 1.763e+03 1.108e+02\r", " 300 3.55e+01 2.052e+05 1.648e+03 1.128e+02\n", " 301 3.56e+01 2.054e+05 1.550e+03 1.144e+02\r", " 302 3.56e+01 2.054e+05 1.474e+03 1.159e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 303 3.57e+01 2.054e+05 1.424e+03 1.172e+02\r", " 304 3.57e+01 2.052e+05 1.402e+03 1.183e+02\r", " 305 3.58e+01 2.049e+05 1.406e+03 1.192e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 306 3.59e+01 2.046e+05 1.433e+03 1.200e+02\r", " 307 3.61e+01 2.042e+05 1.477e+03 1.208e+02\r", " 308 3.62e+01 2.037e+05 1.532e+03 1.212e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 309 3.63e+01 2.031e+05 1.592e+03 1.212e+02\r", " 310 3.63e+01 2.026e+05 1.651e+03 1.211e+02\r", " 311 3.64e+01 2.020e+05 1.705e+03 1.210e+02\r", " 312 3.64e+01 2.014e+05 1.751e+03 1.214e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 313 3.65e+01 2.008e+05 1.787e+03 1.213e+02\r", " 314 3.65e+01 2.002e+05 1.811e+03 1.205e+02\r", " 315 3.66e+01 1.996e+05 1.823e+03 1.194e+02\r", " 316 3.66e+01 1.991e+05 1.822e+03 1.182e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 317 3.67e+01 1.986e+05 1.809e+03 1.171e+02\r", " 318 3.67e+01 1.981e+05 1.785e+03 1.160e+02\r", " 319 3.68e+01 1.976e+05 1.752e+03 1.148e+02\r", " 320 3.68e+01 1.972e+05 1.711e+03 1.141e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 321 3.69e+01 1.968e+05 1.667e+03 1.137e+02\r", " 322 3.69e+01 1.964e+05 1.620e+03 1.137e+02\r", " 323 3.70e+01 1.960e+05 1.576e+03 1.138e+02\r", " 324 3.70e+01 1.956e+05 1.537e+03 1.133e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 325 3.71e+01 1.952e+05 1.507e+03 1.131e+02\r", " 326 3.72e+01 1.948e+05 1.488e+03 1.131e+02\r", " 327 3.72e+01 1.944e+05 1.483e+03 1.134e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 328 3.75e+01 1.940e+05 1.493e+03 1.131e+02\r", " 329 3.77e+01 1.935e+05 1.517e+03 1.133e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 330 3.78e+01 1.931e+05 1.553e+03 1.132e+02\r", " 331 3.78e+01 1.927e+05 1.599e+03 1.126e+02\r", " 332 3.79e+01 1.923e+05 1.653e+03 1.122e+02\r", " 333 3.79e+01 1.919e+05 1.711e+03 1.118e+02\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 334 3.80e+01 1.915e+05 1.769e+03 1.108e+02\r", " 335 3.80e+01 1.912e+05 1.825e+03 1.075e+02\r", " 336 3.81e+01 1.908e+05 1.877e+03 1.036e+02\r", " 337 3.81e+01 1.904e+05 1.923e+03 9.975e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 338 3.82e+01 1.901e+05 1.960e+03 9.642e+01\r", " 339 3.82e+01 1.898e+05 1.988e+03 9.338e+01\r", " 340 3.83e+01 1.896e+05 2.005e+03 9.051e+01\r", " 341 3.83e+01 1.894e+05 2.011e+03 8.802e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 342 3.84e+01 1.892e+05 2.006e+03 8.590e+01\r", " 343 3.84e+01 1.891e+05 1.990e+03 8.400e+01\r", " 344 3.85e+01 1.890e+05 1.962e+03 8.252e+01\r", " 345 3.86e+01 1.890e+05 1.925e+03 8.113e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 346 3.86e+01 1.890e+05 1.879e+03 7.985e+01\r", " 347 3.87e+01 1.890e+05 1.825e+03 7.895e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 348 3.89e+01 1.891e+05 1.764e+03 7.867e+01\r", " 349 3.90e+01 1.891e+05 1.700e+03 7.903e+01\r", " 350 3.90e+01 1.892e+05 1.632e+03 7.943e+01\n", " 351 3.91e+01 1.892e+05 1.564e+03 7.929e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 352 3.91e+01 1.892e+05 1.498e+03 7.924e+01\r", " 353 3.92e+01 1.893e+05 1.435e+03 7.937e+01\r", " 354 3.92e+01 1.892e+05 1.378e+03 7.948e+01\r", " 355 3.93e+01 1.892e+05 1.328e+03 7.996e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 356 3.93e+01 1.891e+05 1.286e+03 8.064e+01\r", " 357 3.94e+01 1.891e+05 1.253e+03 8.109e+01\r", " 358 3.94e+01 1.890e+05 1.228e+03 8.142e+01\r", " 359 3.95e+01 1.889e+05 1.212e+03 8.232e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 360 3.95e+01 1.887e+05 1.201e+03 8.361e+01\r", " 361 3.96e+01 1.886e+05 1.195e+03 8.481e+01\r", " 362 3.96e+01 1.884e+05 1.192e+03 8.597e+01\r", " 363 3.97e+01 1.883e+05 1.190e+03 8.708e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 364 3.98e+01 1.881e+05 1.187e+03 8.780e+01\r", " 365 3.98e+01 1.880e+05 1.181e+03 8.863e+01\r", " 366 3.99e+01 1.878e+05 1.172e+03 8.957e+01\r", " 367 3.99e+01 1.877e+05 1.158e+03 9.051e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 368 4.00e+01 1.876e+05 1.139e+03 9.062e+01\r", " 369 4.00e+01 1.874e+05 1.115e+03 8.953e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 370 4.01e+01 1.873e+05 1.087e+03 8.860e+01\r", " 371 4.03e+01 1.871e+05 1.055e+03 8.817e+01\r", " 372 4.04e+01 1.870e+05 1.020e+03 8.824e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 373 4.05e+01 1.868e+05 9.837e+02 8.844e+01\r", " 374 4.05e+01 1.867e+05 9.480e+02 8.903e+01\r", " 375 4.06e+01 1.865e+05 9.150e+02 8.968e+01\r", " 376 4.06e+01 1.863e+05 8.873e+02 8.978e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 377 4.07e+01 1.861e+05 8.675e+02 8.877e+01\r", " 378 4.07e+01 1.859e+05 8.579e+02 8.762e+01\r", " 379 4.08e+01 1.857e+05 8.598e+02 8.687e+01\r", " 380 4.08e+01 1.854e+05 8.737e+02 8.627e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 381 4.09e+01 1.852e+05 8.991e+02 8.598e+01\r", " 382 4.09e+01 1.849e+05 9.345e+02 8.575e+01\r", " 383 4.10e+01 1.846e+05 9.779e+02 8.545e+01\r", " 384 4.10e+01 1.843e+05 1.027e+03 8.549e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 385 4.11e+01 1.840e+05 1.079e+03 8.587e+01\r", " 386 4.11e+01 1.838e+05 1.133e+03 8.667e+01\r", " 387 4.12e+01 1.835e+05 1.185e+03 8.770e+01\r", " 388 4.12e+01 1.833e+05 1.235e+03 8.829e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 389 4.13e+01 1.831e+05 1.280e+03 8.849e+01\r", " 390 4.13e+01 1.829e+05 1.321e+03 8.996e+01\r", " 391 4.14e+01 1.828e+05 1.356e+03 8.590e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 392 4.15e+01 1.826e+05 1.385e+03 8.433e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 393 4.17e+01 1.825e+05 1.407e+03 8.353e+01\r", " 394 4.18e+01 1.823e+05 1.423e+03 8.093e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 395 4.20e+01 1.822e+05 1.431e+03 7.712e+01\r", " 396 4.20e+01 1.820e+05 1.434e+03 7.374e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 397 4.22e+01 1.819e+05 1.430e+03 7.160e+01\r", " 398 4.23e+01 1.817e+05 1.420e+03 6.969e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 399 4.24e+01 1.816e+05 1.405e+03 6.794e+01\r", " 400 4.25e+01 1.815e+05 1.385e+03 6.610e+01\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 401 4.26e+01 1.814e+05 1.362e+03 6.398e+01\r", " 402 4.27e+01 1.813e+05 1.335e+03 6.240e+01\r", " 403 4.28e+01 1.812e+05 1.307e+03 6.161e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 404 4.31e+01 1.811e+05 1.277e+03 6.100e+01\r", " 405 4.32e+01 1.810e+05 1.247e+03 6.074e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 406 4.33e+01 1.810e+05 1.217e+03 6.058e+01\r", " 407 4.35e+01 1.809e+05 1.188e+03 6.034e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 408 4.36e+01 1.808e+05 1.160e+03 6.034e+01\r", " 409 4.37e+01 1.807e+05 1.134e+03 6.028e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 410 4.38e+01 1.806e+05 1.110e+03 6.009e+01\r", " 411 4.39e+01 1.805e+05 1.088e+03 6.010e+01\r", " 412 4.40e+01 1.804e+05 1.067e+03 6.027e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 413 4.41e+01 1.803e+05 1.047e+03 6.057e+01\r", " 414 4.42e+01 1.802e+05 1.028e+03 6.067e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 415 4.45e+01 1.801e+05 1.009e+03 6.061e+01\r", " 416 4.46e+01 1.800e+05 9.890e+02 6.043e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 417 4.48e+01 1.799e+05 9.682e+02 6.024e+01\r", " 418 4.50e+01 1.797e+05 9.458e+02 6.017e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 419 4.52e+01 1.796e+05 9.213e+02 6.011e+01\r", " 420 4.53e+01 1.794e+05 8.946e+02 6.015e+01\r", " 421 4.54e+01 1.792e+05 8.654e+02 6.011e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 422 4.54e+01 1.790e+05 8.339e+02 6.021e+01\r", " 423 4.55e+01 1.788e+05 8.004e+02 6.033e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 424 4.58e+01 1.786e+05 7.652e+02 6.053e+01\r", " 425 4.60e+01 1.784e+05 7.290e+02 6.063e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 426 4.61e+01 1.782e+05 6.928e+02 6.077e+01\r", " 427 4.62e+01 1.779e+05 6.578e+02 6.091e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 428 4.63e+01 1.777e+05 6.254e+02 6.110e+01\r", " 429 4.65e+01 1.774e+05 5.971e+02 6.138e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 430 4.67e+01 1.771e+05 5.747e+02 6.167e+01\r", " 431 4.68e+01 1.768e+05 5.598e+02 6.196e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 432 4.69e+01 1.765e+05 5.534e+02 6.209e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 433 4.72e+01 1.762e+05 5.562e+02 6.210e+01\r", " 434 4.74e+01 1.759e+05 5.680e+02 6.217e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 435 4.75e+01 1.755e+05 5.877e+02 6.220e+01\r", " 436 4.76e+01 1.752e+05 6.141e+02 6.212e+01\r", " 437 4.77e+01 1.749e+05 6.455e+02 6.207e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 438 4.78e+01 1.746e+05 6.801e+02 6.190e+01\r", " 439 4.80e+01 1.743e+05 7.165e+02 6.178e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 440 4.81e+01 1.739e+05 7.532e+02 6.183e+01\r", " 441 4.82e+01 1.736e+05 7.890e+02 6.185e+01\r", " 442 4.83e+01 1.734e+05 8.231e+02 6.177e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 443 4.84e+01 1.731e+05 8.547e+02 6.164e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 444 4.86e+01 1.729e+05 8.832e+02 6.152e+01\r", " 445 4.88e+01 1.726e+05 9.082e+02 6.158e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 446 4.89e+01 1.724e+05 9.295e+02 6.165e+01\r", " 447 4.90e+01 1.723e+05 9.469e+02 6.188e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 448 4.91e+01 1.721e+05 9.603e+02 6.212e+01\r", " 449 4.92e+01 1.720e+05 9.698e+02 6.257e+01\r", " 450 4.93e+01 1.719e+05 9.755e+02 6.438e+01\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 451 4.95e+01 1.719e+05 9.773e+02 6.674e+01\r", " 452 4.96e+01 1.719e+05 9.761e+02 6.672e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 453 4.97e+01 1.718e+05 9.724e+02 6.352e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 454 4.99e+01 1.718e+05 9.665e+02 6.153e+01\r", " 455 5.01e+01 1.718e+05 9.588e+02 6.024e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 456 5.02e+01 1.718e+05 9.497e+02 5.907e+01\r", " 457 5.03e+01 1.718e+05 9.396e+02 5.824e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 458 5.04e+01 1.718e+05 9.289e+02 5.731e+01\r", " 459 5.06e+01 1.719e+05 9.177e+02 5.659e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 460 5.07e+01 1.719e+05 9.065e+02 5.609e+01\r", " 461 5.08e+01 1.720e+05 8.951e+02 5.596e+01\r", " 462 5.09e+01 1.720e+05 8.838e+02 5.649e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 463 5.10e+01 1.721e+05 8.725e+02 5.669e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 464 5.13e+01 1.722e+05 8.611e+02 5.702e+01\r", " 465 5.14e+01 1.723e+05 8.495e+02 5.715e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 466 5.16e+01 1.724e+05 8.375e+02 5.685e+01\r", " 467 5.17e+01 1.725e+05 8.251e+02 5.642e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 468 5.19e+01 1.727e+05 8.118e+02 5.662e+01\r", " 469 5.20e+01 1.728e+05 7.974e+02 5.683e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 470 5.21e+01 1.729e+05 7.820e+02 5.542e+01\r", " 471 5.22e+01 1.730e+05 7.657e+02 5.114e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 472 5.23e+01 1.731e+05 7.482e+02 4.898e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 473 5.25e+01 1.732e+05 7.293e+02 4.771e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 474 5.28e+01 1.733e+05 7.091e+02 4.660e+01\r", " 475 5.29e+01 1.734e+05 6.877e+02 4.590e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 476 5.31e+01 1.735e+05 6.653e+02 4.519e+01\r", " 477 5.33e+01 1.735e+05 6.421e+02 4.462e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 478 5.35e+01 1.736e+05 6.186e+02 4.417e+01\r", " 479 5.36e+01 1.737e+05 5.952e+02 4.383e+01\r", " 480 5.37e+01 1.737e+05 5.725e+02 4.371e+01\r", " 481 5.37e+01 1.737e+05 5.511e+02 4.371e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 482 5.38e+01 1.738e+05 5.317e+02 4.379e+01\r", " 483 5.38e+01 1.738e+05 5.150e+02 4.384e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 484 5.40e+01 1.738e+05 5.016e+02 4.402e+01\r", " 485 5.41e+01 1.738e+05 4.920e+02 4.413e+01\r", " 486 5.42e+01 1.738e+05 4.867e+02 4.438e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 487 5.43e+01 1.737e+05 4.856e+02 4.458e+01\r", " 488 5.43e+01 1.737e+05 4.887e+02 4.495e+01\r", " 489 5.44e+01 1.737e+05 4.957e+02 4.507e+01\r", " 490 5.44e+01 1.736e+05 5.058e+02 4.528e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 491 5.45e+01 1.736e+05 5.186e+02 4.550e+01\r", " 492 5.45e+01 1.735e+05 5.332e+02 4.563e+01\r", " 493 5.46e+01 1.734e+05 5.490e+02 4.567e+01\r", " 494 5.46e+01 1.733e+05 5.654e+02 4.555e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 495 5.47e+01 1.733e+05 5.817e+02 4.542e+01\r", " 496 5.47e+01 1.732e+05 5.974e+02 4.528e+01\r", " 497 5.48e+01 1.731e+05 6.122e+02 4.521e+01\r", " 498 5.48e+01 1.730e+05 6.256e+02 4.525e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 499 5.49e+01 1.729e+05 6.375e+02 4.545e+01\r", " 500 5.49e+01 1.728e+05 6.477e+02 4.537e+01\n", " 501 5.50e+01 1.727e+05 6.562e+02 4.512e+01\r", " 502 5.50e+01 1.726e+05 6.629e+02 4.479e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 503 5.51e+01 1.725e+05 6.680e+02 4.443e+01\r", " 504 5.52e+01 1.724e+05 6.715e+02 4.406e+01\r", " 505 5.52e+01 1.723e+05 6.735e+02 4.365e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 506 5.54e+01 1.723e+05 6.743e+02 4.329e+01\r", " 507 5.55e+01 1.722e+05 6.740e+02 4.293e+01\r", " 508 5.56e+01 1.721e+05 6.728e+02 4.265e+01\r", " 509 5.56e+01 1.720e+05 6.708e+02 4.245e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 510 5.57e+01 1.719e+05 6.683e+02 4.224e+01\r", " 511 5.57e+01 1.719e+05 6.652e+02 4.211e+01\r", " 512 5.58e+01 1.718e+05 6.617e+02 4.207e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 513 5.59e+01 1.717e+05 6.580e+02 4.205e+01\r", " 514 5.60e+01 1.717e+05 6.540e+02 4.203e+01\r", " 515 5.60e+01 1.716e+05 6.497e+02 4.190e+01\r", " 516 5.61e+01 1.715e+05 6.451e+02 4.160e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 517 5.61e+01 1.715e+05 6.402e+02 4.128e+01\r", " 518 5.62e+01 1.714e+05 6.349e+02 4.086e+01\r", " 519 5.62e+01 1.714e+05 6.291e+02 4.040e+01\r", " 520 5.63e+01 1.713e+05 6.227e+02 3.996e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 521 5.63e+01 1.713e+05 6.155e+02 3.966e+01\r", " 522 5.64e+01 1.712e+05 6.075e+02 3.942e+01\r", " 523 5.64e+01 1.712e+05 5.985e+02 3.914e+01\r", " 524 5.65e+01 1.711e+05 5.885e+02 3.891e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 525 5.66e+01 1.711e+05 5.774e+02 3.872e+01\r", " 526 5.67e+01 1.710e+05 5.653e+02 3.843e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 527 5.69e+01 1.710e+05 5.522e+02 3.800e+01\r", " 528 5.69e+01 1.709e+05 5.381e+02 3.767e+01\r", " 529 5.70e+01 1.709e+05 5.233e+02 3.726e+01\r", " 530 5.71e+01 1.708e+05 5.078e+02 3.708e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 531 5.71e+01 1.708e+05 4.919e+02 3.695e+01\r", " 532 5.72e+01 1.707e+05 4.759e+02 3.687e+01\r", " 533 5.72e+01 1.707e+05 4.600e+02 3.685e+01\r", " 534 5.73e+01 1.706e+05 4.447e+02 3.683e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 535 5.73e+01 1.706e+05 4.302e+02 3.696e+01\r", " 536 5.74e+01 1.705e+05 4.169e+02 3.725e+01\r", " 537 5.74e+01 1.704e+05 4.052e+02 3.763e+01\r", " 538 5.75e+01 1.704e+05 3.954e+02 3.797e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 539 5.76e+01 1.703e+05 3.879e+02 3.763e+01\r", " 540 5.76e+01 1.702e+05 3.829e+02 3.731e+01\r", " 541 5.77e+01 1.702e+05 3.804e+02 3.714e+01\r", " 542 5.77e+01 1.701e+05 3.803e+02 3.691e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 543 5.78e+01 1.700e+05 3.826e+02 3.676e+01\r", " 544 5.78e+01 1.699e+05 3.869e+02 3.662e+01\r", " 545 5.79e+01 1.698e+05 3.930e+02 3.658e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 546 5.80e+01 1.697e+05 4.004e+02 3.641e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 547 5.83e+01 1.696e+05 4.089e+02 3.612e+01\r", " 548 5.85e+01 1.695e+05 4.181e+02 3.563e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 549 5.86e+01 1.694e+05 4.277e+02 3.531e+01\r", " 550 5.88e+01 1.693e+05 4.373e+02 3.503e+01\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 551 5.90e+01 1.692e+05 4.468e+02 3.476e+01\r", " 552 5.91e+01 1.691e+05 4.560e+02 3.446e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 553 5.92e+01 1.690e+05 4.646e+02 3.418e+01\r", " 554 5.93e+01 1.689e+05 4.727e+02 3.392e+01\r", " 555 5.94e+01 1.688e+05 4.801e+02 3.363e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 556 5.96e+01 1.687e+05 4.868e+02 3.339e+01\r", " 557 5.97e+01 1.686e+05 4.928e+02 3.313e+01\r", " 558 5.97e+01 1.685e+05 4.981e+02 3.289e+01\r", " 559 5.98e+01 1.684e+05 5.027e+02 3.274e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 560 5.99e+01 1.684e+05 5.067e+02 3.260e+01\r", " 561 5.99e+01 1.683e+05 5.100e+02 3.249e+01\r", " 562 6.00e+01 1.682e+05 5.127e+02 3.237e+01\r", " 563 6.00e+01 1.682e+05 5.149e+02 3.226e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 564 6.01e+01 1.681e+05 5.164e+02 3.215e+01\r", " 565 6.01e+01 1.680e+05 5.175e+02 3.192e+01\r", " 566 6.02e+01 1.680e+05 5.179e+02 3.182e+01\r", " 567 6.02e+01 1.680e+05 5.177e+02 3.162e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 568 6.03e+01 1.679e+05 5.169e+02 3.147e+01\r", " 569 6.03e+01 1.679e+05 5.155e+02 3.123e+01\r", " 570 6.04e+01 1.679e+05 5.133e+02 3.102e+01\r", " 571 6.04e+01 1.679e+05 5.104e+02 3.085e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 572 6.05e+01 1.679e+05 5.066e+02 3.066e+01\r", " 573 6.05e+01 1.679e+05 5.019e+02 3.046e+01\r", " 574 6.06e+01 1.679e+05 4.963e+02 3.020e+01\r", " 575 6.06e+01 1.679e+05 4.897e+02 2.997e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 576 6.07e+01 1.679e+05 4.821e+02 2.977e+01\r", " 577 6.08e+01 1.679e+05 4.735e+02 2.958e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 578 6.10e+01 1.679e+05 4.639e+02 2.945e+01\r", " 579 6.11e+01 1.679e+05 4.532e+02 2.930e+01\r", " 580 6.12e+01 1.679e+05 4.416e+02 2.922e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 581 6.13e+01 1.680e+05 4.291e+02 2.920e+01\r", " 582 6.13e+01 1.680e+05 4.158e+02 2.917e+01\r", " 583 6.14e+01 1.680e+05 4.019e+02 2.912e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 584 6.15e+01 1.680e+05 3.875e+02 2.905e+01\r", " 585 6.16e+01 1.681e+05 3.728e+02 2.891e+01\r", " 586 6.17e+01 1.681e+05 3.580e+02 2.886e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 587 6.17e+01 1.681e+05 3.433e+02 2.883e+01\r", " 588 6.18e+01 1.681e+05 3.290e+02 2.887e+01\r", " 589 6.19e+01 1.682e+05 3.153e+02 2.893e+01\r", " 590 6.19e+01 1.682e+05 3.026e+02 2.892e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 591 6.20e+01 1.682e+05 2.910e+02 2.892e+01\r", " 592 6.20e+01 1.682e+05 2.810e+02 2.890e+01\r", " 593 6.21e+01 1.682e+05 2.727e+02 2.892e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 594 6.23e+01 1.682e+05 2.663e+02 2.895e+01\r", " 595 6.25e+01 1.682e+05 2.620e+02 2.897e+01\r", " 596 6.25e+01 1.682e+05 2.598e+02 2.903e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 597 6.26e+01 1.683e+05 2.597e+02 2.906e+01\r", " 598 6.27e+01 1.683e+05 2.616e+02 2.909e+01\r", " 599 6.27e+01 1.683e+05 2.652e+02 2.912e+01\r", " 600 6.28e+01 1.682e+05 2.703e+02 2.917e+01\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 601 6.29e+01 1.682e+05 2.766e+02 2.927e+01\r", " 602 6.29e+01 1.682e+05 2.839e+02 2.935e+01\r", " 603 6.30e+01 1.682e+05 2.919e+02 2.946e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 604 6.31e+01 1.682e+05 3.004e+02 2.956e+01\r", " 605 6.31e+01 1.682e+05 3.092e+02 2.968e+01\r", " 606 6.32e+01 1.682e+05 3.180e+02 2.970e+01\r", " 607 6.32e+01 1.682e+05 3.269e+02 2.974e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 608 6.33e+01 1.682e+05 3.356e+02 2.967e+01\r", " 609 6.33e+01 1.682e+05 3.441e+02 2.956e+01\r", " 610 6.34e+01 1.682e+05 3.523e+02 2.932e+01\r", " 611 6.35e+01 1.682e+05 3.602e+02 2.911e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 612 6.35e+01 1.682e+05 3.677e+02 2.887e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 613 6.37e+01 1.682e+05 3.747e+02 2.872e+01\r", " 614 6.38e+01 1.682e+05 3.813e+02 2.856e+01\r", " 615 6.39e+01 1.682e+05 3.873e+02 2.838e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 616 6.40e+01 1.682e+05 3.928e+02 2.818e+01\r", " 617 6.40e+01 1.681e+05 3.977e+02 2.800e+01\r", " 618 6.41e+01 1.681e+05 4.019e+02 2.780e+01\r", " 619 6.41e+01 1.681e+05 4.055e+02 2.763e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 620 6.42e+01 1.681e+05 4.084e+02 2.748e+01\r", " 621 6.42e+01 1.681e+05 4.105e+02 2.738e+01\r", " 622 6.43e+01 1.681e+05 4.118e+02 2.721e+01\r", " 623 6.43e+01 1.681e+05 4.123e+02 2.708e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 624 6.44e+01 1.681e+05 4.120e+02 2.694e+01\r", " 625 6.44e+01 1.681e+05 4.108e+02 2.666e+01\r", " 626 6.45e+01 1.681e+05 4.087e+02 2.640e+01\r", " 627 6.45e+01 1.681e+05 4.057e+02 2.615e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 628 6.46e+01 1.681e+05 4.017e+02 2.594e+01\r", " 629 6.46e+01 1.681e+05 3.968e+02 2.571e+01\r", " 630 6.47e+01 1.681e+05 3.910e+02 2.547e+01\r", " 631 6.47e+01 1.681e+05 3.843e+02 2.530e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 632 6.48e+01 1.681e+05 3.767e+02 2.514e+01\r", " 633 6.49e+01 1.680e+05 3.682e+02 2.504e+01\r", " 634 6.49e+01 1.680e+05 3.590e+02 2.491e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 635 6.50e+01 1.680e+05 3.491e+02 2.484e+01\r", " 636 6.52e+01 1.680e+05 3.385e+02 2.476e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 637 6.53e+01 1.680e+05 3.273e+02 2.478e+01\r", " 638 6.53e+01 1.680e+05 3.157e+02 2.480e+01\r", " 639 6.54e+01 1.680e+05 3.037e+02 2.481e+01\r", " 640 6.54e+01 1.680e+05 2.915e+02 2.485e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 641 6.55e+01 1.680e+05 2.793e+02 2.485e+01\r", " 642 6.55e+01 1.679e+05 2.671e+02 2.492e+01\r", " 643 6.56e+01 1.679e+05 2.553e+02 2.496e+01\r", " 644 6.56e+01 1.679e+05 2.439e+02 2.501e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 645 6.57e+01 1.679e+05 2.332e+02 2.506e+01\r", " 646 6.58e+01 1.678e+05 2.234e+02 2.511e+01\r", " 647 6.58e+01 1.678e+05 2.145e+02 2.518e+01\r", " 648 6.59e+01 1.678e+05 2.070e+02 2.521e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 649 6.59e+01 1.677e+05 2.008e+02 2.528e+01\r", " 650 6.60e+01 1.677e+05 1.962e+02 2.534e+01\n", " 651 6.60e+01 1.676e+05 1.931e+02 2.539e+01\r", " 652 6.61e+01 1.676e+05 1.916e+02 2.545e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 653 6.61e+01 1.675e+05 1.918e+02 2.544e+01\r", " 654 6.62e+01 1.675e+05 1.934e+02 2.542e+01\r", " 655 6.62e+01 1.674e+05 1.963e+02 2.541e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 656 6.63e+01 1.673e+05 2.004e+02 2.539e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 657 6.65e+01 1.673e+05 2.054e+02 2.534e+01\r", " 658 6.66e+01 1.672e+05 2.113e+02 2.529e+01\r", " 659 6.67e+01 1.672e+05 2.177e+02 2.521e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 660 6.68e+01 1.671e+05 2.246e+02 2.515e+01\r", " 661 6.68e+01 1.670e+05 2.318e+02 2.505e+01\r", " 662 6.69e+01 1.670e+05 2.392e+02 2.498e+01\r", " 663 6.69e+01 1.669e+05 2.466e+02 2.489e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 664 6.70e+01 1.669e+05 2.539e+02 2.477e+01\r", " 665 6.71e+01 1.668e+05 2.612e+02 2.467e+01\r", " 666 6.71e+01 1.667e+05 2.682e+02 2.457e+01\r", " 667 6.72e+01 1.667e+05 2.749e+02 2.445e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 668 6.72e+01 1.666e+05 2.813e+02 2.431e+01\r", " 669 6.73e+01 1.666e+05 2.874e+02 2.415e+01\r", " 670 6.74e+01 1.665e+05 2.929e+02 2.396e+01\r", " 671 6.74e+01 1.665e+05 2.980e+02 2.383e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 672 6.75e+01 1.665e+05 3.026e+02 2.363e+01\r", " 673 6.75e+01 1.664e+05 3.066e+02 2.341e+01\r", " 674 6.76e+01 1.664e+05 3.100e+02 2.319e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 675 6.77e+01 1.664e+05 3.128e+02 2.296e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 676 6.79e+01 1.663e+05 3.149e+02 2.272e+01\r", " 677 6.80e+01 1.663e+05 3.163e+02 2.250e+01\r", " 678 6.81e+01 1.663e+05 3.170e+02 2.228e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 679 6.82e+01 1.663e+05 3.170e+02 2.205e+01\r", " 680 6.82e+01 1.663e+05 3.163e+02 2.180e+01\r", " 681 6.83e+01 1.663e+05 3.148e+02 2.156e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 682 6.84e+01 1.662e+05 3.127e+02 2.130e+01\r", " 683 6.85e+01 1.662e+05 3.098e+02 2.109e+01\r", " 684 6.86e+01 1.662e+05 3.062e+02 2.086e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 685 6.86e+01 1.662e+05 3.019e+02 2.064e+01\r", " 686 6.87e+01 1.662e+05 2.970e+02 2.044e+01\r", " 687 6.88e+01 1.662e+05 2.914e+02 2.027e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 688 6.89e+01 1.662e+05 2.853e+02 2.012e+01\r", " 689 6.89e+01 1.663e+05 2.787e+02 1.996e+01\r", " 690 6.90e+01 1.663e+05 2.716e+02 1.982e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 691 6.91e+01 1.663e+05 2.641e+02 1.971e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 692 6.93e+01 1.663e+05 2.563e+02 1.960e+01\r", " 693 6.94e+01 1.663e+05 2.482e+02 1.950e+01\r", " 694 6.95e+01 1.663e+05 2.399e+02 1.942e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 695 6.96e+01 1.663e+05 2.315e+02 1.934e+01\r", " 696 6.97e+01 1.663e+05 2.231e+02 1.927e+01\r", " 697 6.97e+01 1.663e+05 2.148e+02 1.921e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 698 6.98e+01 1.664e+05 2.067e+02 1.916e+01\r", " 699 6.99e+01 1.664e+05 1.989e+02 1.913e+01\r", " 700 6.99e+01 1.664e+05 1.915e+02 1.911e+01\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 701 7.00e+01 1.664e+05 1.846e+02 1.911e+01\r", " 702 7.01e+01 1.664e+05 1.783e+02 1.911e+01\r", " 703 7.02e+01 1.664e+05 1.727e+02 1.911e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 704 7.03e+01 1.664e+05 1.679e+02 1.909e+01\r", " 705 7.03e+01 1.664e+05 1.640e+02 1.907e+01\r", " 706 7.04e+01 1.664e+05 1.610e+02 1.905e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 707 7.05e+01 1.665e+05 1.591e+02 1.905e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 708 7.07e+01 1.665e+05 1.581e+02 1.905e+01\r", " 709 7.08e+01 1.665e+05 1.581e+02 1.908e+01\r", " 710 7.09e+01 1.665e+05 1.590e+02 1.908e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 711 7.09e+01 1.665e+05 1.607e+02 1.909e+01\r", " 712 7.10e+01 1.665e+05 1.633e+02 1.912e+01\r", " 713 7.11e+01 1.665e+05 1.665e+02 1.911e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 714 7.12e+01 1.665e+05 1.703e+02 1.910e+01\r", " 715 7.13e+01 1.665e+05 1.746e+02 1.909e+01\r", " 716 7.13e+01 1.665e+05 1.792e+02 1.907e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 717 7.14e+01 1.665e+05 1.841e+02 1.906e+01\r", " 718 7.15e+01 1.665e+05 1.892e+02 1.904e+01\r", " 719 7.16e+01 1.665e+05 1.945e+02 1.902e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 720 7.17e+01 1.665e+05 1.997e+02 1.896e+01\r", " 721 7.17e+01 1.665e+05 2.049e+02 1.892e+01\r", " 722 7.18e+01 1.665e+05 2.100e+02 1.887e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 723 7.19e+01 1.665e+05 2.149e+02 1.882e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 724 7.21e+01 1.665e+05 2.196e+02 1.879e+01\r", " 725 7.22e+01 1.665e+05 2.241e+02 1.874e+01\r", " 726 7.23e+01 1.665e+05 2.282e+02 1.869e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 727 7.24e+01 1.665e+05 2.319e+02 1.866e+01\r", " 728 7.24e+01 1.665e+05 2.353e+02 1.861e+01\r", " 729 7.25e+01 1.665e+05 2.383e+02 1.855e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 730 7.26e+01 1.665e+05 2.408e+02 1.851e+01\r", " 731 7.27e+01 1.665e+05 2.428e+02 1.847e+01\r", " 732 7.27e+01 1.665e+05 2.443e+02 1.843e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 733 7.28e+01 1.665e+05 2.453e+02 1.838e+01\r", " 734 7.29e+01 1.665e+05 2.458e+02 1.831e+01\r", " 735 7.29e+01 1.665e+05 2.457e+02 1.828e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 736 7.30e+01 1.665e+05 2.451e+02 1.823e+01\r", " 737 7.31e+01 1.665e+05 2.440e+02 1.821e+01\r", " 738 7.31e+01 1.665e+05 2.423e+02 1.815e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 739 7.32e+01 1.665e+05 2.402e+02 1.812e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 740 7.34e+01 1.665e+05 2.375e+02 1.807e+01\r", " 741 7.36e+01 1.665e+05 2.344e+02 1.803e+01\r", " 742 7.36e+01 1.665e+05 2.308e+02 1.800e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 743 7.37e+01 1.665e+05 2.268e+02 1.799e+01\r", " 744 7.38e+01 1.665e+05 2.224e+02 1.800e+01\r", " 745 7.39e+01 1.665e+05 2.176e+02 1.796e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 746 7.39e+01 1.665e+05 2.125e+02 1.795e+01\r", " 747 7.40e+01 1.665e+05 2.071e+02 1.791e+01\r", " 748 7.41e+01 1.665e+05 2.015e+02 1.786e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 749 7.42e+01 1.665e+05 1.956e+02 1.784e+01\r", " 750 7.42e+01 1.665e+05 1.896e+02 1.782e+01\n", " 751 7.43e+01 1.665e+05 1.835e+02 1.781e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 752 7.44e+01 1.665e+05 1.773e+02 1.781e+01\r", " 753 7.44e+01 1.665e+05 1.711e+02 1.779e+01\r", " 754 7.45e+01 1.665e+05 1.650e+02 1.777e+01\r", " 755 7.46e+01 1.664e+05 1.591e+02 1.776e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 756 7.47e+01 1.664e+05 1.533e+02 1.774e+01\r", " 757 7.49e+01 1.664e+05 1.479e+02 1.773e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 758 7.50e+01 1.664e+05 1.427e+02 1.770e+01\r", " 759 7.51e+01 1.664e+05 1.380e+02 1.766e+01\r", " 760 7.51e+01 1.664e+05 1.338e+02 1.760e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 761 7.52e+01 1.663e+05 1.302e+02 1.756e+01\r", " 762 7.53e+01 1.663e+05 1.272e+02 1.754e+01\r", " 763 7.54e+01 1.663e+05 1.249e+02 1.752e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 764 7.54e+01 1.663e+05 1.234e+02 1.748e+01\r", " 765 7.55e+01 1.663e+05 1.226e+02 1.744e+01\r", " 766 7.56e+01 1.662e+05 1.225e+02 1.743e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 767 7.57e+01 1.662e+05 1.232e+02 1.739e+01\r", " 768 7.57e+01 1.662e+05 1.245e+02 1.735e+01\r", " 769 7.58e+01 1.661e+05 1.265e+02 1.731e+01\r", " 770 7.58e+01 1.661e+05 1.291e+02 1.728e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 771 7.59e+01 1.661e+05 1.321e+02 1.724e+01\r", " 772 7.60e+01 1.661e+05 1.356e+02 1.721e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 773 7.61e+01 1.660e+05 1.394e+02 1.717e+01\r", " 774 7.63e+01 1.660e+05 1.434e+02 1.712e+01\r", " 775 7.64e+01 1.660e+05 1.476e+02 1.706e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 776 7.65e+01 1.659e+05 1.519e+02 1.701e+01\r", " 777 7.66e+01 1.659e+05 1.562e+02 1.696e+01\r", " 778 7.66e+01 1.659e+05 1.604e+02 1.690e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 779 7.67e+01 1.659e+05 1.646e+02 1.685e+01\r", " 780 7.68e+01 1.658e+05 1.686e+02 1.677e+01\r", " 781 7.69e+01 1.658e+05 1.724e+02 1.668e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 782 7.70e+01 1.658e+05 1.760e+02 1.659e+01\r", " 783 7.70e+01 1.657e+05 1.793e+02 1.650e+01\r", " 784 7.71e+01 1.657e+05 1.824e+02 1.637e+01\r", " 785 7.72e+01 1.657e+05 1.851e+02 1.626e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 786 7.72e+01 1.657e+05 1.874e+02 1.617e+01\r", " 787 7.73e+01 1.657e+05 1.894e+02 1.608e+01\r", " 788 7.73e+01 1.656e+05 1.909e+02 1.599e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 789 7.75e+01 1.656e+05 1.921e+02 1.590e+01\r", " 790 7.77e+01 1.656e+05 1.928e+02 1.581e+01\r", " 791 7.78e+01 1.656e+05 1.932e+02 1.573e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 792 7.79e+01 1.656e+05 1.931e+02 1.565e+01\r", " 793 7.79e+01 1.656e+05 1.925e+02 1.556e+01\r", " 794 7.80e+01 1.656e+05 1.916e+02 1.548e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 795 7.81e+01 1.656e+05 1.902e+02 1.540e+01\r", " 796 7.82e+01 1.655e+05 1.885e+02 1.533e+01\r", " 797 7.82e+01 1.655e+05 1.863e+02 1.526e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 798 7.83e+01 1.655e+05 1.838e+02 1.518e+01\r", " 799 7.84e+01 1.655e+05 1.809e+02 1.512e+01\r", " 800 7.85e+01 1.655e+05 1.776e+02 1.507e+01\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 801 7.85e+01 1.655e+05 1.741e+02 1.499e+01\r", " 802 7.86e+01 1.655e+05 1.702e+02 1.494e+01\r", " 803 7.87e+01 1.655e+05 1.661e+02 1.489e+01\r", " 804 7.87e+01 1.655e+05 1.617e+02 1.482e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 805 7.88e+01 1.655e+05 1.572e+02 1.474e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 806 7.90e+01 1.655e+05 1.524e+02 1.466e+01\r", " 807 7.91e+01 1.655e+05 1.476e+02 1.458e+01\r", " 808 7.92e+01 1.655e+05 1.426e+02 1.451e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 809 7.93e+01 1.656e+05 1.376e+02 1.445e+01\r", " 810 7.94e+01 1.656e+05 1.326e+02 1.437e+01\r", " 811 7.94e+01 1.656e+05 1.277e+02 1.429e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 812 7.95e+01 1.656e+05 1.228e+02 1.422e+01\r", " 813 7.96e+01 1.656e+05 1.181e+02 1.415e+01\r", " 814 7.96e+01 1.656e+05 1.137e+02 1.408e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 815 7.97e+01 1.656e+05 1.095e+02 1.401e+01\r", " 816 7.98e+01 1.656e+05 1.057e+02 1.394e+01\r", " 817 7.99e+01 1.656e+05 1.023e+02 1.387e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 818 7.99e+01 1.656e+05 9.949e+01 1.381e+01\r", " 819 8.00e+01 1.656e+05 9.721e+01 1.373e+01\r", " 820 8.01e+01 1.656e+05 9.556e+01 1.366e+01\r", " 821 8.01e+01 1.656e+05 9.456e+01 1.359e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 822 8.02e+01 1.656e+05 9.423e+01 1.352e+01\r", " 823 8.04e+01 1.656e+05 9.457e+01 1.344e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 824 8.05e+01 1.656e+05 9.555e+01 1.337e+01\r", " 825 8.06e+01 1.656e+05 9.711e+01 1.330e+01\r", " 826 8.07e+01 1.656e+05 9.920e+01 1.322e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 827 8.07e+01 1.656e+05 1.017e+02 1.314e+01\r", " 828 8.08e+01 1.657e+05 1.046e+02 1.307e+01\r", " 829 8.09e+01 1.657e+05 1.078e+02 1.301e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 830 8.10e+01 1.657e+05 1.113e+02 1.294e+01\r", " 831 8.11e+01 1.657e+05 1.148e+02 1.288e+01\r", " 832 8.11e+01 1.657e+05 1.184e+02 1.282e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 833 8.12e+01 1.657e+05 1.220e+02 1.276e+01\r", " 834 8.13e+01 1.657e+05 1.256e+02 1.269e+01\r", " 835 8.13e+01 1.657e+05 1.291e+02 1.264e+01\r", " 836 8.14e+01 1.657e+05 1.324e+02 1.259e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 837 8.15e+01 1.657e+05 1.355e+02 1.253e+01\r", " 838 8.15e+01 1.657e+05 1.385e+02 1.249e+01\r", " 839 8.16e+01 1.657e+05 1.411e+02 1.244e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 840 8.18e+01 1.657e+05 1.436e+02 1.239e+01\r", " 841 8.19e+01 1.657e+05 1.457e+02 1.234e+01\r", " 842 8.20e+01 1.657e+05 1.475e+02 1.228e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 843 8.21e+01 1.657e+05 1.491e+02 1.223e+01\r", " 844 8.21e+01 1.657e+05 1.502e+02 1.219e+01\r", " 845 8.22e+01 1.658e+05 1.511e+02 1.213e+01\r", " 846 8.22e+01 1.658e+05 1.516e+02 1.210e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 847 8.23e+01 1.658e+05 1.518e+02 1.207e+01\r", " 848 8.24e+01 1.658e+05 1.517e+02 1.205e+01\r", " 849 8.24e+01 1.658e+05 1.512e+02 1.204e+01\r", " 850 8.25e+01 1.658e+05 1.504e+02 1.204e+01\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 851 8.25e+01 1.658e+05 1.493e+02 1.203e+01\r", " 852 8.26e+01 1.658e+05 1.479e+02 1.203e+01\r", " 853 8.26e+01 1.658e+05 1.462e+02 1.203e+01\r", " 854 8.27e+01 1.658e+05 1.441e+02 1.202e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 855 8.28e+01 1.658e+05 1.419e+02 1.201e+01\r", " 856 8.28e+01 1.658e+05 1.393e+02 1.202e+01\r", " 857 8.29e+01 1.659e+05 1.366e+02 1.203e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 858 8.30e+01 1.659e+05 1.336e+02 1.204e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 859 8.32e+01 1.659e+05 1.304e+02 1.205e+01\r", " 860 8.33e+01 1.659e+05 1.271e+02 1.204e+01\r", " 861 8.34e+01 1.659e+05 1.237e+02 1.204e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 862 8.34e+01 1.659e+05 1.201e+02 1.203e+01\r", " 863 8.35e+01 1.659e+05 1.165e+02 1.203e+01\r", " 864 8.36e+01 1.659e+05 1.128e+02 1.202e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 865 8.37e+01 1.659e+05 1.091e+02 1.201e+01\r", " 866 8.37e+01 1.659e+05 1.055e+02 1.200e+01\r", " 867 8.38e+01 1.659e+05 1.020e+02 1.200e+01\r", " 868 8.38e+01 1.659e+05 9.853e+01 1.199e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 869 8.39e+01 1.659e+05 9.528e+01 1.198e+01\r", " 870 8.40e+01 1.659e+05 9.225e+01 1.197e+01\r", " 871 8.40e+01 1.659e+05 8.948e+01 1.196e+01\r", " 872 8.41e+01 1.659e+05 8.703e+01 1.195e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 873 8.41e+01 1.659e+05 8.494e+01 1.194e+01\r", " 874 8.42e+01 1.659e+05 8.324e+01 1.193e+01\r", " 875 8.43e+01 1.659e+05 8.196e+01 1.193e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 876 8.44e+01 1.658e+05 8.113e+01 1.193e+01\r", " 877 8.46e+01 1.658e+05 8.074e+01 1.193e+01\r", " 878 8.47e+01 1.658e+05 8.079e+01 1.192e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 879 8.48e+01 1.658e+05 8.128e+01 1.191e+01\r", " 880 8.49e+01 1.658e+05 8.215e+01 1.190e+01\r", " 881 8.49e+01 1.658e+05 8.339e+01 1.189e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 882 8.50e+01 1.658e+05 8.493e+01 1.188e+01\r", " 883 8.51e+01 1.658e+05 8.673e+01 1.186e+01\r", " 884 8.51e+01 1.658e+05 8.875e+01 1.184e+01\r", " 885 8.52e+01 1.658e+05 9.091e+01 1.183e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 886 8.52e+01 1.657e+05 9.319e+01 1.181e+01\r", " 887 8.53e+01 1.657e+05 9.552e+01 1.180e+01\r", " 888 8.53e+01 1.657e+05 9.787e+01 1.178e+01\r", " 889 8.54e+01 1.657e+05 1.002e+02 1.177e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 890 8.55e+01 1.657e+05 1.025e+02 1.176e+01\r", " 891 8.55e+01 1.657e+05 1.046e+02 1.175e+01\r", " 892 8.56e+01 1.657e+05 1.067e+02 1.174e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 893 8.56e+01 1.656e+05 1.086e+02 1.172e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 894 8.58e+01 1.656e+05 1.103e+02 1.171e+01\r", " 895 8.60e+01 1.656e+05 1.119e+02 1.170e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 896 8.61e+01 1.656e+05 1.133e+02 1.167e+01\r", " 897 8.62e+01 1.656e+05 1.145e+02 1.166e+01\r", " 898 8.62e+01 1.656e+05 1.154e+02 1.164e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 899 8.63e+01 1.656e+05 1.162e+02 1.162e+01\r", " 900 8.64e+01 1.655e+05 1.167e+02 1.160e+01\n", " 901 8.65e+01 1.655e+05 1.169e+02 1.159e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 902 8.65e+01 1.655e+05 1.170e+02 1.158e+01\r", " 903 8.66e+01 1.655e+05 1.168e+02 1.155e+01\r", " 904 8.67e+01 1.655e+05 1.164e+02 1.152e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 905 8.67e+01 1.655e+05 1.158e+02 1.149e+01\r", " 906 8.68e+01 1.655e+05 1.149e+02 1.143e+01\r", " 907 8.69e+01 1.655e+05 1.139e+02 1.138e+01\r", " 908 8.69e+01 1.655e+05 1.127e+02 1.134e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 909 8.70e+01 1.655e+05 1.112e+02 1.129e+01\r", " 910 8.70e+01 1.655e+05 1.096e+02 1.123e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 911 8.72e+01 1.654e+05 1.079e+02 1.118e+01\r", " 912 8.74e+01 1.654e+05 1.060e+02 1.113e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 913 8.75e+01 1.654e+05 1.039e+02 1.108e+01\r", " 914 8.75e+01 1.654e+05 1.018e+02 1.103e+01\r", " 915 8.76e+01 1.654e+05 9.954e+01 1.097e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 916 8.77e+01 1.654e+05 9.722e+01 1.091e+01\r", " 917 8.77e+01 1.654e+05 9.484e+01 1.086e+01\r", " 918 8.78e+01 1.654e+05 9.243e+01 1.081e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 919 8.79e+01 1.654e+05 9.002e+01 1.076e+01\r", " 920 8.80e+01 1.654e+05 8.761e+01 1.071e+01\r", " 921 8.80e+01 1.654e+05 8.525e+01 1.066e+01\r", " 922 8.81e+01 1.654e+05 8.296e+01 1.061e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 923 8.81e+01 1.654e+05 8.077e+01 1.055e+01\r", " 924 8.83e+01 1.654e+05 7.870e+01 1.049e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 925 8.84e+01 1.654e+05 7.678e+01 1.044e+01\r", " 926 8.85e+01 1.654e+05 7.504e+01 1.037e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 927 8.86e+01 1.654e+05 7.351e+01 1.031e+01\r", " 928 8.88e+01 1.654e+05 7.220e+01 1.025e+01\r", " 929 8.89e+01 1.654e+05 7.114e+01 1.019e+01\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 930 8.89e+01 1.653e+05 7.034e+01 1.013e+01\r", " 931 8.90e+01 1.653e+05 6.980e+01 1.007e+01\r", " 932 8.90e+01 1.653e+05 6.952e+01 1.000e+01\r", " 933 8.91e+01 1.653e+05 6.951e+01 9.935e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 934 8.91e+01 1.653e+05 6.974e+01 9.868e+00\r", " 935 8.92e+01 1.653e+05 7.021e+01 9.798e+00\r", " 936 8.92e+01 1.653e+05 7.089e+01 9.733e+00\r", " 937 8.93e+01 1.653e+05 7.176e+01 9.664e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 938 8.93e+01 1.653e+05 7.278e+01 9.601e+00\r", " 939 8.94e+01 1.653e+05 7.394e+01 9.544e+00\r", " 940 8.94e+01 1.653e+05 7.519e+01 9.491e+00\r", " 941 8.95e+01 1.653e+05 7.651e+01 9.440e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 942 8.95e+01 1.653e+05 7.788e+01 9.390e+00\r", " 943 8.96e+01 1.653e+05 7.927e+01 9.337e+00\r", " 944 8.96e+01 1.653e+05 8.065e+01 9.293e+00\r", " 945 8.97e+01 1.653e+05 8.200e+01 9.253e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 946 8.98e+01 1.653e+05 8.331e+01 9.215e+00\r", " 947 8.98e+01 1.653e+05 8.455e+01 9.182e+00\r", " 948 8.99e+01 1.653e+05 8.571e+01 9.148e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 949 9.01e+01 1.654e+05 8.677e+01 9.107e+00\r", " 950 9.02e+01 1.654e+05 8.773e+01 9.069e+00\n", " 951 9.03e+01 1.654e+05 8.858e+01 9.029e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 952 9.03e+01 1.654e+05 8.930e+01 9.001e+00\r", " 953 9.04e+01 1.654e+05 8.988e+01 8.967e+00\r", " 954 9.04e+01 1.654e+05 9.034e+01 8.941e+00\r", " 955 9.05e+01 1.654e+05 9.065e+01 8.918e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 956 9.05e+01 1.654e+05 9.081e+01 8.897e+00\r", " 957 9.06e+01 1.654e+05 9.083e+01 8.878e+00\r", " 958 9.07e+01 1.654e+05 9.070e+01 8.861e+00\r", " 959 9.07e+01 1.654e+05 9.043e+01 8.847e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 960 9.08e+01 1.654e+05 9.001e+01 8.832e+00\r", " 961 9.08e+01 1.654e+05 8.945e+01 8.818e+00\r", " 962 9.09e+01 1.654e+05 8.875e+01 8.800e+00\r", " 963 9.09e+01 1.654e+05 8.792e+01 8.779e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 964 9.10e+01 1.654e+05 8.697e+01 8.761e+00\r", " 965 9.10e+01 1.654e+05 8.590e+01 8.742e+00\r", " 966 9.11e+01 1.654e+05 8.472e+01 8.723e+00\r", " 967 9.12e+01 1.654e+05 8.344e+01 8.703e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 968 9.12e+01 1.655e+05 8.207e+01 8.685e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 969 9.15e+01 1.655e+05 8.062e+01 8.667e+00\r", " 970 9.16e+01 1.655e+05 7.910e+01 8.648e+00\r", " 971 9.17e+01 1.655e+05 7.752e+01 8.633e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 972 9.17e+01 1.655e+05 7.590e+01 8.618e+00\r", " 973 9.18e+01 1.655e+05 7.426e+01 8.603e+00\r", " 974 9.19e+01 1.655e+05 7.260e+01 8.581e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 975 9.20e+01 1.655e+05 7.096e+01 8.563e+00\r", " 976 9.20e+01 1.655e+05 6.933e+01 8.544e+00\r", " 977 9.21e+01 1.655e+05 6.775e+01 8.521e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 978 9.22e+01 1.655e+05 6.623e+01 8.502e+00\r", " 979 9.22e+01 1.655e+05 6.480e+01 8.482e+00\r", " 980 9.23e+01 1.655e+05 6.346e+01 8.467e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 981 9.24e+01 1.655e+05 6.223e+01 8.448e+00\r", " 982 9.24e+01 1.655e+05 6.115e+01 8.428e+00\r", " 983 9.25e+01 1.655e+05 6.021e+01 8.407e+00\r", " 984 9.26e+01 1.655e+05 5.943e+01 8.388e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 985 9.26e+01 1.655e+05 5.882e+01 8.368e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 986 9.28e+01 1.655e+05 5.839e+01 8.350e+00\r", " 987 9.29e+01 1.655e+05 5.813e+01 8.329e+00\r", " 988 9.30e+01 1.655e+05 5.806e+01 8.306e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 989 9.31e+01 1.655e+05 5.815e+01 8.289e+00\r", " 990 9.31e+01 1.655e+05 5.841e+01 8.280e+00\r", " 991 9.32e+01 1.655e+05 5.881e+01 8.267e+00\r", " 992 9.33e+01 1.655e+05 5.935e+01 8.256e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 993 9.33e+01 1.655e+05 6.001e+01 8.243e+00\r", " 994 9.34e+01 1.655e+05 6.076e+01 8.235e+00\r", " 995 9.35e+01 1.655e+05 6.160e+01 8.218e+00\r", " 996 9.35e+01 1.655e+05 6.250e+01 8.206e+00\r" ] }, { "name": "stdout", "output_type": "stream", "text": [ " 997 9.36e+01 1.655e+05 6.345e+01 8.193e+00\r", " 998 9.37e+01 1.655e+05 6.442e+01 8.180e+00\r", " 999 9.37e+01 1.655e+05 6.540e+01 8.168e+00\r\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "TV Restruction\n", "SNR: 22.04 (dB), MAE: 0.009\n" ] } ], "source": [ "print(f\"Solving on {device_info()}\\n\")\n", "tangle_recon = solver.solve()\n", "\n", "print(\n", " \"TV Restruction\\nSNR: %.2f (dB), MAE: %.3f\"\n", " % (metric.snr(tangle, tangle_recon), metric.mae(tangle, tangle_recon))\n", ")" ] }, { "cell_type": "markdown", "id": "15c83b02", "metadata": {}, "source": [ "Show the recovered image." ] }, { "cell_type": "code", "execution_count": 5, "id": "42aa5bce", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2025-08-19T21:06:08.556488Z", "iopub.status.busy": "2025-08-19T21:06:08.556293Z", "iopub.status.idle": "2025-08-19T21:06:09.002617Z", "shell.execute_reply": "2025-08-19T21:06:09.001811Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmsAAAHuCAYAAADA2UsSAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAACP+UlEQVR4nO3dd3xTVf8H8M/N7N6DXUoBQYbI3lPhAURxgYshiCgqorhRUR4VxYUDUH8Pw4GKe4KIgqCCCIqyBdmjrAJt6cw4vz9ubkjapE3SpLlJPu/nladyc3PvyWjut99zzvdIQggBIiIiIlIlTbAbQERERETuMVgjIiIiUjEGa0REREQqxmCNiIiISMUYrBERERGpGIM1IiIiIhVjsEZERESkYgzWiIiIiFSMwRoRERGRijFYI6JaI0mSR7dXXnkFkiThu+++c3us//u//4MkSfjss8/c7vPEE084HVev16NRo0aYMGECjh07FoinGFTPPPMMvvjii6C2Ye7cuVi0aFGl7fv374ckSS7vqy3jxo3Df/7zn6CdX7F9+3Y88cQT2L9/f0COr3zu/aFv377o27ev0zZJkvDEE0/45fgVmUwm5OTkYPbs2QE5fqjSBbsBRBQ51q1b5/Tv//73v1i1ahVWrlzptL1u3bp48MEHsWDBArcX14ULFyI9PR3Dhg2r9rzfffcdEhMTce7cOXz//fd48cUXsXbtWvz111/Q6/W+PyGVeeaZZ3DNNddg+PDhQWvD3LlzkZaWhrFjxzptr1u3LtatW4ecnJygtGvTpk14++23sX79+qCc39H27dvx5JNPom/fvmjcuHGwm+O1devWoUGDBgE5tl6vx+OPP4577rkHo0aNQmpqakDOE2oYrBFRrenatavTv9PT06HRaCptB4ArrrgCX3zxBfLy8ip9Ye/cuRPr1q3D1KlTPQq2OnTogLS0NADAJZdcglOnTmHhwoX45Zdf0K9fvxo8o9BVUlKCqKgov2VgqmM0Gl2+z7Xl2WefRefOndGxY8egtcFXxcXFiImJCXYz7AL9Pl5//fW499578eabb+KRRx4J6LlCBbtBiUiVxo8fj/Lycrz//vuV7lu4cCEAuVvLF8oF+/jx407bf/jhBwwYMAAJCQmIiYlBjx498OOPP1Z6/M6dO3H99dcjMzMTRqMRjRo1wujRo1FWVmbfZ+vWrbjiiiuQnJyMqKgotGvXDm+//bbTcX766SdIkoQPPvgA06ZNQ7169ZCQkIBLLrkE//zzj9O+mzZtwmWXXYaMjAwYjUbUq1cPQ4cOxeHDhwHIXVNFRUV4++237d2+SvfVokWLIEkSvv/+e4wbNw7p6emIiYlBWVkZxo4d6zK746orzWq14rXXXkO7du0QHR2NpKQkdO3aFV999RUAoHHjxti2bRtWr15tb4NybHfdoL/88gsGDBiA+Ph4xMTEoHv37vj222+d9lHav2rVKtx+++1IS0tDamoqrrrqKhw9erRS2ys6fvw4Pv/8c4waNarSfWfPnsXUqVPRpEkTGI1GZGRkYMiQIdi5c6d9n/Lycjz11FNo0aIFjEYj0tPTcfPNN+PkyZNOx2rcuDEuu+wyfPfdd2jfvj2io6PRokULLFiwwOm5XHvttQCAfv362V8n5XXp27cvWrdujTVr1qB79+6IiYmxf86XLFmCgQMHom7duoiOjkbLli3x0EMPoaioqNrXwJW9e/fiuuuuQ7169WA0GpGZmYkBAwbgr7/+qvJxrrpBjxw5gltvvRUNGzaEwWBAvXr1cM011zj9jhUUFOC+++5DdnY2DAYD6tevjylTplRqv8FgwMiRI/HWW29BCOHTcws3zKwRkSpdcsklyMrKwoIFC3DXXXfZt1ssFrz77rvo2rUrLrzwQp+OvW/fPgBA8+bN7dvee+89jB49GldccQXefvtt6PV6vPnmmxg0aBCWL1+OAQMGAAD+/vtv9OzZE2lpaZgxYwaaNWuG3NxcfPXVVygvL4fRaMQ///yD7t27IyMjA6+++ipSU1Px3nvvYezYsTh+/DgeeOABp/Y88sgj6NGjB/73v/+hoKAADz74IIYNG4YdO3ZAq9WiqKgIl156KbKzszFnzhxkZmbi2LFjWLVqFQoLCwHIXVP9+/dHv3798NhjjwEAEhISnM4zbtw4DB06FO+++y6Kioq87gIeO3Ys3nvvPYwfPx4zZsyAwWDAn3/+aR979fnnn+Oaa65BYmIi5s6dC0DOqLmzevVqXHrppWjbti3mz58Po9GIuXPnYtiwYfjggw8wcuRIp/1vueUWDB06FO+//z4OHTqE+++/HzfddFOlbvSKvv/+e5hMpkpZ1MLCQvTs2RP79+/Hgw8+iC5duuDcuXNYs2YNcnNz0aJFC1itVlxxxRX4+eef8cADD6B79+44cOAApk+fjr59+2Ljxo2Ijo62H/Pvv//G1KlT8dBDDyEzMxP/+9//MH78eDRt2hS9e/fG0KFD8cwzz+CRRx7BnDlz0L59ewBw6h7Ozc3FTTfdhAceeADPPPMMNBo5r7J7924MGTIEU6ZMQWxsLHbu3InnnnsOv//+e7WvgStDhgyBxWLBrFmz0KhRI5w6dQpr167F2bNnvTrOkSNH0KlTJ5hMJjzyyCNo27Yt8vLysHz5cpw5cwaZmZkoLi5Gnz59cPjwYfs+27Ztw+OPP44tW7bghx9+cPrjoG/fvpg3bx62bt2KNm3aeP3cwo4gIgqSMWPGiNjYWLf3T58+XQAQf/75p33b119/LQCI//u//6v2+Mrjjx07Jkwmkzhz5oz46KOPRGxsrLj++uvt+xUVFYmUlBQxbNgwp8dbLBZx0UUXic6dO9u39e/fXyQlJYkTJ064Pe91110njEajOHjwoNP2wYMHi5iYGHH27FkhhBCrVq0SAMSQIUOc9vvoo48EALFu3TohhBAbN24UAMQXX3xR5fONjY0VY8aMqbR94cKFAoAYPXp0pfvGjBkjsrKyKm1XXjvFmjVrBAAxbdq0KtvQqlUr0adPn0rb9+3bJwCIhQsX2rd17dpVZGRkiMLCQvs2s9ksWrduLRo0aCCsVqtT+ydNmuR0zFmzZgkAIjc3t8o23X777SI6Otp+PMWMGTMEALFixQq3j/3ggw8EAPHpp586bd+wYYMAIObOnWvflpWVJaKiosSBAwfs20pKSkRKSoqYOHGifdvHH38sAIhVq1ZVOl+fPn0EAPHjjz9W+ZysVqswmUxi9erVAoD4+++/7fdVfO9cOXXqlAAgZs+eXeV+ffr0qfR+AhDTp0+3/3vcuHFCr9eL7du3uz3OzJkzhUajERs2bHDa/sknnwgAYunSpU7bd+/eLQCIefPmVdm+SMFuUCJSrZtvvhkajcapG2nhwoWIjY2tlHWpSp06daDX65GcnIwRI0agQ4cOTl2Sa9euxenTpzFmzBiYzWb7zWq14j//+Q82bNiAoqIiFBcXY/Xq1RgxYgTS09Pdnm/lypUYMGAAGjZs6LR97NixKC4urjTR4vLLL3f6d9u2bQEABw4cAAA0bdoUycnJePDBB/HGG29g+/btHj93R1dffbVPjwOAZcuWAQDuuOMOn4/hqKioCOvXr8c111yDuLg4+3atVotRo0bh8OHDlbqCq3ud3Dl69CjS09MrdesuW7YMzZs3xyWXXOL2sd988w2SkpIwbNgwp89Gu3btUKdOHfz0009O+7dr1w6NGjWy/zsqKgrNmzevto2OkpOT0b9//0rb9+7dixtuuAF16tSBVquFXq9Hnz59AAA7duzw+PgAkJKSgpycHDz//PN46aWXsGnTJlitVq+OoVi2bBn69euHli1but3nm2++QevWrdGuXTun13HQoEGQJKnS65iRkQFAztoRx6wRkYplZWVhwIABeP/991FWVoZTp07hm2++wbXXXov4+HiPj/PDDz9gw4YNWL58Oa6++mqsWbPGqWtVGVdzzTXXQK/XO92ee+45CCFw+vRpnDlzBhaLpdqZcHl5eahbt26l7fXq1bPf76jiBAql67CkpAQAkJiYiNWrV6Ndu3Z45JFH0KpVK9SrVw/Tp0+HyWTy+HVw1SZPnTx5ElqtFnXq1PH5GI7OnDkDIYRfXyd3lMkUFZ08ebLa9/L48eM4e/YsDAZDpc/GsWPHcOrUqSrbqLSzujY6cvWanDt3Dr169cL69evx1FNP4aeffsKGDRvspWu8OT4gjzv78ccfMWjQIMyaNQvt27dHeno6Jk+ebO9a95Snr+PmzZsrvYbx8fEQQlR6HZX3y9vnFa44Zo2IVG38+PFYsWIFvvzySxw9ehTl5eUYP368V8e46KKL7LNBL730UgwaNAhvvfUWxo8fj06dOtnve+2119zOdMvMzITFYoFWq7UP6ncnNTUVubm5lbYrg+GV83mjTZs2+PDDDyGEwObNm7Fo0SLMmDED0dHReOihhzw6hquZn1FRUU4TIxQVL57p6emwWCw4duxYjYI+RXJyMjQajd9fJ1fS0tLw559/Vtqenp5e7XupTGZwV/PPmz8aPOXqfVq5ciWOHj2Kn376yZ5NA+D1+DJHWVlZmD9/PgBg165d+Oijj/DEE0+gvLwcb7zxhsfH8fR1jI6OdsqSV7zf0enTp11uj1TMrBGRqg0fPhypqalYsGABFi5ciObNm6Nnz54+H0+SJMyZMwdarRaPPvooAKBHjx5ISkrC9u3b0bFjR5c3g8GA6Oho9OnTBx9//HGlYMbRgAED7BdXR++88w5iYmJqVPpAkiRcdNFFePnll5GUlOQUhHibwQHkGYwnTpxwmrVXXl6O5cuXO+03ePBgAMC8efOqPJ6nbYiNjUWXLl3w2WefOe1vtVrx3nvvoUGDBk4TQGqiRYsWyMvLQ35+vtP2wYMHY9euXVUOzr/sssuQl5cHi8Xi8nNxwQUXeN0eTzOCjpQAruKEjTfffNPr87vSvHlzPProo2jTpo3LwLYqgwcPxqpVqyp1Wzu67LLLsGfPHqSmprp8HSvOSN67dy8A+DyJKNwws0ZEqmY0GnHjjTfitddegxACzz77bI2P2axZM9x6662YO3cufvnlF/Ts2ROvvfYaxowZg9OnT+Oaa65BRkYGTp48ib///hsnT560BykvvfQSevbsiS5duuChhx5C06ZNcfz4cXz11Vd48803ER8fj+nTp+Obb75Bv3798PjjjyMlJQWLFy/Gt99+i1mzZiExMdGr9n7zzTeYO3cuhg8fjiZNmkAIgc8++wxnz57FpZdeat+vTZs2+Omnn/D111+jbt26iI+PrzaYGDlyJB5//HFcd911uP/++1FaWopXX30VFovFab9evXph1KhReOqpp3D8+HFcdtllMBqN2LRpE2JiYuzdykoGcMmSJWjSpAmioqLczuabOXMmLr30UvTr1w/33XcfDAYD5s6di61bt+KDDz7waxV+IQTWr1+PgQMH2rdPmTIFS5YswRVXXIGHHnoInTt3RklJCVavXo3LLrsM/fr1w3XXXYfFixdjyJAhuPvuu9G5c2fo9XocPnwYq1atwhVXXIErr7zSq/a0bt0aAPDWW28hPj4eUVFRyM7OrrIAbPfu3ZGcnIzbbrsN06dPh16vx+LFi/H333/79Jps3rwZd955J6699lo0a9YMBoMBK1euxObNmz3O1CpmzJiBZcuWoXfv3njkkUfQpk0bnD17Ft999x3uvfdetGjRAlOmTMGnn36K3r1745577kHbtm1htVpx8OBBfP/995g6dSq6dOliP+Zvv/0GrVaL3r17+/T8wk4QJzcQUYSrbjao4u+//xYAhFarFUePHvX4+MqsuJMnT1a67/jx4yIuLk7069fPvm316tVi6NChIiUlRej1elG/fn0xdOhQ8fHHHzs9dvv27eLaa68VqampwmAwiEaNGomxY8eK0tJS+z5btmwRw4YNE4mJicJgMIiLLrrIaSakEOdng1Y8fsWZkzt37hTXX3+9yMnJEdHR0SIxMVF07txZLFq0yOlxf/31l+jRo4eIiYkRAOyz+JTZlBVn4imWLl0q2rVrJ6Kjo0WTJk3E66+/7nJGocViES+//LJo3bq1MBgMIjExUXTr1k18/fXX9n32798vBg4cKOLj4wUA+0xTV7NBhRDi559/Fv379xexsbEiOjpadO3a1el4VbVfef1czaqs2O7GjRtXmk0qhBBnzpwRd999t2jUqJHQ6/UiIyNDDB06VOzcudO+j8lkEi+88IK46KKLRFRUlIiLixMtWrQQEydOFLt377bvl5WVJYYOHVrpHK5mVM6ePVtkZ2cLrVbr9Lr06dNHtGrVyuXzWLt2rejWrZuIiYkR6enp4pZbbhF//vlnpdfVk9mgx48fF2PHjhUtWrQQsbGxIi4uTrRt21a8/PLLwmw2V9l2VJgNKoQQhw4dEuPGjRN16tQRer1e1KtXT4wYMUIcP37cvs+5c+fEo48+Ki644AL756dNmzbinnvuEceOHXM6Xq9evSrNzo5kkhCsOEdEROHtxRdfxNNPP40jR4441UUj9dmzZw+aNWuG5cuXO2WOIxmDNSIiCnulpaVo2bIl7rjjDtx3333Bbg5V4eabb8bhw4exYsWKYDdFNTjBgIiIwl5UVBTefffdKldUoOAzm83IycnBnDlzgt0UVWFmjYiIiEjFmFkjIiIiUjEGa0REREQq5lOwtnnzZowfPx45OTmIjo5GdHQ0mjVrhokTJ2Ljxo3+bmOtkiQJTzzxRJX7HD16FE888QT++uuvgLShquOPHTvWaR09X5hMJrRo0cIv9apqaunSpdW+3jXRt29f9O3b1y/HqvjZ+Omnn1yuaecvu3btgsFg8LpAZSRYv349rrzySjRq1AhGoxGZmZno1q0bpk6d6rRf3759IUkS/vOf/1Q6xv79+yFJEl544QX7NuU9VW5arRbp6ekYNmxYjb7bdu3ahfvuuw8dOnRAUlISUlJS0KNHD3zyySeV9v3ss89w/fXXo2nTpoiOjkbjxo1x4403Yvfu3R6fb+/evbjqqquQlJSEuLg4XHrppdV+jo4fP47U1FRIkuSyXe78/PPPMBqNHq192bhxY4wdO9b+74qvtyRJSE5ORpcuXZzWbgXk762cnBzMnj3b47a5onwmlHp1Fa1Zs8belkWLFrk8xquvvgpJkuz10lyp+Lwcb46vgbdOnDiBsWPHIi0tDTExMejWrRt+/PFHjx/vzWfjww8/RLt27RAVFYV69ephypQpOHfuXKX9fv/9dwwaNAjx8fGIi4tDv3798Ouvv/r8HKkyr4O1N998Ex06dMD69etx991345tvvsG3336LKVOmYNu2bejUqRP27NkTiLaqxtGjR/Hkk08GNFgL5PHnzp2LM2fOOK2NGCxLly7Fk08+Gexm+KR9+/ZYt24d2rdvH5DjN2/eHDfeeCPuueeegBw/VH377bfo3r07CgoKMGvWLHz//fd45ZVX0KNHDyxZssTlY5YvX15llfqKnnnmGaxbtw4//fQTHnvsMaxduxZ9+vTxKmBy9P333+Pbb7/F1VdfjY8//hiLFy9Gs2bNcO2112LGjBlO+z733HMoLi7GtGnT8N133+Gpp57Cpk2b0L59e2zbtq3ac508eRK9evXCrl27sGDBAnz00UcoLS1F3759q6wwf8cdd7hcP7MqQghMmTIFEyZMQFZWllePdaS83uvWrcO7776LrKwsjB07Fq+99pp9H71ej8cffxwzZsyotGaot+Lj47Fv3z6Xn4kFCxYgISGhyscrSyZt27YN69evd7vfNddcY39ejrfHHnvMp3aXlZVhwIAB+PHHH/HKK6/gyy+/RGZmJv7zn/9g9erV1T7em8/G4sWLcf3116NTp05YtmwZpk+fjkWLFuGqq65y2m/Dhg3o3bs3SkpK8O677+Ldd99FaWkpBgwYgHXr1vn0PMkFb4qy/fLLL0Kj0Yhhw4aJsrIyl/t89NFH4siRI1Uep6ioyJvT1iq4KPZX0YYNG1wWd3TH2+db1fE9LSLqjslkEvXr1xcPPfSQz8fwpzvuuKPa4o0Kq9UqiouLvTq+q4KOvvLks+FvGzduFADEr7/+WqvnVbPevXuLnJwcYTKZKt1nsVic/t2nTx/RvHlz0aRJE9GhQwdhtVrt9ylFWp9//nn7NndFat9++20BQDz++OM+tfnkyZNO51YMHTpUxMTEOBXTdSwiqjhy5IjQ6/Vi/Pjx1Z7r/vvvF3q9Xuzfv9++LT8/X6SlpYkRI0a4fMwnn3wi4uLi7M+z4vN3Z+nSpQKAUwHZqmRlZYkxY8bY/+3u9VaK2Hbr1s1pe1lZmUhJSRFPP/20R+dzRSk627VrV3HDDTc43VdQUCBiYmLEhAkT3H4HK9/PQ4cOFQDEhAkTXJ4HgLjjjjt8bqcrc+bMEQDE2rVr7dtMJpO48MILRefOnat9vKefDbPZLOrWrSsGDhzo9PjFixcLAGLp0qX2bYMGDRKZmZlO17mCggKRlpYmunfv7tPzpMq8yqw988wz0Gq1ePPNN2EwGFzuc+2116JevXr2fyvddlu2bMHAgQMRHx+PAQMGAJAXap00aRLq168Pg8GAJk2aYNq0aU6LCitdFa7S0RW7pZ544glIkoRt27bh+uuvR2JiIjIzMzFu3LhKa8IVFBRgwoQJSE1NRVxcHP7zn/9g165d1b4GP/30Ezp16gRArgWjpLWVdlT1fCt2ASgcu+qqO77i33//xZAhQxAXF4eGDRti6tSpLhdjruirr77CkSNHMGrUqEr37dy5E9dffz0yMzNhNBrRqFEjjB492um4x44dw8SJE9GgQQMYDAZkZ2fjySefhNlstu/j2L300ksvITs7G3FxcejWrRt+++03+35jx461T8927CLYv3+/fdudd96JN954Ay1btoTRaLR3jTz55JPo0qULUlJSkJCQgPbt22P+/PkuuzU8sXLlSvTt2xepqamIjo5Go0aNcPXVV6O4uNjtY9x1g65fvx7Dhg1DamoqoqKikJOTgylTpjjts3v3btxwww3IyMiA0WhEy5YtXU5V79ChA1q2bOnVosrhLi8vD2lpadDpKq+Wp9FU/krT6/V4+umn8ccff7jNvFWnY8eOAOC0fqY30tLSXC6d1LlzZxQXF9sXrQaAjIyMSvvVq1cPDRo0wKFDh6o91+eff47+/fs7ZboSEhJw1VVX4euvv3b6XQXk7+E77rgDTz/9NBo1auTN08K8efPQqVOnSktamUwmPPDAA6hTpw5iYmLQs2dP/P777x4fV6PRIC4uDnq93mm7wWDAyJEj8dZbb/n8u64YN26cfckuxYcffggAuO6669w+Tln4/Nlnn0X37t3x4YcfVvk94U+ff/45LrjgAnTr1s2+TafT4aabbsLvv/+OI0eOVPt4Tz4bv/32G3Jzc3HzzTc7Pf7aa69FXFwcPv/8c/u2X3/9FX379kVMTIx9W3x8PHr37o21a9ciNze3Rs+ZZB4HaxaLBatWrULHjh1Rt25dr05SXl6Oyy+/HP3798eXX36JJ598EqWlpejXrx/eeecd3Hvvvfj2229x0003YdasWZXSrN66+uqr0bx5c3z66ad46KGH8P777zt1JQkhMHz4cLz77ruYOnUqPv/8c3Tt2tW+UHFV2rdvj4ULFwIAHn30UXta+5Zbbqny+XrKk+ObTCZcfvnlGDBgAL788kuMGzcOL7/8Mp577rlqj//tt98iIyOj0uK4f//9Nzp16oTffvvNvs7bzJkzUVZWhvLycgByoNa5c2csX74cjz/+OJYtW4bx48dj5syZmDBhQqVzzZkzBytWrMDs2bOxePFiFBUVYciQIfbA+bHHHsM111wDAE5dBI6fry+++ALz5s3D448/juXLl6NXr14A5IBw4sSJ+Oijj/DZZ5/hqquuwl133YX//ve/Hr/Wiv3792Po0KEwGAxYsGABvvvuOzz77LOIjY21P3dPKW08ePAgXnrpJSxbtgyPPvqo00V++/bt6NSpE7Zu3YoXX3wR33zzDYYOHYrJkye7/Kz07dsXy5Ytq/HFKVx069YN69evx+TJk7F+/XqYTKZqHzNy5Eh06NABjz76qEf7V7Rv3z4AqLSwuDL+yVerVq1Cenq6ywDN0d69e3HgwAG0atWqyv1KSkqwZ88etG3bttJ9bdu2RUlJiX2BbMXkyZORnZ2NO++806u2l5eX44cffkC/fv0q3TdhwgS88MILGD16NL788ktcffXVuOqqq3DmzBmXx7JarTCbzTCbzTh+/DieffZZbN26FTfddFOlffv27YsDBw5g69atXrW3ouuuuw5arRYffPCBfdv8+fNxzTXXuO0GLSkpwQcffIBOnTqhdevWGDduHAoLC/Hxxx+73F8IYX9ejjfH32Xlj1tPxrFt3brV7XsLoMpucm8+G8prW3FfvV6PFi1aOL325eXlLmvXKdu2bNlS3dMiT3iagjt27JgAIK677rpK95nNZmEymew3x3T/mDFjBACxYMECp8e88cYbAoD46KOPnLY/99xzAoD4/vvvhRDu15MTonK3lLIe2qxZs5z2mzRpkoiKirK3a9myZQKAeOWVV5z2e/rpp2vcDeru+QpRuQtAUbGrzpPjV3zdhgwZIi644IIq2y2EEC1bthT/+c9/Km3v37+/SEpKEidOnHD72IkTJ4q4uDhx4MABp+0vvPCCACC2bdsmhDj/nrVp08Zpjbnff/9dABAffPCBfVtV3aAARGJiojh9+nSVz8lisQiTySRmzJghUlNTnT5/nnSDfvLJJwKA+Ouvv6rcr+Jnw9W6hDk5OSInJ0eUlJS4Pc6gQYNEgwYNRH5+vtP2O++8U0RFRVV6vv/3f/8nAIgdO3ZU2b5IcerUKdGzZ08BQAAQer1edO/eXcycOVMUFhY67eu4zuIPP/wgAIjXXntNCFF1N+iSJUuEyWQSxcXF4tdffxUXXHCBuPDCC8WZM2ecjt+/f3+h1Wp9eh7K+1rxe6gik8kk+vbtKxISEsTBgwer3PfIkSMCgJg5c2al+95///1KXWjffPON0Ov1YsuWLUII992Srqxfv14AEB9++KHT9h07dggA4p577nHarnShueoGrXjTaDRi2rRpLs+7e/duAUDMmzev2ja64viZGDNmjOjYsaMQQoht27YJAOKnn35y+x38zjvvCADijTfeEEIIUVhYKOLi4kSvXr0qncfV81Ju7777rn2//fv3C61WK8aNG1dt2/V6vZg4cWKl7WvXrhUAxPvvv+/2sd58NpRrYW5ubqV9Bw4cKJo3b27/d7t27UTz5s2dhiCYTCbRpEmTattEnvNL6Y4OHTpAr9fbby+++GKlfa6++mqnf69cuRKxsbH2zIpC+evCm9ktFV1++eVO/27bti1KS0tx4sQJAPJfswBw4403Ou13ww03+HzOiio+X3+SJAnDhg1z2ta2bVuPZmMdPXq00l/xxcXFWL16NUaMGIH09HS3j/3mm2/Qr18/1KtXz+mvRCUjWXGA69ChQ6HVap3aCMCjdir69++P5OTkSttXrlyJSy65BImJidBqtfbBx3l5efb32VPt2rWDwWDArbfeirfffrtS5sFTu3btwp49ezB+/Hi3A7VLS0vx448/4sorr0RMTIzT6zhkyBCUlpY6dRUD57vFquviiBSpqan4+eefsWHDBjz77LO44oorsGvXLjz88MNo06YNTp065fJxAwYMwMCBAzFjxgwUFhZWeY6RI0dCr9cjJiYGPXr0QEFBAb799lskJSU57ffjjz9W6lb0xLJly3DHHXfgmmuuqXKijxAC48ePx88//4x33nkHDRs29Oj4VWX7lPvy8/MxceJEPPjgg1XOanTn6NGjACp327r7fh0xYoTLrmtAnlSxYcMGbNiwAStWrMADDzyAZ599Fvfff3+lff35+zBu3Dhs3LgRW7Zswfz585GTk4PevXu73X/+/PmIjo62d5PGxcXh2muvxc8//+xy8smIESPsz8vxNmTIEPs+WVlZMJvN9u7V6njy3vrr8e72ddx+1113YdeuXbjzzjtx5MgRHDp0CLfddpv9e97V0ATynsevYlpaGqKjo11eaN9//31s2LABX331lcvHxsTEVEor5+XloU6dOpU+DBkZGdDpdDWa7ZOamur0byUdW1JSYj+3TqertF+dOnV8PqcjV8/Xn2JiYioFA0ajEaWlpdU+tqSkpNJjz5w5A4vFggYNGlT52OPHj+Prr792Csz1er29a6biRbK698ETrrrcf//9dwwcOBAA8H//93/49ddfsWHDBkybNs3r4wNATk4OfvjhB2RkZOCOO+5ATk4OcnJy8Morr3h1nJMnTwJAla9jXl4ezGYzXnvttUqvo/IFXvF1VN4vb59XuOvYsSMefPBBfPzxxzh69Cjuuece7N+/H7NmzXL7mOeeew6nTp1yKtfhbr8NGzZg9erVmDZtGo4fP47hw4d7NC60OsuXL8dVV12FSy+9FIsXL3Z7QRRC4JZbbsF7772HRYsW4Yorrqj22MnJyZAkyeX3pzIuLiUlBQAwbdo06PV63HnnnTh79izOnj1rL8tQXFyMs2fPVtn1rnweK36fKOeu+H3q6jtX0aRJE3Ts2BEdO3bEJZdcgpkzZ+KWW27Biy++iJ07dzrt68/fh969e6NZs2Z488038e6772LcuHFu349///0Xa9aswdChQyGEsL9mSsJBmSHqKD093f68HG/Ke+Ct1NRUj95bV7z5bCjvk7t9Hc8zbtw4PPvss3j33XfRoEEDNGrUCNu3b7evv1q/fn1Pnx5VwfWfOS5otVr0798f33//PXJzc50uosr4J2VgeEWuPvypqalYv349hBBO9584cQJmsxlpaWkAzv9iVvySrGkwZzabkZeX5/TlcezYMZ+P6cjdL3tUVJTLL/tTp07Zn2+gpaWlOQ1mBuRfUK1Wi8OHD1f72LZt2+Lpp592eb/jxBJ/cfVafvjhh9Dr9fjmm2+cLhRffPGFz+fp1asXevXqBYvFgo0bN+K1117DlClTkJmZWeVgY0dKVrKq1zE5ORlarRajRo3CHXfc4XKf7Oxsp38r71dtfUZCkV6vx/Tp0/Hyyy9XOZapXbt2uP766/HSSy85ZTcqUoIHQL6gR0dH49FHH8Vrr71Wo0XAly9fjuHDh6NPnz749NNP3U7UUgK1hQsXYv78+S7HbrkSHR2Npk2buhwntGXLFkRHR6NJkyYA5HFJ+/fvd/lH6pgxYwDIf8hVzCYqlM9jxe8T5Tv12LFjThdq5TvXU23btoUQAps3b0aLFi3s2/39+3DzzTfj0UcfhSRJ9uftyoIFCyCEwCeffOKyDt3bb7+Np556yqk3wd/atGnj9r0FUGWG1JvPRps2bezbHcc3m81m+0Q0Rw8++CCmTJmC3bt3Iz4+HllZWZg4cSJiY2PRoUMH758oVeJVfvLhhx+GxWLBbbfd5tMgXUcDBgzAuXPnKl1g33nnHfv9AJCZmYmoqChs3rzZab8vv/zS53MrA2IXL17stP3999/36PG+ZIgAeTZoxeexa9euSvVtfD2+J1q0aFGpDl50dDT69OmDjz/+2G0XEgBcdtll2Lp1K3Jyclz+tehLsObLc5UkCTqdzulLUanxU1NarRZdunSxz8z0piBt8+bNkZOTgwULFrjNwMTExKBfv37YtGkT2rZt6/J1rJh92Lt3LzQaTaUZd5HK3eyyHTt2AKj+j4annnoK5eXlXk38eeCBB9C0aVM8++yz1XahuvP9999j+PDh6NmzJ7744gu3C4oLITBhwgQsXLgQb775ZqUZedW58sorsXLlSqeZo4WFhfjss89w+eWX27siZ8+ejVWrVjndXn75ZQDyzPpVq1ZVWYC7ZcuWAFDp+0SZ2V7x+/Wjjz7yqstYqTNZsZtVGaZQcZKUr8aMGYNhw4bh/vvvd5sFslgsePvtt5GTk1PpNVu1ahWmTp2K3NxcLFu2zC9tcufKK6/Ezp07nWq7mc1mvPfee+jSpUu1n31PPxtdunRB3bp1K1Vh+OSTT3Du3DmXkwCNRiNat26NrKwsHDx4EEuWLMGECRMQHR1dg2dMCo8zawDQo0cPzJkzB3fddRfat2+PW2+9Fa1atYJGo0Fubi4+/fRTAPCoC3D06NGYM2cOxowZg/3796NNmzb45Zdf8Mwzz2DIkCG45JJLAMgX5ptuugkLFixATk4OLrroIvz+++8eB1auDBw4EL1798YDDzyAoqIidOzYEb/++qvHF3tl5YbFixejZcuWiIuLQ7169ar9RRk1ahRuuukmTJo0CVdffTUOHDiAWbNmVRon5uvxPdG3b1/MmDEDxcXFTlOtX3rpJfTs2RNdunTBQw89hKZNm+L48eP46quv8OabbyI+Ph4zZszAihUr0L17d0yePBkXXHABSktLsX//fixduhRvvPFGtV2pFSl/wT333HMYPHgwtFot2rZt6zbjAMhj4V566SXccMMNuPXWW5GXl4cXXnjB7cWvOm+88QZWrlyJoUOHolGjRigtLbV3aSifQ0/NmTMHw4YNQ9euXXHPPfegUaNGOHjwIJYvX26/eL3yyivo2bMnevXqhdtvvx2NGzdGYWEh/v33X3z99deVCnX+9ttvaNeuncuxe5Fo0KBBaNCgAYYNG4YWLVrAarXir7/+wosvvoi4uDjcfffdVT4+Ozsbt99+u1fd3Hq9Hs888wxGjBiBV155BY8++igA+Y/K1atXVxuE/PLLLxg+fDjq1KmDRx55pFLB6wsvvND+vTl58mTMnz8f48aNQ5s2bZzGMBqNRlx88cX2f7s6/3333Yd3330XQ4cOxYwZM2A0GvHss8+itLTUqQRQu3bt3La3VatW1a780aBBAzRp0gS//fYbJk+ebN/esmVL3HTTTZg9ezb0ej0uueQSbN26FS+88ILba8Pu3bvtzzM/Px8//PAD5s+fj44dO9pngCt+++03aLXaSmPLJElCnz59vF5RpF69etVm5ZctW4ajR4/iueeec/m6tG7dGq+//jrmz5+Pyy67zL79+PHjlcagAvI1Ugk2Dxw4gJycHIwZM6bacWvjxo3DnDlzcO211+LZZ59FRkYG5s6di3/++Qc//PCD0741+WxotVrMmjULo0aNwsSJE3H99ddj9+7deOCBB3DppZc6rQiydetWfPrpp+jYsSOMRiP+/vtvPPvss2jWrJlPs/PJDV9mJfz111/i5ptvFtnZ2cJoNIqoqCjRtGlTMXr0aPHjjz867VtVEde8vDxx2223ibp16wqdTieysrLEww8/7FQgUgi5aN8tt9wiMjMzRWxsrBg2bJjYv3+/29mgJ0+edHr8woULBQCxb98++7azZ8+KcePGiaSkJBETEyMuvfRSsXPnTo8Ln37wwQeiRYsWQq/XOz2mqudrtVrFrFmzRJMmTURUVJTo2LGjWLlypcsZi94eX3nu1fn333+FJEmVZpMKIcT27dvFtddeK1JTU4XBYBCNGjUSY8eOdXo/Tp48KSZPniyys7OFXq8XKSkpokOHDmLatGni3LlzQgjXs+wUFV/fsrIyccstt4j09HQhSZLT+4QqikouWLBAXHDBBcJoNIomTZqImTNnivnz51d6nz2ZDbpu3Tpx5ZVXiqysLGE0GkVqaqro06eP+Oqrr6psu6vZoMrxBg8eLBITE4XRaBQ5OTmVZsbt27dPjBs3TtSvX1/o9XqRnp4uunfvLp566imn/QoLC0VMTIx48cUXq3wOkWTJkiXihhtuEM2aNRNxcXFCr9eLRo0aiVGjRont27c77es488/RyZMnRUJCgsdFcRVdunQRycnJ4uzZs/bje/J7p/x+urs5foaysrLc7peVlVXp+bk6/7///iuGDx8uEhISRExMjBgwYID4448/qm2nN7NBhRDiscceE8nJyZW+s8vKysTUqVNFRkaGiIqKEl27dhXr1q1zWxTX8RYbGysuvPBCMX369EozpoUQolevXmLYsGFO2woLC91WK6jI3WfCUcXZoMOHDxcGg6HK2fLXXXed0Ol04tixY0KIqmeD9ujRw/445fvSVaUAV44dOyZGjx4tUlJS7K/tihUrXD7Pmn423n//fdG2bVthMBhEnTp1xOTJkyvNuP7nn39E7969RUpKijAYDKJp06bi0UcftV8PyD8kIVi8KdIMGzYMZrM54Cl7qrn58+fj7rvvxqFDh5hZI9U5evQosrOz8c4772DkyJEBP9+ePXvQrFkzLF++HJdeeql9+9KlS3HZZZfh77//tmfricIJg7UItHXrVlx88cVYu3atfbUEUh+z2YwLL7wQY8aMsc90JVKbBx98EMuWLcNff/0V8DINN998Mw4fPowVK1Y4bb///vtx5MiRGg2PIVIzr8asUXho3bo1Fi5c6LfZrxQYhw4dwk033YSpU6cGuylEbj366KOIiYnBkSNHPK4D5wuz2YycnBw8/PDDle57/vnnA3ZeIjVgZo2IiIhIxVhamIiIiEjFGKwRERERqRiDNSIiIgpLa9aswbBhw1CvXj1IkuTRSjerV69Ghw4dEBUVhSZNmuCNN94IfEOr4fEEg7KyMqeq7FarFadPn0ZqaqpHi8cSkeeEECgsLES9evW4ELIH+P1EVHsC8f1UWlqK8vJyj/Y1GAyV1qR1p6ioCBdddBFuvvlmXH311dXuv2/fPgwZMgQTJkzAe++9h19//RWTJk1Cenq6R48PGE8LslVX1JE33njz/+3QoUOBqK8Ydvj9xBtvtX/z1/dTSUmJgC7G4/PWqVNHlJSUeH0eAOLzzz+vcp8HHnhAtGjRwmnbxIkTRdeuXb0+nz95PBu04l+u+fn5aNSoEf7ddwjxHiwvRUSeKywoQNPshjh79iwSExOD3RzVc/f9tGvvQcTH+/795NGXY4QJlTwl3zv3avoeFhYWoHmTRn77fiooKEBiYiKMrW4GtO6XGgQAWMpRtm0hDh065LR8mdForHbJQUmS8Pnnn2P48OFu9+nduzcuvvhip+XoPv/8c4wYMQLFxcXQ6/UePSd/87gb1N0LEZ+Q4NFaoETkPXbhecbt91N8zb6feMGvLNQ+kXwPK/PXe+j37yetAVI1wZryflas6Td9+nSn9U19dezYMWRmZjpty8zMhNlsxqlTp1C3bt0an8MXLIpLREREwScBqC4AtN3tKrPmt2ZUaIPSARnMP54ZrBEREUUQJTsVallSRwkB6tWrU6dOpdV9Tpw4AZ1Oh9TUVL+fz1OcZkZERETBJ2k8uwVQt27dKq09+/3336Njx45BG68GMFgjIiIiNZAkz25eOHfuHP766y/89ddfAOTSHH/99RcOHjwIAHj44YcxevRo+/633XYbDhw4gHvvvRc7duzAggULMH/+fNx3331+e5q+YDcoEVEFHJReWSh3mVHk2rhxI/r162f/97333gsAGDNmDBYtWoTc3Fx74AYA2dnZWLp0Ke655x7MmTMH9erVw6uvvhrcGmtgsEZERERq4Ek3p5fdoH379kVVFcoWLVpUaVufPn3w559/enWeQGM3KBEREZGKMbNGREREwefJmLQIrT3JYI2IyIZj1SiSCKhtLKInsz0js0MwMp81ERERUYhgZo2IiIiCj92gbjGzRkREFKEE2P0fCphZIyIiouALQOmOcMFgjYgIzC4QBR27Qd2KzBCViIiIKEQws0ZEEY0ZNSKVYDeoW5H5rImIiMiOEw3UjZk1IiIiCj6OWXOLwRoRRSRmEYhUht2gbkXmsyYiIiIKEcysERFRtZRMZKh0QjFz6pugvs+S5EFmLVQ+gf7FYI2IIk44XMiFqPwspAi9kBGFOwZrREREFHwaSb5Vt08EYrBGRGFLVPhJkYHvd80JhE6XdyRgsEZERETBx9mgbjFYI6KwFm5ZFldj1RzvC/S4tVCbaEAhhHXW3IrMEJWIIkJkfq3XDjUGwcGqwi+EqDKIDlVc1UA9mFkjIiKi4GM3qFsM1ogobIVTVkCNmRu1dIkG85VR4/vibxWfYfg/Y/VhsEZERETBxzFrbjFYIyKiGglWho0ZnjDDblC3GKwREYUAZZZndd1uyv3BWM3AsWWBOrvaAjSuGkG1gcEaERERBR+7Qd1isEZERH6ntgwYUShjsEZERETBxzFrbjFYIyIKIdWNXeMYKgpZ7AZ1KzJDVCIiIqIQwcwaERERqYAH3aARmmNisEZEFIIkSQpqmQ4iqj0M1oiIiCj4OGbNLQZrREQhSAjBjBqFF0nyYDZoZH7mI7Pzl4iIiChEMLNGRBRCHEt2VCzfwUwbhTTWWXOLwRoRUQA5hk/eVvV3nEBQ3ZqgpH4msxUAoNOe/1QwwCZPMFgjIiKi4OMEA7cYrBER+VFVlxJPs2xKFu1EQRkAQKfVICVWb7vPdiyHg5WUWwAAUXotNJrIvJiplZJNKyg1w2KV37wYgxYAEGvU2rdp+b6xG7QKkfmsiYiIiEJERGXWkjvdWeNjnNnwuh9aQt7ie0dq5ktORHnMkTMlAIC2I56XN5SXAELOxsBqOf9TyShobV/bxlggKRMAUK9pIwDAwyNaoWejNHk3W6amXnK0/ZwsoltzVlsmzGp7LfccL0KPW+bJd57Lk3+Wl1R+oEYLGGzvhTFW/hmbjAYtmwAA5tzcEQBwccMkxBjl9zjism7sBnUrrIM1f1zgPTkmgwD/CsT7VvG4fM/IH2py2UjtcpdvD7SY5J/lJcC50wCAo8f2AADu2rwV+rQ6AIC6jTIAAE9d3xYX1UkCACTbulKVYIA8d6aoHADQ7OoX5Q2lhfLP8pLzQbXC1WQQixkosT1G+Xn2GA4f2QEAuOLvbQCA+m1b4Y3xnQAATdPjAADpCUY/PQsKVfyNJSIiouDjmDW3wjZYC1R2prpzMWPjm9p8vxzPx/eLgsHnrFpFSnepqVT+eeYoTPnHAQAHDxgAAGP2HEWdRnK27YaBTQEA9/XJAQDotZF54fPUjiMFAIA+D3wG6+F/5I1KF6fy2vvLyf0AgCM/7sd4kxkA8NSodgCAK9vUBxAB3dfsBnWLv6lEREREKhZ2mbXaztBUd35mblwL9vukYIaNfOXL3/d+y6i5Y7WcHz9llsdYiQNbkXt4JwDgxY3rAQBfr+sNALj9Pznoly2PbWuQEo1IZrbImbJ/cs9h9Fu/AQAOfP+NbweTJNfj1jx0Ys1yAMD9JXLplk5PDgYANEyNqfJzF+plkyVJqj57GKGZtbAJ1tRy8a+IXaTnqfU9Ahi0UWAFPEirisXkPCkBwK5VawAA9+04hNbtswEAM4e3BgC0rp8QkAkIjiswqKU7z2oVKCiVuxw7PSQHZmc2rTvfrRxk+Rt+AgAMmilPMJh7a2f0a57hdn9Xr2qoB3AkC5tgjYiIiEIXM2vuhUWwpuaMjaNIzd6EyvsDyG2NtPeHvOPtpSL1kif8dGI/XqSK8wEAlsP/4O/CswCAR2yHv+OSHLTNSAQANE6X64H5u86Xq3VOazPbVmaSu4r/PV6E4bN+BACc2bhGaZx/TqI8Hz8c78TPcrfoI8kx+ObungCA1DiDZ81w+G9m2UJXWARrREREFOIkVP/XUGQm1hisBUMkjGMLpWwaUSBst5V9QNHZoLbDJSXbU14CnD4CAPhreTEA4M7Nh9CybUMAwH8vbwUAaNcwCdG29Sy9P5Vn+RxP9/M1AyeEQFGZnFFrftuHAIDyfzb6dKxg2bV5L1btbQ4AuLK1XM4jnFY3YDeoeyEdrIVDQBBuXaN8T4hkvca9Jv9Hxer2aqDUCBM4H7gVnAAAlBWdwV+H5RUR7i2SJyY8dtWFaJMpd41mJkYBAAy6qis/eRp8ecvbJbOUWZ77Txaj68T58sYT+wLStoDb/zdmfdEYANCviTzRwNPuUAptIR2sERERUXhgZs09BmsqEcpdo+GQTSPyF3v359ljwW2IpxyzbICcabOtObp75U8AgIn7T6BF6wYAgGeukLtGO2WnuD5cgDJqVZ2n4gXescuzxR0fAQBKj+cCeYdqpW1O/DjRAAD2rPkFALD5urYAgL7NMyJ1GFdEYbBGREREQcfMmnshGayFeyYnVMZM8X2gSOSYH3F12eg15YPaakpgCOv5xbJLCuUf23/Hpj3bAACTz8krIzw7oi1y0uIABH/1AyXLVmqSs4RnisrR9vrZ8p22tVLDRmEeAOCxz+T349u7k5EQrQcQ+hMlGay5F5LBWqRQY9douAdoRN5QAjfl8nG22AQc3Bqs5gSOEEDpOQDA7qVypf8ROw+jc095ZuKsYXLXaHZ6LIz62l9y+li+vCzT4Od/AgAcXvldrbehtu1YswEAsPu6dujQODnIraFAY7BGREREwcc6a26FXLAWqZmdYGXZIvX1dsRVDag6SoYtZ/hzgTmBmrp+bF2Olj2bsG7vXwCAvr+0AwC0bp+NaUNbAACauegitQrlEPJ/aCTJ56dWUGJGn6fl1QcO/7Ja3mhb+1SVari4eyW2yRJj3/gNG/47EAAQpZdr4XGN0PATcsEaERERhR+OWXOPwVoIqpjt8nfWh9k0Ih+VFwe7BbVLybLZMmx/H9iGGzbtBwC0bt8YAPDclW0AAI1SYpAQLV9ylCyPRgL02oplNwCzLQWnVOeXAJTbits2GPW2vOPBLf7NVNkbYPX/MQMod38uSsvlMiVKZi1Us2iS5EGx48iM1VD7I0GJiIiIyGPMrIUBf4xnYzaNyA8s5mC3IDiUDJe5HOZ9mwEAfx3bDwC4akcuAKDVRQ0xy7bWqJI9SYzRo35ylNMhjp4txdkiuTxInST5vhZjFwL5tiLDxfnO5wwlfi6QCwBRsVEoLJU/d0mxob30lAQPukEjNLXGYC3MMOgiqn3H80vl/1DzAHdvKF2Bkg+dL0rAalsFoWjLOgDA71vWoe97/mhcgIRY96ei9OQJ7D4ll1VRJnP4utg9qReDNSIiIgo6TjBwL2SCNWaMKJi4mgG5YrUNhO849fMgt4RqTMkihlqG7fge/PfrHQCA7nenAjg/0SDksM6aW5xgQERERKRiIZNZIyJSG6XERMnh/cFtCNVcqGXUHGz9cy8AwGzpIW/QB7ExNeFBN6hgNygREXnDoLN1TlhMwW1IoNRkokGoCOEgTSEscp01syX0nwu5xmCNiIiIgs6TCQaROtOVwRoRUU2Fc+YpHIVBNs1RbFI8AKDUFNrPi8Gae/yGISIiIlIxBmtERDVVUiDfwlW4ZKKENbjPRQj/r7yg0aLobCGKzhYiv9iE/OIQHj8peXjzwdy5c5GdnY2oqCh06NABP//8c5X7L168GBdddBFiYmJQt25d3HzzzcjLy/Pt5H7AYI2IyEdHTpfgyOkSwGqRb0S1zWqBpNVC0mqx7VQ+tp3KD3aLVGfJkiWYMmUKpk2bhk2bNqFXr14YPHgwDh486HL/X375BaNHj8b48eOxbds2fPzxx9iwYQNuueWWWm75eR4Ha2VlZSgoKHC6ERGpAb+fiEKfMmatupu3XnrpJYwfPx633HILWrZsidmzZ6Nhw4aYN2+ey/1/++03NG7cGJMnT0Z2djZ69uyJiRMnYuPGjTV9ij7zOFibOXMmEhMT7beGDRsGsl12yZ3u5OoFpBr8PKpTsL6fLntpDS57aU1kZNaC3YVYU6Hc9mqI3D0QuXvw1Mfb8NTH2+wraxBQXl6OP/74AwMHDnTaPnDgQKxdu9blY7p3747Dhw9j6dKlEELg+PHj+OSTTzB06NDaaLJLHgdrDz/8MPLz8+23Q4cOBbJdREQe4/cTUejzJrNWMZNeVlbm8pinTp2CxWJBZmam0/bMzEwcO3bM5WO6d++OxYsXY+TIkTAYDKhTpw6SkpLw2muv+fcJe8HjYM1oNCIhIcHpRkSkBsH4fhJC4NCeIzi050jAz0U1oJaMoL8nFjgqPQeUnsPBP/7GwT/+hlUICNvNn0SAjqvwJlhr2LChUzZ95syZ1R674nNx16W6fft2TJ48GY8//jj++OMPfPfdd9i3bx9uu+02/zxRH7DOGhGRF5QL1ZkiE1BWWlsnBdRSXyoSVjUINbb3RJdWR/6nvw4byACzhg4dOuT0R5nRaHS5X1paGrRabaUs2okTJypl2xQzZ85Ejx49cP/99wMA2rZti9jYWPTq1QtPPfUU6tat66dn4Tn+thEREVHQeZNZq5hJdxesGQwGdOjQAStWrHDavmLFCnTv3t3lY4qLi6HROIdHWq0WQPACWGbWiChsVeyy8Wf1c4NOA5QV+e141VKeh1oybGqnhq5PILDdnwqdAQCQkOy/7n81Z9W8de+992LUqFHo2LEjunXrhrfeegsHDx60d2s+/PDDOHLkCN555x0AwLBhwzBhwgTMmzcPgwYNQm5uLqZMmYLOnTujXr16QXkODNaIiIgo+DwpeuvD3yojR45EXl4eZsyYgdzcXLRu3RpLly5FVlYWACA3N9ep5trYsWNRWFiI119/HVOnTkVSUhL69++P5557zvuT+wmDNSKKGK6yBd5m25T9JQnywG5APVmc2iSsHLdWldrMTFnMAIDS4spjKJXPvKef82Bm1AK5NuikSZMwadIkl/ctWrSo0ra77roLd911l0/nCgQGa0QU0Xy9mBl0GsAShKV92B1avWAHz7Ud8Nieb3HeafsmpdSaxsOPSTh1e4YjBmtEREQUdIHMrIU6BmtERKi67pJyvyOdRjqfQQlGVkIN5TxYxkMdlNe/iowiM2ehjcEaERERBR0za+4xWCMisqk4fq3abESwLxwcv+Ys0saqKWzvf1RKikNT5LYISI67qFuAZoOGA+auiYiIiFSMmTUiogo8Ht+j1Qe2IZ4KdoYt0seuBXs8mK10h05f+ZKutCwUElLsBnWPwRoRkQ+EQIj0LVGkOHfilP2/7eGjQ3coP66hi8EaERERBR0za+4xWCMiqobbXi6Nyr5Cg13OI1irGgRzYkGwu0CrobTOKgS0Kg90JHgQrIVEh67/RegAAyIiIqLQoLI/C4mIQohWhV+hwZ5sUJuCXapDZeIz093eJ+D9ElS1jd2g7qnwm4aISB2q7eFSeRdYUETKzFAVvvdFBUXyfwi4nP4pIimQDzMM1oiIiCj4WBTXLQZrRES+UmM3qCLcsyjsAq3EVZ01R+rLBTpjN6h7YZ6nJiIiIgptKv6zkIhI5cJ9XFZNBGrsGjNqbknVzBxQe06KmTX3GKwREflKzYGDWi5q/gja1Pw6q0h5aTkAt/MLKIQxWCMiIqKgk6Tq/8ZQy98gtY3BGhGRrzTaYLegMrVezRyzY55m2ZhR84ztPY+KiXK/S221pQbkYK26btBaaozKcMAFERERkYoxs0ZE5AMBqG+CQaikHcIhY6a81moojmv7HJrKTVXvZmuzaqu6eNANGhIpwgBgsEZE5AMJcBg8H+QLt+quulRrJMle789cVFz5bjcPodDCYI2IiIiCjqU73GOwRkTkRlUJM1V0g0bohUtVJCm4XaGS5vzn8MRe97tJUqT2IIYFBmtEREQUdCzd4R6DNSIiLzjmULQJyQAAy0lbZkNYaq8hkXrVUqNgj1lURCcE/BSedFX6SqORoKlmFQZRzf3hisEaEVE1XPV0mcxWNGpaDwCw78A2eWN5SS23jFSlNoM2x4DJVu+vYedO8l04P7FACawk1Dy+j9TxYmrAYI2IiIiCjt2g7qk+WDuz4XUAQHKnO4PcEqLzn0cKDTXtshEOGRL7YWyb9DoN/juyDQDgpj+2yBtP7PP5XB6L1KtVKKnlSQfRzdoCAKZde6Ht/M4ZNaVJVWHWTN1UH6wRERFR+GPpDvcYrBERuXG+4rvrDFuzlDgAgC45DQBgDmRmLUIvUlSNqDjEJ8cDAHIS5Z8ahzIdVX1s1Bb4sBvUPZWtlUJEREREjphZIyKqhmMGQsmySRJgtSXcomOjAQCFgRqrFKnphFAW6JmhSiFcnQFarTwbVKlqUd3MT7Vl1BTsBnWPwRoRkRccLxZmi7w2qN6gD1ZzSO0CFbTZ1gN1LBdjrWKB9kgNcsIFgzUiIiIKOmbW3GOwRkTkI41SHsHe/6Tx7yoGEXphCkv+6iKXHD5rAGC1wFRucr4rRD83nGDgHicYEBEREakYM2tERD4y2wYJ2Ut7+OvP/khNH4S7QIxfi0lESZE8bi3UPzUSPOgGDfln6RsGa0REPrLYgjWNxtZJodUDFnMQW+QDqYoOFmGtvXYES1XPX6HG10EJaqLiYDHLXe/aCF3kPBIwWCMiIqKg45g19xisERH5SOn+tJfu0BmcSin4JBBXI0+yR9U9To3ZJV95+3r4+3WoSXeo8hjlpzEGBqMBQOhn1jgb1D1OMCAiIiJSMWbWiIh8ZNDJf++mpsYAAHKNsUBxfjCbVJmvWTV3x/F3hq1idimQmRN/vBaBeh28ZSoFAGj0eiSmyGuCRum1wWxRjbEb1L2QCdbObHgdAJDc6c4gt4QikfL5I3KUkRgFALisU30AwNYNacCZo74dzN9XIX8Fae6O60uwojzWIC/PBZ0BsTkX2u6Sn/+5vbsAi1w3TAlIaty17Hhufwl20GYLcqNjozGwWxYAICE6ZC7p5CW+s0RERBR0HLPmHoO1CFAxK8TsJJF/xBrlbqeejZIBAProKJgCvYB3dQKVUXN1Hm+ySpIGSJEzkHEN5UxQdtNMTOifDQCwQn693vu1PnbvPAYAyD95Wn7swS01a6da1WRVA70RAGC1WtGpQSyAcJhgwG5Qd1T8KSYiIiIiZtbCjCdjq1ztw2wbkfeMtgkGDZLkCQbxyfE4bR/L5Mc1QtWqYtbKVaZNJ5eV0DVpi/89eAkAoF1dORMZY9AixpadtCXWMLhZHRSVya/dH7lnAACPL26A3I3r5R3OnfatbYHibYaxpjS218soZ9My6qWibqw8BlCndf+cJdhfYtViN6h7DNbCSE0GwZ/Z8DoDNiIfFRTLA+KNUcbg9tMEu8vP8fy21yHqgvYAgPn39UefpukAAJ2tu06n1VTquosxapFiWxkiI0Hu6tPeBNxfVg4AyNvws7xjqK0U4S8VPl9CCPuyZ666QSU3/w2oP3ij8xisERERUfB5MGYtQpcGZbAWDvxVVoLlUYi8o3TJFJTLmbWoKJ28PigQ8ZkfXc7FAIDVz1wOAGiYGoNoQ/V1wCRJgk57PvMGAAOb10HC7d0AALfaBuSf3vhL5L7GgL071Gy2wmSVu2GV143CD4M1IiIiCjqOWXOPwVoIC1ShVmbYiNxzdak4XSqPp4qO1p8v+OqPQq6hRpIQ27orAOD3568AANRNiqrxYWOMWvS2jXf7bGo/AMDN82Oxb/m39vMCCF65lGCISQQAWCxWlNsya1rb61BdOKPWV4mlO9xjsBaCaquaPicdEMmquz7E6+Wuz7oZsdgenyZv9HTWYjiw1fwyNm2L5U8MAXB+coA/MiGSJEGvk4/TNDMOAPD8DRfh5uNnAQCFu7bKOwZrqa/aXM3AFpDqEuUZtampsUi2LeQeqYFMJGCwRkREREHHblD3Qi5Yi/RsT22vUckuUa4LGsk8vSw0TJK7Pm/qXB8bNtQFABTk7gpQq1TEVkMto3MvAMDrEzojJ0Ou/xWoi6oySaF7kzS8ckd3AMB/P0kCAOxb+YO8k8UMWGuxzl1t1lmzveb1G9cBAIzo2QgNbJ+/SA1kIkHIBWtEREQUfjhmzT0GayEi2NkdZtiI3FOyPS1TE5CUGg8AKAj3Qe86A1C/BQBg5li5TEf7Bkkw6gNbmFdjK/wapdegX04GAMB8lfwa37rzkLzT8b1A6bmAtsNvvP2c2Ep2pKfLGcy2GfGI8aAkSihgN6h7XBuUiIiISMWYWSMi8pGSC1H+2pckhwRJMNYIFdbALzmlrE2ZUh/1m9QDALS3rfWZEK2v3CTbC1LTjIhwkXlKjJHP16VBKgDY23OkvBQ4vkfeqTbHrtUG2/Mpti1xFqXVnv/8eXgIZT+15XyZWXOPwZrKBbv7syJ2h1Ik8Paip6x1abI4XP6UoCbcggVbHblGbVvg7QlyTbUGKfI2V2tTKhyDLW8vuK4CNUcNU+XzPzayFQDgha/0+HdVrnxnSaFX51I1SQIS5Hpz586V2Tdz5YLwx2CNiIiIgo4TDNxjsKZiasuqOYr0EipEjpQB3mdKyqHT2bohbYViYS737CBK9kitVyNlzc8GzQEAL426GFnpMQAAjZ8zZb66pJlcziL+Gj1u3H5A3rh/c0DO5TfevBaSBoiTu5xLS+Ru0FKLBQZteAw/Zzeoe+HxDhMRERGFKWbWiIhqSClXcbKkFBaLrUCqrXhprVMKtPp7ooEtA9S2YzYAICc1zl6yRC3JjsRo+ZLWNSsVTS/MAgD8e2i7fKfFHJiT1mZBXEmCpJVf8/w8eWmtgnJTlWMFqzwc1DXJgN2g7jFYUyE1d3864mQDIpnSDZgeHYWGDeUFtg/+nSTfWZjn3cGEUOcVqXFbAMCtfRsDAFLjDNCrrPtNCVqSY/W4c2hTAMC9O1sDAKz7/g79mndWC7S2YK1+Y7nGXN3YaFV+XMi/GKwRERFR0HHMmnshGawxo0O1IVQynOR/9lJpXj6uQVI0RnWpDwDYsDYNAFB+Ym/oZ3Q0WrTrLGeq+jaRMzrRBq3Xr09tGmybbLCk1wUAgHVH//Xvqga12f1pP6dAbIK8csHkYfJEj9RYg9cTPNRKggfdoLXSEvVRVw6biIiIiJyEZGYtXIVqJoeZTgpX3mbYog1aXFxHLq2QWkeuqp+705cT+6GMhx8nGhgvaI/nr5THrMVFyZcNucuqxocOCCGAlDh5gsf0QfL6pcM37EPp9t+D2ayai4pD4xy5KG7DOLlsikGnUe374C2NJFWbJQyXLKK3GKwREfmJXqtBSbm8YoFeWdBco/V9JqK/grYaBmzNL2yAzES5bpwyiF/yU/P8ybG3WRnb1DRD7jZs0boB/vJHsObP7k9Pu8eVF1irQ2mp/FkStj8ljDp2kEUCBmtEREQUdCzd4R5DcpUI1S5QR+HwHIhqIsagRbHJgmKTBfHxRsTHGwFt5cXNQ0aj1kCj1nj9uouRHGtAcqwBOo1kXwtVIURw5lAo53V1fsl2i9JrEaXX4pkrWiPqws6IurCzjyez+i+r5u0LptHKt6S60Ok00Ok0kGz/izZoPZpF6bIZXj8isJTnUd0tEjFYIyIiorA2d+5cZGdnIyoqCh06dMDPP/9c5f5lZWWYNm0asrKyYDQakZOTgwULFlR7nj///BNbtmyx//vLL7/E8OHD8cgjj6C83MOl51zwuBu0rKwMZWVl9n8XFBT4fFIiIn9Sy/eTTiuhoFxes9FiseUttH4YbVLTQrk+Tjbo2rcVACAjwWjPplWV2ajNer7eZvIaJEejc9ccAMCanX/IG60WP7cqwLRaFBbKn/NSi9x2nQ+FidWWUVNoJPlW3T7eWrJkCaZMmYK5c+eiR48eePPNNzF48GBs374djRo1cvmYESNG4Pjx45g/fz6aNm2KEydOwGyufuzpxIkT8dBDD6FNmzbYu3cvrrvuOlx55ZX4+OOPUVxcjNmzZ3v/BOBFsDZz5kw8+eSTPp2E3Au3rkPODKVgCNT3k7fXBY0koW5cNAAgNVX+CWMsUFJY88bU5gzR9MYAgBevkKv/x0froNOen1hQ5Skqdkf6OXjzNkhTlgJLjTPgoQHNAABrfmgp33l4h2cBmz+7P32hkVctMMbFQa+X/7tOTBQAecJHZHYMeu6ll17C+PHjccsttwAAZs+ejeXLl2PevHmYOXNmpf2/++47rF69Gnv37kVKSgoAoHHjxh6da9euXWjXrh0A4OOPP0bv3r3x/vvv49dff8V1113nc7DmcUj+8MMPIz8/3347dOiQTyckIvI3fj8RhQGp+nFrSmRaUFDgdHPMrDsqLy/HH3/8gYEDBzptHzhwINauXevyMV999RU6duyIWbNmoX79+mjevDnuu+8+lJSUVPsUhBCwWuXg/ocffsCQIUMAAA0bNsSpU6c8fSUq8TizZjQaYTQafT5RIDCLQ4EQbtnOSKCm76f0BLkdN9pWMvh1RRpw9lgwm1RZVRk2Ywy6Du4CQM6oAXJ5CF/rW/klIViDfjul3Ua9FtlpchmPdl3l6v9/fXccKDhZxYmDUKbDFbM81ikmPgaDOjcEAKTGy5+zSM2qNWzY0Onf06dPxxNPPFFpv1OnTsFisSAzM9Npe2ZmJo4dc/17uXfvXvzyyy+IiorC559/jlOnTmHSpEk4ffp0tePWOnbsiKeeegqXXHIJVq9ejXnz5gEA9u3bV6kN3mDpDiIiIgo6b0p3HDp0CAkJCfbt1f2xVnGspRDC7fhLq9UKSZKwePFiJCYmApC7Uq+55hrMmTMH0dHRbs/z8ssv46abbsIXX3yBadOmoWlTeZm2Tz75BN27d6/6yVWBwVqQhHv2hllPCge+rBEaa5THFPXNltfQjE1JQtFBfzbKj5VoHQvm2sZFpVzUGTd2lbOCabbsjT+qxjsXrPV8X3+RIE+UAIBb+2cDAB7+9yTyN52qfNJgrPtZFVv5l8y6SRh2gbyCwfmVJILWKr9TypFUtw8AJCQkOAVr7qSlpUGr1VbKop04ccJtpqtu3bqoX7++PVADgJYtW0IIgcOHD6NZs2Zuz3fRRRc5zQZVPP/889DpfA+5WLqDiMiPtJIErSShtNyC0nILomOjPUsZBIvSNkM0YIhG4yapuDgzGRdnJkOvlaDX+r/dFeujVVUvzZ+sQsAqBNpnJqN9ZjJatKoL6Azyzd44Pwdq/nhS5nLAXA69XgOjVgujVgutBATgrQk7BoMBHTp0wIoVK5y2r1ixwm2mq0ePHjh69CjOnTtn37Zr1y5oNBo0aNCgyvM1adIEeXl5lbaXlpaiefPmPjwDGYM1IiIiCjqldEd1N2/de++9+N///ocFCxZgx44duOeee3Dw4EHcdtttAOQJSqNHj7bvf8MNNyA1NRU333wztm/fjjVr1uD+++/HuHHjquwCBYD9+/fDYqk8w7isrAyHDx/2vvE27AYNgnDvAnV0ZsPr7AqlkCfgeVeoxnY1UXIp8QnROKV0NQo/1vXyR1EzSXM+q5RUBwAwtlcWGqTIF6RwWzRbeT7ZtvVC7+rXBOu/tHWlmaqYaBBscXL5CL1eY/8cKuOtvKnor9b6aoE2cuRI5OXlYcaMGcjNzUXr1q2xdOlSZGVlAQByc3Nx8OD5sQpxcXFYsWIF7rrrLnTs2BGpqakYMWIEnnrqKbfn+Oqrr+z/vXz5cqcuVIvFgh9//BHZ2dk+PwcGa0RERBR0niwn5etyU5MmTcKkSZNc3rdo0aJK21q0aFGp67Qqw4cPByC3b8yYMU736fV6NG7cGC+++KLHx6uIwRoRkR8J2/ikcrM89kmScH4VA39XzK/pZANJsk8sQFQMAKBteiKiDVo/NE69lNUYWqYnAInyRBAUKhMN/HQSfw6+M8iZzsLCcpit4ZsfC+WF3JXaatnZ2diwYQPS0tL8enwGa0RERER+sG/fvoAcl8EaEZEHPC3joYxZKyqT1xGsXy8Be/Xy0kAwua6yXvPG+Zhh02jtj01OTwYApMQZfBrEHYpSYvVIayCXbzh1ZKe8URI1nxHq7ymtMUkAgPj487NWAzFLN9g0klTtOEk1jqN89dVXceuttyIqKgqvvvpqlftOnjzZp3OERbAWKoPYI2ligaNQqrkWqe8Rea6qoE0A0NsW1s5MlAO0B/o3xc8r5IHMOLA5wI3zcNKB0vWpMwCxcpA2c3xHAHK71XhBDISEaD0evekiAMCUrbb35lye712hAao7ktWqCQDglr6N7Stk6L1cwD18O0+D7+WXX8aNN96IqKgovPzyy273kyQpsoM1IiIiCm2hOmbNseuT3aBERCpSXaYiWi9nPuonRsMYJw/eD1AnqDNPukSVzJpWD8SnAgCaJcUDOD/4PlJ0riuXxUCSrZp90RkAXk4ECWQlXwBRttUKGsbF2D9XagxaaiqQs0FDHYM1IiIiIj+wWCxYtGgRfvzxR5w4ccI+S1SxcuVKn47LYI2IKACUiQbWAGdd3KrqvEpmzRiL9PryWpPKWKhIUzdJHltYv2lDAMCRY/8CFlP1D6zF97WgQM7JaiUJWhdj1cJlPFqodoM6uvvuu7Fo0SIMHToUrVu39lsmkMFaLeCgdVkoTTQgqimtLViLNmih08tftbXSDeoJ24oKmuQMPHdzewBAUow+mC0Kmvho+Xm/MLodAOD69euA8pLKOwYr6AZgsNW9i9JpuR6oyn344Yf46KOPMGTIEL8el8EaERERBV2olu5wZDAY0LRpU78flwu5ExEFgFaSoJUkFJSYodFooNFoPOvnqQ3CCggr9EY9mibFo2lSPGKMOsQYI+/vd8l2a1cvGe3qJQNRcXIWreItiMrKzCgrM6PUbPFoEH6okjy8qdnUqVPxyiuv2Fcy8ZfI+80kIiIiCoBffvkFq1atwrJly9CqVSvo9c7DCz777DOfjstgjYjIj5S/qJUxayazFcZoefB+obKTJAU3W2ObYCCsAnqdrRRE8FqjCsXltnIdCanAqYPBbYzCtqZsoq3AshUirFeXCIfSHUlJSbjyyiv9flwGa0REfuLY9aGzjQSPMWqhN9j+urYN7K/xckY1ZZGXwopPjkdCNC8DABBnlAPYull1kbt3U3Abo8zWtX2eSkvl9ys5yqD6MVuRbuHChQE5Ln9LiYiIKOg0EqrNHIZzZrEqDNYCjGU7KguVtVyJ/CHGqEPzpvIqAblrlCuNBkCF7Fptdosm1QEAPHpTOyRGaMkOhfKqp8bLXdWPXtcKd6xeLm+0ermSQU0oGTPJYd5fTCIAoF49eXWJKIMW2jCu3REO3aDZ2dlVtnHv3r0+HZfBGhEREZEfTJkyxenfJpMJmzZtwnfffYf777/f5+OGTbDGgqtUE8yAUqCkxRlwX78cAMDqRQ5jkZTEmjJ+TflrvBYybLF16wMAujdIhVHHCk6OrmhVH3foDPI/XBXH9SfHDIySUZOk82PWbJm1kZ3rAQDio3TQ2h6j9gyTr0L9ad19990ut8+ZMwcbN270+bj8LSUi8hPHC6hSnqu43ILmGXFonhHnuKN8QdZo5Yu0VLtfxeVl5SgvK0dclC6s63b5IkqvBQzR8q02KO+/UoNP0pz/bNg+RE0SY9EkMda+hFltC7VaZ2o0ePBgfPrppz4/Pmwya0RERBS6wmHMmjuffPIJUlJSfH48gzUiogCwOJXxULq4XJTuULq87JusAe8KVdYqjdJrA3qekBVnu6ieOx34c1WcWODYDWouBwAkGeVuWZ1GqrVuwmCEROEwG/Tiiy+ukGEXOHbsGE6ePIm5c+f6fFwGa0RERER+MHz4cKd/azQapKeno2/fvmjRooXPx2WwRkRU2yQXpTvsEwwCf/rEVHnguj6My0B4SnkFHF/2lMZZAIDTx/4N0EkdsmmOGTVlm5JZs/1UVsOojaxSMD8R4dANOn369IAcl8FagHB2YdU4e5dqgzIYujYXdlIuJo5xkMXqogUVu0TtF20BiADV97INmn9hbHsAgF6n4WBxG8fX4a07ugMArlm/Ut4QqG5pV7NBlckFAKLqNnDaLVgTDCj4GKwRERFR0Hky0zRSw1UGa0REAaBkQ3QaCTtOnqtiR6Vsh0O3aKBqrtWR671lJcTIp/Hv0cOCANC3eYb8D3v208+ZTsdSLfbuT4efypqgJ44BAEwW+d9alXcBUuAwWCMiIqKg00hStQvVR+pC9iyKS0RhLxiFPJXB0lEGLeokRKFOQlQtt8CpMecLrkoa6LXyLVIvfNXJLzYhv9jkNH6sVgmrfCvOB4rzYdRpYNRpoNFIYV3E2P4xreamZosWLUJxcbHfj8tgjYiIiMgPHn74YdSpUwfjx4/H2rVr/XZcBmtEFPYEandGaEXxUTrERwVx1Imy9pWNXqeRZ4KqPEsRDBKAhGgdEqJ1gNZ28zclc+Z0YvfLjmk0kpxV839LVEXJGlZ3U7PDhw/jvffew5kzZ9CvXz+0aNECzz33HI4dO1aj4zJYIyIKsKBfZJSgw9aPpNfKN7Vf+ILFPnZKq5dvwWIL6oQQEAFe1YL8Q6vV4vLLL8dnn32GQ4cO4dZbb8XixYvRqFEjXH755fjyyy9htVqrP1AFDNaIiIgo6MJhzJqjjIwM9OjRA926dYNGo8GWLVswduxY5OTk4KeffvLqWAzWiCisqSEfYRUCVneZEaVLrEJXpV8pA+U1GkCjgVYj2aviUxX0UfItWGxdoxarcF1YOQCCOWRAyWhWd1O748eP44UXXkCrVq3Qt29fFBQU4JtvvsG+fftw9OhRXHXVVRgzZoxXx2TpDiIiIiI/GDZsGJYvX47mzZtjwoQJGD16NFJSUuz3R0dHY+rUqXj55Ze9Oi6DNSIKW2rIqgHAmSJT9TtVHHDuL5IEGGMBADq9/JVvtqjllVGnojKz/B/GmMCcwNVEAvuyY9pK28rN8s8QSCrViCfdnGp/DTIyMrB69Wp069bN7T5169bFvn37vDougzUiIrXxZ+V8ISp15ZWZAxQYhom8c+Xyf/i7C7RipOFm9qedzujf81NAmUwm7N27F6mpqVXuJ0kSsrKyvDo2gzUiIiIKOk9mTat5BrNer8fWrVsD0kZOMCCisKWWr3WNJN+CTaPRQKPRsBSEp6qofVaz47ro73OcYKKc11U9tgALxmofCo2HNzUbPXo05s+f7/fjMrNGRERE5Afl5eX43//+hxUrVqBjx46IjY11uv+ll17y6bgM1ogorEkI/kQDrwf0+zObIkmAxXmCg0Gn9vxEcKXGGeT/sHgwMcTfhDj//lsia2xhqHeDAsDWrVvRvn17AMCuXbuc7qtJ2xmsEVHYU74igxG0CSEQH13FV629m81hMoG/uyjLigAApnI5+GCNtarFGm3vV0lBYE5g7+502KZcyB0Dddsi8kpwLYT6Z0NGulWrVgXkuAzWiIiIKOgkD8Z2RmqwymCNiCiAJEmCUa/1YEc/lutwJARgLq+0iTxQXuLf47nMqLmquSacftZGJlQNMZAnE3FCISm8YcMGfPzxxzh48CDKy51/9z777DOfjsmBC0RERER+8OGHH6JHjx7Yvn07Pv/8c5hMJmzfvh0rV65EYmKiz8dlsEZEFGA6jQRdMFMCFjNgMUNYBYRVwGSxwmSxsnxHdSym2p9k4Fiyw3Z+qwCsIvgTZQJNmWBQ3U3NnnnmGbz88sv45ptvYDAY8Morr2DHjh0YMWIEGjVq5PNxGawREQVYUZn5/BJGQaTVaaHVaWGyCJgsgt2hbuSXmJBfYrIHuX7nWGfNg1pq5WYrys0MrkPBnj17MHToUACA0WhEUVERJEnCPffcg7feesvn4zJYIyIioqBTxqxVd1OzlJQUFBYWAgDq16+PrVu3AgDOnj2L4uJin4/LCQZEFPaCmY8wWaw4WVAWxBacp3QhmW31u6xCQKuKoeXqEmOwTQix+nmyh8LVRANXYpMBAPll8iB1i1VAq5Efq/buQF+Ew0LuvXr1wooVK9CmTRuMGDECd999N1auXIkVK1ZgwIABPh+XwRoRERGRH7z++usoLS0FADz88MPQ6/X45ZdfcNVVV+Gxxx7z+bgM1ogobAkEL6umjC8qM1nRMDWmih2tVf/bj8yHdgIAdp2Ru2maZsbZXx+VJyxqVe5Z+WIbsMyaQlgByUVZF30UACAxuwkAIEYnX6rNVgFDoJpi+xnMz4FGkqCpJnVW3f3BZDab8fXXX2PQoEEA5LV4H3jgATzwwAM1PjaDtQBJ7nQnAODMhteD3BJ1Ul4fonBldYgSXdbJqsUgza5EDtKe/1peBmdA00zoufRUJRdf80xgT2B/rzWVaqoBVnvttXNnz8l72QIUi1XYP1da9cYsEUun0+H222/Hjh07/H5s/pYSERFR0Gk8vKlZly5dsGnTJr8fl5k1IqIAcEx8mK0VOmMds2gVMyu1UJ7h2NEzcrsibKFwj9nWUq0VymdBKJ8Yjb22m+WcvDapck/Fj1FAmmP7ycSdbyZNmoSpU6fi8OHD6NChA2JjY53ub9u2rU/HZbBGREREQRcOs0FHjhwJAJg8ebJ9myRJEEJAkiRYLL6Ng2SwRkQUAMpFRZKA0vIqvqDtmZXay3KVlcilRIrLLUiMCd9SED4LRCFcR47j05SOPcdtysQGjTz5oKDcZNtFOBTGDez75ZjEq61PhgYeTDBQec5v3759ATkugzUiIiIiPzhw4AC6d+8Onc45vDKbzVi7di2ysrJ8Oq7ax+oREYU0nUbCmaJynCkqd15ayGpxvtUijUYDjUaD/GITLFYBS20MhiLXlM+E42dBWZPU9u9dZ85h15lzTmPWnLNsAW6ii1sgKN2g1d3UrF+/fjh9+nSl7fn5+ejXr5/Pxw2bzBpLQVBNsNQK+UvFC2hRmQU3L9jgvJPF7KJ0R+0FTEVHDwMAfj18Gllpcg04vZZ/uwPA3wfza+9kTu+50h0OwCyvWIDCUwCAbzefAAAMyM5AQnTYXLbDkjI2raK8vLxKkw28wXediIiIgs6TtT/VujboVVddBUAe+zl27FgYjUb7fRaLBZs3b0b37t19Pj6DtQBL7nQnszUVMAtKkUAp12G2CuSdlEsw2AeuC2utZtIqOXUAAPB/y//F4GaZAIAYo3w5UOm1MOBKTXJXdP8HPg1OAxw/D0q3eJFcYuXQEfnzU1JuCerHJtAkqfoVCtTaDZqYmAhAzqzFx8cjOjrafp/BYEDXrl0xYcIEn4/PYI2IiIioBhYuXAgAaNy4Me67774adXm6wmCNiCgArLbMWl5hGaxWlRWfjYoDAOQeOoXC0gCXqQgRxUp5ldNHg9sQ4HyWzZZGio+XVwQtNVnsn6twXG8qHOqsTZ8+PSDHZbBGRORHyuBipRtUq5FQcq5EvtNeU00dfVmmchOKyuVgzerQ3kiUEK2X/6M2Vy/wkMb2npiFgEUlnx06r3379vjxxx+RnJyMiy++uMqahX/++adP52CwRkREREEXqhMMrrjiCvuEguHDhwfkHAzWiIgCQEmAJETrYbWoK6OmVMaXJAlbT8mD11vUjQcAxBkj87JwtshWLsNc7rqvLRjvne2cJpP8+dFAOt8NGoYk2/+q20dtHLs+2Q1KREREFAI2btyIHTt2QJIktGzZEh06dKjR8Ris1QIWXJWxZAdFEqU4bkm5RX0TDGzKCwrw2tf/AAD6NU4HELmZtQuGPi7/h7s1Witm22ox01ZWJk9+KDabPVptQkLgVhkIpFDtBnV0+PBhXH/99fj111+RlJQEADh79iy6d++ODz74AA0bNvTpuCxZTUQUABYh3yQJ0Ol10OlrOQiqar0eIeRb6TkcPXACRw+cwOH8YhzOLw7Ji3xNnC024WyxSe7+NJd7HoTV4vpH0dE6REfrUGK2wGQRMFnOt1FycUOF/6baM27cOJhMJuzYsQOnT5/G6dOnsWPHDgghMH78eJ+PG5l/QhEREZGqhENm7eeff8batWtxwQUX2LddcMEFeO2119CjRw+fj8tgjYjIB47XDFe5GJNZ7k4rLDHDVG6qlTYB8CzTYz4/mL68SC5V8UeuPNGga3ZqRJXvOKNMLKhJt6bymgeoa/Ts2VIAwJFzpbjAttKC0hta3Vul3B0KGVNJkqose6Hso2aNGjWCyVT5991sNqN+/fo+H5fdoERERER+MGvWLNx1113YuHGjfdzqxo0bcffdd+OFF17w+bhhkVkLlYHrkTrRIFTeH4BruVL1XP1d77hNKYZbUCL/dX3be39A7PGtEKbXPM06mErP/3eJnFF7/dM6AIDhLeugQYq8rmF16zSGg44jn5X/w93EAm8EKMN2YNteAMCSzDi0S5fXoMxMjAIgFzEOl3cpVLtBk5OTnTJ+RUVF6NKlC3Q6OcQym83Q6XQYN26cz3XYwiJYIyIKNE+vESZbTTXly7tI6WYLpJoEVcYYAMC5/HMAgLxz5fZAwKhT4ZXRj0wWK2Aq8/+BlUkc/mJb0D0v7/wEEE9mhTo1CaHRFRqKZs+eHfBzMFgjIiKioAvVtUHHjBkDQM6gLV68GIMGDUKdOnX8eg6Px6yVlZWhoKDA6UZEpAZq+n5SSiYYdBoYdBro9doAnqyGpSMsJsBiBixmlJeWo7y0HDtOF6C43ILicgsEwjsbs/dE0fmSHWpmLgPMZUhKioLZKmC2Clhst3CikSSPbr6YO3cusrOzERUVhQ4dOuDnn3/26HG//vordDod2rVrV+2+Op0Ot99+O8rK/J+t9ThYmzlzJhITE+03Xwu7ERH5G7+fiMidJUuWYMqUKZg2bRo2bdqEXr16YfDgwTh48GCVj8vPz8fo0aMxYMAAj8/VpUsXbNq0qaZNrsTjYO3hhx9Gfn6+/Xbo0CG/NyZShNKA+5qKpOdKwRPo7ydv/pZXslF6rQS9VsK5wtLqHhI8VgtQeg4oPQdxdDfE0d345I9jOJxXjMN5xbWSvamqdm+gnCk24UyxCcNfXO3QEI188xd/PonCPKAwDxaLFeVWC8qtFnuGTXgxNk7thXKVCQbV3bz10ksvYfz48bjlllvQsmVLzJ49Gw0bNsS8efOqfNzEiRNxww03oFu3bh6fa9KkSZg6dSpef/11rFu3Dps3b3a6+crjMWtGo9G+qjwRkZoE6vvJlwubcu2MsnV/lhSV+K9BjvwRCDgGJ7aB9rv+zcO2jvkAgJzMOACAzuCfrlxvm+y4vz/H6+8/KdeWO7HvsIuTOrwmapkhqpUv1RqNhFJb/T6zRZ1LmNWWikMd3H0HlJeX448//sBDDz3ktH3gwIFYu3at2+MvXLgQe/bswXvvvYennnrK43aNHDkSADB58mT7NkmSIISAJEmwWCweH8sRJxgQERFR8HmSjLTdX3Gow/Tp0/HEE09U2v3UqVOwWCzIzMx02p6ZmYljx465PMXu3bvx0EMP4eeff7aX3/DUvn37vNrfUwzWgiTca66x+5MildWWQVn0hzwe5tzJ0/49gb8zagqr/Bd/7qY/8U4juZZX1wZpAICGKdE+r2rgr55Af6yjXlIuP8dR89bJG47urOakttfJHxm2mrA92ePHi/CnbaWJpqly1jM5Vp21x3yhgQRNNfls5f5Dhw4hISHBvr26zHrFlQ+UTFdFFosFN9xwA5588kk0b97c06bbZWVlef0YTzBYIyIiopCSkJDgFKy5k5aWBq1WWymLduLEiUrZNgAoLCzExo0bsWnTJtx5p5x0sFqtEEJAp9Ph+++/R//+/as97/bt23Hw4EGUlzvPNL788surfawrIR2sMXtDgRDuWU8KrCOn5TFqH/y0X95gK2iqCp4MoC84ifVr5IzT/kuaAQBS4wyINcrj1jwtnRDoelje1p01WwVOFsrj8nJXL/fyZJqaZ9dqUijXIK8ocepYHpb9KRcxvjRbznrWTYqCsKXWQj3BFog6awaDAR06dMCKFStw5ZVX2revWLECV1xxRaX9ExISsGXLFqdtc+fOxcqVK/HJJ58gOzu7yvPt3bsXV155JbZs2WIfqya3W244x6wRERERVXDvvfdi1KhR6NixI7p164a33noLBw8exG233QZAnk1+5MgRvPPOO9BoNGjdurXT4zMyMhAVFVVpuyt33303srOz8cMPP6BJkyb4/fffkZeXh6lTp3Jt0FAWblkcZjspklmFwKGCYgCAyWT7C/pcnn8OXpul2w/IJQbe+aM9ACA7JRYGnbwElUFXdRanNpvpyTg2pfSIyWLFgP9+XwutqoKvM0Mt8jqzZeeKcfasXArmiG2W8QXWeOi9bYbtp9pK6gZqbdCRI0ciLy8PM2bMQG5uLlq3bo2lS5fax5fl5uZWW3PNU+vWrcPKlSuRnp4OjUYDjUaDnj17YubMmZg8ebLPNdgYrBERuaFczDy9PpgtAgkG+dJpUMpdWH3r9nBS0wjIx/pha9bJM9sOd26AlDgDALl+nNykym0K9lJArmIh5T+PnS3F6T/cl2qo/uBBnGygPCGLBfHx8kB6pTvaHGarGATKpEmTMGnSJJf3LVq0qMrHPvHEEy5nmrpisVgQFydP/khLS8PRo0dxwQUXICsrC//88483TXbCYI2IiIiCzpPlpHxdbqq2tG7dGps3b0aTJk3QpUsXzJo1CwaDAW+99RaaNGni83H9WK6ZaiIcug/D4TkQOfK24nthqRlxBh3iDDqUlJhQUmJbezOYalCV/8yG1TizYTWeX7kH+cUm5BebXK5qEOgVB7zl2BalvR3v/EA964B6+4Ip7S44AZPJApPJYv9slpvlmYrerGRgb4bXjwisqlazqI2VLfzh0UcfhdUqZ1+feuopHDhwAL169cLSpUvx6quv+nxcZtaIiIiI/GDQoEH2/27SpAm2b9+O06dPIzk52eXQAU8xWCMi8hOTxYrEGHnMmsYflUpVkkbYsvkIDvZvCgBIipWfn06rUXUxVqU48aE8ecIH8lwsLeULf45d87ach6kUCQnymDXlUeXm8Fl2SgMPukFVlw+sXkpKSo2PwWBNRUJ1Zii7Pync1ORysP+0vO7k0QPHa9CAAK1S4KOzmzfg5dXy8j7PX94KABCt10JS6nup8PqpDLwf+txKeUPRWf+eIBgTDkoKcfz4OQDAVtvP5qnx5+M9Fb4P5B8M1oiIiCjoAlEUN1yEZLDGTA7VhlDNdFLw5BebcM/ivwAA5rN+XhM0mMpL8MsauezAmosyAABXxddHjK08iVZlV1AhgM2H8gEAp9atDHJrquFN7TWNFsePnAIAfPmbXEqlW4MkNEqVVzjQhnhqTYPqZz1G6qzISH3eRERERCEhJDNr4S5UMjrMcBI5O1RQjDN58lgiFJ/17SAqy1Ipynf9AQB49xe56nv3BqlokCJndKKVAsBBpiSnzFaBG15ZUzsn9ce6oR6fS0KZbeWCA3vkMZG5RU0RLnVxJUmqdsZkTWZUhjIGa0RENaTMyEuLMkKnt32tBquelx8nFriyY4s8q3Jz/yb2ma/eBmveXnA9rSGmzAA9eqYEeb/XUrBWm4SwR6SJqYkAgFi9rlLdO48P57eGUaAxWCMiIqKg86QIdWTm1RisqVpypztV2xXKLlCi8wpK5IW2lUwTAKCsOEitCRBbRqfk0F4AwOxldYHB8l2Xt6oHANBpJLe9uDXpvpIkqcrsmnLXyUI5m3nxre/4Z01WT/mjjIcnEw2EFSgpBABER8uX7zi93p7ZjVFJd7SvwmG5qUDhBAMiIiIiFWNmTeXUNtmAGTWKBN7WGLXYHqDTaqDR2v4Grs3MTm0qKQAA7Nx8AE8XyZmsVmkJAIAmGbEw6JxzAP4aEF7xOI6ZNpNFzixdOPotecPxPX45Z63yZFyeEEChXLojNlYu3aGRzhcADofauKHc9kBisEZEVEPKBcZktvq+zJTKVi1wyzZxwpKXi2NGudv3j2NnAACp8UakxxsC3wYbZWD9nuPyqhEhGaR5y2IGABiN8uX7nMns9QQDTiwIPQzWiIiIKOi4goF7DNZCRLC7Q9n9SZHI024lJbOxK68QRQVhNrGgIqW7ruAkSsvkjNbry+oAAGL1WgxuURcA7N2h/r64Kt2fVgEczy8DAPQYp45hIgEnSfZJDCdOyPX8thw/hwvS4m07GIPUMP9gnTX3OMGAiIiISMVCLrMW6Rme2i7nEemvN6DuEipUO9yN8VH+xs87J4/jevH73cg/sL82mhR8VgtQKmd3/lmzHgBwz5EzMN7ZAwDQq0kaACA+Wu/68T5SsphnikxoPfghvx5b9YSwT1w5ekBeweCLdQZ0b5gEAMhKiwEQutknrg3qXsgFa1R7XaIM1IiqpgRxeSVyd9zZs6X22XoRwd4legIAkP93Hh7/WK6sP//mTgCA1g3lf2ukmgURSpBWUi4HK82ve83nY4U02yQS8zk5UD59uginS+U/FpR5Br7OcSH1YrBGREREQccxa+4xWAthgcqwMaNG5J1onVw5vrzcApSXBLk1QaBk2KwW7Fn1EwDg2nz5dZhzaxcAQK+cNET5WGHfbLEiv0QuWdF04MO2jQ5rr9ZGyRK1sX3OdDotdBr5+Vs9TK0p96qthAeXm3IvAj/hRERERKGDmbUw4K8MGzNqRN4pM8njp5TMWmlJefiuXOBKxbUwBeyTDk6uXw0AuNPWbfX0mItxZZv6AACth4OqzLaVCXLPlqLtkIddn9NxWyRl2GyZRZ1OA63t5bR4sgqCirEb1L0I+mQTERERhR5m1sJITUpMMKtG5D0ljxGllzNrJeciaLyaqwyX43Zb5ufEpg0AgAfKTYi+vRsAoGe2XNYjxqCFvsJaokIImG2Lre47KRfd7XLFI961qbYzbO5ei0BQMrdmeQayEIBkG8nl7bJTasPSHe4xWAszroKuigEcAzMiP7FdG/NLTQCAojNng9cWQA4aAh2oeBuYFMnrhp7ZsBo37dgCAJDSGwIAclo2wuwb2wEAGiTJNcJ6P/otCo7byp8c3el7G8O1S1QJ1koKAQAFZ8/hRLEcuClBbnVCO6SLTAzWiIiIKOg4Zs09BmsRgJk0osBQrht78uVB9Qjxbqgq1bSrT1jtBYOF7ee/ezfhsm+/rGnL3J/PUaAybbXZBerI1s2s0Wiw7bi8Hu2g5lV//tT+6WTpDvfCNE9MREREFB5CJrPG7BAFU20t8UWhRVkT9Okl2+QN5cVBbI1NIAbZByt75E9qHcdWw3IbZ06cwddr9gEAJnSUxwKmxhlq3KxgkKTz2eqq9olEIROsERGpTVGZPNi7zLY2Iyxm3w+mXLT9dTWqSdAWDsGZK/4IZNXy2uhsAVlZMc7ly7NmS00qaRv5HYM1IiIiCjoNJGiqGZVW3f3hisEaEZGPlAr7FrNS+6q8ir2DxDETVFVGSS0Zo9rgS4ZNba+PxrbOqrnc/vlTPo+hit2g7qmwA5+IiIiIFMysERH5SFnjUijjzcqKan5QIQKXPlBbdijYgv161GRygdL2smJ7PYtQX8FAsv2vun0iEYM1IiIfaSoGVeURtNwU+c4fC647/IEQmxALADDZgjUhRMQWjw1XDNaIiIgo6DhmzT0Ga0REPrLashtWfw/s9ncZD1IHf2TU7MeyfeasFhgM8qXc6nB8UeFcoZBpkzyYDRqp3aCcYEBERESkYsysERF5QclYCAHotPLfu1arLcvh7wHrzLCRO8pnIyrOo6yZY6ZNrVk2doO6x2CNiMgDjkGaQpkNWlpcGuiTR+5VKlz4swvU6bhWlJfLK2ecM5ntp6rq4xIKgRs5Y7BGREREQcfMmnsM1oiIqiGEcJkYKSqVMxnm0yeVHQPZCPlnpF6tQlUAM2oAgJJCFJwuAADsOSPX+euDdMDDj4uSZauYYas4QcHT+2qCddbc4wQDIiIiIhVTfWYtudOdwW4CkZ3yeTyz4fUgt4Rqg6txavb7AAyesUz+x8n9tdYmZthCRCCzrI7Ht1pQsnc7AOD5D6MBAGM7ZZ0v2Oxlhi2YNJJ8q26fSKT6YI2IqLZ5cuGSABQdORj4xrjDSQfqVduBj6kMAHB6xzbb+S+3L0Flbwk/LiGNwRoREREFHcesucdgjYjIpmJGrdoEiS2jQQSg9jNq9vPaJhuUngvO+f2Es0Hd4wQDIiIiIhVjZo2IIlqNBlZrg/wVyskG6hDMwfnCev78yfXOb65QkkPY/0+9HxcJ1XdzqrTpAcfMGhEREZGKMbNGRBHHb2UKdAb/HKemmGELnmCXvBACsJjk/46KlTchNDNQLN3hHoM1Igpb8soDAbyYSirrnHB8rqEQuLl7/ZQB82oW7CDNkT1Yd3XX+e5QpcVq/WRwNqh7KvumISIiIiJHzKwREfkqFDJAauJpJrLifmp6ndWUUatAo5Fft1DNPbF0h3vMrBERERGpGDNrRES+slqC3QL31DLpwB/j+iSNurJrKqVk1lwRQtjLeKiVhOqzgup+BoHDYI2IyFdqDtaCzd+TL5TjMWirnpuIRqglgCevMVgjIiKioNNAgqaaQFITobk1BmtERL6ymIPdAvUJdDkTZtjcMhecqfJ+doOGLk4wICIiIlIxZtaIiNxQEhEqrtagLrVZJDhYkw5U/KEwpqRVeb/qs1JMrbnFYI2IyFfKMj9qpPIuL78IZpeoCoO2slPH3d4nSVJEfCTCFYM1IiIiCjouN+UegzUiomo4ZiSURIoAOMhdobY1UiOViyyfMqkgJEIcD1YwCI0n4n/8DSMiIiJSMWbWiCgiuSpjIDwYf6Q8TAIAjda/jfKHSByYFOyxayoZt2ZIywRgG6df4XMQCh8Lzi9wj8EaEYUteVC151/vFfetNnhTyUU6aNj9KQv2ZAOdAQBQfrZynbVIDW7CDYM1IiIiCj6m1txisEZE5IZjps1llk1tfUtqa09tC/bqBsHKsCnPN79y6Y5Q+khwNqh7zGETERERqRgza0REHlCybEqGrbjcAuiMwWxScHG8mvoYYwHIZWVc5Z8qfobVRvKgdEcoZQr9ib9tRERecJq0EB0v3zy5ygS2UbV3fkmj/kAt2O2r7c+CVg9o9dBl1Icuoz4AQCPJN+XzqvZF3ANt7ty5yM7ORlRUFDp06ICff/7Z7b6fffYZLr30UqSnpyMhIQHdunXD8uXLa7G1lXn8iS4rK0NBQYHTjYhIDfj9RBT6JA9v3lqyZAmmTJmCadOmYdOmTejVqxcGDx6MgwcPutx/zZo1uPTSS7F06VL88ccf6NevH4YNG4ZNmzb59Lz8weNgbebMmUhMTLTfGjZsGMh2ERF5rDa/n4QQEELAbBHQxsZDGxsfsHMRVUuSAK0O0OoQExeDmLiYauc3qDbTFqBo7aWXXsL48eNxyy23oGXLlpg9ezYaNmyIefPmudx/9uzZeOCBB9CpUyc0a9YMzzzzDJo1a4avv/7at+flBx4Haw8//DDy8/Ptt0OHDgWyXUREHuP3ExG5Ul5ejj/++AMDBw502j5w4ECsXbvWo2NYrVYUFhYiJSUlEE30iMcTDIxGI4zGCB5MS0SqVZvfT0rWotRkgc4gf4VaauXMLtRmdiTY48C8FUllPKzyJ7C4sNh2SndTDNTNm9IdFYc6uPsOOHXqFCwWCzIzM522Z2Zm4tixYx6168UXX0RRURFGjBjh0f6BEGK/fURERBTpGjZs6DT0YebMmVXu72p1Ek+6gj/44AM88cQTWLJkCTIyMmrU5ppg6Q4iIi8o3+/JsQYYo+S/5MvsWZxazLGpccyRGgU7w1YbbM/RfFouiqvK8Wge8KZ0x6FDh5CQkGDf7i6znpaWBq1WWymLduLEiUrZtoqWLFmC8ePH4+OPP8Yll1xS/RMIIGbWiIi8oAzONug0aNaiDpq1qBP80h2BFmpdoGoS6M+GEPYJBnVat0Kd1q2cSnZU3TR1fWa9mV+QkJDgdHMXrBkMBnTo0AErVqxw2r5ixQp0797dbVs++OADjB07Fu+//z6GDh1aw2dWc8ysERERUdi69957MWrUKHTs2BHdunXDW2+9hYMHD+K2224DIE9QOnLkCN555x0AcqA2evRovPLKK+jatas9KxcdHY3ExMSgPAfVB2tnNrwOAEjudGeQW0J0/vNIBABPXtYSAHDZx3p5g8Uc+JOqLBsSMiRNeHaFarRARjYA4JlRFwEAtJoQ/YwEaCH3kSNHIi8vDzNmzEBubi5at26NpUuXIisrCwCQm5vrVHPtzTffhNlsxh133IE77rjDvn3MmDFYtGiR9w3wA9UHa0REREQ1MWnSJEyaNMnlfRUDsJ9++inwDfISgzUiIh+1qW/rEtEZ5J/lJcFrTCCE21i1cJxsIEmAxQQAaJ2RZNvkefpJTeuFelO6I9IwWCMi8lGs0fYValtAG8X5gTsZuz9DmyQFpuaa1YLEjDQAQHKs3ufDqCloo8oYrBEREVHQeVO6I9IwWCMiqqmouGC3wL/CrfuzonDqDhUC+cdOAgCMutB+3wI0vyAshPY7S0RERBTmmFkjIvKRMr5Hn5gEADDlBrExFLlMpQAAjR/6CIM6do2pNbcYrBER+chilS9oTVs2AADs+GdD7SzcHUiB7AKtGExU91ppbZeoQNWvC5Paay26hXh9NaoWgzUiIiIKOpbucC9kgjWuZEDBxJULyBUlk/Hfq1oBAK75dln41Vrzl6g4xDRtDQAoPlcsbzt1+Pz9JQXyT4sZUracKUpIkRfqzv93F5B/PDDtCvXJBvGpeGBYcwCATuu/QEaSpFrvCuVsUPc4wYCIiIhIxUIms0ae8TUDxIwlkfeUwdh9m2fIG2KTQzez5q+xakrqw7aqg9SgBQDgok45uHdQM3mbbdfmqfGIi5IvQyaLnMU5WVCG/YVFAABlCNaSP7KwfdcpAMDhLTvljfnHQ/e19geNFgAQ26QlemWnA/Bu5QI14vwC9xishQF/dNE5HoOBG5GP4lOAM0cDc2ylS0rNF2RJApLqAgAyWl0IAPji3j4AgMQYPeJtgZnGFoVJOD+DUelwS47VIydTXhHCapvAcXGdZBzrJc94fPbHZADAb7/uRvk/G+UH2QIXWC2BeV5qZHuuxeeKEWPQBuQUXNVAPRisERERUfAxteYWg7UQFOjB7syyEfkmqU46zh4MdiuCoE5TAED7vhdhwaiOAIDEaPnyEmXL+ri6xkrS+bl9wr5Ncz6JaNsWa9QhI9EIAJg3Qp58kDvwAsz+RZ7Y8eMquWu0ZNt6359DoCcaBCg7ZYgyQKeV2+7qNfbHWWsrw8bZoO5xggERERGRioVcZu3MhtcjMtsTrNIRFc/L157IvYWTe+PKe/fL/zixL6htCRhblkVq3BatLm4CAJg/ugMAoG5SFKL0WsfdPD+s7ae7Kvw62x4ZCXKGLTXOgBevkDNrR3vJ7Zi5Sp7AsGbNLpRu/927BoSa6HgAwK/PDa+yGG7FzGVNKBm2QE1kYOkO90IuWIskagwS2EVKVJly/ejVLA0NL5IDiEM/7Jc3hsPgbEkCDNEAgOy+/QAAz11/ERonyxMBGqTI9+k0Uq1dTLUaCSmxegBAgm3iwjNDWgIA9nRtiBlf1gMAbPvuR/kBpefC472wqdOlp/wzKcqjjkEJ/gnYKDgYrBEREVHQcX6BeyEZrIX7agZqzKi5wveBIl3FC4cEYMb1bQAAN69bK288d7pW2+QXFWqlJbTuhJzmmQCAl69uCwBonB4Lo14e9qyrxTUpXXXB6XXyNiXDl55gxMyr5dUSnk2MAgD8u/sUTmz+S37AuTz5Z21k2vx9jrgUAMALY9sDAIy6MBp6zmjNrTB6l4mIiIjCT0hm1sJRKGdxOI6N6LxLmskZqIz2nQAAJ375wb/FWv1dHFcpWaEUlpUkQC9noxp0l8dFTb+uDbo2TAUgD+wHAL1WqnKgub/LPFQ8l6sz620lLPRaDS5uKBfPffEKOcO2+cRZPP+lnHn7d82v8gOK8x0arJTs0Kh3nVBJQsNu3QEAPbPTALifkBGKWLrDPWbWiIiIiFQspDNr4TBmKpQzaq7wPaFI4e7ve2XpnwdHymPXpu45CBzZUUut8pJGez6jZhufhtRGaNHxAgDA3BvlcVFNMmIRbXteyvC06so3eFreIVCFVpX3oUmGPGM1M9GI+jfKmbWHbLNHt/6xBzhlq2KsrDNqMQNKYq0mGbZAPK+6zfHYCHm2cawxMEtMBZUHpTsiNLEW2sFaqIqEYIBdoxTprmlTHwDwettm2BeIYE0I37tCHbs+o+IAAPU7dwUAzB3XEY2SYwCcr2mm12rs63n6W6C7UpUaZAnRerRvJHeNzrvhYgDAb32z8PZPBwAAW3/5S36A4wLxKqt1Ua9FDoa0lNde9aX705811wKB8wvcYzcoERERkYqFRWYtVFY1iISMmiuh1DUaqe8R+V+cUf56ffb6trirqAwAcGLN8mA26XxGTWv76k+uhzZ95NUHFo2VJ0RkJhrt5SACVaneU/44v+MRlOfVtI6cTUxPMOKCFPm/H7TIXZ47ft8BnDkqP6CkwLeT+rkLVHeB/N58dW9v+woRYYmpNbeYWSMiIiJSsbDIrAHqzd4wU3Oemsex8X2iQOnbLANLH5CXaBoZL5fE2LP0q8AVZFWyZ46FbZVz2QbM5wz8DwBgxrWt0aGBPI5LKckR7GxabVBKfKTFG5EWL4/LW3SznL1ae2kTvLViLwBgx4qf5AeUFdnLmaD0XOUDBqrUR52mAIBtc68DIL9HNXl31DpWTcHSHe6FTbCmCHbQxou+Z9SyQDzfL/KVgGc9MjqNhOx0eUbizJFy9f878ktw8pfva9YAx6DKsXtT+W+9HIQgqS4S68m13x4dI3d5XndRQwCwz/CMFFUFKzmZcndok4xYXNZCHsS/aoi8MPxL3+7GoX0nAABFu7fIDzCXyz81WocJCbagTZJqFozXl9c4fe/pqwDAvgaqp2GK2oMy8l7YBWtEREQUeiQPSndEQOLXpbAN1mpz0gGzMzVX212kfM+oNlltWZb2DZIAAO/e2RO327pE9637Xd7p7LHzD1BqnykTASzm86sgKFcrSVNpDU8k1bU/JipFXkNy+ZOX2UtwxNnqi0VaRs0bkiQhOVZ+PQfbMmwd66Vg2wl5tYMZn8krBxw9eAoAUFpcCsvhf+QH2zNsAhAVVq1wWBnCnoEzl5/PwMXLK0TE5rTEd48PBgBk2NY1VYR7xozzC9zjBAMiIiIiFQvbzBoQmHFRzMgEnqvXmO8dqZGS6ag2G2DLgCXGyGOP2mcl4Zt7egEANl8nr3TwwHt/AwAObdgIJGYAAGKTEgAARSdOyMVagfMZNqsFiJOzZ4k58tiqZ27phK4N5MyPUsE/IVpnH1AfqV1IgG9ZqSi9/Lo1SIlGui07eeFtiQCAM0XymLV1R89gzhfpAIBju/fLDzybe/59KiuWfyakI7ZxU6fjF506g7SsegCAd++S12FtVS8BRlt5joh7u5hacyusg7WKeLEOXXzvSM28DdokCci0dXFdmlgHAPD7E3KAZrIMQEm5fKEvMcndZVargMFWI0xv+2nQSudrbtlOrJEkVFxoIBJmdwaK4/uldB03So12+tmmYSLGd8oCcH7FBbNVoNwsv3eltvewqNRs/5woQXtitM6+MoROo5yL7xdVFlHBGhEREakTS3e4x2CNiMhPajIAXO+QOavJBABmZgKr4uurgYBGq2yTfxpwvhvaLtFon0vg6i3i+2brBa1uNmittER9OMGAiIiISMWYWSMiUhlPsizCoegqszLuBbrchavXXgjhcjvfpqpxfoF7DNaIiEJQpAZoofCsI/W9ocBhsEZERERBxxUM3GOwRkREqheh1+gIw45QdzjBgIiIiEjFmFkjIiJViswcSuRiN6h7DNaIiEhVIvR6TOQWgzUiIiIKOo5Yc4/BGhERBU2gL77K8QNdb41qjt2g7nGCAREREZGKMbNGRER+F6EJEKoBLuTuHjNrRERERCrGzBoREflNZOY9yC84w8AtBmtEROQXEXodDXsV39dAvc+M1dxjNygRERGRijGzRkRhK1L/Cq9tfJ3JH1i6wz1m1oiIiIhUjJk1IgprElgQNVAiNMlBAcLSHe4xWCOisMcq9v4XmZfMyFLr7zFnGLjFblAiIiIiFWNmjYgiBjNsROrFxJp7zKwRERERqRgza0QUcTjpgKhqwchgsXSHewzWiCgisUvUNxF6rSQKKgZrREREpALVl+6I1D8XGKwRERERgOCGQuwGdY8TDIiIiIhUjJk1IopoHLsW3vi+UjhgsEZERBThIrR3MWQwWCMiIqKg45g19xisERGBtdeqE6HXyLCnpveVC7m7xwkGRERERCrGzBoREREFHbtB3WOwRkREYYdd2hROGKwRERFR0EmofgxdhCbWGKwRESlYc40iieoCH0ZrbnGCAREREZGKeRyslZWVoaCgwOlGRKQG/H4iR8yMhibJw//5Yu7cucjOzkZUVBQ6dOiAn3/+ucr9V69ejQ4dOiAqKgpNmjTBG2+84dN5/cXjYG3mzJlITEy03xo2bBjIdhEReYzfT0TkzpIlSzBlyhRMmzYNmzZtQq9evTB48GAcPHjQ5f779u3DkCFD0KtXL2zatAmPPPIIJk+ejE8//bSWW36eJITw6I+QsrIylJWV2f9dUFCAhg0b4nhePhISEgLWQKJIVFBQgMzUROTn8/fLE+6+n3JPnvXp9WNmprJQGCrE9807vr6nBQUFqJue5Lfvp4KCAiQmJuKoB7+vBQUFqOflubt06YL27dtj3rx59m0tW7bE8OHDMXPmzEr7P/jgg/jqq6+wY8cO+7bbbrsNf//9N9atW+fhs/IvjycYGI1GGI3GQLaFiMgn/H4i8pxaA29v5hdUHOrg7jugvLwcf/zxBx566CGn7QMHDsTatWtdnmPdunUYOHCg07ZBgwZh/vz5MJlM0Ov11bTS/zjBgIiIiEJKw4YNnYY+uMqQAcCpU6dgsViQmZnptD0zMxPHjh1z+Zhjx4653N9sNuPUqVP+eQJeYukOIqIK/FnCQwgBKVLLrtcSdn+GCS9Sa4cOHXLqBq0us17xd7C630tX+7vaXlsYrBEREVFISUhI8GjMWlpaGrRabaUs2okTJyplzxR16tRxub9Op0Nqaqrvja4BdoMSEQWYEMLpRv4hwKxaOAlE6Q6DwYAOHTpgxYoVTttXrFiB7t27u3xMt27dKu3//fffo2PHjkEZrwYwWCMiIhURLm4UGZSF3Ku7eevee+/F//73PyxYsAA7duzAPffcg4MHD+K2224DADz88MMYPXq0ff/bbrsNBw4cwL333osdO3ZgwYIFmD9/Pu677z5/PVWv+dwNqvx1WMjik0R+p/xeMQvjG/v3U2HNvp/8NWatolAcw1ZbLXb1mrs7N387fFPT91L5vfL395Mnxax9KXg9cuRI5OXlYcaMGcjNzUXr1q2xdOlSZGVlAQByc3Odaq5lZ2dj6dKluOeeezBnzhzUq1cPr776Kq6++mqvz+0vHtdZq+jw4cMsPEkUYIcOHUKDBg2C3YyQw+8nosDz1/dTaWkpsrOz3c7OrKhOnTrYt28foqKianzuUOFzsGa1WvHPP//gwgsvrDQrg3ynFPPka+ofofp6CiFQWFiIevXqQaPhaAVv8fspMEL190mtQvX1DMT3U2lpKcrLyz3a12AwRFSgBtSgG1Sj0aB+/foAPJ+VQZ7ja+pfofh6JiYmBrsJIYvfT4HF19S/QvH19Pf3U1RUVMQFYN7gn+xEREREKsZgjYiIiEjFahSsGY1GTJ8+nWvy+RFfU//i6xm5+N77H19T/+LrSZ7yeYIBEREREQUeu0GJiIiIVIzBGhEREZGKMVgjIiIiUjEGa0REREQqxmCNiIiISMUYrBERERGpGIM1IiIiIhVjsEZERESkYv8PQAHDiFxgRrUAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plot.subplots(nrows=1, ncols=2, figsize=(7, 6))\n", "plot.imview(\n", " tangle[32],\n", " title=\"Ground truth (central slice)\",\n", " cmap=plot.cm.Blues,\n", " cbar=None,\n", " fig=fig,\n", " ax=ax[0],\n", ")\n", "plot.imview(\n", " tangle_recon[32],\n", " title=\"TV Reconstruction (central slice)\\nSNR: %.2f (dB), MAE: %.3f\"\n", " % (metric.snr(tangle, tangle_recon), metric.mae(tangle, tangle_recon)),\n", " cmap=plot.cm.Blues,\n", " fig=fig,\n", " ax=ax[1],\n", ")\n", "divider = make_axes_locatable(ax[1])\n", "cax = divider.append_axes(\"right\", size=\"5%\", pad=0.2)\n", "fig.colorbar(ax[1].get_images()[0], cax=cax, label=\"arbitrary units\")\n", "fig.show()" ] } ], "metadata": { "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.12.2" } }, "nbformat": 4, "nbformat_minor": 5 }