{ "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", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\n", "\n", "\n" ], "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\n", "\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 }