{ "cells": [ { "cell_type": "markdown", "source": [ "# Gross-Pitaevskii equation with external magnetic field" ], "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.105161e+01 8.135250e+00\n", " * time: 0.0018990039825439453\n", " 1 2.314313e+01 4.185435e+00\n", " * time: 0.008825063705444336\n", " 2 1.648258e+01 3.595554e+00\n", " * time: 0.015731096267700195\n", " 3 1.273536e+01 2.382911e+00\n", " * time: 0.021572113037109375\n", " 4 9.961725e+00 1.291876e+00\n", " * time: 0.026746034622192383\n", " 5 9.223828e+00 1.860277e+00\n", " * time: 0.031507015228271484\n", " 6 8.741099e+00 1.175847e+00\n", " * time: 0.03589200973510742\n", " 7 8.472479e+00 7.961030e-01\n", " * time: 0.04087519645690918\n", " 8 8.311689e+00 9.525920e-01\n", " * time: 0.09021902084350586\n", " 9 8.275526e+00 8.296795e-01\n", " * time: 0.0932002067565918\n", " 10 8.246208e+00 7.421042e-01\n", " * time: 0.09605216979980469\n", " 11 8.240311e+00 1.001146e+00\n", " * time: 0.0989840030670166\n", " 12 8.150341e+00 4.959490e-01\n", " * time: 0.10269021987915039\n", " 13 8.126966e+00 7.912038e-01\n", " * time: 0.10555601119995117\n", " 14 8.084134e+00 6.355446e-01\n", " * time: 0.10855603218078613\n", " 15 8.039537e+00 5.498641e-01\n", " * time: 0.11141109466552734\n", " 16 7.980361e+00 5.325098e-01\n", " * time: 0.1150510311126709\n", " 17 7.963144e+00 4.207574e-01\n", " * time: 0.11781501770019531\n", " 18 7.924472e+00 4.457566e-01\n", " * time: 0.12082314491271973\n", " 19 7.893223e+00 3.922636e-01\n", " * time: 0.12359809875488281\n", " 20 7.869237e+00 4.798811e-01\n", " * time: 0.12682509422302246\n", " 21 7.841623e+00 3.568805e-01\n", " * time: 0.1300210952758789\n", " 22 7.818467e+00 3.179261e-01\n", " * time: 0.13436317443847656\n", " 23 7.802398e+00 2.418383e-01\n", " * time: 0.13863611221313477\n", " 24 7.783704e+00 3.603037e-01\n", " * time: 0.14192700386047363\n", " 25 7.777750e+00 3.137056e-01\n", " * time: 0.145064115524292\n", " 26 7.770279e+00 2.679355e-01\n", " * time: 0.14835119247436523\n", " 27 7.759338e+00 2.958856e-01\n", " * time: 0.1517200469970703\n", " 28 7.742685e+00 1.641332e-01\n", " * time: 0.15586400032043457\n", " 29 7.741704e+00 4.026390e-01\n", " * time: 0.15929722785949707\n", " 30 7.735567e+00 2.870783e-01\n", " * time: 0.16289806365966797\n", " 31 7.727467e+00 2.017193e-01\n", " * time: 0.1666111946105957\n", " 32 7.715692e+00 1.545526e-01\n", " * time: 0.170166015625\n", " 33 7.707414e+00 1.633176e-01\n", " * time: 0.17330622673034668\n", " 34 7.701193e+00 1.298772e-01\n", " * time: 0.17709612846374512\n", " 35 7.692984e+00 1.095916e-01\n", " * time: 0.18113422393798828\n", " 36 7.686663e+00 9.252611e-02\n", " * time: 0.18497705459594727\n", " 37 7.682248e+00 7.543221e-02\n", " * time: 0.18881607055664062\n", " 38 7.681415e+00 1.741614e-01\n", " * time: 0.19238615036010742\n", " 39 7.678852e+00 1.164613e-01\n", " * time: 0.1958920955657959\n", " 40 7.674706e+00 9.776009e-02\n", " * time: 0.22667908668518066\n", " 41 7.673233e+00 1.114376e-01\n", " * time: 0.22996902465820312\n", " 42 7.672548e+00 1.340283e-01\n", " * time: 0.2328171730041504\n", " 43 7.669825e+00 1.073752e-01\n", " * time: 0.2353651523590088\n", " 44 7.667984e+00 9.663599e-02\n", " * time: 0.23780608177185059\n", " 45 7.667408e+00 1.086421e-01\n", " * time: 0.24035120010375977\n", " 46 7.666101e+00 1.145175e-01\n", " * time: 0.24303102493286133\n", " 47 7.664064e+00 8.339079e-02\n", " * time: 0.24574804306030273\n", " 48 7.662725e+00 1.231989e-01\n", " * time: 0.2483510971069336\n", " 49 7.659795e+00 9.562690e-02\n", " * time: 0.2526891231536865\n", " 50 7.659078e+00 1.404238e-01\n", " * time: 0.2553441524505615\n", " 51 7.656099e+00 1.103888e-01\n", " * time: 0.25800609588623047\n", " 52 7.653176e+00 8.668524e-02\n", " * time: 0.2617311477661133\n", " 53 7.650511e+00 7.947619e-02\n", " * time: 0.2653062343597412\n", " 54 7.648505e+00 8.803798e-02\n", " * time: 0.26812314987182617\n", " 55 7.647721e+00 1.558565e-01\n", " * time: 0.2708871364593506\n", " 56 7.645242e+00 1.237427e-01\n", " * time: 0.27442407608032227\n", " 57 7.642671e+00 8.562518e-02\n", " * time: 0.27784013748168945\n", " 58 7.639461e+00 1.496980e-01\n", " * time: 0.2806510925292969\n", " 59 7.637211e+00 1.154253e-01\n", " * time: 0.28344106674194336\n", " 60 7.634552e+00 8.822295e-02\n", " * time: 0.28612303733825684\n", " 61 7.632310e+00 1.140565e-01\n", " * time: 0.28902101516723633\n", " 62 7.630435e+00 7.065791e-02\n", " * time: 0.29329919815063477\n", " 63 7.628798e+00 9.121952e-02\n", " * time: 0.29630112648010254\n", " 64 7.627685e+00 1.170035e-01\n", " * time: 0.29877209663391113\n", " 65 7.626295e+00 7.804397e-02\n", " * time: 0.3016390800476074\n", " 66 7.624865e+00 4.720245e-02\n", " * time: 0.30547118186950684\n", " 67 7.622896e+00 6.227335e-02\n", " * time: 0.3084390163421631\n", " 68 7.622133e+00 1.004908e-01\n", " * time: 0.31171417236328125\n", " 69 7.620874e+00 8.568546e-02\n", " * time: 0.31481218338012695\n", " 70 7.619603e+00 5.406312e-02\n", " * time: 0.31878113746643066\n", " 71 7.618457e+00 7.451042e-02\n", " * time: 0.32225823402404785\n", " 72 7.617883e+00 8.275176e-02\n", " * time: 0.32557201385498047\n", " 73 7.617249e+00 4.917774e-02\n", " * time: 0.3506031036376953\n", " 74 7.616998e+00 7.541408e-02\n", " * time: 0.35337114334106445\n", " 75 7.616222e+00 5.886760e-02\n", " * time: 0.3566911220550537\n", " 76 7.616076e+00 5.785629e-02\n", " * time: 0.3590400218963623\n", " 77 7.615873e+00 6.288340e-02\n", " * time: 0.3616490364074707\n", " 78 7.615467e+00 3.940573e-02\n", " * time: 0.3644142150878906\n", " 79 7.615114e+00 4.007669e-02\n", " * time: 0.36685705184936523\n", " 80 7.614864e+00 3.507605e-02\n", " * time: 0.3695840835571289\n", " 81 7.614590e+00 1.946838e-02\n", " * time: 0.37322402000427246\n", " 82 7.614205e+00 3.036793e-02\n", " * time: 0.3761141300201416\n", " 83 7.614005e+00 3.361125e-02\n", " * time: 0.3790090084075928\n", " 84 7.613957e+00 3.774864e-02\n", " * time: 0.3817331790924072\n", " 85 7.613711e+00 3.637523e-02\n", " * time: 0.38542604446411133\n", " 86 7.613461e+00 3.933052e-02\n", " * time: 0.388214111328125\n", " 87 7.613056e+00 3.161323e-02\n", " * time: 0.3915441036224365\n", " 88 7.612684e+00 2.196268e-02\n", " * time: 0.3952171802520752\n", " 89 7.612494e+00 3.286099e-02\n", " * time: 0.39803123474121094\n", " 90 7.612237e+00 3.605917e-02\n", " * time: 0.4013381004333496\n", " 91 7.611882e+00 1.669469e-02\n", " * time: 0.40493321418762207\n", " 92 7.611645e+00 2.610068e-02\n", " * time: 0.4072701930999756\n", " 93 7.611348e+00 2.468908e-02\n", " * time: 0.4098222255706787\n", " 94 7.611160e+00 5.301449e-02\n", " * time: 0.4126100540161133\n", " 95 7.610883e+00 3.554461e-02\n", " * time: 0.4160802364349365\n", " 96 7.610636e+00 3.287128e-02\n", " * time: 0.4184250831604004\n", " 97 7.610599e+00 5.401654e-02\n", " * time: 0.4208970069885254\n", " 98 7.610475e+00 4.326270e-02\n", " * time: 0.4235551357269287\n", " 99 7.610326e+00 3.837223e-02\n", " * time: 0.426530122756958\n", " 100 7.610199e+00 4.102934e-02\n", " * time: 0.4295802116394043\n", " 101 7.610146e+00 4.947994e-02\n", " * time: 0.432833194732666\n", " 102 7.609949e+00 2.699965e-02\n", " * time: 0.4366450309753418\n", " 103 7.609878e+00 5.601682e-02\n", " * time: 0.4394512176513672\n", " 104 7.609659e+00 3.411716e-02\n", " * time: 0.44228410720825195\n", " 105 7.609640e+00 3.025929e-02\n", " * time: 0.45594000816345215\n", " 106 7.609444e+00 2.560685e-02\n", " * time: 0.4584650993347168\n", " 107 7.609402e+00 3.994197e-02\n", " * time: 0.46112608909606934\n", " 108 7.609225e+00 2.016656e-02\n", " * time: 0.4647231101989746\n", " 109 7.609103e+00 1.844439e-02\n", " * time: 0.4679450988769531\n", " 110 7.609102e+00 3.893100e-02\n", " * time: 0.4706861972808838\n", " 111 7.609098e+00 3.559330e-02\n", " * time: 0.4739551544189453\n", " 112 7.609040e+00 1.972980e-02\n", " * time: 0.4774460792541504\n", " 113 7.608979e+00 2.505345e-02\n", " * time: 0.48082518577575684\n", " 114 7.608854e+00 2.358798e-02\n", " * time: 0.4842720031738281\n", " 115 7.608724e+00 1.913385e-02\n", " * time: 0.4877440929412842\n", " 116 7.608677e+00 3.886704e-02\n", " * time: 0.4902992248535156\n", " 117 7.608587e+00 3.548366e-02\n", " * time: 0.49326515197753906\n", " 118 7.608576e+00 2.523140e-02\n", " * time: 0.4960181713104248\n", " 119 7.608479e+00 2.157015e-02\n", " * time: 0.4984281063079834\n", " 120 7.608390e+00 2.717405e-02\n", " * time: 0.5008001327514648\n", " 121 7.608247e+00 2.251121e-02\n", " * time: 0.5044100284576416\n", " 122 7.608109e+00 1.950110e-02\n", " * time: 0.5079312324523926\n", " 123 7.607974e+00 1.798669e-02\n", " * time: 0.5114550590515137\n", " 124 7.607944e+00 3.301787e-02\n", " * time: 0.5146310329437256\n", " 125 7.607922e+00 2.950096e-02\n", " * time: 0.5173652172088623\n", " 126 7.607817e+00 2.240230e-02\n", " * time: 0.5198860168457031\n", " 127 7.607629e+00 2.142822e-02\n", " * time: 0.5232431888580322\n", " 128 7.607541e+00 2.595930e-02\n", " * time: 0.5258312225341797\n", " 129 7.607372e+00 1.950580e-02\n", " * time: 0.5293631553649902\n", " 130 7.607354e+00 3.889747e-02\n", " * time: 0.5318241119384766\n", " 131 7.607168e+00 2.416859e-02\n", " * time: 0.5356321334838867\n", " 132 7.607047e+00 2.548437e-02\n", " * time: 0.5386090278625488\n", " 133 7.606951e+00 3.757998e-02\n", " * time: 0.541600227355957\n", " 134 7.606857e+00 2.454002e-02\n", " * time: 0.5448851585388184\n", " 135 7.606706e+00 3.017237e-02\n", " * time: 0.5483200550079346\n", " 136 7.606549e+00 2.093485e-02\n", " * time: 0.5514931678771973\n", " 137 7.606460e+00 3.634124e-02\n", " * time: 0.5653150081634521\n", " 138 7.606391e+00 4.135379e-02\n", " * time: 0.5678691864013672\n", " 139 7.606287e+00 4.141551e-02\n", " * time: 0.5704801082611084\n", " 140 7.606270e+00 3.488872e-02\n", " * time: 0.5729410648345947\n", " 141 7.606256e+00 3.433710e-02\n", " * time: 0.576120138168335\n", " 142 7.606151e+00 3.133224e-02\n", " * time: 0.5785551071166992\n", " 143 7.606061e+00 3.056618e-02\n", " * time: 0.5812702178955078\n", " 144 7.606030e+00 5.119188e-02\n", " * time: 0.5838971138000488\n", " 145 7.605886e+00 3.245798e-02\n", " * time: 0.5872189998626709\n", " 146 7.605858e+00 3.144503e-02\n", " * time: 0.5896902084350586\n", " 147 7.605730e+00 2.435330e-02\n", " * time: 0.5920751094818115\n", " 148 7.605636e+00 4.338568e-02\n", " * time: 0.5948851108551025\n", " 149 7.605567e+00 2.436046e-02\n", " * time: 0.5974130630493164\n", " 150 7.605482e+00 1.918499e-02\n", " * time: 0.5997841358184814\n", " 151 7.605437e+00 3.934304e-02\n", " * time: 0.6024432182312012\n", " 152 7.605373e+00 2.995829e-02\n", " * time: 0.6050770282745361\n", " 153 7.605259e+00 2.396337e-02\n", " * time: 0.6082181930541992\n", " 154 7.605181e+00 3.531049e-02\n", " * time: 0.6108341217041016\n", " 155 7.605061e+00 2.253787e-02\n", " * time: 0.614328145980835\n", " 156 7.604976e+00 1.682726e-02\n", " * time: 0.6170690059661865\n", " 157 7.604935e+00 2.973387e-02\n", " * time: 0.6197750568389893\n", " 158 7.604873e+00 2.313383e-02\n", " * time: 0.622528076171875\n", " 159 7.604791e+00 2.734255e-02\n", " * time: 0.6253671646118164\n", " 160 7.604682e+00 1.430313e-02\n", " * time: 0.628911018371582\n", " 161 7.604633e+00 2.979951e-02\n", " * time: 0.6312961578369141\n", " 162 7.604546e+00 1.461828e-02\n", " * time: 0.6352560520172119\n", " 163 7.604489e+00 2.104824e-02\n", " * time: 0.6380300521850586\n", " 164 7.604419e+00 1.614867e-02\n", " * time: 0.6414940357208252\n", " 165 7.604369e+00 1.317928e-02\n", " * time: 0.6448941230773926\n", " 166 7.604325e+00 1.880639e-02\n", " * time: 0.647834062576294\n", " 167 7.604262e+00 1.106014e-02\n", " * time: 0.6514310836791992\n", " 168 7.604218e+00 2.193885e-02\n", " * time: 0.65431809425354\n", " 169 7.604185e+00 1.288283e-02\n", " * time: 0.6571221351623535\n", " 170 7.604131e+00 1.445975e-02\n", " * time: 0.6710870265960693\n", " 171 7.604057e+00 1.564109e-02\n", " * time: 0.6736660003662109\n", " 172 7.604052e+00 2.595458e-02\n", " * time: 0.6766831874847412\n", " 173 7.604017e+00 1.512751e-02\n", " * time: 0.6791121959686279\n", " 174 7.603964e+00 1.830220e-02\n", " * time: 0.6817021369934082\n", " 175 7.603954e+00 1.628212e-02\n", " * time: 0.6845521926879883\n", " 176 7.603882e+00 1.315696e-02\n", " * time: 0.6871931552886963\n", " 177 7.603822e+00 2.061549e-02\n", " * time: 0.6895661354064941\n", " 178 7.603800e+00 2.522680e-02\n", " * time: 0.6920011043548584\n", " 179 7.603724e+00 2.240699e-02\n", " * time: 0.6956660747528076\n", " 180 7.603687e+00 1.530921e-02\n", " * time: 0.6981790065765381\n", " 181 7.603647e+00 2.062737e-02\n", " * time: 0.7006330490112305\n", " 182 7.603579e+00 2.079025e-02\n", " * time: 0.7032451629638672\n", " 183 7.603507e+00 1.072511e-02\n", " * time: 0.706855058670044\n", " 184 7.603460e+00 1.567670e-02\n", " * time: 0.7102530002593994\n", " 185 7.603452e+00 2.889062e-02\n", " * time: 0.7130792140960693\n", " 186 7.603349e+00 1.421094e-02\n", " * time: 0.7168900966644287\n", " 187 7.603315e+00 2.091882e-02\n", " * time: 0.7196621894836426\n", " 188 7.603310e+00 2.132842e-02\n", " * time: 0.7223532199859619\n", " 189 7.603251e+00 1.918126e-02\n", " * time: 0.725999116897583\n", " 190 7.603214e+00 1.343097e-02\n", " * time: 0.7289071083068848\n", " 191 7.603174e+00 1.154620e-02\n", " * time: 0.7314870357513428\n", " 192 7.603159e+00 1.611538e-02\n", " * time: 0.7344250679016113\n", " 193 7.603158e+00 1.982741e-02\n", " * time: 0.738037109375\n", " 194 7.603130e+00 1.389574e-02\n", " * time: 0.7407591342926025\n", " 195 7.603101e+00 1.557255e-02\n", " * time: 0.7437372207641602\n", " 196 7.603055e+00 1.140129e-02\n", " * time: 0.7466931343078613\n", " 197 7.603037e+00 1.213112e-02\n", " * time: 0.7495810985565186\n", " 198 7.603014e+00 2.186182e-02\n", " * time: 0.752251148223877\n", " 199 7.602979e+00 9.909920e-03\n", " * time: 0.7558391094207764\n", " 200 7.602963e+00 1.447839e-02\n", " * time: 0.7592170238494873\n", " 201 7.602934e+00 9.358954e-03\n", " * time: 0.7628991603851318\n", " 202 7.602915e+00 7.787919e-03\n", " * time: 0.7667181491851807\n", " 203 7.602894e+00 7.530685e-03\n", " * time: 0.7806620597839355\n", " 204 7.602861e+00 7.078135e-03\n", " * time: 0.7837722301483154\n", " 205 7.602858e+00 1.217531e-02\n", " * time: 0.7863330841064453\n", " 206 7.602856e+00 1.211488e-02\n", " * time: 0.788870096206665\n", " 207 7.602840e+00 1.472989e-02\n", " * time: 0.7913720607757568\n", " 208 7.602821e+00 1.116947e-02\n", " * time: 0.7937140464782715\n", " 209 7.602798e+00 1.428451e-02\n", " * time: 0.7967932224273682\n", " 210 7.602777e+00 1.644955e-02\n", " * time: 0.799314022064209\n", " 211 7.602760e+00 1.817088e-02\n", " * time: 0.8018860816955566\n", " 212 7.602739e+00 1.727356e-02\n", " * time: 0.8048431873321533\n", " 213 7.602687e+00 1.441276e-02\n", " * time: 0.8077700138092041\n", " 214 7.602616e+00 1.111332e-02\n", " * time: 0.8105521202087402\n", " 215 7.602587e+00 1.944929e-02\n", " * time: 0.8133060932159424\n", " 216 7.602519e+00 2.105641e-02\n", " * time: 0.8170931339263916\n", " 217 7.602479e+00 9.963882e-03\n", " * time: 0.8204522132873535\n", " 218 7.602476e+00 1.880738e-02\n", " * time: 0.8231351375579834\n", " 219 7.602438e+00 1.553928e-02\n", " * time: 0.8266701698303223\n", " 220 7.602410e+00 1.637626e-02\n", " * time: 0.8293631076812744\n", " 221 7.602406e+00 2.926086e-02\n", " * time: 0.8318970203399658\n", " 222 7.602404e+00 2.392553e-02\n", " * time: 0.8348062038421631\n", " 223 7.602334e+00 1.606513e-02\n", " * time: 0.838820219039917\n", " 224 7.602282e+00 2.045806e-02\n", " * time: 0.8414180278778076\n", " 225 7.602223e+00 1.742602e-02\n", " * time: 0.8441140651702881\n", " 226 7.602170e+00 9.252147e-03\n", " * time: 0.8476741313934326\n", " 227 7.602138e+00 1.528503e-02\n", " * time: 0.8507680892944336\n", " 228 7.602137e+00 2.104744e-02\n", " * time: 0.8535001277923584\n", " 229 7.602096e+00 1.274027e-02\n", " * time: 0.8565280437469482\n", " 230 7.602052e+00 1.674052e-02\n", " * time: 0.8593661785125732\n", " 231 7.602011e+00 2.009477e-02\n", " * time: 0.8621690273284912\n", " 232 7.601967e+00 1.217068e-02\n", " * time: 0.865976095199585\n", " 233 7.601936e+00 9.377714e-03\n", " * time: 0.869725227355957\n", " 234 7.601900e+00 1.319294e-02\n", " * time: 0.8726131916046143\n", " 235 7.601862e+00 1.221589e-02\n", " * time: 0.8875560760498047\n", " 236 7.601853e+00 1.168590e-02\n", " * time: 0.8903930187225342\n", " 237 7.601832e+00 1.281590e-02\n", " * time: 0.8940742015838623\n", " 238 7.601815e+00 2.229555e-02\n", " * time: 0.8970870971679688\n", " 239 7.601772e+00 1.358515e-02\n", " * time: 0.9007291793823242\n", " 240 7.601767e+00 1.984844e-02\n", " * time: 0.9035701751708984\n", " 241 7.601741e+00 1.714341e-02\n", " * time: 0.9063730239868164\n", " 242 7.601705e+00 1.238754e-02\n", " * time: 0.9100620746612549\n", " 243 7.601675e+00 1.386018e-02\n", " * time: 0.9136781692504883\n", " 244 7.601654e+00 6.810668e-03\n", " * time: 0.9175910949707031\n", " 245 7.601625e+00 1.097817e-02\n", " * time: 0.920457124710083\n", " 246 7.601604e+00 6.149129e-03\n", " * time: 0.924126148223877\n", " 247 7.601595e+00 1.215962e-02\n", " * time: 0.926948070526123\n", " 248 7.601585e+00 1.773146e-02\n", " * time: 0.9299290180206299\n", " 249 7.601582e+00 1.514778e-02\n", " * time: 0.9328811168670654\n", " 250 7.601572e+00 1.556494e-02\n", " * time: 0.9359560012817383\n", " 251 7.601555e+00 1.017987e-02\n", " * time: 0.9389231204986572\n", " 252 7.601520e+00 1.276084e-02\n", " * time: 0.9418160915374756\n", " 253 7.601519e+00 1.391268e-02\n", " * time: 0.9447250366210938\n", " 254 7.601496e+00 1.273370e-02\n", " * time: 0.9476661682128906\n", " 255 7.601477e+00 1.276837e-02\n", " * time: 0.9505460262298584\n", " 256 7.601458e+00 6.797296e-03\n", " * time: 0.9541590213775635\n", " 257 7.601443e+00 1.103274e-02\n", " * time: 0.956963062286377\n", " 258 7.601429e+00 8.010818e-03\n", " * time: 0.9600181579589844\n", " 259 7.601418e+00 6.922436e-03\n", " * time: 0.9637742042541504\n", " 260 7.601404e+00 6.794801e-03\n", " * time: 0.9674580097198486\n", " 261 7.601403e+00 1.501458e-02\n", " * time: 0.9705631732940674\n", " 262 7.601389e+00 1.023094e-02\n", " * time: 0.974175214767456\n", " 263 7.601378e+00 1.124396e-02\n", " * time: 0.9770181179046631\n", " 264 7.601360e+00 7.068256e-03\n", " * time: 0.9809291362762451\n", " 265 7.601350e+00 1.037933e-02\n", " * time: 0.9837570190429688\n", " 266 7.601344e+00 7.878459e-03\n", " * time: 0.9989120960235596\n", " 267 7.601332e+00 7.666572e-03\n", " * time: 1.001777172088623\n", " 268 7.601318e+00 8.702500e-03\n", " * time: 1.0053200721740723\n", " 269 7.601309e+00 1.003456e-02\n", " * time: 1.0082221031188965\n", " 270 7.601306e+00 1.246596e-02\n", " * time: 1.0111031532287598\n", " 271 7.601294e+00 1.117840e-02\n", " * time: 1.0147910118103027\n", " 272 7.601293e+00 1.283813e-02\n", " * time: 1.0175871849060059\n", " 273 7.601271e+00 1.041787e-02\n", " * time: 1.0206432342529297\n", " 274 7.601270e+00 1.438363e-02\n", " * time: 1.023407220840454\n", " 275 7.601247e+00 9.010849e-03\n", " * time: 1.0261600017547607\n", " 276 7.601222e+00 1.039401e-02\n", " * time: 1.0298831462860107\n", " 277 7.601198e+00 6.620045e-03\n", " * time: 1.0335700511932373\n", " 278 7.601188e+00 8.522716e-03\n", " * time: 1.0364551544189453\n", " 279 7.601174e+00 1.115562e-02\n", " * time: 1.0395731925964355\n", " 280 7.601172e+00 1.017156e-02\n", " * time: 1.042565107345581\n", " 281 7.601171e+00 1.221839e-02\n", " * time: 1.0454421043395996\n", " 282 7.601162e+00 9.816724e-03\n", " * time: 1.0489721298217773\n", " 283 7.601151e+00 1.170915e-02\n", " * time: 1.0516722202301025\n", " 284 7.601139e+00 1.014762e-02\n", " * time: 1.0544631481170654\n", " 285 7.601130e+00 7.442361e-03\n", " * time: 1.0573430061340332\n", " 286 7.601126e+00 6.946548e-03\n", " * time: 1.0609040260314941\n", " 287 7.601119e+00 6.924314e-03\n", " * time: 1.0642471313476562\n", " 288 7.601113e+00 6.528324e-03\n", " * time: 1.0668811798095703\n", " 289 7.601107e+00 9.988013e-03\n", " * time: 1.069742202758789\n", " 290 7.601099e+00 5.838370e-03\n", " * time: 1.0735981464385986\n", " 291 7.601097e+00 8.713871e-03\n", " * time: 1.0766210556030273\n", " 292 7.601089e+00 5.986736e-03\n", " * time: 1.080369234085083\n", " 293 7.601085e+00 6.224439e-03\n", " * time: 1.0832850933074951\n", " 294 7.601083e+00 5.922046e-03\n", " * time: 1.0860590934753418\n", " 295 7.601080e+00 4.268503e-03\n", " * time: 1.0896790027618408\n", " 296 7.601078e+00 5.097434e-03\n", " * time: 1.0923621654510498\n", " 297 7.601077e+00 4.374066e-03\n", " * time: 1.0951180458068848\n", " 298 7.601074e+00 3.996528e-03\n", " * time: 1.097919225692749\n", " 299 7.601072e+00 3.471563e-03\n", " * time: 1.1130001544952393\n", " 300 7.601071e+00 4.765078e-03\n", " * time: 1.1158699989318848\n", " 301 7.601069e+00 4.877198e-03\n", " * time: 1.1188631057739258\n", " 302 7.601068e+00 4.694021e-03\n", " * time: 1.1218950748443604\n", " 303 7.601066e+00 3.437848e-03\n", " * time: 1.124805212020874\n", " 304 7.601063e+00 3.363479e-03\n", " * time: 1.1276171207427979\n", " 305 7.601061e+00 1.852106e-03\n", " * time: 1.1312150955200195\n", " 306 7.601060e+00 4.371000e-03\n", " * time: 1.1342811584472656\n", " 307 7.601058e+00 2.832171e-03\n", " * time: 1.1380810737609863\n", " 308 7.601057e+00 1.753196e-03\n", " * time: 1.1419000625610352\n", " 309 7.601056e+00 3.548001e-03\n", " * time: 1.1447391510009766\n", " 310 7.601056e+00 2.687940e-03\n", " * time: 1.1473910808563232\n", " 311 7.601054e+00 2.053018e-03\n", " * time: 1.1510100364685059\n", " 312 7.601054e+00 5.189253e-03\n", " * time: 1.153792142868042\n", " 313 7.601054e+00 3.390722e-03\n", " * time: 1.156494140625\n", " 314 7.601053e+00 2.682925e-03\n", " * time: 1.1591901779174805\n", " 315 7.601051e+00 2.697822e-03\n", " * time: 1.162102222442627\n", " 316 7.601048e+00 2.796926e-03\n", " * time: 1.16493821144104\n", " 317 7.601045e+00 1.421794e-03\n", " * time: 1.1684751510620117\n", " 318 7.601044e+00 3.326643e-03\n", " * time: 1.1712582111358643\n", " 319 7.601043e+00 3.365619e-03\n", " * time: 1.174084186553955\n", " 320 7.601042e+00 2.011107e-03\n", " * time: 1.1774470806121826\n", " 321 7.601041e+00 1.535134e-03\n", " * time: 1.1811082363128662\n", " 322 7.601040e+00 2.149178e-03\n", " * time: 1.1840641498565674\n", " 323 7.601039e+00 5.536894e-03\n", " * time: 1.186805009841919\n", " 324 7.601037e+00 2.966835e-03\n", " * time: 1.1904230117797852\n", " 325 7.601036e+00 3.410147e-03\n", " * time: 1.1932260990142822\n", " 326 7.601035e+00 2.656295e-03\n", " * time: 1.1960890293121338\n", " 327 7.601032e+00 2.456018e-03\n", " * time: 1.198768138885498\n", " 328 7.601032e+00 4.986742e-03\n", " * time: 1.2017600536346436\n", " 329 7.601030e+00 2.632713e-03\n", " * time: 1.2053790092468262\n", " 330 7.601028e+00 2.795960e-03\n", " * time: 1.2088990211486816\n", " 331 7.601027e+00 2.809562e-03\n", " * time: 1.2223761081695557\n", " 332 7.601027e+00 3.655918e-03\n", " * time: 1.2251250743865967\n", " 333 7.601025e+00 3.240803e-03\n", " * time: 1.2279551029205322\n", " 334 7.601023e+00 2.940965e-03\n", " * time: 1.231316089630127\n", " 335 7.601022e+00 2.730641e-03\n", " * time: 1.2340002059936523\n", " 336 7.601020e+00 1.946138e-03\n", " * time: 1.2374930381774902\n", " 337 7.601019e+00 3.175974e-03\n", " * time: 1.240260124206543\n", " 338 7.601019e+00 2.416396e-03\n", " * time: 1.2433140277862549\n", " 339 7.601018e+00 3.363247e-03\n", " * time: 1.2464051246643066\n", " 340 7.601017e+00 2.081432e-03\n", " * time: 1.2492311000823975\n", " 341 7.601017e+00 3.599772e-03\n", " * time: 1.2520630359649658\n", " 342 7.601016e+00 3.784154e-03\n", " * time: 1.2549240589141846\n", " 343 7.601015e+00 4.508210e-03\n", " * time: 1.2579841613769531\n", " 344 7.601015e+00 3.241701e-03\n", " * time: 1.2608191967010498\n", " 345 7.601011e+00 3.603575e-03\n", " * time: 1.2637920379638672\n", " 346 7.601008e+00 3.249021e-03\n", " * time: 1.2672860622406006\n", " 347 7.601006e+00 2.880308e-03\n", " * time: 1.269968032836914\n", " 348 7.601004e+00 3.146072e-03\n", " * time: 1.273529052734375\n", " 349 7.601004e+00 3.987089e-03\n", " * time: 1.2763221263885498\n", " 350 7.601000e+00 2.841271e-03\n", " * time: 1.2796711921691895\n", " 351 7.600997e+00 3.213333e-03\n", " * time: 1.2826220989227295\n", " 352 7.600994e+00 2.829795e-03\n", " * time: 1.2864811420440674\n", " 353 7.600993e+00 3.391271e-03\n", " * time: 1.2901430130004883\n", " 354 7.600988e+00 3.787108e-03\n", " * time: 1.2937891483306885\n", " 355 7.600985e+00 3.825092e-03\n", " * time: 1.2967591285705566\n", " 356 7.600982e+00 3.276053e-03\n", " * time: 1.2994861602783203\n", " 357 7.600978e+00 3.517881e-03\n", " * time: 1.302321195602417\n", " 358 7.600974e+00 3.926654e-03\n", " * time: 1.3051741123199463\n", " 359 7.600969e+00 4.748275e-03\n", " * time: 1.3082921504974365\n", " 360 7.600966e+00 3.527200e-03\n", " * time: 1.3118622303009033\n", " 361 7.600964e+00 2.447500e-03\n", " * time: 1.315331220626831\n", " 362 7.600962e+00 2.084165e-03\n", " * time: 1.3177790641784668\n", " 363 7.600961e+00 3.277194e-03\n", " * time: 1.3311331272125244\n", " 364 7.600959e+00 2.720901e-03\n", " * time: 1.3337352275848389\n", " 365 7.600957e+00 2.783322e-03\n", " * time: 1.3364510536193848\n", " 366 7.600955e+00 3.132031e-03\n", " * time: 1.3389880657196045\n", " 367 7.600953e+00 1.823816e-03\n", " * time: 1.3422431945800781\n", " 368 7.600953e+00 3.347557e-03\n", " * time: 1.3450751304626465\n", " 369 7.600953e+00 4.321989e-03\n", " * time: 1.348405122756958\n", " 370 7.600952e+00 2.848158e-03\n", " * time: 1.3509151935577393\n", " 371 7.600951e+00 2.884284e-03\n", " * time: 1.3539221286773682\n", " 372 7.600951e+00 3.022441e-03\n", " * time: 1.3567030429840088\n", " 373 7.600949e+00 1.969183e-03\n", " * time: 1.3594632148742676\n", " 374 7.600947e+00 1.437397e-03\n", " * time: 1.3630990982055664\n", " 375 7.600947e+00 3.470440e-03\n", " * time: 1.3657991886138916\n", " 376 7.600946e+00 2.106716e-03\n", " * time: 1.3692522048950195\n", " 377 7.600945e+00 1.979013e-03\n", " * time: 1.3717920780181885\n", " 378 7.600944e+00 2.310452e-03\n", " * time: 1.3746371269226074\n", " 379 7.600944e+00 4.517373e-03\n", " * time: 1.3771820068359375\n", " 380 7.600943e+00 3.255648e-03\n", " * time: 1.3803400993347168\n", " 381 7.600942e+00 3.305324e-03\n", " * time: 1.3840620517730713\n", " 382 7.600940e+00 3.389194e-03\n", " * time: 1.3875601291656494\n", " 383 7.600940e+00 3.477687e-03\n", " * time: 1.3900542259216309\n", " 384 7.600937e+00 3.100839e-03\n", " * time: 1.3927481174468994\n", " 385 7.600936e+00 3.305250e-03\n", " * time: 1.3956871032714844\n", " 386 7.600933e+00 2.624755e-03\n", " * time: 1.39925217628479\n", " 387 7.600932e+00 3.658885e-03\n", " * time: 1.4016940593719482\n", " 388 7.600931e+00 4.039580e-03\n", " * time: 1.404569149017334\n", " 389 7.600927e+00 3.218110e-03\n", " * time: 1.4073851108551025\n", " 390 7.600924e+00 2.753180e-03\n", " * time: 1.4109601974487305\n", " 391 7.600921e+00 2.613062e-03\n", " * time: 1.4138641357421875\n", " 392 7.600920e+00 4.613085e-03\n", " * time: 1.41664719581604\n", " 393 7.600918e+00 3.952392e-03\n", " * time: 1.419206142425537\n", " 394 7.600915e+00 2.220950e-03\n", " * time: 1.4223790168762207\n", " 395 7.600913e+00 2.064660e-03\n", " * time: 1.4358410835266113\n", " 396 7.600912e+00 3.418463e-03\n", " * time: 1.4383931159973145\n", " 397 7.600910e+00 3.180001e-03\n", " * time: 1.4408960342407227\n", " 398 7.600909e+00 4.785496e-03\n", " * time: 1.4438412189483643\n", " 399 7.600906e+00 2.974750e-03\n", " * time: 1.4474220275878906\n", " 400 7.600903e+00 3.022774e-03\n", " * time: 1.450120210647583\n", " 401 7.600902e+00 5.086027e-03\n", " * time: 1.4528100490570068\n", " 402 7.600901e+00 4.294586e-03\n", " * time: 1.4556801319122314\n", " 403 7.600898e+00 3.590096e-03\n", " * time: 1.4592971801757812\n", " 404 7.600896e+00 2.443877e-03\n", " * time: 1.4621551036834717\n", " 405 7.600894e+00 2.732149e-03\n", " * time: 1.4661800861358643\n", " 406 7.600891e+00 2.201315e-03\n", " * time: 1.469149112701416\n", " 407 7.600889e+00 2.283525e-03\n", " * time: 1.4726731777191162\n", " 408 7.600887e+00 3.475261e-03\n", " * time: 1.4755771160125732\n", " 409 7.600885e+00 1.884211e-03\n", " * time: 1.4791820049285889\n", " 410 7.600885e+00 2.614432e-03\n", " * time: 1.4826841354370117\n", " 411 7.600884e+00 2.989126e-03\n", " * time: 1.4858450889587402\n", " 412 7.600883e+00 2.551096e-03\n", " * time: 1.4893360137939453\n", " 413 7.600881e+00 1.396814e-03\n", " * time: 1.4929330348968506\n", " 414 7.600881e+00 2.814465e-03\n", " * time: 1.495812177658081\n", " 415 7.600880e+00 1.564532e-03\n", " * time: 1.4982120990753174\n", " 416 7.600879e+00 2.007180e-03\n", " * time: 1.5007801055908203\n", " 417 7.600879e+00 2.331050e-03\n", " * time: 1.503521203994751\n", " 418 7.600879e+00 2.637553e-03\n", " * time: 1.5065982341766357\n", " 419 7.600877e+00 2.272948e-03\n", " * time: 1.5094881057739258\n", " 420 7.600876e+00 1.907767e-03\n", " * time: 1.5120491981506348\n", " 421 7.600876e+00 2.144734e-03\n", " * time: 1.5146090984344482\n", " 422 7.600875e+00 1.528154e-03\n", " * time: 1.5180981159210205\n", " 423 7.600874e+00 2.042848e-03\n", " * time: 1.5207040309906006\n", " 424 7.600874e+00 1.522194e-03\n", " * time: 1.523273229598999\n", " 425 7.600873e+00 1.634656e-03\n", " * time: 1.5261850357055664\n", " 426 7.600873e+00 1.623462e-03\n", " * time: 1.5289912223815918\n", " 427 7.600872e+00 1.050503e-03\n", " * time: 1.5431230068206787\n", " 428 7.600872e+00 2.611472e-03\n", " * time: 1.5460131168365479\n", " 429 7.600871e+00 1.196147e-03\n", " * time: 1.54901123046875\n", " 430 7.600871e+00 1.966287e-03\n", " * time: 1.5513651371002197\n", " 431 7.600871e+00 3.125490e-03\n", " * time: 1.5538511276245117\n", " 432 7.600870e+00 1.895833e-03\n", " * time: 1.557035207748413\n", " 433 7.600870e+00 2.000211e-03\n", " * time: 1.5594151020050049\n", " 434 7.600868e+00 1.321102e-03\n", " * time: 1.5625522136688232\n", " 435 7.600868e+00 2.455830e-03\n", " * time: 1.5654630661010742\n", " 436 7.600867e+00 2.270062e-03\n", " * time: 1.5681521892547607\n", " 437 7.600867e+00 2.249795e-03\n", " * time: 1.5709660053253174\n", " 438 7.600866e+00 2.059805e-03\n", " * time: 1.5737340450286865\n", " 439 7.600865e+00 1.179741e-03\n", " * time: 1.5772771835327148\n", " 440 7.600864e+00 2.101183e-03\n", " * time: 1.5798392295837402\n", " 441 7.600864e+00 2.915310e-03\n", " * time: 1.5822641849517822\n", " 442 7.600863e+00 2.389144e-03\n", " * time: 1.5860822200775146\n", " 443 7.600863e+00 1.953627e-03\n", " * time: 1.5895960330963135\n", " 444 7.600862e+00 2.182866e-03\n", " * time: 1.5920422077178955\n", " 445 7.600861e+00 2.299549e-03\n", " * time: 1.594452142715454\n", " 446 7.600860e+00 1.690470e-03\n", " * time: 1.5982232093811035\n", " 447 7.600859e+00 1.200751e-03\n", " * time: 1.6015682220458984\n", " 448 7.600858e+00 9.687741e-04\n", " * time: 1.60545015335083\n", " 449 7.600858e+00 1.179866e-03\n", " * time: 1.608206033706665\n", " 450 7.600858e+00 1.859813e-03\n", " * time: 1.6107380390167236\n", " 451 7.600857e+00 2.048209e-03\n", " * time: 1.6139240264892578\n", " 452 7.600856e+00 8.036487e-04\n", " * time: 1.617845058441162\n", " 453 7.600856e+00 1.625258e-03\n", " * time: 1.6206481456756592\n", " 454 7.600855e+00 1.780060e-03\n", " * time: 1.6230981349945068\n", " 455 7.600855e+00 1.560579e-03\n", " * time: 1.625866174697876\n", " 456 7.600854e+00 1.407059e-03\n", " * time: 1.6284761428833008\n", " 457 7.600853e+00 1.847938e-03\n", " * time: 1.6312100887298584\n", " 458 7.600853e+00 8.683814e-04\n", " * time: 1.6349401473999023\n", " 459 7.600853e+00 1.696652e-03\n", " * time: 1.6485111713409424\n", " 460 7.600852e+00 2.345996e-03\n", " * time: 1.6513621807098389\n", " 461 7.600852e+00 1.624941e-03\n", " * time: 1.6540400981903076\n", " 462 7.600851e+00 1.142548e-03\n", " * time: 1.656865119934082\n", " 463 7.600851e+00 1.097551e-03\n", " * time: 1.6595802307128906\n", " 464 7.600851e+00 2.231028e-03\n", " * time: 1.6621990203857422\n", " 465 7.600850e+00 1.323450e-03\n", " * time: 1.6659550666809082\n", " 466 7.600850e+00 1.312195e-03\n", " * time: 1.6694831848144531\n", " 467 7.600850e+00 1.932481e-03\n", " * time: 1.67205810546875\n", " 468 7.600850e+00 1.224854e-03\n", " * time: 1.6745431423187256\n", " 469 7.600849e+00 1.272548e-03\n", " * time: 1.6773931980133057\n", " 470 7.600849e+00 9.671457e-04\n", " * time: 1.6807892322540283\n", " 471 7.600849e+00 1.283214e-03\n", " * time: 1.683394193649292\n", " 472 7.600849e+00 8.662530e-04\n", " * time: 1.687068223953247\n", " 473 7.600849e+00 8.675408e-04\n", " * time: 1.6896920204162598\n", " 474 7.600848e+00 9.236489e-04\n", " * time: 1.6924190521240234\n", " 475 7.600848e+00 1.097564e-03\n", " * time: 1.69527006149292\n", " 476 7.600848e+00 1.211453e-03\n", " * time: 1.698173999786377\n", " 477 7.600848e+00 8.564432e-04\n", " * time: 1.700857162475586\n", " 478 7.600848e+00 1.346398e-03\n", " * time: 1.7033510208129883\n", " 479 7.600848e+00 1.088860e-03\n", " * time: 1.7070081233978271\n", " 480 7.600848e+00 1.180062e-03\n", " * time: 1.7097690105438232\n", " 481 7.600847e+00 1.203945e-03\n", " * time: 1.7126522064208984\n", " 482 7.600847e+00 7.799915e-04\n", " * time: 1.715994119644165\n", " 483 7.600847e+00 7.120750e-04\n", " * time: 1.7187330722808838\n", " 484 7.600847e+00 7.958916e-04\n", " * time: 1.7214360237121582\n", " 485 7.600847e+00 7.431726e-04\n", " * time: 1.7241160869598389\n", " 486 7.600847e+00 1.058980e-03\n", " * time: 1.7271392345428467\n", " 487 7.600846e+00 1.307090e-03\n", " * time: 1.7297492027282715\n", " 488 7.600846e+00 1.384443e-03\n", " * time: 1.7329370975494385\n", " 489 7.600846e+00 1.517372e-03\n", " * time: 1.7364330291748047\n", " 490 7.600846e+00 1.535546e-03\n", " * time: 1.7394630908966064\n", " 491 7.600845e+00 8.688543e-04\n", " * time: 1.7535011768341064\n", " 492 7.600845e+00 7.929016e-04\n", " * time: 1.7563791275024414\n", " 493 7.600845e+00 8.689368e-04\n", " * time: 1.7592260837554932\n", " 494 7.600844e+00 5.961645e-04\n", " * time: 1.762578010559082\n", " 495 7.600844e+00 1.772172e-03\n", " * time: 1.7650401592254639\n", " 496 7.600844e+00 1.552537e-03\n", " * time: 1.7678430080413818\n", " 497 7.600844e+00 1.253552e-03\n", " * time: 1.771374225616455\n", " 498 7.600843e+00 6.874548e-04\n", " * time: 1.7750310897827148\n", " 499 7.600843e+00 1.188100e-03\n", " * time: 1.777787208557129\n", " 500 7.600843e+00 9.473894e-04\n", " * time: 1.78041410446167\n", " 501 7.600843e+00 6.384992e-04\n", " * time: 1.7837440967559814\n", " 502 7.600842e+00 1.120220e-03\n", " * time: 1.7862920761108398\n", " 503 7.600842e+00 1.303880e-03\n", " * time: 1.7895021438598633\n", " 504 7.600842e+00 7.521000e-04\n", " * time: 1.7924470901489258\n", " 505 7.600842e+00 1.149497e-03\n", " * time: 1.7952520847320557\n", " 506 7.600842e+00 9.815085e-04\n", " * time: 1.798840045928955\n", " 507 7.600842e+00 1.145894e-03\n", " * time: 1.801774024963379\n", " 508 7.600841e+00 7.756741e-04\n", " * time: 1.8047630786895752\n", " 509 7.600841e+00 7.590623e-04\n", " * time: 1.807784080505371\n", " 510 7.600841e+00 1.071154e-03\n", " * time: 1.810617208480835\n", " 511 7.600841e+00 6.880179e-04\n", " * time: 1.8139400482177734\n", " 512 7.600841e+00 9.612964e-04\n", " * time: 1.8168540000915527\n", " 513 7.600841e+00 7.712132e-04\n", " * time: 1.8207461833953857\n", " 514 7.600840e+00 6.974839e-04\n", " * time: 1.8243160247802734\n", " 515 7.600840e+00 1.018041e-03\n", " * time: 1.8271100521087646\n", " 516 7.600840e+00 1.372314e-03\n", " * time: 1.8301401138305664\n", " 517 7.600840e+00 1.554752e-03\n", " * time: 1.8329720497131348\n", " 518 7.600840e+00 7.440250e-04\n", " * time: 1.8368401527404785\n", " 519 7.600840e+00 1.154033e-03\n", " * time: 1.8398971557617188\n", " 520 7.600840e+00 1.024308e-03\n", " * time: 1.8427140712738037\n", " 521 7.600840e+00 1.378447e-03\n", " * time: 1.8454630374908447\n", " 522 7.600840e+00 9.240342e-04\n", " * time: 1.8485820293426514\n", " 523 7.600839e+00 1.359927e-03\n", " * time: 1.8516011238098145\n", " 524 7.600839e+00 9.373707e-04\n", " * time: 1.865535020828247\n", " 525 7.600839e+00 1.999338e-03\n", " * time: 1.8686110973358154\n", " 526 7.600839e+00 1.239880e-03\n", " * time: 1.8721020221710205\n", " 527 7.600838e+00 1.169403e-03\n", " * time: 1.8754780292510986\n", " 528 7.600838e+00 1.260677e-03\n", " * time: 1.878899097442627\n", " 529 7.600838e+00 1.288191e-03\n", " * time: 1.8814690113067627\n", " 530 7.600837e+00 1.132763e-03\n", " * time: 1.8850951194763184\n", " 531 7.600837e+00 1.667892e-03\n", " * time: 1.8879292011260986\n", " 532 7.600837e+00 2.528681e-03\n", " * time: 1.8909780979156494\n", " 533 7.600837e+00 1.845218e-03\n", " * time: 1.8938381671905518\n", " 534 7.600836e+00 1.630686e-03\n", " * time: 1.8966801166534424\n", " 535 7.600836e+00 1.812725e-03\n", " * time: 1.8996570110321045\n", " 536 7.600835e+00 1.098980e-03\n", " * time: 1.9030210971832275\n", " 537 7.600834e+00 1.595411e-03\n", " * time: 1.9058570861816406\n", " 538 7.600834e+00 1.083912e-03\n", " * time: 1.9093971252441406\n", " 539 7.600833e+00 1.844284e-03\n", " * time: 1.912255048751831\n", " 540 7.600833e+00 1.900472e-03\n", " * time: 1.9154231548309326\n", " 541 7.600833e+00 1.731638e-03\n", " * time: 1.9183480739593506\n", " 542 7.600832e+00 1.682862e-03\n", " * time: 1.9211840629577637\n", " 543 7.600832e+00 1.009589e-03\n", " * time: 1.9245731830596924\n", " 544 7.600831e+00 1.651865e-03\n", " * time: 1.9272441864013672\n", " 545 7.600831e+00 1.774669e-03\n", " * time: 1.9302051067352295\n", " 546 7.600831e+00 1.702425e-03\n", " * time: 1.9332962036132812\n", " 547 7.600831e+00 1.927276e-03\n", " * time: 1.9364430904388428\n", " 548 7.600830e+00 2.718352e-03\n", " * time: 1.9392871856689453\n", " 549 7.600830e+00 1.237155e-03\n", " * time: 1.94199800491333\n", " 550 7.600830e+00 9.758484e-04\n", " * time: 1.9455900192260742\n", " 551 7.600829e+00 8.857924e-04\n", " * time: 1.9492170810699463\n", " 552 7.600829e+00 7.833400e-04\n", " * time: 1.952923059463501\n", " 553 7.600829e+00 9.658098e-04\n", " * time: 1.955801010131836\n", " 554 7.600829e+00 1.258714e-03\n", " * time: 1.9587700366973877\n", " 555 7.600828e+00 1.622762e-03\n", " * time: 1.9618501663208008\n", " 556 7.600828e+00 1.016754e-03\n", " * time: 1.97737717628479\n", " 557 7.600828e+00 8.554805e-04\n", " * time: 1.9810750484466553\n", " 558 7.600828e+00 1.606530e-03\n", " * time: 1.9837441444396973\n", " 559 7.600828e+00 1.403772e-03\n", " * time: 1.9870381355285645\n", " 560 7.600828e+00 1.749462e-03\n", " * time: 1.9900751113891602\n", " 561 7.600827e+00 9.794210e-04\n", " * time: 1.993037223815918\n", " 562 7.600827e+00 9.441680e-04\n", " * time: 1.9958980083465576\n", " 563 7.600827e+00 1.114261e-03\n", " * time: 1.9986932277679443\n", " 564 7.600827e+00 1.232025e-03\n", " * time: 2.001591205596924\n", " 565 7.600827e+00 1.248178e-03\n", " * time: 2.0044941902160645\n", " 566 7.600827e+00 1.600084e-03\n", " * time: 2.0072951316833496\n", " 567 7.600827e+00 1.231672e-03\n", " * time: 2.0103862285614014\n", " 568 7.600826e+00 1.139302e-03\n", " * time: 2.013277053833008\n", " 569 7.600826e+00 6.801595e-04\n", " * time: 2.017000198364258\n", " 570 7.600826e+00 1.265150e-03\n", " * time: 2.020055055618286\n", " 571 7.600826e+00 6.900613e-04\n", " * time: 2.023757219314575\n", " 572 7.600826e+00 5.364120e-04\n", " * time: 2.0265932083129883\n", " 573 7.600826e+00 4.886742e-04\n", " * time: 2.0293710231781006\n", " 574 7.600826e+00 7.357897e-04\n", " * time: 2.033388137817383\n", " 575 7.600826e+00 6.481461e-04\n", " * time: 2.0369691848754883\n", " 576 7.600825e+00 5.096648e-04\n", " * time: 2.04072904586792\n", " 577 7.600825e+00 3.530884e-04\n", " * time: 2.0441431999206543\n", " 578 7.600825e+00 5.511537e-04\n", " * time: 2.046901226043701\n", " 579 7.600825e+00 3.825295e-04\n", " * time: 2.0496160984039307\n", " 580 7.600825e+00 5.898895e-04\n", " * time: 2.053046226501465\n", " 581 7.600825e+00 1.024659e-03\n", " * time: 2.055957078933716\n", " 582 7.600825e+00 8.733688e-04\n", " * time: 2.0589001178741455\n", " 583 7.600825e+00 6.154697e-04\n", " * time: 2.062727212905884\n", " 584 7.600825e+00 8.800110e-04\n", " * time: 2.0657150745391846\n", " 585 7.600824e+00 5.817669e-04\n", " * time: 2.068908214569092\n", " 586 7.600824e+00 9.840813e-04\n", " * time: 2.07204008102417\n", " 587 7.600824e+00 9.398667e-04\n", " * time: 2.074770212173462\n", " 588 7.600824e+00 8.762542e-04\n", " * time: 2.0895321369171143\n", " 589 7.600824e+00 1.128688e-03\n", " * time: 2.0928330421447754\n", " 590 7.600824e+00 8.842296e-04\n", " * time: 2.0955910682678223\n", " 591 7.600824e+00 9.544369e-04\n", " * time: 2.0983800888061523\n", " 592 7.600824e+00 8.024149e-04\n", " * time: 2.102081060409546\n", " 593 7.600824e+00 6.919013e-04\n", " * time: 2.1048390865325928\n", " 594 7.600823e+00 9.116241e-04\n", " * time: 2.107722043991089\n", " 595 7.600823e+00 5.052231e-04\n", " * time: 2.1118972301483154\n", " 596 7.600823e+00 6.182628e-04\n", " * time: 2.1147971153259277\n", " 597 7.600823e+00 7.097004e-04\n", " * time: 2.1185522079467773\n", " 598 7.600823e+00 4.535988e-04\n", " * time: 2.1224582195281982\n", " 599 7.600823e+00 4.792520e-04\n", " * time: 2.1260762214660645\n", " 600 7.600823e+00 9.074247e-04\n", " * time: 2.128958225250244\n", " 601 7.600823e+00 5.084135e-04\n", " * time: 2.1320300102233887\n", " 602 7.600823e+00 4.492744e-04\n", " * time: 2.135094165802002\n", " 603 7.600823e+00 9.010732e-04\n", " * time: 2.1381330490112305\n", " 604 7.600823e+00 8.673113e-04\n", " * time: 2.1410181522369385\n", " 605 7.600823e+00 7.551361e-04\n", " * time: 2.144052028656006\n", " 606 7.600823e+00 6.933345e-04\n", " * time: 2.146827220916748\n", " 607 7.600823e+00 5.448176e-04\n", " * time: 2.1504571437835693\n", " 608 7.600822e+00 8.077135e-04\n", " * time: 2.1535820960998535\n", " 609 7.600822e+00 6.614461e-04\n", " * time: 2.1564102172851562\n", " 610 7.600822e+00 9.118665e-04\n", " * time: 2.159210205078125\n", " 611 7.600822e+00 4.364849e-04\n", " * time: 2.1629791259765625\n", " 612 7.600822e+00 8.501459e-04\n", " * time: 2.165811061859131\n", " 613 7.600822e+00 4.720814e-04\n", " * time: 2.168727159500122\n", " 614 7.600822e+00 1.111033e-03\n", " * time: 2.1717021465301514\n", " 615 7.600822e+00 1.222880e-03\n", " * time: 2.17466402053833\n", " 616 7.600822e+00 1.006596e-03\n", " * time: 2.177612066268921\n", " 617 7.600822e+00 6.853244e-04\n", " * time: 2.1805851459503174\n", " 618 7.600822e+00 8.469811e-04\n", " * time: 2.1835250854492188\n", " 619 7.600821e+00 6.676371e-04\n", " * time: 2.1871800422668457\n", " 620 7.600821e+00 6.417877e-04\n", " * time: 2.201423168182373\n", " 621 7.600821e+00 3.756695e-04\n", " * time: 2.2051000595092773\n", " 622 7.600821e+00 7.742260e-04\n", " * time: 2.20758318901062\n", " 623 7.600821e+00 3.941797e-04\n", " * time: 2.211130142211914\n", " 624 7.600821e+00 6.281014e-04\n", " * time: 2.213970184326172\n", " 625 7.600821e+00 4.733041e-04\n", " * time: 2.2165842056274414\n", " 626 7.600821e+00 4.399552e-04\n", " * time: 2.2199552059173584\n", " 627 7.600821e+00 3.405235e-04\n", " * time: 2.223527193069458\n", " 628 7.600821e+00 3.760416e-04\n", " * time: 2.226391077041626\n", " 629 7.600821e+00 2.798225e-04\n", " * time: 2.229098081588745\n", " 630 7.600821e+00 3.700714e-04\n", " * time: 2.2317922115325928\n", " 631 7.600821e+00 4.337665e-04\n", " * time: 2.234724998474121\n", " 632 7.600821e+00 3.307771e-04\n", " * time: 2.2374110221862793\n", " 633 7.600821e+00 5.110619e-04\n", " * time: 2.240143060684204\n", " 634 7.600821e+00 3.771245e-04\n", " * time: 2.2429821491241455\n", " 635 7.600821e+00 3.158347e-04\n", " * time: 2.246613025665283\n", " 636 7.600821e+00 3.436319e-04\n", " * time: 2.249411106109619\n", " 637 7.600821e+00 4.475719e-04\n", " * time: 2.252239227294922\n", " 638 7.600821e+00 3.462347e-04\n", " * time: 2.2552311420440674\n", " 639 7.600820e+00 2.443671e-04\n", " * time: 2.2590270042419434\n", " 640 7.600820e+00 3.211568e-04\n", " * time: 2.2620351314544678\n", " 641 7.600820e+00 3.326460e-04\n", " * time: 2.2650930881500244\n", " 642 7.600820e+00 2.780081e-04\n", " * time: 2.2686221599578857\n", " 643 7.600820e+00 5.720252e-04\n", " * time: 2.2714040279388428\n", " 644 7.600820e+00 2.976044e-04\n", " * time: 2.2752981185913086\n", " 645 7.600820e+00 2.906335e-04\n", " * time: 2.2788290977478027\n", " 646 7.600820e+00 3.828513e-04\n", " * time: 2.281494140625\n", " 647 7.600820e+00 5.143250e-04\n", " * time: 2.284501075744629\n", " 648 7.600820e+00 3.589438e-04\n", " * time: 2.2880260944366455\n", " 649 7.600820e+00 3.038058e-04\n", " * time: 2.290882110595703\n", " 650 7.600820e+00 4.117595e-04\n", " * time: 2.2940421104431152\n", " 651 7.600820e+00 4.150112e-04\n", " * time: 2.297030210494995\n", " 652 7.600820e+00 4.498958e-04\n", " * time: 2.309988021850586\n", " 653 7.600820e+00 5.365927e-04\n", " * time: 2.3128092288970947\n", " 654 7.600820e+00 2.412987e-04\n", " * time: 2.3165721893310547\n", " 655 7.600820e+00 3.229618e-04\n", " * time: 2.3193061351776123\n", " 656 7.600820e+00 2.782164e-04\n", " * time: 2.322134017944336\n", " 657 7.600820e+00 3.361718e-04\n", " * time: 2.3250861167907715\n", " 658 7.600820e+00 2.440254e-04\n", " * time: 2.3285791873931885\n", " 659 7.600820e+00 3.065326e-04\n", " * time: 2.331369161605835\n", " 660 7.600820e+00 2.636625e-04\n", " * time: 2.3353610038757324\n", " 661 7.600820e+00 2.133982e-04\n", " * time: 2.339081048965454\n", " 662 7.600820e+00 2.874086e-04\n", " * time: 2.3427140712738037\n", " 663 7.600820e+00 3.910003e-04\n", " * time: 2.345720052719116\n", " 664 7.600820e+00 2.995252e-04\n", " * time: 2.3487000465393066\n", " 665 7.600820e+00 3.246986e-04\n", " * time: 2.3517541885375977\n", " 666 7.600820e+00 4.390405e-04\n", " * time: 2.354863166809082\n", " 667 7.600820e+00 2.837366e-04\n", " * time: 2.358424186706543\n", " 668 7.600820e+00 3.110767e-04\n", " * time: 2.361316204071045\n", " 669 7.600820e+00 3.071063e-04\n", " * time: 2.3643691539764404\n", " 670 7.600820e+00 4.064561e-04\n", " * time: 2.367112159729004\n", " 671 7.600820e+00 3.708825e-04\n", " * time: 2.370199203491211\n", " 672 7.600820e+00 3.211281e-04\n", " * time: 2.373373031616211\n", " 673 7.600820e+00 4.254178e-04\n", " * time: 2.376434087753296\n", " 674 7.600820e+00 2.239597e-04\n", " * time: 2.3797450065612793\n", " 675 7.600820e+00 2.145667e-04\n", " * time: 2.3824031352996826\n", " 676 7.600820e+00 3.040484e-04\n", " * time: 2.385507106781006\n", " 677 7.600820e+00 2.388455e-04\n", " * time: 2.3883941173553467\n", " 678 7.600820e+00 2.258516e-04\n", " * time: 2.3911831378936768\n", " 679 7.600820e+00 3.513882e-04\n", " * time: 2.394294023513794\n", " 680 7.600820e+00 2.567893e-04\n", " * time: 2.3979220390319824\n", " 681 7.600820e+00 4.056639e-04\n", " * time: 2.4007930755615234\n", " 682 7.600820e+00 2.454788e-04\n", " * time: 2.4046630859375\n", " 683 7.600820e+00 3.998590e-04\n", " * time: 2.407485008239746\n", " 684 7.600820e+00 3.375827e-04\n", " * time: 2.410388231277466\n", " 685 7.600820e+00 2.874519e-04\n", " * time: 2.425387144088745\n", " 686 7.600820e+00 3.121873e-04\n", " * time: 2.4280381202697754\n", " 687 7.600820e+00 2.981234e-04\n", " * time: 2.4316020011901855\n", " 688 7.600820e+00 3.128235e-04\n", " * time: 2.434864044189453\n", " 689 7.600820e+00 3.455741e-04\n", " * time: 2.437636137008667\n", " 690 7.600820e+00 4.512448e-04\n", " * time: 2.440390110015869\n", " 691 7.600820e+00 3.479818e-04\n", " * time: 2.442782163619995\n", " 692 7.600820e+00 3.982809e-04\n", " * time: 2.445650100708008\n", " 693 7.600820e+00 2.982845e-04\n", " * time: 2.4481401443481445\n" ] }, { "output_type": "execute_result", "data": { "text/plain": "Plot{Plots.GRBackend() n=1}", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVxU5f4H8Gdghh0VQVZXVBBUkEjLJTdccMFdzFwqzCW93ay8/cos615Ls82uXdNy7XazXKJrZSVqoGWl5r4gJG4giLKJbDPDzO+PU3Mn+H5xznhgHObzfvnH8PjMmTML8/Cc8znfR2U0GgUAAICjcrL1DgAAANgSBkIAAHBoGAgBAMChYSAEAACHhoEQAAAcGgZCAABwaBgIAQDAoWEgBAAAh4aBEAAAHBoGQgAAcGjqhnkYg8Fw5MSpjhERDfNwDcxgMDg5Nf4/KfA0GxM8zcakwZ5mU1fNbfucO1dYUaGXtdnS0pLgYKf27dtbu193qoEGwqqqqlWfJY98pEXDPBwAAChuXHjQbftMmvTl8ePXZW741JgxJcnJydbt1Z1r/H8rAQAA1KGBZoQAAOAIVCqhUtl6J2TCjBAAABwaZoQAAKAYlUqlkjklNBptPIXEQAgAAIrCoVEAAAA7ghkhAAAoxw7DMhgIAQBAMVacIxS2PkeIQ6MAAODQMBACAIBDw6FRAABQjBUX1Nv6yCgGQgAAUI4V5whtHq7BoVEAAHBomBECAIByVPZ3QT0GQgAAUIxK/qFOmx8axUAIAACKseY6QluPhDhHCAAADg0zQgAAUI4V5whxaBQAABoTexsHZQ6Ex44dy8zMbNGixf333+/m5iY1Hj9+PCMjIzIysnPnzvWwhwAAAPXI0oGwurr6kUceSUtL69mz57Vr16ZOnfrYY48JIZYuXfqvf/1r0KBBf/3rXxcuXPjEE0/U594CAMBdTaUSsi+ov92cMDc3d8OGDcXFxaNHj+7du3ftDnq9fuPGjenp6V26dJk2bZqzs7MQ4sqVK19//fX58+cDAgKmTp0aGBgodTYYDJ988smxY8c6dOiQlJTk4uJiaVjm/fffP3ny5OnTpz/77LPU1NQZM2YIIQoKCpYsWZKSkrJx48Yvv/zyxRdfvHXrloxnDwAAjYzKqn+8kpKS++6779KlS4GBgaNGjfr6669r90lKStqwYUNoaOiqVavmzZsnNU6aNOnnn38OCgo6efJkZGRkVlaW1L5gwYI333yzXbt2n3322eTJk4XlM8KPP/74ySefLCgoOH36dNeuXT09PYUQKSkpHTt2jIiIEELce++9vr6+aWlpI0aMsPT1AgCAxkUlrCixVlf/jz76qH379mvWrBFCeHl5LV26tMYoc+HCha1bt2ZnZ/v6+o4ZMyY0NHTx4sVBQUF79uxxd3eX+gwePHjz5s0vvPBCQUHB6tWrT5w40aFDh+nTpwcFBZ09e9bSGeFvv/326aefTpo0adGiReHh4ceOHRNC5OTktGrVytSnVatW2dnZ3BaMRqOFjwUAAHchm3yNp6WlDRkyRLo9ePDgAwcOaLVa8w779++Pjo729fUVQgQHB4eHhx84cEAIYRoFhRBardbLy0sIcfDgwaCgoA4dOgghvL2977///rS0NEsHwvLycm9v759//nn37t0PP/zwggULhBA6nU46FCtRq9U6ne5OnjAAANy1LBkIpdUnZP2rW25ubosWLaTbAQEBRqMxLy/PvENeXp6pg9Tn6tWr5h22bt2akZExffp0rrOlh0aDg4Pj4uKkCezgwYM/+OADIURQUND169dNffLz84OCgrgtyK41AAAAdxMnJwvmTlZdR3jo0KHExETztkGDBs2aNUsIodFo9Hq91Cjd0Gg05j3VanV1dbXpR51OZ97hxx9/nDdv3vbt2318fLjOlg6E/fv3v3DhgnT7woULwcHBQog+ffrMmTOnqKjIx8cnNzc3IyOjZ8+elj5zAAAAIYQQwcHBEydONG/p2rWr6b9MM7ycnByNRuPv71/jvjk5OaYfr169GhISIt0+cODAuHHjNm/e/MADD9TemrTB+Ph4SwfCp59+ul+/fs2aNWvSpMmrr776zjvvCCHat28/evToMWPGTJkyZcOGDdOmTZMGSIDGgzoUJPc8Cf33MQ6RQGMkHe+UeRcREhJSYyA0SUhIWLp06YsvvqjRaLZt2zZ8+HDplNyvv/4aEBDQsmXLwYMHz5gx49y5c+Hh4cePH8/Nze3fv78Q4siRI+PGjVu/fn1cXJxpa7169aqqqvrxxx979+596dKlo0ePxsfHW3qOMDIycv/+/eXl5VevXk1OTn7wwQel9n//+9+TJ08+ffr0zJkzV69eLevJAwBAYyP/HGHd583Gjx/fpEmT/v37T5s27Z///OfixYul9rlz527btk0I4evr+/zzzw8ePHjmzJkjRox4+eWXvb29hRAPPfSQk5PTihUrBg8ePHjw4LffflsI4e7uvmTJkgkTJsycOXPAgAFPPfVUcHCwqmFSQBUVFfP+vnzkI7Ma4LEAlIQZIcAfxoWzKRCT3v02nzp1Q9Zm9doTQwbdSE5O5jrodLq9e/cWFhbGxcWZjouePn3a19fXdJn8kSNHzp4926VLl+joaKnl5MmT5vlSPz+/Nm3amO577NixsLCw7t27C9QaBQCAu5xGoxk6dGiNxhpFPe+555577rnHvMV0lrG2zp07m98dAyEAACjGigvqbb4eIQZCsBsyj+IzveU1y8RtnPo1V3EPyX4nKPBlYesvHGj8LLk0sPZdbAsL8wIAgEPDjBAAAJRjxQX1toaBEAAAFKNSKVx0uwFgIAQAAMXY4YQQ5wgBAMCxYUYIDUJWVpOJh3IpS7I7WymCa5bTX/YF9WRqVPbhI+Zhqe2w22ajqjJ2xtbHseDuZlXRbdvCQAgAAMqx4hxhPe2JxXBoFAAAHBpmhAAAoBiV/IPnNj/YjoEQAAAUo1LJvxzC1iMhBkKwlqwUCZd/oZoNBqYz0072NxroHTFwe8I9qCJ13aiTINxa39yXiMqJaaeanbjO7EaIPec6I3EDdbHDsAzOEQIAgEPDjBAAABRjhxNCDIQAAKAca0qs2XooxKFRAABwaJgRAgCAcuzw2CgGQrgNNjZJ/QfXmQuCku3VejrxWa1nNlJN9K+ulvGIQggD059Mk8qMkjIl1phMJhf4ZNvVRLuzM32wx9mZ2QjVLjd6SuZguYNkRubLD2lSe2fFdYQ2HwhxaBQAABwaZoQAAKAcO6w1ioEQAACUg3OEAADgyFQq+6s1inOEAADg0DAjhN+xC9PKCXxymUwuCKrXEe1koxBCp+U2Uk09orw9YSOpdGpUXmyUPGXCpkCZYKczlQ4VQqjVxJ+zag2TGmXayY04U4117AkZVeWeDldqlf/jXM76w2A7KqGSe4G8zS+ox0AIAADKseIcoa3h0CgAADg0zAgBAEAxVoRlbD6DxEAIAACKUQn51xHa+mQvBsLGjI50MEEPJuZCFzATTLqEy7loq4g8ixBCR7VznWVtRKelO3N7yIVouNJrNOb3mQyGODFxETaKwkZXiHaNC91Z4+LMtFOJG64zl7ghNyInniP4cI2K7o46bXcfnCMEAACwLxgIAQBAOarflySU4zabXLduXbt27Xx9fR999NGKioraHY4fP96jRw8fH5/evXufOXNGajxy5MicOXN69eo1ffp08573/tn+/fsxEAIAgGJU8v/V7cSJE88888zWrVuzsrIuX7786quv1uhgMBgmTJgwefLkvLy8kSNHPvjgg1L7jRs32rZte//996enp5s6h4eHb/nDs88+m56e3q1bNwyEAABw99qwYcOECRPuvffepk2bLly4cN26dTU6pKamlpaWPvnkk66urgsWLLhy5covv/wihBgyZMhzzz0XExNj3tnNzS30D3v27ElMTPT29sZACAAAyrFiSljnrDA9PT0qKkq6HRUVlZeXV1JSYt4hIyOjS5cuUvpMo9FEREScO3futrtZUVGxZcuWpKQkgdRo48BWRyPXzmWroNHtfMEzKvBZSWc1Kyv0ZHsV1V7FbERbSW+EfFD5qVEZJdbYX1vmP8iVbLkSa9zaucqkRl3pIKgL1c5FTMnOQggXN6KdbKxjI2xxOOppyq3TZvOYviOw6KRfjbsIlVarLSoqMm/08PBwdXUVQhQWFnp7e0uNTZo0EUIUFBQ0bdrU1LOwsNDLy8v0Y9OmTQsKCm77oFu2bPH39+/du7fAQAgAAAqyYoV6lUqkpKS0b9/evHH27NlLly4VQvj6+t68eVNqlOaCfn5+5j2bN29eWlpq+rGkpKRGB9L69euTkpKkMRsDIQAA2NiIESOSk5PJ/woLCzt16pR0+9SpUwEBAdK80LzDmTNnDAaDk5OTXq9PT08PCwur++GysrJ++umnzZs3Sz/iHCEAAChH6XOEjz766LZt244ePVpaWrp06dIZM2ZI7S+++OI333wjhOjfv7+np+d7772n0+nefvvt4ODg++67TwhRUVGRlZV1/fr1qqqqrKysvLw80zbXrl07bNiw4OBg6UcMhAAAoBj5FxHepjhpdHT0G2+8MXbs2NatWwcHBy9atEhqv3DhQmFhoRDCyclp27Zt//73v/38/JKTkz/77DOpw/HjxxMTEz/55BONRpOYmPj666+btnnq1Km5c+eafsShUQAAuKs99thjjz32WI3Gjz/+2HS7W7duhw4dqtHh/vvvP3z4MLnBHTt2mP+IgdCeyEqHCqZOJldUk1v2lstqktnOynK6M9+uI7ZcITM1StUa5dKhXO1QbklhZmFesm8d/0Ego6TCigV7qeVwFUmNsulQd/p7w82d6O+qpTtXe9DtLm50u8aFeG25iKkTc6zLyYnKACNKqigrVp+w+RuAQ6MAAODQMCMEAADlWHEdoa2nhArMCLVa7YULF3Q64hgXAADAXc7SgXDBggXmIZ/y8nKpfceOHSEhIcOGDWvVqtXu3bvrbT8BAMAOSOcIZf2z+UlCGYdGX3jhhSVLlpi3VFVVzZgxY9OmTcOHD5eKtl24cMHZmT7BDrKQqQtZoRjBBEa45W3JamdCZs6loozpXEYfMCA3zpdY46qmEe3VzGvCVZhjwzLUdrjOMkM0slaglReW4VIkXP6FDNFwnV2Z1JKuSkM0KpRaMhqoLyvmBVe70O0q8k9/th4b3Q63Z28vnbxDozUWgtq1a1eTJk2GDx8uhJgwYYJWq923b5+SewcAAFDPZAyEb775ZosWLQIDA9966y2p5cKFC6ZKNk5OTh06dLh48aLiuwgAAPbCigvqbX4Fi6UD4Zw5c4qKim7duvXFF1+8+uqr27dvF0KUlpa6u7ub+nh4eJhKo9ZmMNBHSAAAwC5UV9OnJ8xZcY7Q1uOgxQNhhw4dpDHv/vvvnzZt2rfffiuE8Pf3N184o6ioKCAggH0kdsUUAACwA5ZEQKwYB21+UtGawamoqMjT01MIER0dffToUa1WK4QoKys7ffp0dHS0wjsIAABQnyxNjS5evPiBBx7w9vZOTU3dunXrDz/8IITo0aNHWFjYU089NXv27BUrVvTs2TMiIqI+97YRklU1TVY6VDA1ybgUKBf4rGACn+W3iPbyUi4dyrRTUVUuHcqttUuuqcu9VlxqVNafpFyZOnZ1X+rd5B6QO1/CpkapFWu51X21zFq7slKj2iqunXhZZL1roq7UKNnM4FZIpv6DKWnHbsXmx/Hudraf4Mlm6YzQaDQuW7Zs/vz5J06cSE1NjY2Nldq/+OKLioqKxx57TKPRmGp+AwCAg7LDc4SWzgj//ve/k+3BwcHr169Xbn8AAAAaFGqNAgCAYlTyF/SwmxkhAADA7dnhOUIMhA1EkaUEuVBMFVPyiszFcHmW8lIt3U6FYoQQFWRYhunMJXTIPdcxReC4FRPJ14p7YdVMWsTTm67K5Uqtscc9nZtFVWS7llphkdtDbp1CLv9Crl/IJWu46IpGRzxNPdUohNBx7dTGyQJ4gg/LcGkmOVXq2BkJ2a5mitrxVwpQG7e3r/76Y816hLZ+9XBtHwAAODTMCAEAQDHWlEyz9ZQQM0IAAHBoGAgBAMCh4dAoAAAoxpqwTP3sieUwENYDMt7GpUOZjJy8NXWZECMZEC27SadDuXY2NUqVXqvkFuZlVv0lc4ZcOlQvp1gXtzJtMz86HRoa6UO2N/f3qN2Yf7WM7PzbyQKyvZx6bbXM03GSWWKNTI1y0VPuZSFfQy7YyVZHIxcxZmvdkc2CXWyXwn3hkq+JYF5DrjP3GjqpyIJ5Nv8yv1uohOxzhHazDBMAAECjhBkhAAAox4oL6m09ncZACAAAyrkLimjLhYEQAAAUY8V1hDhHCAAAYEuYESqPDL1xy8EamAAenRplaopya+qS5UPLmJqiXHvFLa6QKZUaZfKr7Fq71NPkqlPqqHVfhRDV1VS7iv5su3nQ7YFtvMn24LZNqI2TfcXljGKynXyaXNaX4+TMlMRUE3tDrtYr+AKnsjrzG5FR9JXDzQ3IDCe7WDH3WpGpUWYjXGpU1vTF1lMdsAgGQgAAUIxVh0braV8shUOjAADg0DAjBAAAxVhRWea2l0/odLrvv/++uLh4wIABLVq0IPscPXr03LlznTt37tq1q/kdz507p1arO3XqZGo8fvy4Xv/7iQkfH5/Q0FAMhAAAoChFD3Vqtdq4uDidThcaGjpv3ryUlJRu3brV6PPaa6+9//77gwcPfvrpp5999tn58+cLITZt2jR79mw3N7ewsLCDBw+aOsfFxYWHh3t4eAgh+vXrt2jRIgyEd0DOWrtGptCUXk8HQMjaY1XU+q6CqXYmmBJr3MK8fHU0GWvtsmEZpjgcmQmSG5YhX0OusldJYSXZnp9DV00jf6dvMCXWSovphXnJYFEVU3aOy5ZwlcA01FLD7p50Zzd3GRkiLuFVzSwQXUk1G7hSaszG2cJmVKSFywSpuXaqwpwzU3aOXwmZ6K+StYpvo6ZSyb4cou7+n3/+eXFx8ZEjRzQazSuvvPLKK68kJyebdygsLHz11VcPHz4cERFx9OjR/v37z5gxw9vbe8iQITk5OTt37ly5cmWNbW7cuLFjx46mH3GOEAAA7l5ffvnlmDFjNBqNEGLixIlff/11dfWf/lZOSUnp0KFDRESEECImJiYgICA1NVUIERQU5OvrS27z+PHjP/zwQ1FRkfQjBkIAAFCO6vfThLL+1SE7O7tly5bS7VatWul0uvz8fPMOOTk5pg5CiJYtW+bk5NSxwSZNmqxaterZZ59t06bNxo0bBcIyAACgJCtqjQqRkZHx+uuvm7fExsYOGjRICKHX652dfz/0LN3Q6f50rkGn05k6CCHUanWNDjWcO3dOml9+880348aNGz58OAZCAACwMZ1OZzpQKdFqfy/xERQUdP36del2fn6+SqUKDAw07xkUFGQ+R8zPzw8ODq7jsaRRUAgxbNgwb2/v06dPYyAEAADFSAc75d6pc+fOy5YtI/+vb9++X3zxxfPPPy+E2L17d8+ePV1cXIQQ5eXlGo1Go9H06dNn9uzZhYWFzZs3z83NTU9P79mzpyWPeuXKlcLCwpYtW2IgtB5bZor6D7oMGBObFExNMm55WzY1Sq2pW8l05kKMWiaqWkXtIZcOZUus0Qvzcp2ZEmtUalTPdL52+RazEfr99GxCLORbWiQvekomcrm3nitJxi7MS7W7e2nIzq06NiPbfQPcazfeZEKw3GtIft64pac53NMkg6BqF/pD60IlaYUQGleiXeMiL3pK7iEXduXYvJZK/bFmhfo6+0+fPv2NN954/PHHw8PDlyxZsmHDBql9wIABkydPnj9/fmho6Pjx40ePHj1lypSNGzdOnz5dmhFmZmauW7fu9OnT2dnZzz33XKdOnR555JG9e/d++OGH3bt3r6ioWLdu3UMPPdSxY0cMhAAAoChFL6hv1qzZwYMHN2zYkJ2d/cUXX/Tp00dq/9vf/hYWFibd3rhx44YNG86ePTtnzpxp06ZJje7u7qGhoaGhoQkJCUII6YBqdHR0//79MzMz3d3dV65cOXz4cIGwDAAA3OWCgoIWLlxYo3HChAmm22q1eubMmTU6tGzZctasWTUafX19Z8+eXaMRAyEAACjGmqLb9bQrFsNACAAAiqmPWqP1DRfUAwCAQ8OM0AJyaooKIYxUTK6aqzXKJCTp1CgT+Kyg0nqCSZmywU6Z7WS2k6vwyb1WzPrD8lKjZC5RJ5g4LvOCc2VCyYggtydc8pZ8mmw6lPnr1NnA/IcbkYT0aupK9g2N9CHb23Yi2nMvlZKduXK1ZTeJtZ11zOeHmzRwJVWd1cRrywU+K6l0qBDChXqtXKlGwURMBVOb1MmJfjtVTAi2MbPignpbv0gYCAEAQDnWnCO08UiIgRAAABRjhxNCnCMEAADHhhkhAAAoxorLJ2xeaAcD4e2xFaKY/yADI2QZMCGEjqmzVVUpY9lbrvQamVPggh6yCpgJIVyoKAFZkEzw+RdDdUXtRu5pcnsia2FeuXEespnfOBfnIRq5331usVmNK7e2M/FGcK8VF9oi27nqaNzTJ98IrnYd94LzYRniZdG40J8rF+rXRwihrSS+8bg4mAvzG6F2Idqdnel3TcWFaGz91V+/7O3YKA6NAgCAQ8OMEAAAFGNN0e362RPLYSAEAADFqFSyD/za/DgxBkIAAFCOHV4/gXOEAADg0DAjtIS82BuZtWNTo3JWsiXXwq1jI7KWvdXr6HY1VVBKCOHf0qt2Y2Brb7IzuWSrEOK3Ezcs71xRRr/g9NNkXiu9zBAj+eYbmY+EgX5M+iPB/RHMveBGI1Pxy5lISN64Si8RfOZQPtmen02stXuziC47l59DL8zLxX1lcVZzn0Ni41VMiTWNK93u6k4WHaS/BrnfiGod8UYY1PQvuJMz/a5xHzebHyS8c3Y4IcRACAAACrLD6whxaBQAABwaZoQAAKAomx/rlEnejFCn0yUkJDz00EOmlrNnzw4cODAoKCg+Pj4rK0vp3QMAAHsiXUco959tyZsRvvbaa5cvX3b6Y8E0o9E4bty4qVOnbt++/c0335w0adKhQ4fqYSdtjV2PkG5nwjJc0INZNo+sjsaEYsjl7rh27hG5jXBV04LaELmYTrEtyM63iomV6oQQRdeJEmvcMnhcxS/yGXHBomqmChq7EAzZzAZruCUqiUYDV9SNL+pnueu55WR7FfMRcjtHfBXwH1p6I07U95maybNwpdS4zyH5LrOLaDLvPtnOZ81klPTTMKtFsiEslq3HhDtmRa1RmxeckzEjPHny5I4dO55++mlTy759+woKCp577jkfH5+XXnrp3LlzR44cqYedBAAAqC+WDoR6vX7mzJkrV67UaDSmxrNnz0ZHRzs7OwshXF1dIyIizpw5Uy+7CQAAdkEl/5+tWToQvvnmmz179uzVq5d5Y0FBgbf3/46PNWvW7MYN4rIwiYGswA8AAHaiupq5VNaMNScIbT0WWjQQXrx48b333nv00UezsrLy8/OrqqqysrKMRqOPj8+tW/+7tLakpMTX15d9JCdcqgEAYMecmfoA9s6isMyNGzcCAwOTkpKEEIWFhdeuXUtMTNy/f3+HDh3Onj1rNBpVKpVer8/MzGzfvn097zAAANy9pGmezLvYmEUD4b333nv48GHp9ieffPLGG29IPw4cOFClUm3YsCEpKWnVqlUtWrTo2bNnPe5sPZNTY0vqL2PZUjJpJuoKfBKHINiNMOk+JjUqL3rKRQfJwF55KV0djWvXUguosq8Jl+Kj2rl0KBngFEIIlZzAJ5Nf5cjKxHHL3uoEs7gx9Tnk9pArG6ZxIf7S53ZbxQQ+ybWahYr+kuHWH2aXX6b2XKelN8KvPk1shIuecmnSaurzyS74zHzejMwBMvqVtflAIYsd1liTfbjSy8srJCREuq1Wqz/77LM33njD09Nz7dq1mzdvtnkKFgAAQBbZlWVGjRo1atQo0489e/Y0HR1VdMcAAMAeyb+O0NZTQmVKrGEUBAAAYdUK9bYeB1FrFAAAFOQI5wgBAAAaE8wILcClSZk8GBnY40JlXEaOzGoamAAnH7Qj2rnsJRe0u1VCr8564UxR7caifKJ2qODX2r3yWwn1iHRhUi7dR4Zp2foN7HLKdHfyteUymdw5ArKZ7cz8dcrWttUTL4uR+bxxyM+nkzO9h87OMsqHqjXy8tLOzIOS/ckAZx3t5IecjZhy6WXyI8GlRtkwOv00yWZbz5fksabW6O2e4k8//fTmm2+WlJSMGTNm3rx5tbd//fr1l156KT09vUuXLn//+999fHyEEDk5Od98882JEyeCgoKef/55qWdRUdG6dev2799fWVl53333Pf30082aNcOMEAAAFGNFhbW6h8Hs7Oz4+PghQ4a8+OKLK1euXL16de0+48ePr6ysXLZsWWFh4eTJk6XGvXv3fvvtt1euXElOTjb1PHbs2MmTJx9++OH/+7//++mnnyZMmCAwIwQAgLvZunXrBg8ePHv2bCHEa6+9tmjRoscff9y8w9GjR48ePbpr1y43N7c1a9a0aNHi3Llz4eHh06ZNmzZt2r///e+VK1eaOg8YMGDAgAHS7datW4eFhZWWlmJGCAAAylF6SnjkyBFTqZaePXump6eXl5fX6BATE+Pm5iaE8PLy6tq169GjRy3Z0/T0dD8/Py8vLwyEAACgGNUfpwktV/dAeO3aNemcnxCiefPmUot5h/z8fFMHqU9eXt5t97OgoOCJJ55YtmyZSsVUPwJzskus0QvzMmf1uRANVSGMKxsmLyzDdOZKr1VVEFXQhBBlN4lIy5Xf6E90FVVKTTC5GC5Zw5UHM1DPiEuLcDkX7t0kK4F5e2lqNwoh1Br6L0vyGVWW068JU+tNcFXrq8lwFhfnYaIoZC5GzaTKVU4yXkNuTwzMh9lQTT8o/WvFvMuyKhqyxQXZooMy8lNGrp1stbdcDMmK6whVQqSkpNQoVf3QQw/94x//EEJ4e3tXVPyewpPmgk2aNDHv6eXlVVlZafqxrKysRofaSkpK4uPjExMTpRraGAgBAMDGevXqVSMFExAQIN1o06ZNVlaWdDsrK8vLy1RU7M4AACAASURBVEuaF5q0bt36/Pnz0m2j0XjhwoXWrVvX8VilpaXDhg3r1avX66+/LrXg0CgAACjHqnOEnp6eoX/m6ekpbe/BBx/csmVLYWGhEGL16tWTJk2SLp/45JNPjhw5IoQYMmRIcXHx7t27hRBfffWVwWDo168ft3fl5eWjRo2KjIxcsWKFqREzQgAAUI4V1xHW2T8uLi4hISEiIsLPz8/Z2fmbb76R2t99993Jkyffc8897u7ua9asefDBB1u3bp2dnb1+/XqNRiOE+P7778ePH6/Vaquqqpo3bz506NDNmzd/9dVXqampx44d+/zzz6XtHDx4EAMhAADcvVQq1apVq15++eXi4uKOHTuaRs0DBw6Ybo8fPz4+Pv7KlSutW7f28PCQGgcMGCDNI80lJiYmJibWaMRACAAAiqmnotv+/v7+/v7mLc7Of1r/0tPTs1OnTjIf+HcYCC3A5AnZ2klUko3PzjHtVDU1LgvHrZ1LPigXY+P2UFtJZzXJgmdc9JRtr6JKXnELpTKlsMjooKxMr+ArewW28qrdGHGvf+1GIYRnExey/fypgtqNv50kGgWfJpWFTSNzHxWyP/P1oGbW1HVxJe6g0VCr9fJHw/iQtqWNQmZFQy4dyv66kb/gbIk1eg/5/yA1gjDpXQ0DIQAAKEZlh+sRIjUKAAAODTNCAABQjDUX1Nv60C8GQgAAUIw1yzDZeiTEoVEAAHBomBFaj80lUu3cKr6ywo1sxFRuYo3uTTdzBU7JCo1cOpRbDZgM5nG1H7l28ulzf2Jy7S5udLixVYemtRv7JrQjO/sFeZDt+74kftGuXblFdq4oo1OjfJVUopH/C1vGRsg6q0KIJj6uZLtfsGftRrb+aqm8irLcjtOt3K8V9VGRW8KX3IjsWqPyFk62K7ddYPDug4EQAACUUz/XEdYrDIQAAKAYlbjdukrEXWwM5wgBAMChYUYIAADKseIcoa2nhBgI7wBbPIlqY+u0yTiXzlZH4zYiZ8lWrp17muSes2XnuHYl9pBsd3JiKnixe0K3k/W05C6zTG5c1msihDCy8R8Z3yJcZ/Ll8vCm1x9uE+5Dtkd2JyrPcW9E1pma1ZAluZdKyXZlPm9ylg5mA25k5UKZhRhZZH9bjxOyqORfF2jz54dDowAA4NAwIwQAAMWoVPIvkLf1BfUYCAEAQDk4RwgAAI7Mmlqj9bMnlsM5QgAAcGiYEZqzRdUjmbE3mRuREbRjnz3z1xq5ki1XlIsrAqeiIoVOzBq5XP5QkZJXVZV0YbPLGcW1G79PziI7ezahY5ZZZ4pqN5YUVNJ7yLxD7Eq2VH+uM/enOvluenrT6wyHtGtCtne6pwW1f/QjlhZXke0FeeVkO/0Gyfww0y+tEr/3bE7VJt8qNoWi2wAAAHYGAyEAADg0HBoFAAAloeg2AAA4LmvOEdbTrlgMA+Ht2eRkN/lBkpuAIPvLXaiPW0/OxZVYwI9Yj04IIURFOR1FuUUlJrgCZtye0GEZZpE57neOi/NcyyZWDSxmci7ca0guMcit3ejszJywkBOJEkywiH3+VDP3Rty6qSXbr18tq93IpZPIt17wq07K+j3kIip0pEXmb4Ssr22bf8XbgB1eR4hzhAAA4NAwIwQAAMXY4+UTGAgBAEAxVlSWsTkcGgUAAIeGgRAAABwaDo3eATlhMz7wKScIyj0iExEka5KxBcyY6mie7vSHpEUIERFt6utGdubKiV06RxQwI4OaQgi9ls4T6pyI+KWRjqny+TSmnQwxVlXQiUe+OhrRyJWj4w4ryVuXmIj01oUMzXJV0M6fKiDby0uJNCn3NAvzmVJqFfQ7p3YhtsNGoJm3U97vJvdrRYexlciYyu9+F7LHc4SYEQIAgEPDQAgAAEqS8jKy/t3W9evXz58/zx10EUKUlZVlZmZWVFRYsoeVlZWZmZm3bv1+8AkDIQAAKEdl1T+e0Wj8y1/+EhkZOXz48G7dul29erV2n88//7xNmzYTJkxo06bNzp07pcZvv/22R48ebm5uPXr0MO+8Z88eU+fNmzcLywfCXbt2zZ07d9y4cbNmzfrhhx9M7SUlJc8///y4ceMWL15cVkbUlQAAAMehkq/uDe7du/eLL75IT08/d+5cbGzsyy+/XKNDZWXl7NmzP/nkk+PHj3/44YczZ87U6/VCiKCgoCVLlixfvty8s8FgmDVr1ttvv338+PH//ve/c+bMuXXrlqUDYVZWVlRU1KOPPhoWFhYfH79//36pfeLEiZmZmbNmzfr1118ffvhhC7cGAABgiU8//TQxMdHX11cI8fjjj3/66ac1DpDu2rWrSZMmQ4YMEUKMGjVKCJGWliaEiI6OHjJkiI+Pj3nnX375paioaNKkSUKIPn36tG3b9quvvrI0NTpnzhzpRkJCwqFDh1JSUh544IGTJ08eOHAgPz/fw8PjvvvuCwoKysrKCg0NvdPnfZdh02BygqBsvI35U4ResZaJsXEbJwOizk70Q6rV9MH3Js1dyfawaL/aje27+JKd83PoIKiOCoLeLKLDipVMwVKyOKde0BFTdolgObE1tYbLGdKvLffGkbh0KFdsk84Gy/yokC9L+S0d2TnvcinZXkbVIHVxY75kmJM9XKpZ40IFYdnAJ/2Y8nLUTDv5u8n+IisUJrUjVlxQX3f/CxcujB07VrrdoUOH0tLSgoICP7//ff9cunSpQ4cOf2xK1b59+4sXL3Jbu3TpUtu2bdVqtWmDFy9elHf5RHV19enTpw8ePDhz5kwhxKFDh2JjYz08PIQQPj4+kZGRhw8fbnwDIQAAWMqqotv5+fm7d+82b2vXrl379u2FELdu3XJ3d5capeGmtLTUfCA07yD1KS2l/1aT7lu7s4yBcM2aNQsWLLh169aTTz4ZFxcnhMjLy2vevLmpg5+fX15eHnd3g4EpKg8AAPbAYDA4MYc97lBmZubrr79u3jJmzJh58+YJIQICAoqLf7/guLCwUGox7+nv719UVGT6sbCwsEYHcwEBATU6d+/eXcZTmj17dmlpaVZWVmpq6ttvvy2E8PT0rKr631GsiooKT09uHR5RTy8fAAA0DEu+xlVCyL96QtW7d++UP5NGQSFETEzMTz/9JN3++eefw8PDpXmhSbdu3Y4ePSoNRmVlZadOnYqJieF2Lyoq6vz589KAWl1dfejQoZiYGNmDU7t27aZOnSrNYVu2bGl+KPbSpUutWrWSu0EAAGg0VFZcR1jnBpOSknbt2rV27doffvhh4cKFTzzxhNSemJj40UcfCSFiY2OjoqLmzZt36NChxx9/vE+fPp06dRJCXLt2bevWrQcPHiwsLNy6dat0vUPbtm2HDRs2e/bsQ4cOPfnkky1btuzTp4+lh0avXLkiDXJVVVW7du3q2rWrEGLo0KFJSUk///zz/fffv3v37srKyn79+ln98t0FZKxZKurKv5Dn0mWWcZJzVp9bytWZ6u+kZh5Rz2ycKZHl4kaEFzy8NWRnd0+6XUPVzSJ3u6526hlxrxVbqIwL0ciKonChGLKZCb9wISy1i4xMh5qr38bsIZlcMTJnM3RV9H+Q4RruBafDL/znTVYmiK8jSMXHuM5yNi73I8F/9dt/isa6hXn5hZdbt269c+fOt9566z//+c/jjz8+d+5cqT0qKiooKEi6/fnnn7/00kvPPPNMVFTUihUrpMaCggJpzhYXF7d79+6uXbv26dNHCLFp06bFixc/88wzYWFhX3/9tUqlsnQg7N+/v5ubm6+vb3p6elRU1IsvviiE8PLyeuedd0aOHBkdHX38+PF//etfrq50whAAAMA6vXv37t27d43GRYsWmW77+/uvXr26RofIyMg1a9bU3lrTpk1Ng6XE0oEwIyMjIyOjuLi4VatWLVu2NLUnJSWNHDny/PnzYWFh0nUeAADguOQX3RYqVR0zwgZg6UDo7OwcERFB/pe/v7+/v79yuwQAAPbKmusI62dPLIckJwAAODSsRwgAAIqxaj3CetoXS2EgvD22RpKcUBmfhGQCnxoqTskF6rggqJw94T6L5aV0na2L6UW1G7miXMU36IV5r18lCrVrq4iFdutAvoZq6gUUgg2nGdg1dWUUzOOQ+XBZ75qoI8SoRK6VvDyMjOMKIZyZCnOyKslxndkHpdq5Wnfcu0+GaZ2Y30FZpdfk1rSz/Xc/mMFACAAAirHiHKHNTxLiHCEAADg0zAgBAEAxKiH/HKGtp4QYCAEAQDnWVZaxKRwaBQAAh4YZoQXYWqNysppM4JMstsm1a5gsHBc9petkygwlVpTTQdBL54prN+ZcoJcBqyyjN1J0vaJ2o04rLzVKPk0uN8i9a1xdC7I3twqriytdP9O7GVF60MPLheysZZ5+WQmx7K2Q+XJxqVEypcyVA9UwT1NNfWi5N4J8RMH/ppDtXGe1htlz6hnJipgKroSvnFV8RaMOjSq+MG8DwEAIAADKsaLEmq2PjWIgBAAA5eAcIQAAgH3BjBAAABRjhxNCDIRmZL95zB2c5WQ3uHYX6qw+uRauEMLFld6ItpJo13NF3bilXJksRkWZnth4cRW9J0zVNB3VLiuFJJjwgtHIBYiYWAy7TC7RyKVFAlt7k+2R9xLLs/i39CI75168SbafOZxPtt/ILa/dyC2Hy33eyBSJhvlckR9OIYSLK/F94sp8aLnXkF+wV0ZEhcu50Bk0JrDGxnzoOm3yMmhs5TWm3Y6o+F9h9i62fto4NAoAAA4NM0IAAFCOHR4bxUAIAACKsabotq1hIAQAAMVYtR6hjUdOnCMEAACHhhmhGfaPEq54Et1bZok1JoBHZe1c3ej3y4Vp11QQmUyd1kB25kpecdWqqvVELlFPV1JjOVEvi4aoRyYEHy2j185lUnwGPf30DXSzIOOkXCm1FsGeZHvXnoG1G0Mjm5Odz/5Kp0Pzrtwi228WEUndauZpsp836hmRKVAhhLunhmz3bEoUjXNzpzfC1iRj3mU1mWtlS6kxv27U0+ReEzY1qszCvDLb7QjOEQIAgGNTyV1WyebLMOHQKAAAODTMCAEAQDGqu+ACebkwEAIAgHJwjtChcH/1kKfN2dJWTOyCzL+4uBNVzYQQru70RqoqiHZu+bpqPb2HBqbOlsFAbLy6mq7sxUVRyJyLmiqmJYTQsyvVEc9Ir6Mfkms3VHMhGuIZcYvMca9t2U1iKcFSphwdt3ajkamaRgZDZJVSE0w4y9ObXjHRL8iDbPcPIYrGuXnQXzKVzDqXZOk+Dlt0kC1GSFaSkxmWIUussWGZRpyKoeHyCQAAADuDGSEAACjGmhXq62dPLIeBEAAAlGOH5whxaBQAABwaZoQAAKAYK8Iyt50R5uXlbdiwoaSkZNSoUb169ardQa/Xb9q0KT09vUuXLlOnTnV2/j0AlZGR8Z///MdgMEyePDkyMlLa1I4dO8zvO3jwYAyEtyf3PWVKrNFb4WpBkfE2rsQaV8Wqyp0ssUYnG7mlXI3cQrZUytLApEa519BQTT19pjMXmyRX/a2qlLEUsOBrktHPiNnDwmsVZPuxH3JrN146VyRrIyWFlWS7E7XSMv95oxOSZLbTp4Ub2blNWDOyPbQLUTSOq9OWn01XjLt+tYxsJ99lfolgGav+1uvCvDY/6NfwFD8yWlJS0qNHj6FDh3bq1CkhIeGjjz4aMWJEjT5JSUm//fbblClT3nvvvZ9++mn16tVCiMzMzB49esybN8/FxaVnz54//vhjly5dqqqqsrKypHsVFBSsXbv2zJkzGAgBAODu9dFHH4WGhn744YdCCG9v76VLl9YYCC9evLh169YrV674+fmNHTs2NDR08eLFQUFB//znPydNmvTqq68KIW7evPnOO++sW7euTZs2y5Ytk+74zjvv9OzZMyIiAucIAQBAUSqZ/+q0b9++IUOGSLcHDx584MABrVZbo0N0dLSfn58QIjg4ODw8/MCBA0KItLQ00x2HDBmSmppaY8sbN26cMWOGQFgGAACUpLJGHdu7evWqv7+/dDsgIMBoNObl5Zl3yMvLa9GihenHgICAq1evCiFyc3NNd/T398/N/dMZil9++eX8+fOJiYkCYRkAAFCQNdcRqsShQ4ekMclk0KBBs2bNEkJoNBq9/vd6Q9INjeZPC4Gp1erq6v+dQtbpdFIHtVptfsca91q/fn1iYqK3t7fAQAgAADYXHBw8ceJE85auXbua/kua4QkhcnJy1Gq1aZ5Xu4MQ4urVqyEhIUKIkJAQ8zsGBweb+lRUVGzZsuXLL7+UfsRAeCdkLMUpe2FeVyLE6MqkQ9086IVSyaAdV2yTTIEKPjVqoOqHGqrpp8Nl6shwo4cX/XS49YerKojqlCUFdMaSXMZWCKFlUqZkmpRL2HL1My+mEwHRnCz6IyG3Sipda5T5vHFJSFey1mgTutaobyBdazSwlbflj0i+a4KpyyoUKuErb2Fepp380HIH9xyu0qiwKjYqREhISI2B0CQhIWHp0qUvvviiRqPZtm3b8OHDpasjjh496u/vHxISMnjw4KSkpIyMjLCwsBMnTuTm5vbv318IMXLkyK1bt06ZMkUIsW3btoSEBNM2t2zZ4u/v37t3b+lHDIQAAKAYlUp+Ee06+48fP/69997r379/+/btd+7cuWvXLql9zpw5kydPnj9/vq+v78KFCwcNGhQfH79z587FixdLBzznzp3bq1evkSNHuri4/Prrrz/99JNpm+vWrZsxY4ZpPzEQAgDA3cvFxSU1NXXPnj2FhYVvvPFGQECA1L5u3TopKSqEWLRo0bBhw86ePTtv3rzo6Gip0d/f/9ixYykpKQaDYePGjU2aNDFtc8WKFWFhYaYfMRACAIBirFiY97b9NRpNfHx8jcYuXbqY/xgbGxsbG1ujj5eX19ixY2tv8J577jH/EQMhAAAoRiVUKpknCW1+xhQDofVkrbjpzKRF2LP9VHjBVUe/X1zVNJ2WSJ2wK9OyJdaYtXap2mNc4kZjpJ8mGccIbtekdqMQIqAlse6rEKKCWsn2wlm6gFn2byVkO1fATKcl3jiu2Jus14qL53Cr/nKLzZK5GG4j/OeQ2Di3Ea583c0i4jXkYmJcWIb7tdJQC0RzOReuGCG5hDUZFBJ86TVnqqYdtzAv+3xs/t1ff7D6BAAAgH3BjBAAAJQj/4J6m8NACAAAirFiGSaVSsVcrtxALB0ICwsLU1JSLl++3Lp169GjR7u5/b4+i9Fo/PLLL8+ePRsdHV071QMAAHCXs/QcYWxs7ObNmwsKClatWtWtW7fi4mKp/YknnnjhhRfKy8vnz5+/cOHCettPAACwB3KXnrCqEo2yLJ0RHjp0SLp0sbq6ukuXLp9//nlSUlJOTs66det+++23kJCQadOmRUVFPfPMM76+vvW5w/aAOiygYv7kkLVgL5l5E0JU6+iaZGRAlFuBliulxiH7c0+T493MtXZjYGuiUpcQIjSSWPdVCFFJpUYryuhQYvENOh1ayYYYiYQklw6VdUSIjRPKDHyS5eu4PeE+b2R6mfuoFOXTSwc7Ue8+mUcV/AvOpZfJgCi3JLUrtc6wYIoUcqv4coluul6g44VDOdYU3RbCtodGLf3eMl3A7+zs7OLiIpV6+/7777t27SqVN+3QoUPbtm33799fTzsKAAB3P9XvQ6GMfzb/g0H25RNbtmy5cePG6NGjhRBXr14NCgoy/VdgYKB5CfAauD+iAQDALjTWr3F5qdG0tLS//OUv//3vf5s1ayakqI/Z62I0GmXXWgUAgEbG3i6olzEQHjhwIDEx8dNPP+3Zs6fUEhQUZL7mb15envmCTzVgjAQAsGuWfI1btzCvbVl6aPTIkSPjxo1bv379wIEDTY0DBw48depUdna2ECIjI+Py5ct9+/atl90EAACoH5bOCIcPH+7p6blp06ZNmzYJIcaOHTt58uTg4ODZs2cPGTJk7NixW7ZseeaZZ3x8fOpzb+2DrBU6nai6hUIINZUDdXGlD9BXe9DtenJRWT1TJ5NslYlbgJcrzkn211bKW7KVXH+Ye0QuCshV8qw/7B/BzAeIXdyYaucqfMp6+tweksVdhRCF1yzdsuCfDpcyJcuHunlyS1UzqVFqZ7hVfLnXkIz1slMlW891Gp41F9TX065YzNKBcNWqVdXV//u6iYiIkG6sWLHiu+++O3369Jo1a8wniwAA4IjssOi2pQPhuHHjuP8aOnTo0KFDFdofAACwZ1bUGrX1QIjVJwAAwKGh6DYAAChGJf8aAbkL+SoOA2ED4d5nbjlPMgHBLUNKrvsqhKjWE+8v19nAXCoraxFabnlSWasB38gtJztz68FWUxsvvkGXAeOuCObqbJGrsHK4jZOvIXdpsuzqaFQ1Pu6jomHCMmREhVwLV/Dvsqy8FReK4aqmuVO5GHdPurigOxOWcaESN/JKqTG/s7Y+tgd3BIdGAQDAoWFGCAAAirFuPcJ62hkLYSAEAADFWFFZxuZHljEQAgCAcmy/mIRsOEcIAAAODTPChsIX1CJbnahIndpI/+FiZBKPZCZTVgq0jnYSVzFOR1VBE8zSr+RCu0KIynK69BoZSa0spzfCJR7dPOj8If1yMe8m99rSi81yqVFmD7lwIxkQ5QKfaqZsGFlOTG45OnLZW7I0muBXmeaCoGTVNK4zuSdCCBfqteJKqXEfFfrAn73NgeqPVecI62lfLIWBEAAAlNOIV58AAABolDAQAgCAQ8OhUQAAUIxKWLEME64jdGyyPjD8Em70zN7VSLy/crIvdSE/607OdJ5Fy5QH02mJnAtXj02vpRM3ZJyHS0BwFbw0LtzrIuP14l5bI/mE5FTXEzKXGOQ+Elw7uXENVbxNMIXKBLPaH5es4d4IV6Y6GtmfDcUwD+pMPX0uFMOuMGjrE1p3OXu8jhCHRgEAwKFhIAQAAOWorPpXp/Xr17dr187Pzy8pKamigqinf/z48fvuu6958+a9e/c+e/asqf2VV14JDg4OCgpauHCh6dBRXFzcvX94/vnnBQZCAABQkMoqdWzwxIkTTz/99JYtW86fP3/x4sVXX321RgeDwTBhwoRJkybl5uaOGDFi0qRJUvv27ds3bNhw4MCBw4cPb9++/eOPP5bajx8/vmjRojVr1qxZs2bmzJkCAyEAANzNNmzYMGHChO7duzdt2vSFF15Yv359jQ5paWmlpaXz5893dXX929/+duXKlYMHDwoh1q9fP3fu3LZt24aEhDz55JPr1q0z3aVz586xsbGxsbGhoaECAyEAAChICsvI+1fnBtPT06OioqTbUVFRubm5JSUl5h3OnTvXpUsXJycnIYRGo4mIiEhPT5fuGB0dbbrjuXPnTHcZPXp0586dH3300ezsbIHU6F2L/mw40alEJy5PSDW6MW86d3CCy9SREVZZyUbBlF7TMelQHZeEpFKmZPE2IUQ1sy6xvMpzXDqUbmZqcjEvLJca5crXkf3ZN4JpJwOiGq6UGpMa1VCl/rgVj8mIqagjCEpth1t/mEyHCua1UjHTAZtn+u2Y/NSoVqstKioyb/Pw8HB1dRVCFBYWent7S41NmjQRQhQUFDRt2tTUs7Cw0MvLy/Rj06ZNCwoKat9RahRCrFq1qlu3blVVVa+99trQoUOPHj2KgRAAAGwsJSWlffv25i2zZs1atmyZEMLX1/fmzZtSozQX9PPzM+/ZvHnz0tJS048lJSVSh+bNm5vf0XSvxMRE6camTZuaN29+4sQJDIQAAKAkKybTI0aMSE5OJv8rLCzs1KlT0u1Tp04FBARI80LzDmfOnDEYDE5OTnq9Pj09PSwszHTH+Ph46Y4dO3assWW1Wu3s7KzX63GOEAAAFGPNOcI6x81HH31027Ztx44du3Xr1rJly5KSkqT2l1566dtvvxVC9O/f38PD41//+pdOp3vnnXeCg4N79OghhEhKSnr//fezs7OvXbu2cuVK6Y6ZmZnff//9rVu3bty4MX/+fB8fn+joaAyEAACgGKsunqhrJIyOjl6+fPmoUaNatWoVGBi4aNEiqf38+fM3btwQQjg5OW3btm3Tpk2+vr7bt2//7LPPpOsxJkyY8NBDD8XGxnbt2jUhIWHatGlCiLKysmeeeSYoKCgiIuLChQs7d+50d3dXyVptzmoVFRXz/r585COzGuCxGgk6osEtd0dvgwyM6KmqZkKIqkq6OloVswpgObVqYEUZs5EKul1eWIbZc7IkG8IyZGeEZWo3yg7LOHCGZlx40G37rP8uM7+4UtZmzxzaW3o+lTs02gBwjvBuRf6yGZmvSCZNKqhvPX6FYPrDwKZGqfKh7FKuzDqx2kpizNMyq/iyNUjpgZD5o6Ga3ghbJpT6D64zX52SrMsqcyDk1tSl3iB249wbRK/uSw8z3NhGj6ZsZxkbEcxHS+6auuSYhwFPYRZUiiHuYlMYCAEAQDEoug0AAGBnMCMEAADF2ON6hJgRAgCAQ8OM0J6wf2bJCtFwuUGZJdac6bCMvHgFmYsho6RCqbCMrHQol4uRGZYhoyt8WIZbO1dGmpRdw5nZCPnGKbK6L99Z3vrDTOCTW1NX3lq7oCArzhHa/H3BjBAAABwaZoQAAKCc260vSN2lfvbEYhgIAQBAOVZcR2hrODQKAAAODTNCAABQjD2GZTAQNgpciSgqTcrVY1MxhwdUKmY1YCqw5+xMBzs1XGUvN7IaKpMa5cqHUgFRrpQaV5eVLboro9SovMWN5aZGZZVkk71xKsPJdpazVjO3J+xG2CCo5a22/2J1ZCr55whln1NUGgZCAABQlL39IYJzhAAA4NAwIwQAAMXYYWgUAyEAACjHHs8R4tAoAAA4NMwIGzVyYXTuuAXzRxET4uNWXadripKhRCGEmlqFtZopTMoFQclF57kV55lmfrFdWZg/bMnXkKuTyS8qK6O/7Ewm9ZI7cZlMOcFO2WU/EQS1d1iYFwAAHJpK9qFOmx8axUAIAACKseaCTkK1dQAAEBJJREFU+vrZE8vJO0eYl5eXl5dXozE3Nzc1NfXatWvK7RUAAEADsXQg/OCDD4KCgoKDgx977LEa7V27dl22bFnnzp0/+eSTethDAACwKyr5/2zK0kOjffr02bt371dffZWWlmZqLC0tXbBgwe7du3v06LF3795JkyaNHz/e1dW1fnYV6hd/NENGvELFJW4M9H8YqY2o1VzOhdkIlaGRt9CuqKNsGtMuq68SKRJZ/dmTLswbRPZn94RuZj8pzCNyWwH7phIqNpF3t7J0RhgZGRkREeHk9Kf+3333XatWrXr06CGEGDhwoIeHR2pqquK7CAAAUH/uKCxz5cqVtm3bmn5s06bN5cuX73SPAADAbjnc6hOVlZUuLi6mH11dXSsqKrjOBq7sPwAA2AODwVDjuCDhLjjnJ9cdDYSBgYEFBQWmH2/cuBEUFMR1vv3LBwAAdzFLvsYdrsRajx49fv3117KyMiFEYWHh2bNnu3fvrtCOAQAANARLZ4SnT5/+6quv9u/ff/78+ddffz0qKmrYsGGdO3fu27fvlClTkpKS3n///VGjRpmfMoTGQW6alO5KF01jlr3llnKVUx1NkYppSpH39y6XDmX732nnuv+D6GpvB76gIVlxjtDmh1ItnRFqtdqioqIuXbqMHj26qKhImgUKIbZs2RITE/Pxxx/37dt306ZN9bafAABgD+ztIkJh+YwwJiYmJiamdru3t/fixYsV3SUAAICGg1qjAACgIOXDMnq9fu/evTdv3uzXr1+LFi3IPseOHcvIyIiMjOzSpYupsby8fM+ePQaDIS4uzsvLS2osLCz85ZdfKisrY2JipNN5SHICAIBirDgyWvcwqNVqBw4c+OKLL27fvj0yMvLYsWO1+yxdunTkyJHffPPNkCFD3n33Xanxxo0b3bp1e++999auXdu1a9fc3FwhxA8//BAaGrpixYrNmzdHR0evXLlSYEYIypKdrCGbmZyLvL8ybRKWkfl38F3TG0A5Sq9HmJycXFhYeOTIERcXl5dffvmVV15JTk4271BYWLhkyZLDhw9HREQcPXq0f//+M2bM8PLyWrVqVadOnXbs2CGEmDx58rvvvrts2bKOHTueP3/e19dXCPHtt9+OGzduzpw5mBECAMDda8eOHWPHjpWKtyQmJn799dfV1dXmHVJSUtq3bx8RESGEiImJCQgI+P7776U7JiYmSn0SExOlETEgIEAaBYUQrVu31ul0er0eAyEAAChGJV/dG8zOzm7ZsqV0u1WrVjqdLj8/37xDTk5Oq1atTD+2bNkyJydHaje/o9Ro7rXXXktMTHR3d8ehUQAAUIx1tUYzMjJef/1188bY2NhBgwYJIfR6vbPz71ciSze0Wq15T51OZ+oghFCr1TqdrsYd1Wp1jXstW7bsyJEj+/fvFzhHCAAANqfT6YqKisxbTONWUFDQ9evXpdv5+fkqlapGLc+goCDzOWJ+fr7UITAw0HTHa9euBQcHm/q8++6769atS01NlQ6TYiAEAADFqFSya4eqVKrOnTsvW7aM/N++ffv+97//ff7554UQu3fvvv/++6XzhRUVFWq1WqPR9OnTZ/bs2YWFhc2bN8/Ly0tPT+/Vq5d0x927d48bN066Y9++faUNrl279p133klNTQ0JCZFaMBDCXUd+BS9ZRcYAoJ4p+ts3ffr05cuXz507Nzw8/B//+MeGDRuk9v79+0+ePHn+/PmhoaHjxo0bPXr01KlTN2zYMG3aNGny99e//rVHjx4tWrRwcXH54IMPpKOg+/btmzVr1sSJE1evXi1t55lnnsFACAAAd69mzZodPHhww4YNly9fTk5OfuCBB6T2BQsWhIWFSbc3bty4cePGM2fOzJ49e/r06VJjeHj4L7/88vHHH5eVlf3444/ShfYtW7Y0DYESjUajMjZIceKKiop5f18+8pFZDfBYAABQH8aFswvtmXz+65XCMu1tu5k7uj/l6uHvalwd2JAwIwQAAMXY43qEGAgBAEA5SleWaQC4oB4AABwaZoQAAKAY6y6oty0MhAAAoBgrzhHaHA6NAgCAQ8OMEAAAlGOHYRkMhAAAoBiV/HN+Nj+SioEQAAAUoxIqlexlpG08EuIcIQAAODTMCAEAQDlWnCO0NQyEAACgHPnXEdp84MShUQAAcGiYEQIAgGJQdBsAABwbzhECAIAjUwkrZoT1tC+WwjlCAABwaJgRAgCAYuywwhoGQgAAUJAdjoQ4NAoAAA4NM0IAAFCMSmXF5RC4fAIAABoNrFAPAAAODecIAQAA7AtmhAAAoBhrLqi39ZQQAyEAAChGhdUnAAAA7AsGQgAAcGg4NAoAAIqxahmmetoXS2EgBAAAxVhxjtDW4yAOjQIAgGPDjBAAAJRjhxfUYyAEAADlyD9HaPOThDg0CgAADg0zQgAAUIxK/gTP1hNCDIQAAKAcay6fsPVJwgYaCNVqdXrad2lbPmqYh2tgBQUFPj4+Tk6N/DgznmZjcuPGDV9fX/nrxtkZPE1lRXz1VURERN194tr6yd1saEXMT06V1u6UAlRGo7FhHqmwsLC4uLhhHgsAABTXsmVLFxcXW++F8hpuIAQAALgLNfLjPwAAAHXDQAgAAA4NAyEAADg0DIQAAODQcB2hbNeuXTt8+HBOTk5cXFz79u1N7ZcuXdq4cWNZWdnEiRO7d+9uwz1UxJEjR3bt2nX9+vWIiIgpU6a4u7tL7Tdv3vzwww9zcnL69+8/atQo2+7kndu9e/eBAweKi4tbt249bdo0X19fqb24uPjDDz/Mzc0dNGjQ8OHDbbuTCtqyZYurq+vo0aOlH6uqqtatW/fbb7/FxMRMmTLF3i8a2bFjR15ennS7efPmEyZMkG4XFhauXbs2Ly8vPj5+yJAhtttBxeTn52/atOnq1avt2rV7+OGHmzZtKsw+tHFxcSNGjLD1PtoT+/7c28QDDzzw2muv/d///d/hw4dNjXl5ed27dy8uLvb39x88ePC+fftsuId3rri4OCEh4fr1661bt/74448feOCBqqoqIYTBYBg4cOCBAwfat28/f/78d99919Z7eqc+++wzg8EQGhr6ww8/dOvWrbCwUAih1+v79u17+PDh0NDQxx9/fPXq1bbeTWXs2LFj5syZy5cvN7U8+OCD27Zt69ix49tvv71gwQIb7psili9fvmvXrqysrKysrOzsbKlRq9X26dPnxIkT7dq1S0pK2rhxo033UQGZmZlRUVGnT59u27ZtRkaG9Ez1en2/fv0OHz7cvn37efPmvf/++7beTbtiBJmqq6uNRmN0dPSnn35qanzllVfGjBkj3V6+fPnw4cNts3MKqa6urqqqkm5XVFQ0bdp03759RqNx586dbdu21el0RqNxz549wcHB0u1GwGAwtGvXLjk52Wg0Jicnh4WFSW+09JSl23atuLi4c+fOr7zySq9evaSWU6dOeXp63rx502g0nj9/3t3dvaCgwKb7eKd69+69Y8eOGo2bN2/u0qWLwWAwGo3JyckdO3aUbtuv+Pj4hQsX1miUnpr0Qf3mm2/atGmj1+ttsXd2CTNC2cjDR/v27TMdchk8eHBaWlrD7pTCnJycTJfNGgwGrVbr7e0thEhLSxswYIBarRZC9OvXr6CgICMjw5Y7qpyMjIzi4uJOnToJIdLS0uLi4qQ3Oi4u7vLlyxcvXrTx/t2xp5566qmnngoODja17Nu377777pPe2dDQ0JCQkIMHD9puB5XxzTffvPXWWzt37jT+cYX0vn37Bg0aJFVdGTJkSGZm5tWrV226j3dEp9OlpKSMHj16/fr1q1evNk18a3xos7OzG8GHtsFgIFRGbm5uixYtpNv+/v5lZWU3b9607S4p5W9/+1vfvn27desmhMjLyzM9TWdnZ19f39zcXJvunQKeffbZkJCQqKioN998UxoIzd9NFxcXHx8fe3+ae/bsuXDhQlJSknmj+bsphPD397frEUIIERkZ6erqmp+f/8QTTyQkJBgMBvHnd9PDw8PLy8uu380rV64YDIa5c+devHjx5MmT0dHRZ86cEX9+NzUaTSP40DYkhGWUoVar9Xq9dFu6odFobLpHylixYkVKSorplKdara6urjb9r06nawT1lhYvXvz000/v379/zpw5Xbt27d69u0ajaUxPs6ys7Iknnvj8889r1KJsfO/mBx98IN147rnnwsLCdu3aFR8fb/67KYTQ6/V2/TSdnJyMRuPcuXOlP2t0Ot1bb721bt26xvduNiTMCJUREhJi+ms6JyenefPmppil/frnP//53nvv7d27NzAwUGoJCQnJycmRbldUVBQVFZkfarNTnp6egYGBEydOjI+P/+KLL8Sfn2ZpaWlpaaldP820tLScnJypU6fee++9S5YsOXHixL333mswGMyfphAiJyfHrp+mOR8fn8jIyAsXLog//27euHGjsrLSrp9mUFCQk5NTZGSk9GPnzp0vXbokan1ob968addPs4FhIFRGQkLC9u3bpUMxW7duTUhIsPUe3am1a9e+/fbbKSkpLVu2NDUmJCSkpKRIxdOTk5PDw8PNLyCxO3q9XqfTSbe1Wu2JEydat24thEhISPj2229LS0uFENu2bYuJiQkJCbHljt6Z3r177927d82aNWvWrJk2bVpoaOiaNWucnJzi4+OPHj0qfY3++OOPVVVVvXr1svXOWk+n05lmfleuXDl27Fjnzp2FEAkJCTt37iwvLxdCbNu2rWfPnn5+spdHuHu4uroOGzbs559/ln78+eefpUExISHhu+++kz6027dv79atm/lvLtyGrdM69mfevHmxsbHu7u6hoaGxsbGHDx82Go23bt265557+vbtm5iYGBAQcO7cOVvv5h3JyclRqVStW7eO/cPXX38t/deUKVMiIyMfeeQRPz+/L7/80rb7eYcuX74cEBAwZsyYKVOmtGnTZtCgQRUVFdJ/TZgwoUuXLg8//LCfn993331n2/1U0IcffmhKjRqNxoULF7Zt2zYpKSkwMHDNmjU23LE7d/78+aCgoHHjxiUmJvr4+MydO1dqNxgMCQkJ3bp1mz59uq+v7/fff2/T3VTAkSNH/P39p0+fPnz48A4dOly9elVqnzhxoulD++2339p2J+0LVp+QLTMz0zwIExYWJuXuqqqq9u7de+vWrUGDBvn4+NhuBxWg1WpPnjxp3tK2bVvpYnOj0bh///6cnJxevXq1adPGRjuomMuXLx87dqyysrJjx44xMTGmdqPRmJaWlpeX17t371atWtlwD5V148aNgoKC8PBwU8vhw4czMzO7det223Xm7n5nz549e/aswWCIiooKCwsztRsMhtTU1Pz8/AceeMCuJ/cmBQUFe/fubdasWZ8+fUxnYYxG4759+3JzcxvZh7YBYCAEAACHhnOEAADg0DAQAgCAQ8NACAAADg0DIQAAODQMhAAA4NAwEAIAgEPDQAgAAA4NAyEAADg0DIQAAODQMBACAIBDw0AIAAAO7f8B6pyYWgcnnEMAAAAASUVORK5CYII=", "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.7.0" }, "kernelspec": { "name": "julia-1.7", "display_name": "Julia 1.7.0", "language": "julia" } }, "nbformat": 4 }