diff --git a/simpegPF/Dev/GRAV/Intgrl_GRAV_Inv_Driver.py b/simpegPF/Dev/GRAV/Intgrl_GRAV_Inv_Driver.py index a6d9399d..27fa00bf 100644 --- a/simpegPF/Dev/GRAV/Intgrl_GRAV_Inv_Driver.py +++ b/simpegPF/Dev/GRAV/Intgrl_GRAV_Inv_Driver.py @@ -3,16 +3,14 @@ from SimPEG import * import simpegPF as PF import pylab as plt import os -home_dir = 'C:\\Users\\dominiquef.MIRAGEOSCIENCE\\ownCloud\\Research\\CraigModel\\GRAV\\checkerboard_tests' +home_dir = '.' + inpfile = 'PYGRAV3D_inv.inp' -dsep = '\\' -os.chdir(home_dir) + +dsep = os.path.sep plt.close('all') #%% User input -# Initial beta -beta_in = 1e-2 - # Treshold values for compact norm eps_p = 1e-3 # Small model values eps_q = 1e-3 # Small model gradient @@ -122,22 +120,20 @@ reg.mref = mref reg.wght = wr -# Create pre-conditioner -diagA = np.sum(prob.G**2.,axis=0) + beta_in*(reg.W.T*reg.W).diagonal() -PC = Utils.sdiag(diagA**-1.) - # Data misfit function dmis = DataMisfit.l2_DataMisfit(survey) dmis.Wd = 1./wd opt = Optimization.ProjectedGNCG(maxIter=20,lower=bounds[0],upper=bounds[1], maxIterCG= 50, tolCG = 1e-4) -opt.approxHinv = PC +invProb = InvProblem.BaseInvProblem(dmis, reg, opt) -invProb = InvProblem.BaseInvProblem(dmis, reg, opt, beta = beta_in) +# Add directives to the inversion beta = Directives.BetaSchedule(coolingFactor=2, coolingRate=1) +beta_init = Directives.BetaEstimate_ByEig() target = Directives.TargetMisfit() +update_Jacobi = Directives.Update_lin_PreCond(onlyOnStart=True) -inv = Inversion.BaseInversion(invProb, directiveList=[beta,target]) +inv = Inversion.BaseInversion(invProb, directiveList=[beta,target,beta_init,update_Jacobi]) m0 = mstart @@ -148,7 +144,6 @@ m_out = actvMap*mrec # Write result Mesh.TensorMesh.writeModelUBC(mesh,'SimPEG_inv_l2l2.sus',m_out) -#Utils.meshutils.writeUBCTensorModel(home_dir+dsep+'wr.dat',mesh,wr_out) # Plot predicted pred = prob.fields(mrec) @@ -204,25 +199,23 @@ reg.eps_p = eps_p reg.eps_q = eps_q reg.norms = lpnorms -diagA = np.sum(prob.G**2.,axis=0) + beta_in*(reg.W.T*reg.W).diagonal() -PC = Utils.sdiag(diagA**-1.) - -#reg.alpha_s = 1. - dmis = DataMisfit.l2_DataMisfit(survey) dmis.Wd = 1./wd -opt = Optimization.ProjectedGNCG(maxIter=20 , maxIterLS = 20,lower=bounds[0],upper=bounds[1], maxIterCG= 50, tolCG = 1e-4) -opt.approxHinv = PC -#opt.phim_last = reg.eval(mrec) -# opt = Optimization.InexactGaussNewton(maxIter=6) +opt = Optimization.ProjectedGNCG(maxIter=20 , maxIterLS = 20,lower=bounds[0],upper=bounds[1], maxIterCG= 50, tolCG = 1e-4) + invProb = InvProblem.BaseInvProblem(dmis, reg, opt, beta = invProb.beta) + +# Add directives to the inversion beta = Directives.BetaSchedule(coolingFactor=1, coolingRate=1) -#betaest = Directives.BetaEstimate_ByEig() +update_beta = Directives.Scale_Beta(tol = 0.05) target = Directives.TargetMisfit() IRLS =Directives.Update_IRLS( phi_m_last = phim, phi_d_last = phid ) +update_Jacobi = Directives.Update_lin_PreCond(onlyOnStart=False) +save_log = Directives.SaveOutputEveryIteration() +save_log.fileName = 'LogName_blabla' -inv = Inversion.BaseInversion(invProb, directiveList=[beta,IRLS]) +inv = Inversion.BaseInversion(invProb, directiveList=[beta,IRLS,update_beta,update_Jacobi,save_log]) m0 = mrec @@ -242,24 +235,26 @@ print "Final misfit:" + str(np.sum( ((d-pred)/wd)**2. ) ) #%% Plot out a section of the model yslice = midx + m_out[m_out==-100] = np.nan + plt.figure() ax = plt.subplot(221) -mesh.plotSlice(m_out, ax = ax, normal = 'Z', ind=-18, clim = (vmin,vmax) , pcolorOpts = {'cmap':'bwr'}) +mesh.plotSlice(m_out, ax = ax, normal = 'Z', ind=-5, clim = (mrec.min(), mrec.max()), pcolorOpts = {'cmap':'bwr'}) plt.plot(np.array([mesh.vectorCCx[0],mesh.vectorCCx[-1]]), np.array([mesh.vectorCCy[yslice],mesh.vectorCCy[yslice]]),c='w',linestyle = '--') -plt.title('Z: ' + str(mesh.vectorCCz[-18]) + ' m') +plt.title('Z: ' + str(mesh.vectorCCz[-5]) + ' m') plt.xlabel('x');plt.ylabel('z') plt.gca().set_aspect('equal', adjustable='box') ax = plt.subplot(222) -mesh.plotSlice(m_out, ax = ax, normal = 'Z', ind=-23, clim = (vmin,vmax), pcolorOpts = {'cmap':'bwr'}) +mesh.plotSlice(m_out, ax = ax, normal = 'Z', ind=-8, clim = (mrec.min(), mrec.max()), pcolorOpts = {'cmap':'bwr'}) plt.plot(np.array([mesh.vectorCCx[0],mesh.vectorCCx[-1]]), np.array([mesh.vectorCCy[yslice],mesh.vectorCCy[yslice]]),c='w',linestyle = '--') -plt.title('Z: ' + str(mesh.vectorCCz[-23]) + ' m') +plt.title('Z: ' + str(mesh.vectorCCz[-8]) + ' m') plt.xlabel('x');plt.ylabel('z') plt.gca().set_aspect('equal', adjustable='box') ax = plt.subplot(212) -mesh.plotSlice(m_out, ax = ax, normal = 'Y', ind=yslice, clim = (vmin,vmax), pcolorOpts = {'cmap':'bwr'}) +mesh.plotSlice(m_out, ax = ax, normal = 'Y', ind=yslice, clim = (mrec.min(), mrec.max()), pcolorOpts = {'cmap':'bwr'}) plt.title('Cross Section') plt.xlabel('x');plt.ylabel('z') plt.gca().set_aspect('equal', adjustable='box') diff --git a/simpegPF/Dev/GRAV/PYGRAV3D_inv.inp b/simpegPF/Dev/GRAV/PYGRAV3D_inv.inp index 828f95d7..ba5850d6 100644 --- a/simpegPF/Dev/GRAV/PYGRAV3D_inv.inp +++ b/simpegPF/Dev/GRAV/PYGRAV3D_inv.inp @@ -1,7 +1,7 @@ Mesh_10m.msh ! Mesh file Grav.dat ! Obsfile Gaussian.topo ! Topofile | null -VALUE 1e-4 ! Starting model +VALUE 0.0 ! Starting model VALUE 0.0 ! Reference model DEFAULT ! Cell based weight file 1 ! target chi factor | DEFAULT=1 diff --git a/simpegPF/Dev/GRAV/Pred0.dat b/simpegPF/Dev/GRAV/Pred0.dat index 2db7a0eb..4f7214dc 100644 --- a/simpegPF/Dev/GRAV/Pred0.dat +++ b/simpegPF/Dev/GRAV/Pred0.dat @@ -1,192 +1,617 @@ -191 -3.624724e+05 6.008600e+06 2.185513e+03 2.950277e-02 5.000000e-02 -3.630474e+05 6.007742e+06 2.182354e+03 2.958855e-02 5.000000e-02 -3.658851e+05 6.010287e+06 2.181814e+03 2.909199e-02 5.000000e-02 -3.634466e+05 6.005839e+06 2.190788e+03 2.947114e-02 5.000000e-02 -3.641340e+05 6.005157e+06 2.173695e+03 2.936905e-02 5.000000e-02 -3.644170e+05 6.003660e+06 2.189776e+03 2.896865e-02 5.000000e-02 -3.673445e+05 6.005950e+06 2.194062e+03 2.885706e-02 5.000000e-02 -3.605751e+05 6.011702e+06 2.267910e+03 2.855588e-02 5.000000e-02 -3.687754e+05 6.014567e+06 2.230211e+03 2.611820e-02 5.000000e-02 -3.600941e+05 6.011233e+06 2.205465e+03 2.849530e-02 5.000000e-02 -3.662178e+05 6.013153e+06 2.310513e+03 2.835436e-02 5.000000e-02 -3.623110e+05 6.012305e+06 2.287369e+03 2.878394e-02 5.000000e-02 -3.666170e+05 6.013924e+06 2.328716e+03 2.779231e-02 5.000000e-02 -3.633072e+05 6.013055e+06 2.350074e+03 2.838244e-02 5.000000e-02 -3.593416e+05 6.011481e+06 2.188002e+03 2.819680e-02 5.000000e-02 -3.592337e+05 6.010404e+06 2.202372e+03 2.842105e-02 5.000000e-02 -3.670648e+05 6.010916e+06 2.176537e+03 2.865192e-02 5.000000e-02 -3.594689e+05 6.013023e+06 2.187709e+03 2.760522e-02 5.000000e-02 -3.606836e+05 6.009231e+06 2.185043e+03 2.913644e-02 5.000000e-02 -3.634058e+05 6.008465e+06 2.171398e+03 2.959731e-02 5.000000e-02 -3.630400e+05 6.007708e+06 2.181436e+03 2.959071e-02 5.000000e-02 -3.634243e+05 6.005860e+06 2.188994e+03 2.947753e-02 5.000000e-02 -3.619294e+05 6.007571e+06 2.196007e+03 2.945406e-02 5.000000e-02 -3.613222e+05 6.010390e+06 2.161111e+03 2.912337e-02 5.000000e-02 -3.609065e+05 6.005685e+06 2.192528e+03 2.913821e-02 5.000000e-02 -3.600812e+05 6.004739e+06 2.201676e+03 2.870265e-02 5.000000e-02 -3.582056e+05 6.003746e+06 2.328486e+03 2.813477e-02 5.000000e-02 -3.644598e+05 6.010922e+06 2.168886e+03 2.914594e-02 5.000000e-02 -3.660609e+05 6.004171e+06 2.184777e+03 2.884816e-02 5.000000e-02 -3.662338e+05 6.004030e+06 2.186003e+03 2.876446e-02 5.000000e-02 -3.594967e+05 6.014888e+06 2.161400e+03 2.641712e-02 5.000000e-02 -3.596059e+05 6.013750e+06 2.154978e+03 2.730693e-02 5.000000e-02 -3.598438e+05 6.013214e+06 2.185253e+03 2.765375e-02 5.000000e-02 -3.599208e+05 6.012727e+06 2.203143e+03 2.790024e-02 5.000000e-02 -3.600763e+05 6.012250e+06 2.220439e+03 2.817421e-02 5.000000e-02 -3.611366e+05 6.011684e+06 2.297049e+03 2.871383e-02 5.000000e-02 -3.615742e+05 6.011880e+06 2.309140e+03 2.912090e-02 5.000000e-02 -3.620552e+05 6.012062e+06 2.303044e+03 2.912972e-02 5.000000e-02 -3.626264e+05 6.012780e+06 2.289942e+03 2.863952e-02 5.000000e-02 -3.630185e+05 6.012978e+06 2.325430e+03 2.841604e-02 5.000000e-02 -3.634810e+05 6.013343e+06 2.355786e+03 2.840198e-02 5.000000e-02 -3.639229e+05 6.013535e+06 2.357851e+03 2.835332e-02 5.000000e-02 -3.644825e+05 6.012914e+06 2.363058e+03 2.837733e-02 5.000000e-02 -3.646570e+05 6.012452e+06 2.351115e+03 2.838988e-02 5.000000e-02 -3.649237e+05 6.012810e+06 2.359124e+03 2.843872e-02 5.000000e-02 -3.651010e+05 6.013278e+06 2.358159e+03 2.837036e-02 5.000000e-02 -3.656116e+05 6.013387e+06 2.346855e+03 2.825637e-02 5.000000e-02 -3.661738e+05 6.013631e+06 2.317519e+03 2.820977e-02 5.000000e-02 -3.668565e+05 6.014278e+06 2.312226e+03 2.754386e-02 5.000000e-02 -3.674190e+05 6.014362e+06 2.298255e+03 2.724274e-02 5.000000e-02 -3.679217e+05 6.014372e+06 2.261726e+03 2.672379e-02 5.000000e-02 -3.684490e+05 6.014348e+06 2.244935e+03 2.643617e-02 5.000000e-02 -3.691523e+05 6.014957e+06 2.201091e+03 2.568247e-02 5.000000e-02 -3.696587e+05 6.014932e+06 2.201995e+03 2.542427e-02 5.000000e-02 -3.590002e+05 6.012880e+06 2.180407e+03 2.753383e-02 5.000000e-02 -3.588834e+05 6.012462e+06 2.177161e+03 2.767358e-02 5.000000e-02 -3.585091e+05 6.012213e+06 2.167542e+03 2.764672e-02 5.000000e-02 -3.589367e+05 6.011716e+06 2.184280e+03 2.797252e-02 5.000000e-02 -3.588939e+05 6.011250e+06 2.215462e+03 2.804477e-02 5.000000e-02 -3.585927e+05 6.010922e+06 2.213598e+03 2.802561e-02 5.000000e-02 -3.582362e+05 6.010625e+06 2.228342e+03 2.800097e-02 5.000000e-02 -3.577469e+05 6.010694e+06 2.268556e+03 2.831717e-02 5.000000e-02 -3.574252e+05 6.011061e+06 2.307562e+03 2.796307e-02 5.000000e-02 -3.584220e+05 6.011412e+06 2.303356e+03 2.808073e-02 5.000000e-02 -3.593649e+05 6.011057e+06 2.167662e+03 2.838074e-02 5.000000e-02 -3.591553e+05 6.009881e+06 2.166452e+03 2.857786e-02 5.000000e-02 -3.595398e+05 6.009477e+06 2.160818e+03 2.880609e-02 5.000000e-02 -3.600429e+05 6.009469e+06 2.165836e+03 2.896599e-02 5.000000e-02 -3.609785e+05 6.009612e+06 2.160491e+03 2.920621e-02 5.000000e-02 -3.612110e+05 6.009171e+06 2.161913e+03 2.931579e-02 5.000000e-02 -3.615351e+05 6.008806e+06 2.162689e+03 2.941500e-02 5.000000e-02 -3.616756e+05 6.008330e+06 2.161829e+03 2.947662e-02 5.000000e-02 -3.586542e+05 6.009829e+06 2.201714e+03 2.833175e-02 5.000000e-02 -3.583539e+05 6.009427e+06 2.297539e+03 2.866337e-02 5.000000e-02 -3.579539e+05 6.009058e+06 2.333486e+03 2.859184e-02 5.000000e-02 -3.574419e+05 6.008757e+06 2.449519e+03 2.803439e-02 5.000000e-02 -3.571929e+05 6.008326e+06 2.470826e+03 2.787830e-02 5.000000e-02 -3.568360e+05 6.008077e+06 2.500811e+03 2.762602e-02 5.000000e-02 -3.613670e+05 6.006269e+06 2.186354e+03 2.931716e-02 5.000000e-02 -3.609547e+05 6.005975e+06 2.201936e+03 2.916378e-02 5.000000e-02 -3.605453e+05 6.005187e+06 2.192131e+03 2.896361e-02 5.000000e-02 -3.596760e+05 6.004399e+06 2.216005e+03 2.846017e-02 5.000000e-02 -3.591839e+05 6.004206e+06 2.223183e+03 2.822088e-02 5.000000e-02 -3.586900e+05 6.004318e+06 2.240749e+03 2.800929e-02 5.000000e-02 -3.582144e+05 6.004122e+06 2.267646e+03 2.815678e-02 5.000000e-02 -3.574795e+05 6.004361e+06 2.273342e+03 2.785160e-02 5.000000e-02 -3.573075e+05 6.005061e+06 2.257032e+03 2.768883e-02 5.000000e-02 -3.577691e+05 6.005261e+06 2.257991e+03 2.782176e-02 5.000000e-02 -3.582753e+05 6.005377e+06 2.253481e+03 2.824865e-02 5.000000e-02 -3.589901e+05 6.005582e+06 2.225082e+03 2.843000e-02 5.000000e-02 -3.594262e+05 6.005850e+06 2.208702e+03 2.867859e-02 5.000000e-02 -3.599277e+05 6.005832e+06 2.206934e+03 2.884661e-02 5.000000e-02 -3.616705e+05 6.006416e+06 2.159826e+03 2.944395e-02 5.000000e-02 -3.689559e+05 6.013437e+06 2.160188e+03 2.686265e-02 5.000000e-02 -3.689809e+05 6.012858e+06 2.160165e+03 2.713001e-02 5.000000e-02 -3.689575e+05 6.012164e+06 2.166898e+03 2.742852e-02 5.000000e-02 -3.687815e+05 6.011610e+06 2.169930e+03 2.772659e-02 5.000000e-02 -3.685491e+05 6.011021e+06 2.161041e+03 2.805256e-02 5.000000e-02 -3.683284e+05 6.010552e+06 2.162389e+03 2.829110e-02 5.000000e-02 -3.681385e+05 6.009930e+06 2.160917e+03 2.852168e-02 5.000000e-02 -3.683213e+05 6.009170e+06 2.161018e+03 2.855719e-02 5.000000e-02 -3.681582e+05 6.008637e+06 2.160607e+03 2.868224e-02 5.000000e-02 -3.678554e+05 6.008233e+06 2.161901e+03 2.883354e-02 5.000000e-02 -3.676810e+05 6.007757e+06 2.161971e+03 2.891626e-02 5.000000e-02 -3.675351e+05 6.007272e+06 2.161757e+03 2.896734e-02 5.000000e-02 -3.674284e+05 6.006771e+06 2.161644e+03 2.898216e-02 5.000000e-02 -3.672650e+05 6.006275e+06 2.161193e+03 2.899572e-02 5.000000e-02 -3.670451e+05 6.005801e+06 2.161079e+03 2.900593e-02 5.000000e-02 -3.668871e+05 6.005298e+06 2.168517e+03 2.895119e-02 5.000000e-02 -3.669246e+05 6.004782e+06 2.160798e+03 2.884703e-02 5.000000e-02 -3.674886e+05 6.004466e+06 2.160969e+03 2.858275e-02 5.000000e-02 -3.680010e+05 6.007243e+06 2.186598e+03 2.873252e-02 5.000000e-02 -3.685053e+05 6.007066e+06 2.247928e+03 2.844525e-02 5.000000e-02 -3.689689e+05 6.006886e+06 2.294719e+03 2.876087e-02 5.000000e-02 -3.694464e+05 6.006798e+06 2.340813e+03 2.847357e-02 5.000000e-02 -3.699384e+05 6.006629e+06 2.452104e+03 2.790729e-02 5.000000e-02 -3.702331e+05 6.006504e+06 2.478823e+03 2.768065e-02 5.000000e-02 -3.648724e+05 6.003990e+06 2.162273e+03 2.903588e-02 5.000000e-02 -3.615819e+05 6.006914e+06 2.166497e+03 2.945036e-02 5.000000e-02 -3.619111e+05 6.005892e+06 2.161213e+03 2.943030e-02 5.000000e-02 -3.624320e+05 6.006129e+06 2.161200e+03 2.951822e-02 5.000000e-02 -3.628240e+05 6.006286e+06 2.162131e+03 2.955949e-02 5.000000e-02 -3.632946e+05 6.006160e+06 2.161080e+03 2.956698e-02 5.000000e-02 -3.676218e+05 6.013497e+06 2.169281e+03 2.740393e-02 5.000000e-02 -3.679297e+05 6.013094e+06 2.168942e+03 2.751932e-02 5.000000e-02 -3.677949e+05 6.012264e+06 2.167650e+03 2.796881e-02 5.000000e-02 -3.672235e+05 6.011537e+06 2.168687e+03 2.843018e-02 5.000000e-02 -3.666300e+05 6.012199e+06 2.168907e+03 2.836406e-02 5.000000e-02 -3.659040e+05 6.011713e+06 2.165875e+03 2.871230e-02 5.000000e-02 -3.659362e+05 6.011209e+06 2.167131e+03 2.887116e-02 5.000000e-02 -3.657312e+05 6.010915e+06 2.166870e+03 2.899740e-02 5.000000e-02 -3.655597e+05 6.011327e+06 2.167944e+03 2.891109e-02 5.000000e-02 -3.650121e+05 6.011050e+06 2.167057e+03 2.906459e-02 5.000000e-02 -3.644476e+05 6.011501e+06 2.165267e+03 2.897208e-02 5.000000e-02 -3.642410e+05 6.011986e+06 2.164931e+03 2.880292e-02 5.000000e-02 -3.640583e+05 6.012445e+06 2.165348e+03 2.861521e-02 5.000000e-02 -3.637110e+05 6.012686e+06 2.165028e+03 2.851203e-02 5.000000e-02 -3.632662e+05 6.012413e+06 2.165095e+03 2.863726e-02 5.000000e-02 -3.628353e+05 6.012198e+06 2.165437e+03 2.872215e-02 5.000000e-02 -3.626366e+05 6.011643e+06 2.167954e+03 2.892334e-02 5.000000e-02 -3.619752e+05 6.011006e+06 2.169750e+03 2.905700e-02 5.000000e-02 -3.608859e+05 6.012099e+06 2.332313e+03 2.884199e-02 5.000000e-02 -3.607787e+05 6.012623e+06 2.334431e+03 2.864554e-02 5.000000e-02 -3.604558e+05 6.012994e+06 2.340932e+03 2.825683e-02 5.000000e-02 -3.604541e+05 6.013514e+06 2.417598e+03 2.775451e-02 5.000000e-02 -3.605505e+05 6.014052e+06 2.428877e+03 2.743902e-02 5.000000e-02 -3.605665e+05 6.014532e+06 2.431259e+03 2.703843e-02 5.000000e-02 -3.609731e+05 6.014173e+06 2.431723e+03 2.750775e-02 5.000000e-02 -3.614862e+05 6.014207e+06 2.463004e+03 2.751270e-02 5.000000e-02 -3.619207e+05 6.014382e+06 2.454613e+03 2.747475e-02 5.000000e-02 -3.624145e+05 6.014385e+06 2.450971e+03 2.754979e-02 5.000000e-02 -3.640306e+05 6.003363e+06 2.230716e+03 2.912848e-02 5.000000e-02 -3.636639e+05 6.003018e+06 2.254717e+03 2.944835e-02 5.000000e-02 -3.632341e+05 6.002875e+06 2.270788e+03 2.935057e-02 5.000000e-02 -3.630920e+05 6.002515e+06 2.342944e+03 2.897828e-02 5.000000e-02 -3.627738e+05 6.002067e+06 2.449163e+03 2.845754e-02 5.000000e-02 -3.625173e+05 6.007599e+06 2.169206e+03 2.958475e-02 5.000000e-02 -3.635370e+05 6.008139e+06 2.164473e+03 2.963331e-02 5.000000e-02 -3.627917e+05 6.008156e+06 2.169742e+03 2.959324e-02 5.000000e-02 -3.633881e+05 6.009013e+06 2.164473e+03 2.956369e-02 5.000000e-02 -3.663215e+05 6.003654e+06 2.202612e+03 2.858425e-02 5.000000e-02 -3.664925e+05 6.003153e+06 2.185419e+03 2.839115e-02 5.000000e-02 -3.664874e+05 6.002625e+06 2.175218e+03 2.818385e-02 5.000000e-02 -3.664629e+05 6.002138e+06 2.182963e+03 2.793174e-02 5.000000e-02 -3.665085e+05 6.001658e+06 2.183867e+03 2.764954e-02 5.000000e-02 -3.667288e+05 6.001228e+06 2.186336e+03 2.730574e-02 5.000000e-02 -3.669117e+05 6.000761e+06 2.194478e+03 2.690493e-02 5.000000e-02 -3.672005e+05 6.000532e+06 2.193351e+03 2.667164e-02 5.000000e-02 -3.615060e+05 6.010861e+06 2.164254e+03 2.903521e-02 5.000000e-02 -3.608114e+05 6.010798e+06 2.166447e+03 2.890123e-02 5.000000e-02 -3.678002e+05 6.003356e+06 2.180167e+03 2.808437e-02 5.000000e-02 -3.682322e+05 6.003116e+06 2.219280e+03 2.777212e-02 5.000000e-02 -3.685762e+05 6.002753e+06 2.243722e+03 2.751404e-02 5.000000e-02 -3.690196e+05 6.002528e+06 2.262227e+03 2.771937e-02 5.000000e-02 -3.604802e+05 6.004043e+06 2.252471e+03 2.865938e-02 5.000000e-02 -3.609087e+05 6.003672e+06 2.320061e+03 2.916318e-02 5.000000e-02 -3.609935e+05 6.003061e+06 2.341764e+03 2.896186e-02 5.000000e-02 -3.610941e+05 6.002503e+06 2.353045e+03 2.870830e-02 5.000000e-02 -3.613520e+05 6.002090e+06 2.402026e+03 2.841739e-02 5.000000e-02 -3.616648e+05 6.001616e+06 2.460820e+03 2.804598e-02 5.000000e-02 -3.621238e+05 6.001462e+06 2.490509e+03 2.797261e-02 5.000000e-02 -3.657493e+05 6.004426e+06 2.169680e+03 2.901485e-02 5.000000e-02 -3.636824e+05 6.000022e+06 2.865274e+03 2.708391e-02 5.000000e-02 -3.571072e+05 6.007140e+06 2.504079e+03 2.777789e-02 5.000000e-02 -3.687904e+05 6.004471e+06 2.394426e+03 2.822128e-02 5.000000e-02 -3.660711e+05 6.001991e+06 2.247601e+03 2.794515e-02 5.000000e-02 -3.565341e+05 6.013529e+06 2.495766e+03 2.596333e-02 5.000000e-02 -3.572975e+05 6.002759e+06 2.586419e+03 2.739037e-02 5.000000e-02 -3.602847e+05 6.000220e+06 2.741396e+03 2.681741e-02 5.000000e-02 -3.592838e+05 6.002004e+06 2.769861e+03 2.752123e-02 5.000000e-02 -3.624843e+05 5.999891e+06 2.693787e+03 2.696451e-02 5.000000e-02 +616 +1.050000e+02 1.550000e+02 2.408312e+02 0.000000e+00 2.000000e-03 +1.150000e+02 1.550000e+02 2.440469e+02 0.000000e+00 2.000000e-03 +1.250000e+02 1.550000e+02 2.471986e+02 0.000000e+00 2.000000e-03 +1.350000e+02 1.550000e+02 2.502569e+02 0.000000e+00 2.000000e-03 +1.450000e+02 1.550000e+02 2.531923e+02 0.000000e+00 2.000000e-03 +1.550000e+02 1.550000e+02 2.559753e+02 0.000000e+00 2.000000e-03 +1.650000e+02 1.550000e+02 2.585770e+02 0.000000e+00 2.000000e-03 +1.750000e+02 1.550000e+02 2.609695e+02 0.000000e+00 2.000000e-03 +1.850000e+02 1.550000e+02 2.631267e+02 0.000000e+00 2.000000e-03 +1.950000e+02 1.550000e+02 2.650244e+02 0.000000e+00 2.000000e-03 +2.050000e+02 1.550000e+02 2.666411e+02 0.000000e+00 2.000000e-03 +2.150000e+02 1.550000e+02 2.679578e+02 0.000000e+00 2.000000e-03 +2.250000e+02 1.550000e+02 2.689593e+02 0.000000e+00 2.000000e-03 +2.350000e+02 1.550000e+02 2.696337e+02 0.000000e+00 2.000000e-03 +2.450000e+02 1.550000e+02 2.699729e+02 0.000000e+00 2.000000e-03 +2.550000e+02 1.550000e+02 2.699729e+02 0.000000e+00 2.000000e-03 +2.650000e+02 1.550000e+02 2.696337e+02 0.000000e+00 2.000000e-03 +2.750000e+02 1.550000e+02 2.689593e+02 0.000000e+00 2.000000e-03 +2.850000e+02 1.550000e+02 2.679578e+02 0.000000e+00 2.000000e-03 +2.950000e+02 1.550000e+02 2.666411e+02 0.000000e+00 2.000000e-03 +3.050000e+02 1.550000e+02 2.650244e+02 0.000000e+00 2.000000e-03 +3.150000e+02 1.550000e+02 2.631267e+02 0.000000e+00 2.000000e-03 +3.250000e+02 1.550000e+02 2.609695e+02 0.000000e+00 2.000000e-03 +3.350000e+02 1.550000e+02 2.585770e+02 0.000000e+00 2.000000e-03 +3.450000e+02 1.550000e+02 2.559753e+02 0.000000e+00 2.000000e-03 +3.550000e+02 1.550000e+02 2.531923e+02 0.000000e+00 2.000000e-03 +3.650000e+02 1.550000e+02 2.502569e+02 0.000000e+00 2.000000e-03 +3.750000e+02 1.550000e+02 2.471986e+02 0.000000e+00 2.000000e-03 +1.050000e+02 1.650000e+02 2.427601e+02 0.000000e+00 2.000000e-03 +1.150000e+02 1.650000e+02 2.460869e+02 0.000000e+00 2.000000e-03 +1.250000e+02 1.650000e+02 2.493475e+02 0.000000e+00 2.000000e-03 +1.350000e+02 1.650000e+02 2.525115e+02 0.000000e+00 2.000000e-03 +1.450000e+02 1.650000e+02 2.555483e+02 0.000000e+00 2.000000e-03 +1.550000e+02 1.650000e+02 2.584275e+02 0.000000e+00 2.000000e-03 +1.650000e+02 1.650000e+02 2.611190e+02 0.000000e+00 2.000000e-03 +1.750000e+02 1.650000e+02 2.635942e+02 0.000000e+00 2.000000e-03 +1.850000e+02 1.650000e+02 2.658260e+02 0.000000e+00 2.000000e-03 +1.950000e+02 1.650000e+02 2.677892e+02 0.000000e+00 2.000000e-03 +2.050000e+02 1.650000e+02 2.694617e+02 0.000000e+00 2.000000e-03 +2.150000e+02 1.650000e+02 2.708240e+02 0.000000e+00 2.000000e-03 +2.250000e+02 1.650000e+02 2.718600e+02 0.000000e+00 2.000000e-03 +2.350000e+02 1.650000e+02 2.725577e+02 0.000000e+00 2.000000e-03 +2.450000e+02 1.650000e+02 2.729086e+02 0.000000e+00 2.000000e-03 +2.550000e+02 1.650000e+02 2.729086e+02 0.000000e+00 2.000000e-03 +2.650000e+02 1.650000e+02 2.725577e+02 0.000000e+00 2.000000e-03 +2.750000e+02 1.650000e+02 2.718600e+02 0.000000e+00 2.000000e-03 +2.850000e+02 1.650000e+02 2.708240e+02 0.000000e+00 2.000000e-03 +2.950000e+02 1.650000e+02 2.694617e+02 0.000000e+00 2.000000e-03 +3.050000e+02 1.650000e+02 2.677892e+02 0.000000e+00 2.000000e-03 +3.150000e+02 1.650000e+02 2.658260e+02 0.000000e+00 2.000000e-03 +3.250000e+02 1.650000e+02 2.635942e+02 0.000000e+00 2.000000e-03 +3.350000e+02 1.650000e+02 2.611190e+02 0.000000e+00 2.000000e-03 +3.450000e+02 1.650000e+02 2.584275e+02 0.000000e+00 2.000000e-03 +3.550000e+02 1.650000e+02 2.555483e+02 0.000000e+00 2.000000e-03 +3.650000e+02 1.650000e+02 2.525115e+02 0.000000e+00 2.000000e-03 +3.750000e+02 1.650000e+02 2.493475e+02 0.000000e+00 2.000000e-03 +1.050000e+02 1.750000e+02 2.445174e+02 0.000000e+00 2.000000e-03 +1.150000e+02 1.750000e+02 2.479454e+02 0.000000e+00 2.000000e-03 +1.250000e+02 1.750000e+02 2.513052e+02 0.000000e+00 2.000000e-03 +1.350000e+02 1.750000e+02 2.545655e+02 0.000000e+00 2.000000e-03 +1.450000e+02 1.750000e+02 2.576947e+02 0.000000e+00 2.000000e-03 +1.550000e+02 1.750000e+02 2.606614e+02 0.000000e+00 2.000000e-03 +1.650000e+02 1.750000e+02 2.634348e+02 0.000000e+00 2.000000e-03 +1.750000e+02 1.750000e+02 2.659854e+02 0.000000e+00 2.000000e-03 +1.850000e+02 1.750000e+02 2.682850e+02 0.000000e+00 2.000000e-03 +1.950000e+02 1.750000e+02 2.703080e+02 0.000000e+00 2.000000e-03 +2.050000e+02 1.750000e+02 2.720313e+02 0.000000e+00 2.000000e-03 +2.150000e+02 1.750000e+02 2.734350e+02 0.000000e+00 2.000000e-03 +2.250000e+02 1.750000e+02 2.745027e+02 0.000000e+00 2.000000e-03 +2.350000e+02 1.750000e+02 2.752216e+02 0.000000e+00 2.000000e-03 +2.450000e+02 1.750000e+02 2.755831e+02 0.000000e+00 2.000000e-03 +2.550000e+02 1.750000e+02 2.755831e+02 0.000000e+00 2.000000e-03 +2.650000e+02 1.750000e+02 2.752216e+02 0.000000e+00 2.000000e-03 +2.750000e+02 1.750000e+02 2.745027e+02 0.000000e+00 2.000000e-03 +2.850000e+02 1.750000e+02 2.734350e+02 0.000000e+00 2.000000e-03 +2.950000e+02 1.750000e+02 2.720313e+02 0.000000e+00 2.000000e-03 +3.050000e+02 1.750000e+02 2.703080e+02 0.000000e+00 2.000000e-03 +3.150000e+02 1.750000e+02 2.682850e+02 0.000000e+00 2.000000e-03 +3.250000e+02 1.750000e+02 2.659854e+02 0.000000e+00 2.000000e-03 +3.350000e+02 1.750000e+02 2.634348e+02 0.000000e+00 2.000000e-03 +3.450000e+02 1.750000e+02 2.606614e+02 0.000000e+00 2.000000e-03 +3.550000e+02 1.750000e+02 2.576947e+02 0.000000e+00 2.000000e-03 +3.650000e+02 1.750000e+02 2.545655e+02 0.000000e+00 2.000000e-03 +3.750000e+02 1.750000e+02 2.513052e+02 0.000000e+00 2.000000e-03 +1.050000e+02 1.850000e+02 2.460835e+02 0.000000e+00 2.000000e-03 +1.150000e+02 1.850000e+02 2.496018e+02 0.000000e+00 2.000000e-03 +1.250000e+02 1.850000e+02 2.530500e+02 0.000000e+00 2.000000e-03 +1.350000e+02 1.850000e+02 2.563960e+02 0.000000e+00 2.000000e-03 +1.450000e+02 1.850000e+02 2.596076e+02 0.000000e+00 2.000000e-03 +1.550000e+02 1.850000e+02 2.626524e+02 0.000000e+00 2.000000e-03 +1.650000e+02 1.850000e+02 2.654988e+02 0.000000e+00 2.000000e-03 +1.750000e+02 1.850000e+02 2.681164e+02 0.000000e+00 2.000000e-03 +1.850000e+02 1.850000e+02 2.704766e+02 0.000000e+00 2.000000e-03 +1.950000e+02 1.850000e+02 2.725528e+02 0.000000e+00 2.000000e-03 +2.050000e+02 1.850000e+02 2.743215e+02 0.000000e+00 2.000000e-03 +2.150000e+02 1.850000e+02 2.757621e+02 0.000000e+00 2.000000e-03 +2.250000e+02 1.850000e+02 2.768578e+02 0.000000e+00 2.000000e-03 +2.350000e+02 1.850000e+02 2.775956e+02 0.000000e+00 2.000000e-03 +2.450000e+02 1.850000e+02 2.779668e+02 0.000000e+00 2.000000e-03 +2.550000e+02 1.850000e+02 2.779668e+02 0.000000e+00 2.000000e-03 +2.650000e+02 1.850000e+02 2.775956e+02 0.000000e+00 2.000000e-03 +2.750000e+02 1.850000e+02 2.768578e+02 0.000000e+00 2.000000e-03 +2.850000e+02 1.850000e+02 2.757621e+02 0.000000e+00 2.000000e-03 +2.950000e+02 1.850000e+02 2.743215e+02 0.000000e+00 2.000000e-03 +3.050000e+02 1.850000e+02 2.725528e+02 0.000000e+00 2.000000e-03 +3.150000e+02 1.850000e+02 2.704766e+02 0.000000e+00 2.000000e-03 +3.250000e+02 1.850000e+02 2.681164e+02 0.000000e+00 2.000000e-03 +3.350000e+02 1.850000e+02 2.654988e+02 0.000000e+00 2.000000e-03 +3.450000e+02 1.850000e+02 2.626524e+02 0.000000e+00 2.000000e-03 +3.550000e+02 1.850000e+02 2.596076e+02 0.000000e+00 2.000000e-03 +3.650000e+02 1.850000e+02 2.563960e+02 0.000000e+00 2.000000e-03 +3.750000e+02 1.850000e+02 2.530500e+02 0.000000e+00 2.000000e-03 +1.050000e+02 1.950000e+02 2.474409e+02 0.000000e+00 2.000000e-03 +1.150000e+02 1.950000e+02 2.510374e+02 0.000000e+00 2.000000e-03 +1.250000e+02 1.950000e+02 2.545621e+02 0.000000e+00 2.000000e-03 +1.350000e+02 1.950000e+02 2.579825e+02 0.000000e+00 2.000000e-03 +1.450000e+02 1.950000e+02 2.612655e+02 0.000000e+00 2.000000e-03 +1.550000e+02 1.950000e+02 2.643779e+02 0.000000e+00 2.000000e-03 +1.650000e+02 1.950000e+02 2.672876e+02 0.000000e+00 2.000000e-03 +1.750000e+02 1.950000e+02 2.699634e+02 0.000000e+00 2.000000e-03 +1.850000e+02 1.950000e+02 2.723760e+02 0.000000e+00 2.000000e-03 +1.950000e+02 1.950000e+02 2.744984e+02 0.000000e+00 2.000000e-03 +2.050000e+02 1.950000e+02 2.763063e+02 0.000000e+00 2.000000e-03 +2.150000e+02 1.950000e+02 2.777790e+02 0.000000e+00 2.000000e-03 +2.250000e+02 1.950000e+02 2.788990e+02 0.000000e+00 2.000000e-03 +2.350000e+02 1.950000e+02 2.796533e+02 0.000000e+00 2.000000e-03 +2.450000e+02 1.950000e+02 2.800326e+02 0.000000e+00 2.000000e-03 +2.550000e+02 1.950000e+02 2.800326e+02 0.000000e+00 2.000000e-03 +2.650000e+02 1.950000e+02 2.796533e+02 0.000000e+00 2.000000e-03 +2.750000e+02 1.950000e+02 2.788990e+02 0.000000e+00 2.000000e-03 +2.850000e+02 1.950000e+02 2.777790e+02 0.000000e+00 2.000000e-03 +2.950000e+02 1.950000e+02 2.763063e+02 0.000000e+00 2.000000e-03 +3.050000e+02 1.950000e+02 2.744984e+02 0.000000e+00 2.000000e-03 +3.150000e+02 1.950000e+02 2.723760e+02 0.000000e+00 2.000000e-03 +3.250000e+02 1.950000e+02 2.699634e+02 0.000000e+00 2.000000e-03 +3.350000e+02 1.950000e+02 2.672876e+02 0.000000e+00 2.000000e-03 +3.450000e+02 1.950000e+02 2.643779e+02 0.000000e+00 2.000000e-03 +3.550000e+02 1.950000e+02 2.612655e+02 0.000000e+00 2.000000e-03 +3.650000e+02 1.950000e+02 2.579825e+02 0.000000e+00 2.000000e-03 +3.750000e+02 1.950000e+02 2.545621e+02 0.000000e+00 2.000000e-03 +1.050000e+02 2.050000e+02 2.485739e+02 0.000000e+00 2.000000e-03 +1.150000e+02 2.050000e+02 2.522356e+02 0.000000e+00 2.000000e-03 +1.250000e+02 2.050000e+02 2.558243e+02 0.000000e+00 2.000000e-03 +1.350000e+02 2.050000e+02 2.593068e+02 0.000000e+00 2.000000e-03 +1.450000e+02 2.050000e+02 2.626493e+02 0.000000e+00 2.000000e-03 +1.550000e+02 2.050000e+02 2.658182e+02 0.000000e+00 2.000000e-03 +1.650000e+02 2.050000e+02 2.687807e+02 0.000000e+00 2.000000e-03 +1.750000e+02 2.050000e+02 2.715050e+02 0.000000e+00 2.000000e-03 +1.850000e+02 2.050000e+02 2.739614e+02 0.000000e+00 2.000000e-03 +1.950000e+02 2.050000e+02 2.761223e+02 0.000000e+00 2.000000e-03 +2.050000e+02 2.050000e+02 2.779631e+02 0.000000e+00 2.000000e-03 +2.150000e+02 2.050000e+02 2.794624e+02 0.000000e+00 2.000000e-03 +2.250000e+02 2.050000e+02 2.806028e+02 0.000000e+00 2.000000e-03 +2.350000e+02 2.050000e+02 2.813707e+02 0.000000e+00 2.000000e-03 +2.450000e+02 2.050000e+02 2.817570e+02 0.000000e+00 2.000000e-03 +2.550000e+02 2.050000e+02 2.817570e+02 0.000000e+00 2.000000e-03 +2.650000e+02 2.050000e+02 2.813707e+02 0.000000e+00 2.000000e-03 +2.750000e+02 2.050000e+02 2.806028e+02 0.000000e+00 2.000000e-03 +2.850000e+02 2.050000e+02 2.794624e+02 0.000000e+00 2.000000e-03 +2.950000e+02 2.050000e+02 2.779631e+02 0.000000e+00 2.000000e-03 +3.050000e+02 2.050000e+02 2.761223e+02 0.000000e+00 2.000000e-03 +3.150000e+02 2.050000e+02 2.739614e+02 0.000000e+00 2.000000e-03 +3.250000e+02 2.050000e+02 2.715050e+02 0.000000e+00 2.000000e-03 +3.350000e+02 2.050000e+02 2.687807e+02 0.000000e+00 2.000000e-03 +3.450000e+02 2.050000e+02 2.658182e+02 0.000000e+00 2.000000e-03 +3.550000e+02 2.050000e+02 2.626493e+02 0.000000e+00 2.000000e-03 +3.650000e+02 2.050000e+02 2.593068e+02 0.000000e+00 2.000000e-03 +3.750000e+02 2.050000e+02 2.558243e+02 0.000000e+00 2.000000e-03 +1.050000e+02 2.150000e+02 2.494692e+02 0.000000e+00 2.000000e-03 +1.150000e+02 2.150000e+02 2.531825e+02 0.000000e+00 2.000000e-03 +1.250000e+02 2.150000e+02 2.568218e+02 0.000000e+00 2.000000e-03 +1.350000e+02 2.150000e+02 2.603533e+02 0.000000e+00 2.000000e-03 +1.450000e+02 2.150000e+02 2.637429e+02 0.000000e+00 2.000000e-03 +1.550000e+02 2.150000e+02 2.669565e+02 0.000000e+00 2.000000e-03 +1.650000e+02 2.150000e+02 2.699606e+02 0.000000e+00 2.000000e-03 +1.750000e+02 2.150000e+02 2.727234e+02 0.000000e+00 2.000000e-03 +1.850000e+02 2.150000e+02 2.752143e+02 0.000000e+00 2.000000e-03 +1.950000e+02 2.150000e+02 2.774056e+02 0.000000e+00 2.000000e-03 +2.050000e+02 2.150000e+02 2.792724e+02 0.000000e+00 2.000000e-03 +2.150000e+02 2.150000e+02 2.807928e+02 0.000000e+00 2.000000e-03 +2.250000e+02 2.150000e+02 2.819493e+02 0.000000e+00 2.000000e-03 +2.350000e+02 2.150000e+02 2.827280e+02 0.000000e+00 2.000000e-03 +2.450000e+02 2.150000e+02 2.831197e+02 0.000000e+00 2.000000e-03 +2.550000e+02 2.150000e+02 2.831197e+02 0.000000e+00 2.000000e-03 +2.650000e+02 2.150000e+02 2.827280e+02 0.000000e+00 2.000000e-03 +2.750000e+02 2.150000e+02 2.819493e+02 0.000000e+00 2.000000e-03 +2.850000e+02 2.150000e+02 2.807928e+02 0.000000e+00 2.000000e-03 +2.950000e+02 2.150000e+02 2.792724e+02 0.000000e+00 2.000000e-03 +3.050000e+02 2.150000e+02 2.774056e+02 0.000000e+00 2.000000e-03 +3.150000e+02 2.150000e+02 2.752143e+02 0.000000e+00 2.000000e-03 +3.250000e+02 2.150000e+02 2.727234e+02 0.000000e+00 2.000000e-03 +3.350000e+02 2.150000e+02 2.699606e+02 0.000000e+00 2.000000e-03 +3.450000e+02 2.150000e+02 2.669565e+02 0.000000e+00 2.000000e-03 +3.550000e+02 2.150000e+02 2.637429e+02 0.000000e+00 2.000000e-03 +3.650000e+02 2.150000e+02 2.603533e+02 0.000000e+00 2.000000e-03 +3.750000e+02 2.150000e+02 2.568218e+02 0.000000e+00 2.000000e-03 +1.050000e+02 2.250000e+02 2.501165e+02 0.000000e+00 2.000000e-03 +1.150000e+02 2.250000e+02 2.538671e+02 0.000000e+00 2.000000e-03 +1.250000e+02 2.250000e+02 2.575429e+02 0.000000e+00 2.000000e-03 +1.350000e+02 2.250000e+02 2.611099e+02 0.000000e+00 2.000000e-03 +1.450000e+02 2.250000e+02 2.645335e+02 0.000000e+00 2.000000e-03 +1.550000e+02 2.250000e+02 2.677793e+02 0.000000e+00 2.000000e-03 +1.650000e+02 2.250000e+02 2.708137e+02 0.000000e+00 2.000000e-03 +1.750000e+02 2.250000e+02 2.736041e+02 0.000000e+00 2.000000e-03 +1.850000e+02 2.250000e+02 2.761201e+02 0.000000e+00 2.000000e-03 +1.950000e+02 2.250000e+02 2.783334e+02 0.000000e+00 2.000000e-03 +2.050000e+02 2.250000e+02 2.802189e+02 0.000000e+00 2.000000e-03 +2.150000e+02 2.250000e+02 2.817546e+02 0.000000e+00 2.000000e-03 +2.250000e+02 2.250000e+02 2.829227e+02 0.000000e+00 2.000000e-03 +2.350000e+02 2.250000e+02 2.837092e+02 0.000000e+00 2.000000e-03 +2.450000e+02 2.250000e+02 2.841048e+02 0.000000e+00 2.000000e-03 +2.550000e+02 2.250000e+02 2.841048e+02 0.000000e+00 2.000000e-03 +2.650000e+02 2.250000e+02 2.837092e+02 0.000000e+00 2.000000e-03 +2.750000e+02 2.250000e+02 2.829227e+02 0.000000e+00 2.000000e-03 +2.850000e+02 2.250000e+02 2.817546e+02 0.000000e+00 2.000000e-03 +2.950000e+02 2.250000e+02 2.802189e+02 0.000000e+00 2.000000e-03 +3.050000e+02 2.250000e+02 2.783334e+02 0.000000e+00 2.000000e-03 +3.150000e+02 2.250000e+02 2.761201e+02 0.000000e+00 2.000000e-03 +3.250000e+02 2.250000e+02 2.736041e+02 0.000000e+00 2.000000e-03 +3.350000e+02 2.250000e+02 2.708137e+02 0.000000e+00 2.000000e-03 +3.450000e+02 2.250000e+02 2.677793e+02 0.000000e+00 2.000000e-03 +3.550000e+02 2.250000e+02 2.645335e+02 0.000000e+00 2.000000e-03 +3.650000e+02 2.250000e+02 2.611099e+02 0.000000e+00 2.000000e-03 +3.750000e+02 2.250000e+02 2.575429e+02 0.000000e+00 2.000000e-03 +1.050000e+02 2.350000e+02 2.505080e+02 0.000000e+00 2.000000e-03 +1.150000e+02 2.350000e+02 2.542811e+02 0.000000e+00 2.000000e-03 +1.250000e+02 2.350000e+02 2.579790e+02 0.000000e+00 2.000000e-03 +1.350000e+02 2.350000e+02 2.615674e+02 0.000000e+00 2.000000e-03 +1.450000e+02 2.350000e+02 2.650117e+02 0.000000e+00 2.000000e-03 +1.550000e+02 2.350000e+02 2.682770e+02 0.000000e+00 2.000000e-03 +1.650000e+02 2.350000e+02 2.713296e+02 0.000000e+00 2.000000e-03 +1.750000e+02 2.350000e+02 2.741368e+02 0.000000e+00 2.000000e-03 +1.850000e+02 2.350000e+02 2.766679e+02 0.000000e+00 2.000000e-03 +1.950000e+02 2.350000e+02 2.788945e+02 0.000000e+00 2.000000e-03 +2.050000e+02 2.350000e+02 2.807913e+02 0.000000e+00 2.000000e-03 +2.150000e+02 2.350000e+02 2.823363e+02 0.000000e+00 2.000000e-03 +2.250000e+02 2.350000e+02 2.835114e+02 0.000000e+00 2.000000e-03 +2.350000e+02 2.350000e+02 2.843026e+02 0.000000e+00 2.000000e-03 +2.450000e+02 2.350000e+02 2.847007e+02 0.000000e+00 2.000000e-03 +2.550000e+02 2.350000e+02 2.847007e+02 0.000000e+00 2.000000e-03 +2.650000e+02 2.350000e+02 2.843026e+02 0.000000e+00 2.000000e-03 +2.750000e+02 2.350000e+02 2.835114e+02 0.000000e+00 2.000000e-03 +2.850000e+02 2.350000e+02 2.823363e+02 0.000000e+00 2.000000e-03 +2.950000e+02 2.350000e+02 2.807913e+02 0.000000e+00 2.000000e-03 +3.050000e+02 2.350000e+02 2.788945e+02 0.000000e+00 2.000000e-03 +3.150000e+02 2.350000e+02 2.766679e+02 0.000000e+00 2.000000e-03 +3.250000e+02 2.350000e+02 2.741368e+02 0.000000e+00 2.000000e-03 +3.350000e+02 2.350000e+02 2.713296e+02 0.000000e+00 2.000000e-03 +3.450000e+02 2.350000e+02 2.682770e+02 0.000000e+00 2.000000e-03 +3.550000e+02 2.350000e+02 2.650117e+02 0.000000e+00 2.000000e-03 +3.650000e+02 2.350000e+02 2.615674e+02 0.000000e+00 2.000000e-03 +3.750000e+02 2.350000e+02 2.579791e+02 0.000000e+00 2.000000e-03 +1.050000e+02 2.450000e+02 2.506390e+02 0.000000e+00 2.000000e-03 +1.150000e+02 2.450000e+02 2.544197e+02 0.000000e+00 2.000000e-03 +1.250000e+02 2.450000e+02 2.581250e+02 0.000000e+00 2.000000e-03 +1.350000e+02 2.450000e+02 2.617206e+02 0.000000e+00 2.000000e-03 +1.450000e+02 2.450000e+02 2.651717e+02 0.000000e+00 2.000000e-03 +1.550000e+02 2.450000e+02 2.684435e+02 0.000000e+00 2.000000e-03 +1.650000e+02 2.450000e+02 2.715022e+02 0.000000e+00 2.000000e-03 +1.750000e+02 2.450000e+02 2.743151e+02 0.000000e+00 2.000000e-03 +1.850000e+02 2.450000e+02 2.768512e+02 0.000000e+00 2.000000e-03 +1.950000e+02 2.450000e+02 2.790823e+02 0.000000e+00 2.000000e-03 +2.050000e+02 2.450000e+02 2.809829e+02 0.000000e+00 2.000000e-03 +2.150000e+02 2.450000e+02 2.825310e+02 0.000000e+00 2.000000e-03 +2.250000e+02 2.450000e+02 2.837084e+02 0.000000e+00 2.000000e-03 +2.350000e+02 2.450000e+02 2.845013e+02 0.000000e+00 2.000000e-03 +2.450000e+02 2.450000e+02 2.849001e+02 0.000000e+00 2.000000e-03 +2.550000e+02 2.450000e+02 2.849001e+02 0.000000e+00 2.000000e-03 +2.650000e+02 2.450000e+02 2.845013e+02 0.000000e+00 2.000000e-03 +2.750000e+02 2.450000e+02 2.837084e+02 0.000000e+00 2.000000e-03 +2.850000e+02 2.450000e+02 2.825310e+02 0.000000e+00 2.000000e-03 +2.950000e+02 2.450000e+02 2.809829e+02 0.000000e+00 2.000000e-03 +3.050000e+02 2.450000e+02 2.790823e+02 0.000000e+00 2.000000e-03 +3.150000e+02 2.450000e+02 2.768512e+02 0.000000e+00 2.000000e-03 +3.250000e+02 2.450000e+02 2.743151e+02 0.000000e+00 2.000000e-03 +3.350000e+02 2.450000e+02 2.715022e+02 0.000000e+00 2.000000e-03 +3.450000e+02 2.450000e+02 2.684435e+02 0.000000e+00 2.000000e-03 +3.550000e+02 2.450000e+02 2.651717e+02 0.000000e+00 2.000000e-03 +3.650000e+02 2.450000e+02 2.617206e+02 0.000000e+00 2.000000e-03 +3.750000e+02 2.450000e+02 2.581250e+02 0.000000e+00 2.000000e-03 +1.050000e+02 2.550000e+02 2.505080e+02 0.000000e+00 2.000000e-03 +1.150000e+02 2.550000e+02 2.542811e+02 0.000000e+00 2.000000e-03 +1.250000e+02 2.550000e+02 2.579790e+02 0.000000e+00 2.000000e-03 +1.350000e+02 2.550000e+02 2.615674e+02 0.000000e+00 2.000000e-03 +1.450000e+02 2.550000e+02 2.650117e+02 0.000000e+00 2.000000e-03 +1.550000e+02 2.550000e+02 2.682770e+02 0.000000e+00 2.000000e-03 +1.650000e+02 2.550000e+02 2.713296e+02 0.000000e+00 2.000000e-03 +1.750000e+02 2.550000e+02 2.741368e+02 0.000000e+00 2.000000e-03 +1.850000e+02 2.550000e+02 2.766679e+02 0.000000e+00 2.000000e-03 +1.950000e+02 2.550000e+02 2.788945e+02 0.000000e+00 2.000000e-03 +2.050000e+02 2.550000e+02 2.807913e+02 0.000000e+00 2.000000e-03 +2.150000e+02 2.550000e+02 2.823363e+02 0.000000e+00 2.000000e-03 +2.250000e+02 2.550000e+02 2.835114e+02 0.000000e+00 2.000000e-03 +2.350000e+02 2.550000e+02 2.843026e+02 0.000000e+00 2.000000e-03 +2.450000e+02 2.550000e+02 2.847007e+02 0.000000e+00 2.000000e-03 +2.550000e+02 2.550000e+02 2.847007e+02 0.000000e+00 2.000000e-03 +2.650000e+02 2.550000e+02 2.843026e+02 0.000000e+00 2.000000e-03 +2.750000e+02 2.550000e+02 2.835114e+02 0.000000e+00 2.000000e-03 +2.850000e+02 2.550000e+02 2.823363e+02 0.000000e+00 2.000000e-03 +2.950000e+02 2.550000e+02 2.807913e+02 0.000000e+00 2.000000e-03 +3.050000e+02 2.550000e+02 2.788945e+02 0.000000e+00 2.000000e-03 +3.150000e+02 2.550000e+02 2.766679e+02 0.000000e+00 2.000000e-03 +3.250000e+02 2.550000e+02 2.741368e+02 0.000000e+00 2.000000e-03 +3.350000e+02 2.550000e+02 2.713296e+02 0.000000e+00 2.000000e-03 +3.450000e+02 2.550000e+02 2.682770e+02 0.000000e+00 2.000000e-03 +3.550000e+02 2.550000e+02 2.650117e+02 0.000000e+00 2.000000e-03 +3.650000e+02 2.550000e+02 2.615674e+02 0.000000e+00 2.000000e-03 +3.750000e+02 2.550000e+02 2.579791e+02 0.000000e+00 2.000000e-03 +1.050000e+02 2.650000e+02 2.501165e+02 0.000000e+00 2.000000e-03 +1.150000e+02 2.650000e+02 2.538671e+02 0.000000e+00 2.000000e-03 +1.250000e+02 2.650000e+02 2.575429e+02 0.000000e+00 2.000000e-03 +1.350000e+02 2.650000e+02 2.611099e+02 0.000000e+00 2.000000e-03 +1.450000e+02 2.650000e+02 2.645335e+02 0.000000e+00 2.000000e-03 +1.550000e+02 2.650000e+02 2.677793e+02 0.000000e+00 2.000000e-03 +1.650000e+02 2.650000e+02 2.708137e+02 0.000000e+00 2.000000e-03 +1.750000e+02 2.650000e+02 2.736041e+02 0.000000e+00 2.000000e-03 +1.850000e+02 2.650000e+02 2.761201e+02 0.000000e+00 2.000000e-03 +1.950000e+02 2.650000e+02 2.783334e+02 0.000000e+00 2.000000e-03 +2.050000e+02 2.650000e+02 2.802189e+02 0.000000e+00 2.000000e-03 +2.150000e+02 2.650000e+02 2.817546e+02 0.000000e+00 2.000000e-03 +2.250000e+02 2.650000e+02 2.829227e+02 0.000000e+00 2.000000e-03 +2.350000e+02 2.650000e+02 2.837092e+02 0.000000e+00 2.000000e-03 +2.450000e+02 2.650000e+02 2.841048e+02 0.000000e+00 2.000000e-03 +2.550000e+02 2.650000e+02 2.841048e+02 0.000000e+00 2.000000e-03 +2.650000e+02 2.650000e+02 2.837092e+02 0.000000e+00 2.000000e-03 +2.750000e+02 2.650000e+02 2.829227e+02 0.000000e+00 2.000000e-03 +2.850000e+02 2.650000e+02 2.817546e+02 0.000000e+00 2.000000e-03 +2.950000e+02 2.650000e+02 2.802189e+02 0.000000e+00 2.000000e-03 +3.050000e+02 2.650000e+02 2.783334e+02 0.000000e+00 2.000000e-03 +3.150000e+02 2.650000e+02 2.761201e+02 0.000000e+00 2.000000e-03 +3.250000e+02 2.650000e+02 2.736041e+02 0.000000e+00 2.000000e-03 +3.350000e+02 2.650000e+02 2.708137e+02 0.000000e+00 2.000000e-03 +3.450000e+02 2.650000e+02 2.677793e+02 0.000000e+00 2.000000e-03 +3.550000e+02 2.650000e+02 2.645335e+02 0.000000e+00 2.000000e-03 +3.650000e+02 2.650000e+02 2.611099e+02 0.000000e+00 2.000000e-03 +3.750000e+02 2.650000e+02 2.575429e+02 0.000000e+00 2.000000e-03 +1.050000e+02 2.750000e+02 2.494692e+02 0.000000e+00 2.000000e-03 +1.150000e+02 2.750000e+02 2.531825e+02 0.000000e+00 2.000000e-03 +1.250000e+02 2.750000e+02 2.568218e+02 0.000000e+00 2.000000e-03 +1.350000e+02 2.750000e+02 2.603533e+02 0.000000e+00 2.000000e-03 +1.450000e+02 2.750000e+02 2.637429e+02 0.000000e+00 2.000000e-03 +1.550000e+02 2.750000e+02 2.669565e+02 0.000000e+00 2.000000e-03 +1.650000e+02 2.750000e+02 2.699606e+02 0.000000e+00 2.000000e-03 +1.750000e+02 2.750000e+02 2.727234e+02 0.000000e+00 2.000000e-03 +1.850000e+02 2.750000e+02 2.752143e+02 0.000000e+00 2.000000e-03 +1.950000e+02 2.750000e+02 2.774056e+02 0.000000e+00 2.000000e-03 +2.050000e+02 2.750000e+02 2.792724e+02 0.000000e+00 2.000000e-03 +2.150000e+02 2.750000e+02 2.807928e+02 0.000000e+00 2.000000e-03 +2.250000e+02 2.750000e+02 2.819493e+02 0.000000e+00 2.000000e-03 +2.350000e+02 2.750000e+02 2.827280e+02 0.000000e+00 2.000000e-03 +2.450000e+02 2.750000e+02 2.831197e+02 0.000000e+00 2.000000e-03 +2.550000e+02 2.750000e+02 2.831197e+02 0.000000e+00 2.000000e-03 +2.650000e+02 2.750000e+02 2.827280e+02 0.000000e+00 2.000000e-03 +2.750000e+02 2.750000e+02 2.819493e+02 0.000000e+00 2.000000e-03 +2.850000e+02 2.750000e+02 2.807928e+02 0.000000e+00 2.000000e-03 +2.950000e+02 2.750000e+02 2.792724e+02 0.000000e+00 2.000000e-03 +3.050000e+02 2.750000e+02 2.774056e+02 0.000000e+00 2.000000e-03 +3.150000e+02 2.750000e+02 2.752143e+02 0.000000e+00 2.000000e-03 +3.250000e+02 2.750000e+02 2.727234e+02 0.000000e+00 2.000000e-03 +3.350000e+02 2.750000e+02 2.699606e+02 0.000000e+00 2.000000e-03 +3.450000e+02 2.750000e+02 2.669565e+02 0.000000e+00 2.000000e-03 +3.550000e+02 2.750000e+02 2.637429e+02 0.000000e+00 2.000000e-03 +3.650000e+02 2.750000e+02 2.603533e+02 0.000000e+00 2.000000e-03 +3.750000e+02 2.750000e+02 2.568218e+02 0.000000e+00 2.000000e-03 +1.050000e+02 2.850000e+02 2.485739e+02 0.000000e+00 2.000000e-03 +1.150000e+02 2.850000e+02 2.522356e+02 0.000000e+00 2.000000e-03 +1.250000e+02 2.850000e+02 2.558243e+02 0.000000e+00 2.000000e-03 +1.350000e+02 2.850000e+02 2.593068e+02 0.000000e+00 2.000000e-03 +1.450000e+02 2.850000e+02 2.626493e+02 0.000000e+00 2.000000e-03 +1.550000e+02 2.850000e+02 2.658182e+02 0.000000e+00 2.000000e-03 +1.650000e+02 2.850000e+02 2.687807e+02 0.000000e+00 2.000000e-03 +1.750000e+02 2.850000e+02 2.715050e+02 0.000000e+00 2.000000e-03 +1.850000e+02 2.850000e+02 2.739614e+02 0.000000e+00 2.000000e-03 +1.950000e+02 2.850000e+02 2.761223e+02 0.000000e+00 2.000000e-03 +2.050000e+02 2.850000e+02 2.779631e+02 0.000000e+00 2.000000e-03 +2.150000e+02 2.850000e+02 2.794624e+02 0.000000e+00 2.000000e-03 +2.250000e+02 2.850000e+02 2.806028e+02 0.000000e+00 2.000000e-03 +2.350000e+02 2.850000e+02 2.813707e+02 0.000000e+00 2.000000e-03 +2.450000e+02 2.850000e+02 2.817570e+02 0.000000e+00 2.000000e-03 +2.550000e+02 2.850000e+02 2.817570e+02 0.000000e+00 2.000000e-03 +2.650000e+02 2.850000e+02 2.813707e+02 0.000000e+00 2.000000e-03 +2.750000e+02 2.850000e+02 2.806028e+02 0.000000e+00 2.000000e-03 +2.850000e+02 2.850000e+02 2.794624e+02 0.000000e+00 2.000000e-03 +2.950000e+02 2.850000e+02 2.779631e+02 0.000000e+00 2.000000e-03 +3.050000e+02 2.850000e+02 2.761223e+02 0.000000e+00 2.000000e-03 +3.150000e+02 2.850000e+02 2.739614e+02 0.000000e+00 2.000000e-03 +3.250000e+02 2.850000e+02 2.715050e+02 0.000000e+00 2.000000e-03 +3.350000e+02 2.850000e+02 2.687807e+02 0.000000e+00 2.000000e-03 +3.450000e+02 2.850000e+02 2.658182e+02 0.000000e+00 2.000000e-03 +3.550000e+02 2.850000e+02 2.626493e+02 0.000000e+00 2.000000e-03 +3.650000e+02 2.850000e+02 2.593068e+02 0.000000e+00 2.000000e-03 +3.750000e+02 2.850000e+02 2.558243e+02 0.000000e+00 2.000000e-03 +1.050000e+02 2.950000e+02 2.474409e+02 0.000000e+00 2.000000e-03 +1.150000e+02 2.950000e+02 2.510374e+02 0.000000e+00 2.000000e-03 +1.250000e+02 2.950000e+02 2.545621e+02 0.000000e+00 2.000000e-03 +1.350000e+02 2.950000e+02 2.579825e+02 0.000000e+00 2.000000e-03 +1.450000e+02 2.950000e+02 2.612655e+02 0.000000e+00 2.000000e-03 +1.550000e+02 2.950000e+02 2.643779e+02 0.000000e+00 2.000000e-03 +1.650000e+02 2.950000e+02 2.672876e+02 0.000000e+00 2.000000e-03 +1.750000e+02 2.950000e+02 2.699634e+02 0.000000e+00 2.000000e-03 +1.850000e+02 2.950000e+02 2.723760e+02 0.000000e+00 2.000000e-03 +1.950000e+02 2.950000e+02 2.744984e+02 0.000000e+00 2.000000e-03 +2.050000e+02 2.950000e+02 2.763063e+02 0.000000e+00 2.000000e-03 +2.150000e+02 2.950000e+02 2.777790e+02 0.000000e+00 2.000000e-03 +2.250000e+02 2.950000e+02 2.788990e+02 0.000000e+00 2.000000e-03 +2.350000e+02 2.950000e+02 2.796533e+02 0.000000e+00 2.000000e-03 +2.450000e+02 2.950000e+02 2.800326e+02 0.000000e+00 2.000000e-03 +2.550000e+02 2.950000e+02 2.800326e+02 0.000000e+00 2.000000e-03 +2.650000e+02 2.950000e+02 2.796533e+02 0.000000e+00 2.000000e-03 +2.750000e+02 2.950000e+02 2.788990e+02 0.000000e+00 2.000000e-03 +2.850000e+02 2.950000e+02 2.777790e+02 0.000000e+00 2.000000e-03 +2.950000e+02 2.950000e+02 2.763063e+02 0.000000e+00 2.000000e-03 +3.050000e+02 2.950000e+02 2.744984e+02 0.000000e+00 2.000000e-03 +3.150000e+02 2.950000e+02 2.723760e+02 0.000000e+00 2.000000e-03 +3.250000e+02 2.950000e+02 2.699634e+02 0.000000e+00 2.000000e-03 +3.350000e+02 2.950000e+02 2.672876e+02 0.000000e+00 2.000000e-03 +3.450000e+02 2.950000e+02 2.643779e+02 0.000000e+00 2.000000e-03 +3.550000e+02 2.950000e+02 2.612655e+02 0.000000e+00 2.000000e-03 +3.650000e+02 2.950000e+02 2.579825e+02 0.000000e+00 2.000000e-03 +3.750000e+02 2.950000e+02 2.545621e+02 0.000000e+00 2.000000e-03 +1.050000e+02 3.050000e+02 2.460835e+02 0.000000e+00 2.000000e-03 +1.150000e+02 3.050000e+02 2.496018e+02 0.000000e+00 2.000000e-03 +1.250000e+02 3.050000e+02 2.530500e+02 0.000000e+00 2.000000e-03 +1.350000e+02 3.050000e+02 2.563960e+02 0.000000e+00 2.000000e-03 +1.450000e+02 3.050000e+02 2.596076e+02 0.000000e+00 2.000000e-03 +1.550000e+02 3.050000e+02 2.626524e+02 0.000000e+00 2.000000e-03 +1.650000e+02 3.050000e+02 2.654988e+02 0.000000e+00 2.000000e-03 +1.750000e+02 3.050000e+02 2.681164e+02 0.000000e+00 2.000000e-03 +1.850000e+02 3.050000e+02 2.704766e+02 0.000000e+00 2.000000e-03 +1.950000e+02 3.050000e+02 2.725528e+02 0.000000e+00 2.000000e-03 +2.050000e+02 3.050000e+02 2.743215e+02 0.000000e+00 2.000000e-03 +2.150000e+02 3.050000e+02 2.757621e+02 0.000000e+00 2.000000e-03 +2.250000e+02 3.050000e+02 2.768578e+02 0.000000e+00 2.000000e-03 +2.350000e+02 3.050000e+02 2.775956e+02 0.000000e+00 2.000000e-03 +2.450000e+02 3.050000e+02 2.779668e+02 0.000000e+00 2.000000e-03 +2.550000e+02 3.050000e+02 2.779668e+02 0.000000e+00 2.000000e-03 +2.650000e+02 3.050000e+02 2.775956e+02 0.000000e+00 2.000000e-03 +2.750000e+02 3.050000e+02 2.768578e+02 0.000000e+00 2.000000e-03 +2.850000e+02 3.050000e+02 2.757621e+02 0.000000e+00 2.000000e-03 +2.950000e+02 3.050000e+02 2.743215e+02 0.000000e+00 2.000000e-03 +3.050000e+02 3.050000e+02 2.725528e+02 0.000000e+00 2.000000e-03 +3.150000e+02 3.050000e+02 2.704766e+02 0.000000e+00 2.000000e-03 +3.250000e+02 3.050000e+02 2.681164e+02 0.000000e+00 2.000000e-03 +3.350000e+02 3.050000e+02 2.654988e+02 0.000000e+00 2.000000e-03 +3.450000e+02 3.050000e+02 2.626524e+02 0.000000e+00 2.000000e-03 +3.550000e+02 3.050000e+02 2.596076e+02 0.000000e+00 2.000000e-03 +3.650000e+02 3.050000e+02 2.563960e+02 0.000000e+00 2.000000e-03 +3.750000e+02 3.050000e+02 2.530500e+02 0.000000e+00 2.000000e-03 +1.050000e+02 3.150000e+02 2.445174e+02 0.000000e+00 2.000000e-03 +1.150000e+02 3.150000e+02 2.479454e+02 0.000000e+00 2.000000e-03 +1.250000e+02 3.150000e+02 2.513052e+02 0.000000e+00 2.000000e-03 +1.350000e+02 3.150000e+02 2.545655e+02 0.000000e+00 2.000000e-03 +1.450000e+02 3.150000e+02 2.576947e+02 0.000000e+00 2.000000e-03 +1.550000e+02 3.150000e+02 2.606614e+02 0.000000e+00 2.000000e-03 +1.650000e+02 3.150000e+02 2.634348e+02 0.000000e+00 2.000000e-03 +1.750000e+02 3.150000e+02 2.659854e+02 0.000000e+00 2.000000e-03 +1.850000e+02 3.150000e+02 2.682850e+02 0.000000e+00 2.000000e-03 +1.950000e+02 3.150000e+02 2.703080e+02 0.000000e+00 2.000000e-03 +2.050000e+02 3.150000e+02 2.720313e+02 0.000000e+00 2.000000e-03 +2.150000e+02 3.150000e+02 2.734350e+02 0.000000e+00 2.000000e-03 +2.250000e+02 3.150000e+02 2.745027e+02 0.000000e+00 2.000000e-03 +2.350000e+02 3.150000e+02 2.752216e+02 0.000000e+00 2.000000e-03 +2.450000e+02 3.150000e+02 2.755831e+02 0.000000e+00 2.000000e-03 +2.550000e+02 3.150000e+02 2.755831e+02 0.000000e+00 2.000000e-03 +2.650000e+02 3.150000e+02 2.752216e+02 0.000000e+00 2.000000e-03 +2.750000e+02 3.150000e+02 2.745027e+02 0.000000e+00 2.000000e-03 +2.850000e+02 3.150000e+02 2.734350e+02 0.000000e+00 2.000000e-03 +2.950000e+02 3.150000e+02 2.720313e+02 0.000000e+00 2.000000e-03 +3.050000e+02 3.150000e+02 2.703080e+02 0.000000e+00 2.000000e-03 +3.150000e+02 3.150000e+02 2.682850e+02 0.000000e+00 2.000000e-03 +3.250000e+02 3.150000e+02 2.659854e+02 0.000000e+00 2.000000e-03 +3.350000e+02 3.150000e+02 2.634348e+02 0.000000e+00 2.000000e-03 +3.450000e+02 3.150000e+02 2.606614e+02 0.000000e+00 2.000000e-03 +3.550000e+02 3.150000e+02 2.576947e+02 0.000000e+00 2.000000e-03 +3.650000e+02 3.150000e+02 2.545655e+02 0.000000e+00 2.000000e-03 +3.750000e+02 3.150000e+02 2.513052e+02 0.000000e+00 2.000000e-03 +1.050000e+02 3.250000e+02 2.427601e+02 0.000000e+00 2.000000e-03 +1.150000e+02 3.250000e+02 2.460869e+02 0.000000e+00 2.000000e-03 +1.250000e+02 3.250000e+02 2.493475e+02 0.000000e+00 2.000000e-03 +1.350000e+02 3.250000e+02 2.525115e+02 0.000000e+00 2.000000e-03 +1.450000e+02 3.250000e+02 2.555483e+02 0.000000e+00 2.000000e-03 +1.550000e+02 3.250000e+02 2.584275e+02 0.000000e+00 2.000000e-03 +1.650000e+02 3.250000e+02 2.611190e+02 0.000000e+00 2.000000e-03 +1.750000e+02 3.250000e+02 2.635942e+02 0.000000e+00 2.000000e-03 +1.850000e+02 3.250000e+02 2.658260e+02 0.000000e+00 2.000000e-03 +1.950000e+02 3.250000e+02 2.677892e+02 0.000000e+00 2.000000e-03 +2.050000e+02 3.250000e+02 2.694617e+02 0.000000e+00 2.000000e-03 +2.150000e+02 3.250000e+02 2.708240e+02 0.000000e+00 2.000000e-03 +2.250000e+02 3.250000e+02 2.718600e+02 0.000000e+00 2.000000e-03 +2.350000e+02 3.250000e+02 2.725577e+02 0.000000e+00 2.000000e-03 +2.450000e+02 3.250000e+02 2.729086e+02 0.000000e+00 2.000000e-03 +2.550000e+02 3.250000e+02 2.729086e+02 0.000000e+00 2.000000e-03 +2.650000e+02 3.250000e+02 2.725577e+02 0.000000e+00 2.000000e-03 +2.750000e+02 3.250000e+02 2.718600e+02 0.000000e+00 2.000000e-03 +2.850000e+02 3.250000e+02 2.708240e+02 0.000000e+00 2.000000e-03 +2.950000e+02 3.250000e+02 2.694617e+02 0.000000e+00 2.000000e-03 +3.050000e+02 3.250000e+02 2.677892e+02 0.000000e+00 2.000000e-03 +3.150000e+02 3.250000e+02 2.658260e+02 0.000000e+00 2.000000e-03 +3.250000e+02 3.250000e+02 2.635942e+02 0.000000e+00 2.000000e-03 +3.350000e+02 3.250000e+02 2.611190e+02 0.000000e+00 2.000000e-03 +3.450000e+02 3.250000e+02 2.584275e+02 0.000000e+00 2.000000e-03 +3.550000e+02 3.250000e+02 2.555483e+02 0.000000e+00 2.000000e-03 +3.650000e+02 3.250000e+02 2.525115e+02 0.000000e+00 2.000000e-03 +3.750000e+02 3.250000e+02 2.493475e+02 0.000000e+00 2.000000e-03 +1.050000e+02 3.350000e+02 2.408312e+02 0.000000e+00 2.000000e-03 +1.150000e+02 3.350000e+02 2.440469e+02 0.000000e+00 2.000000e-03 +1.250000e+02 3.350000e+02 2.471986e+02 0.000000e+00 2.000000e-03 +1.350000e+02 3.350000e+02 2.502569e+02 0.000000e+00 2.000000e-03 +1.450000e+02 3.350000e+02 2.531923e+02 0.000000e+00 2.000000e-03 +1.550000e+02 3.350000e+02 2.559753e+02 0.000000e+00 2.000000e-03 +1.650000e+02 3.350000e+02 2.585770e+02 0.000000e+00 2.000000e-03 +1.750000e+02 3.350000e+02 2.609695e+02 0.000000e+00 2.000000e-03 +1.850000e+02 3.350000e+02 2.631267e+02 0.000000e+00 2.000000e-03 +1.950000e+02 3.350000e+02 2.650244e+02 0.000000e+00 2.000000e-03 +2.050000e+02 3.350000e+02 2.666411e+02 0.000000e+00 2.000000e-03 +2.150000e+02 3.350000e+02 2.679578e+02 0.000000e+00 2.000000e-03 +2.250000e+02 3.350000e+02 2.689593e+02 0.000000e+00 2.000000e-03 +2.350000e+02 3.350000e+02 2.696337e+02 0.000000e+00 2.000000e-03 +2.450000e+02 3.350000e+02 2.699729e+02 0.000000e+00 2.000000e-03 +2.550000e+02 3.350000e+02 2.699729e+02 0.000000e+00 2.000000e-03 +2.650000e+02 3.350000e+02 2.696337e+02 0.000000e+00 2.000000e-03 +2.750000e+02 3.350000e+02 2.689593e+02 0.000000e+00 2.000000e-03 +2.850000e+02 3.350000e+02 2.679578e+02 0.000000e+00 2.000000e-03 +2.950000e+02 3.350000e+02 2.666411e+02 0.000000e+00 2.000000e-03 +3.050000e+02 3.350000e+02 2.650244e+02 0.000000e+00 2.000000e-03 +3.150000e+02 3.350000e+02 2.631267e+02 0.000000e+00 2.000000e-03 +3.250000e+02 3.350000e+02 2.609695e+02 0.000000e+00 2.000000e-03 +3.350000e+02 3.350000e+02 2.585770e+02 0.000000e+00 2.000000e-03 +3.450000e+02 3.350000e+02 2.559753e+02 0.000000e+00 2.000000e-03 +3.550000e+02 3.350000e+02 2.531923e+02 0.000000e+00 2.000000e-03 +3.650000e+02 3.350000e+02 2.502569e+02 0.000000e+00 2.000000e-03 +3.750000e+02 3.350000e+02 2.471986e+02 0.000000e+00 2.000000e-03 +1.050000e+02 3.450000e+02 2.387514e+02 0.000000e+00 2.000000e-03 +1.150000e+02 3.450000e+02 2.418474e+02 0.000000e+00 2.000000e-03 +1.250000e+02 3.450000e+02 2.448817e+02 0.000000e+00 2.000000e-03 +1.350000e+02 3.450000e+02 2.478261e+02 0.000000e+00 2.000000e-03 +1.450000e+02 3.450000e+02 2.506521e+02 0.000000e+00 2.000000e-03 +1.550000e+02 3.450000e+02 2.533315e+02 0.000000e+00 2.000000e-03 +1.650000e+02 3.450000e+02 2.558362e+02 0.000000e+00 2.000000e-03 +1.750000e+02 3.450000e+02 2.581396e+02 0.000000e+00 2.000000e-03 +1.850000e+02 3.450000e+02 2.602165e+02 0.000000e+00 2.000000e-03 +1.950000e+02 3.450000e+02 2.620435e+02 0.000000e+00 2.000000e-03 +2.050000e+02 3.450000e+02 2.635999e+02 0.000000e+00 2.000000e-03 +2.150000e+02 3.450000e+02 2.648676e+02 0.000000e+00 2.000000e-03 +2.250000e+02 3.450000e+02 2.658318e+02 0.000000e+00 2.000000e-03 +2.350000e+02 3.450000e+02 2.664810e+02 0.000000e+00 2.000000e-03 +2.450000e+02 3.450000e+02 2.668076e+02 0.000000e+00 2.000000e-03 +2.550000e+02 3.450000e+02 2.668076e+02 0.000000e+00 2.000000e-03 +2.650000e+02 3.450000e+02 2.664810e+02 0.000000e+00 2.000000e-03 +2.750000e+02 3.450000e+02 2.658318e+02 0.000000e+00 2.000000e-03 +2.850000e+02 3.450000e+02 2.648676e+02 0.000000e+00 2.000000e-03 +2.950000e+02 3.450000e+02 2.635999e+02 0.000000e+00 2.000000e-03 +3.050000e+02 3.450000e+02 2.620435e+02 0.000000e+00 2.000000e-03 +3.150000e+02 3.450000e+02 2.602165e+02 0.000000e+00 2.000000e-03 +3.250000e+02 3.450000e+02 2.581396e+02 0.000000e+00 2.000000e-03 +3.350000e+02 3.450000e+02 2.558362e+02 0.000000e+00 2.000000e-03 +3.450000e+02 3.450000e+02 2.533315e+02 0.000000e+00 2.000000e-03 +3.550000e+02 3.450000e+02 2.506521e+02 0.000000e+00 2.000000e-03 +3.650000e+02 3.450000e+02 2.478261e+02 0.000000e+00 2.000000e-03 +3.750000e+02 3.450000e+02 2.448817e+02 0.000000e+00 2.000000e-03 +1.050000e+02 3.550000e+02 2.365428e+02 0.000000e+00 2.000000e-03 +1.150000e+02 3.550000e+02 2.395115e+02 0.000000e+00 2.000000e-03 +1.250000e+02 3.550000e+02 2.424211e+02 0.000000e+00 2.000000e-03 +1.350000e+02 3.550000e+02 2.452446e+02 0.000000e+00 2.000000e-03 +1.450000e+02 3.550000e+02 2.479545e+02 0.000000e+00 2.000000e-03 +1.550000e+02 3.550000e+02 2.505237e+02 0.000000e+00 2.000000e-03 +1.650000e+02 3.550000e+02 2.529255e+02 0.000000e+00 2.000000e-03 +1.750000e+02 3.550000e+02 2.551343e+02 0.000000e+00 2.000000e-03 +1.850000e+02 3.550000e+02 2.571258e+02 0.000000e+00 2.000000e-03 +1.950000e+02 3.550000e+02 2.588778e+02 0.000000e+00 2.000000e-03 +2.050000e+02 3.550000e+02 2.603702e+02 0.000000e+00 2.000000e-03 +2.150000e+02 3.550000e+02 2.615858e+02 0.000000e+00 2.000000e-03 +2.250000e+02 3.550000e+02 2.625104e+02 0.000000e+00 2.000000e-03 +2.350000e+02 3.550000e+02 2.631330e+02 0.000000e+00 2.000000e-03 +2.450000e+02 3.550000e+02 2.634461e+02 0.000000e+00 2.000000e-03 +2.550000e+02 3.550000e+02 2.634461e+02 0.000000e+00 2.000000e-03 +2.650000e+02 3.550000e+02 2.631330e+02 0.000000e+00 2.000000e-03 +2.750000e+02 3.550000e+02 2.625104e+02 0.000000e+00 2.000000e-03 +2.850000e+02 3.550000e+02 2.615858e+02 0.000000e+00 2.000000e-03 +2.950000e+02 3.550000e+02 2.603702e+02 0.000000e+00 2.000000e-03 +3.050000e+02 3.550000e+02 2.588778e+02 0.000000e+00 2.000000e-03 +3.150000e+02 3.550000e+02 2.571258e+02 0.000000e+00 2.000000e-03 +3.250000e+02 3.550000e+02 2.551343e+02 0.000000e+00 2.000000e-03 +3.350000e+02 3.550000e+02 2.529255e+02 0.000000e+00 2.000000e-03 +3.450000e+02 3.550000e+02 2.505237e+02 0.000000e+00 2.000000e-03 +3.550000e+02 3.550000e+02 2.479545e+02 0.000000e+00 2.000000e-03 +3.650000e+02 3.550000e+02 2.452446e+02 0.000000e+00 2.000000e-03 +3.750000e+02 3.550000e+02 2.424211e+02 0.000000e+00 2.000000e-03 +1.050000e+02 3.650000e+02 2.342278e+02 0.000000e+00 2.000000e-03 +1.150000e+02 3.650000e+02 2.370632e+02 0.000000e+00 2.000000e-03 +1.250000e+02 3.650000e+02 2.398421e+02 0.000000e+00 2.000000e-03 +1.350000e+02 3.650000e+02 2.425387e+02 0.000000e+00 2.000000e-03 +1.450000e+02 3.650000e+02 2.451270e+02 0.000000e+00 2.000000e-03 +1.550000e+02 3.650000e+02 2.475808e+02 0.000000e+00 2.000000e-03 +1.650000e+02 3.650000e+02 2.498748e+02 0.000000e+00 2.000000e-03 +1.750000e+02 3.650000e+02 2.519843e+02 0.000000e+00 2.000000e-03 +1.850000e+02 3.650000e+02 2.538864e+02 0.000000e+00 2.000000e-03 +1.950000e+02 3.650000e+02 2.555596e+02 0.000000e+00 2.000000e-03 +2.050000e+02 3.650000e+02 2.569850e+02 0.000000e+00 2.000000e-03 +2.150000e+02 3.650000e+02 2.581461e+02 0.000000e+00 2.000000e-03 +2.250000e+02 3.650000e+02 2.590291e+02 0.000000e+00 2.000000e-03 +2.350000e+02 3.650000e+02 2.596237e+02 0.000000e+00 2.000000e-03 +2.450000e+02 3.650000e+02 2.599228e+02 0.000000e+00 2.000000e-03 +2.550000e+02 3.650000e+02 2.599228e+02 0.000000e+00 2.000000e-03 +2.650000e+02 3.650000e+02 2.596237e+02 0.000000e+00 2.000000e-03 +2.750000e+02 3.650000e+02 2.590291e+02 0.000000e+00 2.000000e-03 +2.850000e+02 3.650000e+02 2.581461e+02 0.000000e+00 2.000000e-03 +2.950000e+02 3.650000e+02 2.569850e+02 0.000000e+00 2.000000e-03 +3.050000e+02 3.650000e+02 2.555596e+02 0.000000e+00 2.000000e-03 +3.150000e+02 3.650000e+02 2.538864e+02 0.000000e+00 2.000000e-03 +3.250000e+02 3.650000e+02 2.519843e+02 0.000000e+00 2.000000e-03 +3.350000e+02 3.650000e+02 2.498748e+02 0.000000e+00 2.000000e-03 +3.450000e+02 3.650000e+02 2.475808e+02 0.000000e+00 2.000000e-03 +3.550000e+02 3.650000e+02 2.451270e+02 0.000000e+00 2.000000e-03 +3.650000e+02 3.650000e+02 2.425387e+02 0.000000e+00 2.000000e-03 +3.750000e+02 3.650000e+02 2.398421e+02 0.000000e+00 2.000000e-03 diff --git a/simpegPF/Dev/MAG/Intgrl_MAG_Inv_Driver.py b/simpegPF/Dev/MAG/Intgrl_MAG_Inv_Driver.py index 65e47f77..7f6b64f0 100644 --- a/simpegPF/Dev/MAG/Intgrl_MAG_Inv_Driver.py +++ b/simpegPF/Dev/MAG/Intgrl_MAG_Inv_Driver.py @@ -14,7 +14,6 @@ wd = survey.std ndata = survey.srcField.rxList[0].locs.shape[0] -beta_in = 1e+5 eps_p = 1e-4 eps_q = 1e-4 @@ -24,8 +23,8 @@ nC = len(actv) # Create active map to go from reduce set to full actvMap = Maps.InjectActiveCells(mesh, actv, -100) -# Creat reduced identity map -idenMap = Maps.IdentityMap(nP = nC) +# Create reduced identity map +idenMap = Maps.IdentityMap(nP=nC) # Get magnetization vector for MOF M_xyz = driver.magnetizationModel @@ -34,12 +33,6 @@ M_xyz = driver.magnetizationModel midx = int(mesh.nCx/2) midy = int(mesh.nCy/2) -# Get distance weighting function -#============================================================================== -# wr = PF.Magnetics.get_dist_wgt(mesh,rxLoc,actv,3.,np.min(mesh.hx)/4) -# #wrMap = PF.BaseMag.WeightMap(nC, wr) -#============================================================================== - #%% Plot obs data PF.Magnetics.plot_obs_2D(rxLoc,d, 'Observed Data') @@ -47,9 +40,7 @@ PF.Magnetics.plot_obs_2D(rxLoc,d, 'Observed Data') prob = PF.Magnetics.Problem3D_Integral(mesh, mapping=idenMap, actInd=actv) prob.solverOpts['accuracyTol'] = 1e-4 survey.pair(prob) -#survey.makeSyntheticData(data, std=0.01) -#survey.dobs=d -#survey.mtrue = model + # Write out the predicted pred = prob.fields(driver.m0) PF.Magnetics.writeUBCobs('Pred.dat', survey, pred) @@ -68,24 +59,19 @@ plt.gca().set_aspect('equal', adjustable='box') reg = Regularization.Simple(mesh, indActive=actv, mapping=idenMap) reg.mref = driver.mref reg.wght = wr -#reg.alpha_s = 1. - -# Create pre-conditioner -diagA = np.sum(prob.G**2.,axis=0) + beta_in*(reg.W.T*reg.W).diagonal() -PC = Utils.sdiag(diagA**-1.) dmis = DataMisfit.l2_DataMisfit(survey) dmis.Wd = 1/wd opt = Optimization.ProjectedGNCG(maxIter=10,lower=0.,upper=1., maxIterCG= 20, tolCG=1e-3) -opt.approxHinv = PC -# opt = Optimization.InexactGaussNewton(maxIter=6) -invProb = InvProblem.BaseInvProblem(dmis, reg, opt, beta=beta_in) +invProb = InvProblem.BaseInvProblem(dmis, reg, opt) + beta = Directives.BetaSchedule(coolingFactor=2, coolingRate=1) -#betaest = Directives.BetaEstimate_ByEig() +beta_init = Directives.BetaEstimate_ByEig() target = Directives.TargetMisfit() +update_Jacobi = Directives.Update_lin_PreCond(onlyOnStart=True) -inv = Inversion.BaseInversion(invProb, directiveList=[beta,target]) +inv = Inversion.BaseInversion(invProb, directiveList=[beta,target,beta_init,update_Jacobi]) # Run inversion mrec = inv.run(driver.m0) @@ -94,7 +80,6 @@ m_out = actvMap*mrec # Write result Mesh.TensorMesh.writeModelUBC(mesh,'SimPEG_inv_l2l2.sus',m_out) -#Utils.meshutils.writeUBCTensorModel(home_dir+dsep+'wr.dat',mesh,wr_out) # Plot predicted pred = prob.fields(mrec) @@ -134,6 +119,7 @@ plt.gca().set_aspect('equal', adjustable='box') phim = invProb.phi_m_last phid = invProb.phi_d +# Set parameters for sparsity reg = Regularization.Sparse(mesh, indActive = actv, mapping = idenMap) reg.recModel = mrec reg.mref = driver.mref @@ -143,26 +129,22 @@ reg.eps_q = eps_q reg.norms = driver.lpnorms -diagA = np.sum(prob.G**2.,axis=0) + beta_in*(reg.W.T*reg.W).diagonal() -PC = Utils.sdiag(diagA**-1.) - -#reg.alpha_s = 1. - dmis = DataMisfit.l2_DataMisfit(survey) dmis.Wd = wd opt = Optimization.ProjectedGNCG(maxIter=20 ,lower=0.,upper=1., maxIterCG= 10, tolCG = 1e-4) -opt.approxHinv = PC -#opt.phim_last = reg.eval(mrec) -# opt = Optimization.InexactGaussNewton(maxIter=6) invProb = InvProblem.BaseInvProblem(dmis, reg, opt, beta = invProb.beta) + +# Create inversion directives beta = Directives.BetaSchedule(coolingFactor=1, coolingRate=1) update_beta = Directives.Scale_Beta(tol = 0.05) -#betaest = Directives.BetaEstimate_ByEig() target = Directives.TargetMisfit() IRLS =Directives.Update_IRLS( phi_m_last = phim, phi_d_last = phid ) +update_Jacobi = Directives.Update_lin_PreCond(onlyOnStart=False) +save_log = Directives.SaveOutputEveryIteration() +save_log.fileName = 'LogName_blabla' -inv = Inversion.BaseInversion(invProb, directiveList=[beta,IRLS,update_beta]) +inv = Inversion.BaseInversion(invProb, directiveList=[beta,IRLS,update_beta,update_Jacobi,save_log]) m0 = mrec @@ -171,6 +153,7 @@ mrec = inv.run(m0) m_out = actvMap*mrec +# Write final model out. Mesh.TensorMesh.writeModelUBC(mesh,'SimPEG_inv_l0l2.sus',m_out) pred = prob.fields(mrec) diff --git a/simpegPF/notebooks/SimPEG Tutorial - MAG Linear Problem.ipynb b/simpegPF/notebooks/SimPEG Tutorial - MAG Linear Problem.ipynb index d504509e..b2512692 100644 --- a/simpegPF/notebooks/SimPEG Tutorial - MAG Linear Problem.ipynb +++ b/simpegPF/notebooks/SimPEG Tutorial - MAG Linear Problem.ipynb @@ -43,7 +43,7 @@ }, { "cell_type": "code", - "execution_count": 4, + "execution_count": 1, "metadata": { "collapsed": false }, @@ -53,15 +53,19 @@ "output_type": "stream", "text": [ "Using matplotlib backend: nbAgg\n", - "Populating the interactive namespace from numpy and matplotlib\n" + "Populating the interactive namespace from numpy and matplotlib\n", + "Efficiency Warning: Interpolation will be slow, use setup.py!\n", + "\n", + " python setup.py build_ext --inplace\n", + " \n" ] }, { "name": "stderr", "output_type": "stream", "text": [ - "WARNING: pylab import has clobbered these variables: ['linalg']\n", - "`%matplotlib` prevents importing * from pylab and numpy\n" + "C:\\Users\\dominiquef.MIRAGEOSCIENCE\\AppData\\Local\\Continuum\\Anaconda\\lib\\site-packages\\IPython\\kernel\\__init__.py:13: ShimWarning: The `IPython.kernel` package has been deprecated. You should import from ipykernel or jupyter_client instead.\n", + " \"You should import from ipykernel or jupyter_client instead.\", ShimWarning)\n" ] } ], @@ -75,7 +79,7 @@ }, { "cell_type": "code", - "execution_count": 5, + "execution_count": 2, "metadata": { "collapsed": false, "scrolled": true @@ -176,7 +180,7 @@ }, { "cell_type": "code", - "execution_count": 7, + "execution_count": 3, "metadata": { "collapsed": false }, @@ -227,7 +231,7 @@ }, { "cell_type": "code", - "execution_count": 8, + "execution_count": 4, "metadata": { "collapsed": false }, @@ -971,7 +975,7 @@ { "data": { "text/html": [ - "" + "" ], "text/plain": [ "" @@ -983,10 +987,10 @@ { "data": { "text/plain": [ - "" + "" ] }, - "execution_count": 8, + "execution_count": 4, "metadata": {}, "output_type": "execute_result" } @@ -1022,7 +1026,7 @@ }, { "cell_type": "code", - "execution_count": 25, + "execution_count": 5, "metadata": { "collapsed": false }, @@ -1800,7 +1804,7 @@ }, { "cell_type": "code", - "execution_count": 26, + "execution_count": 10, "metadata": { "collapsed": false }, @@ -1820,7 +1824,7 @@ "beta_in = 1e+4\n", "\n", "# Create a jacobi pre-conditioner\n", - "diagA = np.sum(prob.G**2.,axis=0) + beta_in*(reg.W.T*reg.W).diagonal()*wr\n", + "diagA = np.sum(prob.G**2.,axis=0) + beta_in*(reg.W.T*reg.W).diagonal()\n", "PC = Utils.sdiag(diagA**-1.)\n", "\n", "# Define the misfit function\n", @@ -1843,7 +1847,7 @@ }, { "cell_type": "code", - "execution_count": 27, + "execution_count": 11, "metadata": { "collapsed": false }, @@ -1858,16 +1862,15 @@ "=============================== Projected GNCG ===============================\n", " # beta phi_d phi_m f |proj(x-g)-x| LS Comment \n", "-----------------------------------------------------------------------------\n", - " 0 1.00e+04 6.70e+04 0.00e+00 6.70e+04 8.30e+01 0 \n", - " 1 5.00e+03 2.94e+03 4.20e-02 3.15e+03 6.89e+01 0 \n", - " 2 2.50e+03 3.61e+02 6.08e-02 5.13e+02 7.04e+01 0 \n", - " 3 1.25e+03 2.58e+02 5.84e-02 3.31e+02 8.36e+01 0 Skip BFGS \n", + " 0 1.00e+04 6.73e+04 0.00e+00 6.73e+04 8.30e+01 0 \n", + " 1 5.00e+03 1.33e+03 3.52e-02 1.50e+03 6.45e+01 0 \n", + " 2 2.50e+03 3.14e+02 5.17e-02 4.44e+02 7.06e+01 0 Skip BFGS \n", "------------------------- STOP! -------------------------\n", - "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 6.6998e+03\n", - "1 : |xc-x_last| = 8.5475e-03 <= tolX*(1+|x0|) = 1.0111e-01\n", - "0 : |proj(x-g)-x| = 8.3598e+01 <= tolG = 1.0000e-01\n", - "0 : |proj(x-g)-x| = 8.3598e+01 <= 1e3*eps = 1.0000e-02\n", - "0 : maxIter = 10 <= iter = 4\n", + "1 : |fc-fOld| = 0.0000e+00 <= tolF*(1+|f0|) = 6.7288e+03\n", + "1 : |xc-x_last| = 1.6646e-02 <= tolX*(1+|x0|) = 1.0111e-01\n", + "0 : |proj(x-g)-x| = 7.0593e+01 <= tolG = 1.0000e-01\n", + "0 : |proj(x-g)-x| = 7.0593e+01 <= 1e3*eps = 1.0000e-02\n", + "0 : maxIter = 10 <= iter = 3\n", "------------------------- DONE! -------------------------\n" ] } @@ -1885,7 +1888,7 @@ }, { "cell_type": "code", - "execution_count": 28, + "execution_count": 12, "metadata": { "collapsed": false }, @@ -2629,7 +2632,7 @@ { "data": { "text/html": [ - "" + "" ], "text/plain": [ "" @@ -2669,7 +2672,7 @@ }, { "cell_type": "code", - "execution_count": 29, + "execution_count": 9, "metadata": { "collapsed": false }, @@ -3413,7 +3416,7 @@ { "data": { "text/html": [ - "" + "" ], "text/plain": [ "" @@ -3425,10 +3428,10 @@ { "data": { "text/plain": [ - "" + "" ] }, - "execution_count": 29, + "execution_count": 9, "metadata": {}, "output_type": "execute_result" }