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