{
"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.102290e+01 7.041416e+00\n",
" * time: 0.0018489360809326172\n",
" 1 2.840777e+01 4.949200e+00\n",
" * time: 0.006920814514160156\n",
" 2 2.097527e+01 5.716790e+00\n",
" * time: 0.012439966201782227\n",
" 3 1.370089e+01 2.886502e+00\n",
" * time: 0.016038894653320312\n",
" 4 1.116671e+01 1.133748e+00\n",
" * time: 0.01960587501525879\n",
" 5 1.040005e+01 2.048981e+00\n",
" * time: 0.022694826126098633\n",
" 6 9.888087e+00 1.255660e+00\n",
" * time: 0.025656938552856445\n",
" 7 9.461884e+00 1.011384e+00\n",
" * time: 0.028729915618896484\n",
" 8 9.000984e+00 8.456923e-01\n",
" * time: 0.03168487548828125\n",
" 9 8.655796e+00 9.045388e-01\n",
" * time: 0.03452897071838379\n",
" 10 8.390322e+00 5.943918e-01\n",
" * time: 0.03753781318664551\n",
" 11 8.350848e+00 1.030523e+00\n",
" * time: 0.040007829666137695\n",
" 12 8.294046e+00 8.644571e-01\n",
" * time: 0.04257392883300781\n",
" 13 8.198117e+00 6.974578e-01\n",
" * time: 0.04592585563659668\n",
" 14 8.048967e+00 4.116745e-01\n",
" * time: 0.04960989952087402\n",
" 15 7.957076e+00 3.845398e-01\n",
" * time: 0.05331587791442871\n",
" 16 7.943642e+00 6.021342e-01\n",
" * time: 0.092681884765625\n",
" 17 7.911156e+00 5.087298e-01\n",
" * time: 0.09547591209411621\n",
" 18 7.879390e+00 5.412756e-01\n",
" * time: 0.09806990623474121\n",
" 19 7.842516e+00 3.481325e-01\n",
" * time: 0.10102200508117676\n",
" 20 7.812414e+00 2.676037e-01\n",
" * time: 0.10408997535705566\n",
" 21 7.782916e+00 3.394597e-01\n",
" * time: 0.10641288757324219\n",
" 22 7.765237e+00 3.985629e-01\n",
" * time: 0.10876202583312988\n",
" 23 7.746721e+00 1.605332e-01\n",
" * time: 0.11177182197570801\n",
" 24 7.732907e+00 2.088353e-01\n",
" * time: 0.1141350269317627\n",
" 25 7.728162e+00 2.831760e-01\n",
" * time: 0.11638784408569336\n",
" 26 7.725938e+00 2.903319e-01\n",
" * time: 0.11867284774780273\n",
" 27 7.715058e+00 2.002263e-01\n",
" * time: 0.12154793739318848\n",
" 28 7.705152e+00 2.600964e-01\n",
" * time: 0.12373900413513184\n",
" 29 7.699972e+00 2.184188e-01\n",
" * time: 0.12601494789123535\n",
" 30 7.692585e+00 2.279659e-01\n",
" * time: 0.12842392921447754\n",
" 31 7.689894e+00 2.321909e-01\n",
" * time: 0.1306748390197754\n",
" 32 7.683332e+00 2.787574e-01\n",
" * time: 0.13326787948608398\n",
" 33 7.680395e+00 1.890848e-01\n",
" * time: 0.13609600067138672\n",
" 34 7.672437e+00 1.547245e-01\n",
" * time: 0.13898992538452148\n",
" 35 7.669930e+00 2.164050e-01\n",
" * time: 0.14186501502990723\n",
" 36 7.665492e+00 1.855944e-01\n",
" * time: 0.14486002922058105\n",
" 37 7.661820e+00 1.272048e-01\n",
" * time: 0.14773893356323242\n",
" 38 7.657584e+00 1.362699e-01\n",
" * time: 0.15016889572143555\n",
" 39 7.654934e+00 8.766086e-02\n",
" * time: 0.15271401405334473\n",
" 40 7.654480e+00 1.598830e-01\n",
" * time: 0.15536999702453613\n",
" 41 7.652239e+00 9.934691e-02\n",
" * time: 0.1589670181274414\n",
" 42 7.650453e+00 7.527233e-02\n",
" * time: 0.16277003288269043\n",
" 43 7.649083e+00 1.150946e-01\n",
" * time: 0.165618896484375\n",
" 44 7.648177e+00 8.135321e-02\n",
" * time: 0.16861701011657715\n",
" 45 7.646946e+00 7.720579e-02\n",
" * time: 0.17166686058044434\n",
" 46 7.645058e+00 5.406924e-02\n",
" * time: 0.17465686798095703\n",
" 47 7.643652e+00 5.052017e-02\n",
" * time: 0.17757081985473633\n",
" 48 7.642202e+00 6.189246e-02\n",
" * time: 0.18045401573181152\n",
" 49 7.641855e+00 9.368634e-02\n",
" * time: 0.18334603309631348\n",
" 50 7.641054e+00 6.412553e-02\n",
" * time: 0.18594789505004883\n",
" 51 7.639683e+00 3.917063e-02\n",
" * time: 0.20866894721984863\n",
" 52 7.638715e+00 5.194167e-02\n",
" * time: 0.21099591255187988\n",
" 53 7.638035e+00 6.588134e-02\n",
" * time: 0.21354293823242188\n",
" 54 7.636989e+00 6.142308e-02\n",
" * time: 0.2164309024810791\n",
" 55 7.636072e+00 4.057438e-02\n",
" * time: 0.21939587593078613\n",
" 56 7.635475e+00 4.872530e-02\n",
" * time: 0.222398042678833\n",
" 57 7.634980e+00 3.900013e-02\n",
" * time: 0.22528791427612305\n",
" 58 7.634404e+00 4.278678e-02\n",
" * time: 0.2277059555053711\n",
" 59 7.634009e+00 5.571754e-02\n",
" * time: 0.2299790382385254\n",
" 60 7.633680e+00 5.079089e-02\n",
" * time: 0.23236584663391113\n",
" 61 7.633375e+00 5.332203e-02\n",
" * time: 0.23460984230041504\n",
" 62 7.632928e+00 7.453823e-02\n",
" * time: 0.2369098663330078\n",
" 63 7.632317e+00 5.953921e-02\n",
" * time: 0.2391979694366455\n",
" 64 7.631522e+00 4.226178e-02\n",
" * time: 0.24220800399780273\n",
" 65 7.631328e+00 7.342402e-02\n",
" * time: 0.24451398849487305\n",
" 66 7.630578e+00 4.370234e-02\n",
" * time: 0.2468578815460205\n",
" 67 7.629892e+00 6.013812e-02\n",
" * time: 0.24912190437316895\n",
" 68 7.629020e+00 3.830866e-02\n",
" * time: 0.2514188289642334\n",
" 69 7.628330e+00 4.472795e-02\n",
" * time: 0.25385403633117676\n",
" 70 7.627872e+00 5.994983e-02\n",
" * time: 0.25613903999328613\n",
" 71 7.627804e+00 7.089002e-02\n",
" * time: 0.25854992866516113\n",
" 72 7.627268e+00 6.465028e-02\n",
" * time: 0.261538028717041\n",
" 73 7.627131e+00 9.168751e-02\n",
" * time: 0.2640049457550049\n",
" 74 7.627032e+00 8.133178e-02\n",
" * time: 0.2662949562072754\n",
" 75 7.626524e+00 6.939501e-02\n",
" * time: 0.26857900619506836\n",
" 76 7.625646e+00 8.837095e-02\n",
" * time: 0.2708468437194824\n",
" 77 7.624706e+00 5.111281e-02\n",
" * time: 0.27394986152648926\n",
" 78 7.624605e+00 1.144389e-01\n",
" * time: 0.27618885040283203\n",
" 79 7.623676e+00 6.356848e-02\n",
" * time: 0.27858996391296387\n",
" 80 7.622880e+00 6.605025e-02\n",
" * time: 0.2809908390045166\n",
" 81 7.622372e+00 6.413931e-02\n",
" * time: 0.2839798927307129\n",
" 82 7.621806e+00 4.222750e-02\n",
" * time: 0.286984920501709\n",
" 83 7.621357e+00 5.332367e-02\n",
" * time: 0.2892758846282959\n",
" 84 7.620831e+00 4.491713e-02\n",
" * time: 0.30832791328430176\n",
" 85 7.620395e+00 3.744157e-02\n",
" * time: 0.311262845993042\n",
" 86 7.619932e+00 2.750860e-02\n",
" * time: 0.314162015914917\n",
" 87 7.619625e+00 4.952384e-02\n",
" * time: 0.31640100479125977\n",
" 88 7.619325e+00 4.109088e-02\n",
" * time: 0.3193199634552002\n",
" 89 7.618805e+00 4.077561e-02\n",
" * time: 0.3223879337310791\n",
" 90 7.618440e+00 6.256811e-02\n",
" * time: 0.32460999488830566\n",
" 91 7.618277e+00 5.624196e-02\n",
" * time: 0.3268589973449707\n",
" 92 7.617923e+00 4.926166e-02\n",
" * time: 0.3297419548034668\n",
" 93 7.617529e+00 3.595190e-02\n",
" * time: 0.33199286460876465\n",
" 94 7.617030e+00 4.915827e-02\n",
" * time: 0.33424997329711914\n",
" 95 7.616617e+00 4.009074e-02\n",
" * time: 0.3372659683227539\n",
" 96 7.616325e+00 3.765083e-02\n",
" * time: 0.3400909900665283\n",
" 97 7.615827e+00 4.831876e-02\n",
" * time: 0.3425748348236084\n",
" 98 7.615529e+00 4.149599e-02\n",
" * time: 0.3454148769378662\n",
" 99 7.615277e+00 5.030796e-02\n",
" * time: 0.3477199077606201\n",
" 100 7.614993e+00 2.840701e-02\n",
" * time: 0.35056400299072266\n",
" 101 7.614957e+00 6.053424e-02\n",
" * time: 0.3530080318450928\n",
" 102 7.614873e+00 4.829080e-02\n",
" * time: 0.3552408218383789\n",
" 103 7.614364e+00 4.870874e-02\n",
" * time: 0.3575890064239502\n",
" 104 7.614015e+00 6.855373e-02\n",
" * time: 0.3598289489746094\n",
" 105 7.613472e+00 5.783060e-02\n",
" * time: 0.36226892471313477\n",
" 106 7.613176e+00 5.116245e-02\n",
" * time: 0.36455202102661133\n",
" 107 7.612838e+00 3.983134e-02\n",
" * time: 0.36757683753967285\n",
" 108 7.612509e+00 3.729807e-02\n",
" * time: 0.36981201171875\n",
" 109 7.612240e+00 3.433138e-02\n",
" * time: 0.37215685844421387\n",
" 110 7.611957e+00 4.622372e-02\n",
" * time: 0.37439489364624023\n",
" 111 7.611711e+00 3.217555e-02\n",
" * time: 0.3774528503417969\n",
" 112 7.611399e+00 3.286435e-02\n",
" * time: 0.380511999130249\n",
" 113 7.611322e+00 4.109409e-02\n",
" * time: 0.3829958438873291\n",
" 114 7.611141e+00 3.494357e-02\n",
" * time: 0.3852348327636719\n",
" 115 7.610745e+00 5.965660e-02\n",
" * time: 0.38759398460388184\n",
" 116 7.610352e+00 4.204375e-02\n",
" * time: 0.3898189067840576\n",
" 117 7.610031e+00 3.620723e-02\n",
" * time: 0.3923499584197998\n",
" 118 7.609874e+00 4.097104e-02\n",
" * time: 0.4032559394836426\n",
" 119 7.609354e+00 4.943510e-02\n",
" * time: 0.4055008888244629\n",
" 120 7.609015e+00 2.944610e-02\n",
" * time: 0.40834784507751465\n",
" 121 7.608768e+00 4.623617e-02\n",
" * time: 0.41059398651123047\n",
" 122 7.608608e+00 3.723652e-02\n",
" * time: 0.41300296783447266\n",
" 123 7.608435e+00 2.723995e-02\n",
" * time: 0.4158759117126465\n",
" 124 7.608404e+00 5.236463e-02\n",
" * time: 0.41828489303588867\n",
" 125 7.608051e+00 3.931512e-02\n",
" * time: 0.421205997467041\n",
" 126 7.607862e+00 3.253523e-02\n",
" * time: 0.4241750240325928\n",
" 127 7.607687e+00 3.829422e-02\n",
" * time: 0.4265158176422119\n",
" 128 7.607525e+00 3.936380e-02\n",
" * time: 0.4288468360900879\n",
" 129 7.607300e+00 2.277224e-02\n",
" * time: 0.43172693252563477\n",
" 130 7.607037e+00 2.225021e-02\n",
" * time: 0.43400001525878906\n",
" 131 7.606898e+00 2.747788e-02\n",
" * time: 0.43624091148376465\n",
" 132 7.606719e+00 2.867698e-02\n",
" * time: 0.43854784965515137\n",
" 133 7.606591e+00 1.936495e-02\n",
" * time: 0.44152283668518066\n",
" 134 7.606474e+00 2.327721e-02\n",
" * time: 0.4444398880004883\n",
" 135 7.606350e+00 1.397533e-02\n",
" * time: 0.44738292694091797\n",
" 136 7.606283e+00 2.037536e-02\n",
" * time: 0.4496469497680664\n",
" 137 7.606139e+00 2.166394e-02\n",
" * time: 0.45197582244873047\n",
" 138 7.606028e+00 1.691680e-02\n",
" * time: 0.4548468589782715\n",
" 139 7.605931e+00 3.363934e-02\n",
" * time: 0.4571249485015869\n",
" 140 7.605885e+00 3.210235e-02\n",
" * time: 0.4593830108642578\n",
" 141 7.605735e+00 3.743706e-02\n",
" * time: 0.46168994903564453\n",
" 142 7.605672e+00 2.553263e-02\n",
" * time: 0.4640069007873535\n",
" 143 7.605569e+00 3.283540e-02\n",
" * time: 0.46627283096313477\n",
" 144 7.605387e+00 2.316794e-02\n",
" * time: 0.468552827835083\n",
" 145 7.605260e+00 3.224783e-02\n",
" * time: 0.4708118438720703\n",
" 146 7.605209e+00 3.186753e-02\n",
" * time: 0.47321295738220215\n",
" 147 7.605165e+00 3.512431e-02\n",
" * time: 0.4754478931427002\n",
" 148 7.605161e+00 5.130535e-02\n",
" * time: 0.4778409004211426\n",
" 149 7.605003e+00 2.895404e-02\n",
" * time: 0.4800739288330078\n",
" 150 7.604861e+00 4.363684e-02\n",
" * time: 0.482496976852417\n",
" 151 7.604709e+00 1.558589e-02\n",
" * time: 0.4855179786682129\n",
" 152 7.604545e+00 2.500268e-02\n",
" * time: 0.49652791023254395\n",
" 153 7.604422e+00 2.204940e-02\n",
" * time: 0.4995279312133789\n",
" 154 7.604363e+00 3.075254e-02\n",
" * time: 0.5018858909606934\n",
" 155 7.604280e+00 3.219137e-02\n",
" * time: 0.5042359828948975\n",
" 156 7.604230e+00 3.849645e-02\n",
" * time: 0.5066089630126953\n",
" 157 7.604108e+00 2.877284e-02\n",
" * time: 0.5096139907836914\n",
" 158 7.603934e+00 2.888318e-02\n",
" * time: 0.512617826461792\n",
" 159 7.603731e+00 1.754829e-02\n",
" * time: 0.5156219005584717\n",
" 160 7.603524e+00 2.132989e-02\n",
" * time: 0.5178928375244141\n",
" 161 7.603486e+00 3.557791e-02\n",
" * time: 0.520143985748291\n",
" 162 7.603342e+00 2.550083e-02\n",
" * time: 0.5224509239196777\n",
" 163 7.603185e+00 2.814757e-02\n",
" * time: 0.5247318744659424\n",
" 164 7.603074e+00 4.622987e-02\n",
" * time: 0.5271279811859131\n",
" 165 7.602846e+00 2.494884e-02\n",
" * time: 0.5300228595733643\n",
" 166 7.602685e+00 3.211304e-02\n",
" * time: 0.5323209762573242\n",
" 167 7.602513e+00 2.603963e-02\n",
" * time: 0.535161018371582\n",
" 168 7.602400e+00 4.820417e-02\n",
" * time: 0.5375568866729736\n",
" 169 7.602260e+00 4.101240e-02\n",
" * time: 0.5397899150848389\n",
" 170 7.602105e+00 3.786428e-02\n",
" * time: 0.5428140163421631\n",
" 171 7.602016e+00 4.638605e-02\n",
" * time: 0.5450470447540283\n",
" 172 7.601765e+00 4.340607e-02\n",
" * time: 0.5480949878692627\n",
" 173 7.601575e+00 4.458416e-02\n",
" * time: 0.5511059761047363\n",
" 174 7.601400e+00 5.454856e-02\n",
" * time: 0.5535039901733398\n",
" 175 7.601009e+00 3.113959e-02\n",
" * time: 0.5557270050048828\n",
" 176 7.600793e+00 2.206292e-02\n",
" * time: 0.5587968826293945\n",
" 177 7.600657e+00 2.376004e-02\n",
" * time: 0.5618908405303955\n",
" 178 7.600463e+00 2.238766e-02\n",
" * time: 0.5642938613891602\n",
" 179 7.600371e+00 3.672642e-02\n",
" * time: 0.5665860176086426\n",
" 180 7.600316e+00 4.326935e-02\n",
" * time: 0.5689868927001953\n",
" 181 7.600138e+00 5.547518e-02\n",
" * time: 0.5715160369873047\n",
" 182 7.599836e+00 4.137294e-02\n",
" * time: 0.5740549564361572\n",
" 183 7.599775e+00 4.531431e-02\n",
" * time: 0.5765838623046875\n",
" 184 7.599562e+00 5.338432e-02\n",
" * time: 0.5789139270782471\n",
" 185 7.599180e+00 4.253677e-02\n",
" * time: 0.5901570320129395\n",
" 186 7.599165e+00 4.784022e-02\n",
" * time: 0.5927720069885254\n",
" 187 7.599000e+00 4.425419e-02\n",
" * time: 0.5951118469238281\n",
" 188 7.598760e+00 3.243815e-02\n",
" * time: 0.5982010364532471\n",
" 189 7.598566e+00 4.819262e-02\n",
" * time: 0.6005439758300781\n",
" 190 7.598439e+00 5.266591e-02\n",
" * time: 0.6029269695281982\n",
" 191 7.598423e+00 6.374650e-02\n",
" * time: 0.6052608489990234\n",
" 192 7.598010e+00 6.289213e-02\n",
" * time: 0.6083559989929199\n",
" 193 7.597793e+00 5.683004e-02\n",
" * time: 0.6106388568878174\n",
" 194 7.597546e+00 6.211485e-02\n",
" * time: 0.613090991973877\n",
" 195 7.597473e+00 5.042880e-02\n",
" * time: 0.6154630184173584\n",
" 196 7.597134e+00 4.642671e-02\n",
" * time: 0.6178610324859619\n",
" 197 7.597072e+00 6.263889e-02\n",
" * time: 0.620189905166626\n",
" 198 7.596587e+00 5.182877e-02\n",
" * time: 0.6230578422546387\n",
" 199 7.596243e+00 4.544697e-02\n",
" * time: 0.6260459423065186\n",
" 200 7.596068e+00 3.643112e-02\n",
" * time: 0.6284899711608887\n",
" 201 7.596041e+00 6.219180e-02\n",
" * time: 0.6307530403137207\n",
" 202 7.595757e+00 5.405962e-02\n",
" * time: 0.633836030960083\n",
" 203 7.595583e+00 3.919552e-02\n",
" * time: 0.6368780136108398\n",
" 204 7.595555e+00 5.148544e-02\n",
" * time: 0.6393508911132812\n",
" 205 7.595329e+00 4.462371e-02\n",
" * time: 0.6418328285217285\n",
" 206 7.595030e+00 4.471488e-02\n",
" * time: 0.6442320346832275\n",
" 207 7.594872e+00 3.658050e-02\n",
" * time: 0.6473250389099121\n",
" 208 7.594684e+00 3.072993e-02\n",
" * time: 0.6495640277862549\n",
" 209 7.594620e+00 3.742502e-02\n",
" * time: 0.6520109176635742\n",
" 210 7.594491e+00 3.668803e-02\n",
" * time: 0.6549868583679199\n",
" 211 7.594338e+00 2.719559e-02\n",
" * time: 0.6582369804382324\n",
" 212 7.594192e+00 1.184969e-02\n",
" * time: 0.6614518165588379\n",
" 213 7.594130e+00 1.735984e-02\n",
" * time: 0.6637728214263916\n",
" 214 7.594117e+00 3.105997e-02\n",
" * time: 0.6662440299987793\n",
" 215 7.594115e+00 2.907798e-02\n",
" * time: 0.6686220169067383\n",
" 216 7.594099e+00 2.492571e-02\n",
" * time: 0.6709568500518799\n",
" 217 7.594037e+00 1.858997e-02\n",
" * time: 0.6733658313751221\n",
" 218 7.593967e+00 1.892479e-02\n",
" * time: 0.675743818283081\n",
" 219 7.593841e+00 1.893340e-02\n",
" * time: 0.6869869232177734\n",
" 220 7.593720e+00 7.826913e-03\n",
" * time: 0.6900689601898193\n",
" 221 7.593704e+00 1.684473e-02\n",
" * time: 0.6925289630889893\n",
" 222 7.593659e+00 1.214727e-02\n",
" * time: 0.694943904876709\n",
" 223 7.593654e+00 1.549655e-02\n",
" * time: 0.6974649429321289\n",
" 224 7.593621e+00 1.408484e-02\n",
" * time: 0.6998920440673828\n",
" 225 7.593587e+00 8.559957e-03\n",
" * time: 0.7030520439147949\n",
" 226 7.593542e+00 1.016134e-02\n",
" * time: 0.7062668800354004\n",
" 227 7.593495e+00 1.094938e-02\n",
" * time: 0.7093589305877686\n",
" 228 7.593444e+00 6.977957e-03\n",
" * time: 0.7124788761138916\n",
" 229 7.593424e+00 1.182483e-02\n",
" * time: 0.71486496925354\n",
" 230 7.593403e+00 1.310013e-02\n",
" * time: 0.7173929214477539\n",
" 231 7.593392e+00 1.383707e-02\n",
" * time: 0.7197329998016357\n",
" 232 7.593380e+00 1.695454e-02\n",
" * time: 0.7223069667816162\n",
" 233 7.593365e+00 1.192060e-02\n",
" * time: 0.7246530055999756\n",
" 234 7.593353e+00 1.355933e-02\n",
" * time: 0.7271528244018555\n",
" 235 7.593332e+00 9.094047e-03\n",
" * time: 0.7295999526977539\n",
" 236 7.593288e+00 8.053882e-03\n",
" * time: 0.7320489883422852\n",
" 237 7.593265e+00 1.405932e-02\n",
" * time: 0.73447585105896\n",
" 238 7.593234e+00 1.048732e-02\n",
" * time: 0.7369370460510254\n",
" 239 7.593213e+00 7.938562e-03\n",
" * time: 0.7400178909301758\n",
" 240 7.593193e+00 7.865888e-03\n",
" * time: 0.7432348728179932\n",
" 241 7.593187e+00 1.138788e-02\n",
" * time: 0.7456369400024414\n",
" 242 7.593174e+00 9.469342e-03\n",
" * time: 0.7486598491668701\n",
" 243 7.593160e+00 6.798705e-03\n",
" * time: 0.7516608238220215\n",
" 244 7.593159e+00 1.093102e-02\n",
" * time: 0.7542450428009033\n",
" 245 7.593146e+00 7.775876e-03\n",
" * time: 0.7575018405914307\n",
" 246 7.593136e+00 8.500382e-03\n",
" * time: 0.7603809833526611\n",
" 247 7.593127e+00 6.138444e-03\n",
" * time: 0.7627630233764648\n",
" 248 7.593117e+00 8.954750e-03\n",
" * time: 0.7649998664855957\n",
" 249 7.593111e+00 6.341075e-03\n",
" * time: 0.7674088478088379\n",
" 250 7.593105e+00 6.304217e-03\n",
" * time: 0.7697160243988037\n",
" 251 7.593095e+00 7.786700e-03\n",
" * time: 0.7720248699188232\n",
" 252 7.593082e+00 3.817842e-03\n",
" * time: 0.7751288414001465\n",
" 253 7.593072e+00 7.646007e-03\n",
" * time: 0.7859518527984619\n",
" 254 7.593064e+00 8.282611e-03\n",
" * time: 0.7883968353271484\n",
" 255 7.593050e+00 5.374068e-03\n",
" * time: 0.7912449836730957\n",
" 256 7.593043e+00 5.705911e-03\n",
" * time: 0.7935609817504883\n",
" 257 7.593038e+00 5.497065e-03\n",
" * time: 0.7958269119262695\n",
" 258 7.593029e+00 4.878459e-03\n",
" * time: 0.7981808185577393\n",
" 259 7.593017e+00 4.652318e-03\n",
" * time: 0.8010509014129639\n",
" 260 7.593015e+00 1.104281e-02\n",
" * time: 0.8035018444061279\n",
" 261 7.593013e+00 8.086203e-03\n",
" * time: 0.8057370185852051\n",
" 262 7.593011e+00 7.618414e-03\n",
" * time: 0.808107852935791\n",
" 263 7.592999e+00 4.706563e-03\n",
" * time: 0.8104839324951172\n",
" 264 7.592989e+00 9.396902e-03\n",
" * time: 0.8128619194030762\n",
" 265 7.592985e+00 1.042203e-02\n",
" * time: 0.8151969909667969\n",
" 266 7.592970e+00 8.589162e-03\n",
" * time: 0.8176829814910889\n",
" 267 7.592963e+00 1.237687e-02\n",
" * time: 0.8200428485870361\n",
" 268 7.592962e+00 1.540658e-02\n",
" * time: 0.8223910331726074\n",
" 269 7.592951e+00 8.275880e-03\n",
" * time: 0.8246569633483887\n",
" 270 7.592946e+00 1.262295e-02\n",
" * time: 0.8271238803863525\n",
" 271 7.592933e+00 1.933851e-02\n",
" * time: 0.8294618129730225\n",
" 272 7.592911e+00 1.027468e-02\n",
" * time: 0.8324909210205078\n",
" 273 7.592893e+00 1.180618e-02\n",
" * time: 0.8354220390319824\n",
" 274 7.592883e+00 1.038783e-02\n",
" * time: 0.8378329277038574\n",
" 275 7.592880e+00 1.482674e-02\n",
" * time: 0.8400628566741943\n",
" 276 7.592855e+00 9.504104e-03\n",
" * time: 0.8431220054626465\n",
" 277 7.592834e+00 9.866064e-03\n",
" * time: 0.8453998565673828\n",
" 278 7.592828e+00 1.989648e-02\n",
" * time: 0.8478028774261475\n",
" 279 7.592787e+00 1.205759e-02\n",
" * time: 0.8507568836212158\n",
" 280 7.592785e+00 1.579968e-02\n",
" * time: 0.8531529903411865\n",
" 281 7.592764e+00 1.371804e-02\n",
" * time: 0.8560290336608887\n",
" 282 7.592729e+00 1.070743e-02\n",
" * time: 0.8590638637542725\n",
" 283 7.592704e+00 8.149899e-03\n",
" * time: 0.8612768650054932\n",
" 284 7.592688e+00 1.431637e-02\n",
" * time: 0.8636538982391357\n",
" 285 7.592658e+00 9.374058e-03\n",
" * time: 0.8666868209838867\n",
" 286 7.592642e+00 8.178616e-03\n",
" * time: 0.8699498176574707\n",
" 287 7.592636e+00 1.358512e-02\n",
" * time: 0.8808798789978027\n",
" 288 7.592625e+00 1.516025e-02\n",
" * time: 0.8833749294281006\n",
" 289 7.592611e+00 9.263103e-03\n",
" * time: 0.8863129615783691\n",
" 290 7.592602e+00 9.580045e-03\n",
" * time: 0.8887369632720947\n",
" 291 7.592579e+00 6.095229e-03\n",
" * time: 0.8910820484161377\n",
" 292 7.592565e+00 1.353807e-02\n",
" * time: 0.8934588432312012\n",
" 293 7.592553e+00 6.658948e-03\n",
" * time: 0.896759033203125\n",
" 294 7.592552e+00 7.859525e-03\n",
" * time: 0.8992648124694824\n",
" 295 7.592549e+00 5.818541e-03\n",
" * time: 0.9017248153686523\n",
" 296 7.592542e+00 5.185018e-03\n",
" * time: 0.9041600227355957\n",
" 297 7.592530e+00 7.175007e-03\n",
" * time: 0.906466007232666\n",
" 298 7.592529e+00 1.066231e-02\n",
" * time: 0.9087958335876465\n",
" 299 7.592523e+00 9.963244e-03\n",
" * time: 0.9110579490661621\n",
" 300 7.592514e+00 1.142981e-02\n",
" * time: 0.9133880138397217\n",
" 301 7.592513e+00 1.093670e-02\n",
" * time: 0.9157090187072754\n",
" 302 7.592503e+00 8.528525e-03\n",
" * time: 0.9180879592895508\n",
" 303 7.592503e+00 1.203676e-02\n",
" * time: 0.9204168319702148\n",
" 304 7.592493e+00 1.025464e-02\n",
" * time: 0.9227769374847412\n",
" 305 7.592474e+00 7.746780e-03\n",
" * time: 0.9259998798370361\n",
" 306 7.592463e+00 6.256538e-03\n",
" * time: 0.9289908409118652\n",
" 307 7.592460e+00 6.178687e-03\n",
" * time: 0.931380033493042\n",
" 308 7.592452e+00 7.753947e-03\n",
" * time: 0.9336960315704346\n",
" 309 7.592441e+00 4.178463e-03\n",
" * time: 0.936621904373169\n",
" 310 7.592434e+00 3.672045e-03\n",
" * time: 0.9396660327911377\n",
" 311 7.592430e+00 8.876585e-03\n",
" * time: 0.9418928623199463\n",
" 312 7.592419e+00 5.661403e-03\n",
" * time: 0.9449348449707031\n",
" 313 7.592418e+00 1.081678e-02\n",
" * time: 0.9472808837890625\n",
" 314 7.592414e+00 7.256687e-03\n",
" * time: 0.9497430324554443\n",
" 315 7.592405e+00 5.136565e-03\n",
" * time: 0.9519619941711426\n",
" 316 7.592398e+00 5.228207e-03\n",
" * time: 0.9543240070343018\n",
" 317 7.592386e+00 4.698692e-03\n",
" * time: 0.9565229415893555\n",
" 318 7.592385e+00 9.713478e-03\n",
" * time: 0.9589769840240479\n",
" 319 7.592378e+00 5.990810e-03\n",
" * time: 0.9618508815765381\n",
" 320 7.592371e+00 6.680337e-03\n",
" * time: 0.9642598628997803\n",
" 321 7.592366e+00 5.503588e-03\n",
" * time: 0.9760148525238037\n",
" 322 7.592362e+00 6.541491e-03\n",
" * time: 0.9784619808197021\n",
" 323 7.592358e+00 5.164568e-03\n",
" * time: 0.9808449745178223\n",
" 324 7.592357e+00 1.151333e-02\n",
" * time: 0.9831838607788086\n",
" 325 7.592351e+00 5.503384e-03\n",
" * time: 0.9860849380493164\n",
" 326 7.592347e+00 6.053568e-03\n",
" * time: 0.9884779453277588\n",
" 327 7.592346e+00 5.830202e-03\n",
" * time: 0.9908618927001953\n",
" 328 7.592343e+00 4.895123e-03\n",
" * time: 0.9931788444519043\n",
" 329 7.592343e+00 7.158989e-03\n",
" * time: 0.9955618381500244\n",
" 330 7.592340e+00 3.902016e-03\n",
" * time: 0.9978978633880615\n",
" 331 7.592338e+00 3.387040e-03\n",
" * time: 1.0002338886260986\n",
" 332 7.592336e+00 4.372013e-03\n",
" * time: 1.0024988651275635\n",
" 333 7.592333e+00 2.911150e-03\n",
" * time: 1.0056488513946533\n",
" 334 7.592331e+00 3.097144e-03\n",
" * time: 1.0086169242858887\n",
" 335 7.592329e+00 1.127305e-03\n",
" * time: 1.011512041091919\n",
" 336 7.592329e+00 3.052394e-03\n",
" * time: 1.0137269496917725\n",
" 337 7.592328e+00 2.238078e-03\n",
" * time: 1.015941858291626\n",
" 338 7.592328e+00 2.812973e-03\n",
" * time: 1.018263816833496\n",
" 339 7.592328e+00 3.167805e-03\n",
" * time: 1.0206060409545898\n",
" 340 7.592327e+00 2.458869e-03\n",
" * time: 1.0228359699249268\n",
" 341 7.592327e+00 2.349392e-03\n",
" * time: 1.0251789093017578\n",
" 342 7.592326e+00 1.384441e-03\n",
" * time: 1.0275919437408447\n",
" 343 7.592325e+00 7.736827e-04\n",
" * time: 1.0308029651641846\n",
" 344 7.592324e+00 1.574218e-03\n",
" * time: 1.0332670211791992\n",
" 345 7.592324e+00 2.662153e-03\n",
" * time: 1.035696029663086\n",
" 346 7.592324e+00 2.477905e-03\n",
" * time: 1.038064956665039\n",
" 347 7.592323e+00 1.475692e-03\n",
" * time: 1.0412909984588623\n",
" 348 7.592323e+00 2.079524e-03\n",
" * time: 1.0436749458312988\n",
" 349 7.592322e+00 1.936475e-03\n",
" * time: 1.0466468334197998\n",
" 350 7.592322e+00 1.695085e-03\n",
" * time: 1.0498619079589844\n",
" 351 7.592321e+00 1.225224e-03\n",
" * time: 1.0528738498687744\n",
" 352 7.592321e+00 2.245174e-03\n",
" * time: 1.0553250312805176\n",
" 353 7.592321e+00 1.579957e-03\n",
" * time: 1.05769681930542\n",
" 354 7.592320e+00 1.190595e-03\n",
" * time: 1.0601890087127686\n",
" 355 7.592320e+00 2.164596e-03\n",
" * time: 1.0628750324249268\n",
" 356 7.592319e+00 8.348839e-04\n",
" * time: 1.0747179985046387\n",
" 357 7.592319e+00 1.344716e-03\n",
" * time: 1.076982021331787\n",
" 358 7.592319e+00 1.485746e-03\n",
" * time: 1.0794780254364014\n",
" 359 7.592319e+00 8.743250e-04\n",
" * time: 1.0824038982391357\n",
" 360 7.592318e+00 8.728192e-04\n",
" * time: 1.084758996963501\n",
" 361 7.592318e+00 9.937941e-04\n",
" * time: 1.0870108604431152\n",
" 362 7.592318e+00 1.526824e-03\n",
" * time: 1.0894699096679688\n",
" 363 7.592318e+00 1.139431e-03\n",
" * time: 1.0923569202423096\n",
" 364 7.592317e+00 9.830437e-04\n",
" * time: 1.0948569774627686\n",
" 365 7.592317e+00 1.132396e-03\n",
" * time: 1.0971310138702393\n",
" 366 7.592317e+00 1.005122e-03\n",
" * time: 1.0995948314666748\n",
" 367 7.592317e+00 8.292739e-04\n",
" * time: 1.1018319129943848\n",
" 368 7.592317e+00 4.801830e-04\n",
" * time: 1.1048710346221924\n",
" 369 7.592317e+00 6.446411e-04\n",
" * time: 1.1071078777313232\n",
" 370 7.592317e+00 5.355400e-04\n",
" * time: 1.1095709800720215\n",
" 371 7.592316e+00 6.340734e-04\n",
" * time: 1.1118450164794922\n",
" 372 7.592316e+00 9.140430e-04\n",
" * time: 1.1142048835754395\n",
" 373 7.592316e+00 5.866546e-04\n",
" * time: 1.1170649528503418\n",
" 374 7.592316e+00 8.865160e-04\n",
" * time: 1.1195900440216064\n",
" 375 7.592316e+00 5.263003e-04\n",
" * time: 1.1225619316101074\n",
" 376 7.592316e+00 6.642944e-04\n",
" * time: 1.1250660419464111\n",
" 377 7.592316e+00 9.916131e-04\n",
" * time: 1.1274638175964355\n",
" 378 7.592316e+00 6.791649e-04\n",
" * time: 1.1299118995666504\n",
" 379 7.592316e+00 6.415214e-04\n",
" * time: 1.1323189735412598\n",
" 380 7.592316e+00 6.370638e-04\n",
" * time: 1.1353459358215332\n",
" 381 7.592316e+00 7.560462e-04\n",
" * time: 1.1375980377197266\n",
" 382 7.592316e+00 6.905203e-04\n",
" * time: 1.1400208473205566\n",
" 383 7.592315e+00 9.276054e-04\n",
" * time: 1.142401933670044\n",
" 384 7.592315e+00 5.670935e-04\n",
" * time: 1.145495891571045\n",
" 385 7.592315e+00 4.707314e-04\n",
" * time: 1.148535966873169\n",
" 386 7.592315e+00 7.029753e-04\n",
" * time: 1.1510210037231445\n",
" 387 7.592315e+00 1.050090e-03\n",
" * time: 1.1535580158233643\n",
" 388 7.592315e+00 1.014299e-03\n",
" * time: 1.1559009552001953\n",
" 389 7.592315e+00 4.029314e-04\n",
" * time: 1.1592259407043457\n",
" 390 7.592315e+00 7.978936e-04\n",
" * time: 1.1702518463134766\n",
" 391 7.592315e+00 6.604000e-04\n",
" * time: 1.1725749969482422\n",
" 392 7.592315e+00 5.234645e-04\n",
" * time: 1.1756179332733154\n",
" 393 7.592315e+00 6.564376e-04\n",
" * time: 1.17787504196167\n",
" 394 7.592315e+00 3.842406e-04\n",
" * time: 1.1808419227600098\n",
" 395 7.592315e+00 5.187252e-04\n",
" * time: 1.1831228733062744\n",
" 396 7.592315e+00 4.479018e-04\n",
" * time: 1.1856958866119385\n",
" 397 7.592315e+00 6.042387e-04\n",
" * time: 1.187978982925415\n",
" 398 7.592315e+00 3.461923e-04\n",
" * time: 1.1910479068756104\n",
" 399 7.592315e+00 5.946887e-04\n",
" * time: 1.1933579444885254\n",
" 400 7.592314e+00 6.422876e-04\n",
" * time: 1.1957910060882568\n",
" 401 7.592314e+00 3.233411e-04\n",
" * time: 1.1988060474395752\n",
" 402 7.592314e+00 6.672690e-04\n",
" * time: 1.2010698318481445\n",
" 403 7.592314e+00 5.501135e-04\n",
" * time: 1.2034008502960205\n",
" 404 7.592314e+00 5.179317e-04\n",
" * time: 1.2057700157165527\n",
" 405 7.592314e+00 4.986598e-04\n",
" * time: 1.20806884765625\n",
" 406 7.592314e+00 6.392393e-04\n",
" * time: 1.210526943206787\n",
" 407 7.592314e+00 5.532676e-04\n",
" * time: 1.2128360271453857\n",
" 408 7.592314e+00 5.947960e-04\n",
" * time: 1.2152738571166992\n",
" 409 7.592314e+00 6.107469e-04\n",
" * time: 1.217512845993042\n",
" 410 7.592314e+00 8.784793e-04\n",
" * time: 1.2199790477752686\n",
" 411 7.592314e+00 5.457388e-04\n",
" * time: 1.2224488258361816\n",
" 412 7.592314e+00 5.207802e-04\n",
" * time: 1.2247488498687744\n",
" 413 7.592314e+00 6.685354e-04\n",
" * time: 1.227741003036499\n",
" 414 7.592314e+00 5.542814e-04\n",
" * time: 1.230841875076294\n",
" 415 7.592314e+00 4.433835e-04\n",
" * time: 1.2340190410614014\n",
" 416 7.592314e+00 2.840974e-04\n",
" * time: 1.2370150089263916\n",
" 417 7.592314e+00 2.116143e-04\n",
" * time: 1.2401340007781982\n",
" 418 7.592314e+00 6.686302e-04\n",
" * time: 1.2423639297485352\n",
" 419 7.592314e+00 4.637688e-04\n",
" * time: 1.2446789741516113\n",
" 420 7.592314e+00 7.000199e-04\n",
" * time: 1.2469170093536377\n",
" 421 7.592313e+00 6.319254e-04\n",
" * time: 1.2493369579315186\n",
" 422 7.592313e+00 5.643321e-04\n",
" * time: 1.2516438961029053\n",
" 423 7.592313e+00 5.706201e-04\n",
" * time: 1.2542028427124023\n",
" 424 7.592313e+00 6.002635e-04\n",
" * time: 1.265002965927124\n",
" 425 7.592313e+00 5.652962e-04\n",
" * time: 1.2672548294067383\n",
" 426 7.592313e+00 6.694805e-04\n",
" * time: 1.2696988582611084\n",
" 427 7.592313e+00 7.561422e-04\n",
" * time: 1.2719309329986572\n",
" 428 7.592313e+00 7.700220e-04\n",
" * time: 1.274271011352539\n",
" 429 7.592313e+00 5.886958e-04\n",
" * time: 1.2771458625793457\n",
" 430 7.592313e+00 1.242297e-03\n",
" * time: 1.2795889377593994\n",
" 431 7.592313e+00 1.034358e-03\n",
" * time: 1.2818570137023926\n",
" 432 7.592312e+00 1.205310e-03\n",
" * time: 1.2842848300933838\n",
" 433 7.592312e+00 1.119195e-03\n",
" * time: 1.2865538597106934\n",
" 434 7.592312e+00 6.624688e-04\n",
" * time: 1.2888619899749756\n",
" 435 7.592312e+00 1.445924e-03\n",
" * time: 1.2911510467529297\n",
" 436 7.592312e+00 1.369466e-03\n",
" * time: 1.2935190200805664\n",
" 437 7.592312e+00 9.197517e-04\n",
" * time: 1.296483039855957\n",
" 438 7.592312e+00 6.728694e-04\n",
" * time: 1.2993009090423584\n",
" 439 7.592311e+00 8.488112e-04\n",
" * time: 1.3016059398651123\n",
" 440 7.592311e+00 9.128676e-04\n",
" * time: 1.3039519786834717\n",
" 441 7.592311e+00 5.240053e-04\n",
" * time: 1.3068568706512451\n",
" 442 7.592311e+00 1.034452e-03\n",
" * time: 1.3091819286346436\n",
" 443 7.592311e+00 9.269072e-04\n",
" * time: 1.3121328353881836\n",
" 444 7.592311e+00 8.845126e-04\n",
" * time: 1.3146109580993652\n",
" 445 7.592310e+00 8.025110e-04\n",
" * time: 1.316922903060913\n",
" 446 7.592310e+00 6.202968e-04\n",
" * time: 1.3200538158416748\n",
" 447 7.592310e+00 5.620058e-04\n",
" * time: 1.323056936264038\n",
" 448 7.592310e+00 5.136854e-04\n",
" * time: 1.3255388736724854\n",
" 449 7.592310e+00 1.002645e-03\n",
" * time: 1.3279979228973389\n",
" 450 7.592310e+00 5.838946e-04\n",
" * time: 1.3312909603118896\n",
" 451 7.592310e+00 6.229958e-04\n",
" * time: 1.3338639736175537\n",
" 452 7.592309e+00 1.037268e-03\n",
" * time: 1.3361828327178955\n",
" 453 7.592309e+00 7.810620e-04\n",
" * time: 1.3391430377960205\n",
" 454 7.592309e+00 1.455935e-03\n",
" * time: 1.3415579795837402\n",
" 455 7.592309e+00 1.561557e-03\n",
" * time: 1.3440959453582764\n",
" 456 7.592309e+00 9.561142e-04\n",
" * time: 1.3464748859405518\n",
" 457 7.592309e+00 1.151604e-03\n",
" * time: 1.3487999439239502\n",
" 458 7.592308e+00 8.964022e-04\n",
" * time: 1.3598418235778809\n",
" 459 7.592308e+00 1.087445e-03\n",
" * time: 1.3621859550476074\n",
" 460 7.592308e+00 8.456089e-04\n",
" * time: 1.3644330501556396\n",
" 461 7.592308e+00 7.271225e-04\n",
" * time: 1.36673903465271\n",
" 462 7.592308e+00 9.024107e-04\n",
" * time: 1.369246006011963\n",
" 463 7.592308e+00 5.428255e-04\n",
" * time: 1.3721988201141357\n",
" 464 7.592308e+00 9.906250e-04\n",
" * time: 1.374675989151001\n",
" 465 7.592308e+00 6.878895e-04\n",
" * time: 1.377619981765747\n",
" 466 7.592308e+00 3.721208e-04\n",
" * time: 1.3805079460144043\n",
" 467 7.592308e+00 4.070284e-04\n",
" * time: 1.3829419612884521\n",
" 468 7.592307e+00 3.219611e-04\n",
" * time: 1.38588285446167\n",
" 469 7.592307e+00 4.782827e-04\n",
" * time: 1.3881959915161133\n",
" 470 7.592307e+00 4.801304e-04\n",
" * time: 1.391122817993164\n",
" 471 7.592307e+00 5.861392e-04\n",
" * time: 1.393460988998413\n",
" 472 7.592307e+00 7.772588e-04\n",
" * time: 1.3957159519195557\n",
" 473 7.592307e+00 5.881656e-04\n",
" * time: 1.3980779647827148\n",
" 474 7.592307e+00 5.232730e-04\n",
" * time: 1.4003100395202637\n",
" 475 7.592307e+00 7.462374e-04\n",
" * time: 1.4027659893035889\n",
" 476 7.592307e+00 6.163629e-04\n",
" * time: 1.4049999713897705\n",
" 477 7.592307e+00 6.590480e-04\n",
" * time: 1.4072790145874023\n",
" 478 7.592307e+00 7.110202e-04\n",
" * time: 1.4095098972320557\n",
" 479 7.592307e+00 8.450639e-04\n",
" * time: 1.4119348526000977\n",
" 480 7.592307e+00 5.254595e-04\n",
" * time: 1.4141738414764404\n",
" 481 7.592307e+00 5.403903e-04\n",
" * time: 1.416496992111206\n",
" 482 7.592307e+00 7.955876e-04\n",
" * time: 1.4188239574432373\n",
" 483 7.592307e+00 8.131631e-04\n",
" * time: 1.4218459129333496\n",
" 484 7.592307e+00 8.033843e-04\n",
" * time: 1.4246819019317627\n",
" 485 7.592306e+00 9.001554e-04\n",
" * time: 1.4271318912506104\n",
" 486 7.592306e+00 6.425147e-04\n",
" * time: 1.4300038814544678\n",
" 487 7.592306e+00 6.873360e-04\n",
" * time: 1.4323358535766602\n",
" 488 7.592306e+00 1.030023e-03\n",
" * time: 1.4351859092712402\n",
" 489 7.592306e+00 5.899199e-04\n",
" * time: 1.438060998916626\n",
" 490 7.592306e+00 6.312570e-04\n",
" * time: 1.440269947052002\n",
" 491 7.592306e+00 6.195972e-04\n",
" * time: 1.4428019523620605\n",
" 492 7.592306e+00 4.270171e-04\n",
" * time: 1.4537379741668701\n",
" 493 7.592306e+00 3.880703e-04\n",
" * time: 1.4559988975524902\n",
" 494 7.592306e+00 3.948447e-04\n",
" * time: 1.458225965499878\n",
" 495 7.592306e+00 4.855579e-04\n",
" * time: 1.460433006286621\n",
" 496 7.592306e+00 4.496147e-04\n",
" * time: 1.462751865386963\n",
" 497 7.592306e+00 4.425022e-04\n",
" * time: 1.4656648635864258\n",
" 498 7.592306e+00 7.195816e-04\n",
" * time: 1.468066930770874\n",
" 499 7.592306e+00 4.748474e-04\n",
" * time: 1.4703538417816162\n",
" 500 7.592306e+00 4.960634e-04\n",
" * time: 1.4727048873901367\n",
" 501 7.592306e+00 6.709209e-04\n",
" * time: 1.4750800132751465\n",
" 502 7.592306e+00 3.651232e-04\n",
" * time: 1.478022813796997\n",
" 503 7.592306e+00 5.944828e-04\n",
" * time: 1.4802608489990234\n",
" 504 7.592306e+00 4.468584e-04\n",
" * time: 1.4832520484924316\n",
" 505 7.592306e+00 5.799489e-04\n",
" * time: 1.4854748249053955\n",
" 506 7.592306e+00 8.113803e-04\n",
" * time: 1.4876840114593506\n",
" 507 7.592306e+00 1.043600e-03\n",
" * time: 1.4898948669433594\n",
" 508 7.592306e+00 7.190400e-04\n",
" * time: 1.4922139644622803\n",
" 509 7.592306e+00 5.879789e-04\n",
" * time: 1.4951839447021484\n",
" 510 7.592305e+00 6.571553e-04\n",
" * time: 1.4975030422210693\n",
" 511 7.592305e+00 9.001956e-04\n",
" * time: 1.4998040199279785\n",
" 512 7.592305e+00 7.378372e-04\n",
" * time: 1.502089023590088\n",
" 513 7.592305e+00 3.964742e-04\n",
" * time: 1.5044119358062744\n",
" 514 7.592305e+00 7.574339e-04\n",
" * time: 1.506654977798462\n",
" 515 7.592305e+00 7.486537e-04\n",
" * time: 1.508997917175293\n",
" 516 7.592305e+00 9.885412e-04\n",
" * time: 1.511214017868042\n",
" 517 7.592305e+00 6.927226e-04\n",
" * time: 1.5143029689788818\n",
" 518 7.592305e+00 4.341122e-04\n",
" * time: 1.5172319412231445\n",
" 519 7.592305e+00 4.593987e-04\n",
" * time: 1.51951003074646\n",
" 520 7.592305e+00 5.829484e-04\n",
" * time: 1.5217769145965576\n",
" 521 7.592305e+00 7.709669e-04\n",
" * time: 1.5240809917449951\n",
" 522 7.592305e+00 9.223135e-04\n",
" * time: 1.5269789695739746\n",
" 523 7.592304e+00 4.144980e-04\n",
" * time: 1.5299839973449707\n",
" 524 7.592304e+00 6.367513e-04\n",
" * time: 1.5323128700256348\n",
" 525 7.592304e+00 7.192259e-04\n",
" * time: 1.5345749855041504\n",
" 526 7.592304e+00 5.958315e-04\n",
" * time: 1.5460779666900635\n",
" 527 7.592304e+00 6.007328e-04\n",
" * time: 1.5484249591827393\n",
" 528 7.592304e+00 5.538829e-04\n",
" * time: 1.5507190227508545\n",
" 529 7.592304e+00 6.352578e-04\n",
" * time: 1.5531208515167236\n",
" 530 7.592304e+00 9.156473e-04\n",
" * time: 1.5553369522094727\n",
" 531 7.592304e+00 1.113372e-03\n",
" * time: 1.5577068328857422\n",
" 532 7.592304e+00 1.270301e-03\n",
" * time: 1.5599958896636963\n",
" 533 7.592303e+00 8.576146e-04\n",
" * time: 1.5623078346252441\n",
" 534 7.592303e+00 8.850594e-04\n",
" * time: 1.5646400451660156\n",
" 535 7.592303e+00 1.627696e-03\n",
" * time: 1.5669519901275635\n",
" 536 7.592303e+00 9.423361e-04\n",
" * time: 1.5693318843841553\n",
" 537 7.592303e+00 1.034615e-03\n",
" * time: 1.5715668201446533\n",
" 538 7.592303e+00 8.449986e-04\n",
" * time: 1.5749030113220215\n",
" 539 7.592303e+00 6.295347e-04\n",
" * time: 1.577193021774292\n",
" 540 7.592302e+00 7.176453e-04\n",
" * time: 1.579530954360962\n",
" 541 7.592302e+00 6.586048e-04\n",
" * time: 1.5824990272521973\n",
" 542 7.592302e+00 1.942404e-03\n",
" * time: 1.584751844406128\n",
" 543 7.592302e+00 1.070989e-03\n",
" * time: 1.5877599716186523\n",
" 544 7.592302e+00 9.657404e-04\n",
" * time: 1.5900559425354004\n",
" 545 7.592302e+00 9.402287e-04\n",
" * time: 1.5923488140106201\n",
" 546 7.592302e+00 1.020338e-03\n",
" * time: 1.594667911529541\n",
" 547 7.592302e+00 8.866450e-04\n",
" * time: 1.5969128608703613\n",
" 548 7.592302e+00 6.839219e-04\n",
" * time: 1.5999388694763184\n",
" 549 7.592302e+00 9.526252e-04\n",
" * time: 1.6022298336029053\n",
" 550 7.592302e+00 1.080698e-03\n",
" * time: 1.6045429706573486\n",
" 551 7.592301e+00 8.966271e-04\n",
" * time: 1.6068480014801025\n",
" 552 7.592301e+00 7.540693e-04\n",
" * time: 1.6091740131378174\n",
" 553 7.592301e+00 1.231184e-03\n",
" * time: 1.6114189624786377\n",
" 554 7.592301e+00 1.019640e-03\n",
" * time: 1.6138370037078857\n",
" 555 7.592301e+00 9.502074e-04\n",
" * time: 1.6161129474639893\n",
" 556 7.592301e+00 7.177609e-04\n",
" * time: 1.619182825088501\n",
" 557 7.592301e+00 6.939495e-04\n",
" * time: 1.6221868991851807\n",
" 558 7.592301e+00 8.901981e-04\n",
" * time: 1.6245129108428955\n",
" 559 7.592301e+00 8.014725e-04\n",
" * time: 1.62680983543396\n",
" 560 7.592301e+00 6.701744e-04\n",
" * time: 1.6298069953918457\n",
" 561 7.592301e+00 4.573906e-04\n",
" * time: 1.6411418914794922\n",
" 562 7.592301e+00 5.266435e-04\n",
" * time: 1.6434438228607178\n",
" 563 7.592300e+00 8.559942e-04\n",
" * time: 1.6456730365753174\n",
" 564 7.592300e+00 4.540772e-04\n",
" * time: 1.6487088203430176\n",
" 565 7.592300e+00 5.476328e-04\n",
" * time: 1.6509678363800049\n",
" 566 7.592300e+00 6.307882e-04\n",
" * time: 1.653291940689087\n",
" 567 7.592300e+00 3.943299e-04\n",
" * time: 1.6555328369140625\n",
" 568 7.592300e+00 5.727576e-04\n",
" * time: 1.6577799320220947\n",
" 569 7.592300e+00 3.645821e-04\n",
" * time: 1.6606390476226807\n",
" 570 7.592300e+00 5.189390e-04\n",
" * time: 1.6630849838256836\n",
" 571 7.592300e+00 3.436121e-04\n",
" * time: 1.6653549671173096\n",
" 572 7.592300e+00 5.047516e-04\n",
" * time: 1.6675758361816406\n",
" 573 7.592300e+00 3.588241e-04\n",
" * time: 1.6697778701782227\n",
" 574 7.592300e+00 3.213661e-04\n",
" * time: 1.6720948219299316\n",
" 575 7.592300e+00 3.363915e-04\n",
" * time: 1.6745679378509521\n",
" 576 7.592300e+00 4.599944e-04\n",
" * time: 1.6768219470977783\n",
" 577 7.592300e+00 3.321465e-04\n",
" * time: 1.6796610355377197\n",
" 578 7.592300e+00 3.555291e-04\n",
" * time: 1.6819629669189453\n",
" 579 7.592300e+00 2.982073e-04\n",
" * time: 1.6850378513336182\n",
" 580 7.592300e+00 3.049252e-04\n",
" * time: 1.6872608661651611\n",
" 581 7.592300e+00 2.045465e-04\n",
" * time: 1.690194845199585\n",
" 582 7.592300e+00 1.968502e-04\n",
" * time: 1.6924188137054443\n",
" 583 7.592300e+00 1.599554e-04\n",
" * time: 1.6948168277740479\n",
" 584 7.592300e+00 1.835117e-04\n",
" * time: 1.697026014328003\n",
" 585 7.592300e+00 1.217371e-04\n",
" * time: 1.6999719142913818\n",
" 586 7.592300e+00 1.793300e-04\n",
" * time: 1.7022149562835693\n",
" 587 7.592300e+00 2.289451e-04\n",
" * time: 1.7046480178833008\n",
" 588 7.592300e+00 1.842698e-04\n",
" * time: 1.706885814666748\n",
" 589 7.592300e+00 1.771351e-04\n",
" * time: 1.709928035736084\n",
" 590 7.592300e+00 1.610654e-04\n",
" * time: 1.7121479511260986\n",
" 591 7.592300e+00 1.221196e-04\n",
" * time: 1.7152299880981445\n",
" 592 7.592300e+00 1.481587e-04\n",
" * time: 1.7174479961395264\n",
" 593 7.592300e+00 1.301925e-04\n",
" * time: 1.7198619842529297\n"
]
},
{
"output_type": "execute_result",
"data": {
"text/plain": "Plot{Plots.GRBackend() n=1}",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3dd1xUV/o/8DMMVYoooICKKIKgCBp7iSWKPbZY1hbrmpjYkhg1u9mY2KImP5Os+SZqrGuJscSeoqyxr1Ew1kVRREEBFRh6H+b3x81OJvA8OHe8MA7zeb/4YzicuXNm5g5nzr3PfR6VTqcTAAAA1srG3AMAAAAwJ0yEAABg1TARAgCAVcNECAAAVg0TIQAAWDVMhAAAYNUwEQIAgFXDRAgAAFYNEyEAAFg1TIQAAGDVbKvmYUpLSy9dvR4YElI1D1fFSktLbWyq/1cKPM3qBE+zOqmyp1nTwe6pfW7dSs/PL5G12ezsTF9fm4CAAFPH9ayqaCIsLCz86rt9Ayd6Vc3DAQCA4oY19Xlqn1GjDl258kTmhq8PGZK5b98+00b17Kr/dyUAAIAKVNGKEAAArIFKJVQqcw9CJqwIAQDAqmFFCAAAilGpVCqZS0KdzsxLSEyEAACgKBwaBQAAsCBYEQIAgHIsMFgGEyEAACjGhHOEwtznCHFoFAAArBomQgAAsGo4NAoAAIox4YJ6cx8ZxUQIAADKMeEcodmDa3BoFAAAnmu//vrryJEj+/bt+/XXX+t0uvIdUlNTZ86cGRER8dZbb2k0GqkxMjLy9ddf79Onz/jx48+cOaPvnJWV9e6770ZEREyfPv3Ro0cCEyEAAChJJf+nQg8fPuzdu3f37t3ffffdVatWrVu3rnyf4cOHZ2Zmfvjhh8nJyWPHjpUa169fHxwcPG/evBdeeKF3794XL16U2idPnnz79u0PP/ywtLT05ZdfFjg0CgAAClLJP9RZcf/169f36tXrjTfeEEJ8/PHHH3zwwWuvvWbY4fLly9HR0T/99JOjo2N4eLiXl9etW7eaNm26c+dOqUPPnj1Pnz59+PDhtm3bJiQkHDp06OHDh56enu3bt/f29j579ixWhAAAoBiVfBXPhL/99lunTp2k2506dYqJicnLyzPscOnSpVatWjk6OgohXFxcWrRo8dtvvxl2KC0tvXPnTsOGDYUQV65cady4saenpxDC1ta2ffv2ly5dwkQIAADPr5SUlFq1akm3a9euLYSQTuzpPXr0SN9B6pOSkmLYYfny5UKIcePGcZ1xaBQAAJRjxGm/8nc5duxYQECAYduYMWMWL14shHB1dc3Pz5capbWgq6urYU8XF5eCggL9r3l5eYYd1q9fv3bt2lOnTtnb23OdMRECAICS5M+DolOnTmvWrDFsrFOnjnTDz8/v7t270u27d+86Ozt7eHgY9jTsoNPp4uPjpaOgQoh//etfH3300fHjx/Utfn5+9+/f12q1arVa2uDw4cPlHRq9fPny7t27T5w4YTijXrlyZffu3Tdu3JC1KQAAAImzs3PjP3NxcZH+9Je//GX37t3SRRHr1q0bNWqUdJ3id999J50L7N27d3p6+vHjx4UQP/zwg1ar7datmxBi9+7dCxYs+PnnnwMDA/UP1KFDBzc3tz179gghzp8/HxcXN2DAAGMnQq1WO378+EGDBu3Zs+fDDz/ctm2b1P7xxx8PGDDgyJEjvXr1Wr16tVIvCgAAWCKVSn68TIVryF69evXr1y8kJCQsLOzs2bMfffSR1L5q1aqTJ08KIZycnNasWTNy5Mi2bdtOnDhxzZo1dnZ2Qoh58+ZpNJouXbrUrl27du3as2fPFkLY2NisXbt25syZbdu2HTBgwP/93/+5ubmpyIsTy/vyyy/Xr19/+vRp6dirTqdTqVRpaWl+fn5RUVEhISFRUVG9evV68OCBfho3lJ+f/+ailQMnTpP5kgIAwPNiWFOfp/bp0Hn7tWupsjarLbnWr0/avn37KuiTnJycmZkZFBRkY/P7+q2oqEitVktHOIUQOTk59+7da9SokbOz81MfMT8/Py4urmHDhtKMZuyKcNu2bbNnz05LSzt//nxubq60Mj127FhgYGBISIgQok2bNh4eHtL8DAAA1smEBaHKiAsPfXx8goOD9bOgEMLe3l4/CwohXFxcQkNDjZkFhRBOTk6hoaH6mBpjg2Xu3Lmzc+fONWvWuLq63rx58/Dhwy1btnz48GGDBg30fRo0aPDgwQNuC0YuPQEA4PkkHQs09yiUZ+yKUIoxPX/+fGRk5IQJE+bOnSuEKC4uNpyQbW1ti4uLK2WYAABgbsasZ6Tr42X9mJ2xE6Gvr2/Pnj2l7wIRERFXrlwRQvj4+Dx58kTf5/Hjxz4+7BHkavk9AgDAehgemWSZkGvU3JODsRNh9+7d4+Pjpdvx8fG+vr5CiC5duly6dEmKak1OTo6Nje3YsWMlDRQAAKAyGHuO8O233+7WrZu7u7ubm9vSpUs/++wzIURAQMDgwYOHDBkyduzYTZs2jR8/XpogAQDAOknHO2XexcyMXRE2a9bs9OnTeXl5SUlJ+/bt+8tf/iK1b926dfTo0Tdu3PjrX/9aJi8AAABYHfnnCM1+3kxGirXg4OAlS5aUabSzs3v99dcVHRIAAFis5+Ccn1yoPgEAAFYNSbcBAEAx0gX1Mu9j5iUkJkKwfNSlTVaevoH9v2Jpx6zA4phwzs/c8yAOjQIAgHXDihAAAJRjgcEymAgBAEAxRibRLnOXShqMkTARAgCAYixwQYhzhAAAYN2wIoQqwQRxVmKzvG0I8fyXCWMOH5GtOuY7uYp7ljK2jZBUqJAJS0Jz7zmYCAEAQDkmnCOspJEYDYdGAQDAqmFFCAAAilHJv0De3EGjmAgBAEA5KpX8yyHMPRNiIgRTyUpsxoSikM06OZ25/rpSGZ3lbrxScf9EuP8VZH8Vc9JD1sbZ/05yomUQWWN1LDBYBucIAQDAqmFFCAAAirHABSEmQgAAUI4pKdbMPRXi0CgAAFg1rAgBAEA5FnhsFBMhPAUbNUnGajKdS0vpP5DtpVoZnbn+sjqLCqJJqe0oEkkqKwpUCKGyodtt1ES7jZzOXH9ZnYUQNuQBJplho+YOpIdnZcJ1hGafCHFoFAAArBpWhAAAoBwLzDWKiRAAAJSDc4QAAGDNVCrLyzWKc4QAAPC8y8rKSkpKqqBDcXFxQkJCcXGxCRvHitD6sBVrubBJuj8d8MnEampL6K1oS4j+JcV0Z66d3LhWS3dmQ1K5djodKtmXjSaV9YVXdsAn1a5W099x1bZ0u60d0U42CiHUtvRIyI3LCzEV7Jdz+pprc68koDyVUMm9QP6p/d99993169e7ubl5e3sfOnSoTp06ZTocOXJk0qRJ7u7uWVlZ27Zt69WrlxDi6NGjS5YsuXz5cnBw8IULF6SeZ8+e7dKli+F9f/rpJ6wIAQBAOSr5PxU6efLktm3bbt68ef/+/cDAwA8//LBMh8LCwsmTJ2/cuDE2Nnb16tWTJ0/WarVCCE9Pz7lz537wwQeGnTt37qz7n59//tnDw6N79+6YCAEA4Pm1ffv2UaNG1a1bVwgxY8aM7du3l7nq99ixYzVq1Bg4cKAQYvjw4SUlJadOnRJCvPDCC4MGDZLuSNqwYcOrr77q4OCAiRAAABQjBcvI+ql4URgfHx8YGCjdDgoKysrKSk9PL9MhKCjof4+uCggIiI+Pf+o409LSDhw4MHHiRIFzhAAAoCCVkH8doUr1+PHjyMhIw8YmTZr4+/sLIbKzs52cnKRGZ2dnIURmZqaHh4e+p2EHqU9WVtZTH3Tbtm3h4eFhYWECE2H1RsduMBEdTJiLKGWiTrTFxB2Ki7Vk5+IieiPFRUR/slEIUVwoI4iGi6xhnw6Xeq3qU6wx0SVqNlhGRpwL125HHRyys1fTndl2aiN2dGe1HfN0mGgZGxX5oiNP2/PHhOsIhbh9+/aKFSsMWwYOHDh79mwhRJ06dTIyMqRGjUYjhChztNOwg9SnfDRNeZs3b54+fbp0GxMhAACYWefOnfft20f+KTw8XB/z+euvvzZp0kRaFxp2mDt3blFRkb29fV5e3vXr18PDwyt+uKioqFu3bo0cOVL6FecIAQBAOarfSxLKUdH2pkyZ8sMPP2zfvj0qKur9999/8803pfZx48Zt375dCNG2bdvg4OC33377+vXrs2fPbt++ffPmzYUQjx492r1794ULF9LT03fv3n3mzBn9Njds2DBixAh3d3fpV0yEAACgGKWvnhD+/v4HDhzYunXrrFmzxo8fP2vWLKm9UaNGtWvXlm5///33OTk5EydO1Ol0O3fulBqlifDRo0cvvPBCmYlQq9XOmDFD/ysOjQIAwHOtR48ePXr0KNO4ePFi/W1fX9/NmzeX6RAWFrZr1y5yg+vWrTP8FRMhAAAoB0m3wSxk1c7ls6DR7WwgaCHRXlRAdy5k2sn+RQUldGfqEQUbeiovxZrMqFGuui/ZTAcxyi3Ay0eNEu12bHQoHcNpT7XbO9L/HxwcmY1Q7VpH+o2wK2VCT+3IZiGovG5cnjYdoknN5+kn/crfxdwzISZCAABQjAkV6s3+BQXBMgAAYNWwIgQAAOXgHCEAAFgzE84Rmv3YKA6NAgCAVcOK0JJwQYlsECMV8cgl4eTCLAvz6RjOAqqd61yYJ2MjBUxnLiSVjBplq/jKLcxLRo0ykbfcG0F+QeaiQ+UW5iWjSdmcokyaUDLg07EG/f/B0YmJJqX6OxQzG2GilHVOZLOw08lIqapiv+FTb4S5j8tVM78XlJB1l8oZifGwIgQAAKuGFSEAACjHhOsIzb0kVGBFWFRUFB8fX1xc/OybAgAAqGLGToRz5841TBWel5cntR88eLBevXr9+vVr0KBBmbKKAABgbRSvUF8FZBwa/fvf/75kyRLDlsLCwilTpmzZsqV///67du2aPHlyfHy8Wk2fjQdZyKgLWUExggkYIVOjCSZuRQiRn8u05xDHAPJz6QMD7Eao/lywDDdyujBvibwUa1zmObI/11kWuUExsvrb2sotzCsjWMbJmU6D5kS9QSUyc91xO7nQyfhnxT1N5ps/8rEpzdJeOnmHRvPz8w1/PXr0qJubW//+/YUQw4cPLyoqOnXqlJKjAwAAqGQyJsJPP/3Uy8vL29v7//2//ye1xMfHBwUF/b4hG5smTZrcu3dP8SECAIClkF+VV/4F+EozdiJ8/fXXNRpNTk7O/v37ly5dunfvXiFEdna2k9MfV/3UqFEjKyuL20JpKX2EBAAALIJWS5+eMGTCOUJzz4NGT4RNmjSR5rwOHTqMHz/+p59+EkLUqVNHo9Ho+2g0mrp167KPxFVMAQAAS2BMCIgJ86DZTyqaMjlpNBpnZ2chRHh4+G+//VZUVCSEyM3NvXHjRnh4uMIDBAAAqEzGBmItXLjwxRdfdHV1PXHixO7du8+cOSOEaNeuXVBQ0FtvvfXaa699/vnnHTt2DAkJqczRVkdyaurKig4VTCVbLttZHhPwmZdNt+dmFRGdqVBSwYSYCiHyqcFwhXnZrGlUsi42apRpl5V6jY1slIM7LyIrlZoQwoYKEOWiRtVUeVshhG0B0Z+thMxWSCaiSbl3jYu85aNGjWyrCBlNKisfWwXN8DvzL/BkM3ZFqNPpli9fPmfOnKtXr544caJ169ZS+/79+/Pz86dOnWpnZ/fdd99V2jgBAMASWOA5QmNXhIsWLSLbfX19N27cqNx4AAAAqhRyjQIAgGJU/GF/9i6WsiIEAAB4Ogs8R4iJ0My4s/1kJAEZFSJklhLk4llys4ngFyFEbpacYBlmI2RQjGBKDJL1BQX/9Mm4mBJmI1wQDbdxLhMYjfv8y9kGGyzDxbmQwTJM3UEuiIYMaeFeEy0bcES086nUyGZ5wTJsOAuzxKDb7ZgoJMTKmMSUeoTmfk1xbR8AAFg1rAgBAEAxpqRMM/eSECtCAACwapgIAQDAquHQKAAAKMaUYJnKGYnxMBFWETYUjqsHSwXmceGUZHSoYMrh5jIp03IyuahRpp3aDpdKjRshGe9aUqxAwKeWyezFZfxi87dRkZA8BcJG1Wp5NXVL7ainz8RqlsgJPZUbNcoUMSb78lWm2fR1xMhVTGQnHzVqbKMQQsVE2JLhpOY+yfUcUQnZ5wgtpgwTAABAtYQVIQAAKMeEC+rNvZ7GRAgAAMqphCTaGo1m586dGo1m4MCBYWFh5TvodLq9e/fevHkzNDR08ODB+mOtd+/evXbtmouLS8+ePfWdDxw4IFUPFEI0aNCgQ4cOODQKAACKUZmkgg3m5OS0b9/+xIkTBQUF3bp1+/e//12+zxtvvLF06VIbG5v3339/7ty5UuMXX3zRpk2befPmvffee4adp0yZsnfv3sjIyMjIyKtXrwqsCAEA4Hm2Y8cOT0/PnTt3qlQqLy+vpUuXGi7vhBAPHjzYtGlTfHy8j4/PmDFjQkJCFixY4OXlNW3atNmzZ2/dunX16tVltrl48eLAwED9r5gIlUdGvXGxcFx5UjKIsZDKzCmEKGAyeZIxnFw6UC6alMs1mpdDbIcbCZlTVDBxsFwAJ8fegcir6eBIJ9tkR8gUmyXbuVhfGfGOfMSjgyPdbs88I8caxKeYy/DJPc0SKnzX1o4L06VHQu7M8nKHCjb+0oZ6udhixVw0KdUuq7Ngc5DKiFMFWSIjI/v37y+tGvv37z9nzpzi4mI7uz+qQJ84caJFixY+Pj5CCH9//yZNmpw+fXrYsGFOTk7cNg8dOuTh4dG6devQ0FCBqFEAAFCQ0kdGRVJSkre3t3Tb29u7tLQ0JSWlTAdpFtT3SUpKqmCDzZo1u3nz5r///e/OnTtLpXaxIgQAADO7fPnya6+9Ztjy4osvjhs3TghhY2Oj1f5+JKO0tFQIoVb/6ZiEWq3WdxBCaLXaMh3KOHXqlP5B27RpM2nSJEyEAACgGBMyywiVcHd3b926tWFbs2bNpBu+vr76JWBycrJara5Tp45hTx8fn+TkZP2vKSkpvr6+xjxsy5YtPTw8bt++jYkQAAAUJf/MqL+//7Rp08g/9evXb/Xq1f/4xz9sbGwOHDgQERFha2srhIiNja1Vq5aXl1fPnj2nTJly7949f3//W7du3b9/v1u3btwDGa4Xr127lpaWFhgYiInwGXBn+6noAC68gstWVURFkfCp1Lh4FqKdC4rJY1KpkUExgsnfVpBHb5wLlpEVF1OztiPZ3iCwZvlGl5oOZOdHidlke3yMhmwnX0Mu5ITNpEd9QyZjfAQT/CKEaNCEeJpCiLoNXMs35mQWkp0Tb2eS7ZnpBeUbFUlHx2ZMY3CLCTpYhsmCxhU3Jtu5isfcRlQq4kGptt//Iqu5GlCpZKdMq7j/yJEjP//88379+gUGBu7YsePIkSNS+/jx40ePHj1nzpy6devOmjWrd+/egwYN+v777+fNm+fu7i6EiI6OXrBgQXJycmJiYkRERIcOHRYvXvzzzz8vXLiwbdu2+fn5+/bte//99xs0aICJEAAAnl9OTk7nzp07ePBgRkbGpUuX/P39pfZ//vOf+iCaFStW9O3bNyYmZsuWLS+++KLU6O/vP3/+fP12PD09hRAvvfSSg4NDbGxsjRo15s2bFxISIhAsAwAASjKh+sTT+js5OY0aNapMY/v27Q1/7dGjR48ePQxbPDw8evXqVeZejo6OPXv2LHMlIiZCAABQjgm5Rs0N1xECAIBVw4oQAAAUoxIqlaWVn8BEaDo5QaNsyisuAK+YCrOUHTVKZVPjUqyR4ZGCz0lGDqYwn0lUxkWNUkGz9g70gQp3LzpqNPgFr/KNXr7OZOebl+h9Pvk+HU1KvkFFBfRromNiYMmQQjWTqsvJ2Y5s9w+uRbYHtyae/pOHuWTnLA0dTZqRll++saiQiw7lauoSjbJq5FbQn4zhtGUCPrlAULK4sR1T8ZgsViyY+FU27tHqgkZNqlBv7pcDEyEAACjKwhaEOEcIAADWDStCAABQzNOzaJe/SyUNxWiYCAEAQDGm5Ro1LxwaBQAAq4YV4bPgYueIdi7QjqxMK5gavGwBXirtJ9fObaSACUnlQlXJQNBiJgkn9zTJmEzu2yEXekqOkH2aSpQOLqbK2Ao+oyxZ4lVtK+PpCP4NKiDfZaYz+wZRAaJcZ/ZpUo1s2VsunJLpTwd82tP5Wm2ZdrI/vxF6nUAmOLVRy3hNnvIXS2fCBfXmfjEwEQIAgHJMOUeI6wgBAKC6sMAFIc4RAgCAdcOKEAAAFGPC5RNmTy2DifDpuMqiXHspFUrApVLjsliRERNcKrUCrp0KDOGDX+RFkZDFaYuL5QXLkFEnXOeH8VlkO/mpc3ajE5U9epBDtqem5JHt5MtVwgXLMPsEOcJCFf2CpzEjuXHhEdmemkRkU8thyixzr2EulXiP3/O5mCCyUV6wDFcOt1BOsIwdE+dCtnMVku2Ydls74t1Xc1V8mafPMfeMoBBLOzaKQ6MAAGDVsCIEAADFmJJ0u3JGYjxMhAAAoBiViq/Fwd/FvDARAgCAcizw+gmcIwQAAKuGFaExmFRqTKIpMpsaW4CXyWJVSMVwssm35NTO5aNAubhWBXKPce1kqCqXjo57+plpBeUbueBD7rXKzaTDLMlcZdy7yZVfJgej1dJfQzVPiBq5QogiJpg28U6m8SPhXkNtCdGfK2/L/d8gIyRVKnrYbNQoUw7Xzo5KOmhPPx0uatTekQgEJT9rQgiHQvppllDRpFw+NhsuHR0dkSrMvzh6Zha4IMRECAAACrLA6whxaBQAAKwaVoQAAKAosx/rlEneirC4uPjll18eM2aMviUmJuall17y8fHp27fv3bt3lR4eAABYEuk6Qrk/5iVvRbhs2bKEhAQbm9+nT51ON2zYsHHjxu3du/fTTz8dNWrUxYsXK2GQVUVmKjUdHTAhSqm4gxIm9xgXikLWIyzMl9FZCFFUSEQBsKUB5ZQMFEyaMb6zjBRr3Ea46KTcbCLDHBctoi3hwnnkjJDpTGbXE0xNPq5mHjdCbldhyuMxGb/IVibOxbaUCQCxoUdCPqgN80+Oq1PI7yrEYIqpPVwIUVhAv7bkJ4j/+NDtDtRnubSEfkQdExSjY9Yg8t6255IJuUZln1NUmowV4bVr1w4ePPj222/rW06dOpWWlrZgwYJatWp98MEHt27dunTpUiUMEgAAoLIYOxGWlJT89a9/Xb16tZ3dH7mMY2JiwsPD1Wq1EMLBwSEkJOS///1vpQwTAAAsgkr+j7kZOxF++umnHTt27NSpk2FjWlqaq6ur/ld3d/fU1FRuC6WlzJFEAACwBFotfbjYkCknCM09Fxo1Ed67d+/LL7+cNGnS3bt3Hz9+XFhYePfuXZ1OV6tWrZycP4raZGZmenh4sI9kg0s1AAAsmHT8r/oxKlgmNTXV29t78uTJQoj09PRHjx6NHDny9OnTTZo0iYmJ0el0KpWqpKTk9u3bAQEBlTxgAAB4fknLPJl3MTOjJsI2bdpERUVJt3fs2PHJJ59Iv7700ksqlWrTpk2TJ0/+6quvvLy8OnbsWImDrWRMcCj7h1ImnFSrlREJyQUrkonQZGU7E0KUFBMj5BKYcRGPXBVWcufldmhm2/RrxYVNciMnX1suTJfNjkbF+nIPSg5b8AHGKhXxB+5cgZat4Uw/I1uqYq2tHVdUlj4wQ9aVJQvtCiFKmVBIWRGP3H7F7YdM5kK6M/uxol5DLukg/7GidlpmlyhlwkNtdPTrQjabfZ6QxwJzrMk+XOni4lKvXj3ptq2t7XfffffJJ584OzuvX7/+22+/NXsULAAAVDP79+9v27Zt06ZN58+fX1xMXC4VFxc3cODAgICAIUOGJCQkSI0xMTEffPDBoEGDZsyYoe+ZlJQ0a9asF154oVmzZhMnTnzw4IEwYSIcNGjQ4cOH9b927NgxJiYmJyfn6tWrrVq1kv38AACgWlHJVuGSMDY29tVXX120aNGPP/544sSJTz/9tEwHnU43dOjQ8PBw6YTdyJEj9XfMycmpU6fOhQsX9J0TEhJq1qy5bt26/fv3FxcXDx48WCiVYg0LQQAAECZVqK/40OiGDRsGDx7cr18/IcTChQunT5/+3nvvGXY4d+7cw4cPP/roI1tb26VLl3p6el6+fLlly5aDBw8ePHjw1q1br169qu/coUOHDh06SLeXLFnSuHHjrKwsRHICAIByTLiOsMKJ8Pr1661bt5Zut27dOiEhITs727DDjRs3wsPDbW1thRAODg6hoaHXr183ZqTnz5/38/Nzc3PDRAgAAM+v1NTUmjVrSrfd3d2FEE+ePCnTwc3NTf+ru7t7mQ6ke/fuzZkzZ/Xq1QLVJ4zBRZNyqS/JLJeyIh4FV/aWiRtkIyGpEXJJODmONezIdpeaxLcoLgBPpDLFZqmaqNwLyz1NciPcSLinL+tF4c4FyDoixIVHcu0qIkRACCbpq70D94Toj7wNFSHKpQN1qkFvxN3TiRgJVQtXVJBRlosmpd447rWSVQdbVnSoYD7L7H7FBWOTreYPn1SAKblGhWr//v1l7jV79uzPP/9cCGF4wbq0Fqxdu7ZhT3d399zcXP2vWVlZZTqU9/Dhw169er3//vuDBg0SmAgBAEBBpl09MWTIkH379pF/bdy48c2bN6XbN2/erF27trQuLNNBuqK9tLT09u3bjRs3ruDhUlJSevbs+frrr8+cOVNqwaFRAAB4fr366qu7du26d++eVqtdtWrVq6++KrV/9tlnJ0+eFEL06tVLp9N9++23QohNmza5ubl17txZCFFcXKzRaPLy8kpKSjQajbSsfPz4cc+ePV9++eUpU6ZoNBqNRqPVarEiBAAA5Sh9QX2HDh3eeeedli1bqlSqTp06LVy4UGqPjIx0cXHp1q2bra3tjh07JkyYMGvWrFq1am3fvl3K6HnmzJlXXnlF6hwQENCnT59vv/323LlzycnJGzZs2LBhg/SnCxcuYCIEAADFqEy4oO5p3RcsWPDuu+8WFxc7OjrqG48cOaK/3bVr1/j4+NzcXGdnZ31jjx490tPTy2xqyJAhQ4YMKdOIidAISuSCYtOGMe102jCmM3eiXlZ8AZd8y8vXmfr0K8gAACAASURBVGz39Xcr31iYTxdKvX01jWwvyCP6cxvha+rKeK047EeXevflhRtxj8gPhWzmUvqRyeHIQrtCCLUd87I4EO++gyP9/8Hbz5VsDwwncu5zG0m6l0W2P0nKJdtlRX5x7eRewQbFyPnMys1cyKbjs/xwGROuIzSmu1qtfmrKb8NZUBacIwQAAKuGFSEAACjHApNuYyIEAADlmHAdobmTdOLQKAAAWDWsCAEAQDGKJ92uApgIDXGVaZneXMFeOeFtWqYeLBkKyHaWk6yLy/nEpcLybUREhwohwjv5lG/kAj5zs+n8YCkJ2eUbszRMzVIuWxUXmEexYT6j3CbIbXO5x7h2epdghs1sgx05GU0q/7UiNu7IpFKTtUs4MBvh/vdlphWS7WQiPblp6shPEFuTmfu4yQnGlvvfg9kTzT1RVHeYCAEAQDEqYUqu0UoajJFwjhAAAKwaVoQAAKAYUy6oN/ehX0yEAACgGFPKMJl7JsShUQAAsGpYERqQm0GSiweTEyIoL+BTTs1SIUSplkiKKDe8jYudI2ucFjGlg7m0jeQzYiseyxk59x2Ty9ip43JwOhHBtHXru5Cd3Wo7ku1Z6QXlGx89yCE7F+bTryFVOlcI5pnKCmkWzGvO7W/cu0kWQ+aynnL7FbuT08HY8rLvkhtX5rPJ5hSlm1nVIGjUhMwy5oaJEAAAlIPrCAEAwJqphEru5RDmngdxjhAAAKwbVoQAAKAcVJ+wGEpUVtXJ2QoXi8GFBnAPyWxcRhANF6SQx2RBu39LQ7bn5xSVbyxgAj3ux2aQ7WTqNbn1h2W9hlxf7qxGnXpEXEyXgf5k58AWRGVaIcTta0Rd4tOH75GdH9zJJNtljZyNc5FTsZZLjHf/Fv1ukiPhUqw9eUgX4OX2QxL3dNhdQtanTWZcFfOAyvyvoZl7/iCp5F8XaPbngUOjAABg1ax1RQgAAJVApZJ/gby5L6jHRAgAAMrBOUIAALBmpuQarZyRGA/nCAEAwKpZ64qQ/AYiM7xLmRpacuLBuCRbHFlZqYqK6Pbk+0TtXCHEo0QiQxhXmDc7gwgxFUIU5BH9uQxespJycWcpuO+qXE1d11oO5RsbN6tNdm7R0ZtsJ4MbL59JpkeillHdV7BRo2TfCorQEq95QS79bibfzyLbczKJmroOTvQ/Gbmp/uztib1fbmIz+hPEfY6ViPdU5r+E2VdMciDpNgAAgIXBRAgAAFbNWg+NAgBA5UDSbQAAsF6mnCOspKEYDRPhM2DePXIf4HYMPqaDaOdKu8naj7gIAC7uoJApMUjGxRRSFekEU7xQMDEaXAorvpYe0agiyggKUUGdQmbjmalEKcFbvz0hO5cU00OPu0GkWMtMI7ZcwUhk/XORXdORGriWiavKZ4JoSqhIHAdH+p3ggmjsHej+8iJXuA8K9QmS9RkUMj/g5v8fX/Us8DpCnCMEAACrhhUhAAAoxhIvn8BECAAAijEhs4zZ4dAoAAA8727cuHH27NmCAvrMuhAiJSXl5MmTjx8/LtP+5MmThw8flmlMTU09efJkUlKS9CsmQgAAeH6VlpYOHz580KBBf/vb35o2bRoXF1e+z/r160NDQ5ctW9asWbPt27dLjXv27PHz86tbt+7QoUMNO+/duzc4OHjZsmXh4eFffvmlwKHRP5ETBVpROxWZxuXNktVuIzO8jRwiF3rKbYMrT8pU/ZVXU5eMBeQyydkw7WRAqtzSrNwIHycRmeTO/nCf7HzlXArZnpNB5B5LTaYr03IjUTNxsEyMKd2Zew3p11zma0W++6Va+iG5N0jWx4rrzX0iyE+QIp9NuR+rahxlqvg5wh9++OG33367cuWKi4vLW2+9tXDhwm3bthl2yMnJeeedd44ePdq+ffsTJ06MGDFi+PDhDg4OLVu2PHz4cHR09Ndff63vXFJSMmvWrM2bNw8cOPDq1audOnUaO3YsVoQAAPD82rNnz4gRI1xcXIQQEyZM+P7770tL//SV6+eff/b19W3fvr0Qonv37s7Ozr/88osQokmTJmFhYba2f1rvnTt3rri4uH///kKIsLCwkJCQw4cPYyIEAAAlSfEysn4qkJCQ0LBhQ+l2o0aN8vPzU1NTDTskJib6+/vrf/X3909MTOS2lpiY6OfnZ/O/YyNSZxwaBQAA5Zh0Qf29e/fWrVtn2NaqVau2bdsKIfLz8+3t7aVGBwcHIUReXp5hz4KCAn0HqU+ZDk/tbOxEePTo0f3796ekpHh6er766qtdunSR2jMzM5cvX37r1q0WLVrMmzfP2dnZyA0CAED1Y8I5QiFERkZGdHS0YYuLi4s0EXp7e6enp0uNqampKpXK2/tP9c7q1q2blvZH8qbU1FQfHx/ugcp37tixo7ET4d27d8PCwvr163fr1q2+ffv++OOPL774ohBixIgRbm5u06ZN+/LLLydMmLBnzx4jNwgAACBp2bLl2rVryT+1b9/+1KlT7777rhDi1KlTYWFhjo6Ohh3atWs3Y8aM3NxcZ2fnjIyMGzduSDMo6YUXXrh3715ycrKPj09RUdGvv/66ZMkSYyfC119/Xbrx8ssvX7x48dixYy+++OK1a9fOnTv3+PHjGjVqtG/f3sfH5+7du40bNzZym88feQU62cg0MmqUCSrjItPUZFJELjJNTpglNxJu47L6s3GtbLAi0V9tS/e2k5OFkqviK1R0O/dGlFD1ipMT6GLFKqpYseAiJJkBciPhDzcRf7C1o3tzryH5msve3+REKSuyH8oLguX2N5mfTVkfcPlrI4sPGzXhgvqK+0+aNOmTTz5ZsmRJ06ZN58+fv2zZMqm9Z8+ew4cPnz59evPmzbt16zZ27NgpU6asWbNm4MCBjRo1EkLEx8fv2rXrt99+S05OXrFiRUBAwPDhw319fUeNGjV+/PhZs2Zt3749LCysTZs28oJltFrt1atXL1y40LVrVyHExYsXW7duXaNGDSFErVq1mjVrFhUVJe8FAACA6kRl0g+vbt26Z86cefjw4f79+1etWjV+/Hipffjw4WFhYdLt7777rlWrVlu3bu3UqdO//vUvqbGkpESj0fj7+48dO1aj0eTk/P49dd26dREREdu2bQsKCjp48KCQdR3h2rVr586dm5OTM3v27J49ewohUlJSateure/g6emZkkJfRyWEKBPwCgAAlqW0tNSGW4ZXppCQEMNrASXTp0/X33Z1dV24cGGZDoGBgcuXLy+/NUdHx/nz5xu2yHhKr732WnZ29t27d0+cOLFq1SohhLOzc2HhH1cK5+fnVxAsY5aXDwAAlGLMv3GVEPKvnjDzAWHZk1OjRo3GjRsXGRkphKhfv/69e/f0f7p//36DBg0UHBwAAFgWlQnXEZp7zMYeGk1MTJQmucLCwqNHj7Zo0UII0adPn8mTJ58/f75Dhw6RkZEFBQXdunWrxMGaiaz6nEJmijUuMERtS53V5zaipjdiQ7WzYTsy222pEZYysRjc0ydLv3KPKEsRnXlNaOmashW8m0QjV962lKlkK2uX4HLDcchdhStvy5XDdXAk2u0cmJ2T2d/IndmW2cMV2Q/JPbyiEVKvOfkCCv6zqVCKNbP/8680phXmlbnbK8vYibB79+6Ojo4eHh43b94MCwv7xz/+IYRwcXH57LPPBg4cGB4efuXKlf/7v/+TrnYEAACwFMZOhLGxsbGxsRkZGQ0aNKhfv76+ffLkyQMHDoyLiwsKCvLw8KicQQIAgIUw4YJ6lcoyVoRqtTokJIT8U506derUqaPckAAAwFKZch1h5YzEeIjkBAAAq4ak2wAAoBiT6hFW0liMhYnQCLIqhTLhbVwEGhdTp7Yj2slGUUHYG1ndV2YGL+5p2tkTcYlsUCLTTsqmytgKITLTCsh2MoEZXzyWjibl+pNRo2wgMffako/IhJ7qmMwT3BtH7kL2jnTUaA0XO7K9pqdj+UZXdyb2jTmdU5hPhOTyL6y8/ZD+WLHB2Ey7nI8V+9mk2uVmjDP/0UAwgIkQAAAUY8I5QrN/LcA5QgAAsGpYEQIAgGJUQv45QnMvCTERAgCAckzLLGNWODQKAABWDSvCP7Crea7CK3MHMpKNi0CztWeK0NoRUX92THgbGcAphLCl+rOxcGx6RjqK0daO2Hk86tYgOzds6k62O1FBjAmxGWTn2MupZHtJMTFCrjCvlgliLNbR0aTku8/WfWV2IbouLxMdym2cfDcFU2uXzB0qhKjl5US2B7X0LN/oF0S/a/m5xWT7/ZvEG5f2KI/srGOCZrn9kN6ZZX4iyE8Q+VkT/GeT/ARx8atsgDHZ+hxcSPDsFC/MWwUwEQIAgHJMSLFm7mOjmAgBAEA5OEcIAABgWbAiBAAAxVjgghAToTHYnE9MO1n8kzurz6V3ok7Uk2ERFbST1VmLC+moEC4pF1colnxGdeo5k51D29cl28nYDa6o7KPEHLI9g0q9Zk9F0FSAiyEiE6HxwTJcrV1iI2wqNa4SMhO7Qb5c3C5R04NIpSaEaNKCqKEW1tGb7Kx5kk+25+cQQTS52UVkZy3zBnH7oR319LldRdYnhQ9Yk5H+kEuAx/2XMP///kqjkl922OzBMjg0CgAAVg0rQgAAUI4FHhvFRAgAAIoxJem2uWEiBAAAxZhUj9DMMyfOEQIAgFXDitAYSlQQZSqFygoF5ALq2Ng5KtEUF1DH5R7jUmGRT5/bSEEeUbJVCJFHxRkWFdBxrUz4KvOCMyF/Dsx3Tzt7OoiRjO3kvr9yVVjJ15B7OlxIqg0TgGhD7Vpc6CkZvyqY15x8dwT/bpJp7fhwSpk7Mxkcy6RS4zZCfoK4zlz+NvKzzL3g/FFCSzt6aDycIwQAAOumkltWyexlmHBoFAAArBpWhAAAoBjVc3CBvFyYCAEAQDk4R1gtya0cRkYHcBm8uFpo5Al8rsicvQMdvEDnaWPiC0qZOBcyx5hgUrWlpdDF566dTyHb7alnxKVS45J1kWy4OnDMLq8W9MtCBt04OtFb4QKRyNeqIJ8JOZGZHE5W6HluFh3/EncjvXxjTib9ghcV0iOn330mJkhWXkDB7LeyYs0E8wmS9YiC+SxzMUFsXBXdXB3g8gkAAAALgxUhAAAoxpQK9U/rcPfu3a+//jojI2PIkCEDBgwo3yE/P//LL7+8efNmaGjoG2+84eDgILVHR0dv3ry5tLT01Vdfbd++vRDi4cOH27ZtM7zvkCFDsCIEAADlqEz64aWlpXXs2FGtVnfq1GnKlCm7d+8u32fs2LGRkZG9evU6ePDg1KlTpcbr16/36NHD398/ODi4d+/e0dHRZe6VkpLy3nvv2draYkUIAADPr82bN7dq1Wr58uVCCBsbm5UrV44YMcKww61bt3788cdHjx65ublFRETUr19/6dKlfn5+//znPydNmvTOO+8IIZKTkz/77LNt27bVq1dv/vz50h2XL1/erVu3gIAArAgBAEAxKvkqXhGeO3euR48e0u0ePXpER0cXFhYadjh//nyrVq3c3NyEEJ6ens2aNfv111/L3/HcuXNltrxly5YpU6YInCM0isy4LybFmowCvILJBeXApFhzcJIRI1dSRAclctGhbDsVZZqfSwcl3o/NMH4jXAYvrqQwUwmZy3bG5Lpj3iCykm39Jm5kZ3dPos6wECIjlahk++BOFtk5k6ozLIQoKaHfODIlGxfEWFxEv4YpCdnlG8lhV7BxMkCU2/Pl1tSVFUfNfyJkZC7kPpvkM+JSrCHDmpF3qUBycrKnp6d028vLS6fTpaSkNGzYkOwg9UlOThZCpKSkGN5RatQ7depUcnLysGHDBCZCAAAwu7Nnz0ZERBi2DBgwYM6cOUIIR0fHoqLfr+SRbjg6/unrqaOjY3HxH1/BCwsLpQ4ODg6Gdyxzr40bN44ePbpGjRoCEyEAAChM/no3MDBQf+pO0rRpU+lGvXr1Hjx4IN1OSEhwcHDw8vIy7FmvXr3ExET9r4mJifXr1y/TnpiYWK9ePX2fnJycvXv3Hj9+XPoV5wgBAEA5JpwkVKnq1KnT688aNGggbW/o0KF79+4tKCgQQuzYsWPQoEFSJZaTJ0/GxcUJIXr37n3//v3Lly8LIc6fP5+ent69e3fpjjt27JAqrmzfvn3o0KH6Me7YscPf379t27bSr1gRAgCAYky5jrDC/oMGDVq/fn3btm0bNmx46dKlyMhIqX3evHmjR4+eM2dOzZo1ly5d2rdv3y5dupw+fXrlypXSAc/p06fv3LmzS5cudnZ2jx8/Xrt2rX6bGzdulMJkJJgIAQDg+WVra3v48OGoqKj09PTOnTu7urpK7bt379bfnjVr1oABA27duvXpp5/6+/tLje7u7lFRUefOnSstLe3UqZP+KnshxI4dO3x9ff94iCp6KtWTjDgxLtCOrSxKRo3WoN8vxwI7sp2stlpSTMcNarVMNKmWyU1KRZMW5tMhpoX59IOSgaAlTLJN7jUkQwHVTGeuHC4XOujT0LV8Y6sXfcs3CiH8gtzJ9gQqaFZbzCRxZd4grl4xub9xRWW515CMJS7Q0eG7XGAnGavJvbDcnm/HxGqSG3dyZj4RNehPBPkJkjtC8jXks2VW3/BQjglho09jY2PTrl27Mo1+fn6GvwYEBAQEBJTpY2dn161bt/IbbNy4seGvmAgBAEAxKpX8JNpIug0AAGBGWBECAIBiTCjMa+4FISZCAABQjko8LWcacRczw0RoOllnx7kMTFxQg70D8dY4OtFRJMUudHsJlU+LC0XRyizMq9MRoQQ6mXnauP4kMpeYYF5bO3t5I+EiJmq4EmEXbrWJvGtCCHcqH5sQIoPqT25ZCOHIhERxuxDZzgXFcPubWk20c53tHGRkTeNTqXFBMfTTd3ImXq4aLvRr6MS0kxWVyc+aqCDgiHrB2QK8Zv8fX/UssEI9zhECAIBVw4oQAACUI/+CerPDRAgAAIr5vbKSzLvIOEdSCYydCNPT048dO5aQkODn5zd48GB9Gm+dTnfo0KGYmJjw8PC+fftW2jgBAAAqhbHnCFu3bv3tt9+mpaV99dVXLVu2zMj4PVPGzJkz//73v+fl5c2ZM+dvf/tbpY0TAAAsgcqkH7MydkV48eJFqcKhVqsNDQ39/vvvJ0+e/PDhww0bNty5c6devXrjx48PCwt75513PDw8KnPAloA6LGDDfOVQ2zKp16iYOgct/X5pS+jjClqqlKuWjQKlR8iRnTzC6I2omeKxJcVcoVTqaTJlbLnXyoZ5OrlZReUb4/+bTnbOSqdr6qal5Bm55QpGwmX8InchrhwuGR0qhLCl6hjbMo/IF4gm9k8uCpQLjnV2tafbaxLtNdzozjWoEFPBpFjj4le5zyYdvmtxp8UqjSlJt+mizlXH2BWhvs6vWq22t7dXq9VCiF9++aVFixZSkacmTZr4+/ufPn26kgYKAADPP9XvU6GMH7MvCWVfPrFr167U1NTBgwcLIZKSknx8fPR/8vb2TkpK4u6ok7viAACA50l1/TcuL2r05MmTM2bMOHDggLu7u5BCfQxeF51Op8jhMgAAsGCWdkG9jInw3LlzI0eO3LlzZ8eOHaUWHx+f5ORkfYeUlBTDCk9lYI4EALBoxvwbV7wwbxUw9tDopUuXhg0btnHjxpdeeknf+NJLL12/fv3BgwdCiNjY2ISEhK5du1bKMAEAACqHsSvC/v37Ozs7b9myZcuWLUKIoUOHjh492tfX97XXXuvdu/fQoUN37dr1zjvv1KpVqzJHaxnobzdcnkwdlxOSTNvI5M905vJqEpGT3EF+NnMqk+KSrHzLlcPlUl+SQYxFTMRjMRdNWkJsXFtMP6JWLe8kR15OcfnGu0zU6IO4TLK9iKo/TG5ZVBDwyXxY6ahRJk+mLbNxMiSVSxNqT9XIFUwmTy461MmFDvh0caMDPp1rOhCNTIipI1Owl3xG5GdN8MldyZy3Zl/TPD9MuaC+koZiNGMnwq+++kqr/ePDHBISIt34/PPPf/755xs3bqxdu9ZwsQgAANbIApNuGzsRDhs2jPtTnz59+vTpo9B4AADAkpmQa9TcEyGqTwAAgFVD0m0AAFCMSv41AnIL+SoOE2GVkXHiXQihps7fc0m2dHQ2MbrsbeUGy7CBHlzsRgnRmQv0KKLbS4qI58/VHy7V0u2yQmgK8+mwHa6dxAUQcYEe3BtkI6emrq29jGAZWanUBBMXQxbUFXxdYi7+hezvxLxW3AjJp8lFeHF7vvkP5IHScGgUAACsGlaEAACgGNPqEVbSYIyEiRAAABRjQmYZsx9sxkQIAADKMX8xCdlwjhAAAKwaVoRVhD9WwOUkIxptuS8uTMijTk4sJHeYno9rpUq5MsGKXLyrHRXEaOdAx16SicqEECVU6jUualSrZeoSK1Sv+Nlxu4qs8F0+apTJmkblHuOiRh1r0AGfTNQol2KNiSZlUq+R2yEL7Qoh7JmnSb4sXPgu+4mwtOVOFTPpHGEljcVYmAgBAEA51bj6BAAAQLWEiRAAAKwaDo0CAIBiVMKEMky4jtC6yQqisWEW8FxkhFAR7y8fFCOvlCBZxc3Onq6xZ8dUtiPbHfKJvGtCiEImWKaYameDZUq4dq6mY1VHy3Bl8Mi6g4JJX8ftEmyJQSouxt5RRio1IYQT1e7IpFjjNsK1k1nTuKfDlV0kX1sExSirMq4jXL169eLFi/Pz8/v167dx40YXF5cyHS5cuDBlypTbt2+HhIRs2rSpZcuWQgidTjdv3rxvvvlGCDFhwoTPPvvMxsZGCNG6deuMjAzpjv3791+9ejUOjQIAwPPr0qVLCxcuPHHixOPHj7OzsxcvXlymg1arHTVq1IwZM/Ly8iZMmDB69GidTieE+O677/bv33/r1q24uLh///vfmzdvlvrfv3//m2++OXbs2LFjx95//32Bc4QAAKAklUk/vM2bN48YMaJZs2ZOTk7z58/ftGlTmQ6//PJLYWHhtGnTbGxsZsyY8eTJk/PnzwshNm3aNH369Lp163p4eMycOdPwjg0aNGjcuHHjxo3r1q0rMBECAICCVCapYIOxsbGhoaHS7dDQ0CdPnugPbOo7NG/eXNqIra1t06ZNY2NjhRC3b982vOPt27f1d+nVq1fDhg1HjBgRFxcncI4QAADMLjc39+7du4Ytnp6ebm5uQgiNRqM/Kejq6iqESE9Pd3d31/fMyMhwdnbW/+rm5paenl7+jlKjEGLLli0vvPBCQUHBkiVL+vXrd+3aNUyEAACgGBOCZVRCnDt3LiIiwrBx7NixixYtEkJ4enpmZmZKjdJa0MvLy7Cnh4dHVlaW/teMjAypg4eHh+Ed9fcaMGCAdGPt2rW1atW6fPkyJsLnFL0ncYGdzG5nS6VH4/ZRPliRyZpGxSWSKdMEE5QomFDAIi5qtIBuJ1OvFVPVekUF0aRyokwViSTl3jW2uDGXNU3WG8GEWTpQAaL2XAFeJxnvpiOzEW7j3AjJZ8S9VrKypiE6VHnyo0YjIiL27dtH/jE4OPjq1avS7atXr9arV09aFxp2uH79ularVavVRUVFMTExwcHBQoiQkJArV67069dPuqPUSDy4issjCQAA8ByYPHny3r17z5w5k5qaumTJkqlTp0rt77777qFDh4QQXbt2rV279ieffJKdnb1s2bJGjRq1adNGCDF16tSvvvrq1q1b8fHxX3zxhXTHmzdvHj58ODk5+d69e6+//rq3t3d4eDhWhAAAoCRlL5Bv3rz5119/PW3atIyMjKFDh7733ntSu0ajycvLE0KoVKrvv//+zTffXLVqVXh4+O7du6UOgwcPvnnzZr9+/bRa7dSpU//yl78IIYqLi1esWBEXF+fg4NChQ4effvrJwcFBpauS7Pr5+flvLlo5cOK0Knisaox7r7g3sZSqtMBdUc4dS+RqPpAHMAu4o5pMe0E+sXEcGqXbcWi0HBwarWLDmvo8tc+Gn2IfZxTI2ux/o47nxJ3kDo1WAawIAQBAMaaUYTJ3ijWcIwQAAKuGFaElUaS6LxchJavuqxDClkp9ydZ9ZaNGqUOjBfRuyR4aLVDi0ChzxLiEPDTKVfdljlGTX5DZJK7coVGmXd6hUe6NIA+NMp3ZdgcyHai8Ws3cEWByP+R2WqQPNaenZYqh72JWmAgBAEAxlZF0u7Lh0CgAAFg1rAgBAEAxlliPECtCAACwalgRVgeygmhUdIwCG1/ApV4jYz3UdnTICX9ZG9G/mIqgEUI4FtK7a1ERVZhXoWAZ+lpM7kJCrp16DdUyKyHLCpax5V5wJkSFrJDMxbNw7yb5oHJjf7j9jYvxolvNfc7JmpmSa9Tc7xdWhAAAYNWwIgQAAOXIv6De3KcIMRECAICCTLiO0NxwaBQAAKwaVoQAAKAYSwyWwURYnXF5pujezNEBrkgCmapNbtowO3sizLKkmA5WLHHkAkGpqFE2OpQO7NRq6f46KhCUrT7B1QchU6xxLyyb644Lv6Ry3bGlKujXlg495bKdUY8omHefjwJFdrTqyZSk2+Z+dzERAgCAoiztWwvOEQIAgFXDihAAABRjgUGjmAgBAEA5lniOEIdGAQDAqmFFCL/jvpNxieF1ZNVfLn+mDR1OqaY2wgUrcoGd2hJiK6VcZ6ambikTCaqjNsMV4JWFDY/kwnflVE62kRNiKpiQVJlpP+l3X+5+ZXmH1aAMFOYFAACrppJ9qNPsh0YxEQIAgGJMuaC+ckZiPHnnCFNSUlJSUso0Jicnnzhx4tGjR8qNCgAAoIoYOxGuW7fOx8fH19d36tSpZdpbtGixfPny5s2b79ixoxJGCAAAFkUl/8esjD002qVLl+PHjx8+fPjkyZP6xuzs7Llz50ZGRrZr1+748eOjRo165ZVXHBwcKmeoYCZygx3IzkzqNUEF0djomIgbJtDD1o7YiK6UziXGBsUw8S/PT7CMrMrJsuJZhBA21MbZA1xyjnyZ++wPVDWVUMn6lVqbywAADmFJREFU5/A8MHZF2KxZs5CQEBubP/X/+eefGzRo0K5dOyHESy+9VKNGjRMnTig+RAAAgMrzTMEyiYmJ/v7++l8bNmyYkJDwrCMCAACLZXXVJwoKCuzt7fW/Ojg45Ofnc51LS+nrugAAwCKUlpaWOS5IeA7O+cn1TBOht7d3Wlqa/tfU1FQfHx+u89NfPgAAeI4Z82/c6lKstWvXLjo6Ojc3VwiRnp4eExPTtm1bhQYGAABQFYxdEd64cePw4cOnT5+Oi4tbsWJFWFhYv379mjdv3rVr17Fjx06ePPnrr78eNGiQ4SlDsEKyv9iRd2BCMukwUCF0ZCQksxE23JOLGqWa+aBR7g9yco/JjNWUVX2Zf4PkRI0C8Ew4R2j2Q6nGrgiLioo0Gk1oaOjgwYM1Go20ChRC7Nq1q1WrVtu2bevateuWLVsqbZwAAGAJLO0iQmH8irBVq1atWrUq3+7q6rpw4UJFhwQAAPCHoqKiY8eOZWRk9OzZ09vbm+xz8eLFmJiYFi1aGE5V2dnZR48eLS0t7d27d82aNaXGx48f/+c//8nLy2vTpk1gYKBAGSYAAFCUygQVbK6oqKh79+4ff/zx0aNHmzdvfunSpfJ9Fi1aNHz48LNnz7788suffvqp1Pj48ePw8PDNmzfv3LmzRYsWSUlJQohTp041bdp0/fr1R44cadu27apVqwSSbgMAgIIUr8K0Z8+enJyc6OhoOzu7xYsXf/TRRwcOHDDskJaWtnz58suXLwcFBV29erVLly7Tpk1zc3P76quvwsPD9+3bJ4QYP378559/vnLlypCQkHv37kmrw2PHjg0aNGjmzJlYEcJzRzrZTv3IYUP/2HA/avpHTf7Ycj82zA/Vmdk4NxJu5PQzlf+NvPwPgClMOEdY4c52+PDhIUOG2NnZCSGGDx/+ww8/aLVaww7Hjh0LCgoKCgoSQoSFhfn4+Eg5zg4dOjR8+HCpzyuvvHLo0CEhhJeXl/4Yqa+vb0lJSUlJCSZCAAB4fj18+LBevXrS7fr165eUlJQpdmTYQQhRr169hw8flr/jgwcPymx58eLFY8aMcXJywqFRAABQzFPP+ZFu3LixYMECw5auXbv2799fCFFSUqJW/37xlHSjpKTEsKdWqzW80t/W1lbqoNVqDe9YZh25aNGi//73v1IZCUyEAACgGNNyjdrZ2dWqVcuwUZ+/09fX9/Hjx9LtR48eqVSqMoGjPj4+T5480f/66NEjX19fqd3wjlKj5JNPPtmxY8eJEyekB8VECAAAZhYUFDR//nzyT926dduzZ8/f/vY3IcTRo0c7d+4szZFZWVkODg4ODg4vvvjiX//619TUVE9Pz6SkpNjY2M6dOwshunfvfvTo0VdeeUUIcezYse7du0sb/OKLL9atW3fixAn9hIqJEAAAFCPFtcm8S0X9x48f/8knn0ydOjUoKGjFihVbt26V2iMiIkaPHj1nzhx/f/9Ro0a9/PLLY8aM2bJly8SJE6UZbubMme3atatdu7adnd2GDRvOnj0rhPjll1/eeuutvn37Llq0SNrOokWLMBGCxVMovpEtQGx0VwAQQij8CalZs+bFixc3b96clZX1ww8/tG/fXmp/7733AgICpNsbNmzYvn17TEzMW2+9NXr0aKkxMDAwKipq+/btOp3u119/bdq0qRCiYcOGa9asMdy+g4ODSpFa20+Vn5//5qKVAydOq4LHAgCAyjCsKVtfSO/76MS03CJZm718+lhS1M/SBX9mgRUhAAAoyeIOmWAiBAAAxZhw+YQJl1soCxMhAAAoR/Eca5UPmWUAAMCqYUUIAACKMe2CevPCRAgAAIoxLcWaeeHQKAAAWDWsCAEAQDkWGCyDiRAAABSjkn/Oz+xHUjERAgCAYqS6znLvUkmDMRLOEQIAgFXDihAAAJRjwjlCc8NECAAAypF/HaHZJ04cGgUAAKuGFSEAACgGSbcBAMC64RwhAABYM5UwYUVYSWMxFs4RAgCAVcOKEAAAFGOBGdYwEQIAgIIscCbEoVEAALBqWBECAIBiVCoTLofA5RMAAFBtoEI9AABYNZwjBAAAsCxYEQIAgGJMuaDe3EtCTIQAAKAYFapPAAAAWBZMhAAAYNVwaBQAABRjUhmmShqLsTARAgCAYkw4R2jueRCHRgEAwLphRQgAAMqxwAvqMRECAIBy5J8jNPtJQhwaBQAAq4YVIQAAKEYlf4Fn7gUhJkIAAFCOKZdPmPskYRVNhLa2tjdP/nxy17+q5uGqWFpaWq1atWxsqvlxZjzN6iQ1NdXDw0N+3TgLg6eprJDDh0NCQiru09PfU+5mG+e3+o9NgamDUoBKp9NVzSOlp6dnZGRUzWMBAIDi6tevb29vb+5RKK/qJkIAAIDnUDU//gMAAFAxTIQAAGDVMBECAIBVw0QIAABWDROhbMuWLevTp0+TJk2OHDli2L58+XJPT8+aNWtOmjSpsLDQXMNTRGFh4cSJE+vXr29nZ9eiRYvDhw/r/7Rv376GDRs6Ozv36tUrKSnJjINUxIQJE7y8vNRqdePGjdetW6dv37Vrl5+fn4uLS9++fR89emTGESooOzu7Xbt2I0aM0LecOXOmWbNmTk5O7du3j4mJMePYFDFy5MiA/xkwYIC+ffv27fXr13dxcRk4cGBqaqoZR6iUL774omHDhmq1OjAwMDY2VmrU77R9+vSpNjttFdGBTMuWLdu5c2dgYODOnTv1jcePH69bt25cXFx2dnbnzp2XLVtmxhE+u+zs7MWLF8fHx5eUlGzdutXJyenBgwc6ne7x48cuLi6RkZHFxcWvvfbasGHDzD3SZ3Xq1CmNRlNaWnr8+HEnJ6dLly7pdLqHDx+6uLicPHmyqKho0qRJY8aMMfcwlTF9+vROnTp16tRJ+rWwsNDb23vLli1arXbJkiWtW7c27/CeXefOnb/55pu4uLi4uLjExESp8d69ey4uLv/5z38KCwvHjBkzadIk8w7y2a1du7ZJkybR0dFarTYmJubJkyc6nS4pKUm/006ePHn06NHmHqYlwURoovDwcMOJcNy4cfPmzZNu79u3LzAw0EzjqhR16tQ5duyYTqf74osvevToITUmJCTY2tqmpqaadWhKCgoK2rdvn06nW7lyZb9+/aTGO3fu2NvbZ2ZmmnVoCvjll1+6d+++bt06/US4f//+gIAA6XZ+fr6Li8uVK1fMN0AFdO7c+eDBg2UaFy1aNGTIEOn2jRs3nJyc8vLyqnxoSvLz89u/f3+ZxpUrV/bt21e6HRcXVz122iqDQ6PKuH37dmhoqHQ7NDT07t27Wq3WvENSSnR0dG5ubnh4uBDi9u3bLVq0kNobNGjg7Ox89+5ds45OARcvXty5c+eMGTO8vLz69OkjhIiNjdU/zYCAALVaff/+fbOO8Vnl5eXNmDHj66+/Nkw+YrjTOjo6BgQE6A+yWa4333zT19e3X79+UVFRUovhThscHFxUVJSYmGi+AT6rjIyMhISExMTEsLCw4ODgpUuXlpaWij8/zcaNG6vV6nv37plzoBYFE6EyNBqNi4uLdNvV1VWr1WZmZpp3SIpITU0dPXr0ihUrvLy8hBAajcbZ2Vn/Vzc3t/T0dPONThkXL17cs2fPTz/91Lp1a7VaLco9TVdXV0t/mn//+9/HjRsXHBxs2Fj93s0PP/zw9OnTFy5c6NChQ0REhHQO2/Bp2tjYuLi4WPTTTE5OFkIcPnz4+PHjBw4c2Lhx44YNG0R1fDerEiZCZXh4eOhnvoyMDFtbW3d3d/MO6dlpNJqIiIgRI0a8+eabUouHh0dWVpa+Q0ZGhjRBWrQ33nhjz549N27c+OWXXzZv3iyE8PT0NHyamZmZFv00r127tmfPnu7du0dHR9+/fz8nJyc6OlpUx3ezV69eDRs2rF+//sKFCxs2bBgZGSn+/DS1Wm12drZFP01PT08hxLx58zw9PZs2bTp16lQplq3Mu2npO20Vw0SojODg4KtXr0q3r169GhQUZOlZmzMzM/v27du9e/elS5fqG4ODg69cuSLdjo+PLywsbNSokZkGqDAHB4fQ0FDpoJnh07x586aNjY2fn59ZR/dM8vPzg4OD//GPfyxYsGDv3r33799fsGCBTqeTdlqdTieEyMvLu3PnTtOmTc09WMWoVL/njzR8N69fv16jRo169eqZdWjPxNPTk5zhDJ/mrVu3VCpVw4YNq3ZolszM5ygtUGxsbFRUVGBg4McffxwVFZWVlaXT6c6ePevh4REdHZ2cnNy6detVq1aZe5jPJD8/v2PHjj179rx48WJUVFRUVJQUFJOenu7m5rZ79+7MzMyxY8daejhlWlra1q1bExMTHz16tGvXLldX1zNnzuh0usePH7u6uh44cCAjI2PEiBGTJ08290gV88033+iDZYqLi/38/L744ovs7Ox58+Z17tzZvGN7Rqmpqdu3b3/w4EFSUtLixYtr1qz58OFD3f9igH/44QeNRjN48ODp06ebe6TP6v333+/bt296evqdO3cCAgI2bNig+99Ou3///oyMjJEjR1aD4NiqhIlQtjfffLO1gaioKKl9zZo1TZo08fX1ffvtt4uLi807yGeUlJTU+s+OHDki/SkyMrJVq1ZeXl4jR4609JBRjUYzZMgQPz8/b2/vLl26SCGjkh9//DE8PLxOnTqjR4/WaDRmHKSy9u3bZ/gv8tKlS507d/bw8Ojbt298fLz5xqWAJ0+eRERE+Pr6SsEyv/76q/5PBw8ebNGiRZ06dcaPH18NYikLCgqkmKDg4OCVK1eWlpZK7T/99FO13GmrAKpPAACAVbPs81gAAADPCBMhAABYNUyEAABg1TARAgCAVcNECAAAVg0TIQAAWDVMhAAAYNUwEQIAgFXDRAgAAFYNEyEAAFg1TIQAAGDVMBECAIBV+/8aSFLLGkGZVQAAAABJRU5ErkJggg==",
"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.ρ[:, :, 1, 1], c=:blues)"
],
"metadata": {},
"execution_count": 5
}
],
"nbformat_minor": 3,
"metadata": {
"language_info": {
"file_extension": ".jl",
"mimetype": "application/julia",
"name": "julia",
"version": "1.6.2"
},
"kernelspec": {
"name": "julia-1.6",
"display_name": "Julia 1.6.2",
"language": "julia"
}
},
"nbformat": 4
}