{
"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.178874e+01 8.073477e+00\n",
" * time: 0.0019261837005615234\n",
" 1 2.979741e+01 5.177095e+00\n",
" * time: 0.008530139923095703\n",
" 2 2.136178e+01 5.751207e+00\n",
" * time: 0.01683521270751953\n",
" 3 1.751289e+01 7.333431e+00\n",
" * time: 0.02464914321899414\n",
" 4 1.192422e+01 2.676646e+00\n",
" * time: 0.03176307678222656\n",
" 5 9.798932e+00 1.021552e+00\n",
" * time: 0.03886914253234863\n",
" 6 9.152591e+00 1.245030e+00\n",
" * time: 0.04542708396911621\n",
" 7 8.787339e+00 7.303126e-01\n",
" * time: 0.05136704444885254\n",
" 8 8.488029e+00 5.894573e-01\n",
" * time: 0.05793118476867676\n",
" 9 8.298257e+00 5.768405e-01\n",
" * time: 0.063751220703125\n",
" 10 8.145146e+00 5.158911e-01\n",
" * time: 0.0696101188659668\n",
" 11 8.033811e+00 4.742253e-01\n",
" * time: 0.0761561393737793\n",
" 12 7.951600e+00 4.700619e-01\n",
" * time: 0.11073899269104004\n",
" 13 7.923141e+00 6.973155e-01\n",
" * time: 0.11396503448486328\n",
" 14 7.870060e+00 4.118104e-01\n",
" * time: 0.11706900596618652\n",
" 15 7.857083e+00 5.350892e-01\n",
" * time: 0.12032914161682129\n",
" 16 7.820183e+00 4.187180e-01\n",
" * time: 0.12432599067687988\n",
" 17 7.781929e+00 5.020283e-01\n",
" * time: 0.12743902206420898\n",
" 18 7.755426e+00 4.000637e-01\n",
" * time: 0.13132405281066895\n",
" 19 7.741935e+00 4.450722e-01\n",
" * time: 0.13444304466247559\n",
" 20 7.733792e+00 3.611611e-01\n",
" * time: 0.1376180648803711\n",
" 21 7.722208e+00 2.701878e-01\n",
" * time: 0.1407790184020996\n",
" 22 7.713792e+00 3.308856e-01\n",
" * time: 0.1438741683959961\n",
" 23 7.707877e+00 3.095860e-01\n",
" * time: 0.147003173828125\n",
" 24 7.702640e+00 2.850163e-01\n",
" * time: 0.15026307106018066\n",
" 25 7.693755e+00 2.520440e-01\n",
" * time: 0.15374016761779785\n",
" 26 7.686768e+00 2.699542e-01\n",
" * time: 0.15687012672424316\n",
" 27 7.685180e+00 2.624248e-01\n",
" * time: 0.16009020805358887\n",
" 28 7.681179e+00 1.951314e-01\n",
" * time: 0.16480302810668945\n",
" 29 7.677673e+00 2.632465e-01\n",
" * time: 0.1702280044555664\n",
" 30 7.672823e+00 2.599632e-01\n",
" * time: 0.1748790740966797\n",
" 31 7.664361e+00 1.334921e-01\n",
" * time: 0.18097305297851562\n",
" 32 7.662072e+00 2.750370e-01\n",
" * time: 0.18641018867492676\n",
" 33 7.656434e+00 1.760239e-01\n",
" * time: 0.1929171085357666\n",
" 34 7.656242e+00 2.819530e-01\n",
" * time: 0.19822907447814941\n",
" 35 7.651541e+00 1.863965e-01\n",
" * time: 0.20467400550842285\n",
" 36 7.648331e+00 2.164933e-01\n",
" * time: 0.20975208282470703\n",
" 37 7.647984e+00 2.305950e-01\n",
" * time: 0.21582913398742676\n",
" 38 7.644223e+00 1.377365e-01\n",
" * time: 0.22176098823547363\n",
" 39 7.640991e+00 1.308886e-01\n",
" * time: 0.2278749942779541\n",
" 40 7.638552e+00 1.094372e-01\n",
" * time: 0.25028109550476074\n",
" 41 7.637247e+00 1.559888e-01\n",
" * time: 0.2533860206604004\n",
" 42 7.634933e+00 1.335024e-01\n",
" * time: 0.2563471794128418\n",
" 43 7.633332e+00 1.378269e-01\n",
" * time: 0.2592802047729492\n",
" 44 7.631486e+00 1.471017e-01\n",
" * time: 0.262559175491333\n",
" 45 7.629342e+00 7.429038e-02\n",
" * time: 0.26625514030456543\n",
" 46 7.627941e+00 7.927926e-02\n",
" * time: 0.26926112174987793\n",
" 47 7.627447e+00 1.294184e-01\n",
" * time: 0.2722771167755127\n",
" 48 7.625968e+00 7.551238e-02\n",
" * time: 0.2768821716308594\n",
" 49 7.625001e+00 6.862656e-02\n",
" * time: 0.27999305725097656\n",
" 50 7.624053e+00 6.496034e-02\n",
" * time: 0.2830522060394287\n",
" 51 7.623417e+00 6.819647e-02\n",
" * time: 0.2860841751098633\n",
" 52 7.622818e+00 4.990528e-02\n",
" * time: 0.28994321823120117\n",
" 53 7.622578e+00 6.896675e-02\n",
" * time: 0.2931480407714844\n",
" 54 7.622075e+00 3.735801e-02\n",
" * time: 0.2970302104949951\n",
" 55 7.621869e+00 6.357034e-02\n",
" * time: 0.30001115798950195\n",
" 56 7.621503e+00 4.867156e-02\n",
" * time: 0.3030731678009033\n",
" 57 7.621233e+00 5.502045e-02\n",
" * time: 0.30612802505493164\n",
" 58 7.620985e+00 4.253322e-02\n",
" * time: 0.3092520236968994\n",
" 59 7.620853e+00 5.839866e-02\n",
" * time: 0.31233716011047363\n",
" 60 7.620371e+00 5.027390e-02\n",
" * time: 0.31530308723449707\n",
" 61 7.620139e+00 6.203246e-02\n",
" * time: 0.31838321685791016\n",
" 62 7.619735e+00 5.082239e-02\n",
" * time: 0.3214151859283447\n",
" 63 7.619393e+00 4.064180e-02\n",
" * time: 0.3247561454772949\n",
" 64 7.619333e+00 6.924952e-02\n",
" * time: 0.32776904106140137\n",
" 65 7.619021e+00 4.010354e-02\n",
" * time: 0.3307371139526367\n",
" 66 7.618750e+00 7.635480e-02\n",
" * time: 0.3337230682373047\n",
" 67 7.618612e+00 4.561417e-02\n",
" * time: 0.33678722381591797\n",
" 68 7.618541e+00 7.632950e-02\n",
" * time: 0.3405771255493164\n",
" 69 7.618310e+00 7.969570e-02\n",
" * time: 0.34369611740112305\n",
" 70 7.617821e+00 4.287992e-02\n",
" * time: 0.3630940914154053\n",
" 71 7.617425e+00 6.870406e-02\n",
" * time: 0.3660430908203125\n",
" 72 7.616875e+00 3.366292e-02\n",
" * time: 0.3699030876159668\n",
" 73 7.616566e+00 5.532085e-02\n",
" * time: 0.372866153717041\n",
" 74 7.616297e+00 3.974816e-02\n",
" * time: 0.3768000602722168\n",
" 75 7.616243e+00 5.849349e-02\n",
" * time: 0.37979817390441895\n",
" 76 7.616223e+00 4.705534e-02\n",
" * time: 0.3827931880950928\n",
" 77 7.616054e+00 6.111076e-02\n",
" * time: 0.3867170810699463\n",
" 78 7.615660e+00 4.233811e-02\n",
" * time: 0.39074206352233887\n",
" 79 7.615379e+00 7.147692e-02\n",
" * time: 0.3938741683959961\n",
" 80 7.614933e+00 6.416441e-02\n",
" * time: 0.3968780040740967\n",
" 81 7.614879e+00 4.140046e-02\n",
" * time: 0.39997005462646484\n",
" 82 7.614550e+00 3.417496e-02\n",
" * time: 0.4038112163543701\n",
" 83 7.614260e+00 2.930962e-02\n",
" * time: 0.40680503845214844\n",
" 84 7.614059e+00 6.696491e-02\n",
" * time: 0.4098021984100342\n",
" 85 7.613832e+00 4.517744e-02\n",
" * time: 0.41358208656311035\n",
" 86 7.613485e+00 3.696788e-02\n",
" * time: 0.4175751209259033\n",
" 87 7.613430e+00 4.164094e-02\n",
" * time: 0.4206051826477051\n",
" 88 7.613236e+00 2.950946e-02\n",
" * time: 0.4236412048339844\n",
" 89 7.613060e+00 4.014367e-02\n",
" * time: 0.42660999298095703\n",
" 90 7.612920e+00 4.530450e-02\n",
" * time: 0.4296431541442871\n",
" 91 7.612666e+00 4.512535e-02\n",
" * time: 0.4337332248687744\n",
" 92 7.612589e+00 4.888226e-02\n",
" * time: 0.43712615966796875\n",
" 93 7.612517e+00 5.051830e-02\n",
" * time: 0.4408760070800781\n",
" 94 7.612259e+00 4.780878e-02\n",
" * time: 0.4441821575164795\n",
" 95 7.612153e+00 4.427455e-02\n",
" * time: 0.4472529888153076\n",
" 96 7.612051e+00 5.523071e-02\n",
" * time: 0.4503471851348877\n",
" 97 7.611797e+00 4.794785e-02\n",
" * time: 0.45337820053100586\n",
" 98 7.611531e+00 4.616947e-02\n",
" * time: 0.4676392078399658\n",
" 99 7.611316e+00 3.163070e-02\n",
" * time: 0.4714641571044922\n",
" 100 7.611276e+00 4.972612e-02\n",
" * time: 0.47460412979125977\n",
" 101 7.611120e+00 4.902382e-02\n",
" * time: 0.4784560203552246\n",
" 102 7.610873e+00 4.003331e-02\n",
" * time: 0.48153209686279297\n",
" 103 7.610599e+00 2.962464e-02\n",
" * time: 0.485288143157959\n",
" 104 7.610364e+00 2.879457e-02\n",
" * time: 0.48934412002563477\n",
" 105 7.610209e+00 2.745396e-02\n",
" * time: 0.49321413040161133\n",
" 106 7.610100e+00 3.967563e-02\n",
" * time: 0.4966590404510498\n",
" 107 7.609961e+00 3.071004e-02\n",
" * time: 0.5007731914520264\n",
" 108 7.609774e+00 2.141246e-02\n",
" * time: 0.5046210289001465\n",
" 109 7.609687e+00 1.699303e-02\n",
" * time: 0.5085041522979736\n",
" 110 7.609656e+00 3.522981e-02\n",
" * time: 0.5128481388092041\n",
" 111 7.609654e+00 3.987945e-02\n",
" * time: 0.5158631801605225\n",
" 112 7.609501e+00 2.501769e-02\n",
" * time: 0.519017219543457\n",
" 113 7.609378e+00 3.072030e-02\n",
" * time: 0.5220541954040527\n",
" 114 7.609239e+00 3.377606e-02\n",
" * time: 0.5250520706176758\n",
" 115 7.609011e+00 2.629497e-02\n",
" * time: 0.5283660888671875\n",
" 116 7.608991e+00 4.683545e-02\n",
" * time: 0.5314991474151611\n",
" 117 7.608933e+00 3.356111e-02\n",
" * time: 0.5345752239227295\n",
" 118 7.608736e+00 2.505826e-02\n",
" * time: 0.5377521514892578\n",
" 119 7.608555e+00 3.665424e-02\n",
" * time: 0.5407431125640869\n",
" 120 7.608372e+00 2.363868e-02\n",
" * time: 0.5449399948120117\n",
" 121 7.608194e+00 2.957404e-02\n",
" * time: 0.5493760108947754\n",
" 122 7.608047e+00 1.836890e-02\n",
" * time: 0.5532190799713135\n",
" 123 7.607988e+00 2.857453e-02\n",
" * time: 0.556779146194458\n",
" 124 7.607896e+00 2.995345e-02\n",
" * time: 0.5611422061920166\n",
" 125 7.607725e+00 1.992903e-02\n",
" * time: 0.5763092041015625\n",
" 126 7.607663e+00 2.746698e-02\n",
" * time: 0.5792632102966309\n",
" 127 7.607529e+00 2.106534e-02\n",
" * time: 0.5830950736999512\n",
" 128 7.607421e+00 2.387112e-02\n",
" * time: 0.5860772132873535\n",
" 129 7.607341e+00 2.323160e-02\n",
" * time: 0.5892210006713867\n",
" 130 7.607216e+00 3.027004e-02\n",
" * time: 0.5922009944915771\n",
" 131 7.607112e+00 2.959508e-02\n",
" * time: 0.5951700210571289\n",
" 132 7.607094e+00 2.932824e-02\n",
" * time: 0.5982379913330078\n",
" 133 7.607015e+00 2.870490e-02\n",
" * time: 0.6014730930328369\n",
" 134 7.606992e+00 3.914016e-02\n",
" * time: 0.6047031879425049\n",
" 135 7.606850e+00 2.764842e-02\n",
" * time: 0.6077940464019775\n",
" 136 7.606713e+00 3.164189e-02\n",
" * time: 0.61080002784729\n",
" 137 7.606542e+00 2.572833e-02\n",
" * time: 0.6146290302276611\n",
" 138 7.606435e+00 1.820954e-02\n",
" * time: 0.6183362007141113\n",
" 139 7.606255e+00 2.055410e-02\n",
" * time: 0.6223130226135254\n",
" 140 7.606119e+00 1.972259e-02\n",
" * time: 0.6263031959533691\n",
" 141 7.606002e+00 2.063439e-02\n",
" * time: 0.6293540000915527\n",
" 142 7.605994e+00 2.788884e-02\n",
" * time: 0.6323561668395996\n",
" 143 7.605975e+00 3.637467e-02\n",
" * time: 0.6356191635131836\n",
" 144 7.605871e+00 2.875584e-02\n",
" * time: 0.6395890712738037\n",
" 145 7.605787e+00 3.737618e-02\n",
" * time: 0.6427321434020996\n",
" 146 7.605721e+00 4.622652e-02\n",
" * time: 0.6458430290222168\n",
" 147 7.605668e+00 3.754051e-02\n",
" * time: 0.649224042892456\n",
" 148 7.605544e+00 3.361410e-02\n",
" * time: 0.6524381637573242\n",
" 149 7.605355e+00 2.880915e-02\n",
" * time: 0.656400203704834\n",
" 150 7.605282e+00 3.323432e-02\n",
" * time: 0.6595320701599121\n",
" 151 7.605104e+00 3.567696e-02\n",
" * time: 0.6632671356201172\n",
" 152 7.605087e+00 4.277122e-02\n",
" * time: 0.6663320064544678\n",
" 153 7.604903e+00 3.671312e-02\n",
" * time: 0.6702370643615723\n",
" 154 7.604721e+00 3.423207e-02\n",
" * time: 0.6845030784606934\n",
" 155 7.604599e+00 3.365080e-02\n",
" * time: 0.6875340938568115\n",
" 156 7.604468e+00 2.139750e-02\n",
" * time: 0.690709114074707\n",
" 157 7.604294e+00 2.916278e-02\n",
" * time: 0.6946451663970947\n",
" 158 7.604144e+00 1.680242e-02\n",
" * time: 0.698479175567627\n",
" 159 7.604071e+00 4.157167e-02\n",
" * time: 0.7015590667724609\n",
" 160 7.603942e+00 3.413969e-02\n",
" * time: 0.7047510147094727\n",
" 161 7.603774e+00 2.456552e-02\n",
" * time: 0.7085740566253662\n",
" 162 7.603711e+00 2.305136e-02\n",
" * time: 0.712367057800293\n",
" 163 7.603599e+00 2.211205e-02\n",
" * time: 0.7171900272369385\n",
" 164 7.603512e+00 2.205009e-02\n",
" * time: 0.7213001251220703\n",
" 165 7.603482e+00 2.602874e-02\n",
" * time: 0.7245612144470215\n",
" 166 7.603381e+00 2.400506e-02\n",
" * time: 0.7282421588897705\n",
" 167 7.603311e+00 1.492612e-02\n",
" * time: 0.732370138168335\n",
" 168 7.603247e+00 2.066597e-02\n",
" * time: 0.7353332042694092\n",
" 169 7.603235e+00 2.816003e-02\n",
" * time: 0.7382621765136719\n",
" 170 7.603204e+00 2.462509e-02\n",
" * time: 0.7412362098693848\n",
" 171 7.603131e+00 2.253180e-02\n",
" * time: 0.7443389892578125\n",
" 172 7.602970e+00 2.212456e-02\n",
" * time: 0.747298002243042\n",
" 173 7.602913e+00 2.262561e-02\n",
" * time: 0.7504680156707764\n",
" 174 7.602869e+00 2.376628e-02\n",
" * time: 0.7536020278930664\n",
" 175 7.602775e+00 2.924218e-02\n",
" * time: 0.7567241191864014\n",
" 176 7.602663e+00 1.756108e-02\n",
" * time: 0.760897159576416\n",
" 177 7.602573e+00 1.902948e-02\n",
" * time: 0.7639551162719727\n",
" 178 7.602494e+00 1.892095e-02\n",
" * time: 0.7669742107391357\n",
" 179 7.602451e+00 3.367653e-02\n",
" * time: 0.7715680599212646\n",
" 180 7.602349e+00 2.319436e-02\n",
" * time: 0.775705099105835\n",
" 181 7.602341e+00 2.888374e-02\n",
" * time: 0.7899911403656006\n",
" 182 7.602234e+00 3.170918e-02\n",
" * time: 0.7933192253112793\n",
" 183 7.602180e+00 4.313864e-02\n",
" * time: 0.7963712215423584\n",
" 184 7.601940e+00 3.215992e-02\n",
" * time: 0.7999510765075684\n",
" 185 7.601717e+00 3.566516e-02\n",
" * time: 0.8038690090179443\n",
" 186 7.601505e+00 3.347708e-02\n",
" * time: 0.8069941997528076\n",
" 187 7.601498e+00 4.710457e-02\n",
" * time: 0.8099150657653809\n",
" 188 7.601187e+00 4.480785e-02\n",
" * time: 0.8130011558532715\n",
" 189 7.601025e+00 6.264988e-02\n",
" * time: 0.8160312175750732\n",
" 190 7.600704e+00 4.261948e-02\n",
" * time: 0.8197472095489502\n",
" 191 7.600438e+00 4.055995e-02\n",
" * time: 0.826714038848877\n",
" 192 7.600185e+00 4.764402e-02\n",
" * time: 0.8305211067199707\n",
" 193 7.599901e+00 6.146472e-02\n",
" * time: 0.8335700035095215\n",
" 194 7.599566e+00 3.747784e-02\n",
" * time: 0.8375611305236816\n",
" 195 7.599287e+00 3.859206e-02\n",
" * time: 0.8407211303710938\n",
" 196 7.599091e+00 3.549398e-02\n",
" * time: 0.8445870876312256\n",
" 197 7.598845e+00 3.426347e-02\n",
" * time: 0.8476331233978271\n",
" 198 7.598604e+00 2.526132e-02\n",
" * time: 0.8515481948852539\n",
" 199 7.598401e+00 2.921411e-02\n",
" * time: 0.85465407371521\n",
" 200 7.598266e+00 3.131272e-02\n",
" * time: 0.8577160835266113\n",
" 201 7.598058e+00 3.919834e-02\n",
" * time: 0.8607370853424072\n",
" 202 7.598027e+00 4.014128e-02\n",
" * time: 0.8637611865997314\n",
" 203 7.597894e+00 2.835059e-02\n",
" * time: 0.8667440414428711\n",
" 204 7.597665e+00 1.999775e-02\n",
" * time: 0.8702352046966553\n",
" 205 7.597518e+00 3.115744e-02\n",
" * time: 0.8732750415802002\n",
" 206 7.597366e+00 2.292758e-02\n",
" * time: 0.8762969970703125\n",
" 207 7.597291e+00 2.748302e-02\n",
" * time: 0.8792481422424316\n",
" 208 7.597222e+00 2.824212e-02\n",
" * time: 0.8825821876525879\n",
" 209 7.597101e+00 1.707600e-02\n",
" * time: 0.8856470584869385\n",
" 210 7.597000e+00 2.337972e-02\n",
" * time: 0.901576042175293\n",
" 211 7.596877e+00 2.414328e-02\n",
" * time: 0.9053490161895752\n",
" 212 7.596748e+00 3.249982e-02\n",
" * time: 0.9084970951080322\n",
" 213 7.596599e+00 2.574114e-02\n",
" * time: 0.9126532077789307\n",
" 214 7.596406e+00 2.198796e-02\n",
" * time: 0.9175870418548584\n",
" 215 7.596308e+00 4.039008e-02\n",
" * time: 0.9206202030181885\n",
" 216 7.596129e+00 2.585984e-02\n",
" * time: 0.9245471954345703\n",
" 217 7.596015e+00 4.297080e-02\n",
" * time: 0.927666187286377\n",
" 218 7.595744e+00 2.005980e-02\n",
" * time: 0.9334790706634521\n",
" 219 7.595662e+00 3.387969e-02\n",
" * time: 0.9367051124572754\n",
" 220 7.595448e+00 3.439935e-02\n",
" * time: 0.9401421546936035\n",
" 221 7.595436e+00 4.093395e-02\n",
" * time: 0.9435520172119141\n",
" 222 7.595265e+00 4.132066e-02\n",
" * time: 0.946641206741333\n",
" 223 7.595023e+00 2.564528e-02\n",
" * time: 0.9515261650085449\n",
" 224 7.594813e+00 3.123543e-02\n",
" * time: 0.9546411037445068\n",
" 225 7.594670e+00 3.760023e-02\n",
" * time: 0.9577770233154297\n",
" 226 7.594502e+00 3.632711e-02\n",
" * time: 0.9617199897766113\n",
" 227 7.594477e+00 4.688546e-02\n",
" * time: 0.9649491310119629\n",
" 228 7.594417e+00 3.668339e-02\n",
" * time: 0.9680631160736084\n",
" 229 7.594359e+00 2.934261e-02\n",
" * time: 0.9711260795593262\n",
" 230 7.594220e+00 3.488656e-02\n",
" * time: 0.9741899967193604\n",
" 231 7.594100e+00 2.637627e-02\n",
" * time: 0.9780139923095703\n",
" 232 7.593950e+00 1.568404e-02\n",
" * time: 0.9824151992797852\n",
" 233 7.593949e+00 2.929536e-02\n",
" * time: 0.9854280948638916\n",
" 234 7.593877e+00 2.825020e-02\n",
" * time: 0.988685131072998\n",
" 235 7.593875e+00 3.694605e-02\n",
" * time: 0.9917352199554443\n",
" 236 7.593813e+00 2.540033e-02\n",
" * time: 0.9950160980224609\n",
" 237 7.593732e+00 1.816817e-02\n",
" * time: 0.9981842041015625\n",
" 238 7.593668e+00 2.021470e-02\n",
" * time: 1.001270055770874\n",
" 239 7.593652e+00 1.959412e-02\n",
" * time: 1.0156030654907227\n",
" 240 7.593643e+00 1.528841e-02\n",
" * time: 1.018583059310913\n",
" 241 7.593592e+00 1.293105e-02\n",
" * time: 1.0214831829071045\n",
" 242 7.593561e+00 1.578187e-02\n",
" * time: 1.0244781970977783\n",
" 243 7.593536e+00 1.488465e-02\n",
" * time: 1.027576208114624\n",
" 244 7.593500e+00 1.042508e-02\n",
" * time: 1.0312600135803223\n",
" 245 7.593492e+00 1.393784e-02\n",
" * time: 1.0345280170440674\n",
" 246 7.593490e+00 1.751374e-02\n",
" * time: 1.0375101566314697\n",
" 247 7.593470e+00 1.080887e-02\n",
" * time: 1.0422492027282715\n",
" 248 7.593452e+00 1.220518e-02\n",
" * time: 1.0475430488586426\n",
" 249 7.593433e+00 1.197320e-02\n",
" * time: 1.0522961616516113\n",
" 250 7.593413e+00 1.118678e-02\n",
" * time: 1.0567421913146973\n",
" 251 7.593411e+00 1.307943e-02\n",
" * time: 1.0599250793457031\n",
" 252 7.593396e+00 1.074430e-02\n",
" * time: 1.0628890991210938\n",
" 253 7.593383e+00 1.401209e-02\n",
" * time: 1.0659661293029785\n",
" 254 7.593372e+00 1.474635e-02\n",
" * time: 1.0691230297088623\n",
" 255 7.593362e+00 1.341213e-02\n",
" * time: 1.0723280906677246\n",
" 256 7.593359e+00 2.128065e-02\n",
" * time: 1.0757040977478027\n",
" 257 7.593319e+00 1.562594e-02\n",
" * time: 1.079876184463501\n",
" 258 7.593305e+00 2.016169e-02\n",
" * time: 1.082949161529541\n",
" 259 7.593294e+00 1.591901e-02\n",
" * time: 1.0859460830688477\n",
" 260 7.593269e+00 1.251385e-02\n",
" * time: 1.089076042175293\n",
" 261 7.593259e+00 2.037200e-02\n",
" * time: 1.0921401977539062\n",
" 262 7.593213e+00 1.303774e-02\n",
" * time: 1.096013069152832\n",
" 263 7.593202e+00 1.089937e-02\n",
" * time: 1.0990991592407227\n",
" 264 7.593187e+00 1.346883e-02\n",
" * time: 1.1021602153778076\n",
" 265 7.593168e+00 1.186850e-02\n",
" * time: 1.1060290336608887\n",
" 266 7.593152e+00 1.106988e-02\n",
" * time: 1.1098110675811768\n",
" 267 7.593128e+00 5.417195e-03\n",
" * time: 1.1243290901184082\n",
" 268 7.593104e+00 7.217412e-03\n",
" * time: 1.1273550987243652\n",
" 269 7.593090e+00 1.280546e-02\n",
" * time: 1.1303582191467285\n",
" 270 7.593069e+00 8.677875e-03\n",
" * time: 1.1340100765228271\n",
" 271 7.593063e+00 1.373361e-02\n",
" * time: 1.137495994567871\n",
" 272 7.593046e+00 1.046297e-02\n",
" * time: 1.141235113143921\n",
" 273 7.593022e+00 9.802773e-03\n",
" * time: 1.1442251205444336\n",
" 274 7.593007e+00 1.508448e-02\n",
" * time: 1.1472210884094238\n",
" 275 7.593006e+00 1.019948e-02\n",
" * time: 1.1512970924377441\n",
" 276 7.593002e+00 1.179033e-02\n",
" * time: 1.155245065689087\n",
" 277 7.592979e+00 9.477377e-03\n",
" * time: 1.1585991382598877\n",
" 278 7.592963e+00 8.562695e-03\n",
" * time: 1.1624600887298584\n",
" 279 7.592952e+00 1.152652e-02\n",
" * time: 1.1655972003936768\n",
" 280 7.592939e+00 9.713536e-03\n",
" * time: 1.1688392162322998\n",
" 281 7.592921e+00 6.764368e-03\n",
" * time: 1.1719081401824951\n",
" 282 7.592906e+00 1.019102e-02\n",
" * time: 1.1749610900878906\n",
" 283 7.592896e+00 6.175085e-03\n",
" * time: 1.1778922080993652\n",
" 284 7.592890e+00 1.363663e-02\n",
" * time: 1.180891990661621\n",
" 285 7.592879e+00 9.814314e-03\n",
" * time: 1.1841211318969727\n",
" 286 7.592865e+00 8.229212e-03\n",
" * time: 1.187880039215088\n",
" 287 7.592847e+00 5.698404e-03\n",
" * time: 1.1916141510009766\n",
" 288 7.592838e+00 8.633349e-03\n",
" * time: 1.194653034210205\n",
" 289 7.592823e+00 8.669478e-03\n",
" * time: 1.1977641582489014\n",
" 290 7.592800e+00 6.695463e-03\n",
" * time: 1.2015750408172607\n",
" 291 7.592790e+00 7.781029e-03\n",
" * time: 1.2046151161193848\n",
" 292 7.592776e+00 6.837265e-03\n",
" * time: 1.2076141834259033\n",
" 293 7.592774e+00 9.142954e-03\n",
" * time: 1.2106280326843262\n",
" 294 7.592760e+00 7.336330e-03\n",
" * time: 1.214637041091919\n",
" 295 7.592752e+00 5.954362e-03\n",
" * time: 1.2184321880340576\n",
" 296 7.592744e+00 7.816504e-03\n",
" * time: 1.232335090637207\n",
" 297 7.592738e+00 6.038311e-03\n",
" * time: 1.235262155532837\n",
" 298 7.592731e+00 9.058387e-03\n",
" * time: 1.2382631301879883\n",
" 299 7.592718e+00 4.804684e-03\n",
" * time: 1.241927146911621\n",
" 300 7.592704e+00 5.059608e-03\n",
" * time: 1.2451519966125488\n",
" 301 7.592693e+00 6.592889e-03\n",
" * time: 1.2483441829681396\n",
" 302 7.592687e+00 9.547255e-03\n",
" * time: 1.2514231204986572\n",
" 303 7.592677e+00 9.878041e-03\n",
" * time: 1.2551991939544678\n",
" 304 7.592662e+00 6.916640e-03\n",
" * time: 1.2589170932769775\n",
" 305 7.592659e+00 1.265315e-02\n",
" * time: 1.2630221843719482\n",
" 306 7.592654e+00 1.058305e-02\n",
" * time: 1.2661221027374268\n",
" 307 7.592629e+00 9.488867e-03\n",
" * time: 1.269336223602295\n",
" 308 7.592610e+00 6.029180e-03\n",
" * time: 1.27311110496521\n",
" 309 7.592598e+00 1.202902e-02\n",
" * time: 1.277940034866333\n",
" 310 7.592590e+00 8.424549e-03\n",
" * time: 1.2809841632843018\n",
" 311 7.592571e+00 8.635251e-03\n",
" * time: 1.2839791774749756\n",
" 312 7.592571e+00 1.616907e-02\n",
" * time: 1.2869141101837158\n",
" 313 7.592555e+00 1.177452e-02\n",
" * time: 1.2906970977783203\n",
" 314 7.592545e+00 1.191066e-02\n",
" * time: 1.293992042541504\n",
" 315 7.592543e+00 1.376595e-02\n",
" * time: 1.2969729900360107\n",
" 316 7.592536e+00 1.861995e-02\n",
" * time: 1.2999351024627686\n",
" 317 7.592517e+00 1.171532e-02\n",
" * time: 1.3036072254180908\n",
" 318 7.592504e+00 1.067271e-02\n",
" * time: 1.310554027557373\n",
" 319 7.592483e+00 1.106798e-02\n",
" * time: 1.3136861324310303\n",
" 320 7.592470e+00 7.090305e-03\n",
" * time: 1.3177011013031006\n",
" 321 7.592458e+00 5.308948e-03\n",
" * time: 1.321681022644043\n",
" 322 7.592449e+00 5.219042e-03\n",
" * time: 1.3247640132904053\n",
" 323 7.592441e+00 3.782524e-03\n",
" * time: 1.328613042831421\n",
" 324 7.592440e+00 7.543645e-03\n",
" * time: 1.343257188796997\n",
" 325 7.592434e+00 5.702459e-03\n",
" * time: 1.347254991531372\n",
" 326 7.592428e+00 6.136651e-03\n",
" * time: 1.3511149883270264\n",
" 327 7.592426e+00 8.060949e-03\n",
" * time: 1.3544120788574219\n",
" 328 7.592422e+00 6.809108e-03\n",
" * time: 1.357713222503662\n",
" 329 7.592414e+00 5.171380e-03\n",
" * time: 1.3608860969543457\n",
" 330 7.592407e+00 4.092982e-03\n",
" * time: 1.3642570972442627\n",
" 331 7.592402e+00 4.660772e-03\n",
" * time: 1.3672871589660645\n",
" 332 7.592402e+00 6.164150e-03\n",
" * time: 1.3705930709838867\n",
" 333 7.592398e+00 4.531960e-03\n",
" * time: 1.3737151622772217\n",
" 334 7.592396e+00 4.721444e-03\n",
" * time: 1.3767940998077393\n",
" 335 7.592394e+00 6.525967e-03\n",
" * time: 1.3799610137939453\n",
" 336 7.592390e+00 6.314100e-03\n",
" * time: 1.3829681873321533\n",
" 337 7.592389e+00 6.543146e-03\n",
" * time: 1.386641025543213\n",
" 338 7.592384e+00 4.495488e-03\n",
" * time: 1.3905880451202393\n",
" 339 7.592382e+00 4.804942e-03\n",
" * time: 1.3936800956726074\n",
" 340 7.592381e+00 4.764296e-03\n",
" * time: 1.3967361450195312\n",
" 341 7.592377e+00 3.627922e-03\n",
" * time: 1.399705171585083\n",
" 342 7.592376e+00 5.834643e-03\n",
" * time: 1.4028542041778564\n",
" 343 7.592373e+00 3.985008e-03\n",
" * time: 1.4067580699920654\n",
" 344 7.592370e+00 4.312561e-03\n",
" * time: 1.4105451107025146\n",
" 345 7.592368e+00 2.515741e-03\n",
" * time: 1.4142770767211914\n",
" 346 7.592368e+00 4.809344e-03\n",
" * time: 1.4179151058197021\n",
" 347 7.592367e+00 4.686798e-03\n",
" * time: 1.421138048171997\n",
" 348 7.592364e+00 2.976821e-03\n",
" * time: 1.425159215927124\n",
" 349 7.592363e+00 3.593114e-03\n",
" * time: 1.4282660484313965\n",
" 350 7.592360e+00 2.649794e-03\n",
" * time: 1.432314157485962\n",
" 351 7.592360e+00 3.578413e-03\n",
" * time: 1.435420036315918\n",
" 352 7.592358e+00 4.016736e-03\n",
" * time: 1.44926118850708\n",
" 353 7.592357e+00 3.370890e-03\n",
" * time: 1.4522972106933594\n",
" 354 7.592356e+00 2.969040e-03\n",
" * time: 1.4554600715637207\n",
" 355 7.592356e+00 3.496307e-03\n",
" * time: 1.4586620330810547\n",
" 356 7.592355e+00 3.216857e-03\n",
" * time: 1.461724042892456\n",
" 357 7.592354e+00 1.794594e-03\n",
" * time: 1.4658520221710205\n",
" 358 7.592353e+00 2.042724e-03\n",
" * time: 1.4689981937408447\n",
" 359 7.592352e+00 1.902743e-03\n",
" * time: 1.4723050594329834\n",
" 360 7.592352e+00 2.324789e-03\n",
" * time: 1.4754149913787842\n",
" 361 7.592351e+00 1.705313e-03\n",
" * time: 1.4801630973815918\n",
" 362 7.592351e+00 9.396563e-04\n",
" * time: 1.4841279983520508\n",
" 363 7.592350e+00 1.183471e-03\n",
" * time: 1.4872591495513916\n",
" 364 7.592350e+00 8.055184e-04\n",
" * time: 1.4912211894989014\n",
" 365 7.592350e+00 1.295509e-03\n",
" * time: 1.4947261810302734\n",
" 366 7.592350e+00 8.265826e-04\n",
" * time: 1.498716115951538\n",
" 367 7.592350e+00 1.523596e-03\n",
" * time: 1.5018670558929443\n",
" 368 7.592349e+00 1.144253e-03\n",
" * time: 1.5054430961608887\n",
" 369 7.592349e+00 1.515984e-03\n",
" * time: 1.5084760189056396\n",
" 370 7.592349e+00 9.582551e-04\n",
" * time: 1.5121300220489502\n",
" 371 7.592349e+00 8.825099e-04\n",
" * time: 1.515334129333496\n",
" 372 7.592349e+00 5.697221e-04\n",
" * time: 1.518799066543579\n",
" 373 7.592349e+00 8.375615e-04\n",
" * time: 1.522038221359253\n",
" 374 7.592349e+00 1.304432e-03\n",
" * time: 1.5251061916351318\n",
" 375 7.592348e+00 7.301717e-04\n",
" * time: 1.5292251110076904\n",
" 376 7.592348e+00 1.200565e-03\n",
" * time: 1.5325210094451904\n",
" 377 7.592348e+00 8.754370e-04\n",
" * time: 1.535696029663086\n",
" 378 7.592348e+00 8.733190e-04\n",
" * time: 1.5399761199951172\n",
" 379 7.592348e+00 9.046054e-04\n",
" * time: 1.543424129486084\n",
" 380 7.592348e+00 1.287683e-03\n",
" * time: 1.5467371940612793\n",
" 381 7.592348e+00 9.302263e-04\n",
" * time: 1.5653040409088135\n",
" 382 7.592348e+00 7.341308e-04\n",
" * time: 1.5692510604858398\n",
" 383 7.592348e+00 8.416323e-04\n",
" * time: 1.5724570751190186\n",
" 384 7.592348e+00 1.021006e-03\n",
" * time: 1.5755162239074707\n",
" 385 7.592348e+00 1.131927e-03\n",
" * time: 1.5787580013275146\n",
" 386 7.592347e+00 1.004103e-03\n",
" * time: 1.5818841457366943\n",
" 387 7.592347e+00 6.921731e-04\n",
" * time: 1.586055040359497\n",
" 388 7.592347e+00 6.374985e-04\n",
" * time: 1.5895881652832031\n",
" 389 7.592347e+00 7.999362e-04\n",
" * time: 1.5928301811218262\n",
" 390 7.592347e+00 6.665308e-04\n",
" * time: 1.5969011783599854\n",
" 391 7.592347e+00 7.635281e-04\n",
" * time: 1.6007680892944336\n",
" 392 7.592346e+00 4.271714e-04\n",
" * time: 1.605147123336792\n",
" 393 7.592346e+00 6.001808e-04\n",
" * time: 1.6084561347961426\n",
" 394 7.592346e+00 9.803160e-04\n",
" * time: 1.6119441986083984\n",
" 395 7.592346e+00 1.112591e-03\n",
" * time: 1.615922212600708\n",
" 396 7.592346e+00 1.677715e-03\n",
" * time: 1.619277000427246\n",
" 397 7.592346e+00 1.795487e-03\n",
" * time: 1.6230721473693848\n",
" 398 7.592346e+00 1.217728e-03\n",
" * time: 1.6261351108551025\n",
" 399 7.592345e+00 8.018916e-04\n",
" * time: 1.6292321681976318\n",
" 400 7.592345e+00 5.000411e-04\n",
" * time: 1.6331980228424072\n",
" 401 7.592345e+00 5.054228e-04\n",
" * time: 1.6364622116088867\n",
" 402 7.592345e+00 8.793834e-04\n",
" * time: 1.63950514793396\n",
" 403 7.592345e+00 1.104538e-03\n",
" * time: 1.6425971984863281\n",
" 404 7.592345e+00 9.294481e-04\n",
" * time: 1.645719051361084\n",
" 405 7.592345e+00 7.831222e-04\n",
" * time: 1.6487882137298584\n",
" 406 7.592345e+00 8.385505e-04\n",
" * time: 1.652095079421997\n",
" 407 7.592344e+00 4.870714e-04\n",
" * time: 1.6560101509094238\n",
" 408 7.592344e+00 7.425085e-04\n",
" * time: 1.659121036529541\n",
" 409 7.592344e+00 7.580187e-04\n",
" * time: 1.662142038345337\n",
" 410 7.592344e+00 5.659748e-04\n",
" * time: 1.676516056060791\n",
" 411 7.592344e+00 4.802877e-04\n",
" * time: 1.6796081066131592\n",
" 412 7.592344e+00 5.594980e-04\n",
" * time: 1.683091163635254\n",
" 413 7.592344e+00 3.749001e-04\n",
" * time: 1.6869161128997803\n",
" 414 7.592344e+00 4.705403e-04\n",
" * time: 1.690737009048462\n",
" 415 7.592344e+00 3.800642e-04\n",
" * time: 1.6945600509643555\n",
" 416 7.592344e+00 6.702253e-04\n",
" * time: 1.6985750198364258\n",
" 417 7.592344e+00 8.769411e-04\n",
" * time: 1.7017850875854492\n",
" 418 7.592344e+00 7.031342e-04\n",
" * time: 1.7053380012512207\n",
" 419 7.592344e+00 1.062697e-03\n",
" * time: 1.708470106124878\n",
" 420 7.592344e+00 9.701545e-04\n",
" * time: 1.7115070819854736\n",
" 421 7.592344e+00 6.826293e-04\n",
" * time: 1.714691162109375\n",
" 422 7.592343e+00 9.183249e-04\n",
" * time: 1.7178261280059814\n",
" 423 7.592343e+00 1.384686e-03\n",
" * time: 1.7209692001342773\n",
" 424 7.592343e+00 1.018084e-03\n",
" * time: 1.724133014678955\n",
" 425 7.592343e+00 9.877895e-04\n",
" * time: 1.7273261547088623\n",
" 426 7.592343e+00 9.544149e-04\n",
" * time: 1.731236219406128\n",
" 427 7.592343e+00 6.709234e-04\n",
" * time: 1.735044002532959\n",
" 428 7.592342e+00 5.551956e-04\n",
" * time: 1.7389750480651855\n",
" 429 7.592342e+00 1.294225e-03\n",
" * time: 1.7420220375061035\n",
" 430 7.592342e+00 1.048887e-03\n",
" * time: 1.7459421157836914\n",
" 431 7.592342e+00 6.578802e-04\n",
" * time: 1.7496991157531738\n",
" 432 7.592342e+00 9.223837e-04\n",
" * time: 1.752655029296875\n",
" 433 7.592342e+00 1.299705e-03\n",
" * time: 1.7556390762329102\n",
" 434 7.592341e+00 1.420530e-03\n",
" * time: 1.7585911750793457\n",
" 435 7.592341e+00 8.023613e-04\n",
" * time: 1.7624812126159668\n",
" 436 7.592341e+00 1.098398e-03\n",
" * time: 1.7654402256011963\n",
" 437 7.592340e+00 1.258763e-03\n",
" * time: 1.7684111595153809\n",
" 438 7.592340e+00 1.745657e-03\n",
" * time: 1.7815861701965332\n",
" 439 7.592339e+00 2.592376e-03\n",
" * time: 1.784811019897461\n",
" 440 7.592339e+00 2.049833e-03\n",
" * time: 1.7892119884490967\n",
" 441 7.592338e+00 1.398813e-03\n",
" * time: 1.7933721542358398\n",
" 442 7.592338e+00 1.510427e-03\n",
" * time: 1.7969591617584229\n",
" 443 7.592337e+00 1.240393e-03\n",
" * time: 1.8001160621643066\n",
" 444 7.592336e+00 1.602733e-03\n",
" * time: 1.8035390377044678\n",
" 445 7.592336e+00 2.628793e-03\n",
" * time: 1.8068370819091797\n",
" 446 7.592336e+00 2.290619e-03\n",
" * time: 1.8100759983062744\n",
" 447 7.592336e+00 3.092115e-03\n",
" * time: 1.8133020401000977\n",
" 448 7.592335e+00 2.249988e-03\n",
" * time: 1.8164551258087158\n",
" 449 7.592335e+00 2.264389e-03\n",
" * time: 1.8195490837097168\n",
" 450 7.592334e+00 1.794218e-03\n",
" * time: 1.8235549926757812\n",
" 451 7.592334e+00 1.971992e-03\n",
" * time: 1.826700210571289\n",
" 452 7.592333e+00 2.026081e-03\n",
" * time: 1.8297932147979736\n",
" 453 7.592333e+00 1.377200e-03\n",
" * time: 1.8328380584716797\n",
" 454 7.592332e+00 1.749173e-03\n",
" * time: 1.8358771800994873\n",
" 455 7.592332e+00 1.377882e-03\n",
" * time: 1.8388922214508057\n",
" 456 7.592332e+00 7.152823e-04\n",
" * time: 1.842616081237793\n",
" 457 7.592332e+00 1.173806e-03\n",
" * time: 1.8455910682678223\n",
" 458 7.592331e+00 1.009981e-03\n",
" * time: 1.8492870330810547\n",
" 459 7.592331e+00 8.126221e-04\n",
" * time: 1.8530230522155762\n",
" 460 7.592331e+00 1.061977e-03\n",
" * time: 1.8577570915222168\n",
" 461 7.592331e+00 1.241536e-03\n",
" * time: 1.8607861995697021\n",
" 462 7.592331e+00 9.277153e-04\n",
" * time: 1.863821029663086\n",
" 463 7.592330e+00 1.313096e-03\n",
" * time: 1.8667922019958496\n",
" 464 7.592330e+00 9.075196e-04\n",
" * time: 1.8699851036071777\n",
" 465 7.592330e+00 1.518819e-03\n",
" * time: 1.8731040954589844\n",
" 466 7.592330e+00 1.146519e-03\n",
" * time: 1.8763470649719238\n",
" 467 7.592330e+00 8.918493e-04\n",
" * time: 1.892902135848999\n",
" 468 7.592329e+00 1.495819e-03\n",
" * time: 1.8960530757904053\n",
" 469 7.592329e+00 1.364271e-03\n",
" * time: 1.8991100788116455\n",
" 470 7.592329e+00 1.666263e-03\n",
" * time: 1.9022080898284912\n",
" 471 7.592329e+00 8.560064e-04\n",
" * time: 1.9059371948242188\n",
" 472 7.592328e+00 1.334364e-03\n",
" * time: 1.9090960025787354\n",
" 473 7.592328e+00 8.781824e-04\n",
" * time: 1.9129881858825684\n",
" 474 7.592328e+00 2.185644e-03\n",
" * time: 1.9190821647644043\n",
" 475 7.592327e+00 1.660704e-03\n",
" * time: 1.9230351448059082\n",
" 476 7.592327e+00 1.607191e-03\n",
" * time: 1.9269700050354004\n",
" 477 7.592327e+00 8.625173e-04\n",
" * time: 1.9309871196746826\n",
" 478 7.592327e+00 1.232970e-03\n",
" * time: 1.9341351985931396\n",
" 479 7.592327e+00 1.234049e-03\n",
" * time: 1.9371230602264404\n",
" 480 7.592326e+00 1.704748e-03\n",
" * time: 1.9403181076049805\n",
" 481 7.592326e+00 1.197144e-03\n",
" * time: 1.9435491561889648\n",
" 482 7.592326e+00 1.306136e-03\n",
" * time: 1.9465701580047607\n",
" 483 7.592326e+00 1.878991e-03\n",
" * time: 1.950618028640747\n",
" 484 7.592326e+00 2.050965e-03\n",
" * time: 1.9536521434783936\n",
" 485 7.592325e+00 1.025332e-03\n",
" * time: 1.9574170112609863\n",
" 486 7.592325e+00 1.433895e-03\n",
" * time: 1.9603941440582275\n",
" 487 7.592325e+00 9.639850e-04\n",
" * time: 1.9634981155395508\n",
" 488 7.592325e+00 6.557239e-04\n",
" * time: 1.9666571617126465\n",
" 489 7.592325e+00 1.165949e-03\n",
" * time: 1.9697670936584473\n",
" 490 7.592325e+00 1.391569e-03\n",
" * time: 1.9727270603179932\n",
" 491 7.592324e+00 8.816796e-04\n",
" * time: 1.976546049118042\n",
" 492 7.592324e+00 1.249776e-03\n",
" * time: 1.9799981117248535\n",
" 493 7.592324e+00 1.093522e-03\n",
" * time: 1.982945203781128\n",
" 494 7.592324e+00 1.115956e-03\n",
" * time: 1.9868900775909424\n",
" 495 7.592324e+00 5.608803e-04\n",
" * time: 2.003369092941284\n",
" 496 7.592324e+00 7.588004e-04\n",
" * time: 2.0064401626586914\n",
" 497 7.592324e+00 1.148327e-03\n",
" * time: 2.0094141960144043\n",
" 498 7.592323e+00 7.867448e-04\n",
" * time: 2.0133349895477295\n",
" 499 7.592323e+00 9.756409e-04\n",
" * time: 2.016308069229126\n",
" 500 7.592323e+00 1.103842e-03\n",
" * time: 2.0194032192230225\n",
" 501 7.592323e+00 1.154035e-03\n",
" * time: 2.0224151611328125\n",
" 502 7.592323e+00 1.370253e-03\n",
" * time: 2.033802032470703\n",
" 503 7.592323e+00 1.519947e-03\n",
" * time: 2.036867141723633\n",
" 504 7.592323e+00 1.285826e-03\n",
" * time: 2.039914131164551\n",
" 505 7.592323e+00 1.016815e-03\n",
" * time: 2.0432209968566895\n",
" 506 7.592322e+00 1.322440e-03\n",
" * time: 2.046541213989258\n",
" 507 7.592322e+00 1.036879e-03\n",
" * time: 2.049751043319702\n",
" 508 7.592322e+00 9.490958e-04\n",
" * time: 2.0536410808563232\n",
" 509 7.592322e+00 9.219205e-04\n",
" * time: 2.056910991668701\n",
" 510 7.592322e+00 1.101652e-03\n",
" * time: 2.0604372024536133\n",
" 511 7.592321e+00 7.473969e-04\n",
" * time: 2.064846992492676\n",
" 512 7.592321e+00 1.293986e-03\n",
" * time: 2.068303108215332\n",
" 513 7.592321e+00 1.748231e-03\n",
" * time: 2.071533203125\n",
" 514 7.592321e+00 1.355053e-03\n",
" * time: 2.0755200386047363\n",
" 515 7.592321e+00 9.725563e-04\n",
" * time: 2.0792620182037354\n",
" 516 7.592320e+00 1.136404e-03\n",
" * time: 2.082357168197632\n",
" 517 7.592320e+00 1.850260e-03\n",
" * time: 2.0855281352996826\n",
" 518 7.592320e+00 9.532045e-04\n",
" * time: 2.0888261795043945\n",
" 519 7.592320e+00 9.980898e-04\n",
" * time: 2.092229127883911\n",
" 520 7.592319e+00 1.922664e-03\n",
" * time: 2.095449209213257\n",
" 521 7.592319e+00 1.152683e-03\n",
" * time: 2.0992350578308105\n",
" 522 7.592319e+00 1.528200e-03\n",
" * time: 2.102262020111084\n",
" 523 7.592318e+00 1.496721e-03\n",
" * time: 2.106304168701172\n",
" 524 7.592318e+00 9.032715e-04\n",
" * time: 2.121448040008545\n",
" 525 7.592318e+00 1.200588e-03\n",
" * time: 2.124624013900757\n",
" 526 7.592318e+00 1.744404e-03\n",
" * time: 2.127732992172241\n",
" 527 7.592317e+00 1.267899e-03\n",
" * time: 2.130690097808838\n",
" 528 7.592317e+00 2.444900e-03\n",
" * time: 2.1338090896606445\n",
" 529 7.592317e+00 2.053855e-03\n",
" * time: 2.1369872093200684\n",
" 530 7.592317e+00 1.484579e-03\n",
" * time: 2.140122175216675\n",
" 531 7.592316e+00 2.667103e-03\n",
" * time: 2.143104076385498\n",
" 532 7.592316e+00 2.056917e-03\n",
" * time: 2.1461751461029053\n",
" 533 7.592315e+00 1.032773e-03\n",
" * time: 2.1500470638275146\n",
" 534 7.592315e+00 1.218417e-03\n",
" * time: 2.154355049133301\n",
" 535 7.592315e+00 1.653209e-03\n",
" * time: 2.157539129257202\n",
" 536 7.592314e+00 1.982212e-03\n",
" * time: 2.1605942249298096\n",
" 537 7.592314e+00 9.748209e-04\n",
" * time: 2.164414167404175\n",
" 538 7.592314e+00 7.655541e-04\n",
" * time: 2.1685280799865723\n",
" 539 7.592314e+00 1.742880e-03\n",
" * time: 2.1716132164001465\n",
" 540 7.592314e+00 1.480555e-03\n",
" * time: 2.1755261421203613\n",
" 541 7.592313e+00 1.715307e-03\n",
" * time: 2.1787350177764893\n",
" 542 7.592313e+00 2.342237e-03\n",
" * time: 2.181934118270874\n",
" 543 7.592313e+00 1.529450e-03\n",
" * time: 2.18495512008667\n",
" 544 7.592312e+00 1.789278e-03\n",
" * time: 2.188035011291504\n",
" 545 7.592312e+00 1.992743e-03\n",
" * time: 2.191035032272339\n",
" 546 7.592312e+00 3.251821e-03\n",
" * time: 2.1941351890563965\n",
" 547 7.592311e+00 3.186465e-03\n",
" * time: 2.1973211765289307\n",
" 548 7.592311e+00 3.430129e-03\n",
" * time: 2.2003660202026367\n",
" 549 7.592311e+00 2.919029e-03\n",
" * time: 2.2035059928894043\n",
" 550 7.592310e+00 3.061003e-03\n",
" * time: 2.2073891162872314\n",
" 551 7.592309e+00 2.923366e-03\n",
" * time: 2.2104811668395996\n",
" 552 7.592308e+00 2.602203e-03\n",
" * time: 2.2136940956115723\n",
" 553 7.592307e+00 1.994923e-03\n",
" * time: 2.2288641929626465\n",
" 554 7.592307e+00 1.267185e-03\n",
" * time: 2.232771158218384\n",
" 555 7.592306e+00 1.393623e-03\n",
" * time: 2.2359001636505127\n",
" 556 7.592306e+00 1.968744e-03\n",
" * time: 2.238938093185425\n",
" 557 7.592306e+00 1.485990e-03\n",
" * time: 2.2427830696105957\n",
" 558 7.592305e+00 1.337964e-03\n",
" * time: 2.2468771934509277\n",
" 559 7.592304e+00 1.108946e-03\n",
" * time: 2.2506980895996094\n",
" 560 7.592304e+00 1.494092e-03\n",
" * time: 2.2538201808929443\n",
" 561 7.592304e+00 1.405447e-03\n",
" * time: 2.257740020751953\n",
" 562 7.592304e+00 2.261017e-03\n",
" * time: 2.2619781494140625\n",
" 563 7.592303e+00 1.533952e-03\n",
" * time: 2.265944004058838\n",
" 564 7.592303e+00 1.355182e-03\n",
" * time: 2.269846200942993\n",
" 565 7.592303e+00 9.969089e-04\n",
" * time: 2.272935152053833\n",
" 566 7.592303e+00 1.185422e-03\n",
" * time: 2.276196002960205\n",
" 567 7.592302e+00 7.337092e-04\n",
" * time: 2.279283046722412\n",
" 568 7.592302e+00 8.200066e-04\n",
" * time: 2.2824110984802246\n",
" 569 7.592302e+00 1.215014e-03\n",
" * time: 2.2854421138763428\n",
" 570 7.592302e+00 1.158982e-03\n",
" * time: 2.2885310649871826\n",
" 571 7.592302e+00 1.155542e-03\n",
" * time: 2.2920000553131104\n",
" 572 7.592302e+00 1.092271e-03\n",
" * time: 2.2950611114501953\n",
" 573 7.592301e+00 7.106530e-04\n",
" * time: 2.2989280223846436\n",
" 574 7.592301e+00 8.167775e-04\n",
" * time: 2.3021349906921387\n",
" 575 7.592301e+00 6.901315e-04\n",
" * time: 2.305183172225952\n",
" 576 7.592301e+00 3.683400e-04\n",
" * time: 2.309218168258667\n",
" 577 7.592301e+00 4.379561e-04\n",
" * time: 2.312296152114868\n",
" 578 7.592301e+00 5.364424e-04\n",
" * time: 2.3152620792388916\n",
" 579 7.592301e+00 7.207831e-04\n",
" * time: 2.3183610439300537\n",
" 580 7.592301e+00 7.871642e-04\n",
" * time: 2.321453094482422\n",
" 581 7.592301e+00 8.199022e-04\n",
" * time: 2.3352301120758057\n",
" 582 7.592301e+00 6.618446e-04\n",
" * time: 2.338654041290283\n",
" 583 7.592301e+00 3.338160e-04\n",
" * time: 2.3424901962280273\n",
" 584 7.592301e+00 4.852660e-04\n",
" * time: 2.3455212116241455\n",
" 585 7.592301e+00 6.237890e-04\n",
" * time: 2.3485472202301025\n",
" 586 7.592301e+00 5.299105e-04\n",
" * time: 2.351516008377075\n",
" 587 7.592301e+00 4.818951e-04\n",
" * time: 2.354936122894287\n",
" 588 7.592301e+00 5.012142e-04\n",
" * time: 2.357990026473999\n",
" 589 7.592301e+00 3.445064e-04\n",
" * time: 2.361138105392456\n",
" 590 7.592301e+00 7.206400e-04\n",
" * time: 2.364143133163452\n",
" 591 7.592301e+00 4.871841e-04\n",
" * time: 2.3679890632629395\n",
" 592 7.592300e+00 3.987234e-04\n",
" * time: 2.3804190158843994\n",
" 593 7.592300e+00 5.620618e-04\n",
" * time: 2.3836100101470947\n",
" 594 7.592300e+00 5.805472e-04\n",
" * time: 2.3868260383605957\n",
" 595 7.592300e+00 2.927758e-04\n",
" * time: 2.3906571865081787\n",
" 596 7.592300e+00 3.763173e-04\n",
" * time: 2.394559144973755\n",
" 597 7.592300e+00 2.608761e-04\n",
" * time: 2.398339033126831\n",
" 598 7.592300e+00 2.409925e-04\n",
" * time: 2.4024240970611572\n",
" 599 7.592300e+00 1.766842e-04\n",
" * time: 2.406294107437134\n",
" 600 7.592300e+00 2.981309e-04\n",
" * time: 2.4093780517578125\n",
" 601 7.592300e+00 3.753438e-04\n",
" * time: 2.412440061569214\n",
" 602 7.592300e+00 2.961948e-04\n",
" * time: 2.415574073791504\n",
" 603 7.592300e+00 4.956437e-04\n",
" * time: 2.418837070465088\n",
" 604 7.592300e+00 4.682776e-04\n",
" * time: 2.4219861030578613\n",
" 605 7.592300e+00 4.215651e-04\n",
" * time: 2.425140142440796\n",
" 606 7.592300e+00 4.610572e-04\n",
" * time: 2.4281561374664307\n",
" 607 7.592300e+00 2.893157e-04\n",
" * time: 2.431274175643921\n",
" 608 7.592300e+00 3.654500e-04\n",
" * time: 2.434512138366699\n",
" 609 7.592300e+00 3.085864e-04\n",
" * time: 2.4376330375671387\n",
" 610 7.592300e+00 3.597196e-04\n",
" * time: 2.451763153076172\n",
" 611 7.592300e+00 3.208603e-04\n",
" * time: 2.455587148666382\n",
" 612 7.592300e+00 2.892662e-04\n",
" * time: 2.459484100341797\n",
" 613 7.592300e+00 3.864581e-04\n",
" * time: 2.4625661373138428\n",
" 614 7.592300e+00 3.362437e-04\n",
" * time: 2.467254161834717\n",
" 615 7.592300e+00 1.801129e-04\n",
" * time: 2.4712531566619873\n",
" 616 7.592300e+00 1.580500e-04\n",
" * time: 2.475097179412842\n",
" 617 7.592300e+00 4.511348e-04\n",
" * time: 2.478142023086548\n",
" 618 7.592300e+00 2.609911e-04\n",
" * time: 2.4815640449523926\n",
" 619 7.592300e+00 1.558484e-04\n",
" * time: 2.4855270385742188\n",
" 620 7.592300e+00 2.731299e-04\n",
" * time: 2.488862991333008\n",
" 621 7.592300e+00 1.139913e-04\n",
" * time: 2.492964029312134\n",
" 622 7.592300e+00 1.752264e-04\n",
" * time: 2.4960880279541016\n",
" 623 7.592300e+00 3.147842e-04\n",
" * time: 2.4991111755371094\n",
" 624 7.592300e+00 2.427371e-04\n",
" * time: 2.502286195755005\n",
" 625 7.592300e+00 1.822185e-04\n",
" * time: 2.5053250789642334\n",
" 626 7.592300e+00 2.102145e-04\n",
" * time: 2.5083601474761963\n",
" 627 7.592300e+00 1.884880e-04\n",
" * time: 2.511660099029541\n",
" 628 7.592300e+00 1.921035e-04\n",
" * time: 2.515496015548706\n"
]
},
{
"output_type": "execute_result",
"data": {
"text/plain": "Plot{Plots.GRBackend() n=1}",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3dd3xUVfo/8DM9vZAEUiDFAAFCJ0GKQAihgzSRKEQR/aGCui4oZUEEdAEBKyqiGzGrgqBfpbfQVlRqBAJIMaGF9N4zmfb747KzEc8TuOMNN5P5vF/8MTlz5+ZM4+Tc+7nPUVgsFgYAAOColHJ3AAAAQE4YCAEAwKFhIAQAAIeGgRAAABwaBkIAAHBoGAgBAMChYSAEAACHhoEQAAAcGgZCAABwaBgIAQDAoan/4uOvpKV7+Po5OztL0ptGyGw2KxQKhUIhd0fuB4vFYrFYlEpH+fPIbDbjyTZJeLINxFOnues2ly8XVVcbRe22oqIsKEgVFhZma7/+KsVfrDU6efqMCbNfk6o3AADQaI2PCLjrNl27Jp09my9yx+fHjCnZsmWLbb366xzljyYAAACuv3poFAAAwEqhYHZ3KgkzQgAAcGiYEQIAgGRsSBdaLDJPITEQAgCApHBoFAAAwI5gRggAANKxw7AMBkIAAJCMLRVI5D5HiEOjAADg0DAQAgCAQ8OhUQAAkIwNF9TLfWQUAyEAAEjHhnOEsodrMBACAECjZrFYUlJSysvLe/bs6erqese9aWlp169fj4yMDAi4XRM8KyvrypUr7u7uXbp0UatvD3MZGRmXL19u3bp1aGio9bGpqan5+fk4RwgAANJRiP9XL5PJ9PDDDyckJKxYsaJt27aXL1+ue++yZcseeuihDz74oHPnzt999x1j7OWXX+7ates///nPhISEzp07Z2VlMcbWr1/frVu3NWvW9OzZ86OPPhIe+8wzz4wZM+add97BQAgAAJJR/Pc0oah/9di+ffvly5dTUlL27t07efLkxYsXW+/Kycl58803f/zxx23btq1fv37WrFkmkykhISEjIyM5Ofn8+fNhYWErVqyorq5+5ZVXfvjhh61bt+7Zs2fevHllZWUnTpzYunXrqVOndu7ciUOjAADQeH3//fcTJkxwcXFhjE2ZMqVnz57WtYi3b9/evXv3tm3bMsZGjBhRWVl54sSJ3r17Cw9UKpUdOnTIyck5fPiwu7t7v379GGPdu3cPCQnZu3dvSkrKiBEjfHx8GC6fAAAACSnEq39KmJGRERISItwOCQnR6/X5+bcX/r1165b1LqVS2bJly1u3blkfmJeX9/XXXz/22GN198AYCw4OvnXr1q1bt6wnCzEjBAAA6dzDab8/P+TatWvWU3eCSZMm+fr6Msb0er1GoxEatVotY6y6ulr4saamxpqFYYzpdDrrXVVVVY888sj48eNHjBixZs2aP29W97EYCAEAQGbl5eW//fZb3Ra9Xi/c8Pf3LywsFG4XFBQoFAprOtTf3//ChQvWh+Tn5wcGBjLGampqxo4dGxYW9sEHH9yxB+tm2dnZ1kYMhAAAICXxE0LWuXPnO2aEVr179z506NCcOXMYY4cPH+7WrZtOp7Pe9cYbb9TW1mq12uvXr+fk5HTr1q22tnbixImenp6JiYnCqcSePXtevHixoKDA19e3vLz8119/7dWrl0ajWbVqlbAfJWMsMzPz5ZdfHjx48OTJk0+dOsUYM5lM77777rBhw6ZMmXL+/HnbXgsAAHA0CoX404T1Dp1PPfXUyZMnFyxY8OWXX86dO/fVV19ljM2aNeuJJ57o1atXx44dH3/88c2bN0+ePHnatGk+Pj4vvPDCkSNHunbt+vbbb7/11lsbNmwICQmZMGFCfHz85s2b4+PjBw8e3K5du0ceeaSsrGzmzJkbN25U5uTk9OrVS6FQvPrqqyNHjqyoqGCMrVq1av369XPmzOnUqdPAgQNLSkru00sIAABQh6+v79GjR2tqag4dOrRu3br4+HjG2NChQ8eNG8cY27FjR5cuXXbu3PnYY4+9//77jLFBgwbNnz+/7klBxtjnn38+fPjwnTt39u/ff+PGjYwxnU73008/eXh47NmzRzFz5sySkpKvvvrK+gCTyRQcHJyUlBQXF8cYi42NHT9+/AsvvMDt4uTpMybMfq3BXgEAAGgsxkcE3HWbXn2/PneuQNRuTcZzw4YUbNmyxdZ+/VXKn376acCAAYsWLXr66aeFcTI3NzcrK6tPnz7CFn369ElJSZGrfwAAYEdsODKqkLvYqPLmzZtLly4NDAwcOnTowoULV65cmZub6+TkJFy9yBjz8fHJzc2lHm/NqgIAQNNmMpnuuo20ZWXuD7WLi0t8fPxzzz3HGFOr1bNnzx43blxtba3JZFKpVIyxyspKd3d36vHW9A4AADRtwqDQ9ChDQkKs12QEBAQUFRUFBgYqlcrr168LjdeuXWvVqhX5eCVq0wAAwH/ZUHRb7kmhcvLkyTt27DAajYyxH374oVevXq6urqNGjfr0008ZYzk5Odu2bRNSOgBwdxbiH4BjEH9kVP6Do+pp06bt2bMnIiLC09NTr9f/8MMPjLEVK1aMGDFi3759mZmZU6dOjYqKkrufAAAADULt5OS0bdu2mzdvGgyGsLAw4VBnRETElStXLl265Ofn17x5c7k7CQAAdkL8FE/2KeHtSw6Dg4PvuEOlUkVGRt73/gAAgD1rBOf8xELUBQAAHBqKbgMAgGSEC+pFPkbmKSQGQrBvFskCmcSOpNo/tR/yfwAJ/mtoyH0D8NkQA5V7HMShUQAAcGyYEQIAgHTsMCyDgRAAACRjQxFt2YtuYyAEAADJ2OGEEOcIAQDAsWFGCA1PdB6Td4/YnVDbi4yZSpdKFUehIH4xcRCJ22oh/jKn/2Andm53f+GDjGyYEsr9AcNACAAA0rHhHGED9eSe4dAoAAA4NMwIAQBAMgrxx9JlP/aOgRAAACSjUIi/HELukRCHRgEAwKFhRgii0UFK/h3U9uR+eHeYzdS2xC81i+2MuO3pFCv/DgUVCKCaiXalkgh28tqVVMSU+AOYaqeJiLDWewc0IUiNAgCAI7PDcRCHRgEAwLFhRggAAJKxpdao3HNCDIQAACAdOzw2ioEQxIVW6tmeypuYidyK2SRie2pjk5idMDpEI3Z7qUqvcUMujA7FqFTE9rx2lYp/7kOpFrlzKqFD/f9FnXIhar7JHZ4HKdlwHaHsAyHOEQIAgEPDjBAAAKRjh7VGMRACAIB07PAcIQ6NAgCAQ8OMEAAAJKNQoOg2NAKiS6AR1cvItKfIFKjJyP8FRqLdZOTsx2QQuxN+u6ioaj3bi02NUtXLqEAmNwXKGFOp+TviBj7VGv7Gao2KaKe2F/FLGd15uoQb0qRNh4IpxF4XeNft09LS1q5dW1xcPG7cuNGjR9e9q6ys7P33309LS4uKinr++efVanVNTc2JEyfOnDljNptffvllxlh+fv7nn39e91EjRoxQq9Xbtm0TfsShUQAAaLwKCwv79Omj0+liYmKmT5++efPmuveOHj36zJkzI0eO3LhxozDs/ec//3nhhRd27ty5ePHiP++tuLh43rx5ZrP57NmziYmJQiNmhAAAIB0bwjL1Wr9+fY8ePZYtWyb8uHLlykcffVS4ffz48XPnzmVnZ+t0up49e3bo0GHJkiVDhw4dOnTosWPHhg0bJmzm5+c3d+5c4fZHH33Uo0ePLl26XLx4MSwsTGjHjBAAACQjnCMU9a/+gfOXX36JiYkRbsfExPz66696vd56lzBZZIyFhoa2aNHi9OnT9XcvMTHx6aefFm7fvHlz/vz5q1atwkAIAACSUdyuNipOPTvMycnx9fUVbvv5+VksluzsbOtdPj4+1i39/Pysd3GlpqZevHgxPj6eMebr6ztkyBA/P7+LFy/i0CgAAMjsyJEjffr0qduSlJTUpk0bxpiTk1Ntba3QKMwFnZ2dhR+dnJwMBoP1IXq93snJqZ7f8tlnn02YMMHb25sxFhcXFxcXJ7RjILRjoteYFZkC5aY3GWNGIsBpqDXxt68ltjeI2N5ApUaJnUuWJiVTo8SrTy60K7KmKFUOlEiNcoOdGi0/HarV8ds1RLtWx/+l1P7JlCnReSV1ZIqqcUpsLvt12cCYjecIIyMjV6xYUbclKChIuNGyZcuMjAzhdkZGhk6n8/Pzs9514MAB4bbZbM7MzGzZsiX1K2pqajZs2PDtt9/++S4MhAAAIB0F+Wcf+QgFa9as2R0zQqtx48bNmzfvtddec3Z2/uqrr8aOHatUKs+ePavX60ePHv23v/3t999/b9Omze7du52dnaOjo6nf8sMPP3h6elpPNxYVFTVr1owxZjQaMRACAEDjNXr06MTExKioqODg4NTUVGEKmJiYWFBQsGHDhgULFvTv3//BBx/8+eefP/74Y7VafevWrQEDBuj1+vLy8vDw8LCwsP379wsPmTZtmvK/xx+mTp168+bNoKAgnCMEAAApSX31BFOr1du3b09JSSkpKendu7erqytjbNGiRSaTiTG2YMGCRx999OrVq2vXrg0ICGCM+fv7JycnWx+u0WiEG//6179atGhhbf/uu+9SU1MLCwsDAwMxEAIAgHQaoOi2QqGIioqq22LNkTLG2rRpI8RqBGq1+oEHHvjzTkJDQ+v+qNVqrfvE5RMAAODQMCO0A6JXhBe5mLuRSG8a9PyAZa3eyG+v4e+HbNeLaBcbPaVqk1JRWKlWuheLrDVKpEapQCY3wKl14qc6dWQ7/z8EnTPR7sR/kanfKzZlqqRSpiJjo6hZej/d9bpAzkPkzvtiIAQAAMkoxP/lIftfKjg0CgAADg0zQgAAkI4drlCPgRAAACRjwzlC2Y+NYiBsRMhQjMjSaGJLoFGhFX0VPxRTU020E9tTYRl9jYjtDUQnqdJrZFiGKplGvchUvTqxC/aKLL1GlVhTE/kRbnU0qpSazplq5/+H4Oyq4bYbXPjbOxn47Vpn/qtDhXc01ItJvDhY+LcxsGWF+obpyb3DOUIAAHBoGAgBAMCh4dAoAABIx4brCOU+Nqo2Go3l5eXWn11dXbVaLWPMZDJdvny5efPmdSvZAAAANDHK48eP+/n5hf/Xli1bGGO///57u3bt4uPj27dvP3/+fLk7CQAA9kEIy4j6J3taRs0Ya9u27W+//Va3de7cuWPHjl21alV2dnanTp0mTpzYvXt3mXrYBElVMo0KTFIBSyqlWVPJb6+qNHDbqyv47TXE9nopSqyRqVEiCkuXUuO/aBZ+M70Ar0TpUKVK3IK91PaiFubVVfPbnVyoF1nkUszEJ9NZZGSXeu01FhGrEzNWTxaC82LKfpiuKbC31/D2OcKysjInJyfhoGhlZeX27dsvXrzIGAsICHj44Yc3btyIgRAAAJokJWMsLS2tbdu2np6eDz/8cF5eXmZmptlsDgsLE7YIDw+/efMm9XizmfhDGgAAHI/CJvL2WRkZGZmdnZ2Tk5OdnW00Gl944YWKigqtVqtS3T544uLiUlZWRj1er9ffr64CAICc7mXmY8M5QrnHQab08vLy8fFhjHl5ec2bN2///v0tWrSoqamprq4WtigqKvL396ce7+zsfJ96CgAAslIqm+al5394Vrm5uR4eHv7+/gEBAUePHhUajx49ihOEAABwL2yYEMqerlEnJiYqlcrw8PD09PSFCxe+8MILKpVq5syZs2bNWrNmzcmTJ8+cOfPtt9/K20v7xY29iU2HGo3iFsilaoRWE6nOynKivayW215FpkaJ/hBpVXKtXV4WkYqYUgFFcqFdIjVKHvIhi4fyiU17KlXE2rNi98PbXq3hv/IGPb+SBp0O5bcbRUZ2qU84ufoxVejVIq4SCJUm5dYmtaAw6V8k/7gmmjooKOjLL7/86quv/Pz83n777UmTJjHG5s2b5+TktHTpUl9f34MHD3p7e8vdTwAAsAc2FN2We+BUDxs2bNiwYXe0qlSq2bNnz549W5Y+AQAA3DeoNQoAAJJR0OUjyIfIPiOU+fcDAEBTYofnCJtmFhYAAOAeYUYoETJayLmDqqxoItOh/GwetSI8VQu0spyfAq0o4bdTqVEqfaqvFldTlKpUyS0rSgYXiRdNbECRDC4S7dSfvFQKlFpxnkqHKqjUqJg0qYqImBr04mqHim03UoVeiTeLivhSLz71dRM7D+GmSekr5agysiJ/a1Nnywr1cr+GGAgBAEAytpRMk3skxKFRAABwaBgIAQDAoeHQKAAASMaWc4QN05N7hxkhAAA4NMwIpUGFRrnrnlMZOSoYqa+mVpYXWTu0lEiNEu1UypTqD1mpUsxK9FR7LbHMPbXiPPWWiE2HUqVGqVSnkwv/j0udM/+7Ru2nlijQqideBy4qYqrWEu+UgUiHEgVdqUKvUkV2RcdDyXaqfCivXUPkchVUflj2+UzjomCiwzKyr0eIgRAAAKRjwwX1cv8tgYEQAACk0wgW2hUL5wgBAKBR+/jjj/39/T09PR977LHKysq6d/34448dOnRwcXHp06fP77//zhhLT09/7rnnoqOju3btat1swIAB4f/11FNPCY2rVq1q3ry5l5cXBkIAAJCMwib17PD06dMLFixITk7Ozs4uLi5+4403rHfV1NQ88sgjCxYsKC8vj4uLe/LJJxljVVVVoaGhjz/++PXr161bZmRkvPvuu8nJycnJycuXL2eM/fjjj6tXr/75558zMjJwaFQcKjpBrrXLSwGYiIiB2FJq1AK5VURptAqinQrFUKXUxBboIkusiWmndiK2OpfYvAb1FXV25X93fP1duO3Bbb247dRCu1lXy7jtOTfLue3cDwmVxFFTIRejit8u8kWmF9TlNpN3UC++ZO28iQC3kTGmoI6fUdvb2+HBRuuLL76YOHFip06dGGNz586Nj49fsWKFcNeOHTt8fHwmT54s3LVq1apLly516tSpU6dOx44du2M/QUFBDzzwgPXH9evXJyQktGnThuHQKAAANGaXL18WRkHGWKdOnfLy8kpKSoQfr1y50rFjR+G2q6trWFjYlStXqP1MnDgxJCRkwoQJwjZXrlyx7hYzQgAAkIwNtUYVClZVVVX3SCZjLDAwUKvVMsaKi4vd3NyERnd3d8ZYUVGRl5fXHXcxxjw8PIqKiri/Ys2aNZ07dzYajatWrRo8ePCFCxfqPhYDIQAASMaGyjJMwX766aeYmJi6bbt27erQoQNjzNfXt6zs9pmC0tJSxpifn5/wo4+Pz7Vr16wPKSkpsd51h5EjRwo3Pvzww//7v/87duyYj4+PdbcYCAEAQGZDhgzZsmUL966IiIjU1FThdmpqamBgoDAvZIy1a9cuMTFRuF1eXn79+vW2bdvW/4sUCoVSqbRYLBEREWfPnhUaMRACAICkJA0KTZs27aGHHnr66acjIiLefPPNZ555hjGWmJhYVlY2c+bMGTNmfPLJJ1OmTFm8eHHv3r3btGmj1+vPnz9/6dIlk8mUkpLi5OTk4uJy/vz56Ohok8m0evVqlUr14IMPuru7jxo16oknnggKCsJAKJa4+Bx3JVIqXUmWUqsi0qFUapRopxbspVKpVFUzsjQaVQWNWI6V+raoeMulKqgqX8QvNYpcM5Z6B6kFdV3cNdx2vyBXbnunXi247WqNuKBmfnYlt52qb8ffuch6clK1k6glb8kUKH97UasWM+LNVREr8yqIEmvEvpnDLuSrUIgumVb/9h07dlyzZs1TTz1VVlY2ZsyY+fPnM8ZKS0tLSkq0Wu3WrVtfeumlhQsX9uzZMykpiTFWWFj47LPPMsYiIiKeffbZ4ODglStXrlq16sqVK1qtNioqat++fR4eHr169Vq2bNmkSZMqKysxEAIAgHQaYIX6hISEhISEui2zZs0SbkRHRx89erTuXYGBgadOnbpjDz/++OOfdzt9+vTp06czXD4BAAAODjNCAACQjg1Ft+WGgRAAACSjYArxS1PJPHLi0CgAADg0zAj5xMbhqBiekReYrCVqdVKp0epKamFeantx6VDq91LtVGqUCjo6u/EDlj4t+GU5uRm/vFv8tGRuBr/2ZhVZyJQos0m8tWo1/49FfbW4wrDUm6gh1sglX3wqykukZLmoJ0vWCBVZoJVczVhkOpQIgZIpUDUvb8wYU2v57Rre9iriHVcS+WEqquqooVFbLqiXvS4rBkIAAJCUnR0ZxaFRAABwbJgRAgCAZGwput1AXblnGAgBAEAythXdlhcOjQIAgEPDjFAcMj5HBCaNvAqQtXpxKU0yiEjUIKW2rxGZDqW2p8paarX8spktWrpx2zv18ufvx4mzn/PHc7kbF+dXc9tNxhpuu8HAT1eSeWDinS0v0XPbr18sFrUfcoX6a/wV6ksK+M+L+2ZReUazif9OkSvOE3VMqZKi1CLvKuLJiq8RSqRAiU+gVke087andqLW8F8FqvOOW4PUhgvq5X7uGAgBAEA6tpwjlHkkxEAIAACSscMJIc4RAgCAY8OMEAAAJGPD5ROyl5bBQEgga03xNzcRy71ycyVUlSwqnEKWRhNZMq2WqAqmp6p2EdsbiTVyqWpkOl74hTHm6ePEbXdy4XwsnV3FfVap1YCpd4pcWpnKm5Txt89IL+W2F+VVcdup/zLE1sPjVk1TqfjvCFUajfoPgVwgl8qJEEERg5r/iaLCLyrik1mtJurV6fjxMZ0zEZZx4jxfjRP/l2qIxI1KTcSgHLn0mr0dG8WhUQAAcGiYEQIAgGRsKbrdMD25dxgIAQBAMgoFeSy9nofIC4dGAQDAoWFGCAAA0rHDCwkdfSAUuwCviVqAlwhSGnjLpVIpTWqtV3LhXGo/VAqUWKvWSKzpSj0pqp1aEDgno4Lb/ttJftU0Na/G1c3f+WnMilJ+qTPuK19PO5UypXLCVBUx6s2qLK0l9kP8H0B9BAnc/ViI7zdVuow6QsUtesfoamQU+hNFvVlEKrWGitry+0NFcJ2cOe2GWv5OyLA0sRqwRcV/B+n/8+UeDSRih+MgDo0CAIBjc/QZIQAASAkX1AMAgKOT/VinSP8bCD/99FNXV9fJkycLP27btm3v3r2+vr7PP/+8vz9/uRwAAIC6bLmOUO6B8/Y5wk2bNs2dO3fdunXCj4mJiTNnzuzatWtOTk7fvn2rq/lLvgEAANg7NWOsoKBg6dKls2fP3rdvH2PMYrG89dZbH3zwwbhx4/7f//t/0dHRmzdvfvLJJ+XuagMhFtolMntmk7ggpYEX1KytIWqBUu1E2pNqFxuYNBqIMpvEk6XWttUT/b95pYTbXpBNld/kNFIL0pYU8tvJCC7RSe47xRgzE6lRcjlWtbg1Zqn9kGU8ibKc3OyiUsl/Z6laoF6+ztz2oDAPbrtGx+9MST7/Taks50dnjURkV00UhhX1dWN0dV/uN8ioF/cdpwq3krFfcjljot3e2FB0W/Q5RakpGWMvvvjiggUL/Pz8hKa8vLzff/990KBBwo+DBg366aefZOsgAADYEYX4f3JTbt++vbS09PHHH7c25eTkaLVaD4/bfwP6+fllZ2dTj6+p4f/pBwAATYyZOkJi59SzZ89OTk6u26TVak0mk9lsViqVjLHa2lonJ/5yOYwxtRq5UwAAh3AvxzBtCMvIPilUZ2ZmxsbGMsbKysoqKipat26dkpJiNpuzs7ODgoIYY7du3QoMDCQfj4EQAMAx3NNAyBQKkSOb3OMgU547dy45OTk5Ofnll1/u1KnTvn37PD09Y2JivvnmG8ZYRUXFjh07xo4dK3c/AQAAGoT6gQceEG75+vo6OTkJPy5btmz06NHHjh27ePFiVFTUwIEDZe2kBMTWFCXXKydjbEQsjRfUFJVhY3QWjqwRWkvE3qhF26WIwjK6xmkFUWazKI9/WQ43hkf+UqK2Jxm1JTpJPVmKmUh1mvi7pxd5p9YrJ/av0fKDmgpecQyqDKbOmX8gp2U4Px3aIyaI264lFm2/mJLPbb/6WxG3nfzEaqg0Kf/FEVXyl2o3EDuhqtGaTVRqlGgn5j+KJpMmbYBio0aj8fDhw+Xl5f379/fx8bnj3pSUlKtXr3br1q1169ZCi8ViSUtLq6ys7Nq1q9BSVlZ2/PjxysrKzp07C8NcVlaWNf7yv+/D+PHjY2JihNu9evW6ePHizz//3KJFiwcffFD2bCsAANgJ8ZdP1DsSGgyGwYMHV1ZWhoWFTZ8+ff/+/V26dLHeO2/evI0bNw4YMGDmzJmrVq168sknjxw5MmrUKIvFolQqS0pKGGNnz57t379/jx49fHx8pk6dOn/+/Llz565du/azzz5r2bIlqzsQ+vn5Wa+gYIz5+vqOGTNG1JMBAACQ1g8//JCfn3/69GmtVvv6668vWbLk+++/F+66devWBx98cOnSpeDg4AMHDkyZMuXxxx+PjIw8f/58ZmbmsGHDhM2CgoIuX74slEj76aefYmJiZs6cyRibPHny22+/zbD6BAAASEhIjYr6V/+h0a1bt44fP16r1TLGJk2atGPHDtN/z0Ds2LEjOjo6ODiYMRYbG2symU6cONGsWbNWrVrV3YOvr6+1UGhYWJjJZNLr9YyxwsLCQ4cOXb16FZlPAACQjtTnCG/dutW/f3/hdsuWLQ0GQ15eXkBAAGMsMzPTOuYpFIqgoKBbt27V/6uWL18+atQo4UTjr7/+evPmzdTUVIcfCMmwDL9dbIk1bgrAQEQDDCJrO4ltN0mUAhAbSTARJdyoNAF3/6KTRGQyiFgulSp1JvIrTceyxL0I5P6JGBe3hJvWzA+zUHXaXN213HZffxduOxW6yUjjr6JM1ZOjPoFkHIwIAVFvOv2J5bSTXx/iu0+WWKPSV9QnhIpNEbtptGwpscYUV65c+ec//1m38ZlnnmnRogVjzGg0qlS3P8kajYYxVlt7O39nMBisdzHG1Gq1wcBfHlzwySef7Nq1S6iVtnDhwjfeeIMxVlpa6vADIQAAyM1kMlVV/aH4sPVvxICAgPz82znk3NxchUIhTAeFu1JTU60Psc4UudavX79s2bJDhw4JV8brdDqh3dPTEwMhAABIxrYjo+3bt79jRmjVr1+/bdu2zZ8/nzF24MCBXr16abVaYebXr1+/RURlLYoAACAASURBVIsWVVdXOzs7X7lypaCgoEePHtydbN68eeHChfv37w8PD//zvRgIAQBAOlKfI3zyySdXrlz5wgsvtGvXbsmSJYmJiYyx2bNnFxQUbNiwoU+fPmPHjn3kkUc+/vjjGTNmeHl5FRUVrVy5Misrq6amZt68eb6+vrGxsY8//vioUaOSkpKEfb7wwgsvvfRS165dfXx8Dh8+jIEQAAAaLy8vrxMnTiQmJqanp3/33XcDBgxgjE2YMEFY8uH777//7LPPzp07N2vWLGFheZVK5e3t7e3tvWTJEsaYp6ent7f3HdNNtVo9bdq0Y8eOXbp06aGHHsJACAAAklHYsL7g3TYPCgpatGhR3RZhOGSMOTs7v/TSS3Xv8vT0nDt37h17+HPLqFGjRo0aJdx29IGQSuxRSTAqc0iXXhMRSyODjiIDkFS0leok1S52xVEqJkcFI6nfyw2I1lTxS6mJLY1GElkQi06B8ndP5CWZkvgvg14amvqkcdqpF6emkp+sy8mo4LZTJdN0TvxUai6xH7FvopJY5dgsMucsqp36+ojdOV1ijc/u0qEUG1afkP2544J6AABwaI4+IwQAACk1QNHthoaBEAAApGPDBfVyr+uAgRAAACRjjyvU4xwhAAA4NMeZEYpbmZdKfNFpUhGBTzIFStYCpYoZioulUVS8MpWMMQWx7Cr1IlDLq9K1TEUVaOWXnaSiqtzam6yed5bI/in5rwFZNpPMFpL95LdTB4vInDPvRaZeNGqp5OuXirnt5SV6bjuVGq2u5KdDqVWUqRdTRdQUJfPSxCeTelO4+5EsHUq9VeQdck+LHJjjDIQAANDgFFIvzHsf4NAoAAA4NMwIAQBAMrZcUC/3UWEMhAAAIBlb1iOUeyTEoVEAAHBoDjMjFBcapctmku0iYmlia4qKrflJ0RIZP2c3J247tV55VTk/c5h9s5zbXpTLf161NfxMI//5Ek9WQfwtR/2JSb2DOif+d8EvyJXb7uXDf9HKivkBy7xMfvlNPfEiULVJqSfG/aRxV2BnjNUQ6c2Sgmpue3UFvzaphsgVq4kaoVQKVEvlk6WoHcoYM5GLyPNSo1T0lEqHkl9PsfHtpsKGyjJyc5iBEAAA7gNcUA8AAGBfMCMEAADJKJhC7HWBck8IMRACAICEsPpEoyW22hFdNok6kc7fD/fEO3mqX6KT61QW2cWd/3Y/0KEZtz2svTe3vbSInwc5fSSL215Vzo9aUMvDcvMvdMk0bjMZeaDOXvgFunDb+44I4bZHdPXltqdfKOK2H9lxnduemV7KbRdbhotb5M/ITzUxIxFaMej54RqLhR+uoSoLaonkkRPxJooltgKiqDwLXRqN6gyxvaNSiL8uUO5xEOcIAQDAsTnKjBAAAO4DhUL8BfJyX1CPgRAAAKRjh+cIcWgUAAAcGmaEAAAgGVuKbjdMT+4dBkKJIl9iyi+Jrr0k0cdEq+O/3b4B/Cpioe34qdGSghpu+/WL/GVdNVpqcVtiQWDe14gsmUYWySNSoyr+URA3Lx23PbwjP1LbpW8Af/9ER8/+nM1tzyb6QwUyKdzXgawWJrIqGJm6pPpI7YfYnPyIS/TJF7127l/e2JGh6DYAAICdwUAIAAAODYdGAQBASii6DQAAYE8wIwQAAMnYEpZpoK7cMwyEEr0FxBuv4AUjRUekJIqrGfT8NWALc6u47TeulHDby4haoyWF/GVdDbX830sVVxVXAZLAfeXr2VF5Cf9JpZ/j1w5VEfunao1SL5pZZDqUwv1MKYlOUu3UJ5N6MelVkYnt+Zszuo4n+QBRyP6I+SLKHWy0H3Z4QT0GQgAAkAwunwAAALAzmBECAIBkbKgsIzsMhAAA0Kilpqa+//77xcXF48aNS0hIqHtXfn7+ihUr0tPTo6KiXnnlFScnp/Ly8kOHDp05c8ZsNi9evFjYTK/Xv/POO8ePHw8NDZ0/f36LFi0YY9nZ2StWrLhx4wYOjQIAQOOVk5MTExMTERExderU11577Ysvvqh778iRIwsLC59//vnDhw/PmDGDMXbixInVq1enpqa+99571s1efvnlPXv2PP/881VVVUOHDrVYLGazefDgwXq9/rnnnnOUGSE1U6em8HTMjIjhEfvhxvOoxdapnYhFVVasLOcvWJ5+nh90zLlZzm2nVpzPvsHfXl/NX9+cLhPKaTSbiK3FvPL17Kcgq5Lb/vPuG9z288dzuO1lRPo0P5O/fyI5K/rDoOLVLFUTVV7VxAr1Gh2/XefM/49Cq+PvX6Xm74eM8opEfg3JdKsEVW3Jeqj2dhiwoUkelvn888/79es3Z84cxlhVVdWyZcumTp0q3HXkyJFr16798ssvarW6Y8eO4eHhK1asGDRo0KBBg44dO7Z//35hs6Kioi+++OL8+fPh4eGDBw9u1arVwYMHDQZDWVnZ2rVrFQoq/gwAAGAT4TShqH/1OHny5EMPPSTc7tev37lz56qrq6139e7dW61WM8aCgoKCgoLOnDnz5z2cP3++WbNm4eHhjDGlUtm3b98TJ04IuxXGYAyEAADQeOXk5Pj4+Ai3hRs5ObePx+Tm5jZr9r/1YXx9fa133bGHupv5+Pjk5OTUbbx9xMNgMKjV6jvmp7W1tVqtVrJnAwAATZ5NF9QfPHgwMjKybtv3338fERHBGHNxcampub30mzAXdHNzE350cXHR6/93MqK6utrVlbOonKura93NampqmjVrVltbW1Jyu2aI8ocffmjVqpWnp6erq+vo0aNzc3MZY/n5+XFxcb6+vr6+vp9++qnI5wQAAA5KIR5jLDo6evMfhYSECDsMDg6+fv26cPv69evOzs7WCWKrVq2sd5lMpoyMjFatWv25Sy1btszKyqqtrbXupGXLlq1atbp27ZrQouzatevRo0erqqry8vIYY3PnzmWMzZkzJyAgoLi4+NChQ6+88srvv//eIC8YAAAAY+7u7pF/5OTkJNz1yCOPfPfdd+Xl5Yyx9evXP/LII0ql8ujRo4cPH3744YdTU1PPnz/PGPvhhx+8vb2joqL+vPPOnTu3bNly06ZNjLErV66cOHFi7Nix48ePP3r0aFpaGmNMHRYWJmzq5uYWFxe3Y8eOmpqaTZs2HT9+XKVSderUaeTIkf/+97/feOON+/NyNBSRiS+yiCLZzt+RihcQVROZOiprR1eG5HeGUlvDr/lZXckPNOZn8fejr+KnRiuJNKnRwC+nSfWf/3yp0qFEqU6LyOyfnnhxsq6XcdtzbopLpVIRWSpCLBb3k6bR8j9RTi78FKi3nzO33dffhduuJdKk1RX8T0JNFT8/TNZEpfLVYtvFpEmp77JSbCFWR42T2nBBff3bDx8+vHfv3h07dgwMDMzLyxOyoBs3biwoKNiwYcOKFStiYmI6dux44cKFpKQkpVJ548aNbt26GY3GiooKISNz8uTJNWvWTJ48OTEx8cKFC0uXLg0ICGCMLVq0qHfv3pGRkWrGWHFxcXJy8o0bN7744os1a9ZkZWXV1NS0b99e6ET79u0vXbpEdZFK6gMAgCOyreg2PZIolcqvv/46LS2tpKSkS5cuGo2GMbZ8+XJh9HnxxRcnTpx448aNdu3aeXp6MsZatWqVnp5ufbhKpWKMDR48OD09/bfffgsODhZGQcbY3LlzExISMjIy1Iyx8vLyAwcO3Lx5U6PRuLu7l5SU6HQ6IZDKGPPw8Cgq4l9nxhirquIvXAAAAE2M2WxWKu9yrYGCMYXIkfBetm/dunXdH+uGYvz9/f39/a0/KpVKb2/vP+/B3d39wQcfvKMxMDAwMDBQzRgLDg5et24dY2zdunVTp07dvXt3TU2NXq/X6XSMseLiYqEaDRc3ogMAAE3PXUdBO/WHZxUZGZmdnR0QEODl5WW9LPHMmTPWw6QAAAD1UNhwQb3cfVZv2LDhgQceCA4OzsjI+Mc//vHwww9rNJqnnnpKKOl28uTJQ4cOrV27Vu5+/nVUKoZoFlnDiVvgihH5FyoUoyIKX6nUIs/Sc1tpJiP/8Dy1oC6VK6FCMRQl8aKpeREPjZFfzUvsL6VOy3PDJvU8wEytKkz9XpF/SVNn36nIhpL3IdEQJdbcPPjXB4e24xxQYoxF9uQfE3IiwjLXLhVz2zN+5y/1TMW4qE8+9WaJDd1w90PunGgXvepvkw/RSH2O8D5QV1RU/P3vf8/Ozvb19R0xYoRw+cSbb775yiuv9O3b18/P79tvvw0MDJSzjwAAAA1GPX369OnTp9/R6uLi8vHHH8vSIQAAsGPii27LPkt2lNUnAADgPrDlOsKG6cm9a5oRIAAAgHuEGSEAAEjGpvUIG6gv98rRB0Lq9afToUQ7EW/jroBKLYsqtvQaHW8Tlz5VmfjtJqO4OBxZf44MOvK3p5Z75TLoiQirkZ8mNfGLfDGm4KfWxC7dLLbUktjaTFR2kfuhoj45OqLEWvMgN257RDdfbjtVqq2qgr8EdH5mBbedqktHfSOUVLuYdCi1PfX1EVvXjUyTcltBVjg0CgAADs3RZ4QAACAhG8Iysk+TMRACAIBkFEz8OUK5R0IMhAAAIB3bKsvICucIAQDAoTn8jFBkAJIsj0kFQXnt1HKpGp0EO6+nnaxlSmT2lCp+8JLav8bET3uqVOKCkWYdZ3vql+qJzF4tkSatZUTKtJaqHUqkSUVGZEW3E7ll8sOj4bz4VKCRSmlWlOm57fmZ/KWbdUSt0fISfmrURPxe6pOpJmqlUmlSsp36xPJeTPLrQ6VJyVV/uc1iyx7bH8kX5r0PHH4gBAAACdlQYk3uPwNwaBQAABwaZoQAACAdOwzLYCAEAADJ2OE4iEOjAADg2BxlRig6xUTVGiVrivLjbdyymVonERszxjREO5mp0/KDkWoD/+8eaoV6Kj5nNhPpUCKzR62T7uqu4bZz363SohrutiX51fx9VBj4OydQHxAzPzlLV6MlPiFkpUpqUXWqnSo8y3uzqCelr+ZXXM1ML+O2W4gXgfrElhTy3xSqMKyaiMLS1XdFFFxl9CeQ266hUqMiC5yKrjUq+7RIIgr6uZMPkfu5O8pACAAA94MdHhvFQAgAAJKxpdao3HCOEAAAHBpmhAAAIBmbFuZF0W2ZEZEE4n0RW2KNm3PROvFfdrJdx482UFGFWiIaYKwl1qo1EO3Ek9UQGRe1C397/2D+cq8PdPThtnOfV/q5Qu7Gl07nc9upBBCFegepnAhVQIuKVFC5DyqWxYjSblQzdz/UzqkXh0okGYhPCFXvjfp/jdpeslAM8Y0QFU9TU9k04peSJdaoN9beDhuKZofnCHFoFAAAHBpmhAAAICGF2PUFsR4hAAA0HYpGcF2gWDg0CgAADg0zQgAAkI4dhmUcfSAUm+xSEQkxUbWddESGjVrmlGqvreGnSQ16/v6pdKjRSKRGRQYvnVz5/WzRyp3bHhndnNvu4saJpRpq+dW5Mq/zq4KVl/DXmKXqw1HBRfKTQGzv5MJ/EVzdtfzdE7+2por/5lLV0cxmzptFrSVLpTEtxBtOlUajUJ9weilpESXQGGM6Il9NvfhUOzenTUZbqdSoyBJr8v+v38Ds8fIJHBoFAIDGzmQyVVfza9gyxioqKkRtfwcMhAAAIBmhxJq4f3fb5/Lly5s1axYYGDh8+PCSkpK6d+3atSswMDA0NLRDhw6pqalC47Jly4TtR4wYUVpayhjr1q1bszqmT59eU1Nj/REDIQAASEdh0z/a8ePH33nnnbNnzxYUFDg5OS1ZssR6V1VV1ZQpUz777LOCgoJp06Y99dRTjLFjx469++67qamp+fn5Wq1W2P7w4cPp6enp6em///67s7PzkCFDLBZLcXHxhQsX0tPTMRACAIBkFOLVPxAmJSU9+uijoaGhKpVq9uzZ//73v613bd26NSgoaOTIkYyxmTNnXrp06fz580lJSZMmTQoJCVGr1dbtPT09vb29vb29U1JS9Hr96NGjhT14eXl5e3tjIAQAgMYrPT29ffv2wu327dsXFRUVFxdb7+rQoYNw29nZOTQ0ND09PS0tre72hYWFdY+mfv7550888YROpxN+DAkJad68uaOnRknUoppKfq6OiuFxE2hU5UMnZ367nsi8GfT8op8GoqYomQ418dstRIiQXLWYiM/VEpnD0kJ+ZcvaGs72NZX8tCRVe5MOTBJP1ixuQV1nV/6LHxDKj8iGtPXitlNZxMxr/DRszo1ybjs3ZUpl8aglncUGJqmorfgFcvnbk98U4htBvSnU9tx0q5boJLnksshao008M2rr1ROlpaVnzpyp2xgREeHs7MwYKykpcXV1FRrd3NwYY8XFxd7e3oyx0tJSFxcX60Pc3d2Li4tLS0v/vL2XlxdjrKioaOvWrceOHWOMaTSa//znP9HR0UVFRRgIAQBAUuJH+19//XXq1Kl1WzZt2hQREcEY8/Pzs07phLlg8+a3r7zy8/NLT0+3PqS4uLh58+b1bP/VV1917ty5S5cujDG1Wt2/f3/GWFBQEAZCAACQ2cCBA7ds2cK9KzIy8vTp08Lt06dPBwcHC/M8xliHDh0+/vhj4XZJScn169fbt2/foUOHutuHhIRYJ4jr169/7rnn/vwrMBACAIB0pL6g/umnn46Kitq9e3dERMTixYufffZZxtiHH35YWlo6d+5cxthbb701ZcqUpUuXxsbGhoWFPfPMMz179tyzZ0/btm2XLFkibM8YS0lJuXz58qRJk4Qff/7554KCgo4dOxYXF2MgBAAAyQiXBop9SD3atm375Zdfvv7662VlZWPHjp0zZw5jTKlUqlQqtVq9c+fO2bNnf/bZZ9HR0V988QVjLCIiIikpadGiRcL2r776qrCf06dPz5kzRzhZKPjwww9v3Ljh5uaGgRAAABq1MWPGjBkzpm7LjBkzhBuRkZF79uy5Y/uxY8eOHTv2jsZnnnmm7o99+/ZNTk4WbmMgFIcKTFLZQjU3Naqjaory05XOLvx2qgIkVZbTRKRGzSZxNUWVSv7+qf3k3bqz+pHgzE/8/nDLeOZl8neiJwquqqjgoomfCTQTEVmqpqi7F792aGgEPx3aM64Vt50KRp79OYfbTtUgLcnnVJOi3lc61Um0U7VAxa4IL7LKrjNRvdaZV422nu2p/XP7w/3OMvqTQKVGERu1IxgIAQBAMgqF+CLaKLoNAAAgI8wIAQBAMjasUC/3hBADIQAASOfe1pO44yHk+ez7AwMhn9i/UKgT5twaVFod/4g0dUqfKplGhWKMxAK8VFjGJDIsQyWGuKXRGGPVlQZue+bVUm47N3RTTYREzMTqwVShLOqdsvAWtmV0DEorchVlai1ZMj9C5U2IfAo/t0K8sVQehCqBRv1SseEXsjQaEX5xceMnkrhLN7N6SqyJCctQySCUUrtXdrhCPc4RAgCAQ8OMEAAApCP+gnrZYSAEAADJKGwqsSbzOUKLxXLy5MmTJ08qFIqYmBjr2k6pqanJycl+fn4TJ04U1sIAAABoepQfffRRQkLChQsXzp0717t3b6FW2/bt2wcOHJifn//111/HxMQYjcRScAAAAHUpbPonK/Vjjz02Y8YMpVLJGIuOjl6+fPnUqVOXLFmyevXqp556ymg0durUafv27ePGjZO5p40DFQtWEKkjblkmKpunc+YfHqBSoEYDPyNHpUDFtpML8xLHPahYHbUwL1UijruAMJXqpEqpUWlPjZm/vZnYP/Vkqf7kZ1Vy288f55dMowp3ZV3nL8xLfRi41dGozlMl1ho6HUqlQF3d+Z9kVw9ie6KdSp9S/eRGbVUacR9vxEbvYEvR7Ybpyb1T+vj4CKMgY8zFxUWhUBQWFqakpIwePZoxplarhw8fvm/fPlk7CQAA0FD+94dhVVXV0qVLX3rppaysLI1G4+PjI7QHBARkZWVRj6+trW3wPgIAQCNAHSWqS3F7Tijin+xzwtsDYW1t7cSJE7t27frcc88pFIq6z9ZisYiuoAoAAA7L3s4RKhljBoPhscce0+l0SUlJSqUyICDAaDQWFhYKW+Tk5AQEBFCP12r5B+sBAKCJaaqTIqXJZJo6dWpVVdXGjRs1Gg1jzMfHp0ePHtu3b2eMGY3G3bt3DxkyRO5+AgCAHRB9YLQRXICvXrNmzcaNG0ePHp2QkMAYUygUmzZtev3115966qlLly6dOXPG09NTCM4AY+QUnkqTKnnpPCroSGX2TEQGj1oI18LPFZJBR6qdolTyL6ehgppU3I4s3shv529sMlHPlt9MpUPJF4fYD9VelMtZIJcxVlPJf9GovDEVtaVO0nCDkWQJXLKmKPHJJAqlUjU8qVqgLlQKlEyNEvshtqdqjZKpUV5ZUaWK/yJQ75Ts/4k3NrZcUN9AXbln6iFDhmzatOmO1tGjRx8+fDg5OfmJJ54YP368Wo0CNAAAcA/ssOi2ukOHDtZqMnV17NixY8eO979DAAAA9xOmegAAIB0bzvnJPiOU+fcDAEATohAfLhW7kK/ksB4hAAA4NMwIpUH+AcRNjVJlMIl1w81m/ttEBx2JdiJJSWVGqSdF9V9fzQ86kqlRYj/c/auJmpzcwqSsvkgtv51Kk5I1SLmtNH0NEbUVF50lV7rnvpj0J61ha4q6uhMryxNpTxdieypNSqVSnZz5/dcQaViVmvP6kCVF5Z61QMPBQAgAAJKxbT3CBurMPcKhUQAAcGiYEQIAgGRsqRQj91FnDIQAACCdRlBEWywMhA2N84kQtYovY0yr429/L0ui/OH3kguLiiuBRpWIU2v4eRA1kRMx1PKjDUYdJ3RDlRyjFqqlVhs2EyXZzESlNrKWGoEKuVC5Faqgl9iEEXd7NbHGrFbH/+KLDcU4uxIl1qiwDBFyoRbUpfYvaqFdRiet+C8yFtr9a2w6R9hAfblXOEcIAAAODTNCAACQjvhzhJgRAgAAyAkDIQAAODQcGgUAAMkomA3rEd5l+5qamh07dpSVlQ0ePLhVq1Z33HvgwIH09PSoqKju3btT21++fDkjI8P6kNjYWKVSyRg7ePBgWloaBsKGRXweqEAgfydqCzFxJ7J/1KeKCjQqucsH0ylWqkCXzsnAbddX8/tJBUG57QZiY0Mtv91k5Kc9qYV8ydWMyZV5+c0U8sUXu5oxsb2aF+Wl3yl+O5UOdSIWvKVSndQCuVQ79XvJdChRjJD7IjD6xed+P2U/X2XvJL+OUK/X9+vXz93dvXXr1q+88sqePXt69uxpvXfmzJkHDx6Mi4tbsmTJggULZsyYUVNT89BDD3l6eoaHh7/yyit79+6Njo5eu3bttm3bwsPDhUfFxMQolcqXXnpp3759Q4YMwUAIAACN1+bNm41GY3JyskqlCgsLW7p06Y4dO4S7rl27tn79+qtXr/r7+8fHx48dO/bpp5/evHmz2Wzet2+fSqUKDQ1dunTp9u3bGWMJCQlLliyx7vbGjRufffZZenp6YGAgzhECAIB0FDb9o+3atWvMmDEqlYoxNn78+L179xqNty9N3r17d69evfz9/Rljffr0UalUx48f37lzZ93t9+zZYzKZGGMZGRlbtmw5c+aM8FhhZhkYGMgQlgEAAAkpbFLPDjMzM4XhijEWFBRkNBrz8vKEH7OysoKCgqy/NzAwMCsrKysr68/bazSaq1evJiUljRgxYtiwYXq9vu5jcWgUAABkduHChTlz5tRt+fvf/x4QEMAYM5vN1pFSCDQIM7w77hLuNZlM3O3feust4XZFRUV0dPS6devqboaBEAAAJGNDWEbBmE6n8/X1rdsoHNtkjAUEBFingDk5OUqlUjgWyhjz9/dPSUmxPiQ3NzcwMPDP27do0cIaCXRzcxs+fPjZs2e7d+9+9OhRoREDoQzoT4nINClVtJSKyBGbkwvk8pYtZXRmj1rWtbaGH+yk1qrlbk+mRqlaoyJrkJIL9lKrHFOrIlM1S0UWeqVefKpsJrfMppaovakjFrClUppi055OxH60RDuVYqVSr9SLQ0VtyfQ2AqINRHxqtHXr1nfMCK1iY2M3bNiwYMECxtiePXv69eun0WjKysrMZnNsbOyCBQvKyso8PDzOnz9fWloaFRU1cODATZs2/eMf/xC279+/v0bzvw+wxWI5depUXFxcbGzsnDlzSktLPT09MRACAICU7npdoChTpkxZvXr1lClT2rVr984773zzzTeMsYULFxYUFGzYsGHo0KEjRowYM2ZMYmLirFmz3N3dExIS3n777YSEhIiIiHfeeWfTpk2MsX79+vXq1cvT0/PAgQMZGRkzZszw9fUdNWrU8OHDx40bh4EQAAAaL3d39xMnTnz11VclJSUHDhzo1q0bY2zKlCl6vZ4x9s0332zcuPHKlSurVq0aPXo0Y8zDw+PkyZN3bL948eKUlJSKiopp06ZNmDDBxcWFMfb1119/8803ly9fxkAIAACSseUc4d229/Hx+dvf/la3xXpNvVqtTkhIuOv2gwYNGjRo0B2bqdXqKVOmMJwjBAAACdmyHqHciz3iOkIAAHBomBE2IuRfUVTZSWJ7NfH3jdhylypi8XQ1lRolMopUarRWL6IGqbGWH8c0GIhao1RqlKpBSqVGqZSpRDVIqaAjFYxUEeU0uQFLKsdLpTR1TuJqflL7oX4vlQKlaoQqqZXlqW8K0qGNwd0qxfAfIisMhAAAIBnJi27fBzg0CgAADg0zQgAAkExDrEfY0DAjBAAAh4aBEAAAHBoOjdoBqWqTUvtREg8gU6NkcJEf4NQ5EYFPYnF5Iy/wSdUUpdKkRiOVGuW3U+lQ0alRkTVFxaZGqVqj3DQpldKk8r3cgqWMMa1ORFSV0Z8QFZUCJV4EKucs8gsB91VDXFDf0DAQAgCAdMRfUC/7XzAYCAEAQDo2XEcoN5wjBAAAh4YZIQAASAbnCOG+EhuiUZAhGqowFT8PQpVkoyISJh2VZ+F3iFsFjSqZJlkoRnSJNW5zPWEZfjuVE6FXzEMvPgAAEN5JREFUSyZyKLztqXdErREXciHDL0Qnqc6T4ReURmtCbCm6Lfc7jUOjAADg0DAjBAAASdnbVB4DIQAASMYOQ6M4NAoAAI4NM0IAAJCMPYZlMBA2QWLTpNRxASL6RwYdLSp+kpKsFmbi/2JugJNKb5qoVCeVAhXZTu2HTI2KRK6WTJVeExPUpIrkKalVf6mdS5X2RDrUEdjhwrw4NAoAAA4NM0IAAJCOQvShThwaBQCApsOWyjIN05N7pywvL1++fPnEiRMHDx5cUlIitFZVVT377LMhISHR0dG7d++Wt4sAAGBPFOL/yUpZVlZ28+bNPn367N+/v7a2VmhduHDhtWvXjh49+tprr02aNCkrK0veXgIAADQQdVBQ0Nq1awsLC2fNmiU0GQyG9evX79mzJzAw8OGHHx44cGBSUtL8+fPl7Sg0HLEZPyZuSVo6GEnUPuUmMs380qHkArli057kfqh0KLUf4g4FWYSUaBYZ2eW+yEoqpUmE5OhaoPztyaq24jaHJkXBFOSnvbHinCPMysoqKSnp1q2b8GO3bt0uXrx4f3sFAAB2yR5Xn+D8ZVhYWOjk5KTVaoUfvby88vLyqMdXVlY2VNcAAKAxMVNHZuwcZyD09vauqakxGAzCj6WlpT4+PtTjXVxcGqprAADQmCiV93Dpub0lZRh3IAwICHB2dr506ZLw46VLl8LDw6nHy379BwAANB4Km8jbZyVj7Nq1azdu3GCM3bhx4+rVq05OTpMmTVq5cqXJZDp37tyuXbsSEhLk7SUAAEADUTPGhg0bZjQaH3jggfj4eIVCkZaWtnLlysmTJzdr1kyr1a5evbpNmzZy9xMaEakyhBYxfwVSh2Qs1O5FpkMpUtUUFUt0uU4RreLfQRz3gXtmQ1hG9qOjasbY5cuX72j18/Pbt2+fwWDQaDRy9AoAAOxT4zjtJ0p9JdYwCgIAgEjSL8N09OjR1atXl5aWjh07dubMmXW3v3Xr1uuvv3716tWoqKhFixa5u7szxn755Ze333677vbff//91q1bMzIygoODX3rppe7du1+/fn358uXCTrD6BAAANF63bt0aNmzYkCFDXnvttTVr1nzyySfWuywWy4gRI1xcXJYvX37lypXp06czxjIyMoYPHz506NDXXnvtgw8+WLduHWMsMTGxb9++ixcvDg8P79+/f1paWl5e3rZt2+Li4uLi4lB0GwAAJCP5coT/+te/Bg8e/OyzzzLGli1btnDhwueff1646+DBg/n5+e+//75Sqfz000+Dg4OzsrL+9a9/DR06VBgUly1btmjRoueee27nzp3CQ/r3779r1679+/d3797dzc1t4sSJDKtPwP0grrgYcQ8ZuRH7S0XuqKHDMg1ZjgylzkAGUo+Ep0+f7t+/v3C7d+/ely5dqqqqEi5hP336dM+ePYWrG1u0aBEcHHzu3LnTp0/HxMRYt7948WJ1dbWzs7PQUlNTc/369dDQUMZYYWHhlClTXF1dcWgUAAAar9zcXG9vb+F2s2bNhBbhx7y8POtdwr05OTn1bM8YmzVrVmRk5NChQ729vWfNmjV+/PiIiAjMCAEAQDK2XSC/d+/eoKCgui3JyckdOnRgjLm7u1dXVwuNVVVVjDEPDw/hRzc3t5s3b1ofUllZ6eHh4ebmRm2/ZMmSn3766dChQwqFok2bNgsXLhTaMRACAIBkbCu6PWDAgMTExLqNzZs3F26EhIRcvXpVuH316lU3NzdhnscYCw4O3r59u3DbYDAIodA7tnd3dxcmiCtWrNi0adOhQ4f+XDQUh0YBAEBmTk5OQX9kvX4vPj5+8+bNRUVFjLFPPvlEKPxy8ODBnTt3jhkz5uLFiydPnmSMbdy4MSAgoHv37sL2xcXFdbd/7733Pv/88/3797do0ULYbVpamrAEr16vx4wQAAAko1CIrkFd//aDBg0aPXp0+/btfX19VSrV7t27GWNbtmwpKCgYOXLkhx9+OGzYsJCQkKysrE2bNikUisGDB48cObJdu3Y+Pj5qtXrPnj2Msblz56rV6o4dOwr7fOWVV/R6/XvvvdeqVausrCyF2IpTd5g8fcaE2a/9lT0ANCkiU6wAdmR8RMBdt/k+JaOwslbUbs8cSc48uWfLli31bJOXl1dSUtKmTRth1DSbzRaLRaVSMcbKy8szMzPDwsJ0Oh21PVdxcXFOTk7z5s0xIwQAgMauefPm1rOG7I8LQrm7u7dr167+7bm8vb2F04cYCAEAQEp2d/gDAyEAAEjGhssnZF+PEAMhAABIR/Iaaw0Pl08AAIBDw4wQQFJy/20LIC/bLqiXFwZCAACQjG0l1uSFQ6MAAODQMCMEAADp2GFYBgMhAABIRiH+nJ/sR1JxaBQAABwaZoQAACAZBVMoRB7rFLu95DAQAgCAdGw4Ryg3HBoFAACHhhkhAABIR/wF9bLPIDEQAgCAZFB0GwAAHBvOEQIAANgXzAgBAEAyCmbDodEG6su9wkAIAACSscMKazg0CgAAjg0zQgAAkI4dTgkxEAIAgGQUChsuh5B5JMShUQAAcGiYEQIAgHTEV5ZBahQAAJoQOzxHiEOjAADg0DAjBAAAydhyQb3cU0IMhAAAIBmFHa4+gUOjAADg0DAQAgCAQ8OhUQAAkIxN6xE2UF/uFQZCAACQjA3nCOUeB3FoFAAAHBtmhAAAIB07vKAeAyEAAEhH/DlC2U8S4tAoAAA4NMwIAQBAMgrxEzy5J4R/eSA0mUyS9AMAABq5yspKV1fX+rex5fIJuU8SKiwWy195/MaNGz/88EO1usnOLKurq5VKpU6nk7sj94Nerzebzc7OznJ35H4wGAx6vd7NzU3ujtwPJpOpsrLSw8ND7o7cDxaLpbS01MvLS+6O3CelpaUeHh7i18K1RWJiYuvWrSXfbUFBQU5OTseOHSXf8z36qwMhAACAXUNYBgAAHBoGQgAAcGgYCAEAwKFhIAQAAIfWZNOef8WtW7dyc3PbtWtnDQoXFRUdOHDAzc0tNja2KSVI9Xr9sWPHcnJyWrdu3aNHD6HRaDQePHiwpKQkJiamefPm8vZQQkVFRadOnSoqKvL39+/Xr59KpRLaT5w4cfXq1a5du7Zr107eHjaE9PT06upqax6voKDg0KFDHh4esbGxGo1G3r5Jpaio6Nq1a9YfIyIihDBwWlraqVOngoOD+/TpI1/vGsT169ePHTvm5ubWu3dvHx8fxlhZWdn+/fs1Gk1cXJyDBL8lhNToH9TW1gYHB1dVVVVWVh47diw6Opoxdvny5f79+w8YMCAnJ6eqqurHH390cXGRu6cSMJvNPj4+7du3DwsLO3LkSM+ePTdv3mwymQYPHlxdXd26deu9e/fu3bvXOkDau+HDhzPG/Pz8zpw5o1arDx8+7OHhMWfOnG+//TYmJmbXrl3Lly+fNm2a3N2U0rVr17p169a8efMrV64wxs6fPz9w4MBBgwZlZGRYLJZDhw41jb/qvvrqq7///e9du3YVfnzvvfciIyO/+eabF198ceTIkb/88ktcXNzHH38sbycl9Pbbby9fvjwmJkav14eGhq5ZsyYjI6NPnz49evSorKzMzMz8+eefvb295e6mXbFAHSaT6eLFiyaTyd3d/cSJE0LjE0888fLLLwv39urVa926dbL2UTJms/nixYvC7cLCQk9Pz4MHD27evLl9+/Z6vd5isbzxxhujRo2StY8NwmAwREREJCUl3bx509nZ+caNGxaL5eDBgy1atBCeeNNgNpuHDBny0ksvtWnTRmh59NFH582bZ7FYDAZD9+7dk5KSZO2gZL788ssRI0bUbTEajSEhIVu2bLFYLLm5ue7u7pcuXZKpdxI7fvy4h4dHenp63ca//e1vTz75pMViMZvNQ4cOXb58uTyds1s4R/gHSqWyXbt2SuUfXpZt27Y9+uijwr2PPPLI9u3bZeqdxBQKhfVgYLNmzdzc3Kqrq7dv3z527FitVssYmzhx4u7duw0Gg6zdlJ7FYjGbzd7e3rt27YqKigoODmaMxcTEWCyW48ePy907yXz66aehoaGxsbHCjxaLZfv27cInWa1Wjx8/ftu2bbJ2UEoVFRWHDh1KTU01m82MsdTU1IKCgpEjRzLGmjdvPmDAgCbztd2wYUN8fLybm9uRI0fy8vKERuv/UQqFYuLEiU3myd43GAjvorKysqSkpGXLlsKPLVu2zMzMlLdLDeGLL77QarUxMTGZmZnWJ9uqVSuTyZSTkyNv3yT01VdfjRs3rn379pMmTRo1alRmZmarVq2EuxQKRVBQUJN5c7Oyst57770VK1ZYW4qLi6urq5vqJzk3N3fFihWjRo2Kjo7Ozs7Oysry9/e3Vrxq2bJlVlaWvD2USnp6+uXLl2NjY99+++327dtv2LDBYrFkZ2fX/do2pXf2/kBY5i6E+ZA1WKFWq2tra2XtkfR+/PHHOXPmbN261cXFxWg0Wp+scKMpzQgffPBBb2/vkydPrlu3Lj4+3mAwWJ8sY0ytVjeZJ/v888//85//rHuiqAl/kidNmjRlyhTGmMFgmDBhwj/+8Y8xY8bc8c42mSdbVVV169at8+fPOzk57dy5c/LkyePHj6/7tW1KT/a+wYzwLry8vFxcXPLz84Ufc3NzAwMD5e2StI4dOzZx4sSNGzf27t2bMebv7299ssKBl4CAADn7J6k2bdqMHDly8eLFI0eO/PjjjwMCAqwHlxhjeXl5TePJnj179tChQ3v37n322WfXrl2bl5f37LPP6nQ6tVrdJD/J1virRqOJj4//9ddfAwIC8vPzLf9NAubl5TWZJxsYGNinTx8nJyfGWFxcXGlpaUZGRosWLZrkO3vfYCC8u/79+ycnJwu39+/f379/f3n7I6HTp0+PHTs2MTFx0KBBQssdT7Znz55NMopdWlrq5ubWv3//X375pbKykjF25cqVgoKCqKgoubsmgcDAwNWrV/fo0aNHjx7h4eE6na5Hjx46ne6hhx5qqp9kq9TU1KCgoE6dOjHGUlJSGGO1tbX/+c9/msyTjYmJsV4rcvXqVZVK5e/vf8fXdsCAAfJ10C7h8ok7rVq1qrCw8N13342Pjw8ICJg3b97p06fHjRu3cOHCrKysr7/+OjU1tUWLFnJ3UwJVVVUtW7Zs27ZtTEyM0DJmzJjIyMhOnToNGzasQ4cOb7755meffTZ27FhZuymNmpqaYcOGxcXFeXp6Hj16dN++fcePHw8PDx81apRerx8/fvwnn3wydOjQlStXyt1TiW3duvXVV18VLp/Ys2fP448/vnDhwmvXrv3f//3fuXPnhEvQ7N2MGTM8PDyCgoLOnj37zTff7N27t2/fvkuXLt24ceOLL764a9euioqKw4cPy91NaVRXV3fp0mXQoEHdu3f/6KOPYmNj33nnnZSUlNjY2Llz51ZWVq5du/bUqVMPPPCA3D21JxgI7/TNN9+UlZVZf3z88cfd3NxOnjz57bffurm5TZ06VQgZNgF6vT4pKaluS9++fSMjI3NyctavX19YWDhmzJh+/frJ1T1pWSyWHTt2nDp1qry8vHXr1vHx8c2aNWOM6fX6xMTEtLS0qKioxx577P6sZXM/3bhx49ixY5MmTRJ+PHr06Pfff+/l5TV16tSgoCB5+yaVY8eOJScnFxcXBwYGPvroo9Zv6HfffffLL7+EhoY+88wzTePaX0FRUVFiYmJBQUHv3r2tf6empqZu3LhRo9E88cQTDbFSUtOGgRAAABwazhECAIBDw0AIAAAODQMhAAA4tP8PC7swR2LPjVcAAAAASUVORK5CYII=",
"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.3"
},
"kernelspec": {
"name": "julia-1.5",
"display_name": "Julia 1.5.3",
"language": "julia"
}
},
"nbformat": 4
}