{
"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.122979e+01 8.041349e+00\n",
" * time: 0.0014331340789794922\n",
" 1 2.962691e+01 5.194341e+00\n",
" * time: 0.00603795051574707\n",
" 2 2.182118e+01 6.600144e+00\n",
" * time: 0.012025117874145508\n",
" 3 1.563171e+01 4.917492e+00\n",
" * time: 0.01752305030822754\n",
" 4 1.186117e+01 2.149640e+00\n",
" * time: 0.022284984588623047\n",
" 5 1.137053e+01 2.929259e+00\n",
" * time: 0.026194095611572266\n",
" 6 1.053403e+01 1.750955e+00\n",
" * time: 0.031018972396850586\n",
" 7 9.747891e+00 1.552575e+00\n",
" * time: 0.034928083419799805\n",
" 8 9.090514e+00 1.679429e+00\n",
" * time: 0.038886070251464844\n",
" 9 8.779532e+00 1.127845e+00\n",
" * time: 0.04311013221740723\n",
" 10 8.582690e+00 9.808984e-01\n",
" * time: 0.0481419563293457\n",
" 11 8.381509e+00 8.715960e-01\n",
" * time: 0.052304983139038086\n",
" 12 8.223375e+00 6.743384e-01\n",
" * time: 0.05624508857727051\n",
" 13 8.197315e+00 8.439144e-01\n",
" * time: 0.05989503860473633\n",
" 14 8.139930e+00 7.023806e-01\n",
" * time: 0.06377601623535156\n",
" 15 8.074470e+00 6.752596e-01\n",
" * time: 0.09373712539672852\n",
" 16 8.036873e+00 5.538479e-01\n",
" * time: 0.09615302085876465\n",
" 17 7.995625e+00 5.610557e-01\n",
" * time: 0.0984959602355957\n",
" 18 7.942597e+00 4.213812e-01\n",
" * time: 0.10098791122436523\n",
" 19 7.923442e+00 8.019301e-01\n",
" * time: 0.10300397872924805\n",
" 20 7.907011e+00 8.261580e-01\n",
" * time: 0.10527896881103516\n",
" 21 7.870397e+00 4.219163e-01\n",
" * time: 0.10834407806396484\n",
" 22 7.855176e+00 5.505108e-01\n",
" * time: 0.11079597473144531\n",
" 23 7.815045e+00 4.594445e-01\n",
" * time: 0.11294412612915039\n",
" 24 7.799625e+00 4.662650e-01\n",
" * time: 0.1157069206237793\n",
" 25 7.777751e+00 3.029118e-01\n",
" * time: 0.11850595474243164\n",
" 26 7.756934e+00 3.360724e-01\n",
" * time: 0.12151503562927246\n",
" 27 7.749661e+00 3.981501e-01\n",
" * time: 0.12383699417114258\n",
" 28 7.748521e+00 3.674228e-01\n",
" * time: 0.12640595436096191\n",
" 29 7.743930e+00 4.164525e-01\n",
" * time: 0.12885308265686035\n",
" 30 7.724561e+00 2.842537e-01\n",
" * time: 0.1319410800933838\n",
" 31 7.721135e+00 5.095137e-01\n",
" * time: 0.13440299034118652\n",
" 32 7.708546e+00 4.804166e-01\n",
" * time: 0.13692307472229004\n",
" 33 7.691247e+00 2.350169e-01\n",
" * time: 0.13994407653808594\n",
" 34 7.682490e+00 3.041951e-01\n",
" * time: 0.14243292808532715\n",
" 35 7.679673e+00 2.418369e-01\n",
" * time: 0.14464807510375977\n",
" 36 7.674242e+00 1.786314e-01\n",
" * time: 0.14735102653503418\n",
" 37 7.664293e+00 1.583925e-01\n",
" * time: 0.1502389907836914\n",
" 38 7.661325e+00 2.861231e-01\n",
" * time: 0.15380597114562988\n",
" 39 7.652962e+00 2.294544e-01\n",
" * time: 0.1591949462890625\n",
" 40 7.647011e+00 1.843652e-01\n",
" * time: 0.16402602195739746\n",
" 41 7.640637e+00 1.971624e-01\n",
" * time: 0.16797494888305664\n",
" 42 7.636014e+00 1.372626e-01\n",
" * time: 0.17357897758483887\n",
" 43 7.633162e+00 1.062178e-01\n",
" * time: 0.1783759593963623\n",
" 44 7.631505e+00 1.006395e-01\n",
" * time: 0.1817479133605957\n",
" 45 7.631221e+00 1.342639e-01\n",
" * time: 0.1857929229736328\n",
" 46 7.629368e+00 8.469784e-02\n",
" * time: 0.20447993278503418\n",
" 47 7.628685e+00 9.651779e-02\n",
" * time: 0.20651507377624512\n",
" 48 7.627694e+00 8.109871e-02\n",
" * time: 0.20859003067016602\n",
" 49 7.626311e+00 4.839437e-02\n",
" * time: 0.21058011054992676\n",
" 50 7.625166e+00 5.113374e-02\n",
" * time: 0.21263599395751953\n",
" 51 7.624780e+00 7.262042e-02\n",
" * time: 0.21499013900756836\n",
" 52 7.624503e+00 6.358642e-02\n",
" * time: 0.2175159454345703\n",
" 53 7.624236e+00 8.465746e-02\n",
" * time: 0.21996712684631348\n",
" 54 7.624051e+00 7.048943e-02\n",
" * time: 0.22196507453918457\n",
" 55 7.623163e+00 7.060252e-02\n",
" * time: 0.22400403022766113\n",
" 56 7.622416e+00 3.661848e-02\n",
" * time: 0.22653603553771973\n",
" 57 7.621963e+00 7.270198e-02\n",
" * time: 0.22852802276611328\n",
" 58 7.621737e+00 8.689250e-02\n",
" * time: 0.23093008995056152\n",
" 59 7.621224e+00 5.613937e-02\n",
" * time: 0.23328304290771484\n",
" 60 7.620919e+00 6.026545e-02\n",
" * time: 0.2357640266418457\n",
" 61 7.620240e+00 5.781363e-02\n",
" * time: 0.23807501792907715\n",
" 62 7.619508e+00 4.033903e-02\n",
" * time: 0.24078893661499023\n",
" 63 7.619377e+00 6.979910e-02\n",
" * time: 0.24297404289245605\n",
" 64 7.618908e+00 6.716091e-02\n",
" * time: 0.24570703506469727\n",
" 65 7.618355e+00 3.383349e-02\n",
" * time: 0.24873113632202148\n",
" 66 7.617652e+00 3.791594e-02\n",
" * time: 0.251068115234375\n",
" 67 7.617143e+00 4.593312e-02\n",
" * time: 0.253277063369751\n",
" 68 7.616893e+00 8.083180e-02\n",
" * time: 0.255648136138916\n",
" 69 7.616190e+00 4.484075e-02\n",
" * time: 0.2586350440979004\n",
" 70 7.616032e+00 7.532133e-02\n",
" * time: 0.2624959945678711\n",
" 71 7.615566e+00 4.377557e-02\n",
" * time: 0.2677431106567383\n",
" 72 7.615047e+00 5.773214e-02\n",
" * time: 0.27170705795288086\n",
" 73 7.614853e+00 5.599879e-02\n",
" * time: 0.275130033493042\n",
" 74 7.614841e+00 7.548428e-02\n",
" * time: 0.27919602394104004\n",
" 75 7.614441e+00 7.990188e-02\n",
" * time: 0.28328704833984375\n",
" 76 7.614056e+00 7.432898e-02\n",
" * time: 0.2869751453399658\n",
" 77 7.613328e+00 6.679307e-02\n",
" * time: 0.29041504859924316\n",
" 78 7.612781e+00 5.976445e-02\n",
" * time: 0.30577611923217773\n",
" 79 7.612486e+00 5.926440e-02\n",
" * time: 0.3080129623413086\n",
" 80 7.612055e+00 6.092114e-02\n",
" * time: 0.31112003326416016\n",
" 81 7.611601e+00 4.362699e-02\n",
" * time: 0.3140439987182617\n",
" 82 7.611309e+00 3.585826e-02\n",
" * time: 0.31603407859802246\n",
" 83 7.610903e+00 5.130359e-02\n",
" * time: 0.31803107261657715\n",
" 84 7.610518e+00 3.544783e-02\n",
" * time: 0.3205270767211914\n",
" 85 7.610339e+00 3.879526e-02\n",
" * time: 0.3227670192718506\n",
" 86 7.609995e+00 3.580910e-02\n",
" * time: 0.325747013092041\n",
" 87 7.609913e+00 7.756345e-02\n",
" * time: 0.3281710147857666\n",
" 88 7.609824e+00 5.500128e-02\n",
" * time: 0.3305840492248535\n",
" 89 7.609462e+00 5.539547e-02\n",
" * time: 0.33307504653930664\n",
" 90 7.609008e+00 5.098483e-02\n",
" * time: 0.3355889320373535\n",
" 91 7.608538e+00 3.716341e-02\n",
" * time: 0.33866405487060547\n",
" 92 7.608480e+00 5.906393e-02\n",
" * time: 0.3411281108856201\n",
" 93 7.608180e+00 4.713079e-02\n",
" * time: 0.3442881107330322\n",
" 94 7.607972e+00 5.505148e-02\n",
" * time: 0.3466460704803467\n",
" 95 7.607676e+00 5.182515e-02\n",
" * time: 0.3488340377807617\n",
" 96 7.607194e+00 3.502500e-02\n",
" * time: 0.35094308853149414\n",
" 97 7.607184e+00 5.021192e-02\n",
" * time: 0.35303807258605957\n",
" 98 7.606871e+00 5.456627e-02\n",
" * time: 0.3560910224914551\n",
" 99 7.606568e+00 3.041682e-02\n",
" * time: 0.3591940402984619\n",
" 100 7.606377e+00 5.002653e-02\n",
" * time: 0.3623509407043457\n",
" 101 7.606121e+00 3.774674e-02\n",
" * time: 0.3657550811767578\n",
" 102 7.606108e+00 5.366885e-02\n",
" * time: 0.3691270351409912\n",
" 103 7.605826e+00 3.891657e-02\n",
" * time: 0.3742949962615967\n",
" 104 7.605723e+00 4.869222e-02\n",
" * time: 0.37764906883239746\n",
" 105 7.605583e+00 7.508928e-02\n",
" * time: 0.38048505783081055\n",
" 106 7.605452e+00 5.287392e-02\n",
" * time: 0.3834221363067627\n",
" 107 7.605135e+00 3.823169e-02\n",
" * time: 0.38662004470825195\n",
" 108 7.604837e+00 2.049159e-02\n",
" * time: 0.39923095703125\n",
" 109 7.604801e+00 5.590117e-02\n",
" * time: 0.4014449119567871\n",
" 110 7.604780e+00 4.867948e-02\n",
" * time: 0.40383100509643555\n",
" 111 7.604596e+00 4.126217e-02\n",
" * time: 0.40619611740112305\n",
" 112 7.604281e+00 4.214552e-02\n",
" * time: 0.4082491397857666\n",
" 113 7.604201e+00 5.422412e-02\n",
" * time: 0.4102311134338379\n",
" 114 7.604004e+00 6.177909e-02\n",
" * time: 0.4121971130371094\n",
" 115 7.603796e+00 6.492113e-02\n",
" * time: 0.4141850471496582\n",
" 116 7.603532e+00 4.775545e-02\n",
" * time: 0.416748046875\n",
" 117 7.603377e+00 5.655964e-02\n",
" * time: 0.41919398307800293\n",
" 118 7.603273e+00 5.289166e-02\n",
" * time: 0.4216439723968506\n",
" 119 7.603077e+00 5.213989e-02\n",
" * time: 0.4237360954284668\n",
" 120 7.603010e+00 6.035143e-02\n",
" * time: 0.42572808265686035\n",
" 121 7.602684e+00 4.559793e-02\n",
" * time: 0.42824792861938477\n",
" 122 7.602468e+00 4.612345e-02\n",
" * time: 0.4308459758758545\n",
" 123 7.602223e+00 4.130215e-02\n",
" * time: 0.4330711364746094\n",
" 124 7.601883e+00 4.218979e-02\n",
" * time: 0.43625593185424805\n",
" 125 7.601691e+00 2.352147e-02\n",
" * time: 0.4395010471343994\n",
" 126 7.601578e+00 2.887224e-02\n",
" * time: 0.4419851303100586\n",
" 127 7.601397e+00 2.530797e-02\n",
" * time: 0.44446301460266113\n",
" 128 7.601246e+00 2.756150e-02\n",
" * time: 0.4475901126861572\n",
" 129 7.601098e+00 3.283176e-02\n",
" * time: 0.45073699951171875\n",
" 130 7.600982e+00 3.638977e-02\n",
" * time: 0.4540529251098633\n",
" 131 7.600751e+00 2.043097e-02\n",
" * time: 0.45861005783081055\n",
" 132 7.600673e+00 3.205637e-02\n",
" * time: 0.4622640609741211\n",
" 133 7.600524e+00 3.257767e-02\n",
" * time: 0.4670069217681885\n",
" 134 7.600368e+00 4.152013e-02\n",
" * time: 0.4708230495452881\n",
" 135 7.600159e+00 3.068077e-02\n",
" * time: 0.4757080078125\n",
" 136 7.599995e+00 2.019210e-02\n",
" * time: 0.47872400283813477\n",
" 137 7.599949e+00 3.422529e-02\n",
" * time: 0.48124098777770996\n",
" 138 7.599757e+00 2.716672e-02\n",
" * time: 0.48394107818603516\n",
" 139 7.599620e+00 3.177704e-02\n",
" * time: 0.49555110931396484\n",
" 140 7.599512e+00 2.554000e-02\n",
" * time: 0.4980599880218506\n",
" 141 7.599491e+00 3.421710e-02\n",
" * time: 0.5007600784301758\n",
" 142 7.599362e+00 3.093265e-02\n",
" * time: 0.5038719177246094\n",
" 143 7.599191e+00 4.511819e-02\n",
" * time: 0.506385087966919\n",
" 144 7.599054e+00 3.398795e-02\n",
" * time: 0.5094509124755859\n",
" 145 7.598900e+00 2.125321e-02\n",
" * time: 0.5125961303710938\n",
" 146 7.598814e+00 4.055044e-02\n",
" * time: 0.5152161121368408\n",
" 147 7.598655e+00 3.750311e-02\n",
" * time: 0.5175681114196777\n",
" 148 7.598470e+00 2.800670e-02\n",
" * time: 0.5205399990081787\n",
" 149 7.598397e+00 3.532578e-02\n",
" * time: 0.5229840278625488\n",
" 150 7.598195e+00 2.672998e-02\n",
" * time: 0.5253291130065918\n",
" 151 7.598021e+00 3.399095e-02\n",
" * time: 0.5278270244598389\n",
" 152 7.597696e+00 2.404643e-02\n",
" * time: 0.5311551094055176\n",
" 153 7.597640e+00 4.313169e-02\n",
" * time: 0.5336430072784424\n",
" 154 7.597565e+00 3.324135e-02\n",
" * time: 0.5361270904541016\n",
" 155 7.597423e+00 3.367817e-02\n",
" * time: 0.5386481285095215\n",
" 156 7.597405e+00 4.493830e-02\n",
" * time: 0.5411150455474854\n",
" 157 7.597230e+00 4.822960e-02\n",
" * time: 0.5436179637908936\n",
" 158 7.597003e+00 3.339382e-02\n",
" * time: 0.5469570159912109\n",
" 159 7.596835e+00 3.311032e-02\n",
" * time: 0.5498371124267578\n",
" 160 7.596641e+00 3.784571e-02\n",
" * time: 0.5519981384277344\n",
" 161 7.596476e+00 2.975651e-02\n",
" * time: 0.5548510551452637\n",
" 162 7.596466e+00 3.859681e-02\n",
" * time: 0.557002067565918\n",
" 163 7.596405e+00 3.852983e-02\n",
" * time: 0.5595309734344482\n",
" 164 7.596259e+00 3.761790e-02\n",
" * time: 0.561967134475708\n",
" 165 7.596253e+00 4.994300e-02\n",
" * time: 0.5644609928131104\n",
" 166 7.596093e+00 4.381489e-02\n",
" * time: 0.5665950775146484\n",
" 167 7.595878e+00 2.992704e-02\n",
" * time: 0.5700831413269043\n",
" 168 7.595690e+00 3.570843e-02\n",
" * time: 0.5722470283508301\n",
" 169 7.595489e+00 3.992696e-02\n",
" * time: 0.5833749771118164\n",
" 170 7.595335e+00 3.227723e-02\n",
" * time: 0.5861029624938965\n",
" 171 7.595157e+00 3.180085e-02\n",
" * time: 0.588231086730957\n",
" 172 7.595055e+00 2.444863e-02\n",
" * time: 0.5906410217285156\n",
" 173 7.594920e+00 1.698692e-02\n",
" * time: 0.5940279960632324\n",
" 174 7.594808e+00 2.074471e-02\n",
" * time: 0.5965960025787354\n",
" 175 7.594718e+00 2.181492e-02\n",
" * time: 0.5991120338439941\n",
" 176 7.594632e+00 2.167301e-02\n",
" * time: 0.6022260189056396\n",
" 177 7.594516e+00 9.238920e-03\n",
" * time: 0.6053309440612793\n",
" 178 7.594486e+00 2.221289e-02\n",
" * time: 0.6078379154205322\n",
" 179 7.594391e+00 2.262643e-02\n",
" * time: 0.610522985458374\n",
" 180 7.594366e+00 2.374135e-02\n",
" * time: 0.6130239963531494\n",
" 181 7.594310e+00 1.990355e-02\n",
" * time: 0.6155319213867188\n",
" 182 7.594264e+00 2.243086e-02\n",
" * time: 0.6180441379547119\n",
" 183 7.594248e+00 2.789877e-02\n",
" * time: 0.6205379962921143\n",
" 184 7.594231e+00 2.883251e-02\n",
" * time: 0.6230831146240234\n",
" 185 7.594153e+00 2.499256e-02\n",
" * time: 0.6264820098876953\n",
" 186 7.594082e+00 2.225329e-02\n",
" * time: 0.629209041595459\n",
" 187 7.594033e+00 2.538817e-02\n",
" * time: 0.6318800449371338\n",
" 188 7.594014e+00 1.940278e-02\n",
" * time: 0.6344740390777588\n",
" 189 7.593959e+00 2.115205e-02\n",
" * time: 0.6379349231719971\n",
" 190 7.593910e+00 1.127434e-02\n",
" * time: 0.6413710117340088\n",
" 191 7.593904e+00 1.684426e-02\n",
" * time: 0.6440620422363281\n",
" 192 7.593871e+00 1.607887e-02\n",
" * time: 0.6465799808502197\n",
" 193 7.593842e+00 1.681870e-02\n",
" * time: 0.649183988571167\n",
" 194 7.593808e+00 1.131391e-02\n",
" * time: 0.6524429321289062\n",
" 195 7.593778e+00 9.706294e-03\n",
" * time: 0.6557459831237793\n",
" 196 7.593753e+00 9.213728e-03\n",
" * time: 0.6583030223846436\n",
" 197 7.593729e+00 1.640549e-02\n",
" * time: 0.6610541343688965\n",
" 198 7.593727e+00 1.756481e-02\n",
" * time: 0.6640870571136475\n",
" 199 7.593699e+00 1.183837e-02\n",
" * time: 0.6672670841217041\n",
" 200 7.593675e+00 1.020887e-02\n",
" * time: 0.6791789531707764\n",
" 201 7.593641e+00 1.286421e-02\n",
" * time: 0.6817100048065186\n",
" 202 7.593602e+00 1.505342e-02\n",
" * time: 0.6842019557952881\n",
" 203 7.593580e+00 1.100117e-02\n",
" * time: 0.6875169277191162\n",
" 204 7.593543e+00 9.048546e-03\n",
" * time: 0.6906700134277344\n",
" 205 7.593537e+00 2.092083e-02\n",
" * time: 0.6931331157684326\n",
" 206 7.593531e+00 1.725994e-02\n",
" * time: 0.6956090927124023\n",
" 207 7.593504e+00 1.538315e-02\n",
" * time: 0.6981029510498047\n",
" 208 7.593502e+00 1.934268e-02\n",
" * time: 0.7005951404571533\n",
" 209 7.593463e+00 1.756513e-02\n",
" * time: 0.703963041305542\n",
" 210 7.593426e+00 1.343708e-02\n",
" * time: 0.7064769268035889\n",
" 211 7.593404e+00 1.335173e-02\n",
" * time: 0.7089779376983643\n",
" 212 7.593382e+00 1.235428e-02\n",
" * time: 0.7114889621734619\n",
" 213 7.593377e+00 1.183915e-02\n",
" * time: 0.7139899730682373\n",
" 214 7.593353e+00 1.006407e-02\n",
" * time: 0.7165710926055908\n",
" 215 7.593324e+00 1.223770e-02\n",
" * time: 0.7192609310150146\n",
" 216 7.593300e+00 1.016437e-02\n",
" * time: 0.7224149703979492\n",
" 217 7.593283e+00 1.331974e-02\n",
" * time: 0.7255489826202393\n",
" 218 7.593274e+00 1.318231e-02\n",
" * time: 0.7280299663543701\n",
" 219 7.593251e+00 9.042707e-03\n",
" * time: 0.7312679290771484\n",
" 220 7.593242e+00 6.715157e-03\n",
" * time: 0.7346489429473877\n",
" 221 7.593226e+00 4.483438e-03\n",
" * time: 0.7372019290924072\n",
" 222 7.593214e+00 9.274588e-03\n",
" * time: 0.7397110462188721\n",
" 223 7.593200e+00 8.244464e-03\n",
" * time: 0.7422099113464355\n",
" 224 7.593197e+00 1.429226e-02\n",
" * time: 0.7447259426116943\n",
" 225 7.593190e+00 9.599830e-03\n",
" * time: 0.7472310066223145\n",
" 226 7.593164e+00 7.083517e-03\n",
" * time: 0.7499289512634277\n",
" 227 7.593149e+00 1.274406e-02\n",
" * time: 0.7524240016937256\n",
" 228 7.593140e+00 1.613032e-02\n",
" * time: 0.7549431324005127\n",
" 229 7.593116e+00 8.387166e-03\n",
" * time: 0.7585489749908447\n",
" 230 7.593098e+00 7.035766e-03\n",
" * time: 0.7610471248626709\n",
" 231 7.593092e+00 1.133280e-02\n",
" * time: 0.7728760242462158\n",
" 232 7.593074e+00 9.500377e-03\n",
" * time: 0.7755060195922852\n",
" 233 7.593072e+00 1.284885e-02\n",
" * time: 0.7780959606170654\n",
" 234 7.593052e+00 9.812644e-03\n",
" * time: 0.7807259559631348\n",
" 235 7.593028e+00 9.473956e-03\n",
" * time: 0.7839131355285645\n",
" 236 7.593008e+00 6.991877e-03\n",
" * time: 0.7870900630950928\n",
" 237 7.592992e+00 7.683299e-03\n",
" * time: 0.7896599769592285\n",
" 238 7.592979e+00 1.506843e-02\n",
" * time: 0.7922711372375488\n",
" 239 7.592959e+00 7.767108e-03\n",
" * time: 0.7956540584564209\n",
" 240 7.592941e+00 8.363681e-03\n",
" * time: 0.798285961151123\n",
" 241 7.592928e+00 1.692358e-02\n",
" * time: 0.8007810115814209\n",
" 242 7.592906e+00 8.651789e-03\n",
" * time: 0.8039300441741943\n",
" 243 7.592901e+00 1.205663e-02\n",
" * time: 0.8064329624176025\n",
" 244 7.592883e+00 9.318605e-03\n",
" * time: 0.8095600605010986\n",
" 245 7.592863e+00 9.330215e-03\n",
" * time: 0.8122439384460449\n",
" 246 7.592843e+00 9.139900e-03\n",
" * time: 0.8147399425506592\n",
" 247 7.592828e+00 1.274444e-02\n",
" * time: 0.8172080516815186\n",
" 248 7.592808e+00 1.301942e-02\n",
" * time: 0.8197169303894043\n",
" 249 7.592787e+00 6.569008e-03\n",
" * time: 0.8229320049285889\n",
" 250 7.592771e+00 1.179640e-02\n",
" * time: 0.8254570960998535\n",
" 251 7.592757e+00 1.293004e-02\n",
" * time: 0.8281240463256836\n",
" 252 7.592753e+00 1.154983e-02\n",
" * time: 0.8306069374084473\n",
" 253 7.592736e+00 1.110358e-02\n",
" * time: 0.8330719470977783\n",
" 254 7.592723e+00 1.204782e-02\n",
" * time: 0.8355491161346436\n",
" 255 7.592718e+00 1.395266e-02\n",
" * time: 0.8380270004272461\n",
" 256 7.592702e+00 1.817329e-02\n",
" * time: 0.8404891490936279\n",
" 257 7.592700e+00 1.595070e-02\n",
" * time: 0.8430891036987305\n",
" 258 7.592680e+00 9.572685e-03\n",
" * time: 0.8461079597473145\n",
" 259 7.592661e+00 1.268456e-02\n",
" * time: 0.8484909534454346\n",
" 260 7.592642e+00 1.087397e-02\n",
" * time: 0.8509681224822998\n",
" 261 7.592633e+00 1.201636e-02\n",
" * time: 0.8534460067749023\n",
" 262 7.592621e+00 1.444585e-02\n",
" * time: 0.855942964553833\n",
" 263 7.592606e+00 7.205131e-03\n",
" * time: 0.8680529594421387\n",
" 264 7.592592e+00 1.108725e-02\n",
" * time: 0.8705160617828369\n",
" 265 7.592578e+00 1.022238e-02\n",
" * time: 0.873661994934082\n",
" 266 7.592563e+00 7.992463e-03\n",
" * time: 0.8769099712371826\n",
" 267 7.592560e+00 1.313093e-02\n",
" * time: 0.879418134689331\n",
" 268 7.592543e+00 9.243734e-03\n",
" * time: 0.88254714012146\n",
" 269 7.592535e+00 1.114048e-02\n",
" * time: 0.8850550651550293\n",
" 270 7.592523e+00 8.722124e-03\n",
" * time: 0.887566089630127\n",
" 271 7.592520e+00 1.131925e-02\n",
" * time: 0.890265941619873\n",
" 272 7.592520e+00 1.123008e-02\n",
" * time: 0.8927750587463379\n",
" 273 7.592508e+00 7.344187e-03\n",
" * time: 0.8952910900115967\n",
" 274 7.592504e+00 9.714913e-03\n",
" * time: 0.8977820873260498\n",
" 275 7.592492e+00 1.299734e-02\n",
" * time: 0.9003260135650635\n",
" 276 7.592487e+00 7.806256e-03\n",
" * time: 0.9028329849243164\n",
" 277 7.592482e+00 8.780856e-03\n",
" * time: 0.9054291248321533\n",
" 278 7.592467e+00 6.286919e-03\n",
" * time: 0.908668041229248\n",
" 279 7.592461e+00 6.726336e-03\n",
" * time: 0.9118211269378662\n",
" 280 7.592455e+00 7.431509e-03\n",
" * time: 0.9143400192260742\n",
" 281 7.592449e+00 6.141836e-03\n",
" * time: 0.916823148727417\n",
" 282 7.592442e+00 6.426480e-03\n",
" * time: 0.9199850559234619\n",
" 283 7.592431e+00 3.286473e-03\n",
" * time: 0.9232759475708008\n",
" 284 7.592426e+00 3.485459e-03\n",
" * time: 0.9264049530029297\n",
" 285 7.592425e+00 8.584793e-03\n",
" * time: 0.9289059638977051\n",
" 286 7.592421e+00 6.630196e-03\n",
" * time: 0.9314169883728027\n",
" 287 7.592414e+00 6.464129e-03\n",
" * time: 0.9345629215240479\n",
" 288 7.592413e+00 8.853763e-03\n",
" * time: 0.9372251033782959\n",
" 289 7.592408e+00 5.543389e-03\n",
" * time: 0.9397470951080322\n",
" 290 7.592397e+00 5.780869e-03\n",
" * time: 0.9422569274902344\n",
" 291 7.592391e+00 3.792298e-03\n",
" * time: 0.9453740119934082\n",
" 292 7.592389e+00 4.982756e-03\n",
" * time: 0.9479129314422607\n",
" 293 7.592383e+00 3.892845e-03\n",
" * time: 0.9595239162445068\n",
" 294 7.592380e+00 5.351848e-03\n",
" * time: 0.9620330333709717\n",
" 295 7.592379e+00 5.014534e-03\n",
" * time: 0.9645059108734131\n",
" 296 7.592376e+00 3.359004e-03\n",
" * time: 0.9677109718322754\n",
" 297 7.592373e+00 2.971451e-03\n",
" * time: 0.9703159332275391\n",
" 298 7.592370e+00 4.519380e-03\n",
" * time: 0.9727880954742432\n",
" 299 7.592366e+00 3.067496e-03\n",
" * time: 0.9759190082550049\n",
" 300 7.592365e+00 2.264150e-03\n",
" * time: 0.9791109561920166\n",
" 301 7.592364e+00 3.396637e-03\n",
" * time: 0.9818410873413086\n",
" 302 7.592362e+00 2.864486e-03\n",
" * time: 0.98465895652771\n",
" 303 7.592362e+00 3.623735e-03\n",
" * time: 0.9871630668640137\n",
" 304 7.592361e+00 3.182557e-03\n",
" * time: 0.9896650314331055\n",
" 305 7.592359e+00 4.812012e-03\n",
" * time: 0.9921650886535645\n",
" 306 7.592358e+00 2.192644e-03\n",
" * time: 0.9953370094299316\n",
" 307 7.592357e+00 2.777160e-03\n",
" * time: 0.9978220462799072\n",
" 308 7.592355e+00 2.433774e-03\n",
" * time: 1.0005090236663818\n",
" 309 7.592354e+00 3.600501e-03\n",
" * time: 1.0029981136322021\n",
" 310 7.592353e+00 2.668612e-03\n",
" * time: 1.0055041313171387\n",
" 311 7.592352e+00 2.998975e-03\n",
" * time: 1.007997989654541\n",
" 312 7.592351e+00 1.791824e-03\n",
" * time: 1.0112099647521973\n",
" 313 7.592350e+00 3.367376e-03\n",
" * time: 1.013746976852417\n",
" 314 7.592349e+00 1.922486e-03\n",
" * time: 1.017091989517212\n",
" 315 7.592348e+00 2.153743e-03\n",
" * time: 1.0202760696411133\n",
" 316 7.592347e+00 1.436843e-03\n",
" * time: 1.023442029953003\n",
" 317 7.592347e+00 2.600490e-03\n",
" * time: 1.0259771347045898\n",
" 318 7.592346e+00 2.584244e-03\n",
" * time: 1.0285530090332031\n",
" 319 7.592346e+00 1.942423e-03\n",
" * time: 1.0312299728393555\n",
" 320 7.592345e+00 2.375374e-03\n",
" * time: 1.0337319374084473\n",
" 321 7.592345e+00 2.430385e-03\n",
" * time: 1.0363221168518066\n",
" 322 7.592345e+00 2.365220e-03\n",
" * time: 1.038841962814331\n",
" 323 7.592344e+00 2.608496e-03\n",
" * time: 1.041327953338623\n",
" 324 7.592343e+00 1.141218e-03\n",
" * time: 1.0539300441741943\n",
" 325 7.592342e+00 1.268935e-03\n",
" * time: 1.0565690994262695\n",
" 326 7.592342e+00 1.910160e-03\n",
" * time: 1.0592501163482666\n",
" 327 7.592342e+00 2.041865e-03\n",
" * time: 1.0619499683380127\n",
" 328 7.592341e+00 1.164659e-03\n",
" * time: 1.0651249885559082\n",
" 329 7.592341e+00 2.831403e-03\n",
" * time: 1.0676429271697998\n",
" 330 7.592341e+00 2.235623e-03\n",
" * time: 1.0707969665527344\n",
" 331 7.592341e+00 1.857724e-03\n",
" * time: 1.0733470916748047\n",
" 332 7.592341e+00 2.221466e-03\n",
" * time: 1.0759799480438232\n",
" 333 7.592341e+00 2.085523e-03\n",
" * time: 1.0787150859832764\n",
" 334 7.592340e+00 1.091504e-03\n",
" * time: 1.081860065460205\n",
" 335 7.592340e+00 1.782067e-03\n",
" * time: 1.084378957748413\n",
" 336 7.592339e+00 1.398026e-03\n",
" * time: 1.0869081020355225\n",
" 337 7.592339e+00 1.263535e-03\n",
" * time: 1.0901720523834229\n",
" 338 7.592339e+00 8.668438e-04\n",
" * time: 1.0934779644012451\n",
" 339 7.592338e+00 1.195870e-03\n",
" * time: 1.0959889888763428\n",
" 340 7.592338e+00 7.669958e-04\n",
" * time: 1.099133014678955\n",
" 341 7.592338e+00 1.121864e-03\n",
" * time: 1.1016099452972412\n",
" 342 7.592338e+00 1.307612e-03\n",
" * time: 1.1041460037231445\n",
" 343 7.592338e+00 1.041482e-03\n",
" * time: 1.1073849201202393\n",
" 344 7.592338e+00 6.369906e-04\n",
" * time: 1.1106879711151123\n",
" 345 7.592338e+00 8.703047e-04\n",
" * time: 1.113171100616455\n",
" 346 7.592338e+00 1.051278e-03\n",
" * time: 1.1156530380249023\n",
" 347 7.592338e+00 7.105433e-04\n",
" * time: 1.1181590557098389\n",
" 348 7.592337e+00 5.603412e-04\n",
" * time: 1.1206979751586914\n",
" 349 7.592337e+00 8.054896e-04\n",
" * time: 1.1233630180358887\n",
" 350 7.592337e+00 7.510595e-04\n",
" * time: 1.1261470317840576\n",
" 351 7.592337e+00 1.051305e-03\n",
" * time: 1.1288061141967773\n",
" 352 7.592337e+00 7.131994e-04\n",
" * time: 1.1321649551391602\n",
" 353 7.592337e+00 1.006199e-03\n",
" * time: 1.1348590850830078\n",
" 354 7.592337e+00 7.986403e-04\n",
" * time: 1.1374831199645996\n",
" 355 7.592337e+00 7.457208e-04\n",
" * time: 1.149277925491333\n",
" 356 7.592337e+00 7.778563e-04\n",
" * time: 1.151871919631958\n",
" 357 7.592337e+00 3.907001e-04\n",
" * time: 1.1551151275634766\n",
" 358 7.592337e+00 7.493911e-04\n",
" * time: 1.157745122909546\n",
" 359 7.592337e+00 5.842344e-04\n",
" * time: 1.1603729724884033\n",
" 360 7.592337e+00 6.562276e-04\n",
" * time: 1.1630771160125732\n",
" 361 7.592337e+00 7.506908e-04\n",
" * time: 1.1663539409637451\n",
" 362 7.592337e+00 8.352800e-04\n",
" * time: 1.1689810752868652\n",
" 363 7.592337e+00 9.472751e-04\n",
" * time: 1.1717050075531006\n",
" 364 7.592336e+00 5.720017e-04\n",
" * time: 1.174854040145874\n",
" 365 7.592336e+00 6.304469e-04\n",
" * time: 1.1773531436920166\n",
" 366 7.592336e+00 5.805973e-04\n",
" * time: 1.1798839569091797\n",
" 367 7.592336e+00 6.345697e-04\n",
" * time: 1.182507038116455\n",
" 368 7.592336e+00 7.522445e-04\n",
" * time: 1.185157060623169\n",
" 369 7.592336e+00 6.243511e-04\n",
" * time: 1.1886589527130127\n",
" 370 7.592336e+00 1.197960e-03\n",
" * time: 1.1913409233093262\n",
" 371 7.592336e+00 7.431643e-04\n",
" * time: 1.1945970058441162\n",
" 372 7.592336e+00 7.106787e-04\n",
" * time: 1.1970889568328857\n",
" 373 7.592335e+00 7.888861e-04\n",
" * time: 1.1997339725494385\n",
" 374 7.592335e+00 5.318967e-04\n",
" * time: 1.2030761241912842\n",
" 375 7.592335e+00 7.995488e-04\n",
" * time: 1.2055590152740479\n",
" 376 7.592335e+00 8.509128e-04\n",
" * time: 1.208078145980835\n",
" 377 7.592335e+00 9.756695e-04\n",
" * time: 1.2106060981750488\n",
" 378 7.592335e+00 9.066597e-04\n",
" * time: 1.2130990028381348\n",
" 379 7.592335e+00 9.898951e-04\n",
" * time: 1.2157230377197266\n",
" 380 7.592335e+00 9.154441e-04\n",
" * time: 1.2184100151062012\n",
" 381 7.592335e+00 6.177933e-04\n",
" * time: 1.2215609550476074\n",
" 382 7.592335e+00 7.644241e-04\n",
" * time: 1.2240569591522217\n",
" 383 7.592335e+00 7.735228e-04\n",
" * time: 1.2266380786895752\n",
" 384 7.592334e+00 6.502734e-04\n",
" * time: 1.2298400402069092\n",
" 385 7.592334e+00 1.077262e-03\n",
" * time: 1.2324090003967285\n",
" 386 7.592334e+00 1.493294e-03\n",
" * time: 1.2444329261779785\n",
" 387 7.592334e+00 7.828781e-04\n",
" * time: 1.2477269172668457\n",
" 388 7.592334e+00 1.047423e-03\n",
" * time: 1.2504401206970215\n",
" 389 7.592334e+00 1.151099e-03\n",
" * time: 1.2529981136322021\n",
" 390 7.592334e+00 6.647514e-04\n",
" * time: 1.2561349868774414\n",
" 391 7.592334e+00 9.258147e-04\n",
" * time: 1.2586381435394287\n",
" 392 7.592333e+00 1.467664e-03\n",
" * time: 1.261228084564209\n",
" 393 7.592333e+00 1.117447e-03\n",
" * time: 1.264456033706665\n",
" 394 7.592333e+00 1.002816e-03\n",
" * time: 1.2680439949035645\n",
" 395 7.592333e+00 9.770372e-04\n",
" * time: 1.270678997039795\n",
" 396 7.592333e+00 1.433352e-03\n",
" * time: 1.2733500003814697\n",
" 397 7.592333e+00 1.068466e-03\n",
" * time: 1.2767040729522705\n",
" 398 7.592333e+00 1.444333e-03\n",
" * time: 1.2793560028076172\n",
" 399 7.592332e+00 1.232882e-03\n",
" * time: 1.2827229499816895\n",
" 400 7.592332e+00 8.380539e-04\n",
" * time: 1.28603196144104\n",
" 401 7.592332e+00 1.055473e-03\n",
" * time: 1.2885429859161377\n",
" 402 7.592332e+00 9.731639e-04\n",
" * time: 1.2911529541015625\n",
" 403 7.592332e+00 1.863298e-03\n",
" * time: 1.2937560081481934\n",
" 404 7.592331e+00 1.681105e-03\n",
" * time: 1.2964191436767578\n",
" 405 7.592331e+00 1.911270e-03\n",
" * time: 1.298996925354004\n",
" 406 7.592331e+00 1.853985e-03\n",
" * time: 1.3014941215515137\n",
" 407 7.592330e+00 2.457718e-03\n",
" * time: 1.303995132446289\n",
" 408 7.592330e+00 2.942402e-03\n",
" * time: 1.3065130710601807\n",
" 409 7.592330e+00 3.361880e-03\n",
" * time: 1.3091521263122559\n",
" 410 7.592329e+00 2.767590e-03\n",
" * time: 1.3117690086364746\n",
" 411 7.592329e+00 3.703688e-03\n",
" * time: 1.3143870830535889\n",
" 412 7.592328e+00 2.786403e-03\n",
" * time: 1.3168599605560303\n",
" 413 7.592328e+00 2.792871e-03\n",
" * time: 1.3193731307983398\n",
" 414 7.592326e+00 1.512939e-03\n",
" * time: 1.3225090503692627\n",
" 415 7.592326e+00 1.408154e-03\n",
" * time: 1.3249740600585938\n",
" 416 7.592325e+00 1.366084e-03\n",
" * time: 1.3283140659332275\n",
" 417 7.592324e+00 1.533627e-03\n",
" * time: 1.3406391143798828\n",
" 418 7.592324e+00 1.312936e-03\n",
" * time: 1.3439459800720215\n",
" 419 7.592324e+00 2.541340e-03\n",
" * time: 1.3464200496673584\n",
" 420 7.592323e+00 1.340366e-03\n",
" * time: 1.3495259284973145\n",
" 421 7.592323e+00 1.901944e-03\n",
" * time: 1.3520350456237793\n",
" 422 7.592323e+00 1.669735e-03\n",
" * time: 1.3552250862121582\n",
" 423 7.592322e+00 2.349817e-03\n",
" * time: 1.3577120304107666\n",
" 424 7.592322e+00 1.591539e-03\n",
" * time: 1.3604021072387695\n",
" 425 7.592321e+00 1.115719e-03\n",
" * time: 1.363508939743042\n",
" 426 7.592321e+00 1.260761e-03\n",
" * time: 1.3659911155700684\n",
" 427 7.592321e+00 1.413226e-03\n",
" * time: 1.3684749603271484\n",
" 428 7.592320e+00 1.146712e-03\n",
" * time: 1.3710739612579346\n",
" 429 7.592320e+00 1.210393e-03\n",
" * time: 1.3736140727996826\n",
" 430 7.592320e+00 1.099374e-03\n",
" * time: 1.376986026763916\n",
" 431 7.592319e+00 1.825105e-03\n",
" * time: 1.3794970512390137\n",
" 432 7.592319e+00 1.076177e-03\n",
" * time: 1.3826360702514648\n",
" 433 7.592319e+00 1.374964e-03\n",
" * time: 1.3852241039276123\n",
" 434 7.592319e+00 8.315456e-04\n",
" * time: 1.388519048690796\n",
" 435 7.592318e+00 1.413261e-03\n",
" * time: 1.3912811279296875\n",
" 436 7.592318e+00 1.220718e-03\n",
" * time: 1.3937859535217285\n",
" 437 7.592318e+00 1.689296e-03\n",
" * time: 1.3962929248809814\n",
" 438 7.592317e+00 1.386973e-03\n",
" * time: 1.3994309902191162\n",
" 439 7.592317e+00 1.359832e-03\n",
" * time: 1.4019320011138916\n",
" 440 7.592317e+00 9.856739e-04\n",
" * time: 1.404400110244751\n",
" 441 7.592316e+00 1.252170e-03\n",
" * time: 1.407094955444336\n",
" 442 7.592316e+00 1.045679e-03\n",
" * time: 1.410236120223999\n",
" 443 7.592316e+00 2.139240e-03\n",
" * time: 1.412708044052124\n",
" 444 7.592316e+00 1.061380e-03\n",
" * time: 1.4158539772033691\n",
" 445 7.592315e+00 9.521794e-04\n",
" * time: 1.4184820652008057\n",
" 446 7.592315e+00 1.260532e-03\n",
" * time: 1.4217870235443115\n",
" 447 7.592315e+00 1.434296e-03\n",
" * time: 1.4334900379180908\n",
" 448 7.592315e+00 1.395546e-03\n",
" * time: 1.436094045639038\n",
" 449 7.592315e+00 1.023787e-03\n",
" * time: 1.438730001449585\n",
" 450 7.592315e+00 8.422402e-04\n",
" * time: 1.441188097000122\n",
" 451 7.592315e+00 8.789574e-04\n",
" * time: 1.443681001663208\n",
" 452 7.592315e+00 8.822153e-04\n",
" * time: 1.4461770057678223\n",
" 453 7.592315e+00 9.217619e-04\n",
" * time: 1.448734998703003\n",
" 454 7.592314e+00 6.674760e-04\n",
" * time: 1.45194411277771\n",
" 455 7.592314e+00 8.688512e-04\n",
" * time: 1.4545409679412842\n",
" 456 7.592314e+00 8.099588e-04\n",
" * time: 1.4570159912109375\n",
" 457 7.592314e+00 8.769903e-04\n",
" * time: 1.4595179557800293\n",
" 458 7.592314e+00 6.599960e-04\n",
" * time: 1.4626719951629639\n",
" 459 7.592314e+00 5.023326e-04\n",
" * time: 1.465951919555664\n",
" 460 7.592314e+00 9.451352e-04\n",
" * time: 1.468644142150879\n",
" 461 7.592314e+00 6.629701e-04\n",
" * time: 1.4711670875549316\n",
" 462 7.592314e+00 6.815671e-04\n",
" * time: 1.4736671447753906\n",
" 463 7.592314e+00 1.197415e-03\n",
" * time: 1.4761900901794434\n",
" 464 7.592314e+00 8.774134e-04\n",
" * time: 1.4793930053710938\n",
" 465 7.592314e+00 1.506129e-03\n",
" * time: 1.482116937637329\n",
" 466 7.592313e+00 7.846167e-04\n",
" * time: 1.4855120182037354\n",
" 467 7.592313e+00 5.622763e-04\n",
" * time: 1.4886629581451416\n",
" 468 7.592313e+00 1.004569e-03\n",
" * time: 1.4911730289459229\n",
" 469 7.592313e+00 6.267797e-04\n",
" * time: 1.4936730861663818\n",
" 470 7.592313e+00 6.173425e-04\n",
" * time: 1.4961750507354736\n",
" 471 7.592313e+00 8.763243e-04\n",
" * time: 1.4987709522247314\n",
" 472 7.592313e+00 7.011771e-04\n",
" * time: 1.5020229816436768\n",
" 473 7.592313e+00 1.022897e-03\n",
" * time: 1.5045020580291748\n",
" 474 7.592313e+00 1.062944e-03\n",
" * time: 1.5069921016693115\n",
" 475 7.592312e+00 6.552263e-04\n",
" * time: 1.5094640254974365\n",
" 476 7.592312e+00 1.344156e-03\n",
" * time: 1.511965036392212\n",
" 477 7.592312e+00 7.983643e-04\n",
" * time: 1.5145149230957031\n",
" 478 7.592312e+00 1.206719e-03\n",
" * time: 1.5170960426330566\n",
" 479 7.592312e+00 1.091780e-03\n",
" * time: 1.5287930965423584\n",
" 480 7.592312e+00 9.127935e-04\n",
" * time: 1.5314481258392334\n",
" 481 7.592312e+00 8.307703e-04\n",
" * time: 1.5339291095733643\n",
" 482 7.592312e+00 1.670341e-03\n",
" * time: 1.5364179611206055\n",
" 483 7.592311e+00 1.393137e-03\n",
" * time: 1.5389249324798584\n",
" 484 7.592311e+00 1.084014e-03\n",
" * time: 1.5421299934387207\n",
" 485 7.592311e+00 6.979647e-04\n",
" * time: 1.5453450679779053\n",
" 486 7.592311e+00 9.840887e-04\n",
" * time: 1.5480151176452637\n",
" 487 7.592311e+00 8.268222e-04\n",
" * time: 1.5511541366577148\n",
" 488 7.592311e+00 5.040654e-04\n",
" * time: 1.5543699264526367\n",
" 489 7.592311e+00 1.622945e-03\n",
" * time: 1.5569190979003906\n",
" 490 7.592310e+00 1.203215e-03\n",
" * time: 1.5596251487731934\n",
" 491 7.592310e+00 1.095909e-03\n",
" * time: 1.562330961227417\n",
" 492 7.592310e+00 1.383502e-03\n",
" * time: 1.5648219585418701\n",
" 493 7.592310e+00 1.868239e-03\n",
" * time: 1.567331075668335\n",
" 494 7.592310e+00 1.071701e-03\n",
" * time: 1.5699360370635986\n",
" 495 7.592309e+00 1.228139e-03\n",
" * time: 1.5725600719451904\n",
" 496 7.592309e+00 1.447945e-03\n",
" * time: 1.575186014175415\n",
" 497 7.592309e+00 1.472284e-03\n",
" * time: 1.5778329372406006\n",
" 498 7.592308e+00 9.794771e-04\n",
" * time: 1.5804810523986816\n",
" 499 7.592308e+00 8.326715e-04\n",
" * time: 1.5830590724945068\n",
" 500 7.592308e+00 1.244760e-03\n",
" * time: 1.5862109661102295\n",
" 501 7.592308e+00 1.803163e-03\n",
" * time: 1.5887501239776611\n",
" 502 7.592308e+00 1.425418e-03\n",
" * time: 1.591357946395874\n",
" 503 7.592307e+00 1.291239e-03\n",
" * time: 1.5940399169921875\n",
" 504 7.592307e+00 1.051159e-03\n",
" * time: 1.5971579551696777\n",
" 505 7.592307e+00 9.475327e-04\n",
" * time: 1.6003789901733398\n",
" 506 7.592306e+00 6.802385e-04\n",
" * time: 1.6035120487213135\n",
" 507 7.592306e+00 1.711117e-03\n",
" * time: 1.6060740947723389\n",
" 508 7.592306e+00 9.968570e-04\n",
" * time: 1.609386920928955\n",
" 509 7.592306e+00 1.051930e-03\n",
" * time: 1.611889123916626\n",
" 510 7.592306e+00 8.129281e-04\n",
" * time: 1.62391996383667\n",
" 511 7.592306e+00 1.062359e-03\n",
" * time: 1.6267590522766113\n",
" 512 7.592306e+00 1.312330e-03\n",
" * time: 1.6294050216674805\n",
" 513 7.592306e+00 7.440486e-04\n",
" * time: 1.6320641040802002\n",
" 514 7.592305e+00 8.143635e-04\n",
" * time: 1.6346640586853027\n",
" 515 7.592305e+00 1.374633e-03\n",
" * time: 1.6372430324554443\n",
" 516 7.592305e+00 1.303545e-03\n",
" * time: 1.6405699253082275\n",
" 517 7.592305e+00 1.334270e-03\n",
" * time: 1.6431519985198975\n",
" 518 7.592305e+00 7.748057e-04\n",
" * time: 1.6462829113006592\n",
" 519 7.592305e+00 1.047324e-03\n",
" * time: 1.648771047592163\n",
" 520 7.592305e+00 9.067938e-04\n",
" * time: 1.6513340473175049\n",
" 521 7.592305e+00 7.493153e-04\n",
" * time: 1.6538851261138916\n",
" 522 7.592305e+00 1.057233e-03\n",
" * time: 1.6566050052642822\n",
" 523 7.592305e+00 9.580617e-04\n",
" * time: 1.659924030303955\n",
" 524 7.592305e+00 6.010324e-04\n",
" * time: 1.6632399559020996\n",
" 525 7.592305e+00 1.082501e-03\n",
" * time: 1.6658079624176025\n",
" 526 7.592305e+00 1.143050e-03\n",
" * time: 1.668442964553833\n",
" 527 7.592305e+00 8.288707e-04\n",
" * time: 1.6710829734802246\n",
" 528 7.592304e+00 6.077063e-04\n",
" * time: 1.6736769676208496\n",
" 529 7.592304e+00 4.346339e-04\n",
" * time: 1.6761720180511475\n",
" 530 7.592304e+00 5.527290e-04\n",
" * time: 1.6786811351776123\n",
" 531 7.592304e+00 6.257253e-04\n",
" * time: 1.6811821460723877\n",
" 532 7.592304e+00 4.734995e-04\n",
" * time: 1.6845099925994873\n",
" 533 7.592304e+00 3.386831e-04\n",
" * time: 1.6878421306610107\n",
" 534 7.592304e+00 5.460467e-04\n",
" * time: 1.6903760433197021\n",
" 535 7.592304e+00 4.113025e-04\n",
" * time: 1.6929590702056885\n",
" 536 7.592304e+00 1.013812e-03\n",
" * time: 1.6954660415649414\n",
" 537 7.592304e+00 8.324561e-04\n",
" * time: 1.6987159252166748\n",
" 538 7.592304e+00 5.762140e-04\n",
" * time: 1.7020759582519531\n",
" 539 7.592304e+00 8.159976e-04\n",
" * time: 1.7046599388122559\n",
" 540 7.592304e+00 6.416366e-04\n",
" * time: 1.7071759700775146\n",
" 541 7.592304e+00 1.174705e-03\n",
" * time: 1.7188971042633057\n",
" 542 7.592304e+00 5.963562e-04\n",
" * time: 1.7213780879974365\n",
" 543 7.592304e+00 5.636157e-04\n",
" * time: 1.7238540649414062\n",
" 544 7.592304e+00 7.771219e-04\n",
" * time: 1.726341962814331\n",
" 545 7.592304e+00 1.097493e-03\n",
" * time: 1.7288200855255127\n",
" 546 7.592303e+00 1.174788e-03\n",
" * time: 1.7314491271972656\n",
" 547 7.592303e+00 6.213433e-04\n",
" * time: 1.7347290515899658\n",
" 548 7.592303e+00 7.352676e-04\n",
" * time: 1.7371830940246582\n",
" 549 7.592303e+00 7.122682e-04\n",
" * time: 1.739673137664795\n",
" 550 7.592303e+00 3.879021e-04\n",
" * time: 1.742811918258667\n",
" 551 7.592303e+00 6.751406e-04\n",
" * time: 1.745361089706421\n",
" 552 7.592303e+00 5.760892e-04\n",
" * time: 1.7487330436706543\n",
" 553 7.592303e+00 3.866982e-04\n",
" * time: 1.75205397605896\n",
" 554 7.592303e+00 3.531694e-04\n",
" * time: 1.7552099227905273\n",
" 555 7.592303e+00 4.800337e-04\n",
" * time: 1.757699966430664\n",
" 556 7.592303e+00 6.236399e-04\n",
" * time: 1.7602031230926514\n",
" 557 7.592303e+00 4.359017e-04\n",
" * time: 1.7628891468048096\n",
" 558 7.592303e+00 4.449010e-04\n",
" * time: 1.7656080722808838\n",
" 559 7.592303e+00 5.581119e-04\n",
" * time: 1.768226146697998\n",
" 560 7.592303e+00 4.148960e-04\n",
" * time: 1.7707419395446777\n",
" 561 7.592303e+00 3.383739e-04\n",
" * time: 1.7732279300689697\n",
" 562 7.592303e+00 4.381269e-04\n",
" * time: 1.7757360935211182\n",
" 563 7.592303e+00 3.415827e-04\n",
" * time: 1.7790210247039795\n",
" 564 7.592303e+00 5.836848e-04\n",
" * time: 1.7816951274871826\n",
" 565 7.592303e+00 3.883394e-04\n",
" * time: 1.7842481136322021\n",
" 566 7.592303e+00 5.997273e-04\n",
" * time: 1.786790132522583\n",
" 567 7.592303e+00 4.828633e-04\n",
" * time: 1.7899670600891113\n",
" 568 7.592303e+00 6.059820e-04\n",
" * time: 1.7925410270690918\n",
" 569 7.592303e+00 6.807460e-04\n",
" * time: 1.7951021194458008\n",
" 570 7.592302e+00 6.052967e-04\n",
" * time: 1.7984230518341064\n",
" 571 7.592302e+00 8.221753e-04\n",
" * time: 1.8008971214294434\n",
" 572 7.592302e+00 7.858349e-04\n",
" * time: 1.812593936920166\n",
" 573 7.592302e+00 7.776014e-04\n",
" * time: 1.8151681423187256\n",
" 574 7.592302e+00 7.448790e-04\n",
" * time: 1.8176419734954834\n",
" 575 7.592302e+00 5.515709e-04\n",
" * time: 1.820786952972412\n",
" 576 7.592302e+00 5.781691e-04\n",
" * time: 1.8233790397644043\n",
" 577 7.592302e+00 3.993094e-04\n",
" * time: 1.826643943786621\n",
" 578 7.592302e+00 6.794209e-04\n",
" * time: 1.8292970657348633\n",
" 579 7.592302e+00 8.990385e-04\n",
" * time: 1.8318049907684326\n",
" 580 7.592302e+00 8.154840e-04\n",
" * time: 1.83430814743042\n",
" 581 7.592302e+00 8.498985e-04\n",
" * time: 1.8367979526519775\n",
" 582 7.592302e+00 6.293984e-04\n",
" * time: 1.840116024017334\n",
" 583 7.592302e+00 3.552261e-04\n",
" * time: 1.8434350490570068\n",
" 584 7.592302e+00 5.479099e-04\n",
" * time: 1.8459219932556152\n",
" 585 7.592302e+00 3.402267e-04\n",
" * time: 1.8484070301055908\n",
" 586 7.592302e+00 3.828023e-04\n",
" * time: 1.8509209156036377\n",
" 587 7.592302e+00 2.221537e-04\n",
" * time: 1.8540871143341064\n",
" 588 7.592302e+00 6.386296e-04\n",
" * time: 1.8566789627075195\n",
" 589 7.592302e+00 4.635184e-04\n",
" * time: 1.859377145767212\n",
" 590 7.592301e+00 2.196817e-04\n",
" * time: 1.8625340461730957\n",
" 591 7.592301e+00 5.021191e-04\n",
" * time: 1.8650119304656982\n",
" 592 7.592301e+00 3.984708e-04\n",
" * time: 1.8675010204315186\n",
" 593 7.592301e+00 2.675558e-04\n",
" * time: 1.8700709342956543\n",
" 594 7.592301e+00 3.526317e-04\n",
" * time: 1.872736930847168\n",
" 595 7.592301e+00 6.571966e-04\n",
" * time: 1.8755300045013428\n",
" 596 7.592301e+00 5.496870e-04\n",
" * time: 1.8781640529632568\n",
" 597 7.592301e+00 4.318822e-04\n",
" * time: 1.8807640075683594\n",
" 598 7.592301e+00 3.808529e-04\n",
" * time: 1.8834640979766846\n"
]
},
{
"output_type": "execute_result",
"data": {
"text/plain": "Plot{Plots.GRBackend() n=1}",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVxU5f4H8GdWFhFEFtnEBfclUdBcEVFzTREzSVvIyureNFNvWOmvtFua6a17rUxTizK1bDE0N3IrF1QUF3IpQUR22QWG2ZjfH6c7l/T5omc6eBjm837xx8wzZw5nzszwcJ7zOd9HYbFYGAAAgKNSyr0BAAAAckJHCAAADg0dIQAAODR0hAAA4NDQEQIAgENDRwgAAA4NHSEAADg0dIQAAODQ0BECAIBDQ0cIAAAOTf0Xn//blXR3bx8XFxdJtqYRqq2tVSgUCoVC7g25FywWi8ViUSod5d+j2tpavNgmCS+2gXg4ae64zOXLJTqdSdRqKysrAgNV7dq1s3W7/irFX6w1On3m3ybPWyTV1gAAQKMV09n/jsuEhiacPXtD5IrTJk4s27Ztm21b9dc5yj9NAAAAXH91aBQAAMBKoWB2dyoJR4QAAODQcEQIAACSsSFdaLHIfAiJjhAAACSFoVEAAAA7giNCAACQjh2GZdARAgCAZGypQCL3OUIMjQIAgENDRwgAAA4NQ6MAACAZGy6ol3tkFB0hAABIx4ZzhLKHazA0CgAAjd2VK1dOnz5tNBpvf6igoODkyZPl5eXWloqKitTU1PT09LqzSpSUlJw8ebKoqKjuc69du5aSkoKOEAAApKMQ/1Ov2traRx99NCoq6vnnn+/SpUtmZmbdRz/44IPu3bvPmzevQ4cOu3btYowtXLgwODj4mWeeiYyMHDBggNDzbd26tWPHjvPnz+/cufPnn38uPHf+/Pn9+vWbM2cOOkIAAJCM4r+nCUX91GPPnj1HjhxJS0s7fvz4mDFj3njjDetDRUVF8fHx+/fv//nnnz/88MNZs2bV1taOHDkyKysrJSUlIyOjWbNmb7/9tsFgeOGFFzZt2nTo0KHExMQXX3yxqqrq7NmzGzZsOHPmzOHDh9ERAgBA4/X1119PmTLF3d2dMfbUU09t3brVOuCZmJjYs2fP++67jzEWExNTXFx86tSpoUOHCgtrNJrw8PCCgoKDBw9qtdpRo0YxxgYNGtSqVau9e/d+/fXXY8eO9ff3ZzhHCAAAElKIV/8hYVZWlnXy+nbt2lVXV9+4ccP6UPv27YXbarW6devW165dsz6xrKxsy5YtMTExddcgrOTatWt/eq6Erx8AABzdXZz2u/0p2dnZX375Zd22sWPHenp6Msaqq6udnJyERmdnZ8ZYVVWVcLe6ulqr1Vqf4uLiUl1dLdw2GAzTpk0bOnTo5MmT//Of/9y+WN3noiMEAACZFRQU7Ny5s27LkCFDhI7Qz8+vpKREaCwuLhZahLt+fn6XLl2yPqW4uFh4yGg0xsbGOjs7b9iw4ZY1MMaKior8/PzqNqIjBAAAKYk/IGRhYWG3HBFahYeHHz58eP78+Yyxw4cP9+jRw8XFxfrQsmXLzGazSqXKzc3Nzs7u1auX2WyOi4vT6XTbtm1Tq9WMsbCwsAsXLpSVlbVo0UKn050+ffrDDz9UKBSrV68W1qNkjOXk5MyZM2fkyJHTp09PSUlhjJnN5vfee2/06NGPPvpoWlqabfsCAAAcjUIh/jRhvV3njBkzDh06tGLFih07dixYsOCll15ijL322mvPP/98REREmzZtnnnmmaSkpLi4uEceeaRVq1Zz587duXPn2LFjExIS1q5du3379pCQkNGjR8fFxSUlJT355JMDBgzo2bPn1KlTc3JyXnvttd27dyvz8/P79++vUCj+8Y9/jBs3rrKykjH27rvvfvrppy+//HLPnj2HDRtWVlZ2j3YhAABAHf7+/gcPHrxw4cL69esXL148Y8YMxlh4ePigQYMYYzt37mzRosWqVauGDBkiHOHdd999zz77bE5OTkZGRkZGRm5uLmPsiy++6NGjx6pVq9q2bfvNN98wxlxdXX/++efi4uKPPvpI8fe//72srGzjxo3W32o2m4ODgxMSEkaMGMEYi4qKiomJeeGFF7ibOH3m3ybPW9TwuwIAAGQW09n/jsv0H/Tl+fNFd1ysLrPp/OgHirZt22brdv1VysOHDw8dOvT//u//nnrqqc2bNzPGCgoKcnNzBw4cKCwxcODAU6dOybV9AABgR2wYGVXIXWxUmZWVtWTJkoCAgFGjRi1cuHD58uUFBQXOzs6urq7CEl5eXgUFBdTzdTrdvdpUAACQk9lsvuMy0paVuTfUrq6usbGxzz33HGNMrVbPmzdv0qRJBoNByOEwxqqqqpo3b04933p5BwAANG1Cp9D0KNu0aSPUmGGM+fv7l5SUBAQEKJVKa2HTq1evtm7dmny+ErVpAADgv2woui33QaFy+vTpO3bsMJlMjLHvv/++f//+zZo1Gz9+/Nq1axlj+fn5iYmJsbGxMm8mAADYA/Ejo/IPjqpnzJixe/fuzp07e3h46PX677//njG2bNmysWPH7t27NycnJy4uLjw8XO7tBAAAaBBqZ2fnxMTErKwso9HYrl07Yaizc+fOv/3226VLl3x8fHx9feXeSAAAsBPiD/FkPyT8o8RacHDwLQ+oVKru3bvf8+0BAAB71gjO+YmFqAsAADg0FN0GAADJCBfUi3yOzIeQ6AjBPvx3SmoJ1iRFs2jkF138AyJ+qb2NUEETYEMMVPYPKoZGAQDAoeGIEAAApGOHYRl0hAAAIBkbimjLXnQbHSEAAEjGDg8IcY4QAAAcG44IQTZ0EJTzALmw2LQnsSJq/aI2kjFG/TdMjf1Ilq8TORYlycYAcNhwSCj3Bw8dIQAASMeGc4QNtCV3DUOjAADg0HBECAAAklFIOOB/r6AjBAAAySgU4i+HkLsnxNAoAAA4NBwRgmiiy36KDGrW8h6w1BIrqSVWQi0vOjUqTRFS6n9kBfG/KLW8UiliPfRK+L+UDC2Ivj5a1OLQtCA1CgAAjswO+0EMjQIAgGPDESEAAEjGllqjch8ToiMEAADp2OHYKDpCEF1FjIyPSJRnMZs5TzCb+AvX8hZmjJnNxPIiN5IMy4gkNixDhWJUKiL/ouKsSKUWsTBjTEUmcfg7QUEsTyH/8Zf77yBIyIbrCGX/AOAcIQAAODQcEQIAgHTssNYoOkIAAJCOHZ4jxNAoAAA4NBwRAgCAZBQKFN2GxkBkUTCy6hgRsKSCl7VEUJMKfJqM/KSmycRpNxnNolZiNhIbU8tf3kKlTOm9xm+nqppReUkyHUoEO4kgqFrDWZ7byBhTa1R3v5J6fqmSirCKqQPH6DSp7H8fwQYKphB7XeAdl8/Jyfnkk0/Kysqio6MjIyPrPqTT6dauXXvlypW+ffs++uijSqXSZDKlpaWlpKRUVlbOmTOHMXbjxo0NGzbUfdbYsWPVanViYqJwF0OjAADQeJWVlfXr16+oqKhDhw5Tpkyx9l6CmJiYnTt33nfffe+99158fDxjbM+ePdHR0QkJCW+88YawjNlsLv2va9euLViwoKam5uzZs2vWrBEacUQIAADSsSEsU6+EhIQuXbp88MEHjDFnZ+dly5ZNmDBBeOj06dNHjx7Ny8tzdXUdOnRonz59XnvttXHjxo0bNy45OXn06NHCYn5+fsuWLRNur127tkePHn379k1PT+/YsaPQjiNCAACQjHCOUNRP/R3nL7/8MmLECOH2iBEjkpOTDQaD9aFBgwa5uroyxjp16uTl5XXq1Kn6N2/Dhg1PP/20cDsnJ+ftt99eu3YtOkIAAJCM4o9qo+LUs8K8vDwfHx/htq+vr8ViycvLu/0h4VHrQ1y//vpramrqtGnTGGMeHh79+vUzGo179+7F0CgAAMjs2LFjY8aMqdvy4Ycftm/fnjGm0WhMJpPQaDQaGWNOTk7CXa1Wa31IeFSr1dbzW9avXx8dHS30nWPGjLH+RnSEdoysgil2Ilyq5icvvcmIVCdjzGTgtxvJdn4Q1KjntJMLEysn06TExpORVyo2SrUTgUmxtUNVairwKaJd48RPh2qJdo2Waid+KdVObDz1osiJgskRK6RMGzGbzhGGhITMnj27bov1UC8wMDAnJ0e4nZOTo9ForA8FBAQcOnRIuG2xWHJzcwMDA6lfYTAYNm7c+MUXX9z+EDpCAACQjoIsMU8+Q8F8fX1vOSK0mjhx4pIlSxYtWqTVar/66qvx48erVKrLly8bDIbx48fPnTs3KysrODh4//79KpWqX79+1G9JTEx0cnKynm6srq4WTi5aLBZ0hAAA0HhNmjRp9erVgwcPbtu27cGDB5OSkhhjH374YVFR0aZNm1588cWIiIiIiIg9e/a8++67Go0mNzd3woQJVVVVlZWV4eHhwcHB3333HWNsw4YNM2bMUKn+GPaIjY0tKysLDAxMS0tT/MVZZqbP/NvkeYv++ksFG2BotJ6VY2iU2dPQqLgL8DE0KpeYzv53XGbMxO8uXioRtdrqytSBfXO2bdtGLWA2mw8dOlRRUREREdGyZUvGWHZ2tslkatu2LWMsNTU1PT29d+/eISEhjDG9Xp+WlmZ9rrOzc/fu3RljZ86cad++vbu7u9B+8+bNkydPFhcXBwQE4IgQAACk0wBFt1UqVVRUVN2WoKAg6+3evXv37t3betfJySksLOz2lYSGhta927x5c+s6cfkEAAA4NBwR2gFyRniRNUKpSdvNxCgiNRpp0POXN+hN/PYaYj1Eu76Gsx4Db7yUMWYkNkbs0CjZTuw0audTMQGqLKfYUURqaJQ7eql1JoZGnflffCdyeWqIlb8erRO1kfz1qKgap8ROw5BpY3bH6wI5T5F7HiZ0hAAAIBmF+P88ZP9PBUOjAADg0HBECAAA0rHDGerREQIAgGRsOEco+9goOsJGRGwohp4IV9z1fFQORa/jh1/EttcQ7QadFGEZsdcRUu1m4vpCfrNo1MS8SrHXCxLLa3j5FOp6QTIs48Jvd3bhr8fJhf+mOBPr0brwP7FUGEdDThRMXI/IXw1DiOZesmWG+obZkruHc4QAAODQ0BECAIBDw9AoAABIx4brCOUeG1WbTKabN29a7zdr1kyYz8lsNl++fNnX19fb21u+zQMAAGhYyuPHj/v4+IT8l1D29Pfff+/SpUtsbGzXrl1feeUVuTcSAADsgxCWEfUje1pGzRjr1KnThQsX6rbGx8dHR0e/++67eXl5PXv2nDJlSp8+fWTawiaITIcSpdGodCg1C4TYFGhNNdVu5LbrqsSuhyq9RqRPeaXXuFNSMFtKqRE7mZpNgkJ9dUWuhpqTlgpGiiq9RqZGnfh7XkukPfWu/HYXV3FRXmfyzeKvv5ZKnzpxm5naQqRJyVNAnHdR9mG6psDe9uEfH5CKigpnZ2dhULSqqmr79u0XL15kjPn7+0+YMGHz5s3oCAEAoElSMsauXLnSqVMnDw+PCRMmFBYW5uTk1NbWtmvXTlgiJCQkKyuLen6tVJdZAQCA/VPYRN5tVnbv3j0vLy8/Pz8vL89kMr3wwguVlZVardY6ja+rq2tFRQX1fL1ef682FQAA5HQ3Rz42nCOUux9kyhYtWnh5eTHGWrRosWDBgp9++qlVq1Y1NTU6nU5YoqSkxM/Pj3q+i4vLPdpSAACQlZI6p23n/vSqCgoK3N3d/fz8/P39jx07JjQeO3YMJwgBAOBu2HBAKHu6Rr1+/XqlUhkSEpKenr5w4cIXXnhBpVL9/e9/nzt37qpVq06ePHnmzJmtW7fKu5X2ixsQJWuHEsFFI5G1oybCJdObVfwUaHWluPYaIjWqI1KmemJ7uOlQRqReTcQswWRqlJpQl9jJVJS3oVGDQkqiOCk5wS8vNaomJsLVEu1ORI7XqNcQ7VRqlHqz+OuhctH0m0W8W8REwWpiJl+iNikKk/418vdroqkDAwO/+OKLjRs3+vj4rFy5curUqYyxBQsWODs7L1myxNvbe//+/Z6ennJvJwAA2AMbim7L3XGqR48ePXr06FtaVSrVvHnz5s2bJ8s2AQAA3DOoNQoAAJJRMGZ/tUZl/v0AANCU2OE5wqaZhQUAALhLOCKUhqjJ5cnaoUQAksrmUWlMHZH2rLpJpEMrDPx2Yj06In1KpVWpGqcGKjXKyxyaxM5ET8w4T9YCJf6HlazmBfERITeH+L1UmpRIjRJ7mEiNGg1UO/HJlCzKy22m06EE8s0imtW8IwFMc/8X2TJDvdz7EB0hAABIxpaSaXL3hBgaBQAAh4aOEAAAHBqGRgEAQDK2nCNsmC25ezgiBAAAh4YjQqmIKJZIZeeoDJ6eqABJpTerbvJToJXl/PYqIjVKpU91YtOhOhE1RRmxH0zUjPNEQNFEzERP0Trx/yl0ba7ltjs587871JtYTbwpVMFVKk6qUvP/e1ZxU6NG/ouiIrj0zqRmlhfXLrqmKLeVjlYoiEitqOXJhUWmSR2WgokOy8g+HyE6QgAAkI4NF9TL/b8EOkIAAJBOI5hoVyycIwQAgEYtISGhQ4cOvr6+zzzzTE1NTd2Hjh8/3rdv35YtW44aNeratWuMsWvXrs2fPz8qKioyMtK62IQJE8L/a9asWULj6tWr27Vr5+fnh44QAAAko7BJPStMS0t78cUXN27ceOnSpd9///3tt9+2PmQwGCZNmvTMM89kZ2f37NkzLi6OMVZaWurs7Dx8+PAzZ87UXcmcOXPWrFmzZs2a2bNnM8aSk5MXLVr0/fffnz9/HkOj4pCl1IgaUdy5dk1GceERaiJcqmQaFYqpLCPCMkSOQ2zJNLKdyINQuRJuZINKGNF5Df7yVIkyZ08nbnvbzvzJOFu1duO2lxXpuO3paSXcdl1VFbedqlKmMhNhGd7+oXaC6FAMUa+OCr9Qbxa1PDUltdiwjFJJTcArYpZjJTVdL/lb+c12NzzYaG3YsGHy5Mn9+/dnjC1cuPDxxx9fsmSJ8NCPP/7YrFmzmTNnMsZef/11Hx+f33//PTQ0NDQ0NDk5+d133627nq5du4aFhVnvrlu37rHHHgsNDWUYGgUAgMbs0qVLvXr1Em736tUrLy+voqLC+pDQkzHGmjdv3q5du8uXL1PrefTRR3v06PHEE09kZmbeslocEQIAgGRsqDWqUDCDwVBUVFS30dPTU6VSMcaKi4ubN28uNLq7uzPGioqKhBslJSVubv8bm/Hw8LhlJVbLli3r0aOH2WxeuXLliBEjzp07V1JSYl0tOkIAAJCMDZVlmILt27evS5cuddsOHjzYo0cPxpiXl5f1EFC44eXlJdxt2bLl1atXrU8pLy+3PnSLhx9+WLixbt06f3//5OTkli1bWleLjhAAAGQ2ZsyYbdu2cR/q2LHjr7/+KtxOS0vz9fX18PCwPpSQkCDcrqqqyszM7NixY/2/SKVSabVao9FYd7U4RwgAAJJSiPyp15NPPrl169bz589XVVUtW7ZsxowZjLGNGzeuXr36wQcfLCkp+fLLL81m8zvvvNOnT58uXbqYTKaMjIzc3Nza2tqMjIzr169nZWXt27evurq6oqJi0aJFZrP5/vvvf/LJJ7/44ovffvutoqICR4QiEbFRKvbGje0Z9fwMHjXRbjVRSo2aULeKKqVGpEOpUm3UxLlSpUbJQl+87KKJCi4SwUjqHdEQc9K6NON/F9p0bsFt73gffwQmN/Mmt70wm58OLcyp5LaTkVoTFYDk7B/ubL2MMbOKaCd2JlUCjSRmnmp68XrSoUR0Vs1/XWRdOt5+4DYyxhQK/mYSiVTHjY0qFKJLptW/fGho6NKlS0ePHl1VVTV+/PhFixYxxrKyssrKypycnL799tvnn3/+2WefDQ0NFY4O8/PzR44cyRjz8fEZOXJku3btPvjgg3/84x+XL1/WarVhYWE7d+5s0aJFREREfHx8RESEXq9HRwgAANJpgBnqn3vuueeee65uy6uvvircGDJkSFpaWt2HgoKC0tPTb1nD6dOnb1/t3Llz586dyzA0CgAADg5HhAAAIB0bim7LDR0hAABIRsEUCnubfgJDowAA4NBwRMhH1hQl2qmMIrdipEHPT1fWEKlLaoLcamriXKI2KZXqpFKsVO1TekJdooKlyMqW3MAktZFU4Vaq+itVm5TamRUlem57WVENt728mN8utqAr1U7FCrgB0VozPyKr1opMbxIpTaWKv/OpECC5HiX/xaqIRCZ3FmJGR4I1Wv7yal67WiOiMCmjI6zU3mzyYVIbLqiXfZ+gIwQAAEnZ2cgohkYBAMCx4YgQAAAkY0vR7QbalLuGjhAAACRjW9FteWFoFAAAHBqOCClErI6qNUpMvc0NQFI1POkAoUTpUCLVSZW1pLaTerFaZyqzx28XVePUUMN/UQZi46nUqNHA/+ez4Dq/5ufZo3nc9uz0cm57RRk/ZZqTUcFtryIKxhqJqC31v7aGFxBVEmUznYianM4u/HeKquFJfRKoXLGomeIZY3odfz1qDf/D4OTE/0Q5EZ9MrROnnfq4qjXE3wQVv50+yJH78Keh2XBBvdy7BB0hAABIx5ZzhDL3hOgIAQBAMnZ4QIhzhAAA4NhwRAgAAJKx4fIJ2UvLOHxHKFEpNapwF3fuWSqEQk+ESy1P5EeoEmgiS6ZRL8rFTcNt9w9uzm13Jua8Lcji51MyL5Xe3lhZLmLqY0bPJWvm7zOyZBqVSLrmVMZtp+rJUYknqs4cXaCLKvQlYmJe38Bm3PZWQfx3kNrIolz+bMNkDIoK0ZCl1/jtauKTr6/hf9KoqaGdeO1OzsR80U78dmrWX6qeHEXuvkBS9jY2iqFRAABwaA5/RAgAANKxpeh2w2zJ3UNHCAAAklEoyAH8ep4iLwyNAgCAQ8MRIQAASMcOLyR09I6QmomUTI0SNaXIOWZ508ZS2Tm92BSo2JJpItOklJa+Ltz27v18+cu3cuW2Xzx1g9t+g5dFLC6o5i5Mp0OpYnj8nUYFC6kAJJkJJD46VMk3BTEoo1TyH1Dy63/xuRL53vbdWnLbu/drxW2n6sCdPZrPbafqz1GfNDI1SpReMxr4e0H8N4vzB5D6+phMROk1s7h3SvZhwIZmh/0ghkYBAMCxOfoRIQAASAkX1AMAgKOTfaxTpP91hGvXrm3WrNn06dOFu4mJiXv27PH29n7++ef9/Pxk2jwAALAntlxHKHfH+cc5wq+++io+Pn7NmjXC3fXr1//9738PDQ3Nz88fNGiQTqeTbwsBAAAakJoxVlRUtGTJknnz5u3du5cxZrFY3nnnnf/85z+TJk165pln+vbt+/XXXz/xxBNyb2oDEVdslK41SqRJecFOsuYnVYOUTI2KrCkqcv3UdKzUCQBqebWGn5+jsoLcXUnteaq2J/WiuO9IPeun/lelUqPU1M1UDllJ/AK1lr8zuXPJMsYstZx2aiOdXfnnRNxbOnHbqXdc68TfSCrKS+WcFUrqE8j/vSaRU0mL+uRTUyKTOWTixZJvOUnuwyKJ2FB0W/Q5RakpGWOzZs167bXXfHx8hKbCwsLff/99+PDhwt3hw4cfPnxYtg0EAAA7ohD/Izfl9u3by8vLp02bZm3Kz8/XarXu7u7CXR8fn7y8POr5NTX8mv0AANDE1NYS18PaOfW8efOSkpLqNmm1WrPZXFtbK1zSazAYnJ2dyeerkTsFAHAIdzOGaUNYRvaDQnVOTk5UVBRjrKKiorKyskOHDqdOnaqtrc3LywsMDGSMZWdnBwQEkM9HRwgA4BjuqiNkCoXInk3ufpApz58/n5SUlJSUNGfOnJ49e+7du9fDwyMyMnLLli2MscrKyh07dkRHR8u9nQAAAA1C3b59e+GWt7e3s7OzcPftt99+8MEHk5OTL168GB4ePmzYMFk3siGRM9ET7UStUTORNOMm0MgaoRLNLG8g1i+21iiVSi3I5s8sf/4Yv/Kka3Mttz07g1+RsuwG53IdsZHXmmp+eUnqRVGfBKrcJfWvMZUaJbOFBA1R2ZLCjezeLNVzF756sZTbTr0ovY6/Mwuz+TPUUzvfTHx9qMMMNZFDNmilyVFzI8RUDpn6jlN/E6i/IUQR2UZwWCSVBig2Wltbe/z48dLS0sGDB1vzK1YXL15MT0/v1atX69atrY3Xr18vLi4ODQ0V7up0utTU1IqKip49ewqDnbm5udb4y/8GNmNiYiIjI4Xb/fv3v3jx4pEjR1q1anX//ffLnm0FAAA7If7yiXp7QpPJNH78+Nzc3DZt2jz11FP79u3r1q2b9dHFixevXbt24MCBBw8e/OCDD6ZOnXr06NGYmJjS0lIXF5eysjLG2Pnz5yMiIjp37uzj4/Pzzz+/+eabs2fPXr169SeffBIUFMTqdoQ+Pj7WKygYY97e3hMnThT1YgAAAKSVmJiYmZl55swZZ2fnV155ZfHixV999ZXwUG5u7jvvvHP+/PmQkJBdu3bNnDlz8uTJHTt2PHToUHFx8dixY4XF/Pz8zpw506ZNG8bYgQMHHnjggaeeeooxNn369JUrVzLMPgEAABISUqOifuofGv3+++8nT54sXLwwffr0bdu2Wa/i2LFjR1hYWEhICGNs1KhROp3uxIkTPj4+nTt3rjuRmY+Pj9ALMsY6d+5sMpmqq6sZYxUVFSkpKfn5+ch8AgCAdKQ+R5idnT1o0CDhdnBwsMFgKCwsFCpgZ2dnBwcHCw8plcqgoKDs7Oz6f9XKlSsfeOABYfjzwIEDaWlpaWlpjtIRUtWO6Il5iVCMyIl5+SXWeLP1Mrq2EzlTKFkLijqrTyxPvShi/YVEWIYbcmGMzEJQ072WF3OqNFB5DTIHQWy82AlyyU+O6AJaxHqI7aG2n3q93GpkZUX8d+TymSJue/61m/ytId7Bmmr+rMXUxlNF+EwmopQasR5y51BhGTLLxlmeKsJHfceprw9ZbI/6q098oOwuoWFLiTWmSE9Pf++99+o2Pvroo0J3pdfrNZo/ppgWbuj1f6TA9Hp93Uv4tFqt9SGuhISErVu3CrXSFixY8OabbzLGioqKHKUjBACARkun02VlZdVtMRr/+DfL39+/qOiPf91u3LihUCisEyL5+x3ymlgAACAASURBVPunpaVZn3Ljxg1/f3/qV2zZsuXVV1/dt2+fcBDZrFkzod3b2xsdIQAASMa2kdEePXrcckRoNXDgwKSkpPj4eMbYgQMHwsPDnZycrA8tXrxYr9c7OTllZGQUFBT06dOHu5Lvvvtu7ty5e/fu7dKly+2PoiMEAADpSH2OMC4ubvny5fHx8Z07d3711VdXrVrFGHvxxReLioq+/PLL3r17T506dcqUKatWrXr66adbtmxZUlKyfPny3NzcmpqaBQsWeHt7jxw5MjY29oEHHti4caOwzhdeeCE+Pr5Pnz4tW7Y8ePAgOkIAAGi8vLy8kpOTP/744xMnTnz22WejR49mjI0fP16YKDcxMfGDDz44ePDgk08++fTTTzPGVCqVp6enp6fn4sWLGWMeHh7u7u7C6UArtVo9adKkI0eOXL58uUuXLugIAQBAMgob5he80+Lt2rV755136raMHDlSuOHm5rZgwYK6D3l4eAjjqHXd3vLQQw899NBDwm2H7wipTCBREIuemFdE0oyq1US2EzOCUnE1quYTGYCk5rylpj8lgovl1PYT20Nl/Aw6Tjs50S6x5ykqomQa9VWk3nFqp1Fz4ZKzEIv8i0G9Xm4JOoWCn+qk3hGqJBs1Qa5KxU+BaojZg6mdQFc6FPfJJ78pYr5xZDqUaKc+CVQ7GVW3t3QoxYbZJ2R/6bigHgAAHJrDHxECAICEGqDodkNDRwgAANKx4YJ6uasGoCMEAADJ2OMM9ThHCAAADs3RjwjJeoBi06RkjI2XGhWZAjUTNUKp2qFU1k5su4nYTqNeROVGRr9eKjWqr+GUFaUm5qX2PKWWemupGqRUyJQITFIfHXI7xYVYmZlfcpUZGDHhsBjURlKpUY2Wnw6lorMaLX+niQ1e0lVzxeWWuZ/MWrJUr7ivj9gqtfXPyQcNytE7QgAAkJBC6ol57wEMjQIAgEPDESEAAEjGlgvq5R4VRkcIAACSsWU+Qrl7QgyNAgCAQ3OcI0Jx2UIq2SU6kMlLmlHZNjo7x21mtVQ7VWtU5MaTy0sUtaWKN3InYacmJacCimLTm04u/O+CT0AzbnvzFlpuewVRrrMot4rbrifSsGRNVOKTbDJy2qkMglrD35lqYgZ5Sy1/PWKDlNK1c5vJbwr5CeS1kyV8RX5NqL8htKZShNSGyjJyc5yOEAAAGh4uqAcAALAvOCIEAADJKJhC7HWBch8QoiMEAAAJYfaJxossdySyXWSJNW672GgAedadDK2IC7NQ2WWVhj9yrjbx2+lZi/m/V2wIiIs6G2EhvltKflEw5hvID8UMGtuG296+e0tu+5Xzxdz2Izuvcdvzrt3kbxDxwsjXy9v5dL6DvxLq75GSSO7QE/YSG0/8AmoKW2r7xX7yyXiXqK+nZCXW+O1NhkL8dYFy94M4RwgAAI7NYY4IAQCg4SkU4i+Ql/uCenSEAAAgHTs8R4ihUQAAcGg4IgQAAMnYUnS7Ybbk7qEjJIitIiYmxia2EJTI8nDkp5DK/rm6abjtVBWxmmp+VbCK0hpueyURt6Oqo3G3s5ao8kWh3hEq0Ojh5cxt7xzqw22/b6Aft11J7P205AJue8H1Sm47lUWkzr5wm5XEiI+CaFcRsw07u/L/ULi58z8hWqJeXYMnJsWGvflfT4lSoA5bYQ1FtwEAAOwLOkIAAHBoGBoFAAApoeg2AACAPcERIQAASMaWsEwDbcpda3IdoUSpNPGrId5KXrOFWD0ZP6NSoGTqkn+g7+LKb/cJdOO2twrit+uqjdz29LQSbvv138u47XodvwgpN9hJFiwVmeWjpl2lJtS9ksavHUq9iem/8nfCzTL++slZlLmt9RSM5X0YVMREu1Q7lQ71C27ObQ/u6MFt1zrz11Ocz5+duKyInzem/pySszETyNqq3J1JrVv2v9b2wg4vqG9yHSEAAMgHl08AAADYGRwRAgCAZGyoLCM7dIQAANCoXbp0adWqVeXl5dHR0Q899FDdh0pLS1euXHnlypW+ffvOmjVLq9VWV1cfPnz49OnTNTU1b7zxhrCY0Wj88MMPk5OT27dvP2/ePC8vL8bYjRs3Vq5cmZmZiaFRAABovAoLCwcPHuzr6/vggw/OmTPnyy+/rPvo+PHjr1y5Ehsb+8MPP8yaNYsxduTIkYULFx49evT999+3LjZv3ryvvvpq6tSp2dnZo0ePtlgsFotl5MiR+fn5U6dObXJHhNQhudhynSJjTFTxRu5JYNErJwYaxM4b7uzKrylKZf969ueX06yu5KdGDTX8GqRFefys4M1yA7ed+7qonUDVFKVyhbVEgLAwm1/z85cdmdz2s0fyuO3lJfwA5I1c/k6gU6D8duqjrOA9gfqEqInUaDOidmi7rp7c9r5RQdx2jRN//WnHC7nt+poibrtRz/9EUQVjyTSpmCCoqGqucDvJwzKffvpp//79X3/9dcaYyWRasWLF9OnThYeOHj166dKlgwcPajSaPn36dO7c+a233ho5cuTIkSOTk5N//vlnYbGysrL169enpqZ26tRpwoQJQUFBhw4dMhqNRUVFGzZsUCqporwAAAA2EU4Tivqpx/HjxyMiIoTbERERZ8+erampsT40cOBAjUbDGAsODvb3909NTb19DefPn/fw8OjUqRNjTKVSDR48ODk5OTk5efDgwUIniI4QAAAar7y8POGUHmPM29vbYrHk5+cLd/Pz81u2bGld0tvbOy+PM1STn59vXYN1sbqNfwyNGo1GtVp9y/GpwWDQavlDJQAAABw2XVB/8ODBsLCwum2bN28WjuFcXV31+j/qUQg3XF1dhbsuLi4Gw//OsNTU1FgfqsvFxcW6BmExT09Po9FYVvZHrQ/1999/P3v27OLiYsbY8OHD161b16pVqxs3bjzyyCMnTpzQarVvv/32zJkzRb4sAABwRDacI2SM9e7de8WKFXVbWrduLdwICgrKysoSbl+7ds3Z2dnb29v6UFJSknDbbDbn5OQEBXHOWwcFBeXl5RmNRmEQNSsrq2/fvkaj8ezZs8ICytDQ0GPHjlVXVxcWFjLG4uPjGWMvv/yyv79/aWnpgQMH5s+f//vvv4t9VQAAAHfJw8Mj7M9cXFyEhyZPnvztt99WVVUxxj7//POYmBilUpmSknLkyJEJEyakpqZeunSJMbZjx47mzZv37dv39pX36tXLz8/v22+/ZYxlZGQkJydPnDgxOjr6yJEjmZmZjDF1u3bthEXd3NxGjBixY8eOmpqar7766vjx4yqVqmfPnuPGjfv888/ffPPNe7I3GoxEte/ISdXFtFMrETVjO6NrilKTpItN1Eo1mzj13yEV1eLvTGIldGqU2mn89eiJyGtORgW3PY94s2rN/FgqMe05+eZSqNfLn6FeZLqSWt7Jmb/XqJQplRql1kOlQM0iq+mS3xQxX09JVsKY+JqlTSWVasMF9fUvP27cuISEhF69erVu3To9PX3fvn2Msc8//7yoqGjTpk1LliyJiIgICws7efLk+vXrVSpVdnb20KFD9Xr9zZs3Q0JC2rVr99NPP/373/9+/PHHP/vss9TU1Ndeey0wMJAx9vLLL/fv3z80NFTNGCstLU1KSrp27dpnn322atWq3Nzcmpqarl27ChvRtWtXob/lor6WAADgiGwruk33JCqV6ttvv01LSysrKwsPD3d2dmaMvfnmm2azmTE2f/782NjYK1eu9OzZUwi/+Pn5WcdLGWPCcOjYsWN/++23c+fOhYSEWAddX3/99bi4uKtXr6oZYzdv3ty3b19WVpZGo2nevHlZWZmTk5Na/UeOxt3dvaSEX1CfMVZdXS3yFQMAgF2qra2940V3Chsulb6L5Xv06FH3rofH/65+DgoKqntqUK1Wt2/f/vY1tGzZMjIy8pbGNm3atGnTRs0YCw4OXrNmDWNszZo1cXFxu3btqqmp0ev1Tk5OjLHS0tJWrVpRG9esWbM7vgAAAGgCmuql5396Vd27d8/Ly/P392/RosWZM2eExjNnzliHSQEAAOqhsOGCerm3Wb1p06b27dsHBwdfv3791VdfnTBhgkajefLJJxctWvTZZ5+dPHnywIEDq1evlns7/zqqbJKI6AGrby5cEdXOqBJoVLvYs/RUFIKaCDf7Sjm3nSpwVVPNXw+VK6GWFzXHLPVi+ZvIWC0ZKiF2psgwDjXBL0XkVLKiz76L2mkUPfFO5Wbe5LanHc/ntmu0/FBMfhZ/PQbikyZJNo3V8/XktVPJHbF13chZhfnN9T5iX6Q+R3gPqCsrK1966aW8vDxvb++xY8cKl0/885//nD9//qBBg3x8fLZu3RoQECDnNgIAADQY9cyZM2+/Xt7V1fWjjz6SZYMAAMCO2XBBvdwVzZvc7BMAACAfW64jbJgtuXtNMwIEAABwl3BECAAAkrFpPsIG2pa75fAdIZUhJJJgdKKMf2yt4s2AqtbwF6amS6VrPnGbyY3XE9m8guv8OWmpuWSpCXjLivlz0lLVy6gPP3enabT8V0ulK2uJVCcZxhRZl66hUX9KqE8gd/9QpcioiF4NkSvOzeTngXVV/CmatU781CiVDq0lCtA5OfP/QNFFB8VNRKzifRPJryG1EpEl2eT/qw+3wdAoAAA4NIc/IgQAAOnYEJaRPS2DjhAAACSjYOLPEcrdE6IjBAAA6dhWWUZWOEcIAAAOzdGPCMk5MqWoKcqIgKiWqMSoIbJ2VMpU7PSkFCoraDTwM35GPX/uWUMNUVOUyARSO5n7ehXEZ5XaCWYjfyPNJqKdmFCXqilKFkol9j2dKxYZgKQixxrOh4damBq5oqK21ZX8dCj1ztKfZH67kwu/XewcudTrpbaHWxOVGyVl9aRDxdYg5bY2IZJPzHsPOHpHCAAAUrKhxJrc/x5gaBQAABwajggBAEA6dhiWQUcIAACSscN+EEOjAADg2BzliJD6j8NCzURPVnqkMntEXI0XS9M48zNsWqKdmu9breGnOsXG28iZ7qndRu00qkArUaSRCl4qVZwHNMQe1hLBRRORGiWjsAZxy1OBSSrFKuoTwujaqtTO5O98ke8UGYDkr4Z8By38fUmGA0XNIM/oqC2106jap9xvHP11o9KkxDsicop62Y+KpKKgY8nkU+R+8Y7SEQIAwL1gh2Oj6AgBAEAyttQalRvOEQIAgEPDESEAAEjGpol5UXT73hB5glpBzXlLllIjzsbzztI7EaEYahpSJ2eiBJqevx6qXa3lZxjURK6Emi61lpwLl/97a4noBBmF4KUSnF3EzdFKVQWrLDdw2/VEnTkqIkHtHCqHQuU4nIjX5eauFbW8iVc6zqATNxEuWSmQyI9Q1cjoOnBEYoj4+lC5FS35TSHaiRJu3G8ilayhvuPU3wTqbwj5N9/ehhNJdniOEEOjAADg0BzmiBAAAO4Fhdj5BTEfIQAANB2KRnBdoFgYGgUAAIeGI0IAAJCOHYZl0BGKS3xRQUc1kQnkzghKZ9v4wUVnV/7yYquCGY1EhtDEj8PRBbSIuWqJdopKxX9dnr4utze2DvHgLuzcjL+S3Ks3ue2Zl0q57WXFNdx2lYq/ky3E3qGmY6WyiC15L5Yx1qaLJ7fd28+V215Rwtn+nKsV3IVvlvGjsxQq8kq2UyXNqAlyxZRAY3SEmPqmUFFb7jdR7PzYKqpCIdEu/1/9BmaPl09gaBQAAOyA0ci/Jop6qJ7lb4GOEAAAJCOUWBP3c6d1rlixwtPT08vLa+LEiRUVfxrnSEpKatOmjZeXV2ho6IULF25ZPjo6+ubNm4yxoUOHhtTx0ksv1dTUWO+iIwQAAOkobPqhnTx5cunSpSdPniwuLrZYLEuWLLE+pNPpYmNj//3vf5eXl0+dOjUuLo4xduLEiWXLlqWkpBQXF5vNZmH5xMTElJSUlJSUEydO6PX6/v37WyyWjIyMw4cPp6SkoCMEAADJKMSrvyNMSEh4+OGHO3TooNFo5s+fn5CQYH0oMTHRz88vOjpaoVC8+OKLaWlpFy5c+Oyzzx5++OGQkJC6y3t4eHh6enp6eqamptbU1ERHRwtraNGihaenJzpCAABovK5cudKtWzfhdrdu3YqKikpLS60Pde/eXbjt6urarl2733///Zblb9y4UVZWZl3bhg0bHn30UScnJ+Fux44dW7du7eipUbFhJaqSpJqYLpVfa5TIsLk003Db6XSoNHPP1pqJFKjY1CixPBWfo7J/QbyAaPiwIO7C7l5O3PbzyQXc9nJeupIxpqu62/PqAnJiXpEv1iewGbf9vgF+3PaQ7i257bmZnIComdhIs4mfJqVmMxY7q7CoiXDraae+KVQ61Jn4BlHfLO76tU7iCqWSsxlT9Y2beGjUxqsnbt68aT3DJwgJCRG6q9LSUjc3N6GxefPmQounpydjrKysrFmz/32DmjdvXlpaWlZWdvvyLVq0EJbftm3b0aNHGWMajSYpKSksLKy0tNTRO0IAAJCY+M7+5MmTDz30UN2W7777rkuXLowxHx+f8vJyoVE4tvP19RXuent7Z2RkWJ9SWlrq6+vr7e19y/I+Pj7C3Y0bN3bv3j00NJQxplarR4wYwRjz9PRERwgAADKLioratm0b96GuXbumpqYKt8+cORMUFGQ94OvateuaNWuE2+Xl5ZmZmV26dOnWrduZM2esy7du3dq6/IYNG5555pnbfwXOEQIAgHRsSMvUO1781FNPbdu2bd++fTk5OUuWLJk5cyZj7OOPP37nnXfGjBljNBrff//9oqKi1157LSIion379jNmzPjuu+/2799fd3nG2Llz5y5duhQbGyvcTU5O/vHHH3Nycs6dO4cjQgAAkIxwaaDYp9SjS5cu69atmz9/fmlp6aRJk+Lj4xljRqPRaDRqNJodO3bMmTPnX//6V3h4uBAQ7dat27p16+bNm1daWhoTEyMszxhLSUmZP3++cHJRWMPy5cszMjJatGiBjhAAABq1KVOmTJkypW7LrFmzhBu9evU6cODAHZdnjM2YMaPu3SFDhhw6dEi4jY6QIi4JRtUg5RZjpGaoN5I1RfkpUBNRO5ReXlzZTCoFSs1vTlUaFTtpO7edrOYqNrgostyl2EgtFZHViAxGUkFHVzd+u7Mrp93Jib9yas9TH29qhnp6BnmRKVCqnSgk6yKynaxByttO6kVR33F6WK+px0MpNsRG5YaOEAAAJKNQiC+iLfc1JQjLAACAQ8MRIQAASMaGGerlPiBERwgAANK5u/kkbnkKmTC4N9AR8pFvIzGWTJVZ4pZlIqcbNfLfDirkYjbxPzwmI7+dDLmIDMvQORF+u4nYOdRHvzhfd3vjryf4JdOoUEneNf7EvHodf/ZjqkiehchxUHuBOjtCRS2qb/JLu2VcKOG23yzll4grK+K0VxALU/ElLTFlNBWukao0GpUAauam5S/fnN8uqpQaIxJMkpVS4zc7ADucoR7nCAEAwKHhiBAAAKQj/oJ62aEjBAAAydyxZBr3KTKfI7RYLCdPnjx58qRCoYiMjLRO43Tu3LmkpCQfH58pU6a4uLjIupEAAAANRfnhhx8+9thjv/766/nz5wcMGPDZZ58xxrZv3z5s2LAbN258+eWXkZGRJhM/XwAAAPAnCpt+ZKV+5JFH/va3vymVSsZY3759ly5dGhcXt3jx4hUrVjz55JMmk6lnz57bt2+fNGmSzFt6j1FvjIVIjhGpIxUvi0gWpnLhDw+4EFW+qOpfZhMxAS+VGuUvThJbZ86oFlfyrbJCf3tjxq/8FCW1NTXERLsGPX9jqBJuKrU031HhK3a76koDtz3zUim3nUrDct90I/FiqXpyoifgpfLPZDqUn/Zs5k6kQ93F1ZkTmxrV8l4X9zvL6OJ5iI3ewpai2w2zJXdP6eXlZf2Kurq6KhSK4uLiU6dOPfjgg4wxtVo9ZsyYvXv3yrqRAAAADeV///tUV1cvWbJk9uzZubm5Go3Gy8tLaPf398/NzaWebzDw/58FAIAmhrrmuC7FH8eEIn5kPyb8oyM0GAxTpkwJDQ197rnnFApF3VdrsVhEV1AFAACHZW/nCJWMMaPR+Mgjjzg5OSUkJCiVSn9/f5PJVFxcLCyRn5/v7+9PPV+r5Q/uAwBAE9NUD4qUZrM5Li6uurp68+bNGo2GMebl5RUWFrZ9+3bGmMlk2rVr1wMPPCD3dgIAgB0QPTDaCC7AV69atWrz5s0PPvjgY489xhhTKBRfffXV66+//uSTT166dOnMmTMeHh5CcAZYfVXS+Q9wo4JUNo+asNdCpD2pFKiZSpkSy9PFQ4lmYi9QxRipGB41sbCJG4A08COm3IUZHUlVqoiymU7c5nqQwWJxayF2JlUTlYoEc3cyFYWlCquSqVEiZUqlMV2pWYWbEzVFiXRoM6qmKDk7sYgJeBnxeqlCrNQ7Jfsf8cbGlgvqG2hT7pr6gQce+Oqrr25pffDBBw8ePJiUlPT444/HxMSo1ShAAwAAd8EOi26ru3XrZq0mU1ePHj169Ohx7zcIAADgXsKhHgAASMeGc36yHxHK/PsBAKAJUYgPl4qdyFdymI8QAAAcGo4IpUH+A8RLoCmJzJ5GZKiTnFleZDqUrJVIpUOJ5cnUqIqfDjUa+curjZztNKmJtCQRPaUKn1KBSTJSSyFDo9R6xEVtqe1UaYjlNbzJ1qmVEAVUxdYUpWp7kjPOEzVFqdQoXVOUSIdSM9ETr4v7TaRqiiId2oShIwQAAMnYNh9hA23MXcLQKAAAODQcEQIAgGRsqRQj97AzOkIAAJBOIyiiLRY6wobF/TxQYRNGRRgs/BFsi4X/9omMa9STDpCmlJpOza8Wpqrhr4ebi1ETIRGTif9LydmJRc5yfDfzztwN0TuTaKeSVtxcDBWKoYv88T9RzkQ4hQzLUKXUiJJp1PJUWIYOxYgLB3E/+Pb2N7zRsekcYQNty93COUIAAHBoOCIEAADpiD9HiCNCAAAAOaEjBAAAh4ahUQAAkIyC2TAf4R2W1+v1u3fvLisrGzlyZEBAwC2P/vLLL7/99lvfvn3vu+8+6/K7du2qqKgYMWKEsPzly5evX79ufUpUVJRSqWSMHT58+PLly+gIGxiRx+S3EsfnVMZP9LaILaVGpEnJ6VuJQlYaJyO3vaaav7xRz6maRk20S6VDzdRsxiJTo2JLr1FfaCqaS6VGqZ0vqiSb2Il2XYiJbZ2pUmfE8lSalJpQl1qemmiXSodSr5famfyvhNznq+yd5NcR6vX6iIgIJyenDh06zJ07d8+ePeHh4dZHZ8+evXv37hEjRixcuPCNN9549tlnheWdnZ1DQkJeeumlvXv3hoWFrV69OjExMSQkRHhWZGSkUql86aWXduzYMXLkSHSEAADQeH3zzTc1NTVHjhxRq9VvvfXWkiVLEhMThYcyMzPXrVt35cqVgICARx555KGHHoqLi/v6668NBoOw/D//+c8lS5b88MMPjLHHHnts8eLF1tVmZWV9/PHHV65cCQwMxDlCAACQjsKmH9qPP/4YHR2tVqsZY5MnT961a5fZ/Meg0a5du+6//35h8HPw4MEKheLEiRM//vjjxIkTb18+Jyfnxx9//PXXX4Xn7t69u2/fvoGBgQxhGQAAkJDCJvWsMDs7W+iuGGOBgYEmk6mgoEC4m5OTExQUZP29AQEBOTk5OTk5dZc3Go2FhYUqlerSpUtr166NiooaP368wWCo+1wMjQIAgMwuXrz46quv1m2ZNWuWv78/Y6y2tlYItjDGVCoVY8xk+qNYldlsrtuJqlQqk8nEXX758uXC7YqKin79+q1Zs6buc3FECAAAkhHCMuJ+GFOpVG5/JvRbjDF/f//CwkLhdkFBgVKp9PPzu/0hxlh+fn5AQICfn9/ty1vX5u7uPnbs2DNnztR9Lo4IZUAPA1Az3vKbxaZJqfEHsWUw1URmT0tkEZ1c+O011fwapIYaTmrUaOBPwEulScWmQM1UrVFieaoGKbmTqRSoyHQoOeEw78MgdqJdKqVJtVM1P52JdidqPUSNU42TiMKqTGw6tBEUNGmyxKdGO3XqdMsRoVVkZOTXX38tPLpnz55BgwZptdqqqqra2trIyMjXXnutsrLSzc3twoUL5eXl4eHhkZGR33777SuvvCIsP3jwYI3mf7Fki8Vy+vTpYcOGRUZGxsfHV1RUuLu7oyMEAAAp3fG6QFEee+yxFStWxMXFde7c+d133920aRNj7JVXXikqKtq0adPw4cPHjRsXHR29bt26F1980d3d/fHHH1+5cuWMGTM6duz47rvvbt68mTE2bNiwAQMGuLm5HThw4OrVq1u3bvXx8Rk1atS4ceNiYmIwNAoAAI2Xu7v7iRMnevbsqdPp9uzZM3r0aMbYtGnTZs6cyRjbunVrXFzcjRs33nrrrTfffJMx5uHhcfLkye7du+t0ur17944aNYox9sorrzRr1kyn002bNu3ChQs+Pj6MsS1btjz11FM3btxQ/MVZZqbP/NvkeYskeK1AT59EvUfUqB01Wmg08Nv1Ov4QZQ3RTi2vJ4Y6qfVgaJRhaFRox9ConYjp7H/HZdbv/q2wrEbUai+k7L955eC2bdts3a6/CkOjAAAgGVvmI5S7nA+GRgEAwKHhiLARaeg0KT1qx18PNa+3RmStUS0177kLMTTKqzVKDo0S473UUKfYIVNyXJoYxxZb0JUaMiVnqCfaueU3yXeEyPdS7xSV+xW7HmqoU6Phr4ecWR5DoI3ZnSrF8J8iK3SEAAAgGcmLbt8DGBoFAACHhiNCAACQTEPMR9jQcEQIAAAODR0hAAA4NAyN2gGxaVIqBco04pa3VnC/hYq4nFlNpkb5gU+jgb+8Sc8JghqJC+fpC+qludD+rxWc+B8yTSo6NSpicnYq36um0p5EFVkqfUpVnaWWJ1OgxIuiwqFIhzZmNoRlZH/j0BECAIB0xF9QL/cpQnSEAAAgIRuuI5QbzhECAIBDwxEh1CC+jAAAESdJREFUAABIBucI4Z4SG6JRKogJE4iBAQWxPBV5UGv4+RSqEBeVczGZOO1mYmGz2FAM1U7O/UE0UyXW+M10mTyRs0/Qs1Jw3kQy1kSFaIh26h2n1k9uPPFiySQRtdfk/rsJ9bCl6LbcPSGGRgEAwKHhiBAAACRlb4fs6AgBAEAydhgaxdAoAAA4NhwRAgCAZOwxLIOOsCkSm8GjJuYVG3SkqoWp+QlLKqPIrXYmtjQatTw50a5cJdZETtgrKpApNr0penmRKVCURnMIdjgxL4ZGAQDAoeGIEAAApKMQPdSJoVEAAGg6bKks0zBbcveUN2/eXLp06ZQpU0aOHFlWVia0VldXP/vss23atOnbt++uXbvk3UQAALAnCvE/slJWVFRkZWUNHDjwp59+MhgMQuvChQuvXr167NixRYsWTZ06NTc3V96tBAAAaCDqwMDA1atXFxcXz507V2gyGo2ffvrp7t27AwICJkyYMGzYsISEhFdeeUXeDYWGI7ZmKdUsNmVqUYkIatJpT/4vJZcna4dKFA8ViUyNik6ZNuTKxaZAicVl/8cf7gEFU5AB9caKc44wNze3rKysd+/ewt3evXtfvHjx3m4VAADYJXucfYLzb2RxcbGzs7NWqxXutmjRorCwkHp+VVVVQ20aAAA0JrW1xAiMneN0hJ6enjU1NUajUbhbXl7u5eVFPd/V1bWhNg0AABoTpfIuLj23t6QM43aE/v7+Li4uly5dEu5eunQpJCSEer7s138AAEDjobCJvNusZIxdvXr12rVrjLFr165lZGQ4OztPnTp1+fLlZrP5/PnzO3fufOyxx+TdSgAAgAaiZoyNHj3aZDK1b98+NjZWoVBcuXJl+fLl06dPb9mypVarXbFiRceOHeXeTpCD6OnBxT2BXA83NUotLm5ieemKhzY0KQKZ4tObSIHCX2VDWEb2D5KaMXb58uVbWn18fPbu3Ws0GjUajRxbBQAA9qlxnPYTpb4Sa+gFAQBAJOmnYUpOTl65cmVFRcXEiROff/75usvn5OS88cYbV69eDQsLW7RokZubm7D8ihUrbt68GR0d/dxzzykUim3btiUmJmZlZQUHB8+ePTs0NDQzM3Pp0qXCSjD7BAAANF45OTmjRo2KioqKj49///3316xZY33IYrGMHTtWq9UuXrz44sWLzz77LGMsOzt71KhRI0aMiI+P/9e//rV27VrG2Jo1a/r167dw4cI2bdoMGTIkPT29sLBw27ZtYWFhYWFhir9YUGP6zL9NnrfoL75OaCLEnq4T8wT6nJ/IX4pzhOIfwDlCEMR09r/jMluPZxVX6kWt9uyRn/JO79m2bRv30cWLF58/f/6bb75hjH377beLFi26cOGC8ND+/funTZuWm5urVCrz8/PbtGmTmZn58ccfX7hwYevWrYyxb7755vXXX//111/rrnDAgAFPPPFEnz59pk+f/vvvvzPMPgFSashwjYLq8KT6pXZOZA8G0GCknpg3NTV1yJAhwu0BAwZcvHixurpauIQ9NTX1/vvvF65u9PPzCw4OPnfuXGpqamRkpHX5Cxcu6HQ6FxcXoUWv12dmZrZt25YxVlxc/Nhjj7m5uWFoFAAAGq/8/HxPT0/hdsuWLRljBQUFwt2CggLrQ8Kj+fn5dRtvWZ4xNnfu3K5du44aNcrT03POnDkTJkwICQnBESEAAEjGtgvk9+7d27p167ote/bs6datG2OsefPmOp1OaKyurhZahLtubm5ZWVnWp1RVVTVv3tzNzY1afsmSJT///PPBgwcVCkXHjh3/7//+T2hHRwgAAJKxrej2kCFD1q1bV7fRz89PuBEcHJyRkSHczsjIcHNzs1b9DA4O3r59u3DbaDRmZ2e3adOmTZs2tywvHBcuX7588+bNBw4cuL1oKIZGAQBAZi4uLq3/zHr9Xmxs7NatW0tKShhja9eunTp1qkKhOHDgwK5duyZOnHjx4sWUlBTG2JYtW1q1atWnTx9h+dLSUmF5oVDM+++//8knn+zbt8/av6anpwsltQ0GA44IAQBAMgqF6BrU9S8/YsSIsWPHdu3a1dfX12Kx7NmzhzH2/fffFxUVjRkz5j//+c+oUaPatWuXnZ29ZcsWhUIxcuTI0aNHd+3a1dvbW6FQ7N69mzEWHx+vVqt79OghrHP+/Pk6nW7VqlXBwcE5OTm4fAIAAO7K3Vw+8d2p68VVBlGrPfNLUs7J3dTlE4L8/PyysrJOnToJGVGz2WyxWNRqNWOsoqLi+vXrISEhzs7OdZcvLy/v2LFjPTNmFBcX5+fn+/r64ogQAAAaOz8/P+uoJmNMpVJZb7u7u3fv3r3+5bm8vLyE84XoCAEAQEp2d/0qOkIAAJCMDZdPyD4fITpCAACQjtSVZe4BXD4BAAAODUeEAAAgGdsuqJcXOkIAAJCMbSXW5IWhUQAAcGg4IgQAAOnYYVgGHSEAAEhGIf6cn+wjqRgaBQAAh4YjQgAAkIyCKRQixzrFLi85dIQAACAdG84Ryg1DowAA4NBwRAgAANIRf0G97EeQ6AgBAEAyKLoNAACODecIAQAA7AuOCAEAQDIKZsPQaANty91CRwgAAJKxwwprGBoFAADHhiNCAACQjh0eEqIjBAAAySgUNlwOIXNPiKFRAABwaDgiBAAA6YivLIPUKAAANCF2eI4QQ6MAAODQcEQIAACSseWCerkPCdERAgCAZBR2OPsEhkYBAMChoSMEAACHhqFRAACQjE3zETbQttwtdIQAACAZG84Ryt0PYmgUAAAcG44IAQBAOnZ4QT06QgAAkI74c4SynyTE0CgAADg0HBECAIBkFOIP8OQ+IPzLHaHZbJZkOwAAoJGrqqpq1qxZ/cvYcvmE3CcJFRaL5a88f/PmzR988IFa3WSPLHU6nVKpdHJykntD7gW9Xl9bW+vi4iL3htwLRqNRr9e7ubnJvSH3gtlsrqqqcnd3l3tD7gWLxVJeXt6iRQu5N+QeKS8vd3d3Fz8Xri3Wr1/foUMHyVdbVFSUn5/fo0cPydd8l/5qRwgAAGDXEJYBAACHho4QAAAcGjpCAABwaOgIAQDAoTXZtOdfkZ2dXVBQ0KVLF2tQuKSkZN++fW5ublFRUU0pQarX65OTk/Pz8zt06BAWFiY0mkym/fv3l5WVRUZG+vr6yruFEiopKUlJSSkpKfHz8xsyZIhKpRLaT5w4kZGRERoa2qVLF3m3sCGkp6frdDprHq+oqOjAgQPu7u5RUVEajUbebZNKSUnJ1atXrXc7d+4shIGvXLmSkpISHBw8cOBA+bauQWRmZiYnJ7u5uQ0YMMDLy4sxVlFR8dNPP2k0mhEjRjhI8FtCSI3+icFgCA4Orq6urqqqSk5O7tu3L2Ps8uXLERERQ4cOzc/Pr66u/vnnn11dXeXeUgnU1tZ6eXl17dq1Xbt2v/zyS79+/b7++muz2Txy5EidTtehQ4c9e/bs2bPH2kHauzFjxjDGfHx8zpw5o1arDx486O7u/vLLL2/dujUyMnLnzp1Lly6dMWOG3JsppatXr/bu3dvX1/e3335jjKWlpQ0bNmz48OHXr1+3WCwHDhxoGv/Vbdy48aWXXgoNDRXuvv/++927d9+yZcusWbPGjRt39OjRESNGfPTRR/JupIRWrly5dOnSyMhIvV7ftm3bVatWXb9+feDAgWFhYVVVVTk5OUeOHPH09JR7M+2KBeowm80XL140m83Nmzc/ceKE0Pj444/PmTNHeLR///5r1qyRdRslU1tbe/HiReF2cXGxh4fH/v37v/76665du+r1eovF8uabb44fP17WbWwQRqOxc+fOCQkJWVlZLi4u165ds1gs+/fvb9WqlfDCm4ba2toHHnhg9uzZHTt2FFoefvjhBQsWWCwWo9HYp0+fhIQEWTdQMl988cXYsWPrtphMpjZt2mzbts1isRQUFDRv3vzSpUsybZ3Ejh8/7u7unp6eXrfxxRdffOKJJywWS21t7ahRo5YuXSrPxtktnCP8E6VS2aVLF6XyT7slMTHx4YcfFh596KGHtm/fLtPWSUyhUFgHA1u2bOnm5qbT6bZv3x4dHa3VahljU6ZM2bVrl9FolHUzpWexWGpraz09PXfu3BkeHh4cHMwYi4yMtFgsx48fl3vrJLN27dq2bdtGRUUJdy0Wy/bt24VPslqtjomJSUxMlHUDpVRZWXngwIFz587V1tYyxs6dO1dUVDRu3DjGmK+v79ChQ5vM13bTpk2xsbFubm6//PJLYWGh0Gj9G6VQKKZMmdJkXuw9g47wDqqqqsrKyoKCgoS7QUFBOTk58m5SQ/jss8+0Wm1kZGROTo71xbZu3dpsNufn58u7bRLauHHjpEmTunbtOnXq1PHjx+fk5LRu3Vp4SKFQBAYGNpk3Nzc39/3331+2bJm1pbS0VKfTNdVPckFBwbJly8aPH9+3b9+8vLzc3Fw/Pz9rxaugoKDc3Fx5t1Aq6enply9fjoqKWrlyZdeuXTdt2mSxWPLy8up+bZvSO3tvICxzB8LxkDVYoVarDQaDrFskvZ9//vnll1/+4YcfXF1dTSaT9cUKN5rSEeH999/v6el58uTJNWvWxMbGGo1G64tljKnV6ibzYp9//vm33nqr7omiJvxJnjp16qOPPsoYMxqNkydPfvXVVydOnHjLO9tkXmx1dXV2dnZaWpqzs/OPP/44ffr0mJiYul/bpvRi7xkcEd5BixYtXF1db9y4IdwtKCgICAiQd5OklZycPGXKlM2bNw8YMIAx5ufnZ32xwsCLv7+/nNsnqY4dO44bN+6NN94YN27cRx995O/vbx1cYowVFhY2jRd79uzZAwcO7Nmz59lnn129enVhYeGzzz7r5OSkVqub5CfZGn/VaDSxsbGnT5/29/e/ceOG5b9JwMLCwibzYgMCAgYOHOjs7MwYGzFiRHl5+fXr11u1atUk39l7Bh3hnUVERCQlJQm3f/rpp4iICHm3R0KpqanR0dHr168fPny40HLLi+3Xr1+TjGKXl5e7ublFREQcPXq0qqqKMfbbb78VFRWFh4fLvWkSCAgIWLFiRVhYWFhYWEhIiJOTU1hYmJOT0+DBg5vqJ9nq3LlzgYGBPXv2ZIydOnWKMWYwGA4dOtRkXmxkZKT1WpGMjAyVSuXn53fL13bo0KHybaBdwuUTt3r33XeLi4vfe++92NhYf3//BQsWpKamTpo0aeHChbm5uV9++eW5c+datWol92ZKoLq6OigoqFOnTpGRkULLxIkTu3fv3rNnz9GjR3fr1u2f//znJ598Eh0dLetmSqOmpmb06NEjRozw8PA4duzY3r17jx8/HhISMn78eL1eHxMT8/HHH48aNWr58uVyb6nEfvjhh3/84x/C5RO7d++eNm3awoULr169+u23354/f164BM3e/e1vf3N3dw8MDDx79uyWLVv27NkzaNCgJUuWbN68edasWTt37qysrDx48KDcmykNnU7Xq1ev4cOH9+nT58MPP4yKivrXv/516tSpqKio+Pj4qqqq1atXp6SktG/fXu4ttSfoCG+1ZcuWiooK691p06a5ubmdPHly69atbm5ucXFxQsiwCdDr9QkJCXVbBg0a1L179/z8/E8//bS4uHjixIlDhgyRa/OkZbFYduzYkZKScvPmzQ4dOsTGxrZs2ZIxptfr169ff+XKlfDw8EceeeTezGVzL127di05OXnq1KnC3WPHjn333XctWrSIi4sLDAyUd9ukkpycnJSUVFpaGhAQ8PDDD1u/od98883Ro0fbtm379NNPN41rfwUlJSXr168vKioaMGCA9f/Uc+fObd68WaPRPP744w0xU1LTho4QAAAcGs4RAgCAQ0NHCAAADg0dIQAAOLT/BwxhIugXKrqYAAAAAElFTkSuQmCC",
"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
}