From ba16e5e19477a2dc1ce082ddbbd9e6a01e172521 Mon Sep 17 00:00:00 2001 From: caes Date: Thu, 10 Aug 2017 04:05:56 -0400 Subject: [PATCH] sed success! --- ...rz1997.tab => magdziarz1997.tab.extracted} | 0 sed/magdziarz1997/magdziarz1997_samples.tab | 5 +- sed/mehdipour2013/mehdipour2013.tab | 2000 ++++++++--------- .../mehdipour2013_powlaw_coords.tab | 1 - src/agn.hpp | 2 +- src/sed.hpp | 44 +- src/table_powerlaw_spline_sed.cpp | 12 +- 7 files changed, 1027 insertions(+), 1037 deletions(-) rename sed/magdziarz1997/{magdziarz1997.tab => magdziarz1997.tab.extracted} (100%) diff --git a/sed/magdziarz1997/magdziarz1997.tab b/sed/magdziarz1997/magdziarz1997.tab.extracted similarity index 100% rename from sed/magdziarz1997/magdziarz1997.tab rename to sed/magdziarz1997/magdziarz1997.tab.extracted diff --git a/sed/magdziarz1997/magdziarz1997_samples.tab b/sed/magdziarz1997/magdziarz1997_samples.tab index 701a8bd..0322d0d 100644 --- a/sed/magdziarz1997/magdziarz1997_samples.tab +++ b/sed/magdziarz1997/magdziarz1997_samples.tab @@ -1,8 +1,9 @@ -0.0003 0.002 +#0.001 0.007 +#0.0003 0.002 0.001 0.0072 0.008 0.04 0.012 0.044 -0.03 0.025 +#0.03 0.025 0.44 0.024 2.1 0.012 60.0 0.029 diff --git a/sed/mehdipour2013/mehdipour2013.tab b/sed/mehdipour2013/mehdipour2013.tab index 058cbc5..c08a190 100644 --- a/sed/mehdipour2013/mehdipour2013.tab +++ b/sed/mehdipour2013/mehdipour2013.tab @@ -1,1000 +1,1000 @@ -1.36193e-07 4.87909e+13 -1.40936e-07 4.87915e+13 -1.45843e-07 4.87922e+13 -1.50922e-07 4.87929e+13 -1.56178e-07 4.87935e+13 -1.61616e-07 4.87942e+13 -1.67244e-07 4.87949e+13 -1.73068e-07 4.87955e+13 -1.79095e-07 4.87962e+13 -1.85332e-07 4.87969e+13 -1.91785e-07 4.87975e+13 -1.98464e-07 4.87982e+13 -2.05375e-07 4.87989e+13 -2.12527e-07 4.87995e+13 -2.19928e-07 4.88002e+13 -2.27586e-07 4.88009e+13 -2.35512e-07 4.88016e+13 -2.43713e-07 4.88022e+13 -2.52200e-07 4.88029e+13 -2.60982e-07 4.88036e+13 -2.70070e-07 4.88042e+13 -2.79475e-07 4.88049e+13 -2.89207e-07 4.88056e+13 -2.99278e-07 4.88062e+13 -3.09700e-07 4.88069e+13 -3.20484e-07 4.88076e+13 -3.31645e-07 4.88082e+13 -3.43194e-07 4.88089e+13 -3.55145e-07 4.88096e+13 -3.67512e-07 4.88102e+13 -3.80310e-07 4.88109e+13 -3.93553e-07 4.88116e+13 -4.07258e-07 4.88122e+13 -4.21440e-07 4.88129e+13 -4.36116e-07 4.88136e+13 -4.51303e-07 4.88142e+13 -4.67018e-07 4.88149e+13 -4.83281e-07 4.88156e+13 -5.00111e-07 4.88163e+13 -5.17526e-07 4.88169e+13 -5.35548e-07 4.88176e+13 -5.54197e-07 4.88183e+13 -5.73496e-07 4.88189e+13 -5.93467e-07 4.88196e+13 -6.14133e-07 4.88203e+13 -6.35519e-07 4.88209e+13 -6.57650e-07 4.88216e+13 -6.80551e-07 4.88223e+13 -7.04250e-07 4.88229e+13 -7.28774e-07 4.88236e+13 -7.54153e-07 4.88243e+13 -7.80414e-07 4.88249e+13 -8.07591e-07 4.88256e+13 -8.35714e-07 4.88263e+13 -8.64816e-07 4.88269e+13 -8.94931e-07 4.88276e+13 -9.26095e-07 4.88283e+13 -9.58345e-07 4.88290e+13 -9.91717e-07 4.88296e+13 -1.02625e-06 4.88303e+13 -1.06199e-06 4.88310e+13 -1.09897e-06 4.88316e+13 -1.13724e-06 4.88323e+13 -1.17684e-06 4.88330e+13 -1.21782e-06 4.88336e+13 -1.26023e-06 4.88343e+13 -1.30412e-06 4.88350e+13 -1.34953e-06 4.88356e+13 -1.39653e-06 4.88363e+13 -1.44516e-06 4.88370e+13 -1.49548e-06 4.88376e+13 -1.54756e-06 4.88383e+13 -1.60145e-06 4.88390e+13 -1.65722e-06 4.88397e+13 -1.71493e-06 4.88403e+13 -1.77464e-06 4.88410e+13 -1.83644e-06 4.88417e+13 -1.90039e-06 4.88423e+13 -1.96657e-06 4.88430e+13 -2.03505e-06 4.88437e+13 -2.10592e-06 4.88443e+13 -2.17925e-06 4.88450e+13 -2.25514e-06 4.88457e+13 -2.33367e-06 4.88463e+13 -2.41494e-06 4.88470e+13 -2.49904e-06 4.88477e+13 -2.58606e-06 4.88483e+13 -2.67611e-06 4.88490e+13 -2.76930e-06 4.88497e+13 -2.86574e-06 4.88504e+13 -2.96553e-06 4.88510e+13 -3.06880e-06 4.88517e+13 -3.17567e-06 4.88524e+13 -3.28625e-06 4.88530e+13 -3.40069e-06 4.88537e+13 -3.51911e-06 4.88544e+13 -3.64166e-06 4.88550e+13 -3.76847e-06 4.88557e+13 -3.89970e-06 4.88564e+13 -4.03550e-06 4.88570e+13 -4.17603e-06 4.88577e+13 -4.32145e-06 4.88584e+13 -4.47194e-06 4.88590e+13 -4.62767e-06 4.88597e+13 -4.78881e-06 4.88604e+13 -4.95558e-06 4.88611e+13 -5.12814e-06 4.88617e+13 -5.30672e-06 4.88624e+13 -5.49152e-06 4.88631e+13 -5.68275e-06 4.88637e+13 -5.88064e-06 4.88644e+13 -6.08542e-06 4.88651e+13 -6.29733e-06 4.88657e+13 -6.51663e-06 4.88664e+13 -6.74356e-06 4.88671e+13 -6.97839e-06 4.88677e+13 -7.22140e-06 4.88684e+13 -7.47287e-06 4.88691e+13 -7.73309e-06 4.88698e+13 -8.00238e-06 4.88704e+13 -8.28105e-06 4.88711e+13 -8.56942e-06 4.88718e+13 -8.86784e-06 4.88724e+13 -9.17664e-06 4.88731e+13 -9.49620e-06 4.88738e+13 -9.82689e-06 4.88744e+13 -1.01691e-05 4.88751e+13 -1.05232e-05 4.88758e+13 -1.08897e-05 4.88764e+13 -1.12689e-05 4.88771e+13 -1.16613e-05 4.88778e+13 -1.20674e-05 4.88785e+13 -1.24876e-05 4.88791e+13 -1.29224e-05 4.88798e+13 -1.33724e-05 4.88805e+13 -1.38381e-05 4.88811e+13 -1.43200e-05 4.88818e+13 -1.48187e-05 4.88825e+13 -1.53347e-05 4.88831e+13 -1.58687e-05 4.88838e+13 -1.64213e-05 4.88845e+13 -1.69931e-05 4.88851e+13 -1.75849e-05 4.88858e+13 -1.81972e-05 4.88865e+13 -1.88309e-05 4.88872e+13 -1.94867e-05 4.88878e+13 -2.01653e-05 4.88885e+13 -2.08675e-05 4.88892e+13 -2.15941e-05 4.88898e+13 -2.23461e-05 4.88905e+13 -2.31243e-05 4.88912e+13 -2.39295e-05 4.88918e+13 -2.47628e-05 4.88925e+13 -2.56252e-05 4.88932e+13 -2.65175e-05 4.88938e+13 -2.74409e-05 4.88945e+13 -2.83965e-05 4.88952e+13 -2.93853e-05 4.88959e+13 -3.04086e-05 4.88965e+13 -3.14676e-05 4.88972e+13 -3.25633e-05 4.88979e+13 -3.36973e-05 4.88985e+13 -3.48707e-05 4.88992e+13 -3.60851e-05 4.88999e+13 -3.73416e-05 4.89005e+13 -3.86420e-05 4.89012e+13 -3.99876e-05 4.89019e+13 -4.13801e-05 4.89026e+13 -4.28211e-05 4.89032e+13 -4.43123e-05 4.89039e+13 -4.58553e-05 4.89046e+13 -4.74522e-05 4.89052e+13 -4.91046e-05 4.89059e+13 -5.08146e-05 4.89066e+13 -5.25841e-05 4.89072e+13 -5.44152e-05 4.89079e+13 -5.63101e-05 4.89086e+13 -5.82710e-05 4.89092e+13 -6.03002e-05 4.89099e+13 -6.24000e-05 4.89106e+13 -6.45730e-05 4.89113e+13 -6.68216e-05 4.89119e+13 -6.91485e-05 4.89126e+13 -7.15565e-05 4.89133e+13 -7.40483e-05 4.89139e+13 -7.66269e-05 4.89146e+13 -7.92953e-05 4.89153e+13 -8.20566e-05 4.89159e+13 -8.49141e-05 4.89166e+13 -8.78710e-05 4.89173e+13 -9.09310e-05 4.89180e+13 -9.40974e-05 4.89186e+13 -9.73742e-05 4.89193e+13 -1.00765e-04 4.89200e+13 -1.04274e-04 4.89206e+13 -1.07905e-04 4.89213e+13 -1.11663e-04 4.89220e+13 -1.15551e-04 4.89226e+13 -1.19575e-04 4.89233e+13 -1.23739e-04 4.89240e+13 -1.28048e-04 4.89247e+13 -1.32507e-04 4.89253e+13 -1.37121e-04 4.89260e+13 -1.41896e-04 4.89267e+13 -1.46838e-04 4.89273e+13 -1.51951e-04 4.89280e+13 -1.57242e-04 4.89287e+13 -1.62718e-04 4.89293e+13 -1.68384e-04 4.89300e+13 -1.74248e-04 4.89307e+13 -1.80316e-04 4.89314e+13 -1.86595e-04 4.89320e+13 -1.93093e-04 4.89327e+13 -1.99817e-04 4.89334e+13 -2.06775e-04 4.89340e+13 -2.13975e-04 4.89347e+13 -2.21427e-04 4.89354e+13 -2.29137e-04 4.89360e+13 -2.37117e-04 4.89367e+13 -2.45374e-04 4.89374e+13 -2.53919e-04 4.89381e+13 -2.62761e-04 4.89387e+13 -2.71911e-04 4.89394e+13 -2.81380e-04 4.89401e+13 -2.91178e-04 4.89407e+13 -3.01318e-04 4.89414e+13 -3.11811e-04 4.89421e+13 -3.22669e-04 4.89427e+13 -3.33905e-04 4.89434e+13 -3.45533e-04 4.89441e+13 -3.57565e-04 4.89448e+13 -3.70017e-04 4.89454e+13 -3.82902e-04 4.89461e+13 -3.96236e-04 4.89468e+13 -4.10034e-04 4.89474e+13 -4.24312e-04 4.89481e+13 -4.39088e-04 4.89488e+13 -4.54379e-04 4.89494e+13 -4.70202e-04 4.89501e+13 -4.86575e-04 4.89508e+13 -5.03519e-04 4.89515e+13 -5.21053e-04 4.89521e+13 -5.39198e-04 4.89528e+13 -5.57975e-04 4.89535e+13 -5.77405e-04 4.89541e+13 -5.97512e-04 4.89548e+13 -6.18319e-04 4.89555e+13 -6.39851e-04 4.89561e+13 -6.62133e-04 4.89568e+13 -6.85190e-04 4.89575e+13 -7.09050e-04 4.89582e+13 -7.33742e-04 4.89588e+13 -7.59293e-04 4.89595e+13 -7.85734e-04 4.89602e+13 -8.13095e-04 4.89608e+13 -8.41410e-04 4.89615e+13 -8.70710e-04 4.89622e+13 -9.01031e-04 4.89629e+13 -9.32408e-04 4.89635e+13 -9.64877e-04 4.89642e+13 -9.98477e-04 4.89649e+13 -1.03325e-03 4.89655e+13 -1.06923e-03 4.89662e+13 -1.10646e-03 4.89669e+13 -1.14499e-03 4.89675e+13 -1.18486e-03 4.89682e+13 -1.22612e-03 4.89689e+13 -1.26882e-03 4.89696e+13 -1.31301e-03 4.89702e+13 -1.35873e-03 4.89709e+13 -1.40604e-03 4.89716e+13 -1.45501e-03 4.89722e+13 -1.50567e-03 4.89729e+13 -1.55811e-03 4.89736e+13 -1.61236e-03 4.89743e+13 -1.66851e-03 4.89749e+13 -1.72661e-03 4.89756e+13 -1.78674e-03 4.89763e+13 -1.84896e-03 4.89769e+13 -1.91335e-03 4.89776e+13 -1.97998e-03 4.89783e+13 -2.04892e-03 4.89789e+13 -2.12027e-03 4.89796e+13 -2.19411e-03 4.89803e+13 -2.27051e-03 4.89810e+13 -2.34958e-03 4.89816e+13 -2.43140e-03 4.89823e+13 -2.51607e-03 4.89830e+13 -2.60369e-03 4.89836e+13 -2.69435e-03 4.89843e+13 -2.78818e-03 4.89850e+13 -2.88527e-03 4.89857e+13 -2.98575e-03 4.89863e+13 -3.08972e-03 4.89870e+13 -3.19731e-03 4.89877e+13 -3.30865e-03 4.89883e+13 -3.42387e-03 4.89890e+13 -3.54310e-03 4.89897e+13 -3.66648e-03 4.89903e+13 -3.79416e-03 4.89910e+13 -3.92628e-03 4.89917e+13 -4.06301e-03 4.89924e+13 -4.20449e-03 4.89930e+13 -4.35091e-03 4.89937e+13 -4.50242e-03 4.89944e+13 -4.65921e-03 4.89950e+13 -4.82146e-03 4.89957e+13 -4.98935e-03 4.89964e+13 -5.16310e-03 4.89971e+13 -5.34289e-03 4.89977e+13 -5.52895e-03 4.89984e+13 -5.72148e-03 4.89991e+13 -5.92072e-03 4.89997e+13 -6.12690e-03 nan -6.34026e-03 nan -6.56104e-03 nan -6.78952e-03 nan -7.02595e-03 nan -7.27062e-03 nan -7.52380e-03 nan -7.78580e-03 nan -8.05693e-03 nan -8.33749e-03 nan -8.62783e-03 nan -8.92828e-03 nan -9.23919e-03 nan -9.56093e-03 nan -9.89387e-03 nan -1.02384e-02 nan -1.05949e-02 nan -1.09639e-02 nan -1.13457e-02 nan -1.17408e-02 nan -1.21496e-02 nan -1.25727e-02 nan -1.30105e-02 nan -1.34636e-02 nan -1.39324e-02 nan -1.44176e-02 nan -1.49197e-02 nan -1.54392e-02 nan -1.59769e-02 nan -1.65332e-02 nan -1.71090e-02 nan -1.77047e-02 nan -1.83213e-02 nan -1.89593e-02 nan -1.96195e-02 nan -2.03027e-02 nan -2.10097e-02 nan -2.17413e-02 nan -2.24984e-02 nan -2.32819e-02 nan -2.40926e-02 nan -2.49316e-02 nan -2.57998e-02 nan -2.66982e-02 nan -2.76280e-02 nan -2.85900e-02 nan -2.95856e-02 nan -3.06159e-02 nan -3.16820e-02 nan -3.27853e-02 nan -3.39270e-02 nan -3.51084e-02 nan -3.63310e-02 nan -3.75962e-02 nan -3.89054e-02 nan -4.02602e-02 nan -4.16622e-02 nan -4.31130e-02 nan -4.46143e-02 nan -4.61679e-02 nan -4.77756e-02 nan -4.94393e-02 nan -5.11609e-02 nan -5.29425e-02 nan -5.47861e-02 nan -5.66939e-02 nan -5.86682e-02 nan -6.07112e-02 nan -6.28253e-02 nan -6.50131e-02 nan -6.72771e-02 nan -6.96199e-02 nan -7.20442e-02 nan -7.45530e-02 nan -7.71492e-02 nan -7.98358e-02 nan -8.26159e-02 nan -8.54928e-02 nan -8.84699e-02 nan -9.15507e-02 nan -9.47388e-02 nan -9.80379e-02 nan -1.01452e-01 nan -1.04985e-01 nan -1.08641e-01 nan -1.12424e-01 nan -1.16339e-01 nan -1.20390e-01 nan -1.24582e-01 nan -1.28921e-01 nan -1.33410e-01 nan -1.38056e-01 nan -1.42863e-01 nan -1.47838e-01 nan -1.52987e-01 nan -1.58314e-01 nan -1.63827e-01 nan -1.69532e-01 nan -1.75436e-01 nan -1.81545e-01 nan -1.87867e-01 nan -1.94409e-01 nan -2.01179e-01 nan -2.08184e-01 nan -2.15434e-01 nan -2.22936e-01 nan -2.30699e-01 nan -2.38733e-01 nan -2.47046e-01 nan -2.55649e-01 nan -2.64552e-01 nan -2.73764e-01 nan -2.83298e-01 nan -2.93163e-01 nan -3.03372e-01 nan -3.13936e-01 nan -3.24868e-01 nan -3.36181e-01 nan -3.47888e-01 nan -3.60002e-01 nan -3.72539e-01 nan -3.85512e-01 nan -3.98936e-01 nan -4.12829e-01 2.18058e+13 -4.27205e-01 2.15972e+13 -4.42081e-01 2.13906e+13 -4.57476e-01 2.11859e+13 -4.73406e-01 2.09833e+13 -4.89892e-01 2.07825e+13 -5.06951e-01 2.05837e+13 -5.24605e-01 2.03868e+13 -5.42873e-01 2.01917e+13 -5.61778e-01 1.99985e+13 -5.81341e-01 1.98072e+13 -6.01585e-01 1.96177e+13 -6.22534e-01 1.94300e+13 -6.44212e-01 1.92442e+13 -6.66646e-01 1.90600e+13 -6.89860e-01 1.88777e+13 -7.13883e-01 1.86971e+13 -7.38743e-01 1.85182e+13 -7.64468e-01 1.83411e+13 -7.91089e-01 1.81656e+13 -8.18637e-01 1.79918e+13 -8.47145e-01 1.78197e+13 -8.76645e-01 1.76492e+13 -9.07172e-01 1.74803e+13 -9.38763e-01 1.73131e+13 -9.71453e-01 1.71475e+13 -1.00528e+00 1.69834e+13 -1.04029e+00 1.68209e+13 -1.07652e+00 1.66600e+13 -1.11400e+00 1.65006e+13 -1.15280e+00 1.63428e+13 -1.19294e+00 1.61864e+13 -1.23448e+00 1.60316e+13 -1.27747e+00 1.58782e+13 -1.32196e+00 1.57263e+13 -1.36799e+00 1.55758e+13 -1.41563e+00 1.54268e+13 -1.46492e+00 1.52792e+13 -1.51594e+00 1.51330e+13 -1.56873e+00 1.49883e+13 -1.62335e+00 1.48449e+13 -1.67988e+00 1.47028e+13 -1.73838e+00 1.45622e+13 -1.79892e+00 1.44229e+13 -1.86156e+00 1.42849e+13 -1.92639e+00 1.41482e+13 -1.99347e+00 1.40129e+13 -2.06289e+00 nan -2.13473e+00 nan -2.20906e+00 nan -2.28599e+00 nan -2.36559e+00 nan -2.44797e+00 nan -2.53322e+00 nan -2.62143e+00 nan -2.71272e+00 nan -2.80718e+00 nan -2.90494e+00 nan -3.00610e+00 nan -3.11078e+00 nan -3.21911e+00 nan -3.33120e+00 nan -3.44721e+00 nan -3.56725e+00 nan -3.69147e+00 nan -3.82002e+00 nan -3.95304e+00 nan -4.09070e+00 nan -4.23315e+00 nan -4.38056e+00 nan -4.53311e+00 nan -4.69096e+00 nan -4.85432e+00 nan -5.02336e+00 nan -5.19829e+00 nan -5.37931e+00 nan -5.56663e+00 nan -5.76048e+00 nan -5.96108e+00 nan -6.16866e+00 nan -6.38347e+00 nan -6.60576e+00 nan -6.83580e+00 nan -7.07384e+00 nan -7.32017e+00 nan -7.57508e+00 nan -7.83887e+00 nan -8.11184e+00 nan -8.39432e+00 nan -8.68664e+00 nan -8.98913e+00 nan -9.30216e+00 nan -9.62609e+00 nan -9.96130e+00 nan -1.03082e+01 nan -1.06671e+01 nan -1.10386e+01 nan -1.14230e+01 nan -1.18208e+01 nan -1.22324e+01 nan -1.26584e+01 nan -1.30992e+01 nan -1.35554e+01 nan -1.40274e+01 nan -1.45159e+01 nan -1.50214e+01 nan -1.55444e+01 nan -1.60858e+01 nan -1.66459e+01 nan -1.72256e+01 nan -1.78254e+01 nan -1.84461e+01 nan -1.90885e+01 nan -1.97532e+01 nan -2.04411e+01 nan -2.11529e+01 nan -2.18895e+01 nan -2.26518e+01 nan -2.34406e+01 nan -2.42569e+01 nan -2.51015e+01 nan -2.59757e+01 nan -2.68802e+01 nan -2.78163e+01 nan -2.87849e+01 nan -2.97873e+01 nan -3.08246e+01 nan -3.18980e+01 nan -3.30088e+01 nan -3.41582e+01 nan -3.53477e+01 nan -3.65786e+01 nan -3.78524e+01 nan -3.91706e+01 nan -4.05346e+01 nan -4.19461e+01 nan -4.34068e+01 nan -4.49184e+01 nan -4.64826e+01 nan -4.81012e+01 nan -4.97763e+01 nan -5.15096e+01 nan -5.33033e+01 nan -5.51595e+01 nan -5.70804e+01 nan -5.90681e+01 nan -6.11250e+01 nan -6.32536e+01 nan -6.54562e+01 nan -6.77356e+01 nan -7.00944e+01 nan -7.25353e+01 nan -7.50612e+01 nan -7.76750e+01 nan -8.03799e+01 nan -8.31790e+01 nan -8.60755e+01 nan -8.90729e+01 nan -9.21747e+01 nan -9.53845e+01 nan -9.87061e+01 nan -1.02143e+02 nan -1.05700e+02 nan -1.09381e+02 nan -1.13190e+02 nan -1.17132e+02 nan -1.21211e+02 nan -1.25432e+02 nan -1.29799e+02 nan -1.34319e+02 nan -1.38997e+02 nan -1.43837e+02 nan -1.48846e+02 nan -1.54029e+02 nan -1.59393e+02 nan -1.64944e+02 nan -1.70687e+02 nan -1.76631e+02 nan -1.82782e+02 nan -1.89147e+02 nan -1.95734e+02 nan -2.02550e+02 nan -2.09603e+02 nan -2.16902e+02 nan -2.24455e+02 nan -2.32272e+02 nan -2.40360e+02 -nan -2.48730e+02 -nan -2.57392e+02 -nan -2.66355e+02 -nan -2.75630e+02 -nan -2.85228e+02 -nan -2.95161e+02 -nan -3.05439e+02 -nan -3.16076e+02 -nan -3.27082e+02 -nan -3.38472e+02 -nan -3.50259e+02 -nan -3.62456e+02 -nan -3.75078e+02 -nan -3.88139e+02 -nan -4.01656e+02 -nan -4.15642e+02 -nan -4.30116e+02 -nan -4.45094e+02 -nan -4.60594e+02 -nan -4.76633e+02 -nan -4.93231e+02 -nan -5.10407e+02 -nan -5.28181e+02 -nan -5.46574e+02 -nan -5.65607e+02 -nan -5.85303e+02 -nan -6.05685e+02 -nan -6.26777e+02 -nan -6.48603e+02 -nan -6.71189e+02 -nan -6.94562e+02 -nan -7.18749e+02 -nan -7.43778e+02 -nan -7.69679e+02 -nan -7.96481e+02 -nan -8.24217e+02 -nan -8.52919e+02 -nan -8.82620e+02 -nan -9.13356e+02 -nan -9.45161e+02 -nan -9.78075e+02 -nan -1.01213e+03 -nan -1.04738e+03 -nan -1.08385e+03 -nan -1.12160e+03 -nan -1.16065e+03 -nan -1.20107e+03 -nan -1.24290e+03 -nan -1.28618e+03 -nan -1.33097e+03 -nan -1.37731e+03 -nan -1.42528e+03 -nan -1.47491e+03 -nan -1.52627e+03 -nan -1.57942e+03 -nan -1.63442e+03 -nan -1.69133e+03 -nan -1.75023e+03 -nan -1.81118e+03 -nan -1.87425e+03 -nan -1.93952e+03 -nan -2.00706e+03 -nan -2.07695e+03 -nan -2.14928e+03 -nan -2.22412e+03 -nan -2.30157e+03 -nan -2.38172e+03 -nan -2.46466e+03 -nan -2.55048e+03 -nan -2.63930e+03 -nan -2.73121e+03 -nan -2.82632e+03 -nan -2.92474e+03 -nan -3.02659e+03 -nan -3.13198e+03 -nan -3.24105e+03 -nan -3.35391e+03 -nan -3.47070e+03 -nan -3.59156e+03 -nan -3.71663e+03 -nan -3.84606e+03 -nan -3.97999e+03 -nan -4.11858e+03 -nan -4.26200e+03 -nan -4.41042e+03 -nan -4.56401e+03 -nan -4.72294e+03 -nan -4.88740e+03 -nan -5.05760e+03 -nan -5.23372e+03 -nan -5.41597e+03 -nan -5.60457e+03 -nan -5.79974e+03 -nan -6.00171e+03 -nan -6.21071e+03 -nan -6.42698e+03 -nan -6.65079e+03 -nan -6.88239e+03 -nan -7.12205e+03 -nan -7.37007e+03 -nan -7.62671e+03 -nan -7.89230e+03 -nan -8.16713e+03 -nan -8.45154e+03 -nan -8.74585e+03 -nan -9.05040e+03 -nan -9.36557e+03 -nan -9.69170e+03 -nan -1.00292e+04 -nan -1.03784e+04 -nan -1.07399e+04 -nan -1.11138e+04 -nan -1.15009e+04 -nan -1.19014e+04 -nan -1.23158e+04 -nan -1.27447e+04 -nan -1.31885e+04 -nan -1.36477e+04 -nan -1.41230e+04 -nan -1.46148e+04 -nan -1.51237e+04 -nan -1.56504e+04 -nan -1.61954e+04 -nan -1.67594e+04 -nan -1.73430e+04 -nan -1.79469e+04 -nan -1.85719e+04 -nan -1.92186e+04 -nan -1.98879e+04 -nan -2.05804e+04 -nan -2.12971e+04 -nan -2.20387e+04 -nan -2.28062e+04 -nan -2.36003e+04 -nan -2.44222e+04 -nan -2.52726e+04 -nan -2.61527e+04 -nan -2.70634e+04 -nan -2.80059e+04 -nan -2.89811e+04 -nan -2.99903e+04 -nan -3.10347e+04 -nan -3.21154e+04 -nan -3.32337e+04 -nan -3.43910e+04 -nan -3.55886e+04 -nan -3.68280e+04 -nan -3.81104e+04 -nan -3.94375e+04 -nan -4.08109e+04 -nan -4.22320e+04 -nan -4.37027e+04 -nan -4.52245e+04 -nan -4.67994e+04 -nan -4.84291e+04 -nan -5.01155e+04 -nan -5.18607e+04 -nan -5.36667e+04 -nan -5.55355e+04 -nan -5.74694e+04 -nan -5.94707e+04 -nan -6.15416e+04 -nan -6.36847e+04 -nan -6.59024e+04 -nan -6.81973e+04 -nan -7.05721e+04 -nan -7.30297e+04 -nan -7.55728e+04 -nan -7.82045e+04 -nan -8.09278e+04 -nan -8.37459e+04 -nan -8.66622e+04 -nan -8.96801e+04 -nan -9.28030e+04 -nan -9.60347e+04 -nan -9.93789e+04 -nan -1.02840e+05 -nan -1.06421e+05 -nan -1.10127e+05 -nan -1.13962e+05 -nan -1.17930e+05 -nan -1.22037e+05 -nan -1.26286e+05 -nan -1.30684e+05 -nan -1.35235e+05 -nan -1.39944e+05 -nan -1.44818e+05 -nan -1.49861e+05 -nan -1.55079e+05 -nan -1.60479e+05 -nan -1.66068e+05 -nan -1.71851e+05 -nan -1.77835e+05 -nan -1.84028e+05 -nan -1.90436e+05 -nan -1.97068e+05 -nan -2.03930e+05 -nan -2.11032e+05 -nan -2.18381e+05 -nan -2.25985e+05 -nan -2.33855e+05 -nan -2.41998e+05 -nan -2.50426e+05 -nan -2.59146e+05 -nan -2.68170e+05 -nan -2.77509e+05 -nan -2.87173e+05 -nan -2.97173e+05 -nan -3.07521e+05 -nan -3.18230e+05 -nan -3.29312e+05 -nan -3.40779e+05 -nan -3.52646e+05 -nan -3.64927e+05 -nan -3.77635e+05 -nan -3.90785e+05 -nan -4.04393e+05 -nan -4.18475e+05 -nan -4.33048e+05 -nan -4.48128e+05 -nan -4.63733e+05 -nan -4.79882e+05 -nan -4.96593e+05 -nan -5.13886e+05 -nan -5.31781e+05 -nan -5.50299e+05 -nan -5.69462e+05 -nan -5.89292e+05 -nan -6.09813e+05 -nan -6.31049e+05 -nan -6.53024e+05 -nan -6.75764e+05 -nan -6.99296e+05 -nan -7.23648e+05 -nan -7.48848e+05 -nan -7.74925e+05 -nan -8.01910e+05 -nan -8.29835e+05 -nan -8.58732e+05 -nan -8.88636e+05 -nan -9.19581e+05 -nan -9.51604e+05 -nan -9.84741e+05 -nan -1.01903e+06 -nan -1.05452e+06 -nan -1.09124e+06 -nan -1.12924e+06 -nan -1.16856e+06 -nan -1.20926e+06 -nan -1.25137e+06 -nan -1.29494e+06 -nan -1.34004e+06 -nan -1.38670e+06 -nan -1.43499e+06 -nan -1.48496e+06 -nan -1.53667e+06 -nan -1.59018e+06 -nan -1.64556e+06 -nan -1.70286e+06 -nan -1.76216e+06 -nan -1.82353e+06 -nan -1.88703e+06 -nan -1.95274e+06 -nan -2.02074e+06 -nan -2.09111e+06 -nan -2.16393e+06 -nan -2.23928e+06 -nan -2.31726e+06 -nan -2.39795e+06 -nan -2.48146e+06 -nan -2.56787e+06 -nan -2.65729e+06 -nan -2.74982e+06 -nan -2.84558e+06 -nan -2.94467e+06 -nan -3.04722e+06 -nan -3.15333e+06 -nan -3.26314e+06 -nan -3.37677e+06 -nan -3.49436e+06 -nan -3.61604e+06 -nan -3.74196e+06 -nan -3.87227e+06 -nan -4.00712e+06 -nan -4.14666e+06 -nan -4.29105e+06 -nan -4.44048e+06 -nan -4.59511e+06 -nan -4.75513e+06 -nan -4.92072e+06 -nan -5.09207e+06 -nan -5.26939e+06 -nan -5.45289e+06 -nan -5.64278e+06 -nan -5.83927e+06 -nan -6.04262e+06 -nan -6.25304e+06 -nan -6.47079e+06 -nan -6.69612e+06 -nan -6.92930e+06 -nan -7.17060e+06 -nan -7.42030e+06 -nan -7.67870e+06 -nan -7.94609e+06 -nan -8.22280e+06 -nan -8.50914e+06 -nan -8.80546e+06 -nan -9.11209e+06 -nan -9.42940e+06 -nan -9.75776e+06 -nan -1.00976e+07 -nan -1.04492e+07 -nan -1.08131e+07 -nan -1.11896e+07 -nan -1.15793e+07 -nan -1.19825e+07 -nan -1.23997e+07 -nan -1.28315e+07 -nan -1.32784e+07 -nan -1.37408e+07 -nan -1.42193e+07 -nan -1.47144e+07 -nan -1.52268e+07 -nan -1.57571e+07 -nan -1.63058e+07 -nan -1.68736e+07 -nan -1.74612e+07 -nan -1.80692e+07 -nan -1.86985e+07 -nan -1.93496e+07 -nan -2.00234e+07 -nan -2.07207e+07 -nan -2.14422e+07 -nan -2.21889e+07 -nan -2.29616e+07 -nan -2.37612e+07 -nan -2.45886e+07 -nan -2.54449e+07 -nan -2.63310e+07 -nan -2.72479e+07 -nan -2.81967e+07 -nan -2.91786e+07 -nan -3.01947e+07 -nan -3.12462e+07 -nan -3.23343e+07 -nan -3.34603e+07 -nan -3.46255e+07 -nan -3.58312e+07 -nan -3.70790e+07 -nan -3.83702e+07 -nan -3.97063e+07 -nan -4.10890e+07 -nan -4.25199e+07 -nan -4.40006e+07 -nan -4.55328e+07 -nan -4.71184e+07 -nan -4.87592e+07 -nan -5.04571e+07 -nan -5.22142e+07 -nan -5.40325e+07 -nan -5.59140e+07 -nan -5.78611e+07 -nan -5.98760e+07 -nan -6.19611e+07 -nan -6.41188e+07 -nan -6.63516e+07 -nan -6.86621e+07 -nan -7.10532e+07 -nan -7.35274e+07 -nan -7.60879e+07 -nan -7.87375e+07 -nan -8.14794e+07 -nan -8.43167e+07 -nan -8.72529e+07 -nan -9.02913e+07 -nan -9.34355e+07 -nan -9.66892e+07 -nan +1.36193e-07 1.26309e+00 +1.40936e-07 1.39969e+00 +1.45843e-07 1.55107e+00 +1.50922e-07 1.71881e+00 +1.56178e-07 1.90470e+00 +1.61616e-07 2.11070e+00 +1.67244e-07 2.33897e+00 +1.73068e-07 2.59192e+00 +1.79095e-07 2.87224e+00 +1.85332e-07 3.18287e+00 +1.91785e-07 3.52709e+00 +1.98464e-07 3.90855e+00 +2.05375e-07 4.33125e+00 +2.12527e-07 4.79967e+00 +2.19928e-07 5.31875e+00 +2.27586e-07 5.89397e+00 +2.35512e-07 6.53140e+00 +2.43713e-07 7.23777e+00 +2.52200e-07 8.02053e+00 +2.60982e-07 8.88794e+00 +2.70070e-07 9.84917e+00 +2.79475e-07 1.09143e+01 +2.89207e-07 1.20947e+01 +2.99278e-07 1.34028e+01 +3.09700e-07 1.48523e+01 +3.20484e-07 1.64585e+01 +3.31645e-07 1.82385e+01 +3.43194e-07 2.02110e+01 +3.55145e-07 2.23968e+01 +3.67512e-07 2.48190e+01 +3.80310e-07 2.75031e+01 +3.93553e-07 3.04776e+01 +4.07258e-07 3.37737e+01 +4.21440e-07 3.74263e+01 +4.36116e-07 4.14739e+01 +4.51303e-07 4.59593e+01 +4.67018e-07 5.09298e+01 +4.83281e-07 5.64378e+01 +5.00111e-07 6.25415e+01 +5.17526e-07 6.93053e+01 +5.35548e-07 7.68007e+01 +5.54197e-07 8.51066e+01 +5.73496e-07 9.43108e+01 +5.93467e-07 1.04510e+02 +6.14133e-07 1.15813e+02 +6.35519e-07 1.28338e+02 +6.57650e-07 1.42218e+02 +6.80551e-07 1.57599e+02 +7.04250e-07 1.74643e+02 +7.28774e-07 1.93530e+02 +7.54153e-07 2.14461e+02 +7.80414e-07 2.37654e+02 +8.07591e-07 2.63357e+02 +8.35714e-07 2.91838e+02 +8.64816e-07 3.23401e+02 +8.94931e-07 3.58376e+02 +9.26095e-07 3.97134e+02 +9.58345e-07 4.40084e+02 +9.91717e-07 4.87679e+02 +1.02625e-06 5.40421e+02 +1.06199e-06 5.98867e+02 +1.09897e-06 6.63634e+02 +1.13724e-06 7.35405e+02 +1.17684e-06 8.14939e+02 +1.21782e-06 9.03074e+02 +1.26023e-06 1.00074e+03 +1.30412e-06 1.10897e+03 +1.34953e-06 1.22890e+03 +1.39653e-06 1.36181e+03 +1.44516e-06 1.50909e+03 +1.49548e-06 1.67230e+03 +1.54756e-06 1.85315e+03 +1.60145e-06 2.05357e+03 +1.65722e-06 2.27566e+03 +1.71493e-06 2.52177e+03 +1.77464e-06 2.79450e+03 +1.83644e-06 3.09672e+03 +1.90039e-06 3.43163e+03 +1.96657e-06 3.80276e+03 +2.03505e-06 4.21403e+03 +2.10592e-06 4.66977e+03 +2.17925e-06 5.17480e+03 +2.25514e-06 5.73446e+03 +2.33367e-06 6.35463e+03 +2.41494e-06 7.04188e+03 +2.49904e-06 7.80346e+03 +2.58606e-06 8.64739e+03 +2.67611e-06 9.58260e+03 +2.76930e-06 1.06190e+04 +2.86574e-06 1.17674e+04 +2.96553e-06 1.30400e+04 +3.06880e-06 1.44503e+04 +3.17567e-06 1.60131e+04 +3.28625e-06 1.77449e+04 +3.40069e-06 1.96640e+04 +3.51911e-06 2.17906e+04 +3.64166e-06 2.41473e+04 +3.76847e-06 2.67588e+04 +3.89970e-06 2.96527e+04 +4.03550e-06 3.28596e+04 +4.17603e-06 3.64134e+04 +4.32145e-06 4.03515e+04 +4.47194e-06 4.47154e+04 +4.62767e-06 4.95514e+04 +4.78881e-06 5.49103e+04 +4.95558e-06 6.08489e+04 +5.12814e-06 6.74296e+04 +5.30672e-06 7.47221e+04 +5.49152e-06 8.28032e+04 +5.68275e-06 9.17583e+04 +5.88064e-06 1.01682e+05 +6.08542e-06 1.12679e+05 +6.29733e-06 1.24865e+05 +6.51663e-06 1.38369e+05 +6.74356e-06 1.53333e+05 +6.97839e-06 1.69916e+05 +7.22140e-06 1.88293e+05 +7.47287e-06 2.08656e+05 +7.73309e-06 2.31222e+05 +8.00238e-06 2.56229e+05 +8.28105e-06 2.83940e+05 +8.56942e-06 3.14648e+05 +8.86784e-06 3.48677e+05 +9.17664e-06 3.86386e+05 +9.49620e-06 4.28173e+05 +9.82689e-06 4.74480e+05 +1.01691e-05 5.25795e+05 +1.05232e-05 5.82659e+05 +1.08897e-05 6.45673e+05 +1.12689e-05 7.15502e+05 +1.16613e-05 7.92883e+05 +1.20674e-05 8.78633e+05 +1.24876e-05 9.73656e+05 +1.29224e-05 1.07896e+06 +1.33724e-05 1.19564e+06 +1.38381e-05 1.32495e+06 +1.43200e-05 1.46825e+06 +1.48187e-05 1.62704e+06 +1.53347e-05 1.80300e+06 +1.58687e-05 1.99799e+06 +1.64213e-05 2.21407e+06 +1.69931e-05 2.45352e+06 +1.75849e-05 2.71887e+06 +1.81972e-05 3.01291e+06 +1.88309e-05 3.33876e+06 +1.94867e-05 3.69984e+06 +2.01653e-05 4.09998e+06 +2.08675e-05 4.54339e+06 +2.15941e-05 5.03475e+06 +2.23461e-05 5.57926e+06 +2.31243e-05 6.18265e+06 +2.39295e-05 6.85130e+06 +2.47628e-05 7.59226e+06 +2.56252e-05 8.41336e+06 +2.65175e-05 9.32326e+06 +2.74409e-05 1.03316e+07 +2.83965e-05 1.14489e+07 +2.93853e-05 1.26871e+07 +3.04086e-05 1.40592e+07 +3.14676e-05 1.55797e+07 +3.25633e-05 1.72646e+07 +3.36973e-05 1.91318e+07 +3.48707e-05 2.12009e+07 +3.60851e-05 2.34937e+07 +3.73416e-05 2.60346e+07 +3.86420e-05 2.88502e+07 +3.99876e-05 3.19703e+07 +4.13801e-05 3.54279e+07 +4.28211e-05 3.92594e+07 +4.43123e-05 4.35052e+07 +4.58553e-05 4.82103e+07 +4.74522e-05 5.34242e+07 +4.91046e-05 5.92020e+07 +5.08146e-05 6.56047e+07 +5.25841e-05 7.26998e+07 +5.44152e-05 8.05622e+07 +5.63101e-05 8.92749e+07 +5.82710e-05 9.89299e+07 +6.03002e-05 1.09629e+08 +6.24000e-05 1.21485e+08 +6.45730e-05 1.34624e+08 +6.68216e-05 1.49184e+08 +6.91485e-05 1.65318e+08 +7.15565e-05 1.83197e+08 +7.40483e-05 2.03009e+08 +7.66269e-05 2.24964e+08 +7.92953e-05 2.49294e+08 +8.20566e-05 2.76255e+08 +8.49141e-05 3.06132e+08 +8.78710e-05 3.39240e+08 +9.09310e-05 3.75928e+08 +9.40974e-05 4.16585e+08 +9.73742e-05 4.61638e+08 +1.00765e-04 5.11564e+08 +1.04274e-04 5.66889e+08 +1.07905e-04 6.28198e+08 +1.11663e-04 6.96137e+08 +1.15551e-04 7.71424e+08 +1.19575e-04 8.54853e+08 +1.23739e-04 9.47305e+08 +1.28048e-04 1.04975e+09 +1.32507e-04 1.16329e+09 +1.37121e-04 1.28909e+09 +1.41896e-04 1.42851e+09 +1.46838e-04 1.58300e+09 +1.51951e-04 1.75420e+09 +1.57242e-04 1.94392e+09 +1.62718e-04 2.15415e+09 +1.68384e-04 2.38712e+09 +1.74248e-04 2.64528e+09 +1.80316e-04 2.93137e+09 +1.86595e-04 3.24840e+09 +1.93093e-04 3.59971e+09 +1.99817e-04 3.98901e+09 +2.06775e-04 4.42042e+09 +2.13975e-04 4.89849e+09 +2.21427e-04 5.42825e+09 +2.29137e-04 6.01532e+09 +2.37117e-04 6.66587e+09 +2.45374e-04 7.38678e+09 +2.53919e-04 8.18565e+09 +2.62761e-04 9.07092e+09 +2.71911e-04 1.00519e+10 +2.81380e-04 1.11390e+10 +2.91178e-04 1.23437e+10 +3.01318e-04 1.36787e+10 +3.11811e-04 1.51580e+10 +3.22669e-04 1.67974e+10 +3.33905e-04 1.86140e+10 +3.45533e-04 2.06271e+10 +3.57565e-04 2.28579e+10 +3.70017e-04 2.53299e+10 +3.82902e-04 2.80694e+10 +3.96236e-04 3.11050e+10 +4.10034e-04 3.63140e+10 +4.24312e-04 4.77862e+10 +4.39088e-04 6.61641e+10 +4.54379e-04 9.20409e+10 +4.70202e-04 1.26044e+11 +4.86575e-04 1.68835e+11 +5.03519e-04 2.21112e+11 +5.21053e-04 2.83608e+11 +5.39198e-04 3.57094e+11 +5.57975e-04 4.42377e+11 +5.77405e-04 5.40298e+11 +5.97512e-04 6.51739e+11 +6.18319e-04 7.77612e+11 +6.39851e-04 9.18866e+11 +6.62133e-04 1.07648e+12 +6.85190e-04 1.25147e+12 +7.09050e-04 1.44486e+12 +7.33742e-04 1.65771e+12 +7.59293e-04 1.89111e+12 +7.85734e-04 2.14612e+12 +8.13095e-04 2.42385e+12 +8.41410e-04 2.72538e+12 +8.70710e-04 3.05177e+12 +9.01031e-04 3.40406e+12 +9.32408e-04 3.78326e+12 +9.64877e-04 4.19029e+12 +9.98477e-04 4.62602e+12 +1.03325e-03 5.09120e+12 +1.06923e-03 5.58646e+12 +1.10646e-03 6.11226e+12 +1.14499e-03 6.66889e+12 +1.18486e-03 7.25637e+12 +1.22612e-03 7.87448e+12 +1.26882e-03 8.52264e+12 +1.31301e-03 9.19989e+12 +1.35873e-03 9.90568e+12 +1.40604e-03 1.06406e+13 +1.45501e-03 1.14053e+13 +1.50567e-03 1.22004e+13 +1.55811e-03 1.30264e+13 +1.61236e-03 1.38837e+13 +1.66851e-03 1.47726e+13 +1.72661e-03 1.56932e+13 +1.78674e-03 1.66454e+13 +1.84896e-03 1.76290e+13 +1.91335e-03 1.86433e+13 +1.97998e-03 1.96875e+13 +2.04892e-03 2.07605e+13 +2.12027e-03 2.18617e+13 +2.19411e-03 2.29906e+13 +2.27051e-03 2.41465e+13 +2.34958e-03 2.53287e+13 +2.43140e-03 2.65359e+13 +2.51607e-03 2.77668e+13 +2.60369e-03 2.90197e+13 +2.69435e-03 3.02925e+13 +2.78818e-03 3.15825e+13 +2.88527e-03 3.28868e+13 +2.98575e-03 3.42017e+13 +3.08972e-03 3.55231e+13 +3.19731e-03 3.68461e+13 +3.30865e-03 3.81649e+13 +3.42387e-03 3.94731e+13 +3.54310e-03 4.07629e+13 +3.66648e-03 4.20258e+13 +3.79416e-03 4.32518e+13 +3.92628e-03 4.44295e+13 +4.06301e-03 4.55460e+13 +4.20449e-03 4.65866e+13 +4.35091e-03 4.75346e+13 +4.50242e-03 4.83712e+13 +4.65921e-03 4.90749e+13 +4.82146e-03 4.96216e+13 +4.98935e-03 4.99842e+13 +5.16310e-03 5.01394e+13 +5.34289e-03 5.00998e+13 +5.52895e-03 4.98958e+13 +5.72148e-03 4.95670e+13 +5.92072e-03 4.91637e+13 +6.12690e-03 4.88048e+13 +6.34026e-03 4.84873e+13 +6.56104e-03 4.81719e+13 +6.78952e-03 4.78585e+13 +7.02595e-03 4.75472e+13 +7.27062e-03 4.72378e+13 +7.52380e-03 4.69305e+13 +7.78580e-03 4.66252e+13 +8.05693e-03 4.63219e+13 +8.33749e-03 4.60205e+13 +8.62783e-03 4.57211e+13 +8.92828e-03 4.54237e+13 +9.23919e-03 4.51282e+13 +9.56093e-03 4.48346e+13 +9.89387e-03 4.45429e+13 +1.02384e-02 4.42531e+13 +1.05949e-02 4.39652e+13 +1.09639e-02 4.36792e+13 +1.13457e-02 4.33951e+13 +1.17408e-02 4.31128e+13 +1.21496e-02 4.28323e+13 +1.25727e-02 4.25536e+13 +1.30105e-02 4.22768e+13 +1.34636e-02 4.20018e+13 +1.39324e-02 4.17285e+13 +1.44176e-02 4.14570e+13 +1.49197e-02 4.11873e+13 +1.54392e-02 4.09194e+13 +1.59769e-02 4.06532e+13 +1.65332e-02 4.03887e+13 +1.71090e-02 4.01259e+13 +1.77047e-02 3.98649e+13 +1.83213e-02 3.96056e+13 +1.89593e-02 3.93479e+13 +1.96195e-02 3.90919e+13 +2.03027e-02 3.88376e+13 +2.10097e-02 3.85849e+13 +2.17413e-02 3.83339e+13 +2.24984e-02 3.80845e+13 +2.32819e-02 3.78368e+13 +2.40926e-02 3.75906e+13 +2.49316e-02 3.73461e+13 +2.57998e-02 3.71031e+13 +2.66982e-02 3.68617e+13 +2.76280e-02 3.66219e+13 +2.85900e-02 3.63837e+13 +2.95856e-02 3.61470e+13 +3.06159e-02 3.59118e+13 +3.16820e-02 3.56782e+13 +3.27853e-02 3.54461e+13 +3.39270e-02 3.52155e+13 +3.51084e-02 3.49864e+13 +3.63310e-02 3.47588e+13 +3.75962e-02 3.45326e+13 +3.89054e-02 3.43080e+13 +4.02602e-02 3.40848e+13 +4.16622e-02 3.38630e+13 +4.31130e-02 3.36427e+13 +4.46143e-02 3.34239e+13 +4.61679e-02 3.32064e+13 +4.77756e-02 3.29904e+13 +4.94393e-02 3.27758e+13 +5.11609e-02 3.25626e+13 +5.29425e-02 3.23507e+13 +5.47861e-02 3.21403e+13 +5.66939e-02 3.19312e+13 +5.86682e-02 3.17234e+13 +6.07112e-02 3.15170e+13 +6.28253e-02 3.13120e+13 +6.50131e-02 3.11083e+13 +6.72771e-02 3.09059e+13 +6.96199e-02 3.07049e+13 +7.20442e-02 3.05051e+13 +7.45530e-02 3.03067e+13 +7.71492e-02 3.01095e+13 +7.98358e-02 2.99136e+13 +8.26159e-02 2.97190e+13 +8.54928e-02 2.95257e+13 +8.84699e-02 2.93336e+13 +9.15507e-02 2.91427e+13 +9.47388e-02 2.89532e+13 +9.80379e-02 2.87648e+13 +1.01452e-01 2.85777e+13 +1.04985e-01 2.83917e+13 +1.08641e-01 2.82070e+13 +1.12424e-01 2.80235e+13 +1.16339e-01 2.78412e+13 +1.20390e-01 2.76601e+13 +1.24582e-01 2.74802e+13 +1.28921e-01 2.73014e+13 +1.33410e-01 2.71238e+13 +1.38056e-01 2.69473e+13 +1.42863e-01 2.67720e+13 +1.47838e-01 2.65978e+13 +1.52987e-01 2.64248e+13 +1.58314e-01 2.62529e+13 +1.63827e-01 2.60821e+13 +1.69532e-01 2.59124e+13 +1.75436e-01 2.57438e+13 +1.81545e-01 2.55764e+13 +1.87867e-01 2.54100e+13 +1.94409e-01 2.52447e+13 +2.01179e-01 2.50804e+13 +2.08184e-01 2.49173e+13 +2.15434e-01 2.47552e+13 +2.22936e-01 2.45941e+13 +2.30699e-01 2.44341e+13 +2.38733e-01 2.42752e+13 +2.47046e-01 2.41172e+13 +2.55649e-01 2.39603e+13 +2.64552e-01 2.38045e+13 +2.73764e-01 2.36496e+13 +2.83298e-01 2.34957e+13 +2.93163e-01 2.33429e+13 +3.03372e-01 2.31910e+13 +3.13936e-01 2.30401e+13 +3.24868e-01 2.28903e+13 +3.36181e-01 2.27413e+13 +3.47888e-01 2.25934e+13 +3.60002e-01 2.24464e+13 +3.72539e-01 2.23004e+13 +3.85512e-01 2.21553e+13 +3.98936e-01 2.20112e+13 +4.12829e-01 2.18640e+13 +4.27205e-01 2.17096e+13 +4.42081e-01 2.15476e+13 +4.57476e-01 2.13780e+13 +4.73406e-01 2.12003e+13 +4.89892e-01 2.10143e+13 +5.06951e-01 2.08198e+13 +5.24605e-01 2.06165e+13 +5.42873e-01 2.04043e+13 +5.61778e-01 2.01829e+13 +5.81341e-01 1.99521e+13 +6.01585e-01 1.97119e+13 +6.22534e-01 1.94621e+13 +6.44212e-01 1.92028e+13 +6.66646e-01 1.89339e+13 +6.89860e-01 1.86557e+13 +7.13883e-01 1.83681e+13 +7.38743e-01 1.80716e+13 +7.64468e-01 1.77666e+13 +7.91089e-01 1.74536e+13 +8.18637e-01 1.71333e+13 +8.47145e-01 1.68066e+13 +8.76645e-01 1.64745e+13 +9.07172e-01 1.61384e+13 +9.38763e-01 1.58000e+13 +9.71453e-01 1.54610e+13 +1.00528e+00 1.51237e+13 +1.04029e+00 1.47908e+13 +1.07652e+00 1.44653e+13 +1.11400e+00 1.41510e+13 +1.15280e+00 1.38519e+13 +1.19294e+00 1.35728e+13 +1.23448e+00 1.33195e+13 +1.27747e+00 1.30983e+13 +1.32196e+00 1.29163e+13 +1.36799e+00 1.27784e+13 +1.41563e+00 1.26854e+13 +1.46492e+00 1.26382e+13 +1.51594e+00 1.26365e+13 +1.56873e+00 1.26797e+13 +1.62335e+00 1.27665e+13 +1.67988e+00 1.28942e+13 +1.73838e+00 1.30596e+13 +1.79892e+00 1.32577e+13 +1.86156e+00 1.34822e+13 +1.92639e+00 1.37250e+13 +1.99347e+00 1.39760e+13 +2.06289e+00 1.41902e+13 +2.13473e+00 1.44034e+13 +2.20906e+00 1.46199e+13 +2.28599e+00 1.48396e+13 +2.36559e+00 1.50626e+13 +2.44797e+00 1.52890e+13 +2.53322e+00 1.55187e+13 +2.62143e+00 1.57520e+13 +2.71272e+00 1.59887e+13 +2.80718e+00 1.62290e+13 +2.90494e+00 1.64729e+13 +3.00610e+00 1.67204e+13 +3.11078e+00 1.69717e+13 +3.21911e+00 1.72268e+13 +3.33120e+00 1.74857e+13 +3.44721e+00 1.77485e+13 +3.56725e+00 1.80152e+13 +3.69147e+00 1.82859e+13 +3.82002e+00 1.85607e+13 +3.95304e+00 1.88397e+13 +4.09070e+00 1.91228e+13 +4.23315e+00 1.94102e+13 +4.38056e+00 1.97019e+13 +4.53311e+00 1.99980e+13 +4.69096e+00 2.02985e+13 +4.85432e+00 2.06036e+13 +5.02336e+00 2.09132e+13 +5.19829e+00 2.12275e+13 +5.37931e+00 2.15465e+13 +5.56663e+00 2.18704e+13 +5.76048e+00 2.21990e+13 +5.96108e+00 2.25326e+13 +6.16866e+00 2.28713e+13 +6.38347e+00 2.32150e+13 +6.60576e+00 2.35639e+13 +6.83580e+00 2.39180e+13 +7.07384e+00 2.42775e+13 +7.32017e+00 2.46423e+13 +7.57508e+00 2.50126e+13 +7.83887e+00 2.53886e+13 +8.11184e+00 2.57701e+13 +8.39432e+00 2.61574e+13 +8.68664e+00 2.65505e+13 +8.98913e+00 2.69495e+13 +9.30216e+00 2.73545e+13 +9.62609e+00 2.77656e+13 +9.96130e+00 2.81829e+13 +1.03082e+01 2.86064e+13 +1.06671e+01 2.90363e+13 +1.10386e+01 2.94727e+13 +1.14230e+01 2.99156e+13 +1.18208e+01 3.03652e+13 +1.22324e+01 3.08216e+13 +1.26584e+01 3.12848e+13 +1.30992e+01 3.17549e+13 +1.35554e+01 3.22322e+13 +1.40274e+01 3.27166e+13 +1.45159e+01 3.32082e+13 +1.50214e+01 3.37073e+13 +1.55444e+01 3.42139e+13 +1.60858e+01 3.47281e+13 +1.66459e+01 3.52500e+13 +1.72256e+01 3.57797e+13 +1.78254e+01 3.63174e+13 +1.84461e+01 3.68632e+13 +1.90885e+01 3.74172e+13 +1.97532e+01 3.79795e+13 +2.04411e+01 3.85503e+13 +2.11529e+01 3.91297e+13 +2.18895e+01 3.97177e+13 +2.26518e+01 4.03146e+13 +2.34406e+01 4.09205e+13 +2.42569e+01 4.15355e+13 +2.51015e+01 4.21597e+13 +2.59757e+01 4.27933e+13 +2.68802e+01 4.34364e+13 +2.78163e+01 4.40892e+13 +2.87849e+01 4.47518e+13 +2.97873e+01 4.54243e+13 +3.08246e+01 4.61070e+13 +3.18980e+01 4.67999e+13 +3.30088e+01 4.75032e+13 +3.41582e+01 4.82171e+13 +3.53477e+01 4.89418e+13 +3.65786e+01 4.96773e+13 +3.78524e+01 5.04238e+13 +3.91706e+01 5.11816e+13 +4.05346e+01 5.19508e+13 +4.19461e+01 5.27316e+13 +4.34068e+01 5.35240e+13 +4.49184e+01 5.43284e+13 +4.64826e+01 5.51449e+13 +4.81012e+01 5.59736e+13 +4.97763e+01 5.68148e+13 +5.15096e+01 5.76687e+13 +5.33033e+01 5.85353e+13 +5.51595e+01 5.94150e+13 +5.70804e+01 6.03079e+13 +5.90681e+01 6.12143e+13 +6.11250e+01 6.21342e+13 +6.32536e+01 6.30680e+13 +6.54562e+01 6.40158e+13 +6.77356e+01 6.49779e+13 +7.00944e+01 6.59544e+13 +7.25353e+01 6.69456e+13 +7.50612e+01 6.79517e+13 +7.76750e+01 6.89729e+13 +8.03799e+01 7.00095e+13 +8.31790e+01 7.10616e+13 +8.60755e+01 7.21295e+13 +8.90729e+01 7.32135e+13 +9.21747e+01 7.43138e+13 +9.53845e+01 7.54306e+13 +9.87061e+01 7.65643e+13 +1.02143e+02 7.76988e+13 +1.05700e+02 7.87849e+13 +1.09381e+02 7.98147e+13 +1.13190e+02 8.07834e+13 +1.17132e+02 8.16859e+13 +1.21211e+02 8.25171e+13 +1.25432e+02 8.32723e+13 +1.29799e+02 8.39467e+13 +1.34319e+02 8.45359e+13 +1.38997e+02 8.50360e+13 +1.43837e+02 8.54433e+13 +1.48846e+02 8.57551e+13 +1.54029e+02 8.59691e+13 +1.59393e+02 8.60845e+13 +1.64944e+02 8.61013e+13 +1.70687e+02 8.60211e+13 +1.76631e+02 8.58472e+13 +1.82782e+02 8.55850e+13 +1.89147e+02 8.52425e+13 +1.95734e+02 8.48302e+13 +2.02550e+02 8.43622e+13 +2.09603e+02 8.38564e+13 +2.16902e+02 8.33352e+13 +2.24455e+02 8.28262e+13 +2.32272e+02 8.23630e+13 +2.40360e+02 8.19859e+13 +2.48730e+02 8.17252e+13 +2.57392e+02 8.15492e+13 +2.66355e+02 8.14021e+13 +2.75630e+02 8.12151e+13 +2.85228e+02 8.09042e+13 +2.95161e+02 8.03680e+13 +3.05439e+02 7.94880e+13 +3.16076e+02 7.81892e+13 +3.27082e+02 7.64501e+13 +3.38472e+02 7.42527e+13 +3.50259e+02 7.15810e+13 +3.62456e+02 6.84213e+13 +3.75078e+02 6.47632e+13 +3.88139e+02 6.06001e+13 +4.01656e+02 5.59301e+13 +4.15642e+02 5.07567e+13 +4.30116e+02 4.50903e+13 +4.45094e+02 3.89494e+13 +4.60594e+02 3.23615e+13 +4.76633e+02 2.53654e+13 +4.93231e+02 1.80128e+13 +5.10407e+02 1.43946e+13 +5.28181e+02 1.34421e+13 +5.46574e+02 1.25526e+13 +5.65607e+02 1.17220e+13 +5.85303e+02 1.09464e+13 +6.05685e+02 1.02220e+13 +6.26777e+02 9.54565e+12 +6.48603e+02 8.91401e+12 +6.71189e+02 8.32417e+12 +6.94562e+02 7.77336e+12 +7.18749e+02 7.25900e+12 +7.43778e+02 6.77867e+12 +7.69679e+02 6.33013e+12 +7.96481e+02 5.91126e+12 +8.24217e+02 5.52011e+12 +8.52919e+02 5.15485e+12 +8.82620e+02 4.81375e+12 +9.13356e+02 4.49522e+12 +9.45161e+02 4.19778e+12 +9.78075e+02 3.92001e+12 +1.01213e+03 3.66062e+12 +1.04738e+03 3.41840e+12 +1.08385e+03 3.19220e+12 +1.12160e+03 2.98098e+12 +1.16065e+03 2.78372e+12 +1.20107e+03 2.59952e+12 +1.24290e+03 2.42751e+12 +1.28618e+03 2.26689e+12 +1.33097e+03 2.11689e+12 +1.37731e+03 1.97681e+12 +1.42528e+03 1.84601e+12 +1.47491e+03 1.72386e+12 +1.52627e+03 1.60979e+12 +1.57942e+03 1.50327e+12 +1.63442e+03 1.40380e+12 +1.69133e+03 1.31091e+12 +1.75023e+03 1.22417e+12 +1.81118e+03 1.14316e+12 +1.87425e+03 1.06752e+12 +1.93952e+03 9.96881e+11 +2.00706e+03 9.30918e+11 +2.07695e+03 8.69319e+11 +2.14928e+03 8.11796e+11 +2.22412e+03 7.58080e+11 +2.30157e+03 7.07917e+11 +2.38172e+03 6.61075e+11 +2.46466e+03 6.17331e+11 +2.55048e+03 5.76482e+11 +2.63930e+03 5.38337e+11 +2.73121e+03 5.02715e+11 +2.82632e+03 4.69450e+11 +2.92474e+03 4.38387e+11 +3.02659e+03 4.09379e+11 +3.13198e+03 3.82290e+11 +3.24105e+03 3.56994e+11 +3.35391e+03 3.33372e+11 +3.47070e+03 3.11312e+11 +3.59156e+03 2.90713e+11 +3.71663e+03 2.71476e+11 +3.84606e+03 2.53513e+11 +3.97999e+03 2.36738e+11 +4.11858e+03 2.21073e+11 +4.26200e+03 2.06445e+11 +4.41042e+03 1.92784e+11 +4.56401e+03 1.80028e+11 +4.72294e+03 1.68115e+11 +4.88740e+03 1.56991e+11 +5.05760e+03 1.46603e+11 +5.23372e+03 1.36902e+11 +5.41597e+03 1.27843e+11 +5.60457e+03 1.19384e+11 +5.79974e+03 1.11484e+11 +6.00171e+03 1.04107e+11 +6.21071e+03 9.72186e+10 +6.42698e+03 9.07857e+10 +6.65079e+03 8.47784e+10 +6.88239e+03 7.91686e+10 +7.12205e+03 7.39300e+10 +7.37007e+03 6.90380e+10 +7.62671e+03 6.44698e+10 +7.89230e+03 6.02038e+10 +8.16713e+03 5.62201e+10 +8.45154e+03 5.25001e+10 +8.74585e+03 4.90261e+10 +9.05040e+03 4.57821e+10 +9.36557e+03 4.27527e+10 +9.69170e+03 3.99237e+10 +1.00292e+04 3.72820e+10 +1.03784e+04 3.48150e+10 +1.07399e+04 3.25113e+10 +1.11138e+04 3.03600e+10 +1.15009e+04 2.83511e+10 +1.19014e+04 2.64751e+10 +1.23158e+04 2.47233e+10 +1.27447e+04 2.30873e+10 +1.31885e+04 2.15596e+10 +1.36477e+04 2.01330e+10 +1.41230e+04 1.88008e+10 +1.46148e+04 1.75568e+10 +1.51237e+04 1.63950e+10 +1.56504e+04 1.53102e+10 +1.61954e+04 1.42971e+10 +1.67594e+04 1.33511e+10 +1.73430e+04 1.24676e+10 +1.79469e+04 1.16427e+10 +1.85719e+04 1.08723e+10 +1.92186e+04 1.01528e+10 +1.98879e+04 9.48103e+09 +2.05804e+04 8.85367e+09 +2.12971e+04 8.26782e+09 +2.20387e+04 7.72074e+09 +2.28062e+04 7.20986e+09 +2.36003e+04 6.73278e+09 +2.44222e+04 6.28727e+09 +2.52726e+04 5.87124e+09 +2.61527e+04 5.48274e+09 +2.70634e+04 5.11995e+09 +2.80059e+04 4.78116e+09 +2.89811e+04 4.46479e+09 +2.99903e+04 4.16936e+09 +3.10347e+04 3.89347e+09 +3.21154e+04 3.63584e+09 +3.32337e+04 3.39526e+09 +3.43910e+04 3.17059e+09 +3.55886e+04 2.96079e+09 +3.68280e+04 2.76488e+09 +3.81104e+04 2.58193e+09 +3.94375e+04 2.41108e+09 +4.08109e+04 2.25154e+09 +4.22320e+04 2.10256e+09 +4.37027e+04 1.96343e+09 +4.52245e+04 1.83351e+09 +4.67994e+04 1.71219e+09 +4.84291e+04 1.59889e+09 +5.01155e+04 1.49309e+09 +5.18607e+04 1.39429e+09 +5.36667e+04 1.30203e+09 +5.55355e+04 1.21588e+09 +5.74694e+04 1.13542e+09 +5.94707e+04 1.06029e+09 +6.15416e+04 9.90133e+08 +6.36847e+04 9.24616e+08 +6.59024e+04 8.63434e+08 +6.81973e+04 8.06300e+08 +7.05721e+04 7.52947e+08 +7.30297e+04 7.03125e+08 +7.55728e+04 6.56599e+08 +7.82045e+04 6.13152e+08 +8.09278e+04 5.72580e+08 +8.37459e+04 5.34692e+08 +8.66622e+04 4.99312e+08 +8.96801e+04 4.66272e+08 +9.28030e+04 4.35419e+08 +9.60347e+04 4.06607e+08 +9.93789e+04 3.79702e+08 +1.02840e+05 3.54577e+08 +1.06421e+05 3.31115e+08 +1.10127e+05 3.09205e+08 +1.13962e+05 2.88745e+08 +1.17930e+05 2.69639e+08 +1.22037e+05 2.51797e+08 +1.26286e+05 2.35135e+08 +1.30684e+05 2.19576e+08 +1.35235e+05 2.05047e+08 +1.39944e+05 1.91479e+08 +1.44818e+05 1.78809e+08 +1.49861e+05 1.66977e+08 +1.55079e+05 1.55928e+08 +1.60479e+05 1.45610e+08 +1.66068e+05 1.35975e+08 +1.71851e+05 1.26978e+08 +1.77835e+05 1.18576e+08 +1.84028e+05 1.10730e+08 +1.90436e+05 1.03403e+08 +1.97068e+05 9.65605e+07 +2.03930e+05 9.01711e+07 +2.11032e+05 8.42044e+07 +2.18381e+05 7.86326e+07 +2.25985e+05 7.34295e+07 +2.33855e+05 6.85707e+07 +2.41998e+05 6.40334e+07 +2.50426e+05 5.97963e+07 +2.59146e+05 5.58395e+07 +2.68170e+05 5.21446e+07 +2.77509e+05 4.86942e+07 +2.87173e+05 4.54721e+07 +2.97173e+05 4.24632e+07 +3.07521e+05 3.96535e+07 +3.18230e+05 3.70296e+07 +3.29312e+05 3.45793e+07 +3.40779e+05 3.22912e+07 +3.52646e+05 3.01545e+07 +3.64927e+05 2.81592e+07 +3.77635e+05 2.62959e+07 +3.90785e+05 2.45559e+07 +4.04393e+05 2.29310e+07 +4.18475e+05 2.14137e+07 +4.33048e+05 1.99967e+07 +4.48128e+05 1.86736e+07 +4.63733e+05 1.74379e+07 +4.79882e+05 1.62841e+07 +4.96593e+05 1.52065e+07 +5.13886e+05 1.42003e+07 +5.31781e+05 1.32607e+07 +5.50299e+05 1.23832e+07 +5.69462e+05 1.15638e+07 +5.89292e+05 1.07987e+07 +6.09813e+05 1.00841e+07 +6.31049e+05 9.41684e+06 +6.53024e+05 8.79373e+06 +6.75764e+05 8.21185e+06 +6.99296e+05 7.66847e+06 +7.23648e+05 7.16105e+06 +7.48848e+05 6.68720e+06 +7.74925e+05 6.24471e+06 +8.01910e+05 5.83150e+06 +8.29835e+05 5.44563e+06 +8.58732e+05 5.08529e+06 +8.88636e+05 4.74880e+06 +9.19581e+05 4.43457e+06 +9.51604e+05 4.14113e+06 +9.84741e+05 3.86711e+06 +1.01903e+06 3.61123e+06 +1.05452e+06 3.37227e+06 +1.09124e+06 3.14913e+06 +1.12924e+06 2.94075e+06 +1.16856e+06 2.74616e+06 +1.20926e+06 2.56445e+06 +1.25137e+06 2.39476e+06 +1.29494e+06 2.23630e+06 +1.34004e+06 2.08832e+06 +1.38670e+06 1.95014e+06 +1.43499e+06 1.82110e+06 +1.48496e+06 1.70059e+06 +1.53667e+06 1.58807e+06 +1.59018e+06 1.48298e+06 +1.64556e+06 1.38485e+06 +1.70286e+06 1.29322e+06 +1.76216e+06 1.20765e+06 +1.82353e+06 1.12774e+06 +1.88703e+06 1.05311e+06 +1.95274e+06 9.83430e+05 +2.02074e+06 9.18356e+05 +2.09111e+06 8.57589e+05 +2.16393e+06 8.00842e+05 +2.23928e+06 7.47850e+05 +2.31726e+06 6.98365e+05 +2.39795e+06 6.52154e+05 +2.48146e+06 6.09001e+05 +2.56787e+06 5.68703e+05 +2.65729e+06 5.31072e+05 +2.74982e+06 4.95931e+05 +2.84558e+06 4.63116e+05 +2.94467e+06 4.32471e+05 +3.04722e+06 4.03855e+05 +3.15333e+06 3.77131e+05 +3.26314e+06 3.52177e+05 +3.37677e+06 3.28873e+05 +3.49436e+06 3.07112e+05 +3.61604e+06 2.86790e+05 +3.74196e+06 2.67813e+05 +3.87227e+06 2.50092e+05 +4.00712e+06 2.33543e+05 +4.14666e+06 2.18090e+05 +4.29105e+06 2.03659e+05 +4.44048e+06 1.90183e+05 +4.59511e+06 1.77598e+05 +4.75513e+06 1.65847e+05 +4.92072e+06 1.54873e+05 +5.09207e+06 1.44625e+05 +5.26939e+06 1.35055e+05 +5.45289e+06 1.26118e+05 +5.64278e+06 1.17773e+05 +5.83927e+06 1.09980e+05 +6.04262e+06 1.02703e+05 +6.25304e+06 9.59068e+04 +6.47079e+06 8.95606e+04 +6.69612e+06 8.36344e+04 +6.92930e+06 7.81003e+04 +7.17060e+06 7.29324e+04 +7.42030e+06 6.81065e+04 +7.67870e+06 6.35999e+04 +7.94609e+06 5.93915e+04 +8.22280e+06 5.54615e+04 +8.50914e+06 5.17916e+04 +8.80546e+06 4.83646e+04 +9.11209e+06 4.51643e+04 +9.42940e+06 4.21758e+04 +9.75776e+06 3.93850e+04 +1.00976e+07 3.67789e+04 +1.04492e+07 3.43452e+04 +1.08131e+07 3.20726e+04 +1.11896e+07 2.99504e+04 +1.15793e+07 2.79686e+04 +1.19825e+07 2.61179e+04 +1.23997e+07 2.43897e+04 +1.28315e+07 2.27758e+04 +1.32784e+07 2.12687e+04 +1.37408e+07 1.98614e+04 +1.42193e+07 1.85471e+04 +1.47144e+07 1.73199e+04 +1.52268e+07 1.61738e+04 +1.57571e+07 1.51036e+04 +1.63058e+07 1.41042e+04 +1.68736e+07 1.31709e+04 +1.74612e+07 1.22994e+04 +1.80692e+07 1.14855e+04 +1.86985e+07 1.07255e+04 +1.93496e+07 1.00158e+04 +2.00234e+07 9.35309e+03 +2.07207e+07 8.73420e+03 +2.14422e+07 8.15625e+03 +2.21889e+07 7.61656e+03 +2.29616e+07 7.11257e+03 +2.37612e+07 6.64193e+03 +2.45886e+07 6.20243e+03 +2.54449e+07 5.79202e+03 +2.63310e+07 5.40876e+03 +2.72479e+07 5.05086e+03 +2.81967e+07 4.71665e+03 +2.91786e+07 4.40455e+03 +3.01947e+07 4.11310e+03 +3.12462e+07 3.84093e+03 +3.23343e+07 3.58678e+03 +3.34603e+07 3.34944e+03 +3.46255e+07 3.12781e+03 +3.58312e+07 2.92084e+03 +3.70790e+07 2.72757e+03 +3.83702e+07 2.54709e+03 +3.97063e+07 2.37855e+03 +4.10890e+07 2.22116e+03 +4.25199e+07 2.07418e+03 +4.40006e+07 1.93693e+03 +4.55328e+07 1.80877e+03 +4.71184e+07 1.68908e+03 +4.87592e+07 1.57732e+03 +5.04571e+07 1.47294e+03 +5.22142e+07 1.37548e+03 +5.40325e+07 1.28446e+03 +5.59140e+07 1.19947e+03 +5.78611e+07 1.12010e+03 +5.98760e+07 1.04598e+03 +6.19611e+07 9.76772e+02 +6.41188e+07 9.12139e+02 +6.63516e+07 8.51783e+02 +6.86621e+07 7.95420e+02 +7.10532e+07 7.42787e+02 +7.35274e+07 6.93637e+02 +7.60879e+07 6.47739e+02 +7.87375e+07 6.04878e+02 +8.14794e+07 5.64853e+02 +8.43167e+07 5.27477e+02 +8.72529e+07 4.92574e+02 +9.02913e+07 4.59980e+02 +9.34355e+07 4.29543e+02 +9.66892e+07 4.01121e+02 diff --git a/sed/mehdipour2013/mehdipour2013_powlaw_coords.tab b/sed/mehdipour2013/mehdipour2013_powlaw_coords.tab index 37166e3..c07a870 100644 --- a/sed/mehdipour2013/mehdipour2013_powlaw_coords.tab +++ b/sed/mehdipour2013/mehdipour2013_powlaw_coords.tab @@ -1,4 +1,3 @@ -12.00 -9 0.0004 3.2e10 0.006 4.9e13 .4 2.2e13 diff --git a/src/agn.hpp b/src/agn.hpp index 31c01c7..e61ad2f 100644 --- a/src/agn.hpp +++ b/src/agn.hpp @@ -17,7 +17,7 @@ namespace agn { -const bool debug = true; +const bool debug = false; const bool line_debug = false; const bool verbose=true; diff --git a/src/sed.hpp b/src/sed.hpp index ba5ed6e..8166db9 100644 --- a/src/sed.hpp +++ b/src/sed.hpp @@ -37,6 +37,11 @@ struct sed_table { // To account for the four main powerlaws in a typical // AGN SED. + +// Hardcoded infrared and gamma ray power laws. +const double IR_POWER = 3; +const double GAMMA_POWER = -2; + struct powerlaw_bounds { double ir_min; double ir_max; @@ -183,10 +188,8 @@ agn::sed_powerlaw_spline::sed_powerlaw_spline( // powerlaws are evaluated across four regions of the sed, first // we construct the powerlaws, here, and locate them iterator1d table_it = powerlaw_coords.table.begin(); - double ir_power = (*table_it).first; - double gamma_power = (*table_it).second; - std::cout << ir_power << ", " << gamma_power; - table_it++; + double ir_power = agn::IR_POWER; + double gamma_power = agn::GAMMA_POWER; coord2d ir_high_point = *table_it; table_it++; coord2d uv_low_point = *table_it; table_it++; coord2d uv_high_point = *table_it; table_it++; @@ -385,34 +388,17 @@ double agn::sed_pow_law::SED_at_2KeV() { agn::sed_table agn::read_sed_table(std::ifstream& table_file) { sed_table resultant; - std::string scratch; + std::string line; double hnu; - std::getline(table_file,scratch); - if(!isdigit(scratch[0])) { - resultant.header = scratch; + if(!isdigit(table_file.peek()) && table_file.peek() != '#') { + std::getline(table_file,resultant.header); } - else - table_file.seekg(0); while(!table_file.eof()) { - char next; - next = table_file.peek(); - if (next == '\n'){ - table_file.get(); - next = table_file.peek(); - } - if (next == '#') { - char commented; - while (commented != '\n') { - if(table_file.eof() || table_file.fail()) - break; - table_file.get(commented); - } - if(table_file.eof() || table_file.fail()) - break; - next = table_file.peek(); - } - table_file >> hnu; - table_file >> resultant.table[hnu]; + std::getline(table_file,line); + if (line[0] == '#') continue; + std::stringstream scratch(line); + scratch >> hnu; + scratch >> resultant.table[hnu]; } return resultant; } diff --git a/src/table_powerlaw_spline_sed.cpp b/src/table_powerlaw_spline_sed.cpp index 9f8381f..63347ba 100644 --- a/src/table_powerlaw_spline_sed.cpp +++ b/src/table_powerlaw_spline_sed.cpp @@ -3,7 +3,9 @@ // Syntax: table_powerlaw_spline - +namespace agn { + bool sed_debug =true; +} int main(int argc, char const *argv[]) { @@ -46,12 +48,14 @@ int main(int argc, char const *argv[]) // Read in sampling table and construct a spline model. samples = agn::read_sed_table(sample_table); - if(agn::debug) debug_file - << "Read samples:\n" + if(agn::sed_debug) debug_file + << "Read " + << samples.table.size() + << " samples:\n" << format_sed_table(samples); powerlaw_coords = agn::read_sed_table(powerlaw_table); - if(agn::debug) debug_file + if(agn::sed_debug) debug_file << "Read power coords:\n" << format_sed_table(powerlaw_coords);