Adding EOSs from master to lalinferenceO2

lalsimulation/src/LALSimNeutronStar.h

@@ -55,6 +55,9 @@
 /** Incomplete type for the neutron star Equation of State (EOS). */
 typedef struct tagLALSimNeutronStarEOS LALSimNeutronStarEOS;
 
+/** Recognised names of equations of state */
+extern const char * const lalSimNeutronStarEOSNames[61];
+
 /** Incomplete type for a neutron star family having a particular EOS. */
 typedef struct tagLALSimNeutronStarFamily LALSimNeutronStarFamily;
 

lalsimulation/src/LALSimNeutronStarEOS.c

@@ -77,6 +77,30 @@ struct tagLALSimNeutronStarEOS {
 
 /** @endcond */
 
+/** Recognised equations of state names */
+const char * const lalSimNeutronStarEOSNames[61] = {
+        "ALF1", "ALF2", "ALF3", "ALF4",
+        "AP1", "AP2", "AP3", "AP4",
+        "BBB2", "BGN1H1", "BPAL12", 
+        "BSK19", "BSK20", "BSK21",
+        "ENG", "FPS", "GNH3",
+        "GS1", "GS2",
+        "H1", "H2", "H3", "H4", "H5", "H6", "H7",
+        "MPA1", "MS1B", "MS1", "MS2",
+        "PAL6", "PCL2", "PS",
+        "QMC700",
+        "SLY4", "SLY",
+        "SQM1", "SQM2", "SQM3",
+        "WFF1", "WFF2", "WFF3",
+        /* From here, EOSs are coming from CompOSE via Ian */
+        "APR", "BHF_BBB2",
+        "KDE0V", "KDE0V1", "RS", "SK255", "SK272",
+        "SKA", "SKB", "SKI2", "SKI3", "SKI4", "SKI5", "SKI6",
+        "SKMP", "SKOP",
+        "SLY2", "SLY230A", "SLY9"
+    };
+
+
 /**
  * @name Destruction routine
  * @{

lalsimulation/src/LALSimNeutronStarEOSTabular.c

@@ -26,6 +26,7 @@
  */
 
 #include <lal/LALSimReadData.h>
+#include <lal/LALSimNeutronStar.h>
 #include <gsl/gsl_interp.h>
 
 void GLBoundConversion(double a, double b, double abcissae[],int nEval);
@@ -304,17 +305,6 @@ static LALSimNeutronStarEOS *eos_alloc_tabular(double *pdat, double *edat,
     return eos;
 }
 
-static int mystrcasecmp(const char *s1, const char *s2)
-{
-    while (*s1) {
-        int c1 = toupper(*s1++);
-        int c2 = toupper(*s2++);
-        if (c1 != c2)
-            return (c1 > c2) - (c1 < c2);
-    }
-    return 0;
-}
-
 /** @endcond */
 
 /**
@@ -357,9 +347,69 @@ LALSimNeutronStarEOS *XLALSimNeutronStarEOSFromFile(const char *fname)
  * @details A known, installed, named tabulated equation of state data file is
  * read and the used to create the equation of state structure.  Presently
  * the known equations of state are:
+ * - ALF1
+ * - ALF2
+ * - ALF3
+ * - ALF4
+ * - AP1
+ * - AP2
+ * - AP3
  * - AP4
+ * - BBB2
+ * - BGN1H1
+ * - BPAL12
+ * - BSK19
+ * - BSK20
+ * - BSK21
+ * - ENG
  * - FPS
+ * - GNH3
+ * - GS1
+ * - GS2
+ * - H1
+ * - H2
+ * - H3
+ * - H4
+ * - H5
+ * - H6
+ * - H7
+ * - MPA1
+ * - MS1B
+ * - MS1
+ * - MS2
+ * - PAL6
+ * - PCL2
+ * - PS
+ * - QMC700
  * - SLY4
+ * - SLY
+ * - SQM1
+ * - SQM2
+ * - SQM3
+ * - WFF1
+ * - WFF2
+ * - WFF3
+ * We also include more modern equations from the CompOSE website
+ * https://compose.obspm.fr/ downloaded on 18 June 2018. These EOSs are:
+ * - APR
+ * - BHF_BBB2
+ * - KDE0V
+ * - KDE0V1
+ * - RS
+ * - SK255
+ * - SK272
+ * - SKA
+ * - SKB
+ * - SKI2
+ * - SKI3
+ * - SKI4
+ * - SKI5
+ * - SKI6
+ * - SKMP
+ * - SKOP
+ * - SLY2
+ * - SLY230A
+ * - SLY9
  * @param[in] name The name of the equation of state.
  * @return A pointer to neutron star equation of state structure.
  */
@@ -367,31 +417,26 @@ LALSimNeutronStarEOS *XLALSimNeutronStarEOSByName(const char *name)
 {
     static const char fname_base[] = "LALSimNeutronStarEOS_";
     static const char fname_extn[] = ".dat";
-    static const char *eos_names[] = {
-        "FPS",
-        "SLY4",
-        "AP4"
-    };
-    size_t n = XLAL_NUM_ELEM(eos_names);
+    size_t n = XLAL_NUM_ELEM(lalSimNeutronStarEOSNames);
     size_t i;
     char fname[FILENAME_MAX];
 
     for (i = 0; i < n; ++i)
-        if (mystrcasecmp(name, eos_names[i]) == 0) {
+        if (XLALStringCaseCompare(name, lalSimNeutronStarEOSNames[i]) == 0) {
             LALSimNeutronStarEOS *eos;
-            snprintf(fname, sizeof(fname), "%s%s%s", fname_base, eos_names[i],
+            snprintf(fname, sizeof(fname), "%s%s%s", fname_base, lalSimNeutronStarEOSNames[i],
                 fname_extn);
             eos = XLALSimNeutronStarEOSFromFile(fname);
             if (!eos)
                 XLAL_ERROR_NULL(XLAL_EFUNC);
-            snprintf(eos->name, sizeof(eos->name), "%s", eos_names[i]);
+            snprintf(eos->name, sizeof(eos->name), "%s", lalSimNeutronStarEOSNames[i]);
             return eos;
         }
 
     XLAL_PRINT_ERROR("Unrecognized EOS name %s...", name);
-    XLALPrintError("\tKnown EOS names are: %s", eos_names[0]);
+    XLALPrintError("\tKnown EOS names are: %s", lalSimNeutronStarEOSNames[0]);
     for (i = 1; i < n; ++i)
-        XLALPrintError(", %s", eos_names[i]);
+        XLALPrintError(", %s", lalSimNeutronStarEOSNames[i]);
     XLALPrintError("\n");
     XLAL_ERROR_NULL(XLAL_ENAME);
 }

lalsimulation/src/LALSimNeutronStarEOS_AP1.dat

+8.343631758107906939e-38	5.760646279117794535e-24
+8.343631580857245259e-37	5.804959210497264011e-24
+8.343631580857246261e-36	5.834500972728693026e-24
+8.343632038754789638e-35	6.019136798346299425e-24
+9.995903227680446623e-34	8.567115626089990214e-24
+1.156545764779204987e-32	1.211212891068081416e-23
+1.404386579748772230e-31	3.330835445636685513e-23
+4.807924266210808738e-30	1.565714166351965477e-22
+1.569554700758486418e-28	8.493260962020868104e-22
+8.049543353278943375e-27	7.710404093022768032e-21
+4.104091370047027957e-26	1.936463467731056542e-20
+2.008249911992246865e-25	4.864791994918719978e-20
+9.508463995767528335e-25	1.221552490862560189e-19
+4.350321992736018382e-24	3.069390689624394171e-19
+1.914897940977587221e-23	7.710403472645451774e-19
+8.058701481400324577e-23	1.936463526814610684e-18
+3.230836530483839324e-22	4.865530760907105231e-18
+4.344561346220565862e-22	6.124748887463802070e-18
+1.185467223548028312e-21	1.222290965864442148e-17
+3.166435353222566329e-21	2.438673691541350992e-17
+8.310692455186995649e-21	4.866269349644829091e-17
+2.151232156996131378e-20	9.711859588326032257e-17
+5.515154112228532168e-20	1.937940497497618577e-16
+7.218533151432849524e-20	2.440150869016797848e-16
+1.345731328335471147e-19	3.867756947236045186e-16
+2.502336294835788531e-19	6.130657183780712154e-16
+3.411041192947317021e-19	7.717789064604915480e-16
+4.160220775893210012e-19	9.719244456511960121e-16
+5.667146953719447417e-19	1.223768099027223581e-15
+1.050801028651695126e-18	1.939417719285730906e-15
+1.946359918714167429e-18	3.075298853003307365e-15
+3.603505884888307650e-18	4.875131897516413091e-15
+4.677349567547897744e-18	6.138781246310792288e-15
+6.362707709195922889e-18	7.725174671335268573e-15
+8.657661113518540592e-18	9.734015669972665520e-15
+1.177195562914600537e-17	1.225245236621271518e-14
+1.601016575098847113e-17	1.543557848371115445e-14
+2.067776666506359893e-17	1.943110352785689595e-14
+2.812081797163315919e-17	2.446797680216988522e-14
+3.823149159380315039e-17	3.081207326570878991e-14
+4.914569977792957113e-17	3.880312335400026662e-14
+6.682349839757224197e-17	4.886948667400894403e-14
+9.087197812587819867e-17	6.153552400687944153e-14
+1.235023519019284382e-16	7.747331269963000231e-14
+1.679450037352667011e-16	9.763557314036988077e-14
+2.145471643418675691e-16	1.228937920342251212e-13
+2.389043447093903019e-16	1.361875875494194413e-13
+2.717843425129702121e-16	1.547989144463240331e-13
+3.695233117251393371e-16	1.949757237840322556e-13
+4.804674567344333288e-16	2.455660272401238043e-13
+6.447788056504376238e-16	3.175002571949458142e-13
+6.517949865510213109e-16	3.293908150529345195e-13
+6.899629978881484835e-16	3.861110313286516761e-13
+7.515870871085186872e-16	4.881778664403938197e-13
+8.121477846155580033e-16	5.881767815932402550e-13
+8.946727017443121825e-16	7.184559915858633373e-13
+1.006192902542618365e-15	8.832991116900644424e-13
+1.155718643767072154e-15	1.086398862554044061e-12
+1.041020880690464571e-14	6.226173798963438979e-12
+2.143842591667660228e-14	1.246447412811411696e-11
+3.880280312113756540e-14	1.871067899484997050e-11
+5.992176115942451360e-14	2.496343966934196134e-11
+8.412503846984265287e-14	3.122191608684835906e-11
+1.111429062208144745e-13	3.748553942045298565e-11
+1.408637360812523970e-13	4.375390480010202386e-11
+1.732619623019208968e-13	5.002671562635437316e-11
+2.083647557845088832e-13	5.630375048359139163e-11
+2.462337999071061373e-13	6.258484260848191244e-11
+2.869544408634723667e-13	6.886986644847367973e-11
+3.306288577519030857e-13	7.515872761758915987e-11
+3.773716254003242311e-13	8.145135713577904988e-11
+4.273067380322587241e-13	8.774770537285798841e-11
+4.805655841780131941e-13	9.404773924203574316e-11
+5.364373251680912945e-13	1.003514389503217084e-10
+5.860202303949822085e-13	1.054715293425258458e-10
+6.397862645392009531e-13	1.117048442900149422e-10
+6.922967731608212056e-13	1.179352050598038139e-10
+7.440244246918799919e-13	1.241626116518924353e-10
+7.955305128954214507e-13	1.303870640662808579e-10
+8.473320188689860128e-13	1.366085623029690042e-10
+8.999311528216132770e-13	1.428300605396571764e-10
+9.538744376278458433e-13	1.490530358651954739e-10
+1.009678854385224647e-12	1.552715799241834215e-10
+1.067846613302788682e-12	1.614901239831713433e-10
+1.128924237255081254e-12	1.677086680421592909e-10
+1.193428707339642593e-12	1.739272121011472385e-10
+1.261877004654013330e-12	1.801457561601351603e-10
+1.334800881184235336e-12	1.863790711076242566e-10
+1.412702547139348006e-12	1.925976151666122043e-10
+1.496143296282395076e-12	1.988309301141013006e-10
+1.585654880599418616e-12	2.050642450615903969e-10
+1.681665655856952247e-12	2.112975600090794933e-10
+1.784766457595041974e-12	2.175308749565685896e-10
+1.895400412468723058e-12	2.237789607925588605e-10
+2.014306064903052854e-12	2.300270466285491055e-10
+2.141778832668055070e-12	2.362899033530405510e-10
+2.278409551303775914e-12	2.425527600775319447e-10
+2.424789056350261592e-12	2.488156168020233384e-10
+2.581360474462546758e-12	2.550932444150159326e-10
+2.748714641180678429e-12	2.613708720280085267e-10
+2.927294683159690854e-12	2.676632705295022696e-10
+3.117691435939630240e-12	2.739704399194971612e-10
+3.320348026175531647e-12	2.802776093094920527e-10
+3.535855289407441684e-12	2.865995495879880930e-10
+3.996263883988567813e-12	2.992582010334813739e-10
+4.497292421947881154e-12	3.117839144824641614e-10
+5.040122574365474121e-12	3.244425659279574423e-10
+5.626231430091463039e-12	3.371012173734507233e-10
+6.260641091216242279e-12	3.497598688189439525e-10
+6.946010317670019772e-12	3.624185202644372334e-10
+7.682486818337807878e-12	3.750771717099305144e-10
+8.475240404195010905e-12	3.878687611519341853e-10
+9.325452746321723690e-12	4.006455797054367076e-10
+1.023725969349826806e-11	4.134223982589392816e-10
+1.121051353683963571e-11	4.261992168124418039e-10
+1.225038408732122931e-11	4.391089733624547679e-10
+1.335583738274796802e-11	4.518857919159572901e-10
+1.453307719629034023e-11	4.647955484659702542e-10
+1.578860271888885582e-11	4.777053050159832699e-10
+1.712389103939362872e-11	4.907332286740054235e-10
+1.852269418045338821e-11	5.037611523320276806e-10
+2.000126011941940177e-11	5.167743051015488406e-10
+2.157140556709260322e-11	5.298022287595709943e-10
+2.321835963497183416e-11	5.429630904141037448e-10
+2.494507650075731916e-11	5.559910140721258985e-10
+2.676337287524999205e-11	5.692552719461667416e-10
+2.865995495879880994e-11	5.824013627121981883e-10
+3.066141035070585215e-11	5.956803914747401285e-10
+3.274262854051915487e-11	6.090775873452913616e-10
+3.492724294984056963e-11	6.223566161078333018e-10
+3.720195977901905998e-11	6.357390410898833346e-10
+3.956973320575484404e-11	6.492691749569450612e-10
+4.204237994084886534e-11	6.626663708274961910e-10
+4.460660618465006807e-11	6.762999009140660102e-10
+4.727570573680949510e-11	6.898152638926264331e-10
+5.004820150847704062e-11	7.034635648676973493e-10
+5.292557058850281691e-11	7.172300329507775586e-10
+5.590929006573692979e-11	7.308635630373473779e-10
+5.900822247328009176e-11	7.447481982284367769e-10
+6.221350527803159678e-11	7.585146663115169862e-10
+6.552366139114133257e-11	7.725322394991168786e-10
+6.895050752341021035e-11	7.864168746902063809e-10
+7.249552076368836823e-11	8.004344478778062733e-10
+7.615870111197580624e-11	8.145554172849142584e-10
+7.993857147942239915e-11	8.286911575805234438e-10
+8.383660895487825924e-11	8.428268978761325258e-10
+8.786610733799443588e-11	8.570808052797510043e-10
+9.199900194061871807e-11	8.714676506798798728e-10
+9.627665125055436615e-11	8.858397251915075409e-10
+1.006724676684992814e-10	9.003447376996457025e-10
+1.051982679052544009e-10	9.148497502077838641e-10
+1.098422352500187876e-10	9.294729298239312153e-10
+1.146294802132444279e-10	9.440961094400786699e-10
+1.245850590630255829e-10	9.735935737768931588e-10
+1.351270421863027758e-10	1.003327372329726234e-09
+1.461165832311650640e-10	1.033445213983589485e-09
+1.577530891923781359e-10	1.063799389853471116e-09
+1.699243013173331477e-10	1.094522837935881930e-09
+1.827454325363391419e-10	1.125497391122812323e-09
+1.960835448528857286e-10	1.156708278525761303e-09
+2.100863471519844722e-10	1.188406605249248448e-09
+2.247390685451341984e-10	1.220370807965756580e-09
+2.399087710358244653e-10	1.252320239793763512e-09
+2.557283926205657922e-10	1.285495655367366562e-09
+2.722127041878592245e-10	1.318685841829471018e-09
+2.893469348492036650e-10	1.351861257403074068e-09
+3.072787934896107044e-10	1.385051443865178318e-09
+3.257128623390571099e-10	1.419452843184377693e-09
+3.448116211710557241e-10	1.453869013392078061e-09
+3.647080079821168854e-10	1.488314725376780830e-09
+3.852395429987279063e-10	1.523912566664574129e-09
+4.064357679978910326e-10	1.559510407952367220e-09
+4.282819120911052189e-10	1.595255958125172110e-09
+4.509109132748808035e-10	1.632035470493057925e-09
+4.741898335527073965e-10	1.668962691745955538e-09
+4.982663818095965365e-10	1.705742204113841353e-09
+5.229780782720355362e-10	1.743851096446831689e-09
+5.484874027135371863e-10	1.781959988779821819e-09
+5.746466462490897931e-10	1.820068881112811948e-09
+6.014705797671944534e-10	1.859359444525895008e-09
+6.290773703758606673e-10	1.898797716823989451e-09
+6.573193091900766373e-10	1.938088280237072511e-09
+6.863588759833552579e-10	1.978708223615259885e-09
+7.160631327591860354e-10	2.019180458108435255e-09
+7.465354757370770627e-10	2.060982072566715352e-09
+7.776725086975202469e-10	2.102783687024995450e-09
+8.094594607520143878e-10	2.144585301483275134e-09
+8.420292698970699788e-10	2.187568587021647955e-09
+8.752489981361765263e-10	2.230551872560020776e-09
+9.091334163578353343e-10	2.273535158098393184e-09
+9.436677536735450989e-10	2.317847823601870319e-09
+9.789849480798162102e-10	2.362012780220335864e-09
+1.014833894472129088e-09	2.407507116803905309e-09
+1.051332759958492923e-09	2.453001453387475167e-09
+1.088481544538907922e-09	2.498495789971045026e-09
+1.126413186209884267e-09	2.545171797634707195e-09
+1.164876579866902172e-09	2.591847805298369777e-09
+1.204108059725980637e-09	2.639853192927136674e-09
+1.244137167564121442e-09	2.687710871670891566e-09
+1.284830965384814404e-09	2.735716259299658463e-09
+1.326425787404078317e-09	2.784903318008518083e-09
+1.368286485416362802e-09	2.833942667832365699e-09
+1.398374785293223653e-09	2.883129726541225319e-09
+1.429762923358186624e-09	2.933498456330177662e-09
+1.461165832311650588e-09	2.983867186119129592e-09
+1.493927663007221267e-09	3.034383624793093526e-09
+1.527900706559887071e-09	3.085934025662138180e-09
+1.564532510042761296e-09	3.137632135416194838e-09
+1.603675364570832765e-09	3.190364207365332629e-09
+1.646954067879228560e-09	3.244425659279574320e-09
+1.692743822232821393e-09	3.297305440113723701e-09
+1.743703387561820099e-09	3.352696271993070120e-09

lalsimulation/src/LALSimNeutronStarEOS_AP2.dat

+8.343631758107906939e-38	5.760646279117794535e-24
+8.343631580857245259e-37	5.804959210497264011e-24
+8.343631580857246261e-36	5.834500972728693026e-24
+8.343632038754789638e-35	6.019136798346299425e-24
+9.995903227680446623e-34	8.567115626089990214e-24
+1.156545764779204987e-32	1.211212891068081416e-23
+1.404386579748772230e-31	3.330835445636685513e-23
+4.807924266210808738e-30	1.565714166351965477e-22
+1.569554700758486418e-28	8.493260962020868104e-22
+8.049543353278943375e-27	7.710404093022768032e-21
+4.104091370047027957e-26	1.936463467731056542e-20
+2.008249911992246865e-25	4.864791994918719978e-20
+9.508463995767528335e-25	1.221552490862560189e-19
+4.350321992736018382e-24	3.069390689624394171e-19
+1.914897940977587221e-23	7.710403472645451774e-19
+8.058701481400324577e-23	1.936463526814610684e-18
+3.230836530483839324e-22	4.865530760907105231e-18
+4.344561346220565862e-22	6.124748887463802070e-18
+1.185467223548028312e-21	1.222290965864442148e-17
+3.166435353222566329e-21	2.438673691541350992e-17
+8.310692455186995649e-21	4.866269349644829091e-17
+2.151232156996131378e-20	9.711859588326032257e-17
+5.515154112228532168e-20	1.937940497497618577e-16
+7.218533151432849524e-20	2.440150869016797848e-16
+1.345731328335471147e-19	3.867756947236045186e-16
+2.502336294835788531e-19	6.130657183780712154e-16
+3.411041192947317021e-19	7.717789064604915480e-16
+4.160220775893210012e-19	9.719244456511960121e-16
+5.667146953719447417e-19	1.223768099027223581e-15
+1.050801028651695126e-18	1.939417719285730906e-15
+1.946359918714167429e-18	3.075298853003307365e-15
+3.603505884888307650e-18	4.875131897516413091e-15
+4.677349567547897744e-18	6.138781246310792288e-15
+6.362707709195922889e-18	7.725174671335268573e-15
+8.657661113518540592e-18	9.734015669972665520e-15
+1.177195562914600537e-17	1.225245236621271518e-14
+1.601016575098847113e-17	1.543557848371115445e-14
+2.067776666506359893e-17	1.943110352785689595e-14
+2.812081797163315919e-17	2.446797680216988522e-14
+3.823149159380315039e-17	3.081207326570878991e-14
+4.914569977792957113e-17	3.880312335400026662e-14
+6.682349839757224197e-17	4.886948667400894403e-14
+9.087197812587819867e-17	6.153552400687944153e-14
+1.235023519019284382e-16	7.747331269963000231e-14
+1.679450037352667011e-16	9.763557314036988077e-14
+2.145471643418675691e-16	1.228937920342251212e-13
+2.389043447093903019e-16	1.361875875494194413e-13
+2.717843425129702121e-16	1.547989144463240331e-13
+3.695233117251393371e-16	1.949757237840322556e-13
+4.804674567344333288e-16	2.455660272401238043e-13
+6.447788056504376238e-16	3.175002571949458142e-13
+6.517949865510213109e-16	3.293908150529345195e-13
+6.899629978881484835e-16	3.861110313286516761e-13
+7.515870871085186872e-16	4.881778664403938197e-13
+8.121477846155580033e-16	5.881767815932402550e-13
+8.946727017443121825e-16	7.184559915858633373e-13
+1.006192902542618365e-15	8.832991116900644424e-13
+1.155718643767072154e-15	1.086398862554044061e-12
+1.041020880690464571e-14	6.226173798963438979e-12
+2.143842591667660228e-14	1.246447412811411696e-11
+3.880280312113756540e-14	1.871067899484997050e-11
+5.992176115942451360e-14	2.496343966934196134e-11
+8.412503846984265287e-14	3.122191608684835906e-11
+1.111429062208144745e-13	3.748553942045298565e-11
+1.408637360812523970e-13	4.375390480010202386e-11
+1.732619623019208968e-13	5.002671562635437316e-11
+2.083647557845088832e-13	5.630375048359139163e-11
+2.462337999071061373e-13	6.258484260848191244e-11
+2.869544408634723667e-13	6.886986644847367973e-11
+3.306288577519030857e-13	7.515872761758915987e-11
+3.773716254003242311e-13	8.145135713577904988e-11
+4.273067380322587241e-13	8.774770537285798841e-11
+4.805655841780131941e-13	9.404773924203574316e-11
+5.364373251680912945e-13	1.003514389503217084e-10
+5.899197449592882027e-13	1.054759606090761956e-10
+6.516472880056316224e-13	1.117210922673662161e-10
+7.161222163131906909e-13	1.179677010145063618e-10
+7.838171983140023331e-13	1.242187410281968572e-10
+8.552196733286050129e-13	1.304712581307374520e-10
+9.308318515660379253e-13	1.367267294109782973e-10
+1.011096859681335826e-12	1.429851548689193671e-10
+1.096531678772039651e-12	1.492450574157105363e-10
+1.187608977270185922e-12	1.555079141402019559e-10
+1.284816194496313404e-12	1.617707708646933755e-10
+1.388625998882458822e-12	1.680483984776859438e-10
+1.499540600637661774e-12	1.743112552021773634e-10
+1.618003126416956914e-12	1.805888828151699317e-10
+1.744589640871889619e-12	1.868665104281625258e-10
+1.879595561772482990e-12	1.931441380411550941e-10
+2.023759433543794585e-12	1.994365365426488369e-10
+2.177376673955848314e-12	2.057141641556414052e-10
+2.341038118548689981e-12	2.120065626571351223e-10
+2.515039185092343092e-12	2.182989611586288651e-10
+2.700118418011865406e-12	2.245913596601226080e-10
+2.896571235077280027e-12	2.308837581616163250e-10
+3.104988471828633565e-12	2.371909275516112424e-10
+3.325960963805972226e-12	2.434980969416061339e-10
+3.559636419894307563e-12	2.498052663316010255e-10
+3.806753384518697737e-12	2.561124357215959170e-10
+4.067754984334177404e-12	2.624196051115908086e-10
+4.343084345995780815e-12	2.687415453900868489e-10
+4.633184596158543432e-12	2.750487147800817404e-10
+4.938646570362511463e-12	2.813706550585777807e-10
+5.259913395262720369e-12	2.877073662255750214e-10
+5.941442190706196711e-12	3.003660176710682506e-10
+6.676589311408880044e-12	3.131428362245708246e-10
+7.469342897266083879e-12	3.259196547780733468e-10
+8.320884619357899014e-12	3.387112442200770178e-10
+9.234020946499548969e-12	3.514880627735795918e-10
+1.021111522085121696e-11	3.643830484350914071e-10
+1.125364453126301849e-11	3.772928049851043711e-10
+1.236692639759537267e-11	3.902025615351173351e-10
+1.355199478204336034e-11	4.032304851931395405e-10
+1.480781572241189743e-11	4.161254708546513558e-10
+1.615492075371759769e-11	4.291533945126736129e-10
+1.756849478327851558e-11	4.422994852787051113e-10
+1.908546503234754550e-11	4.553274089367273683e-10
+2.066742719082167432e-11	4.684882705912600671e-10
+2.235573974650414619e-11	4.817525284653008069e-10
+2.413415472204369043e-11	4.950315572278427470e-10
+2.600562629514053485e-11	5.082958151018834867e-10
+2.798049408774549452e-11	5.215748438644254269e-10
+3.004694138905763562e-11	5.349720397349766601e-10
+3.221826199872801395e-11	5.483692356055278933e-10
+3.449445591675661658e-11	5.618845985840884195e-10
+3.688881694279448641e-11	5.753999615626488423e-10
+3.938657418834047473e-11	5.889153245412093685e-10
+4.199068183109480610e-11	6.025636255162802848e-10
+4.472477329265933526e-11	6.163300935993604941e-10
+4.756226097373198290e-11	6.300965616824406000e-10
+5.051939285166401648e-11	6.438630297655208093e-10
+5.359321474875519851e-11	6.577476649566103117e-10
+5.679849755350670353e-11	6.716323001476997106e-10
+6.012047037741736992e-11	6.856351024467985060e-10
+6.357390410898833346e-11	6.996526756343983984e-10
+6.715879874821962645e-11	7.137884159300075839e-10
+7.086038340661007435e-11	7.280423233336259589e-10
+7.470672277231187521e-11	7.421780636292351444e-10
+7.868304595682388032e-11	7.565501381408628125e-10
+8.280264675979712609e-11	7.709369835409916810e-10
+8.705518555928080069e-11	7.854272251606286422e-10
+9.144952488837560367e-11	7.999322376687668037e-10
+9.597532512513072318e-11	8.145554172849142584e-10
+1.006591738688482450e-10	8.291785969010616096e-10
+1.054863002310270075e-10	8.439347145137194542e-10
+1.104699980113180470e-10	8.588089992343865919e-10
+1.155969734100703402e-10	8.736685130665525292e-10
+1.208687035161339734e-10	8.886609648952289599e-10
+1.263117759288110712e-10	9.036681876124064876e-10
+1.318567674721461622e-10	9.187788065490921080e-10
+1.376129827210478400e-10	9.340223634822883252e-10
+1.436306426964199998e-10	9.493840875234935253e-10
+1.496438714052418357e-10	9.647458115646989321e-10
+1.623320646277374140e-10	9.957055938631281251e-10
+1.756849478327851494e-10	1.027167586370596891e-09
+1.896729792433827185e-10	1.058880683982585237e-09
+2.043257006365324447e-10	1.090963053807102559e-09
+2.198942171167540368e-10	1.123399924955647863e-09
+2.359797146945161697e-10	1.156339006313232327e-09
+2.529957782478513689e-10	1.189514421886835376e-09
+2.707799280032468178e-10	1.223192047669477999e-09
+2.892287677411943720e-10	1.257238945664649391e-09
+3.085786316777126693e-10	1.291640344983848766e-09
+3.287261235932935654e-10	1.327282498937145459e-09
+3.497894105959463016e-10	1.361698669144846034e-09
+3.716355546891604362e-10	1.397340823098142520e-09
+3.943974938694464626e-10	1.434208960797035538e-09
+4.180752281368043291e-10	1.471077098495928349e-09
+4.426687574912340874e-10	1.507960007083322360e-09
+4.681780819327356858e-10	1.546068899416312489e-09
+4.946032014613091761e-10	1.584177791749302826e-09
+5.219293451884533060e-10	1.622286684082292955e-09
+5.503189928876808665e-10	1.661577247495376015e-09
+5.797573736704907606e-10	1.700867810908458868e-09
+6.100967786518713978e-10	1.741487754286646242e-09
+6.414996876053353621e-10	1.781959988779821819e-09
+6.739365587538805630e-10	1.823761603238101709e-09
+7.075551009825184875e-10	1.865563217696381600e-09
+7.422223762947386422e-10	1.908546503234754421e-09
+7.779383846905411304e-10	1.951529788773127242e-09
+8.148212932779351418e-10	1.994660783196511240e-09
+8.528858729454218769e-10	2.038825739814976785e-09
+8.919991856964908421e-10	2.083138405318453920e-09
+9.324271075241630243e-10	2.128632741902023365e-09
+9.738889915469163397e-10	2.175308749565685948e-09
+1.016665484646272769e-09	2.221984757229348530e-09
+1.060741815933731318e-09	2.268660764893010699e-09
+1.105999818301282384e-09	2.316666152521777595e-09
+1.152572429745436770e-09	2.364523831265532901e-09
+1.200326712269683880e-09	2.413710889974392108e-09
+1.249513770978543293e-09	2.464079619763344451e-09
+1.300000667875504826e-09	2.514448349552296795e-09
+1.351270421863027706e-09	2.566146459306353039e-09
+1.404903518010736034e-09	2.617696860175398106e-09
+1.459851223235047474e-09	2.670576641009547488e-09
+1.516083995758959835e-09	2.723456421843696455e-09
+1.572361080948375590e-09	2.777517873757938560e-09
+1.631149217182988382e-09	2.832760996752272975e-09
+1.692743822232821393e-09	2.888004119746607804e-09
+1.754190718397642814e-09	2.944576622706046946e-09