{
"cells": [
{
"cell_type": "markdown",
"source": [
"# Gross-Pitaevskii equation with magnetism"
],
"metadata": {}
},
{
"cell_type": "markdown",
"source": [
"We solve the 2D Gross-Pitaevskii equation with a magnetic field.\n",
"This is similar to the\n",
"previous example (Gross-Pitaevskii equation in one dimension),\n",
"but with an extra term for the magnetic field.\n",
"We reproduce here the results of https://arxiv.org/pdf/1611.02045.pdf Fig. 10"
],
"metadata": {}
},
{
"outputs": [],
"cell_type": "code",
"source": [
"using DFTK\n",
"using StaticArrays\n",
"using Plots"
],
"metadata": {},
"execution_count": 1
},
{
"cell_type": "markdown",
"source": [
"Unit cell. Having one of the lattice vectors as zero means a 2D system"
],
"metadata": {}
},
{
"outputs": [],
"cell_type": "code",
"source": [
"a = 15\n",
"lattice = a .* [[1 0 0.]; [0 1 0]; [0 0 0]];"
],
"metadata": {},
"execution_count": 2
},
{
"cell_type": "markdown",
"source": [
"Confining scalar potential, and magnetic vector potential"
],
"metadata": {}
},
{
"outputs": [],
"cell_type": "code",
"source": [
"pot(x, y, z) = ((x - a/2)^2 + (y - a/2)^2)/2\n",
"ω = .6\n",
"Apot(x, y, z) = ω * @SVector [y - a/2, -(x - a/2), 0]\n",
"Apot(X) = Apot(X...);"
],
"metadata": {},
"execution_count": 3
},
{
"cell_type": "markdown",
"source": [
"Parameters"
],
"metadata": {}
},
{
"outputs": [],
"cell_type": "code",
"source": [
"Ecut = 20 # Increase this for production\n",
"η = 500\n",
"C = η/2\n",
"α = 2\n",
"n_electrons = 1; # Increase this for fun"
],
"metadata": {},
"execution_count": 4
},
{
"cell_type": "markdown",
"source": [
"Collect all the terms, build and run the model"
],
"metadata": {}
},
{
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Iter Function value Gradient norm \n",
" 0 3.148120e+01 7.976782e+00\n",
" * time: 0.0021741390228271484\n",
" 1 2.860931e+01 4.834721e+00\n",
" * time: 0.006232023239135742\n",
" 2 1.878532e+01 4.231944e+00\n",
" * time: 0.010436058044433594\n",
" 3 1.305985e+01 3.290575e+00\n",
" * time: 0.014642000198364258\n",
" 4 1.154349e+01 3.492442e+00\n",
" * time: 0.018026113510131836\n",
" 5 1.133521e+01 3.110641e+00\n",
" * time: 0.021433115005493164\n",
" 6 9.973911e+00 1.461407e+00\n",
" * time: 0.025832176208496094\n",
" 7 8.825201e+00 1.696388e+00\n",
" * time: 0.029942035675048828\n",
" 8 8.340630e+00 1.257247e+00\n",
" * time: 0.0341341495513916\n",
" 9 8.188388e+00 7.369264e-01\n",
" * time: 0.08735418319702148\n",
" 10 8.181798e+00 1.369906e+00\n",
" * time: 0.08998298645019531\n",
" 11 8.143823e+00 1.008025e+00\n",
" * time: 0.09273409843444824\n",
" 12 8.070573e+00 5.562346e-01\n",
" * time: 0.09541201591491699\n",
" 13 8.019478e+00 5.871768e-01\n",
" * time: 0.09789013862609863\n",
" 14 7.987270e+00 5.179631e-01\n",
" * time: 0.10050106048583984\n",
" 15 7.949139e+00 4.413989e-01\n",
" * time: 0.10306119918823242\n",
" 16 7.921160e+00 2.473208e-01\n",
" * time: 0.10559701919555664\n",
" 17 7.903649e+00 4.275275e-01\n",
" * time: 0.1083061695098877\n",
" 18 7.887734e+00 2.483476e-01\n",
" * time: 0.11151003837585449\n",
" 19 7.877275e+00 2.716477e-01\n",
" * time: 0.11405205726623535\n",
" 20 7.863731e+00 3.190057e-01\n",
" * time: 0.11655211448669434\n",
" 21 7.851747e+00 1.992245e-01\n",
" * time: 0.12076020240783691\n",
" 22 7.845183e+00 2.990783e-01\n",
" * time: 0.1255180835723877\n",
" 23 7.831563e+00 2.079099e-01\n",
" * time: 0.12980914115905762\n",
" 24 7.819596e+00 2.182718e-01\n",
" * time: 0.13413715362548828\n",
" 25 7.813778e+00 2.517873e-01\n",
" * time: 0.13860416412353516\n",
" 26 7.803472e+00 2.572012e-01\n",
" * time: 0.14226698875427246\n",
" 27 7.791507e+00 1.526725e-01\n",
" * time: 0.14589715003967285\n",
" 28 7.783590e+00 2.070923e-01\n",
" * time: 0.14956998825073242\n",
" 29 7.774926e+00 1.744098e-01\n",
" * time: 0.1532430648803711\n",
" 30 7.771718e+00 2.111901e-01\n",
" * time: 0.1575150489807129\n",
" 31 7.764113e+00 2.157635e-01\n",
" * time: 0.16212916374206543\n",
" 32 7.757313e+00 1.946915e-01\n",
" * time: 0.16663503646850586\n",
" 33 7.749997e+00 2.537610e-01\n",
" * time: 0.17126011848449707\n",
" 34 7.746078e+00 2.994609e-01\n",
" * time: 0.17576217651367188\n",
" 35 7.736647e+00 1.872230e-01\n",
" * time: 0.18019604682922363\n",
" 36 7.733542e+00 2.962243e-01\n",
" * time: 0.18457603454589844\n",
" 37 7.723686e+00 1.600278e-01\n",
" * time: 0.18925213813781738\n",
" 38 7.716906e+00 1.818841e-01\n",
" * time: 0.19369006156921387\n",
" 39 7.709998e+00 1.974226e-01\n",
" * time: 0.19814014434814453\n",
" 40 7.703846e+00 2.382428e-01\n",
" * time: 0.20284318923950195\n",
" 41 7.696816e+00 1.051032e-01\n",
" * time: 0.20855212211608887\n",
" 42 7.694500e+00 2.571276e-01\n",
" * time: 0.2129380702972412\n",
" 43 7.690573e+00 2.034471e-01\n",
" * time: 0.23691201210021973\n",
" 44 7.685807e+00 2.182316e-01\n",
" * time: 0.23952603340148926\n",
" 45 7.681004e+00 9.912048e-02\n",
" * time: 0.24279117584228516\n",
" 46 7.676336e+00 1.797902e-01\n",
" * time: 0.24538016319274902\n",
" 47 7.671407e+00 1.006924e-01\n",
" * time: 0.24869799613952637\n",
" 48 7.666702e+00 1.548438e-01\n",
" * time: 0.25131702423095703\n",
" 49 7.663790e+00 1.384494e-01\n",
" * time: 0.2539050579071045\n",
" 50 7.660523e+00 1.402333e-01\n",
" * time: 0.2564401626586914\n",
" 51 7.658122e+00 1.531875e-01\n",
" * time: 0.258991003036499\n",
" 52 7.653058e+00 1.387979e-01\n",
" * time: 0.2615470886230469\n",
" 53 7.649242e+00 1.176318e-01\n",
" * time: 0.26500821113586426\n",
" 54 7.646228e+00 7.116394e-02\n",
" * time: 0.268233060836792\n",
" 55 7.642844e+00 9.009411e-02\n",
" * time: 0.2708010673522949\n",
" 56 7.641380e+00 1.094158e-01\n",
" * time: 0.27336812019348145\n",
" 57 7.641047e+00 1.180147e-01\n",
" * time: 0.2758660316467285\n",
" 58 7.640899e+00 8.176400e-02\n",
" * time: 0.27842020988464355\n",
" 59 7.640045e+00 1.683780e-01\n",
" * time: 0.28124117851257324\n",
" 60 7.638125e+00 8.808581e-02\n",
" * time: 0.28443312644958496\n",
" 61 7.636426e+00 1.047422e-01\n",
" * time: 0.2869300842285156\n",
" 62 7.635526e+00 8.817420e-02\n",
" * time: 0.28956103324890137\n",
" 63 7.633994e+00 5.655445e-02\n",
" * time: 0.29207301139831543\n",
" 64 7.633591e+00 8.076289e-02\n",
" * time: 0.2963531017303467\n",
" 65 7.632499e+00 5.568526e-02\n",
" * time: 0.30078601837158203\n",
" 66 7.631722e+00 4.675207e-02\n",
" * time: 0.30649304389953613\n",
" 67 7.631436e+00 9.169248e-02\n",
" * time: 0.31104016304016113\n",
" 68 7.630670e+00 7.257136e-02\n",
" * time: 0.31676602363586426\n",
" 69 7.629986e+00 6.682596e-02\n",
" * time: 0.32233405113220215\n",
" 70 7.629326e+00 5.659019e-02\n",
" * time: 0.3282191753387451\n",
" 71 7.629080e+00 5.506472e-02\n",
" * time: 0.3325660228729248\n",
" 72 7.628601e+00 6.706456e-02\n",
" * time: 0.3370530605316162\n",
" 73 7.627928e+00 3.153190e-02\n",
" * time: 0.34305405616760254\n",
" 74 7.627421e+00 4.915726e-02\n",
" * time: 0.3477320671081543\n",
" 75 7.627369e+00 5.874008e-02\n",
" * time: 0.3653261661529541\n",
" 76 7.626745e+00 5.288012e-02\n",
" * time: 0.3679921627044678\n",
" 77 7.626244e+00 6.666101e-02\n",
" * time: 0.37061214447021484\n",
" 78 7.625692e+00 3.658142e-02\n",
" * time: 0.3741281032562256\n",
" 79 7.625420e+00 5.861060e-02\n",
" * time: 0.37681007385253906\n",
" 80 7.625175e+00 4.771474e-02\n",
" * time: 0.37949514389038086\n",
" 81 7.624876e+00 3.696046e-02\n",
" * time: 0.38205718994140625\n",
" 82 7.624821e+00 6.071859e-02\n",
" * time: 0.38471198081970215\n",
" 83 7.624447e+00 3.546631e-02\n",
" * time: 0.38793206214904785\n",
" 84 7.624115e+00 5.440242e-02\n",
" * time: 0.3907740116119385\n",
" 85 7.624091e+00 7.003944e-02\n",
" * time: 0.39345407485961914\n",
" 86 7.624034e+00 5.958998e-02\n",
" * time: 0.3959810733795166\n",
" 87 7.623693e+00 5.046968e-02\n",
" * time: 0.39850616455078125\n",
" 88 7.623268e+00 4.214160e-02\n",
" * time: 0.4010441303253174\n",
" 89 7.622815e+00 4.709536e-02\n",
" * time: 0.4035520553588867\n",
" 90 7.622685e+00 6.505086e-02\n",
" * time: 0.40639305114746094\n",
" 91 7.622286e+00 4.334952e-02\n",
" * time: 0.40963315963745117\n",
" 92 7.622043e+00 5.810219e-02\n",
" * time: 0.41219210624694824\n",
" 93 7.621672e+00 4.556278e-02\n",
" * time: 0.41547513008117676\n",
" 94 7.621499e+00 4.837781e-02\n",
" * time: 0.4181091785430908\n",
" 95 7.621024e+00 4.616776e-02\n",
" * time: 0.42082810401916504\n",
" 96 7.620702e+00 3.447455e-02\n",
" * time: 0.42513203620910645\n",
" 97 7.620201e+00 2.457585e-02\n",
" * time: 0.43077898025512695\n",
" 98 7.620025e+00 3.474768e-02\n",
" * time: 0.4352850914001465\n",
" 99 7.619919e+00 3.446584e-02\n",
" * time: 0.4399731159210205\n",
" 100 7.619796e+00 4.066059e-02\n",
" * time: 0.44463419914245605\n",
" 101 7.619711e+00 3.385564e-02\n",
" * time: 0.4492671489715576\n",
" 102 7.619422e+00 3.535886e-02\n",
" * time: 0.4538850784301758\n",
" 103 7.619145e+00 2.939781e-02\n",
" * time: 0.4594600200653076\n",
" 104 7.618927e+00 5.123870e-02\n",
" * time: 0.4640212059020996\n",
" 105 7.618819e+00 3.919593e-02\n",
" * time: 0.4678230285644531\n",
" 106 7.618702e+00 2.448087e-02\n",
" * time: 0.47057008743286133\n",
" 107 7.618486e+00 2.572167e-02\n",
" * time: 0.4736900329589844\n",
" 108 7.618325e+00 2.927185e-02\n",
" * time: 0.48492908477783203\n",
" 109 7.618110e+00 4.595612e-02\n",
" * time: 0.4875810146331787\n",
" 110 7.617931e+00 2.446269e-02\n",
" * time: 0.49077320098876953\n",
" 111 7.617770e+00 2.808868e-02\n",
" * time: 0.49405908584594727\n",
" 112 7.617557e+00 1.996933e-02\n",
" * time: 0.4972701072692871\n",
" 113 7.617396e+00 1.993439e-02\n",
" * time: 0.5000290870666504\n",
" 114 7.617339e+00 3.330451e-02\n",
" * time: 0.5026090145111084\n",
" 115 7.617299e+00 4.919919e-02\n",
" * time: 0.5051090717315674\n",
" 116 7.617074e+00 3.974564e-02\n",
" * time: 0.5086400508880615\n",
" 117 7.616943e+00 5.094035e-02\n",
" * time: 0.5113182067871094\n",
" 118 7.616784e+00 5.996913e-02\n",
" * time: 0.5139870643615723\n",
" 119 7.616474e+00 3.445397e-02\n",
" * time: 0.5172979831695557\n",
" 120 7.616265e+00 3.202310e-02\n",
" * time: 0.5205211639404297\n",
" 121 7.616112e+00 2.351918e-02\n",
" * time: 0.5237431526184082\n",
" 122 7.615987e+00 1.983246e-02\n",
" * time: 0.5270512104034424\n",
" 123 7.615848e+00 2.571351e-02\n",
" * time: 0.5297491550445557\n",
" 124 7.615661e+00 2.813978e-02\n",
" * time: 0.5325310230255127\n",
" 125 7.615485e+00 2.614831e-02\n",
" * time: 0.5357010364532471\n",
" 126 7.615394e+00 3.707523e-02\n",
" * time: 0.5383419990539551\n",
" 127 7.615213e+00 2.412467e-02\n",
" * time: 0.5418450832366943\n",
" 128 7.615074e+00 2.228796e-02\n",
" * time: 0.5460531711578369\n",
" 129 7.614939e+00 1.987080e-02\n",
" * time: 0.5510921478271484\n",
" 130 7.614826e+00 2.483946e-02\n",
" * time: 0.5547940731048584\n",
" 131 7.614779e+00 3.045500e-02\n",
" * time: 0.5574660301208496\n",
" 132 7.614654e+00 3.239148e-02\n",
" * time: 0.5600371360778809\n",
" 133 7.614500e+00 1.972008e-02\n",
" * time: 0.5637071132659912\n",
" 134 7.614399e+00 3.007318e-02\n",
" * time: 0.5663430690765381\n",
" 135 7.614329e+00 2.113071e-02\n",
" * time: 0.5690150260925293\n",
" 136 7.614264e+00 2.539635e-02\n",
" * time: 0.5715830326080322\n",
" 137 7.614124e+00 1.723944e-02\n",
" * time: 0.5748600959777832\n",
" 138 7.614026e+00 3.091121e-02\n",
" * time: 0.5776140689849854\n",
" 139 7.613874e+00 3.204256e-02\n",
" * time: 0.5889320373535156\n",
" 140 7.613865e+00 2.837453e-02\n",
" * time: 0.591529130935669\n",
" 141 7.613794e+00 2.913205e-02\n",
" * time: 0.5942111015319824\n",
" 142 7.613635e+00 1.961048e-02\n",
" * time: 0.597520112991333\n",
" 143 7.613449e+00 2.332186e-02\n",
" * time: 0.6000330448150635\n",
" 144 7.613411e+00 2.707925e-02\n",
" * time: 0.6027050018310547\n",
" 145 7.613316e+00 2.222793e-02\n",
" * time: 0.6053330898284912\n",
" 146 7.613127e+00 2.787146e-02\n",
" * time: 0.6087100505828857\n",
" 147 7.613109e+00 2.205062e-02\n",
" * time: 0.6113481521606445\n",
" 148 7.613059e+00 2.701923e-02\n",
" * time: 0.6147921085357666\n",
" 149 7.613004e+00 2.499921e-02\n",
" * time: 0.6174790859222412\n",
" 150 7.612939e+00 2.929434e-02\n",
" * time: 0.6201090812683105\n",
" 151 7.612856e+00 2.152010e-02\n",
" * time: 0.6235010623931885\n",
" 152 7.612782e+00 1.519923e-02\n",
" * time: 0.6269180774688721\n",
" 153 7.612778e+00 2.687760e-02\n",
" * time: 0.6295640468597412\n",
" 154 7.612689e+00 2.270307e-02\n",
" * time: 0.6322581768035889\n",
" 155 7.612674e+00 2.574998e-02\n",
" * time: 0.6348800659179688\n",
" 156 7.612614e+00 2.244497e-02\n",
" * time: 0.6375851631164551\n",
" 157 7.612581e+00 2.725223e-02\n",
" * time: 0.64034104347229\n",
" 158 7.612580e+00 2.602264e-02\n",
" * time: 0.642963171005249\n",
" 159 7.612546e+00 3.133240e-02\n",
" * time: 0.6456341743469238\n",
" 160 7.612474e+00 2.197814e-02\n",
" * time: 0.6486630439758301\n",
" 161 7.612432e+00 1.938327e-02\n",
" * time: 0.6528830528259277\n",
" 162 7.612367e+00 1.817147e-02\n",
" * time: 0.6571981906890869\n",
" 163 7.612353e+00 2.780268e-02\n",
" * time: 0.660865068435669\n",
" 164 7.612277e+00 1.841349e-02\n",
" * time: 0.6642630100250244\n",
" 165 7.612221e+00 2.012426e-02\n",
" * time: 0.6676101684570312\n",
" 166 7.612184e+00 1.866585e-02\n",
" * time: 0.6711091995239258\n",
" 167 7.612158e+00 1.634046e-02\n",
" * time: 0.673821210861206\n",
" 168 7.612153e+00 2.320378e-02\n",
" * time: 0.6765201091766357\n",
" 169 7.612107e+00 1.624230e-02\n",
" * time: 0.6798629760742188\n",
" 170 7.612064e+00 1.172636e-02\n",
" * time: 0.6832690238952637\n",
" 171 7.612035e+00 2.047546e-02\n",
" * time: 0.6946229934692383\n",
" 172 7.612019e+00 1.578750e-02\n",
" * time: 0.6973211765289307\n",
" 173 7.611981e+00 1.632244e-02\n",
" * time: 0.6999650001525879\n",
" 174 7.611929e+00 1.023703e-02\n",
" * time: 0.7036080360412598\n",
" 175 7.611913e+00 1.406426e-02\n",
" * time: 0.7062761783599854\n",
" 176 7.611877e+00 1.575528e-02\n",
" * time: 0.7089290618896484\n",
" 177 7.611849e+00 1.741216e-02\n",
" * time: 0.7115120887756348\n",
" 178 7.611818e+00 1.809885e-02\n",
" * time: 0.7142801284790039\n",
" 179 7.611769e+00 1.009976e-02\n",
" * time: 0.717764139175415\n",
" 180 7.611752e+00 1.870613e-02\n",
" * time: 0.7204470634460449\n",
" 181 7.611719e+00 1.357782e-02\n",
" * time: 0.7231659889221191\n",
" 182 7.611683e+00 1.890593e-02\n",
" * time: 0.725836992263794\n",
" 183 7.611647e+00 1.296869e-02\n",
" * time: 0.7291650772094727\n",
" 184 7.611603e+00 1.016939e-02\n",
" * time: 0.7318680286407471\n",
" 185 7.611594e+00 2.030730e-02\n",
" * time: 0.734745979309082\n",
" 186 7.611565e+00 1.921945e-02\n",
" * time: 0.7374281883239746\n",
" 187 7.611562e+00 2.235582e-02\n",
" * time: 0.740070104598999\n",
" 188 7.611529e+00 2.030852e-02\n",
" * time: 0.7427651882171631\n",
" 189 7.611513e+00 2.291107e-02\n",
" * time: 0.7454700469970703\n",
" 190 7.611467e+00 2.460136e-02\n",
" * time: 0.748237133026123\n",
" 191 7.611432e+00 3.465308e-02\n",
" * time: 0.7509911060333252\n",
" 192 7.611358e+00 2.909797e-02\n",
" * time: 0.7543480396270752\n",
" 193 7.611278e+00 2.695200e-02\n",
" * time: 0.7583701610565186\n",
" 194 7.611234e+00 3.069625e-02\n",
" * time: 0.7625370025634766\n",
" 195 7.611189e+00 2.421526e-02\n",
" * time: 0.7662761211395264\n",
" 196 7.611102e+00 1.703635e-02\n",
" * time: 0.7697300910949707\n",
" 197 7.610992e+00 2.395691e-02\n",
" * time: 0.7724390029907227\n",
" 198 7.610960e+00 3.960296e-02\n",
" * time: 0.7751190662384033\n",
" 199 7.610933e+00 4.963861e-02\n",
" * time: 0.7778241634368896\n",
" 200 7.610807e+00 3.442356e-02\n",
" * time: 0.7813971042633057\n",
" 201 7.610700e+00 2.987762e-02\n",
" * time: 0.784080982208252\n",
" 202 7.610486e+00 3.391001e-02\n",
" * time: 0.7874701023101807\n",
" 203 7.610249e+00 2.099676e-02\n",
" * time: 0.7908420562744141\n",
" 204 7.610182e+00 2.452088e-02\n",
" * time: 0.8022661209106445\n",
" 205 7.610052e+00 2.823041e-02\n",
" * time: 0.8050520420074463\n",
" 206 7.609936e+00 2.002126e-02\n",
" * time: 0.8076980113983154\n",
" 207 7.609828e+00 2.505622e-02\n",
" * time: 0.8112761974334717\n",
" 208 7.609678e+00 4.460734e-02\n",
" * time: 0.8140101432800293\n",
" 209 7.609450e+00 3.118112e-02\n",
" * time: 0.8174221515655518\n",
" 210 7.609296e+00 3.916590e-02\n",
" * time: 0.8208181858062744\n",
" 211 7.609232e+00 4.030161e-02\n",
" * time: 0.823523998260498\n",
" 212 7.609069e+00 4.536255e-02\n",
" * time: 0.82625412940979\n",
" 213 7.608879e+00 4.431682e-02\n",
" * time: 0.8297741413116455\n",
" 214 7.608644e+00 3.192863e-02\n",
" * time: 0.8324601650238037\n",
" 215 7.608344e+00 3.273918e-02\n",
" * time: 0.8350741863250732\n",
" 216 7.608053e+00 2.972049e-02\n",
" * time: 0.8384401798248291\n",
" 217 7.607853e+00 3.100870e-02\n",
" * time: 0.8411750793457031\n",
" 218 7.607648e+00 3.496789e-02\n",
" * time: 0.8447620868682861\n",
" 219 7.607497e+00 2.572620e-02\n",
" * time: 0.8481550216674805\n",
" 220 7.607365e+00 2.669537e-02\n",
" * time: 0.851499080657959\n",
" 221 7.607231e+00 4.274631e-02\n",
" * time: 0.8541522026062012\n",
" 222 7.607020e+00 2.798904e-02\n",
" * time: 0.8574941158294678\n",
" 223 7.606865e+00 2.969221e-02\n",
" * time: 0.8602981567382812\n",
" 224 7.606680e+00 2.520878e-02\n",
" * time: 0.8638100624084473\n",
" 225 7.606672e+00 4.955680e-02\n",
" * time: 0.8670310974121094\n",
" 226 7.606559e+00 3.827245e-02\n",
" * time: 0.8704161643981934\n",
" 227 7.606450e+00 3.428008e-02\n",
" * time: 0.8734099864959717\n",
" 228 7.606253e+00 2.355543e-02\n",
" * time: 0.876230001449585\n",
" 229 7.606079e+00 1.852106e-02\n",
" * time: 0.8788330554962158\n",
" 230 7.606008e+00 3.101181e-02\n",
" * time: 0.8814911842346191\n",
" 231 7.605953e+00 2.070094e-02\n",
" * time: 0.8840770721435547\n",
" 232 7.605917e+00 2.917829e-02\n",
" * time: 0.8866641521453857\n",
" 233 7.605821e+00 3.039691e-02\n",
" * time: 0.8894500732421875\n",
" 234 7.605783e+00 2.392698e-02\n",
" * time: 0.8921670913696289\n",
" 235 7.605700e+00 2.815950e-02\n",
" * time: 0.8947710990905762\n",
" 236 7.605576e+00 1.935715e-02\n",
" * time: 0.9069721698760986\n",
" 237 7.605481e+00 1.931860e-02\n",
" * time: 0.9096810817718506\n",
" 238 7.605422e+00 1.531388e-02\n",
" * time: 0.9129509925842285\n",
" 239 7.605376e+00 1.544694e-02\n",
" * time: 0.9156031608581543\n",
" 240 7.605325e+00 1.717375e-02\n",
" * time: 0.9181911945343018\n",
" 241 7.605256e+00 1.002269e-02\n",
" * time: 0.9217259883880615\n",
" 242 7.605209e+00 1.544462e-02\n",
" * time: 0.9243340492248535\n",
" 243 7.605201e+00 2.207913e-02\n",
" * time: 0.926887035369873\n",
" 244 7.605157e+00 2.511728e-02\n",
" * time: 0.9295141696929932\n",
" 245 7.605137e+00 3.142990e-02\n",
" * time: 0.9321541786193848\n",
" 246 7.605048e+00 1.896229e-02\n",
" * time: 0.9347090721130371\n",
" 247 7.604929e+00 2.352782e-02\n",
" * time: 0.9374721050262451\n",
" 248 7.604840e+00 1.471586e-02\n",
" * time: 0.9406251907348633\n",
" 249 7.604778e+00 2.181153e-02\n",
" * time: 0.943162202835083\n",
" 250 7.604728e+00 2.156717e-02\n",
" * time: 0.945789098739624\n",
" 251 7.604681e+00 1.003909e-02\n",
" * time: 0.9490721225738525\n",
" 252 7.604616e+00 1.723120e-02\n",
" * time: 0.951786994934082\n",
" 253 7.604574e+00 1.797897e-02\n",
" * time: 0.9545040130615234\n",
" 254 7.604519e+00 1.787236e-02\n",
" * time: 0.9577701091766357\n",
" 255 7.604501e+00 2.624247e-02\n",
" * time: 0.9604151248931885\n",
" 256 7.604442e+00 1.489259e-02\n",
" * time: 0.9637119770050049\n",
" 257 7.604428e+00 1.925118e-02\n",
" * time: 0.9672281742095947\n",
" 258 7.604379e+00 1.614183e-02\n",
" * time: 0.9700720310211182\n",
" 259 7.604348e+00 1.747636e-02\n",
" * time: 0.9728331565856934\n",
" 260 7.604309e+00 2.199294e-02\n",
" * time: 0.9756121635437012\n",
" 261 7.604263e+00 1.450830e-02\n",
" * time: 0.9781861305236816\n",
" 262 7.604231e+00 1.182219e-02\n",
" * time: 0.981482982635498\n",
" 263 7.604179e+00 1.238737e-02\n",
" * time: 0.9843151569366455\n",
" 264 7.604136e+00 1.673160e-02\n",
" * time: 0.9869480133056641\n",
" 265 7.604099e+00 1.727775e-02\n",
" * time: 0.9895381927490234\n",
" 266 7.604054e+00 1.332617e-02\n",
" * time: 0.9927091598510742\n",
" 267 7.604032e+00 1.277240e-02\n",
" * time: 0.9952821731567383\n",
" 268 7.603983e+00 1.201359e-02\n",
" * time: 1.0072951316833496\n",
" 269 7.603955e+00 7.990877e-03\n",
" * time: 1.0105111598968506\n",
" 270 7.603943e+00 2.035873e-02\n",
" * time: 1.0129961967468262\n",
" 271 7.603909e+00 1.907264e-02\n",
" * time: 1.0164999961853027\n",
" 272 7.603893e+00 1.543917e-02\n",
" * time: 1.0190041065216064\n",
" 273 7.603865e+00 1.527497e-02\n",
" * time: 1.0216140747070312\n",
" 274 7.603835e+00 1.597506e-02\n",
" * time: 1.0241169929504395\n",
" 275 7.603826e+00 1.806827e-02\n",
" * time: 1.0266761779785156\n",
" 276 7.603779e+00 1.101678e-02\n",
" * time: 1.0292870998382568\n",
" 277 7.603747e+00 1.746900e-02\n",
" * time: 1.0320520401000977\n",
" 278 7.603698e+00 1.475117e-02\n",
" * time: 1.0352811813354492\n",
" 279 7.603667e+00 1.357222e-02\n",
" * time: 1.0384430885314941\n",
" 280 7.603623e+00 1.433077e-02\n",
" * time: 1.0410470962524414\n",
" 281 7.603589e+00 1.042037e-02\n",
" * time: 1.0443601608276367\n",
" 282 7.603570e+00 1.984943e-02\n",
" * time: 1.0471429824829102\n",
" 283 7.603517e+00 1.032937e-02\n",
" * time: 1.050469160079956\n",
" 284 7.603500e+00 1.234476e-02\n",
" * time: 1.0530810356140137\n",
" 285 7.603499e+00 2.251760e-02\n",
" * time: 1.0556550025939941\n",
" 286 7.603447e+00 1.354240e-02\n",
" * time: 1.05892014503479\n",
" 287 7.603420e+00 1.093271e-02\n",
" * time: 1.0616800785064697\n",
" 288 7.603385e+00 1.147930e-02\n",
" * time: 1.0642740726470947\n",
" 289 7.603346e+00 1.181847e-02\n",
" * time: 1.0673491954803467\n",
" 290 7.603314e+00 1.004018e-02\n",
" * time: 1.0705370903015137\n",
" 291 7.603300e+00 1.213438e-02\n",
" * time: 1.0730950832366943\n",
" 292 7.603283e+00 1.092828e-02\n",
" * time: 1.0765011310577393\n",
" 293 7.603248e+00 1.411630e-02\n",
" * time: 1.0792272090911865\n",
" 294 7.603247e+00 1.625853e-02\n",
" * time: 1.081840991973877\n",
" 295 7.603244e+00 2.115780e-02\n",
" * time: 1.0843970775604248\n",
" 296 7.603214e+00 1.537254e-02\n",
" * time: 1.087576150894165\n",
" 297 7.603179e+00 1.520047e-02\n",
" * time: 1.0901110172271729\n",
" 298 7.603149e+00 1.616339e-02\n",
" * time: 1.0928690433502197\n",
" 299 7.603134e+00 1.839369e-02\n",
" * time: 1.0954861640930176\n",
" 300 7.603108e+00 1.229016e-02\n",
" * time: 1.1067261695861816\n",
" 301 7.603089e+00 1.841628e-02\n",
" * time: 1.1095521450042725\n",
" 302 7.603057e+00 1.309476e-02\n",
" * time: 1.112725019454956\n",
" 303 7.603037e+00 1.137991e-02\n",
" * time: 1.1153161525726318\n",
" 304 7.603014e+00 8.580897e-03\n",
" * time: 1.1179661750793457\n",
" 305 7.603006e+00 1.426745e-02\n",
" * time: 1.1204440593719482\n",
" 306 7.602984e+00 1.107401e-02\n",
" * time: 1.1237890720367432\n",
" 307 7.602975e+00 1.145663e-02\n",
" * time: 1.126540184020996\n",
" 308 7.602960e+00 1.026057e-02\n",
" * time: 1.1297721862792969\n",
" 309 7.602948e+00 1.464455e-02\n",
" * time: 1.1323151588439941\n",
" 310 7.602929e+00 1.255633e-02\n",
" * time: 1.134916067123413\n",
" 311 7.602921e+00 1.987387e-02\n",
" * time: 1.1374061107635498\n",
" 312 7.602892e+00 1.225873e-02\n",
" * time: 1.1401610374450684\n",
" 313 7.602874e+00 1.389211e-02\n",
" * time: 1.142765998840332\n",
" 314 7.602848e+00 1.528448e-02\n",
" * time: 1.1454031467437744\n",
" 315 7.602822e+00 1.208840e-02\n",
" * time: 1.1486430168151855\n",
" 316 7.602794e+00 1.366943e-02\n",
" * time: 1.151794195175171\n",
" 317 7.602770e+00 9.392088e-03\n",
" * time: 1.1543970108032227\n",
" 318 7.602734e+00 7.288803e-03\n",
" * time: 1.1577951908111572\n",
" 319 7.602719e+00 7.570639e-03\n",
" * time: 1.1610100269317627\n",
" 320 7.602703e+00 4.957912e-03\n",
" * time: 1.1642131805419922\n",
" 321 7.602690e+00 1.247195e-02\n",
" * time: 1.1668531894683838\n",
" 322 7.602660e+00 1.204497e-02\n",
" * time: 1.1701161861419678\n",
" 323 7.602659e+00 1.607293e-02\n",
" * time: 1.1728360652923584\n",
" 324 7.602638e+00 1.730080e-02\n",
" * time: 1.175356149673462\n",
" 325 7.602604e+00 1.602050e-02\n",
" * time: 1.1779379844665527\n",
" 326 7.602575e+00 1.348129e-02\n",
" * time: 1.1812310218811035\n",
" 327 7.602546e+00 1.154319e-02\n",
" * time: 1.1844501495361328\n",
" 328 7.602537e+00 1.303211e-02\n",
" * time: 1.187195062637329\n",
" 329 7.602508e+00 1.061523e-02\n",
" * time: 1.1898200511932373\n",
" 330 7.602479e+00 1.445873e-02\n",
" * time: 1.1923420429229736\n",
" 331 7.602429e+00 9.369235e-03\n",
" * time: 1.1948561668395996\n",
" 332 7.602404e+00 8.101842e-03\n",
" * time: 1.2064909934997559\n",
" 333 7.602375e+00 6.271999e-03\n",
" * time: 1.2091021537780762\n",
" 334 7.602357e+00 4.978418e-03\n",
" * time: 1.2122361660003662\n",
" 335 7.602352e+00 1.431695e-02\n",
" * time: 1.214827060699463\n",
" 336 7.602350e+00 1.657193e-02\n",
" * time: 1.2174901962280273\n",
" 337 7.602335e+00 1.431146e-02\n",
" * time: 1.2201590538024902\n",
" 338 7.602312e+00 9.876441e-03\n",
" * time: 1.2232930660247803\n",
" 339 7.602307e+00 1.229028e-02\n",
" * time: 1.2259571552276611\n",
" 340 7.602281e+00 9.363639e-03\n",
" * time: 1.2284810543060303\n",
" 341 7.602257e+00 8.939221e-03\n",
" * time: 1.2309889793395996\n",
" 342 7.602247e+00 1.091453e-02\n",
" * time: 1.2338571548461914\n",
" 343 7.602240e+00 1.052333e-02\n",
" * time: 1.2364501953125\n",
" 344 7.602233e+00 1.087575e-02\n",
" * time: 1.23905611038208\n",
" 345 7.602232e+00 8.966033e-03\n",
" * time: 1.241657018661499\n",
" 346 7.602231e+00 1.402991e-02\n",
" * time: 1.2442920207977295\n",
" 347 7.602230e+00 1.536465e-02\n",
" * time: 1.2468061447143555\n",
" 348 7.602216e+00 9.178715e-03\n",
" * time: 1.2495701313018799\n",
" 349 7.602209e+00 9.114835e-03\n",
" * time: 1.2521581649780273\n",
" 350 7.602195e+00 9.065201e-03\n",
" * time: 1.254676103591919\n",
" 351 7.602188e+00 7.022926e-03\n",
" * time: 1.2572300434112549\n",
" 352 7.602175e+00 6.419884e-03\n",
" * time: 1.2597401142120361\n",
" 353 7.602172e+00 9.807445e-03\n",
" * time: 1.2622971534729004\n",
" 354 7.602161e+00 7.302537e-03\n",
" * time: 1.2658610343933105\n",
" 355 7.602153e+00 6.264240e-03\n",
" * time: 1.2691781520843506\n",
" 356 7.602146e+00 5.755593e-03\n",
" * time: 1.2725701332092285\n",
" 357 7.602140e+00 7.100578e-03\n",
" * time: 1.2758691310882568\n",
" 358 7.602132e+00 4.822609e-03\n",
" * time: 1.2790629863739014\n",
" 359 7.602132e+00 9.622222e-03\n",
" * time: 1.2818419933319092\n",
" 360 7.602130e+00 7.285010e-03\n",
" * time: 1.2843680381774902\n",
" 361 7.602126e+00 6.651959e-03\n",
" * time: 1.2868680953979492\n",
" 362 7.602118e+00 5.271510e-03\n",
" * time: 1.2894160747528076\n",
" 363 7.602111e+00 7.457443e-03\n",
" * time: 1.291917085647583\n",
" 364 7.602104e+00 6.304689e-03\n",
" * time: 1.2947111129760742\n",
" 365 7.602103e+00 6.531506e-03\n",
" * time: 1.306042194366455\n",
" 366 7.602096e+00 5.755674e-03\n",
" * time: 1.3087470531463623\n",
" 367 7.602085e+00 8.112194e-03\n",
" * time: 1.3114161491394043\n",
" 368 7.602078e+00 1.030216e-02\n",
" * time: 1.3144311904907227\n",
" 369 7.602074e+00 9.635798e-03\n",
" * time: 1.317397117614746\n",
" 370 7.602064e+00 6.103323e-03\n",
" * time: 1.3209140300750732\n",
" 371 7.602051e+00 7.238742e-03\n",
" * time: 1.3243930339813232\n",
" 372 7.602039e+00 7.585756e-03\n",
" * time: 1.3281550407409668\n",
" 373 7.602038e+00 1.031615e-02\n",
" * time: 1.330937147140503\n",
" 374 7.602024e+00 9.093322e-03\n",
" * time: 1.333923101425171\n",
" 375 7.602012e+00 6.186042e-03\n",
" * time: 1.3377130031585693\n",
" 376 7.602000e+00 1.205058e-02\n",
" * time: 1.3405981063842773\n",
" 377 7.601989e+00 5.999396e-03\n",
" * time: 1.3444840908050537\n",
" 378 7.601988e+00 1.435872e-02\n",
" * time: 1.347398042678833\n",
" 379 7.601976e+00 1.021123e-02\n",
" * time: 1.3510401248931885\n",
" 380 7.601967e+00 8.358129e-03\n",
" * time: 1.354018211364746\n",
" 381 7.601959e+00 9.493374e-03\n",
" * time: 1.356909990310669\n",
" 382 7.601952e+00 9.069549e-03\n",
" * time: 1.3598790168762207\n",
" 383 7.601939e+00 6.853795e-03\n",
" * time: 1.362732172012329\n",
" 384 7.601926e+00 1.062256e-02\n",
" * time: 1.3656129837036133\n",
" 385 7.601915e+00 9.519083e-03\n",
" * time: 1.368636131286621\n",
" 386 7.601910e+00 1.077891e-02\n",
" * time: 1.3716461658477783\n",
" 387 7.601903e+00 1.110914e-02\n",
" * time: 1.374634027481079\n",
" 388 7.601892e+00 1.021958e-02\n",
" * time: 1.3775811195373535\n",
" 389 7.601874e+00 7.029073e-03\n",
" * time: 1.3811171054840088\n",
" 390 7.601866e+00 8.121670e-03\n",
" * time: 1.3839471340179443\n",
" 391 7.601850e+00 7.658265e-03\n",
" * time: 1.3868680000305176\n",
" 392 7.601835e+00 5.896006e-03\n",
" * time: 1.3905301094055176\n",
" 393 7.601825e+00 3.696215e-03\n",
" * time: 1.3941330909729004\n",
" 394 7.601818e+00 5.447625e-03\n",
" * time: 1.3977761268615723\n",
" 395 7.601814e+00 4.229950e-03\n",
" * time: 1.4016261100769043\n",
" 396 7.601811e+00 5.424498e-03\n",
" * time: 1.4047160148620605\n",
" 397 7.601806e+00 5.343858e-03\n",
" * time: 1.4175431728363037\n",
" 398 7.601800e+00 2.822826e-03\n",
" * time: 1.4213612079620361\n",
" 399 7.601798e+00 3.404277e-03\n",
" * time: 1.424375057220459\n",
" 400 7.601794e+00 4.596025e-03\n",
" * time: 1.4273450374603271\n",
" 401 7.601793e+00 5.817366e-03\n",
" * time: 1.4302830696105957\n",
" 402 7.601793e+00 5.660567e-03\n",
" * time: 1.4332501888275146\n",
" 403 7.601790e+00 5.711754e-03\n",
" * time: 1.436156988143921\n",
" 404 7.601786e+00 4.330384e-03\n",
" * time: 1.439073085784912\n",
" 405 7.601783e+00 6.533761e-03\n",
" * time: 1.441946029663086\n",
" 406 7.601783e+00 6.520422e-03\n",
" * time: 1.4449450969696045\n",
" 407 7.601779e+00 3.552697e-03\n",
" * time: 1.448700189590454\n",
" 408 7.601777e+00 4.317950e-03\n",
" * time: 1.4515511989593506\n",
" 409 7.601773e+00 3.658759e-03\n",
" * time: 1.4545130729675293\n",
" 410 7.601773e+00 4.169959e-03\n",
" * time: 1.4573512077331543\n",
" 411 7.601770e+00 3.541989e-03\n",
" * time: 1.460245132446289\n",
" 412 7.601767e+00 2.919978e-03\n",
" * time: 1.4638621807098389\n",
" 413 7.601765e+00 2.306088e-03\n",
" * time: 1.4675991535186768\n",
" 414 7.601765e+00 4.603876e-03\n",
" * time: 1.4706101417541504\n",
" 415 7.601764e+00 2.592792e-03\n",
" * time: 1.4734821319580078\n",
" 416 7.601764e+00 5.820079e-03\n",
" * time: 1.4763050079345703\n",
" 417 7.601762e+00 3.656667e-03\n",
" * time: 1.479858160018921\n",
" 418 7.601761e+00 3.741023e-03\n",
" * time: 1.4827890396118164\n",
" 419 7.601759e+00 2.929388e-03\n",
" * time: 1.4866089820861816\n",
" 420 7.601758e+00 4.013142e-03\n",
" * time: 1.4897840023040771\n",
" 421 7.601758e+00 2.793870e-03\n",
" * time: 1.4927890300750732\n",
" 422 7.601757e+00 2.360518e-03\n",
" * time: 1.495711088180542\n",
" 423 7.601756e+00 2.365338e-03\n",
" * time: 1.4986560344696045\n",
" 424 7.601755e+00 3.192826e-03\n",
" * time: 1.5016820430755615\n",
" 425 7.601755e+00 2.412923e-03\n",
" * time: 1.5052530765533447\n",
" 426 7.601754e+00 3.148039e-03\n",
" * time: 1.50816011428833\n",
" 427 7.601754e+00 2.012344e-03\n",
" * time: 1.511864185333252\n",
" 428 7.601753e+00 1.856823e-03\n",
" * time: 1.5147860050201416\n",
" 429 7.601753e+00 2.644970e-03\n",
" * time: 1.5263500213623047\n",
" 430 7.601752e+00 2.145633e-03\n",
" * time: 1.5293211936950684\n",
" 431 7.601752e+00 2.197014e-03\n",
" * time: 1.5322060585021973\n",
" 432 7.601751e+00 2.249695e-03\n",
" * time: 1.535818099975586\n",
" 433 7.601751e+00 1.901828e-03\n",
" * time: 1.538712978363037\n",
" 434 7.601751e+00 2.038875e-03\n",
" * time: 1.5418040752410889\n",
" 435 7.601750e+00 1.804242e-03\n",
" * time: 1.5455701351165771\n",
" 436 7.601750e+00 1.348346e-03\n",
" * time: 1.548511028289795\n",
" 437 7.601749e+00 1.335443e-03\n",
" * time: 1.5521981716156006\n",
" 438 7.601749e+00 1.622453e-03\n",
" * time: 1.5551090240478516\n",
" 439 7.601749e+00 1.555947e-03\n",
" * time: 1.5581719875335693\n",
" 440 7.601748e+00 1.622047e-03\n",
" * time: 1.56193208694458\n",
" 441 7.601748e+00 3.275000e-03\n",
" * time: 1.564866065979004\n",
" 442 7.601747e+00 2.398138e-03\n",
" * time: 1.5684731006622314\n",
" 443 7.601747e+00 1.801811e-03\n",
" * time: 1.5713770389556885\n",
" 444 7.601746e+00 2.165522e-03\n",
" * time: 1.5744051933288574\n",
" 445 7.601746e+00 1.465532e-03\n",
" * time: 1.5774500370025635\n",
" 446 7.601745e+00 1.569624e-03\n",
" * time: 1.580368995666504\n",
" 447 7.601745e+00 2.723922e-03\n",
" * time: 1.5832490921020508\n",
" 448 7.601744e+00 2.231168e-03\n",
" * time: 1.5861999988555908\n",
" 449 7.601744e+00 1.731660e-03\n",
" * time: 1.5898940563201904\n",
" 450 7.601744e+00 2.431292e-03\n",
" * time: 1.5928571224212646\n",
" 451 7.601743e+00 1.892679e-03\n",
" * time: 1.5957520008087158\n",
" 452 7.601742e+00 1.331262e-03\n",
" * time: 1.5994021892547607\n",
" 453 7.601742e+00 1.163571e-03\n",
" * time: 1.6030900478363037\n",
" 454 7.601742e+00 1.453840e-03\n",
" * time: 1.6060760021209717\n",
" 455 7.601742e+00 3.181580e-03\n",
" * time: 1.609076976776123\n",
" 456 7.601742e+00 2.186529e-03\n",
" * time: 1.6119351387023926\n",
" 457 7.601741e+00 2.224637e-03\n",
" * time: 1.6148111820220947\n",
" 458 7.601740e+00 9.701961e-04\n",
" * time: 1.6184840202331543\n",
" 459 7.601738e+00 1.012902e-03\n",
" * time: 1.6215410232543945\n",
" 460 7.601738e+00 2.499266e-03\n",
" * time: 1.6246321201324463\n",
" 461 7.601738e+00 1.716088e-03\n",
" * time: 1.6283621788024902\n",
" 462 7.601737e+00 1.466418e-03\n",
" * time: 1.641982078552246\n",
" 463 7.601737e+00 1.983545e-03\n",
" * time: 1.6451749801635742\n",
" 464 7.601736e+00 1.854853e-03\n",
" * time: 1.6483571529388428\n",
" 465 7.601736e+00 2.178618e-03\n",
" * time: 1.6515769958496094\n",
" 466 7.601735e+00 1.857955e-03\n",
" * time: 1.6555399894714355\n",
" 467 7.601734e+00 1.228932e-03\n",
" * time: 1.65938401222229\n",
" 468 7.601733e+00 2.009172e-03\n",
" * time: 1.6625690460205078\n",
" 469 7.601733e+00 1.551899e-03\n",
" * time: 1.665701150894165\n",
" 470 7.601733e+00 1.563505e-03\n",
" * time: 1.6688640117645264\n",
" 471 7.601732e+00 8.060359e-04\n",
" * time: 1.672922134399414\n",
" 472 7.601732e+00 2.423650e-03\n",
" * time: 1.675995111465454\n",
" 473 7.601731e+00 1.352108e-03\n",
" * time: 1.679063081741333\n",
" 474 7.601731e+00 1.177188e-03\n",
" * time: 1.6821691989898682\n",
" 475 7.601731e+00 1.393213e-03\n",
" * time: 1.6852049827575684\n",
" 476 7.601730e+00 1.647053e-03\n",
" * time: 1.6882750988006592\n",
" 477 7.601730e+00 1.274533e-03\n",
" * time: 1.6911721229553223\n",
" 478 7.601730e+00 1.798294e-03\n",
" * time: 1.6941301822662354\n",
" 479 7.601729e+00 1.183953e-03\n",
" * time: 1.698038101196289\n",
" 480 7.601729e+00 9.083597e-04\n",
" * time: 1.7019670009613037\n",
" 481 7.601729e+00 1.291266e-03\n",
" * time: 1.705254077911377\n",
" 482 7.601728e+00 1.306728e-03\n",
" * time: 1.708724021911621\n",
" 483 7.601728e+00 2.311934e-03\n",
" * time: 1.713150978088379\n",
" 484 7.601728e+00 1.605010e-03\n",
" * time: 1.7189760208129883\n",
" 485 7.601728e+00 1.830582e-03\n",
" * time: 1.7229721546173096\n",
" 486 7.601727e+00 2.085009e-03\n",
" * time: 1.7288191318511963\n",
" 487 7.601727e+00 2.599085e-03\n",
" * time: 1.7322590351104736\n",
" 488 7.601727e+00 3.076085e-03\n",
" * time: 1.735727071762085\n",
" 489 7.601726e+00 1.539651e-03\n",
" * time: 1.740394115447998\n",
" 490 7.601725e+00 2.111897e-03\n",
" * time: 1.7435951232910156\n",
" 491 7.601724e+00 1.474152e-03\n",
" * time: 1.7466180324554443\n",
" 492 7.601724e+00 2.630003e-03\n",
" * time: 1.749725103378296\n",
" 493 7.601724e+00 2.044035e-03\n",
" * time: 1.7528769969940186\n",
" 494 7.601723e+00 1.217518e-03\n",
" * time: 1.7667920589447021\n",
" 495 7.601723e+00 1.349260e-03\n",
" * time: 1.769975185394287\n",
" 496 7.601722e+00 1.986227e-03\n",
" * time: 1.7729251384735107\n",
" 497 7.601722e+00 2.142998e-03\n",
" * time: 1.7759201526641846\n",
" 498 7.601722e+00 1.561417e-03\n",
" * time: 1.779520034790039\n",
" 499 7.601721e+00 1.344068e-03\n",
" * time: 1.7832400798797607\n",
" 500 7.601721e+00 1.850573e-03\n",
" * time: 1.7860422134399414\n",
" 501 7.601720e+00 1.757191e-03\n",
" * time: 1.78895902633667\n",
" 502 7.601720e+00 1.479737e-03\n",
" * time: 1.7919230461120605\n",
" 503 7.601720e+00 2.041182e-03\n",
" * time: 1.7948870658874512\n",
" 504 7.601719e+00 1.625627e-03\n",
" * time: 1.7979090213775635\n",
" 505 7.601719e+00 1.672737e-03\n",
" * time: 1.8007540702819824\n",
" 506 7.601719e+00 2.710896e-03\n",
" * time: 1.8035860061645508\n",
" 507 7.601718e+00 1.636767e-03\n",
" * time: 1.8074150085449219\n",
" 508 7.601718e+00 1.227500e-03\n",
" * time: 1.8111541271209717\n",
" 509 7.601718e+00 1.591151e-03\n",
" * time: 1.8143761157989502\n",
" 510 7.601717e+00 2.011660e-03\n",
" * time: 1.8174312114715576\n",
" 511 7.601717e+00 2.628653e-03\n",
" * time: 1.8206801414489746\n",
" 512 7.601717e+00 1.219903e-03\n",
" * time: 1.8245580196380615\n",
" 513 7.601717e+00 1.890937e-03\n",
" * time: 1.8275961875915527\n",
" 514 7.601716e+00 2.578137e-03\n",
" * time: 1.830653190612793\n",
" 515 7.601716e+00 1.643552e-03\n",
" * time: 1.8337130546569824\n",
" 516 7.601716e+00 2.750010e-03\n",
" * time: 1.8367671966552734\n",
" 517 7.601715e+00 1.421810e-03\n",
" * time: 1.8406810760498047\n",
" 518 7.601715e+00 1.283363e-03\n",
" * time: 1.843731164932251\n",
" 519 7.601714e+00 1.581835e-03\n",
" * time: 1.8468279838562012\n",
" 520 7.601714e+00 2.081754e-03\n",
" * time: 1.8499090671539307\n",
" 521 7.601714e+00 1.322130e-03\n",
" * time: 1.853734016418457\n",
" 522 7.601713e+00 1.413080e-03\n",
" * time: 1.8575329780578613\n",
" 523 7.601713e+00 1.448607e-03\n",
" * time: 1.8605492115020752\n",
" 524 7.601713e+00 1.958482e-03\n",
" * time: 1.8634309768676758\n",
" 525 7.601713e+00 1.315426e-03\n",
" * time: 1.8669641017913818\n",
" 526 7.601713e+00 1.624884e-03\n",
" * time: 1.8789701461791992\n",
" 527 7.601713e+00 1.421138e-03\n",
" * time: 1.8822510242462158\n",
" 528 7.601713e+00 8.507989e-04\n",
" * time: 1.8854119777679443\n",
" 529 7.601712e+00 1.785049e-03\n",
" * time: 1.8885231018066406\n",
" 530 7.601712e+00 1.547410e-03\n",
" * time: 1.8918161392211914\n",
" 531 7.601712e+00 1.186840e-03\n",
" * time: 1.8949711322784424\n",
" 532 7.601712e+00 9.836154e-04\n",
" * time: 1.898186206817627\n",
" 533 7.601712e+00 1.116964e-03\n",
" * time: 1.901252031326294\n",
" 534 7.601711e+00 1.137647e-03\n",
" * time: 1.905271053314209\n",
" 535 7.601711e+00 9.823750e-04\n",
" * time: 1.9093341827392578\n",
" 536 7.601711e+00 1.103677e-03\n",
" * time: 1.912431001663208\n",
" 537 7.601711e+00 1.453290e-03\n",
" * time: 1.9155771732330322\n",
" 538 7.601711e+00 1.030640e-03\n",
" * time: 1.9193542003631592\n",
" 539 7.601711e+00 1.030530e-03\n",
" * time: 1.9225780963897705\n",
" 540 7.601711e+00 1.301679e-03\n",
" * time: 1.9255690574645996\n",
" 541 7.601710e+00 7.883370e-04\n",
" * time: 1.9293510913848877\n",
" 542 7.601710e+00 1.146433e-03\n",
" * time: 1.9323761463165283\n",
" 543 7.601710e+00 1.488085e-03\n",
" * time: 1.9353580474853516\n",
" 544 7.601710e+00 1.196471e-03\n",
" * time: 1.9383680820465088\n",
" 545 7.601710e+00 1.343770e-03\n",
" * time: 1.941390037536621\n",
" 546 7.601710e+00 1.700005e-03\n",
" * time: 1.9443700313568115\n",
" 547 7.601710e+00 1.291226e-03\n",
" * time: 1.9482860565185547\n",
" 548 7.601709e+00 7.129353e-04\n",
" * time: 1.9519951343536377\n",
" 549 7.601709e+00 1.020169e-03\n",
" * time: 1.9549610614776611\n",
" 550 7.601709e+00 9.157539e-04\n",
" * time: 1.9578571319580078\n",
" 551 7.601709e+00 5.296825e-04\n",
" * time: 1.9615252017974854\n",
" 552 7.601709e+00 8.500480e-04\n",
" * time: 1.9646430015563965\n",
" 553 7.601708e+00 8.133626e-04\n",
" * time: 1.9684550762176514\n",
" 554 7.601708e+00 9.901256e-04\n",
" * time: 1.9723341464996338\n",
" 555 7.601708e+00 1.382346e-03\n",
" * time: 1.9752631187438965\n",
" 556 7.601708e+00 1.158238e-03\n",
" * time: 1.9784011840820312\n",
" 557 7.601708e+00 1.885929e-03\n",
" * time: 1.9814660549163818\n",
" 558 7.601707e+00 1.586508e-03\n",
" * time: 1.9845399856567383\n",
" 559 7.601707e+00 1.515507e-03\n",
" * time: 1.9972620010375977\n",
" 560 7.601706e+00 8.078577e-04\n",
" * time: 2.0013680458068848\n",
" 561 7.601706e+00 1.001990e-03\n",
" * time: 2.0047531127929688\n",
" 562 7.601706e+00 8.909407e-04\n",
" * time: 2.0079591274261475\n",
" 563 7.601706e+00 1.116623e-03\n",
" * time: 2.0110480785369873\n",
" 564 7.601706e+00 1.309177e-03\n",
" * time: 2.014341115951538\n",
" 565 7.601706e+00 1.434728e-03\n",
" * time: 2.01739501953125\n",
" 566 7.601705e+00 6.304497e-04\n",
" * time: 2.0212860107421875\n",
" 567 7.601705e+00 1.514880e-03\n",
" * time: 2.024260997772217\n",
" 568 7.601705e+00 1.817257e-03\n",
" * time: 2.0273079872131348\n",
" 569 7.601705e+00 8.884819e-04\n",
" * time: 2.031238079071045\n",
" 570 7.601705e+00 1.197027e-03\n",
" * time: 2.034196138381958\n",
" 571 7.601705e+00 1.294803e-03\n",
" * time: 2.0372021198272705\n",
" 572 7.601705e+00 1.556670e-03\n",
" * time: 2.0402181148529053\n",
" 573 7.601705e+00 9.812680e-04\n",
" * time: 2.0431342124938965\n",
" 574 7.601705e+00 1.084959e-03\n",
" * time: 2.0462541580200195\n",
" 575 7.601704e+00 7.528757e-04\n",
" * time: 2.0500872135162354\n",
" 576 7.601704e+00 1.414651e-03\n",
" * time: 2.0531210899353027\n",
" 577 7.601704e+00 1.069341e-03\n",
" * time: 2.0560760498046875\n",
" 578 7.601704e+00 9.009442e-04\n",
" * time: 2.059882164001465\n",
" 579 7.601704e+00 4.919636e-04\n",
" * time: 2.0636730194091797\n",
" 580 7.601704e+00 7.756459e-04\n",
" * time: 2.066587209701538\n",
" 581 7.601703e+00 5.640505e-04\n",
" * time: 2.07027006149292\n",
" 582 7.601703e+00 8.453276e-04\n",
" * time: 2.073103189468384\n",
" 583 7.601703e+00 7.959033e-04\n",
" * time: 2.076719045639038\n",
" 584 7.601703e+00 8.019355e-04\n",
" * time: 2.079759120941162\n",
" 585 7.601703e+00 6.254290e-04\n",
" * time: 2.083583116531372\n",
" 586 7.601703e+00 9.470776e-04\n",
" * time: 2.086496114730835\n",
" 587 7.601703e+00 5.079227e-04\n",
" * time: 2.0901410579681396\n",
" 588 7.601703e+00 6.435347e-04\n",
" * time: 2.0931200981140137\n",
" 589 7.601703e+00 5.389022e-04\n",
" * time: 2.096198081970215\n",
" 590 7.601703e+00 6.847118e-04\n",
" * time: 2.099159002304077\n",
" 591 7.601703e+00 3.683523e-04\n",
" * time: 2.113013982772827\n",
" 592 7.601702e+00 8.441647e-04\n",
" * time: 2.116180181503296\n",
" 593 7.601702e+00 8.817415e-04\n",
" * time: 2.119197130203247\n",
" 594 7.601702e+00 6.288657e-04\n",
" * time: 2.122380018234253\n",
" 595 7.601702e+00 4.039963e-04\n",
" * time: 2.126444101333618\n",
" 596 7.601702e+00 3.598258e-04\n",
" * time: 2.129556179046631\n",
" 597 7.601702e+00 4.636102e-04\n",
" * time: 2.1325831413269043\n",
" 598 7.601702e+00 4.791832e-04\n",
" * time: 2.1355772018432617\n",
" 599 7.601702e+00 6.146933e-04\n",
" * time: 2.1386032104492188\n",
" 600 7.601702e+00 4.682963e-04\n",
" * time: 2.142468214035034\n",
" 601 7.601702e+00 5.294664e-04\n",
" * time: 2.14554500579834\n",
" 602 7.601702e+00 8.116131e-04\n",
" * time: 2.1486032009124756\n",
" 603 7.601702e+00 1.174954e-03\n",
" * time: 2.1516151428222656\n",
" 604 7.601702e+00 6.426663e-04\n",
" * time: 2.155407190322876\n",
" 605 7.601702e+00 1.306553e-03\n",
" * time: 2.1583449840545654\n",
" 606 7.601702e+00 7.561872e-04\n",
" * time: 2.1612861156463623\n",
" 607 7.601701e+00 9.179627e-04\n",
" * time: 2.1642889976501465\n",
" 608 7.601701e+00 1.105930e-03\n",
" * time: 2.1671900749206543\n",
" 609 7.601701e+00 7.836292e-04\n",
" * time: 2.170220136642456\n",
" 610 7.601701e+00 5.323647e-04\n",
" * time: 2.1740009784698486\n",
" 611 7.601701e+00 9.100983e-04\n",
" * time: 2.1769490242004395\n",
" 612 7.601701e+00 6.745528e-04\n",
" * time: 2.1798200607299805\n",
" 613 7.601701e+00 1.057477e-03\n",
" * time: 2.1827712059020996\n",
" 614 7.601701e+00 4.638895e-04\n",
" * time: 2.186608076095581\n",
" 615 7.601701e+00 7.871339e-04\n",
" * time: 2.18959903717041\n",
" 616 7.601701e+00 5.060631e-04\n",
" * time: 2.193356990814209\n",
" 617 7.601701e+00 6.418214e-04\n",
" * time: 2.196364164352417\n",
" 618 7.601700e+00 7.404004e-04\n",
" * time: 2.1993601322174072\n",
" 619 7.601700e+00 4.495868e-04\n",
" * time: 2.203050136566162\n",
" 620 7.601700e+00 4.744222e-04\n",
" * time: 2.2059080600738525\n",
" 621 7.601700e+00 6.793283e-04\n",
" * time: 2.208829164505005\n",
" 622 7.601700e+00 7.381390e-04\n",
" * time: 2.2117021083831787\n",
" 623 7.601700e+00 7.786096e-04\n",
" * time: 2.214726209640503\n",
" 624 7.601700e+00 6.792298e-04\n",
" * time: 2.227195978164673\n",
" 625 7.601700e+00 5.415716e-04\n",
" * time: 2.2304940223693848\n",
" 626 7.601700e+00 6.761057e-04\n",
" * time: 2.23362398147583\n",
" 627 7.601700e+00 6.021939e-04\n",
" * time: 2.2368390560150146\n",
" 628 7.601700e+00 3.713081e-04\n",
" * time: 2.2405600547790527\n",
" 629 7.601700e+00 2.312683e-04\n",
" * time: 2.244556188583374\n",
" 630 7.601700e+00 3.595985e-04\n",
" * time: 2.2474091053009033\n",
" 631 7.601700e+00 6.287301e-04\n",
" * time: 2.250326156616211\n",
" 632 7.601700e+00 5.789522e-04\n",
" * time: 2.25337815284729\n",
" 633 7.601700e+00 3.759291e-04\n",
" * time: 2.257052183151245\n",
" 634 7.601700e+00 5.982972e-04\n",
" * time: 2.2601230144500732\n",
" 635 7.601700e+00 4.256669e-04\n",
" * time: 2.2637062072753906\n",
" 636 7.601700e+00 4.665285e-04\n",
" * time: 2.266757011413574\n",
" 637 7.601700e+00 5.853447e-04\n",
" * time: 2.2697842121124268\n",
" 638 7.601700e+00 4.959765e-04\n",
" * time: 2.272855043411255\n",
" 639 7.601700e+00 3.967438e-04\n",
" * time: 2.275834083557129\n",
" 640 7.601700e+00 4.022015e-04\n",
" * time: 2.278771162033081\n",
" 641 7.601700e+00 6.005547e-04\n",
" * time: 2.281776189804077\n",
" 642 7.601700e+00 7.672319e-04\n",
" * time: 2.2847371101379395\n",
" 643 7.601700e+00 5.603547e-04\n",
" * time: 2.288323163986206\n",
" 644 7.601700e+00 6.216656e-04\n",
" * time: 2.291114091873169\n",
" 645 7.601700e+00 6.833152e-04\n",
" * time: 2.2939000129699707\n",
" 646 7.601700e+00 6.028079e-04\n",
" * time: 2.2968060970306396\n",
" 647 7.601700e+00 4.794747e-04\n",
" * time: 2.299669027328491\n",
" 648 7.601699e+00 4.754046e-04\n",
" * time: 2.3033649921417236\n",
" 649 7.601699e+00 4.903187e-04\n",
" * time: 2.3071470260620117\n",
" 650 7.601699e+00 4.276299e-04\n",
" * time: 2.310965061187744\n",
" 651 7.601699e+00 6.280609e-04\n",
" * time: 2.3139162063598633\n",
" 652 7.601699e+00 4.767447e-04\n",
" * time: 2.317676067352295\n",
" 653 7.601699e+00 2.351078e-04\n",
" * time: 2.3216052055358887\n",
" 654 7.601699e+00 3.490193e-04\n",
" * time: 2.3246982097625732\n",
" 655 7.601699e+00 3.438115e-04\n",
" * time: 2.3277502059936523\n",
" 656 7.601699e+00 4.664146e-04\n",
" * time: 2.340193033218384\n",
" 657 7.601699e+00 3.725306e-04\n",
" * time: 2.3434181213378906\n",
" 658 7.601699e+00 4.084074e-04\n",
" * time: 2.346482992172241\n",
" 659 7.601699e+00 3.346699e-04\n",
" * time: 2.349614143371582\n",
" 660 7.601699e+00 3.383938e-04\n",
" * time: 2.3525919914245605\n",
" 661 7.601699e+00 3.068648e-04\n",
" * time: 2.356295108795166\n",
" 662 7.601699e+00 1.728582e-04\n",
" * time: 2.3601272106170654\n",
" 663 7.601699e+00 3.387437e-04\n",
" * time: 2.3630990982055664\n",
" 664 7.601699e+00 4.019420e-04\n",
" * time: 2.3660600185394287\n",
" 665 7.601699e+00 2.850965e-04\n",
" * time: 2.3696770668029785\n",
" 666 7.601699e+00 3.366773e-04\n",
" * time: 2.372636079788208\n",
" 667 7.601699e+00 3.801911e-04\n",
" * time: 2.376376152038574\n",
" 668 7.601699e+00 3.906726e-04\n",
" * time: 2.3791861534118652\n",
" 669 7.601699e+00 2.108485e-04\n",
" * time: 2.382791042327881\n",
" 670 7.601699e+00 3.180458e-04\n",
" * time: 2.3857262134552\n",
" 671 7.601699e+00 2.666382e-04\n",
" * time: 2.388669013977051\n",
" 672 7.601699e+00 1.899630e-04\n",
" * time: 2.391676187515259\n",
" 673 7.601699e+00 3.278806e-04\n",
" * time: 2.3945751190185547\n",
" 674 7.601699e+00 3.942458e-04\n",
" * time: 2.397562026977539\n",
" 675 7.601699e+00 3.937956e-04\n",
" * time: 2.400435209274292\n",
" 676 7.601699e+00 6.402291e-04\n",
" * time: 2.403179168701172\n",
" 677 7.601699e+00 3.456600e-04\n",
" * time: 2.4069151878356934\n",
" 678 7.601699e+00 3.334172e-04\n",
" * time: 2.409675121307373\n",
" 679 7.601699e+00 2.773409e-04\n",
" * time: 2.412428140640259\n",
" 680 7.601699e+00 4.005457e-04\n",
" * time: 2.4153330326080322\n",
" 681 7.601699e+00 2.724667e-04\n",
" * time: 2.4190151691436768\n",
" 682 7.601699e+00 3.768761e-04\n",
" * time: 2.421985149383545\n",
" 683 7.601699e+00 2.822033e-04\n",
" * time: 2.424889087677002\n",
" 684 7.601699e+00 1.898422e-04\n",
" * time: 2.4279141426086426\n",
" 685 7.601699e+00 1.904448e-04\n",
" * time: 2.4309589862823486\n"
]
},
{
"output_type": "execute_result",
"data": {
"text/plain": "Plot{Plots.GRBackend() n=1}",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3dZ3wU1f4G8DOb3Wx6Dyl0CCUEFEgBqSFEBWNCkR6RYsGCFQv6R65gQUEUr6KCgsBFQlOBeEUIxSsqJSD9UhIIAdJ7z/b/i9G9Ec4vMOuEYbPP98OL5OzZyezshpMz88zvCBaLhQEAADgqldI7AAAAoCQMhAAA4NAwEAIAgEPDQAgAAA4NAyEAADg0DIQAAODQMBACAIBDw0AIAAAODQMhAAA4NAyEAADg0NR/8/nnMy94BQS6urrKsje3IbPZLAiCIAhK78itYLFYLBaLSuUofx6ZzWa82GYJL7aJeGs1N+xz7lxpXZ1R0marqytbtnRq3769rfv1dwl/s9Zo8mNPPjDrdbn2BgAAbluju4TcsE/PnquPHy+SuOFTI0aUb9myxba9+vsc5Y8mAAAArr97ahQAAMBKEJjdXUrCjBAAABwaZoQAACAbG9KFFovCU0gMhAAAICucGgUAALAjmBECAIB87DAsg4EQAABkY0sFEqWvEeLUKAAAODQMhAAA4NBwahQAAGRjww31Sp8ZxUAIAADyseEaoeLhGpwaBQAAh4YZIQAAyEewvxvqMRACAIBsBOmnOnFqFAAAQEmYEQIAgGxsuaFe6SkhBkIAAJCPDdcIcWoUAABAQRgIAQBAToL0f43797//3b9//+7du8+ZM8doNDZ86OzZs4mJiV27dp08eXJRURFjbOfOnRMmTIiIiBgwYMDSpUstFgtj7OrVq+PGjevatesDDzxw6dIl8bm//vprXFxceHi4ijGWk5Pz3HPP3X333cnJyYcPH2aMmUymDz/8cNiwYQ8++OCpU6dkP0wAANAsCcIflwklaHQozMzMnDhx4ksvvbR58+bt27d/8MEH1odMJlNCQkJ0dPT27dvVavW0adMYYwcOHLj33ns3btw4d+7ct99++4svvmCMjR8/PigoaPv27R07dhw1ahRjrLS0NCEhYfLkyampqer8/Py+ffuOGTPmpZdeKi4urq6uZowtWrRo3bp1S5YsSU9PHzJkSEZGho+PT9MePAAAgOusWLEiKSlp5MiRjLH58+fPnDnz5ZdfFh/asWOHwWB4/fXXBUFYvHhxSEhIdnb23LlzxUcjIiKmTp26Z8+evn37Hjt2bNeuXa6uru+8886KFSv2799/+PDh3r17i2On+q233ho8ePCHH35o/akmk+njjz9evXp1XFxcXFzcjh071q5dO3PmzFv+8gEAwN7IHZY5depUXFyc+HVUVNSlS5eqq6s9PDzEh6KiosSQqp+fX7t27U6fPt22bVvrc9PT0wcOHHjq1KmIiAhXV1fGmFqtvvPOO0+ePHnq1KnIyEixm+qXX34ZPHjw3LlzH3744ZSUFMZYQUFBbm5uv379xB79+vU7cuSIxJcFAACOyIYzo43fblFYWGg9JSl+UVhYKH5bVFTk7e1t7enr62t9iDG2dOnSCxcuPP/889xuRUVF1s2qLl++PH/+/NDQ0HvvvXfOnDkLFy4sKChwcXFxc3MTe/j7+xcUFFC7WFdXJ+kYAQCAnTKZTDfsI64+IekfY2zbtm3qvzp+/Li4QR8fH/GaHWOsqqqKMebr62t9qKamxvqjq6qqrA+tW7funXfe+fHHHz09PbndGjaq3dzcJkyY8PjjjzPG1Gr1rFmzRo0apdfrTSaTk5MTY6ympsbT05N6zVqt9uYOIAAA2DdxUGgKiYmJ33zzTcMWtfqP29w7dOhw/vx58evz58/7+vpaR7v27dtv3LhR/Lq+vj47O7tDhw6Msc2bN7/44otpaWmdO3cWu2VmZloHtfPnz3fo0KG0tPTYsWPic1Vt27YNCQkRvwkJCSktLQ0NDVWpVNaAaVZWVuvWram9V6lwAwYAAPzJhpsnBCYIwjUzQuv2HnzwwQ0bNly9etVsNi9ZsuTBBx9kjG3dunXdunVJSUnZ2dlpaWmMseXLl4eFhfXo0WPLli0zZsxISUkJDQ0tKyurqqoaMGCAu7v72rVrGWObN2/W6XTx8fGTJk3auXPn6dOnGWPq5OTkzZs3P/PMM2q1+rvvvuvbt6+7u/v999+/fPny9957Lz8/f9u2bdu3b1foiALYyqL0DlxD6doZALeGeL5T4lMa079//6eeeqp79+5OTk69e/f+/PPPGWMHDhwoLS2dNGnSmjVrJk+eLAiCp6fnhg0bGGNr1qyxWCziPRKMsX79+n3//fdr165NTk5+9dVXVSrV119/rdFoOnbsuHDhwoEDB6rVaqGurm7cuHGnT5/29vbW6XTfffdd586dz507d99993l5eeXk5EyZMmXRokXULiY/9uQDs16X9JoBbgUMhAByG90l5IZ9+g9MOXmqWNJmjYYT98QXbdmypZE+er1er9eLYdFrmM3mysrKm7nHr6yszMfHp2E2x2g01tbWql1cXLZt23b58mWDwdC+fXvxVGeXLl3Onz9/9uzZwMDAFi1aSHpJAADguIQmWYbJ2dnZ2dmZ+5BKpbrJO92tFxet1Gq1l5fXH+dh27Rpc83DTk5OERERN7NpAACAP9jhwryIugAAgEPDMkwAACAb8YZ6ic/BeoTgsIg8S1M2NznyF5raIfIJEnN39nYyCporoWmuETYpnBoFAACHhhkhAADIxw7DMhgIAQBANjcsos19ShPtzE3CQAgAALKxwwkhrhECAIBjw4wQpJMa1LTwH6H687tTGyG3QjXLkyclqynKlZeTuHkL7wH6hBMiqdBk5F6Y9xbAQAgAAPKx4RphE+3JTcOpUQAAcGiYEQIAgGwE6efSFT/3joEQAABkIwjSb4dQeiTEqVEAAHBomBECiQ5kSgtwms1Uf6LdLMNGyP68jduA+hOW+luYTofyH1ARf6OS2+f1l9SZ3pdGENuX1h2aF6RGAQDAkdnhOIhTowAA4NgwIwQAANnYUmtU6TkhBkIAAJCPHZ4bxUDoSCRWHaPCLxYih2I28dtNZDs/uMLdjslI/VD+RqgfSu08mQyiEL+6KhWVT6FCMRLbnST0pzo7SdlII9uRWh+ODNco/f8gyMiG+wgVHwhxjRAAABwaZoQAACAfO6w1ioEQAADkY4fXCHFqFAAAHBpmhAAAIBtBQNFtuA3IVRpNatrTZOC3G4l2g57qb7r5jVDtklOjVDu3la5eRpVGUznxH5Aa7HRS87fDbXdSS9uImty4xJ2hUqaSz0AhZWp/BCZIvS9Q8fsIcWoUAAAcGmaEAAAgHxvCMkrDQAgAALKx4Rqh4gMnBkIAAJCNwKTfR6j0VV9cIwQAAIeGGaEdI2uBkgveSqsFSqc9OalOxpi+Xmq7kd+u4/SnIqYmI9UucTVgibVGJS/MS/zNKTnAqSH689q5jY20a6h2rRO/3ZnfTm2fSp+SkVqBivIiTXoba4JrhB999NGbb75ZX18/bNiwVatWeXh4WB/as2fPY489lpOT06NHjzVr1nTt2jUzM3PBggW///67xWI5duwYY+zs2bMJCQkNN7hkyRJBEJ599lnxW8wIAQBAPsIfKzFJ0dj2jhw5Mm/evF9++aW4uLiurm7+/PnWh+rq6saPH79gwYKampr7779/2rRpjDGdTtetW7dp06ZdunRJ7NahQ4e0P33wwQe5ubn9+/evrq4ODQ0VGzEQAgDA7WvVqlXjxo3r2rWri4vLyy+/vGrVKutDqampgYGBY8eOValUL7744vHjx8+cORMRETFr1qyYmBhrN2dn5w5/2rt378iRI/38/Bhjbm5uYiMGQgAAkI0g/V/jMjIyIiIixK8jIiKKiorKy8utD3Xv3l382s3NrV27dhkZGY1sSq/Xr1u3bvr06eK3Bw4cCAoK6tGjB64RAgCAfGwqul1TU3PNGNa2bVtnZ2fGWFlZmfWioKenJ2OstLTUx8eHMVZeXu7u7m59ipeXV2lpaSM/Z8uWLS4uLnFxcYyxqKio3bt3d+jQ4eTJkxgIAQBAYfv377/vvvsatmzbti08PJwxFhAQUFFRITaKc8HAwEDx24CAgIsXL1qfUl5ebn2Ia+XKldOnT3dycmKMhYWFiY2DBw/GQGgHpKZDyZXiiYAlFcjkpjcZY7paftqzvo5oJ/rriP7cn0vWFCVTo/x2C7+ZWaiqomSxUaqZfICLrCkqNfDJC3BqnCV0Zow5E+lQZxd+u9aF/x8I1Z9OmXKbGSNSpnTNUqRJlXfj9Mv1T2FCfHz8li1buI+Gh4cfP35c/Pr48eMtW7YU54WMsa5duy5fvlz8urKy8tKlS127dqV+ytWrV3fv3v3pp59e/xCuEQIAgGyEP4vLSPrXiOnTp3/77bc///xzYWHhm2+++eijjzLGli9f/v7779933306ne6TTz6pqKiYO3du//79O3bsqNPpjhw5cvbsWZPJdOTIkdOnT4vbWbly5ZAhQzp06CB+u2nTpqNHjxYVFe3ZswczQgAAuH1169bt888/f/LJJ8vLy0eOHDl79mzGWG1tbU1NjUaj2bZt27PPPvvWW29FR0evXr2aMVZSUjJjxgzGWJcuXWbMmNGmTZtvv/2WMXby5MmZM2daN5uVlTVv3ryioqKWLVsK1Om1m5T82JMPzHr9b71KuBGcGmU4NfqnZnpqVNoN+IKKOgWKU6NNa3SXkBv2SRj93ZlzjSVWrldbdbRv76vUqdFbADNCAACQjQ3XCBX/UwUD4W1ErpkfWRqNmuERpc7qaogZXo1BWn9q5kf8XO4MVa4ZoeQFe7mt8tWQoiY3VIk1smoab/JHzfCokmlUfxc3/n8UBjf+QdYa+P21rkR/iftJ1aVT8bszhGhuJVtWqG+aPbl5CMsAAIBDw0AIAAAODadGAQBAPjbcR6j0uVG10Wisqqqyfu/u7i5WtTGZTOfOnWvRokVAQIByuwcAANC0VAcPHgwMDOz4JzHAmpGR0bVr1wkTJoSHh7/66qtK7yQAANgHW+6mV3xGyBjr3Lnzf//734atr7zyysiRIxctWpSXl9ejR4+xY8f27t1boT1shpRKh1L389VWUylQfjvVv76W366r4+8PtZ/chX+NxEK7cqVGqYNMpUlJEu8vpKqFqaj7C6Us5Evdt6fR8jdCpUZ1dfwaaNSSy64e/P4mI7+dOvjUsXfW8tvVROhBUppU8dN0zYG9HcM/rhFWVla6uLiIJ0VrampSU1PPnDnDGAsJCUlKSkpJScFACAAAzZKKMZaZmdm5c2dvb++kpKTCwsKcnByz2dy+fXuxR8eOHS9fvkw932wmCnUAAIDjkb48vfQb8OWmioiIyMvLy8/Pz8vLMxqNM2fOrK6udnZ2FheqYIy5ublVVlZSz9fpdLdqVwEAQEk3M/Ox4Rqh0uMgU/n4+Pj7+zPGfHx8Zs+evWvXrqCgoPr6+rq6OrFHaWlpcHAw9XxXV9dbtKcAAKAolap53nr+l1dVUFDg5eUVHBwcEhKyf/9+sXH//v24QAgAADfDhgmh4uka9YoVK1QqVceOHS9cuDBnzpyZM2c6OTk99dRTL7zwwscff5yenn7s2LFNmzYpu5f2i5sDbfJ0KFHbk0p71lTqJfWn0qRUKpWqKarX8V+X0cB7XUSAkKrVSUVzqTQpFYCkDj69P0S7xIURqNdFp0Y5/dUaaas9UGlS6uAY9PwUqJGM8hKffKrQq9TELnUwiXeFaia2Im1nHJfy45pk6pYtW/7rX/9au3ZtYGDg4sWLx48fzxibPXu2i4vL/PnzAwIC9uzZ4+vrq/R+AgCAPbCh6LbSA6d62LBhw4YNu6bVyclp1qxZs2bNUmSfAAAAbhnUGgUAANkI9IUA8imKzwgV/vkAANCc2OE1wuaZhQUAALhJmBHKQ1L5UKnpUD21sjyVDq3ipzqriXQolRqto1KjVE1RIjWqkxjI5J4ncfdy5nb29nfhtquIsy1lxXXc9pL8Wm47dZCpmqXUSR4VkQKl/nam9p9MjfKCoBpn/hGm0qRGI3/jVNpTeuFWbnMjvz/8ZjodSmyfeEDN+/+PKkxKVYuFa9iyQr3ShxYDIQAAyMaWkmlKj4Q4NQoAAA4NAyEAADg0nBoFAADZ2HKNsGn25OZhRggAAA4NM0KJqOKHROyNu745Ve5S6sryVKqTSodWV/DXzKqtlJYOpfaHSo1SqVeqwqSrO6eCpX+QG7dzx+5+3HaqnObF06Xc9oqSem479aZQL4osVEnWGuW3U1X+qdQodzF6E5ECVTvzj7zGxO9PriBPpUOlFgmlkAdHnnZB4Bw0MuhBVpHltzssgUkOyyi+HiEGQgAAkI8NN9Qr/ccEBkIAAJDPbbDQrlS4RggAAA4NM0IAAJCNDTfU4xqhnbEQaRki9sGvQWXQ88MyOmIBWyoUQ5VGIxfaJUIxtdX8/lQohmwnqpFReROKswsnwqDR8itfefpque0aoooYtR0q90G9WVRYhhuPYkxy7oMKxThT1cJ4+RRygVyyyB9VYo1aaJffTsWgpCJjK07EweStTtxIf+7qx9R/ysQ2GHWBS+n/20ECnBoFAACHhhkhAADIxqZTozfuU1hYWF5e3qlTp+s3XllZmZub265dOxcXfgn+RuTn51dXV2NGCAAAshEry0j61/jtExaL5cknn+zevfuIESPuuOOO3Nzcho+uWrWqXbt2EydObNeu3X/+8x/G2LFjx/r37+/u7u7j42Pt1qFDB29vbz8/Pz8/v7FjxzLGzGbzlClTevbsmZCQgIEQAABuX7t3705NTT1z5syZM2diYmL+8Y9/WB8qLy+fOXPmjz/+ePTo0ffee2/GjBkWi8XHx+e1115bt27dNdvZs2dPaWlpaWnppk2bGGNbtmz57bffMjIyzp07h4EQAABkJUj816j169ePHTvW39+fMfb444+vX7/eWrpo69at4eHhMTExjLFJkybl5eX9/vvv7dq1S0hICAoKumY7er2+urq64WaTk5M9PT0ZrhFS6IV2+e2S1tqlAod1NURqtIZKexKl0SQuqCtXOpQssaYnUqNUvpKX5Su8Wn19I2PMzYNTj43RaczcrEpue0Upv8Qa9WKNxIsyU2vPEq+WLJlGpFt9iIWIvQM47SZi6eOKUn6xPb2OWIVY4oK91K+J5FrMxB/qVApUTaRGnYj+3HYykkpUvaN20mFjo4Ig+XaIxvtfunRp9OjR4tdhYWHV1dUlJSUBAQHiQ2FhYeJDGo2mTZs22dnZkZGR3O0kJCQYDIaQkJBPP/00Pj7+0qVLSUlJ4kMYCAEAQD42rVCfn5+fmprasHHw4MFeXl6MserqamsKxtXVlTFWVVUlDoQ1NTVa7f9uoHJzc6uqquL+iJ07d4aFhVkslk8//XTMmDEZGRk1NTXWzWIgBAAAhWVnZy9fvrxhS/fu3cWBMCgoqKysTGwsLS0VW8RvW7Roce7cOetTSktLrz8jKhInjoIgPPXUUwsXLjx06FDDzWIgBAAA+dhQdJuxPn36bNmyhftQr169Dhw4IH594MCBrl27urm5WR9avHix2WxWqVQFBQWXL1++4447Gv9B9fX1lZWVXl5e4mYfffRRhoEQAABkJN4TIflJtIcffrh79+7Lly/v0qXLq6+++uyzzzLG3nzzzbKyssWLFwcGBj733HOTJk169913R40aFRoaWlVVlZKSkpWVpdfrly9f7u3t3b179y1btvTt29doNP7zn/9s27Ztnz59WrRoER0dPXjw4JYtW2IgBACA21fr1q1//PHHxYsXb9y48ZlnnnniiSfYn6kZQRC2b9/+xhtvvPbaa1FRUXPnzmWM6XS6I0eOMMYmT5585MiRoKCggQMH5ubmvvvuuxqNJjo6eu3atc7Ozl26dNm2bds///nP6upq4W8uoZn82JMPzHpdllerDDIdStQUJYooGoh4HjdIWVvFT29WEQvnVldIrClaRfQnfi6VMqUCk/p6Yq1aqp0IWFKZRm7g092Lnw719OHXGqXSm2WFtfz2Yn5qlHpR1CeEDEzym8l0aItW7tz2iBj+JZCOEZwFiqlPyLljxdz2vGx+0IAbfmaMaYjVj7Wu/L+wXdyIdt5SzIwxd0/iTScKzHoR7dSHxMPL+fpGVyKH7Ey8U1Tulxt+ZnYeJh3dJeSGfR6c8WPGhXJJmy0rOtytfSZ1avQWwIwQAABkZW8L8+KGegAAcGiYEQIAgGxsKbrdRLty0zAQAgCAbATpN9QrPhLi1CgAADg0R58RUpFZKvtHLe1NBSC5ZUWpNGY9UWu0XqYaoVQtUKm1Q6lF26lsIXVwdPVEGU9e/2oiUluUy0+Bmk38Hyo1Cit1sXXqpBBVkZLKFqqJLKJ/kBu3vU0nn+sby0v4UdirF/gFV6nfBwNRHZcqsUsdBLqMp7QaoWoirarR8IOdVLvamdOuduZvnE6HElVkqVmOXcdGb4YNN9QrfUgcfSAEAAA52XKNUOGREAMhAADIxg4nhLhGCAAAjg0zQgAAkI0Nt08oft3U4QdCiaXUJC3Ay4gIhuQFb+uokmZEf6rUGRF5oKIQVGk0o8SwDJVDIV8Xb3lY6shTqHfQZCDajfydp7ZD1yaUlhOhro5UE9XRCq7wFyi+6F16fWNVOT9hVFLATxhRS0CTYRni4DhRC+TqiRyKE3/7Tjp+f+oTXl/H33+NVsLqx85EWEZNJG6kJoDI9agVPz8oI3s7N4pTowAA4NAcfkYIAADyseGGeqUnhBgIAQBAPoJA3kjayFOUhVOjAADg0DAjBAAA+djhjYQOMxBSC/BS3YlMoIko3EVVHeMGO3VklS+ipBmV6tRJS29KrRZGHR1q+1TJNCodSpZY4x1MC5ne5DYzM5UHNlJ5YKI+HPVzJe4PmSEklOTzg53/TS/gtude4lRNo3LC5cQqxFQxP4qGKsmm5h9MJzV/f6jgJdWf/I2QmE/WuXA+gc4uRJ02LX8jVL03sxMRqSXjxkqPBjKxw3EQp0YBAMCxOcyMEAAAbgHcUA8AAI5O8XOdEv1vIFy+fLm7u3tycrL47bZt23bs2BEQEPDEE08EBwcrtHsAAGBPbLmPUOmB849rhBs2bHjllVeWLVsmfrtixYqnnnqqZ8+e+fn5/fv3r6urU24PAQAAmpCaMVZcXDx//vxZs2bt3LmTMWaxWN57771//vOfo0aNevTRR6Ojozdu3DhlyhSld/VvkboAr4Wfemuk1qiE8p5kzU8iemqgan4S6U3qzytXDw23nYrJUWvGUtk8utaohHQoIw4yWTuUqBEquXYo8Y5biM8O+cmhUqwSa6VSZUKpg1ZwlV+DlEtqtJVak5b6vSJToFTBVeITS/XXkDlqCdVrGWO6es4n30XHv2BEHXmTlorI8nfe0txrkNpQdFvyNUW5qRhjTz/99P/93/8FBgaKTYWFhRkZGUOHDhW/HTp06C+//KLYDgIAgB0RpP9Tmio1NbWiomLSpEnWpvz8fGdnZy8vL/HbwMDAvLw86vn19fwbkgAAoJkxm4kzJ3ZOPWvWrLS0tIZNzs7OJpPJbDarVCrGmF6vd3FxIZ+vRu4UAMAh3Mw5TBvCMopPCtU5OTlxcXGMscrKyurq6rCwsCNHjpjN5ry8vJYtWzLGrl69GhoaSj4fAyEAgGO4qYGQCdQqm/RTFKY6efJkWlpaWlrac88916NHj507d3p7e8fGxq5fv54xVl1d/f33348cOVLp/QQAAGgS6g4dOohfBQQEuLi4iN++8847iYmJBw4cOHPmTFRU1JAhQxTdSVlIC/nRGUVpK9Rzg6BUZ7Kd+KHUi3L3cua2+we7cds1zvzUaCERRKRSrNUV/KAjHdSUcPDpI8zfGWrj1B+fVGqUqi5LBpEp1M8lPoFktVviw8DNKJIpUKI8poZop3ZGRYRJqYNmJKK8Tgb+fkrNUdNVeSX010sMaVPvCHU1ja41SrTbHTssNvq/E5ujR4+OjY0Vv+7bt++ZM2d+/fXXoKCgPn36KJ5tBQAAOyH99gmlR8L/DYSBgYHWOygYYwEBASNGjFBilwAAAG4dRF0AAEA2dpkaVfjnAwBAc9IE1wjNZvO+ffsqKioGDhzo6+t7zaOnTp06f/58r1692rdvb228ePFiZWVlz549xW9ramrS09Orqqp69uzZunVrxlhubq71FnmHGQglFsSi0gHU8q10zoXTLrUqGNVOnYj38eff99mxuz+33dOHH65xded/PEoK+GvGFufy2ynUwee+XiofQeVBvAL4B8HTR8ttr6/lV+GiFsitq+GvYUu9KIpJYijGRJRqU/OqoFFrxlL/7zgRB9NC/D9BfQKpkmxUGEflRJVwkyezRlVH42fZiM5U0odKBtFLN/Pb6bFA6emSRLaUWGv0NRoMhuHDh5eUlLRv3/7RRx/dvXt39+7drY++/vrrX3311aBBgx5//PEPP/wwOTn5l19+GTFiRF1dnbOzc3l5OWPsxIkTgwYNuuOOO3x9fX/66ae5c+fOmjXrs88+++KLL1q1asUcaCAEAAA7tHXr1tzc3KNHj2q12jlz5sybN2/Tpk3iQzk5OYsXLz59+nT79u137tw5bdq08ePHd+3a9dChQ0VFRcOGDRO7hYSEnDp1Shzzfv7556FDhz7++OOMseTk5MWLFzMMhAAAICPZz4xu2bJl9OjRWq2WMTZx4sTevXubTCYnJyfG2Pfffx8VFSWeEY2Pj9fr9YcOHerXr19AQEBRUZF1Cw1zoGFhYUajUVxSqays7LfffmvdujUGQgAAkI/cI+HVq1cHDhwoft26dWu9Xl9UVCSukpuTkyNe8GOMqVSqli1bXr16tfEftXDhwuHDhwcEBDDG9u/ff+HChePHj2MgBAAAhWVmZi5cuLBhy9SpU1u0aMEY0+v11lqeGo2GMabT/VG7o+FD4qPWh7i+/PLL77777tdff2WMvfbaa2+++SZjrLS0FAMhAADIRrBhfUGB6XS64uLihm0m0x/5uJCQEOt5zsLCQkEQxOkgYyw4OPjkyZPWpxQVFTVSGUt4OHUAACAASURBVHvt2rVvvPHGnj17xIuFrq6uYrufn1+zGwipFKjE/mSJNSKzR7WbefE2uVKj5EqeRPaPGyxkjGk0/BJr1PKqVByOis9RB9NMZgI526F+qFeAK7e954AQbntYD350tii3htt+eC//TEv2uXJuO1Xli1ywl6rEZeQffLJenbOEuCr5zkpcaNeFyBX7teAX8/P05ueTqU8OtQQ0/QmUtnozN2VKl1KTsI40a+SdklYw0e5Co7bcRygwFhERcc2M0GrAgAHbt2+fPXs2Y2zPnj3R0dFarVb8xRkwYMAbb7xRX1/v4uKSmZlZVFTUu3dv7kY2b9788ssvp6Wlde7c+fpHm91ACAAAzciUKVMWLlz4/PPPd+nSZe7cuZ9//jlj7LnnnisuLl63bl10dPQDDzwwZsyYpUuXPvbYY76+vqWlpQsXLszNza2vr589e3ZAQEB8fPzEiROHDx/+r3/9S9zmzJkzZ82a1bt3b29v73379mEgBAAA+cgdlvHz8zt06NCyZctOnz6dkpIydOhQxtjIkSPF5OfWrVs///zzQ4cOPfXUU1OmTGGMOTk5+fr6+vr6zps3jzHm7e3t7e391ltvNdymWq2eOHHi/v37L1++3LNnTwyEAAAgHxtuqL9R/9atW18zkonL6DLG3NzcXnjhhYYPeXt7v/LKK9ds4fqWkSNHWlcYxEAIAACyscdao0RlIwAAAMfg6DNCCxHVopZpJesHSqkrSK8BKzWSyt9QWVEdt/3i6RJuu0bLT43mZFVy20sL+dvX1fHLdZoMEiK1jIrhEX9kevnya4d2uoOfDo0Z2orbfvk8PwWa9d9SbvuVjApuO5UVpFChUepTYiL+dlXxjj1ZC5RIgZo1xMK5RN7YL4ifDu3SK4DbHkAsDU19oqhPIPVJk7TUMyPKh0pfgFeG/xMYYxZiWiQ0lzTp7czRB0IAAJCRYIcL8+LUKAAAODTMCAEAQDa23FCv9GleDIQAACAbW9YjVHokxKlRAABwaJgRSkOuaC+lriAdVZWSomRMR1Ri1BGLqleU1EvafkUpv395Mb+dqgxJF2kklvbmHjTi4NTX8leKL87jH4Qrmfy0JxVQLC+WFpGlDiZVrpP8U5h4QEW0cw+a5JK5RH8qZepLFHrt0M2P2x7azpPbfpmI4JYT+WeyBin16ynl15b6WFJVYSX97jPWSOFjad2VPptIs6GyjNIwEAIAgHxwQz0AAIB9wYwQAABkIzBB6n2BSk8IMRACAICM5F594hbAQCgPade/Ja4eTF2Np2pBUTkOas1Yg56/nXpiO1RUgVyAl4psUCkD3pGg8h0lRDLo8B7+grpUyTSqylfWmTJuu544mGoN/3IDFRAXiLedWl2ZCt1w2yVfqpGYBTMY+AehtpqfYKqu0HPbqcQTWVxQclUzbjMfnayR1t9hCdI/eEqPg7hGCAAAjg0zQgAAkI0gSL9BXukb6jEQAgCAfOzwGiFOjQIAgEPDjBAAAGRjS9HtptmTm4eBUGrhK6KZegK3WUrfxhBxNSprRwUdqXajnipkJTEkS61EKinjR2yECiheOsdPe2Zn8BfgpaK2VESWqjqmItKe1JtCHUtqO2T6lNdObYSKnlKoT0hRTg23/fShAm57zgV+KTXqTaws03HbTURe2kxEdsmjzDsM1K+yUuFQxQcJqVB0GwAAwM5gIAQAAIeGU6MAACAnFN0GAACwJ5gRAgCAbGwJyzTRrtw0RxkIJa5+SqaYyBge1c7bjtTS7GTKlJjPS428kgsCSyyuSH/6Jawly5i0fB61EckFUYlUJ1U7VGo6lOLmqeG2u7jx243EKsf8iK/EfDKV76WWgC4p4Bd6NRB5YzcPZ267k1rab4STmv+mSC37yf/1JPO6xFbI/0Ok9W8+7PCGekcZCAEA4BbA7RMAAAB2BjNCAACQjQ2VZRSHGSEAADg0DIQAAODQmt2pUTIPJ892pK4bzg92Sg6NSis7KbVMJX0ig1g8nehNpVjJRdWJ/eQHOCUeNKlvOPkOEgdHajo0MNSN2x7WI4Db7unDD1jmZFXy2y9y2o16fsSU2nUqJ0xtp76GX6CVQtUs1bry/yPSuvDbnZz4HzXyA079RvDapX5cudFTxhqJ7NpbRWqJ7DEs0+wGQgAAUJTs6/JWVFSsX7++qqpq+PDhERERDR8ymUzffPNNVlZWZGRkfHy82FhYWPj777/X19ePHDnS2nPv3r2HDh1q27btmDFj1Go1Y8xoNG7cuPHKlSs4NQoAALevmpqaPn367Ny5s6SkZMCAAT/99FPDRx966KFFixbV19fPmDFjwYIFjLEdO3a0b9/++eefnzp1qrXbBx98MHXq1Lq6uo8++mjixIli4wMPPLB06dK6uro/ZoQGg0GtVl8zP9Xr9c7O/PMzAAAAHHLfUL9u3TofH5/NmzcLghAcHPzWW2/FxsaKD509e3bLli1Xr1719fVNSkqKi4t75plnBg8eXFFRcfjw4WHDhond6uvrFyxYsHXr1n79+s2aNat169YnTpwwGAw///zzlStXPDw8VN99913r1q29vb3d3d0TExMLCgoYY0VFRfHx8QEBAQEBAcuXL5d8IAAAwCEJ0jW+wbS0tISEBLFbQkLC3r17DQaD9aH+/fv7+voyxnr16uXm5paenu7i4iKe+bQ6evSoxWK56667GGOenp6DBg1KS0tLS0uLjY318PBgjKl69uy5f//+2trawsJCxtgrr7zCGHv55ZdDQkLKysr27t374osvZmRkyH+0AAAAbiQ3Nzc4OFj8OiQkxGw2ixM2xlheXp71IcZYcHBwbm4udwtBQUHWETc4ODgvLy83NzckJERsUbdv3178ysPDIz4+/vvvv6+vr9+wYcPBgwednJx69OiRkJCwZs2aN998s4le5C0iMS8pudYokTTjFkWkFjens23c5kZ2nmqX9qKoio7U0ST3U47cG1UG00S0m03EIuZEMJI+aMTK8sR23L34VxM63cFPh949Lozb7hPgym0/uOsKt726Qn99Y1lRHbez1JK5khdtJx6gDhrVTpErn8z9TaQiqdSvA/lDJadJmwkbbqgXBHb06NHp06c3bJw3b17r1q0ZYyqVyvrrbzabWYP/ygRBaPg/g8Vi4f4v13AL1m4NG9WMsbKysrS0tOzs7FWrVn388ce5ubn19fXh4eFij/Dw8LNnz1IvgPrvCQAAHJFN1wj9/f0HDBjQsM3d3V38IjQ0ND8/X/w6Pz/fyckpKCjI+tDBgwetT8nPzw8NDb1+8yEhIQUFBdZhMi8vLzw83GAw7N+/X+ygZoxVVVXt3r378uXLGo3G09OzvLxcq9Vaz7F6eXmVlpZS+19byy8/DwAAzYzZbFapbnCvgSD9JJDAhDZt2lwzI7S69957P/300zlz5qhUqm3btsXHx6vV6ry8PKPReM8997z88svFxcUBAQHp6ek6nS4mJub6LfTu3VutVu/bt2/QoEEVFRU///zzwoULDQbDggULKisrvby81IyxNm3aLFu2jDG2bNmyqVOnbt++vb6+XqfTabVaxlhZWZl1+L2eddAGAIDm7YajYFOYMGHChx9+mJiY2Llz59WrV6empjLGFixYUFxcvG7dugkTJgwdOvS+++5bt27dnDlzXF1dCwsLZ86cWVJSUltbO27cuJCQkI8++mju3LkTJ06cPHlyWlpaUlJSt27dGGP33ntvXFxcfHz8X6I1EREReXl5ISEhPj4+x44d69OnD2Ps2LFjYtgGAACgcYINN9Q3+qirq+v+/fu3bNlSVlaWnp7esWNHxtiMGTP0ej1j7Msvv/zhhx8yMzNTUlL69evHGPP09HzssccYY6+++ipjzM3NjTH21FNPRUZGHjx48I033khISBC3nJKS8v3332dlZanXrVvXoUOHNm3aXLly5bXXXktKStJoNNOmTXv99ddXrVqVnp6+d+/ezz77TPLBuO1IvHBNXY2XctWdMabmhmWItV65nVlj9dukll7jNpMpAI3Widvu6k4UviK2U12h47aXFvCjHFXlnP5kWMZIXKU2UgEPallhopk6mBp+u4sb/+B4+Wq57X5B/NJrPgEu3HZPH/52nHlvlvRsl7TYlJroT61mTLXTP1fa/lA5FzLPwg3LSNwZ6b+e3ObmlJax6T7CRtMm7u7uycnJDVus9WUEQbAObCJXV1driZmG+vbt27dv34YtKpUqKSmJMaaurq5+/vnn8/LyAgIC7rvvPvH2ibfeeuvFF1/s379/YGDgpk2buJcfAQAAmgH1Y489Js4iG3Jzc/v0008V2SEAALBj0otuK76AIYpuAwCAbGy5j7Bp9uTmoeg2AAA4NMwIAQBANjatR9hE+3KzHH0gpI4/VR5JajUybkCUzM6RmbqmbdcSL9bbnx9cbB/uy233bcGvClZwpZrb/t/0Qm67Qc9ZvtVk5JfholKjZjORMlXxt0O+4xIX7DUa+NunFtQ9RJRM8/Dml2rLPFnCba+p5JRYo9DhZ6mpThlSmo1sh8pRS06fSvm5ZORVYlSVSp47QGzU/uDUKAAAODRHnxECAICMbAjLKD4ZxkAIAACyEZj0a4RKj4QYCAEAQD5yr1B/C+AaIQAAODRHmRGSK4tKrTUqNfbGS6BpNPwanhpnidk5iZUeNc78n0uV8fQJ5KdGO3b357a3CvPmtl88zQ865l7iBykLcjgpUyqbp9YQC+cSL0pl5DaTJUilpka5kVfG2NUL/BfLXVCX0W9WdSW/cGtdDeeFkQFI4pNGfrwlrlVL1vyUuCCwLL9urJHqvvxQN//IUxuRujCv4rcKNDXbFuZVlqMMhAAAcCvYUGJN6XOjODUKAAAODTNCAACQjx2GZTAQAgCAbOxwHMSpUQAAcGyOPiOkrumS65JLXNqbnxrVSkt1kmlSMn1KleXkb8dIlfE08JOUNVX8oGNlST23vbqc39+g4/9c7qGkjjwVFKRQ7yxVm5T8JJCVJPnNBh0/TVpaUMvfDBGkpHCPg9aV/wtOfaKotKdGy/+kubrxt+9M/FyyXKfEWqZ0Gpa/n87kbxanXWqtUckr1HNbb4PkpFwE+reGfIrSr93RB0IAAJCTHZ4bxUAIAACysaXWqNJwjRAAABwaZoQAACAbmxbmRdFthVEXuvm9Ja3wyaiwDHVJ34W4pE9EFdQafv6CCtGonYm1ak389opSfvjl3NEibvvVzApue0k+Pw9SToRruNXOqAJX1NUFqWvDmomDQCEjEkR/autSq4tR/blV06hQjDMVfnHXcNs9fbXcdjdP/urBVNCH+qRRuRI6Pkb8BhGvi/oN4raTSSIqLEOFYhx2/V07vEaIU6MAAODQMCMEAAAZCVLXF8R6hAAA0HwIt8F9gVLh1CgAADg0zAgBAEA+dhiWwUDIR03tycwekSjT8LKOVLZNajsVbzMSpdQ0JqK/gd+uq+cvYpt3qYrbTgUj62oM3Pb6Wn479yBzjyRjzIl4R6iFdtUafl03qalRCrUgsIX/Y8nAJFl1jFrzlpcypTpT6dCgNp7c9tYd+Usue3jzU6NUEb6KYn5OWE/Un6PSns5EkUJZfoPI1X0l5pDpWwKU/l+/idnj7RM4NQoAAA4NM0IAAJCNDSXWFJ8jYyAEAAD54BohAAA4MhuuEd5wIMzMzPzss8/KyspGjRqVmJjY8KHKysqPPvooMzMzKirqiSeeUKvVjLGMjIzPP/+8Yf/PP/+8ouJ/pa86d+6cmJi4ePFi8VtcIwQAgNtXcXFxv379tFptbGzsY489tmHDhoaPJiYmHjt2LCEhISUl5dlnn7X2d3FxEftv3Ljxmg0uWbLk1KlTBoNh9uzZBoOBYUYo9Q8XqQvzcmNvVE1RF2KZU10dv11fz8/aURUdqWCkiVqYl0ifUhk/au1Zqp3CLepIpUAtav4D5DslJXXJ6BQo9aJ0dfx2o4FYhZj4BFKVLSWlSanUJVU7tHVHL257REwLSdspuFzNbc8WyrntVeU6bjsV1HTWEgsOy5HHpkr1ylVrVOmAZJOT/czoqlWrIiMj33nnHfHbRYsWjR8/Xvz64MGDJ0+ezMvL02q1MTEx4eHh8+fPX7lyZUxMzNtvv80Ys1gsixYtGjdu3OOPPy4+5cqVK6+//vqUKVPEb2fNmuXqSiwkDQAAYCNZB/vffvttyJAh4tdDhgyZOnWqTqfTarXiQ+JkkTHWrl274ODgo0eP7t+/PzY21tp/2rRp1v6MsRUrVgwdOrRNmzZ1dXWMsbffflulUuHUKAAA3L7y8/P9/f3FrwMCAiwWS15e3vUPMcYCAwPz8vKu75+fny9+azabV61aNX36dMaYIAjjx4/38/MzmUyYEQIAgHxsuqF+3759/fr1a9i4evXqTp06McZcXFz0+j9KNOh0OsaYq6ur+K2Li4t4kc/6qIuLi1arpfrv2rWruro6KSlJfO769evFdgyEAAAgG1vuIxRYRETEu+++27CxZcuW4hetWrW6cuWK+PWVK1e0Wm1gYKD1od27d4tfm83mnJycVq1aXdPfxcUlICBA/HblypWTJ0+2nia1wqlRAABQmJ+fX7+/cnNzEx8aNWrUN998I17SW7t27ciRI1Uq1bFjxw4ePJiYmPj7779nZGQwxrZv3+7q6hodHT169Ojr+zPGSktLt27dOnXqVHGz5eXlZvMfETbMCCnSkmBUooy7lLazC/+wu7jyg4V6dyKNqSfSnrKlRqmymdK2byb6a6jqnrx26nI2dRLGxZ1/kP1auHHbfQJcuO3UQS7KreG2lxbUctvravmFW6kSrVT1WjpNymmncshuHvxao+7e/BSopw/RTvSv9OCnQLVEOo+qaktFfKnctdR2SSvUq4i8rqOWFKXZEBttVGJi4sqVKyMjI9u0aXPy5ElxCrhy5cri4uJ169bNmTNn0KBBMTExv/3222effaZWq5OSklauXBkVFdW6deuTJ0/u2bNH3M7XX3/do0ePO++8U/x206ZN8+bN69GjR0lJCQZCAACQjSBIL6LdaH+1Wr1t27bff/+9oqKiT58+7u7ujLF//OMfJpOJMfbaa6+NGzcuKytr2bJlwcHBYv/U1NRr+jPGxo8fP3HiROtmH3300djY2KysLG9vbwyEAABwWxMEITIysmFLw7BoWFhYWFhY4/0ZYy1aXHs7bKdOncQ8DgZCAACQjQ0r1CteZAADIQAAyEZggiDxIqFAr2Z6a2Ag5CPfRiosQ9SC4i7ySV26Nxr47a4GfrTBSIRcqNAKVS2MCrNQVc0o1IUBqiqYQc3ff261M5OJ35kKj/gToZjOvQK47R3C/bjt1Bqzpw8VcNuNeqKeHPWRIg6yk8TlYbm5Ei3xSaNWA66t5L/Y/Mv8pZgrS/mhmNIifmLIQBwccl1rojSa1pVoJ2JoVDtRYk3aEtDUwXTctIwdrj6B2ycAAMChYUYIAADykX5DveIwEAIAgGxsWI9QEASFrxFaLJb09PT09HRBEGJjY7t16yY+cOLEibS0tMDAwLFjx1oLtQEAADQzqqVLl06ePPn06dMnT5686667Vq1axRhLTU0dMmRIUVHR119/HRsbazRSRTEAAAAaEGz6pyj1xIkTn3zySbEUW3R09IIFC6ZOnTpv3rz3339/2rRpRqOxR48eqampo0aNUnhPbzGJsVGBSB1xA5POvLprjDETUYCKKnVGl0Dj7wyJOishcWVRci1cPVUVjJ8h5AYjqYishlhG1cPHmdseGOLObQ/twF+TtqqMH4zMuVjJbXf15Ed89URgknqzqOwidZC5y8NSa8xSueLyknpuu5BRwW2n0ptUXTpqdWJuJULGmAvxG0GVaqNKylH9ualRcgFe4siTmVGl/3NXii1Ft5tmT26eyt/f37oqoZubmyAIJSUlR44cSUxMZIyp1erhw4fv3LlT0Z0EAABoKv/726e2tnb+/PnPPPNMbm6uRqOxFrAJCQnJzc2lnm9d9gkAAJo36l7khoQ/5oQS/ik+J/xjINTr9WPHju3Zs+fjjz8uCELDV2uxWCRXUAUAAIdlb9cIVYwxg8EwceJErVa7evVqlUoVEhJiNBpLSkrEHvn5+SEhIdTznZ35V2IAAKCZaa6TIpXJZJo6dWptbW1KSopGo2GM+fv7R0ZGpqamMsaMRuP27dvvuecepfcTAADsgOQTo7fBDfjqjz/+OCUlJTExcfLkyYwxQRA2bNjwj3/8Y9q0aWfPnj127Ji3t7cYnAEmPQmm4lWAVPNzhfwMG2tkwVuinS4SKrV4KNFMHAVucJEx5uTED0ySMTxes5OB+KFExo+K1JYW1XHbr2SUc9vriQV162oN3HaqbCYVXKRSo1SBVjLTyDv41Eaod9Cg479TVeX86Kyunv+JJWuHUulQIu3pSiwg7OrOb6dTo8TCvLz9kXrQlP9f/DZjyw31TbQrN019zz33bNiw4ZrWxMTEn376KS0t7aGHHho9erRajQI0AABwE+yw6La6W7du1moyDXXv3r179+63focAAABuJUz1AABAPjZc81N8RqjwzwcAgGZEkB4ulbqQr+ywHiEAADg0zAjlQmQauauLELE0KjVKhUClriAvNQVK1hQl06FEu06GWqNU+UrqKOh1/LRnzkV+2cxyIk1K1TilapBSqEXSKVR2kQpkciO4VMSUSm9SK8KT7ywvFN3IdqjaoVLToVR/KjXqTBx8bkFX6uNN1RNGaLQZwEAIAACysW09wibamZuEU6MAAODQMCMEAADZ2FIpRunTyxgIAQBAPrdBEW2pMBDKg/oLyML7RNCr+PIfcNbKszMq4udy68CxRhbaJfaTWktWV88PxVAFvTS8dgOxsC29OjG/nSqZpqujFs7lb4daY5Y6OFSuROqHgTrI3P5UZyosQ8W1yHYXfruWaqdKqUkMxZAl1ogwjsZZwidcIMIy9ve/u0JsukbYRPtys3CNEAAAHBpmhAAAIB/p1wgxIwQAAFASBkIAAHBoODUKAACyEZgN6xEqfG4UA2HTooqXcVtV/KgdUxMTd6kLhVK1o6i1bdVUlS8yi8hvr6/jBzX1RJpUz02NEhFTKr1pIlYtplKgVL06M9GfOgjUdsg3RWL1MrL0Gm9/1DKlQ6mSaXRqlP8fC1UCTXI7lQ6lSsRJqUsnteIgXMMe7yPEqVEAALit/fTTT8OHD+/Xr9/ChQvN5r/87ZuVlZWcnBwTE/PUU0+Vl5eLjXv37hX7L1q0SOy/atWqcQ0YjUbGWHZ29oMPPhgTE4OBEAAA5CPY9I+WnZ2dlJQ0adKkJUuWrFmz5pNPPrE+ZDabExISWrZs+eWXXxYXFz/yyCOMsUuXLo0YMSI5OXnJkiWrVq1aunQpY+zYsWMWi2Xsn1QqlcViSUhICAoK+vLLL3FqFAAAZCN70e0vv/xy+PDhkydPZoy9/fbbL7300jPPPCM+tGvXroqKivfee08QhE8++aRVq1ZXr1794osvEhISHnzwQbH/K6+88vTTTzPGunXrNnbsWOtmd+/eXVJSsmjRIpWKqjUCAABwGzh+/HifPn3Er2NiYjIyMmpqaqwPRUdHi+NoYGBg27ZtT506deLEiZiYGGv/8+fP19bWMsa2bduWmJj49NNPnzt3TnxuVFSUOAhiRggAALKxISzTePfCwkIfHx/xaz8/P8ZYQUFBhw4dxId8fX2tPX19fQsKCrj9BwwYEBUV5efnl5aWFhkZeeTIkYbPxUCoAPpTIi1NKkhcKpRMjZLBRam1RqnMoYHbTpX35K6py42SMsaMxIK9JmJBXSoFKnn1Y4nLIlMVLGV7U3ipUWeJC/DKV2uU/x+Llkh7al2p7UurHSopHcqIM3JIh8pAemr0hx9+aDikMcb27dvXvXt3xpiXl5c4pWOMVVdXM8as45yXl9fly5etT6mpqfH29vb09Ly+/5gxY8SW++677/LlyytXrvT29s7MzBQbMRACAICcbLgv8O677167dm3DFi8vL/GLdu3aWUeszMxMLy8v65DZrl277777Tvxar9dfvny5bdu21/T39va2DpyiVq1alZeX33nnnRs2bBBbcI0QAAAUptFofP/KyemP8wSTJk3auHFjYWEhY2zp0qUTJ04UBGHHjh3fffddUlJSRkbGr7/+yhhbs2ZN69ate/XqJfYvKipq2P/IkSPi1s6ePbthw4bBgwcnJiZmZWX9/PPPDDNCAACQkS3XCBvtHxsbO2nSpK5du3p5eQUFBW3bto0x9uOPP5aUlIwaNeqLL75ISkoKCAiora3duHEjYywuLm7ChAldunQR+6empjLGpk6dmpub6+HhUVpa+swzz4ij45dffjlq1KjAwEDBIvFSxzWSH3vygVmv/50tgBV9mUpakRQjsUSfkVjSj7r8piMqwlBL+tXVUP1xjRDXCBnDNUL7N7pLyA37rNyRUVheL2mz/03fU5m5d8uWLY30qaioqKmpCQ0Nvf4hnU6Xn5/fqlUr6ySS27+0tFRsadhNr9fn5eVhRggAALc7b29vb29v7kNarbZt27Y37O/n5yeGSBtydnZu27YtBsLbiNQ0KRUaVRPvqiDw/+4mJyvE39FqKjVK/J3urCVWrnchZqL1nP2kao0aiBkhVYPUSMwUqek1iaodSlWqJA4yNcMjZ4RSVq6nVqKnZoTUTEtyDVKptUyJ/aQKulIHR9LMj2Hy10RuVCmG/xRFYSAEAADZoOg2AACAncGMEAAAZGOP6xFiRggAAA4NAyEAADg0nBq1A3LVJqU2RJ3HkC01Sq2TTqRGXXScj6WeuAmSuo+QTI0S7WYTkSaVuOK81AiuLCvRM6LQq4bqTLxTkjbO6JQptR0q8kq9WOomS+ogIx16O5D9hvpbAAMhAADIR/p6hEpfIsRACAAAMrLhPkKl4RohAAA4NMwIAQBANrhGCLcWldegQjTS1vFlKuIJTlS1MMkL9krIvxiIkIuJaifC92HftQAAENRJREFULyaiIrnZTBXp5jaT6IMp7aCRYRmqPy+fQnaW+EOpRI+K7C8x/ELWpeM2M7LooNL/nwL74xKhxPsIlX7ncGoUAAAcGmaEAAAgK3ubmmMgBAAA2dhhaBSnRgEAwLFhRggAALKxx7AMBkIHItfCv3TGjx/IlFqqzaTlJDWptCeVDjWbqHQov51amJcqsUaSuGCv1DVmyQAnrz+5cYl14Oi0J7e5kRJo/P5IgTYrdrgwL06NAgCAQ8OMEAAA5CNIPtWJU6MAANB82FJZpmn25OapqqqqFixYMHbs2Lvvvru8vFxsra2tnTFjRtu2baOjo7dv367sLgIAgD0RpP9TlKqysvLy5cv9+vXbtWuXXq8XW+fMmZOVlbV///7XX399/Pjxubm5yu4lAABAE1G3bNnys88+KykpeeGFF8Qmg8Hw1Vdf/fjjj6GhoUlJSUOGDFm9evWrr76q7I5C0yHPY1APEEFKInJIZg7NTvwNOZk5/aman1Tas8nToQT6YPKbqTQpGbwk8m3c7UhNb0pe2BYpULiOwASq3PFti3ONMDc3t7y8vFevXuK3vXr1OnPmzK3dKwAAsEv2uPoE58/LkpISFxcXZ2dn8VsfH5/CwkLq+TU1NU21awAAcDuhFmmxd5yB0NfXt76+3mAwiN9WVFT4+/tTz3dzc2uqXQMAgNsJtTTbX9hbUoZxB8KQkBBXV9ezZ8+K3549e7Zjx47U8xW//wMAAG4fgk2U3WcVYywrKys7O5sxlp2dffHiRRcXl/Hjxy9cuNBkMp08efKHH36YPHmysnsJAADQRNSMsWHDhhmNxg4dOkyYMEEQhMzMzIULFyYnJ/v5+Tk7O7///vudOnVSej/hNiJbypRot/CyjkQdUzLuSaVALeQDxPalogKZ9AP8ZokHmdssMdRJPqD03+tgT2wIyyh+dlTNGDt37tw1rYGBgTt37jQYDBqNRom9AgAA+3R7XPaTpLESaxgFAQBAIvtbhgmrTwAAwG1t5cqV7du3DwgImD59el1dXcOHDhw40KtXLz8/v6FDh2ZlZYmNK1asEPs//PDD9fX1jLE33nijR48ePj4+d9xxx7p16xhjp0+fjvoTBkIAAJCNDXdPND4fPHHixAsvvLBx48YLFy5cunTp7bfftj6k1+tHjRo1c+bM3NzcqKioqVOniv1ffPHFTZs2ZWZmXrx4Ueyv0+m++uqrq1evLliw4JFHHklPT6+pqSksLFy2bNmyZcuw+gQ0OVnCNZLPnZCZmNvr8oVceRZ+19vrtYJjkHth3q+++mrMmDHR0dGMsf/7v/+bPHnyW2+9JT70/fffe3p6Pvzww4yxOXPmBAYGZmRkrFixYsyYMVFRUWL/KVOmvPnmmwsWLBCfkpCQ0Lt37/T09KioKK1WGxkZyXBqFAAAbmdnz5694447xK/vuOOOvLy8iooK8dtz587deeed4teenp7t27c/e/bsuXPnGvbPzc2trKy0bq2wsPDkyZPiMJmbm9utW7fo6GjMCAEAQDa23SCv0+kKCgoatvj7+6vVasZYaWmpp6en2Ojl5cUYKykp8fb2Fh/y8PCwPsXb27ukpITbX/xCr9cnJydPmDAhJiYmNzd306ZN4eHh2dnZGAgBAEA2thXd3rt3r3UaJ9qzZ09ERARjzN/f3zqlE+eCAQEB4rd+fn7WgIz4aEBAgJ+fH7e/0WicNGmSu7v7J598whgLDQ0NDQ1ljLVv3x4DIQAAKGzYsGFbtmzhPtS5c+dTp06JX586dSooKEic3okPrVq1Svy6urr60qVLnTp16ty58+nTp639g4ODPT09TSbTtGnTqqqqtm7dev2dgbhGCAAAshEEW8qNNrLBadOmbd68+ejRo1VVVQsWLJg+fTpjbM2aNUuXLr3//vvLy8tXr16t1+vfeeedyMjILl26TJ8+fePGjceOHauqqnr33XfF/tOnTz9x4sSHH36Ym5t78eLFsrKyvXv3nj59WqfTnT17FjNCuO1IvRlXWrNMFdZky2Mi2AnNj6yf6jvvvHPhwoUjRoyoqqpKTEycM2cO+3PdXK1W++233z7xxBNPP/10r1691qxZwxjr2bPne++9l5SUVFVVlZSUJPY/f/68RqN56KGHxG0+8sgjTk5ODz/8cH5+fosWLQSy9OLNSX7syQdmvf63XynArYKBEMBWo7uE3LDPt0eulNToJW322L60nPQfqVOjtwBmhAAAICe7++sOAyEAAMjGhtsnFK81ioEQAADkI3dlmVsAqVEAAHBomBGCg5EWMgUAaWy7oV5ZGAgBAEA2tpVYUxZOjQIAgEPDjBAAAORjh2EZDIQAACAbQfo1P8XPpOLUKAAAODTMCAEAQDYCEwSJ5zql9pcdBkIAAJCPDdcIlYZTowAA4NAwIwQAAPlIv6Fe8RkkBkIAAJANim4DAIBjwzVCAAAA+4IZIQAAyEZgNpwabaJ9uVkYCAEAQDZ2WGENp0YBAMCxYUYIAADyscMpIQZCAACQjSDYcDuEwiMhTo0CAIBDw4wQAADkI72yDFKjAADQjNjhNUKcGgUAAIeGGSEAAMjGlhvqlZ4SYiAEAADZCHa4+gROjQIAgEPDQAgAAA4Np0YBAEA2Nq1H2ET7crMwEAIAgGxsuEao9DiIU6MAAODYMCMEAAD52OEN9RgIAQBAPtKvESp+kRCnRgEAwKFhRggAALIRpE/wlJ4Q/u2B0GQyybIfAABwm6upqXF3d2+8jy23Tyh9kVCwWCx/5/kpKSmffPKJWt1sZ5Z1dXUqlUqr1Sq9I7eCTqczm82urq5K78itYDAYdDqdh4eH0jtyK5hMppqaGi8vL6V35FawWCwVFRU+Pj5K78gtUlFR4eXlJX0tXFusWLEiLCxM9s0WFxfn5+d3795d9i3fpL87EAIAANg1hGUAAMChYSAEAACHhoEQAAAcGgZCAABwaM027fl3XL16taCgoGvXrtagcGlp6e7duz08POLi4ppTglSn0x04cCA/Pz8sLCwyMlJsNBqNe/bsKS8vj42NbdGihbJ7KKPS0tLDhw+XlpYGBwcPHDjQyclJbD906NDFixd79uzZtWtXZfewKVy4cKGurs6axysuLt67d6+Xl1dcXJxGo1F23+RSWlqalZVl/bZLly5iGDgzM/Pw4cNt2rTp16+fcnvXJC5dunTgwAEPD4+77rrL39+fMVZZWblr1y6NRhMfH+8gwW8ZITX6F3q9vk2bNrW1tTU1NQcOHIiOjmaMnTt3btCgQYMHD87Pz6+trf3555/d3NyU3lMZmM1mf3//8PDw9u3b79u3LyYmZuPGjSaT6e67766rqwsLC9uxY8eOHTusA6S9Gz58OGMsMDDw2LFjarX6p59+8vLyevnllzdt2hQbG/vDDz8sWLBg+vTpSu+mnLKysnr16tWiRYvz588zxk6dOjVkyJChQ4deuXLFYrHs3bu3efxVt3bt2ueff75nz57it0uWLImIiFi/fv3TTz+dkJDw22+/xcfHf/rpp8rupIwWL168YMGC2NhYnU7Xrl27jz/++MqVK/369YuMjKypqcnJyfn11199fX2V3k27YoEGTCbTmTNnTCaTp6fnoUOHxMaHHnroueeeEx/t27fvsmXLFN1H2ZjN5jNnzohfl5SUeHt779mzZ+PGjeHh4TqdzmKxvPnmm/fff7+i+9gkDAZDly5dVq9effnyZVdX1+zsbIvFsmfPnqCgIPGFNw9ms/mee+555plnOnXqJLaMGzdu9uzZFovFYDD07t179erViu6gbP71r3/dd999DVuMRmPbtm23bNlisVgKCgo8PT3Pnj2r0N7J7ODBg15eXhcuXGjY+Oyzz06ZMsVisZjN5nvvvXfBggXK7JzdwjXCv1CpVF27dlWp/nJYtm3bNm7cOPHRMWPGpKamKrR3MhMEwXoy0M/Pz8PDo66uLjU1deTIkc7OzoyxsWPHbt++3WAwKLqb8rNYLGaz2dfX94cffoiKimrTpg1jLDY21mKxHDx4UOm9k83y5cvbtWsXFxcnfmuxWFJTU8VPslqtHj169LZt2xTdQTlVV1fv3bv3xIkTZrOZMXbixIni4uKEhATGWIsWLQYPHtxsfm3XrVs3YcIEDw+Pffv2FRYWio3W/6MEQRg7dmyzebG3DAbCG6ipqSkvL2/VqpX4batWrXJycpTdpaawatUqZ2fn2NjYnJwc64tt3bq1yWTKz89Xdt9ktHbt2lGjRoWHh48fP/7+++/Pyclp3bq1+JAgCC1btmw2b25ubu6SJUveffdda0tZWVldXV1z/SQXFBS8++67999/f3R0dF5eXm5ubnBwsLXiVatWrXJzc5XdQ7lcuHDh3LlzcXFxixcvDg8PX7duncViycvLa/hr25ze2VsDYZkbEOdD1mCFWq3W6/WK7pH8fv7555dffnnr1q1ubm5Go9H6YsUvmtOMsE+fPr6+vunp6cuWLZswYYLBYLC+WMaYWq1uNi/2iSeeePvttxteKGrGn+Tx48c/+OCDjDGDwfDAAw+89tprI0aMuOadbTYvtra29urVq6dOnXJxcfn3v/+dnJw8evTohr+2zenF3jKYEd6Aj4+Pm5tbUVGR+G1BQUFoaKiyuySvAwcOjB07NiUl5a677mKMBQcHW1+seOIlJCREyf2TVadOnRISEt54442EhIRPP/00JCTEenKJMVZYWNg8Xuzx48f37t27Y8eOGTNmfPbZZ4WFhTNmzNBqtWq1ull+kq3xV41GM2HChN9//z0kJKSoqMjyZxKwsLCw2bzY0NDQfv36ubi4MMbi4+MrKiquXLkSFBTULN/ZWwYD4Y0NGjQoLS1N/HrXrl2DBg1Sdn9kdPTo0ZEjR65YsWLo0KFiyzUvNiYmpllGsSsqKjw8PAYNGvTbb7/V1NQwxs6fP19cXBwVFaX0rskgNDT0/fffj4yMjIyM7Nixo1arjYyM1Gq1AwYMaK6fZKsTJ060bNmyR48ejLEjR44wxvR6/X/+859m82JjY2Ot94pcvHjRyckpODj4ml/bwYMHK7eDdgm3T1xr0aJFJSUlH3744YQJE0JCQmbPnn306NFRo0bNmTMnNzf366+/PnHiRFBQkNK7KYPa2tpWrVp17tw5NjZWbBkxYkRERESPHj2GDRvWrVu3t95664svvhg5cqSiuymP+vr6YcOGxcfHe3t779+/f+fOnQcPHuzYseP999+v0+lGjx79+eef33vvvQsXLlR6T2W2devWl156Sbx94scff5w0adKcOXOysrK++eabkydPireg2bsnn3zSy8urZcuWx48fX79+/Y4dO/r37z9//vyUlJSnn376hx9+qK6u/umnn5TeTXnU1dXdeeedQ4cO7d2799KlS+Pi4j744IMjR47ExcW98sorNTU1n3322eHDhzt06KD0ntoTDITXWr9+fWVlpfXbSZMmeXh4pKenb9q0ycPDY+rUqWLIsBnQ6XSrV69u2NK/f/+IiIj8/PyvvvqqpKRkxIgRAwcOVGr35GWxWL7//vvDhw9XVVWFhYVNmDDBz8+PMabT6VasWJGZmRkVFTVx4sRbs5bNrZSdnX3gwIHx48eL3+7fv//bb7/18fGZOnVqy5Ytld03uRw4cCAtLa2srCw0NHTcuHHW39DNmzf/9ttv7dq1e+SRR5rHvb+i0tLSFStWFBcX33XXXda/U0+cOJGSkqLRaB566KGmWCmpecNACAAADg3XCAEAwKFhIAQAAIeGgRAAABza/wNZryc5g0JXWwAAAABJRU5ErkJggg==",
"text/html": [
"\n",
"\n"
],
"image/svg+xml": [
"\n",
"\n"
]
},
"metadata": {},
"execution_count": 5
}
],
"cell_type": "code",
"source": [
"terms = [Kinetic(),\n",
" ExternalFromReal(X -> pot(X...)),\n",
" PowerNonlinearity(C, α),\n",
" Magnetic(Apot),\n",
"]\n",
"model = Model(lattice; n_electrons=n_electrons,\n",
" terms=terms, spin_polarization=:spinless) # \"spinless electrons\"\n",
"basis = PlaneWaveBasis(model, Ecut, kgrid=(1, 1, 1))\n",
"scfres = direct_minimization(basis, tol=1e-5) # Reduce tol for production\n",
"heatmap(scfres.ρ.real[:, :, 1], c=:blues)"
],
"metadata": {},
"execution_count": 5
}
],
"nbformat_minor": 3,
"metadata": {
"language_info": {
"file_extension": ".jl",
"mimetype": "application/julia",
"name": "julia",
"version": "1.5.4"
},
"kernelspec": {
"name": "julia-1.5",
"display_name": "Julia 1.5.4",
"language": "julia"
}
},
"nbformat": 4
}