star-schiff

Library for estimating radiative properties
git clone git://git.meso-star.com/star-schiff.git
Log | Files | Refs | README | LICENSE

commit 41ae55cab1343cb6d364f18dbebe4dc5ae11234c
parent e57913c49dac1b95bacf71cb3c3d71d2723a6471
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Thu,  5 Nov 2015 15:51:41 +0100

Push further the test with cylinder distributions

Compare the schiff estimation with the new results sent by Julien. Test
the extinction, absorption, scattering and average projected area
values.

Diffstat:
Acylindres_schiff_part.txt | 57+++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Aoptical_properties.txt | 552+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Msrc/test_sschiff_estimator_cylinder.c | 1031++++++++++++++++++++++++++++++-------------------------------------------------
3 files changed, 998 insertions(+), 642 deletions(-)

diff --git a/cylindres_schiff_part.txt b/cylindres_schiff_part.txt @@ -0,0 +1,57 @@ +#Lambda (m) section_efficace_absorption (m2/MO) section_efficace_diffusion (m2/MO) section_efficace_extinction (m2/MO) +4.00000000000000e-07 2.81148276917168e-13 2.48636076883100e-12 2.76750904574817e-12 +4.10000000000000e-07 2.24009382471161e-13 2.51271968811763e-12 2.73672907058880e-12 +4.20000000000000e-07 1.61344711105939e-13 2.50503191756686e-12 2.66637662867280e-12 +4.30000000000000e-07 1.39243171825275e-13 2.45153550914157e-12 2.59077868096684e-12 +4.40000000000000e-07 1.41848754459858e-13 2.40013784024617e-12 2.54198659470604e-12 +4.50000000000000e-07 1.56193174385628e-13 2.34480377552262e-12 2.50099694990826e-12 +4.60000000000000e-07 1.68136059476062e-13 2.30128761581120e-12 2.46942367528726e-12 +4.70000000000000e-07 1.78670999686088e-13 2.26691165240638e-12 2.44558265209247e-12 +4.80000000000000e-07 1.79662384192516e-13 2.23145247595974e-12 2.41111486015226e-12 +4.90000000000000e-07 1.78635587341230e-13 2.19404479608476e-12 2.37268038342598e-12 +5.00000000000000e-07 1.79927880782966e-13 2.14789647002885e-12 2.32782435081181e-12 +5.10000000000000e-07 2.16339389824781e-13 2.09844594620772e-12 2.31478533603249e-12 +5.20000000000000e-07 2.19354122273430e-13 2.08913310933374e-12 2.30848723160716e-12 +5.30000000000000e-07 1.91474811285185e-13 2.07121119378590e-12 2.26268600507108e-12 +5.40000000000000e-07 1.91816899012804e-13 2.02456238073205e-12 2.21637927974485e-12 +5.50000000000000e-07 2.13209672239836e-13 1.99701046675190e-12 2.21022013899174e-12 +5.60000000000000e-07 1.94418370736651e-13 1.99544678641674e-12 2.18986515715339e-12 +5.70000000000000e-07 1.53288768186043e-13 1.98477219066945e-12 2.13806095885549e-12 +5.80000000000000e-07 1.13058916188841e-13 1.95669352230221e-12 2.06975243849105e-12 +5.90000000000000e-07 9.75058786937038e-14 1.92104033377755e-12 2.01854621247124e-12 +6.00000000000000e-07 7.49967618855186e-14 1.89161985016549e-12 1.96661661205101e-12 +6.10000000000000e-07 5.09040827759909e-14 1.85618232050112e-12 1.90708640327711e-12 +6.20000000000000e-07 3.60890970275237e-14 1.81435711154512e-12 1.85044620857265e-12 +6.30000000000000e-07 2.49641915615705e-14 1.77226056295584e-12 1.79722475451741e-12 +6.40000000000000e-07 2.21770410805031e-14 1.72931103497851e-12 1.75148807605901e-12 +6.50000000000000e-07 2.23815127740530e-14 1.68911545088789e-12 1.71149696366194e-12 +6.60000000000000e-07 2.26417310060833e-14 1.65123794910088e-12 1.67387968010696e-12 +6.70000000000000e-07 2.42954918879993e-14 1.61489226328345e-12 1.63918775517145e-12 +6.80000000000000e-07 2.50199103661435e-14 1.58047335712953e-12 1.60549326749567e-12 +6.90000000000000e-07 2.61526536235239e-14 1.54734484662905e-12 1.57349750025259e-12 +7.00000000000000e-07 2.68952232034447e-14 1.51539949791961e-12 1.54229472112306e-12 +7.10000000000000e-07 2.82870868717291e-14 1.48426511567564e-12 1.51255220254738e-12 +7.20000000000000e-07 2.87138063804889e-14 1.45404603873700e-12 1.48275984511749e-12 +7.30000000000000e-07 2.95114369492585e-14 1.42423598716962e-12 1.45374742411889e-12 +7.40000000000000e-07 3.04942952611187e-14 1.39449735598055e-12 1.42499165124167e-12 +7.50000000000000e-07 3.40894078115008e-14 1.36468712183250e-12 1.39877652964399e-12 +7.60000000000000e-07 3.71254908035759e-14 1.33633922603518e-12 1.37346471683876e-12 +7.70000000000000e-07 3.95486268987460e-14 1.30412696169147e-12 1.34367558859022e-12 +7.80000000000000e-07 6.10373124310176e-14 1.27174247249898e-12 1.33277978493000e-12 +7.90000000000000e-07 7.53698495261906e-14 1.24768526960746e-12 1.32305511913364e-12 +8.00000000000000e-07 8.54123891952407e-14 1.22708141229494e-12 1.31249380149019e-12 +8.10000000000000e-07 8.71620723663670e-14 1.21147153754423e-12 1.29863360991060e-12 +8.20000000000000e-07 7.86416664314719e-14 1.18852783696394e-12 1.26716950339541e-12 +8.30000000000000e-07 8.34793448938697e-14 1.15587575564595e-12 1.23935510053982e-12 +8.40000000000000e-07 1.18350471704465e-13 1.11902174441848e-12 1.23737221612295e-12 +8.50000000000000e-07 1.57804794271177e-13 1.09286032549472e-12 1.25066511976589e-12 +8.60000000000000e-07 2.16325393577261e-13 1.07084592343472e-12 1.28717131701198e-12 +8.70000000000000e-07 2.74648117116301e-13 1.07025425708770e-12 1.34490237420400e-12 +8.80000000000000e-07 3.03581061880574e-13 1.09625202110578e-12 1.39983308298635e-12 +8.90000000000000e-07 2.78148862947325e-13 1.13934695275331e-12 1.41749581570064e-12 +9.00000000000000e-07 1.81456706704593e-13 1.17697051457863e-12 1.35842722128322e-12 +9.10000000000000e-07 8.07784402161031e-14 1.17682749489327e-12 1.25760593510938e-12 +9.20000000000000e-07 2.38496161775180e-14 1.14791980686658e-12 1.17176942304411e-12 +9.30000000000000e-07 7.01062772989699e-15 1.11390219834299e-12 1.12091282607288e-12 +9.40000000000000e-07 2.88084036248655e-15 1.08557878254289e-12 1.08845962290538e-12 +9.50000000000000e-07 1.98318674610597e-15 1.06177669883265e-12 1.06375988557876e-12 diff --git a/optical_properties.txt b/optical_properties.txt @@ -0,0 +1,552 @@ +lambda_vacuum (m) / n_r / kappa_r / n_medium + 4.000000000000000E-07 1.035286987991962E+00 1.642464685700659E-03 1.341779843605380E+00 + 4.010000000000000E-07 1.035394484151721E+00 1.624259338990003E-03 1.341687493674650E+00 + 4.020000000000000E-07 1.035498485175403E+00 1.594876393429627E-03 1.341595857289120E+00 + 4.030000000000000E-07 1.035593745454516E+00 1.563339730364777E-03 1.341504925814500E+00 + 4.040000000000000E-07 1.035689619601582E+00 1.544600372647700E-03 1.341414690748000E+00 + 4.050000000000000E-07 1.035787642729216E+00 1.511779431504190E-03 1.341325143715940E+00 + 4.060000000000000E-07 1.035873513905482E+00 1.471484460290955E-03 1.341236276471320E+00 + 4.070000000000000E-07 1.035966600915103E+00 1.462298106531343E-03 1.341148080891490E+00 + 4.080000000000000E-07 1.036060494953379E+00 1.407755470265253E-03 1.341060548975860E+00 + 4.090000000000000E-07 1.036137924348892E+00 1.368311430748199E-03 1.340973672843670E+00 + 4.100000000000000E-07 1.036210685221302E+00 1.330046390562492E-03 1.340887444731820E+00 + 4.110000000000000E-07 1.036280450985456E+00 1.300974328047466E-03 1.340801856992660E+00 + 4.120000000000000E-07 1.036362608081751E+00 1.286212700588483E-03 1.340716902091980E+00 + 4.130000000000000E-07 1.036447639684917E+00 1.236376528960547E-03 1.340632572606910E+00 + 4.140000000000000E-07 1.036513150647102E+00 1.179708889220116E-03 1.340548861223910E+00 + 4.150000000000000E-07 1.036573630808206E+00 1.146669303767855E-03 1.340465760736840E+00 + 4.160000000000000E-07 1.036625512534661E+00 1.088836702310341E-03 1.340383264045010E+00 + 4.170000000000000E-07 1.036667015511619E+00 1.062695995814180E-03 1.340301364151300E+00 + 4.180000000000000E-07 1.036722768320762E+00 1.052568689767748E-03 1.340220054160350E+00 + 4.190000000000000E-07 1.036779280227207E+00 1.006697219981148E-03 1.340139327276700E+00 + 4.200000000000000E-07 1.036822939800069E+00 9.722694769944399E-04 1.340059176803090E+00 + 4.210000000000000E-07 1.036865833600609E+00 9.469232619753497E-04 1.339979596138690E+00 + 4.220000000000000E-07 1.036897653219950E+00 9.040007897511271E-04 1.339900578777420E+00 + 4.230000000000000E-07 1.036927468207731E+00 9.031823993726944E-04 1.339822118306300E+00 + 4.240000000000000E-07 1.036971750881576E+00 8.924646106519200E-04 1.339744208403820E+00 + 4.250000000000000E-07 1.037006936112468E+00 8.594846944392758E-04 1.339666842838340E+00 + 4.260000000000000E-07 1.037033040892262E+00 8.523605328290717E-04 1.339590015466570E+00 + 4.270000000000000E-07 1.037064992312203E+00 8.534479686759623E-04 1.339513720232010E+00 + 4.280000000000000E-07 1.037100162121758E+00 8.506283100332938E-04 1.339437951163470E+00 + 4.290000000000000E-07 1.037135678599242E+00 8.522590014316155E-04 1.339362702373630E+00 + 4.300000000000000E-07 1.037173395465157E+00 8.565016952270480E-04 1.339287968057550E+00 + 4.310000000000000E-07 1.037214737528130E+00 8.639171883406666E-04 1.339213742491340E+00 + 4.320000000000000E-07 1.037259758643682E+00 8.691165576390881E-04 1.339140020030710E+00 + 4.330000000000000E-07 1.037303754455989E+00 8.678624596520924E-04 1.339066795109690E+00 + 4.340000000000000E-07 1.037346655520205E+00 8.725074419904887E-04 1.338994062239260E+00 + 4.350000000000000E-07 1.037390939381917E+00 8.756139553518682E-04 1.338921816006080E+00 + 4.360000000000000E-07 1.037435399005391E+00 8.788395944917436E-04 1.338850051071230E+00 + 4.370000000000000E-07 1.037478748753387E+00 8.799399526191331E-04 1.338778762168930E+00 + 4.380000000000000E-07 1.037521676917342E+00 8.850314449839643E-04 1.338707944105350E+00 + 4.390000000000000E-07 1.037565219615768E+00 8.888250046471429E-04 1.338637591757410E+00 + 4.400000000000000E-07 1.037608696500256E+00 8.936647271660125E-04 1.338567700071630E+00 + 4.410000000000000E-07 1.037652289517359E+00 8.980316117027969E-04 1.338498264062920E+00 + 4.420000000000000E-07 1.037695839416368E+00 9.024090285591408E-04 1.338429278813520E+00 + 4.430000000000000E-07 1.037737726850386E+00 9.038232640135763E-04 1.338360739471860E+00 + 4.440000000000000E-07 1.037774845995166E+00 9.041864444660936E-04 1.338292641251510E+00 + 4.450000000000000E-07 1.037810077590406E+00 9.145542051765026E-04 1.338224979430060E+00 + 4.460000000000000E-07 1.037846610963693E+00 9.257787115899892E-04 1.338157749348150E+00 + 4.470000000000000E-07 1.037878259630217E+00 9.335983043538121E-04 1.338090946408400E+00 + 4.480000000000000E-07 1.037914617228649E+00 9.702316913975443E-04 1.338024566074410E+00 + 4.490000000000000E-07 1.037970208539347E+00 1.007116632500707E-03 1.337958603869820E+00 + 4.500000000000000E-07 1.038027021857111E+00 1.009173031428449E-03 1.337893055377290E+00 + 4.510000000000000E-07 1.038079358354862E+00 1.027814707250541E-03 1.337827916237600E+00 + 4.520000000000000E-07 1.038129195790007E+00 1.028696923944767E-03 1.337763182148700E+00 + 4.530000000000000E-07 1.038171798198206E+00 1.042122191740612E-03 1.337698848864820E+00 + 4.540000000000000E-07 1.038225010269563E+00 1.074869388950115E-03 1.337634912195530E+00 + 4.550000000000000E-07 1.038283232720702E+00 1.078816153504866E-03 1.337571368004930E+00 + 4.560000000000000E-07 1.038338570451783E+00 1.092245635394565E-03 1.337508212210750E+00 + 4.570000000000000E-07 1.038393288814794E+00 1.094812200306328E-03 1.337445440783500E+00 + 4.580000000000000E-07 1.038440433963812E+00 1.093485331340192E-03 1.337383049745670E+00 + 4.590000000000000E-07 1.038482662253094E+00 1.103439733475396E-03 1.337321035170910E+00 + 4.600000000000000E-07 1.038522598315466E+00 1.113089490368987E-03 1.337259393183210E+00 + 4.610000000000000E-07 1.038570039357183E+00 1.148516858552576E-03 1.337198119956160E+00 + 4.620000000000000E-07 1.038624907459639E+00 1.162434238059191E-03 1.337137211712120E+00 + 4.630000000000000E-07 1.038687104626991E+00 1.195966780613389E-03 1.337076664721550E+00 + 4.640000000000000E-07 1.038755514562735E+00 1.199699923934052E-03 1.337016475302180E+00 + 4.650000000000000E-07 1.038817135628012E+00 1.200459607388018E-03 1.336956639818350E+00 + 4.660000000000000E-07 1.038874078123107E+00 1.201186493947777E-03 1.336897154680260E+00 + 4.670000000000000E-07 1.038928224588712E+00 1.201880424343069E-03 1.336838016343290E+00 + 4.680000000000000E-07 1.038980178159473E+00 1.204040146415200E-03 1.336779221307310E+00 + 4.690000000000000E-07 1.039031437175215E+00 1.207552177169499E-03 1.336720766115970E+00 + 4.700000000000000E-07 1.039082221078632E+00 1.211092670367667E-03 1.336662647356080E+00 + 4.710000000000000E-07 1.039132490661054E+00 1.215106981436297E-03 1.336604861656970E+00 + 4.720000000000000E-07 1.039182661408432E+00 1.219870817381094E-03 1.336547405689770E+00 + 4.730000000000000E-07 1.039232954048947E+00 1.224643796577837E-03 1.336490276166880E+00 + 4.740000000000000E-07 1.039283376885394E+00 1.229326344827146E-03 1.336433469841280E+00 + 4.750000000000000E-07 1.039334064443123E+00 1.234117398218758E-03 1.336376983505970E+00 + 4.760000000000000E-07 1.039385526833898E+00 1.238492582229907E-03 1.336320813993330E+00 + 4.770000000000000E-07 1.039436291084956E+00 1.239066817439063E-03 1.336264958174580E+00 + 4.780000000000000E-07 1.039484948734477E+00 1.239657200984551E-03 1.336209412959190E+00 + 4.790000000000000E-07 1.039531959467774E+00 1.240540672112777E-03 1.336154175294280E+00 + 4.800000000000000E-07 1.039578615206627E+00 1.244442394725602E-03 1.336099242164120E+00 + 4.810000000000000E-07 1.039626229298974E+00 1.248349409793080E-03 1.336044610589550E+00 + 4.820000000000000E-07 1.039673858625137E+00 1.250235937286709E-03 1.335990277627470E+00 + 4.830000000000000E-07 1.039720595019222E+00 1.251865430581906E-03 1.335936240370270E+00 + 4.840000000000000E-07 1.039766629301611E+00 1.253490577728770E-03 1.335882495945360E+00 + 4.850000000000000E-07 1.039811982223968E+00 1.255111379800570E-03 1.335829041514630E+00 + 4.860000000000000E-07 1.039856756652771E+00 1.256727837858235E-03 1.335775874273970E+00 + 4.870000000000000E-07 1.039900966005627E+00 1.258467903243629E-03 1.335722991452750E+00 + 4.880000000000000E-07 1.039944611565546E+00 1.260075944188013E-03 1.335670390313400E+00 + 4.890000000000000E-07 1.039987632682271E+00 1.261679644222744E-03 1.335618068150840E+00 + 4.900000000000000E-07 1.040030012263802E+00 1.263407757979589E-03 1.335566022292100E+00 + 4.910000000000000E-07 1.040071688175533E+00 1.264745004221800E-03 1.335514250095810E+00 + 4.920000000000000E-07 1.040112378686561E+00 1.266076841520105E-03 1.335462748951790E+00 + 4.930000000000000E-07 1.040151955742799E+00 1.267403270885834E-03 1.335411516280550E+00 + 4.940000000000000E-07 1.040190287896555E+00 1.268724293318496E-03 1.335360549532910E+00 + 4.950000000000000E-07 1.040226903080167E+00 1.270352131527932E-03 1.335309846189550E+00 + 4.960000000000000E-07 1.040262146697474E+00 1.273305523147441E-03 1.335259403760560E+00 + 4.970000000000000E-07 1.040295747313031E+00 1.276077378609423E-03 1.335209219785090E+00 + 4.980000000000000E-07 1.040325648091718E+00 1.278823648922353E-03 1.335159291830880E+00 + 4.990000000000000E-07 1.040353759919826E+00 1.289072947859305E-03 1.335109617493900E+00 + 5.000000000000000E-07 1.040379481152520E+00 1.299274503512561E-03 1.335060194397920E+00 + 5.010000000000001E-07 1.040408126491521E+00 1.325811156480386E-03 1.335011020194160E+00 + 5.020000000000000E-07 1.040441036092354E+00 1.344764441067616E-03 1.334962092560890E+00 + 5.030000000000000E-07 1.040477445253276E+00 1.374045011533224E-03 1.334913409203050E+00 + 5.040000000000000E-07 1.040512648034075E+00 1.387903827917666E-03 1.334864967851890E+00 + 5.050000000000000E-07 1.040545900251519E+00 1.418801858069103E-03 1.334816766264610E+00 + 5.060000000000000E-07 1.040590906700738E+00 1.468464949020522E-03 1.334768802224000E+00 + 5.070000000000000E-07 1.040639195455336E+00 1.482563959735157E-03 1.334721073538070E+00 + 5.080000000000000E-07 1.040686458775797E+00 1.525559361876046E-03 1.334673578039750E+00 + 5.090000000000000E-07 1.040741671425000E+00 1.559105378939951E-03 1.334626313586500E+00 + 5.100000000000000E-07 1.040807713567804E+00 1.603233722101158E-03 1.334579278060000E+00 + 5.110000000000000E-07 1.040875584410848E+00 1.613985571121084E-03 1.334532469365830E+00 + 5.120000000000000E-07 1.040943220348469E+00 1.649380251826485E-03 1.334485885433120E+00 + 5.130000000000000E-07 1.041018230706438E+00 1.669091265929868E-03 1.334439524214260E+00 + 5.140000000000000E-07 1.041092713713321E+00 1.676566182400091E-03 1.334393383684570E+00 + 5.150000000000000E-07 1.041165314107454E+00 1.687389518793340E-03 1.334347461841990E+00 + 5.160000000000000E-07 1.041238597163079E+00 1.692922480695940E-03 1.334301756706790E+00 + 5.170000000000000E-07 1.041311956633380E+00 1.695363690601763E-03 1.334256266321260E+00 + 5.180000000000000E-07 1.041389428013160E+00 1.700417308885367E-03 1.334210988749400E+00 + 5.190000000000000E-07 1.041464755955924E+00 1.682652012825744E-03 1.334165922076680E+00 + 5.200000000000000E-07 1.041530069541219E+00 1.658812761280325E-03 1.334121064409700E+00 + 5.210000000000000E-07 1.041594181696574E+00 1.657124535655681E-03 1.334076413875950E+00 + 5.220000000000000E-07 1.041657845088931E+00 1.630122704418439E-03 1.334031968623510E+00 + 5.230000000000000E-07 1.041711910952757E+00 1.607543355228116E-03 1.333987726820790E+00 + 5.240000000000000E-07 1.041759947842067E+00 1.582244525828878E-03 1.333943686656260E+00 + 5.249999999999999E-07 1.041810066348116E+00 1.583879702931137E-03 1.333899846338170E+00 + 5.260000000000000E-07 1.041862095522510E+00 1.557720519804618E-03 1.333856204094340E+00 + 5.270000000000000E-07 1.041913531175414E+00 1.546395579258294E-03 1.333812758171830E+00 + 5.280000000000000E-07 1.041954667678475E+00 1.499811551277397E-03 1.333769506836780E+00 + 5.290000000000000E-07 1.041989797310988E+00 1.498636177737314E-03 1.333726448374050E+00 + 5.300000000000000E-07 1.042024857461786E+00 1.469804039124607E-03 1.333683581087100E+00 + 5.310000000000000E-07 1.042050856148998E+00 1.445880383759585E-03 1.333640903297660E+00 + 5.320000000000001E-07 1.042077186838255E+00 1.447977496199134E-03 1.333598413345510E+00 + 5.330000000000000E-07 1.042102162421808E+00 1.427209064601831E-03 1.333556109588290E+00 + 5.340000000000000E-07 1.042120822250241E+00 1.426699703280908E-03 1.333513990401220E+00 + 5.350000000000000E-07 1.042141900864721E+00 1.433471379762237E-03 1.333472054176910E+00 + 5.360000000000000E-07 1.042166620404672E+00 1.444941845941278E-03 1.333430299325120E+00 + 5.370000000000000E-07 1.042194724365646E+00 1.457717474653508E-03 1.333388724272540E+00 + 5.380000000000000E-07 1.042225414187512E+00 1.471493967735956E-03 1.333347327462570E+00 + 5.390000000000000E-07 1.042256744777569E+00 1.484092000868983E-03 1.333306107355130E+00 + 5.400000000000000E-07 1.042288345839686E+00 1.500763725452133E-03 1.333265062426450E+00 + 5.410000000000000E-07 1.042324503859437E+00 1.528281898750163E-03 1.333224191168830E+00 + 5.420000000000000E-07 1.042364912605382E+00 1.549843563702498E-03 1.333183492090460E+00 + 5.430000000000000E-07 1.042406562887524E+00 1.568902404377878E-03 1.333142963715240E+00 + 5.440000000000000E-07 1.042454714949787E+00 1.608815343234295E-03 1.333102604582520E+00 + 5.450000000000000E-07 1.042513875213928E+00 1.637618035281262E-03 1.333062413247000E+00 + 5.460000000000000E-07 1.042576056800957E+00 1.653580138686441E-03 1.333022388278450E+00 + 5.470000000000000E-07 1.042640409898161E+00 1.675120248233485E-03 1.332982528261570E+00 + 5.480000000000000E-07 1.042707136014031E+00 1.685157909842487E-03 1.332942831795790E+00 + 5.490000000000000E-07 1.042774026206271E+00 1.695250147639022E-03 1.332903297495090E+00 + 5.500000000000000E-07 1.042842766321411E+00 1.705310235117845E-03 1.332863923987830E+00 + 5.510000000000000E-07 1.042912300283037E+00 1.705530449406497E-03 1.332824709916540E+00 + 5.520000000000000E-07 1.042981728226964E+00 1.708356091134469E-03 1.332785653937790E+00 + 5.530000000000000E-07 1.043051509491291E+00 1.704016278805774E-03 1.332746754721980E+00 + 5.540000000000000E-07 1.043124661120986E+00 1.708897901535750E-03 1.332708010953190E+00 + 5.550000000000000E-07 1.043199194602040E+00 1.692798767844739E-03 1.332669421329010E+00 + 5.559999999999999E-07 1.043265008027065E+00 1.670274669430007E-03 1.332630984560360E+00 + 5.570000000000000E-07 1.043336588355702E+00 1.676935291139613E-03 1.332592699371340E+00 + 5.580000000000000E-07 1.043412220918312E+00 1.644255097722219E-03 1.332554564499090E+00 + 5.590000000000000E-07 1.043477770799564E+00 1.615743588209385E-03 1.332516578693600E+00 + 5.600000000000000E-07 1.043536695477078E+00 1.579035215844224E-03 1.332478740717540E+00 + 5.610000000000000E-07 1.043600994183831E+00 1.581811014957712E-03 1.332441049346180E+00 + 5.620000000000000E-07 1.043666231315476E+00 1.535235605622221E-03 1.332403503367140E+00 + 5.630000000000001E-07 1.043718864830297E+00 1.473563431577708E-03 1.332366101580320E+00 + 5.640000000000000E-07 1.043771791283996E+00 1.484273905192923E-03 1.332328842797720E+00 + 5.650000000000000E-07 1.043825841694447E+00 1.417653789592510E-03 1.332291725843300E+00 + 5.660000000000000E-07 1.043867370062960E+00 1.377114388401143E-03 1.332254749552820E+00 + 5.670000000000000E-07 1.043898295606015E+00 1.324202628263995E-03 1.332217912773740E+00 + 5.680000000000000E-07 1.043933664352148E+00 1.332080947209094E-03 1.332181214365030E+00 + 5.690000000000000E-07 1.043977438602228E+00 1.304393676650687E-03 1.332144653197070E+00 + 5.700000000000000E-07 1.044018794045142E+00 1.259390337125543E-03 1.332108228151510E+00 + 5.710000000000000E-07 1.044062030345026E+00 1.259437872346123E-03 1.332071938121120E+00 + 5.720000000000000E-07 1.044109863822610E+00 1.216136708576497E-03 1.332035782009670E+00 + 5.730000000000000E-07 1.044157258612584E+00 1.187015609995578E-03 1.331999758731800E+00 + 5.740000000000000E-07 1.044202028283407E+00 1.143113983172242E-03 1.331963867212890E+00 + 5.750000000000000E-07 1.044232131241066E+00 1.063840725432586E-03 1.331928106388940E+00 + 5.760000000000000E-07 1.044246895319073E+00 1.023527710863081E-03 1.331892475206440E+00 + 5.770000000000000E-07 1.044252610112750E+00 9.697779794138100E-04 1.331856972622240E+00 + 5.780000000000000E-07 1.044253976775146E+00 9.497659966550409E-04 1.331821597603440E+00 + 5.790000000000000E-07 1.044259957725371E+00 9.352618179707669E-04 1.331786349127290E+00 + 5.800000000000000E-07 1.044274809362838E+00 9.395000916553773E-04 1.331751226181010E+00 + 5.810000000000000E-07 1.044292271809736E+00 9.117482360018617E-04 1.331716227761740E+00 + 5.820000000000000E-07 1.044308146494882E+00 9.039549474109366E-04 1.331681352876410E+00 + 5.830000000000000E-07 1.044328059952875E+00 9.029813184193805E-04 1.331646600541600E+00 + 5.840000000000000E-07 1.044352735493407E+00 8.971352837393257E-04 1.331611969783430E+00 + 5.850000000000000E-07 1.044378685632916E+00 8.831537746544346E-04 1.331577459637510E+00 + 5.860000000000000E-07 1.044405836424728E+00 8.786320819678296E-04 1.331543069148750E+00 + 5.869999999999999E-07 1.044435280426248E+00 8.692005288892336E-04 1.331508797371310E+00 + 5.880000000000000E-07 1.044465717884262E+00 8.558876109011710E-04 1.331474643368460E+00 + 5.890000000000000E-07 1.044495153084782E+00 8.403500870864586E-04 1.331440606212510E+00 + 5.900000000000000E-07 1.044523085786671E+00 8.224523096408955E-04 1.331406684984670E+00 + 5.910000000000000E-07 1.044550545262079E+00 8.104699876011172E-04 1.331372878774990E+00 + 5.920000000000000E-07 1.044578615724940E+00 7.960072814810816E-04 1.331339186682210E+00 + 5.930000000000000E-07 1.044607054549288E+00 7.816152660959925E-04 1.331305607813710E+00 + 5.940000000000001E-07 1.044635270377074E+00 7.646952835940411E-04 1.331272141285380E+00 + 5.950000000000000E-07 1.044663935653333E+00 7.513518759906050E-04 1.331238786221540E+00 + 5.960000000000000E-07 1.044693946287332E+00 7.357902455879724E-04 1.331205541754850E+00 + 5.970000000000000E-07 1.044724372554313E+00 7.167716979679335E-04 1.331172407026180E+00 + 5.980000000000000E-07 1.044753615647646E+00 6.946257364021612E-04 1.331139381184580E+00 + 5.990000000000000E-07 1.044780691526674E+00 6.712950232982840E-04 1.331106463387120E+00 + 6.000000000000000E-07 1.044803444015409E+00 6.412373047096696E-04 1.331073652798870E+00 + 6.010000000000000E-07 1.044824539307607E+00 6.249866054301384E-04 1.331040948592740E+00 + 6.020000000000000E-07 1.044845766620673E+00 6.034077108891513E-04 1.331008349949450E+00 + 6.030000000000000E-07 1.044866538700980E+00 5.784136786976474E-04 1.330975856057410E+00 + 6.040000000000000E-07 1.044887221935235E+00 5.676336347762863E-04 1.330943466112670E+00 + 6.050000000000000E-07 1.044906979832908E+00 5.365222304755139E-04 1.330911179318770E+00 + 6.060000000000000E-07 1.044925229429125E+00 5.166167506612064E-04 1.330878994886740E+00 + 6.070000000000000E-07 1.044942099825731E+00 4.987864403064095E-04 1.330846912034950E+00 + 6.080000000000000E-07 1.044955182314457E+00 4.646041110856508E-04 1.330814929989080E+00 + 6.090000000000000E-07 1.044965016968371E+00 4.527564487138666E-04 1.330783047981990E+00 + 6.100000000000000E-07 1.044979238163159E+00 4.409632490707041E-04 1.330751265253680E+00 + 6.110000000000000E-07 1.044994223279349E+00 4.219740270792474E-04 1.330719581051210E+00 + 6.120000000000000E-07 1.045007044211732E+00 4.020154703506148E-04 1.330687994628590E+00 + 6.130000000000000E-07 1.045017963621394E+00 3.842194369970159E-04 1.330656505246750E+00 + 6.140000000000000E-07 1.045027441486893E+00 3.679492624268454E-04 1.330625112173430E+00 + 6.150000000000000E-07 1.045035617615133E+00 3.519800711874577E-04 1.330593814683110E+00 + 6.160000000000000E-07 1.045043007111847E+00 3.395630359374884E-04 1.330562612056960E+00 + 6.170000000000000E-07 1.045051084952903E+00 3.326688305061111E-04 1.330531503582750E+00 + 6.179999999999999E-07 1.045061222694419E+00 3.259457256833445E-04 1.330500488554770E+00 + 6.190000000000000E-07 1.045073433158974E+00 3.187415925020136E-04 1.330469566273780E+00 + 6.200000000000000E-07 1.045088360723040E+00 3.171051848434899E-04 1.330438736046940E+00 + 6.210000000000000E-07 1.045103671053048E+00 3.017345937277039E-04 1.330407997187720E+00 + 6.220000000000000E-07 1.045116656866957E+00 2.864761737697034E-04 1.330377349015850E+00 + 6.230000000000000E-07 1.045126365845574E+00 2.711663840968767E-04 1.330346790857260E+00 + 6.240000000000000E-07 1.045133870688024E+00 2.592620517622472E-04 1.330316322043970E+00 + 6.250000000000001E-07 1.045140983506716E+00 2.527585737208309E-04 1.330285941914100E+00 + 6.260000000000000E-07 1.045149305253753E+00 2.462326221730596E-04 1.330255649811740E+00 + 6.270000000000000E-07 1.045158155208996E+00 2.400150240452965E-04 1.330225445086910E+00 + 6.280000000000000E-07 1.045167543308349E+00 2.351014711807756E-04 1.330195327095500E+00 + 6.290000000000000E-07 1.045177106856348E+00 2.295070011777853E-04 1.330165295199200E+00 + 6.300000000000000E-07 1.045186321622641E+00 2.225635381454606E-04 1.330135348765460E+00 + 6.310000000000000E-07 1.045194792170897E+00 2.172613687535073E-04 1.330105487167400E+00 + 6.320000000000000E-07 1.045203500771309E+00 2.149425311305904E-04 1.330075709783760E+00 + 6.330000000000000E-07 1.045213256484557E+00 2.126151330176184E-04 1.330046015998860E+00 + 6.340000000000000E-07 1.045223553101800E+00 2.102791743765743E-04 1.330016405202500E+00 + 6.350000000000000E-07 1.045234013122680E+00 2.072644387124910E-04 1.329986876789960E+00 + 6.360000000000000E-07 1.045244216794131E+00 2.040711801010856E-04 1.329957430161900E+00 + 6.370000000000000E-07 1.045253740403829E+00 2.006985957335607E-04 1.329928064724300E+00 + 6.380000000000000E-07 1.045263382110148E+00 2.008497337040627E-04 1.329898779888440E+00 + 6.390000000000000E-07 1.045273947322694E+00 2.008317140156538E-04 1.329869575070820E+00 + 6.400000000000000E-07 1.045285057551241E+00 2.008126259646895E-04 1.329840449693100E+00 + 6.410000000000000E-07 1.045296551000445E+00 2.009616291914741E-04 1.329811403182080E+00 + 6.420000000000000E-07 1.045308363549875E+00 2.012795264863135E-04 1.329782434969600E+00 + 6.430000000000000E-07 1.045320492129501E+00 2.015974258466275E-04 1.329753544492540E+00 + 6.440000000000000E-07 1.045332967563536E+00 2.022552520949957E-04 1.329724731192700E+00 + 6.450000000000000E-07 1.045345759899224E+00 2.030843808416504E-04 1.329695994516830E+00 + 6.460000000000000E-07 1.045358824317712E+00 2.035741221492084E-04 1.329667333916530E+00 + 6.470000000000000E-07 1.045372099342502E+00 2.042351660342150E-04 1.329638748848190E+00 + 6.480000000000000E-07 1.045385551462448E+00 2.047262497537979E-04 1.329610238772990E+00 + 6.489999999999999E-07 1.045399124910904E+00 2.052178709316164E-04 1.329581803156790E+00 + 6.500000000000000E-07 1.045412800938988E+00 2.058815975994723E-04 1.329553441470150E+00 + 6.510000000000000E-07 1.045426582720485E+00 2.063745614018979E-04 1.329525153188210E+00 + 6.520000000000000E-07 1.045440471398717E+00 2.070401659792191E-04 1.329496937790700E+00 + 6.530000000000000E-07 1.045454454497395E+00 2.075344725484678E-04 1.329468794761880E+00 + 6.540000000000000E-07 1.045468464599443E+00 2.080293167479006E-04 1.329440723590460E+00 + 6.550000000000000E-07 1.045482498057374E+00 2.086976046627449E-04 1.329412723769620E+00 + 6.560000000000001E-07 1.045496553339788E+00 2.091937918392058E-04 1.329384794796890E+00 + 6.570000000000000E-07 1.045510635129543E+00 2.098639580351773E-04 1.329356936174180E+00 + 6.580000000000000E-07 1.045524666346868E+00 2.103614883294010E-04 1.329329147407660E+00 + 6.590000000000000E-07 1.045538559400521E+00 2.108595564235613E-04 1.329301428007780E+00 + 6.600000000000000E-07 1.045552348654994E+00 2.115324065074228E-04 1.329273777489210E+00 + 6.610000000000000E-07 1.045566110625959E+00 2.129043768487645E-04 1.329246195370780E+00 + 6.620000000000000E-07 1.045580172412483E+00 2.148038996275277E-04 1.329218681175440E+00 + 6.630000000000000E-07 1.045594198084655E+00 2.156579601845352E-04 1.329191234430250E+00 + 6.640000000000000E-07 1.045608130833515E+00 2.168642586307475E-04 1.329163854666310E+00 + 6.650000000000000E-07 1.045621924109802E+00 2.191267300925995E-04 1.329136541418720E+00 + 6.660000000000000E-07 1.045636494884311E+00 2.226260420764523E-04 1.329109294226560E+00 + 6.670000000000000E-07 1.045651816979515E+00 2.257827567840655E-04 1.329082112632820E+00 + 6.680000000000000E-07 1.045667654294284E+00 2.284188832691462E-04 1.329054996184400E+00 + 6.690000000000000E-07 1.045683780252616E+00 2.298254013270867E-04 1.329027944432050E+00 + 6.700000000000000E-07 1.045699609361488E+00 2.305274520667250E-04 1.329000956930320E+00 + 6.710000000000000E-07 1.045715078138228E+00 2.314077646618846E-04 1.328974033237530E+00 + 6.720000000000000E-07 1.045730267886692E+00 2.321122304430071E-04 1.328947172915770E+00 + 6.730000000000000E-07 1.045745266858880E+00 2.335286639357697E-04 1.328920375530800E+00 + 6.740000000000000E-07 1.045760788195037E+00 2.363722416512695E-04 1.328893640652060E+00 + 6.750000000000000E-07 1.045776625515359E+00 2.379755054093924E-04 1.328866967852620E+00 + 6.760000000000000E-07 1.045792600957640E+00 2.386898873686543E-04 1.328840356709140E+00 + 6.770000000000000E-07 1.045808247514113E+00 2.390477550340525E-04 1.328813806801840E+00 + 6.780000000000000E-07 1.045823119905555E+00 2.386893788440745E-04 1.328787317714450E+00 + 6.790000000000000E-07 1.045837555771692E+00 2.388668530728809E-04 1.328760889034230E+00 + 6.799999999999999E-07 1.045851677897483E+00 2.410193638974664E-04 1.328734520351850E+00 + 6.810000000000000E-07 1.045866374227087E+00 2.440765557054830E-04 1.328708211261430E+00 + 6.820000000000000E-07 1.045881785624778E+00 2.467815163039944E-04 1.328681961360470E+00 + 6.830000000000000E-07 1.045897803967143E+00 2.502150407965931E-04 1.328655770249830E+00 + 6.840000000000000E-07 1.045914494666557E+00 2.523929950830751E-04 1.328629637533700E+00 + 6.850000000000000E-07 1.045930993921758E+00 2.525860157679779E-04 1.328603562819540E+00 + 6.860000000000000E-07 1.045946946210705E+00 2.524160983200198E-04 1.328577545718090E+00 + 6.870000000000001E-07 1.045962222578169E+00 2.520631075330667E-04 1.328551585843320E+00 + 6.880000000000000E-07 1.045976831730619E+00 2.529801492320116E-04 1.328525682812380E+00 + 6.890000000000000E-07 1.045991376191637E+00 2.537167945523808E-04 1.328499836245610E+00 + 6.900000000000000E-07 1.046005982980502E+00 2.557304325072635E-04 1.328474045766460E+00 + 6.910000000000000E-07 1.046021138681571E+00 2.586616014485225E-04 1.328448311001520E+00 + 6.920000000000000E-07 1.046036945389053E+00 2.619658887977861E-04 1.328422631580440E+00 + 6.930000000000000E-07 1.046052941616904E+00 2.625325898993536E-04 1.328397007135900E+00 + 6.940000000000000E-07 1.046068836917074E+00 2.630998304190630E-04 1.328371437303630E+00 + 6.950000000000000E-07 1.046084299017366E+00 2.636676103898158E-04 1.328345921722330E+00 + 6.960000000000000E-07 1.046099553526171E+00 2.644198101017839E-04 1.328320460033680E+00 + 6.970000000000000E-07 1.046114617760768E+00 2.649889367910856E-04 1.328295051882260E+00 + 6.980000000000000E-07 1.046129367196580E+00 2.655586030323879E-04 1.328269696915600E+00 + 6.990000000000000E-07 1.046143900278031E+00 2.661288088582095E-04 1.328244394784070E+00 + 7.000000000000000E-07 1.046158216699328E+00 2.668845054481957E-04 1.328219145140910E+00 + 7.010000000000000E-07 1.046171853334797E+00 2.674560582699037E-04 1.328193947642190E+00 + 7.020000000000000E-07 1.046185486369156E+00 2.683991239950843E-04 1.328168801946760E+00 + 7.030000000000000E-07 1.046199164931322E+00 2.725016242455469E-04 1.328143707716250E+00 + 7.040000000000000E-07 1.046214030872363E+00 2.768008603975542E-04 1.328118664615030E+00 + 7.050000000000000E-07 1.046229607267121E+00 2.794347376261044E-04 1.328093672310210E+00 + 7.060000000000000E-07 1.046245758846082E+00 2.828212754461518E-04 1.328068730471570E+00 + 7.070000000000000E-07 1.046261938294243E+00 2.832271806183120E-04 1.328043838771560E+00 + 7.080000000000000E-07 1.046277640167732E+00 2.838201834351253E-04 1.328018996885300E+00 + 7.090000000000000E-07 1.046292931827498E+00 2.842263655385839E-04 1.327994204490490E+00 + 7.100000000000000E-07 1.046307752529465E+00 2.848201808467138E-04 1.327969461267460E+00 + 7.109999999999999E-07 1.046322261845741E+00 2.854145363243779E-04 1.327944766899080E+00 + 7.120000000000000E-07 1.046336495160005E+00 2.858212679038598E-04 1.327920121070800E+00 + 7.130000000000000E-07 1.046350900922971E+00 2.881123228439462E-04 1.327895523470560E+00 + 7.140000000000000E-07 1.046365517830787E+00 2.900313386182835E-04 1.327870973788830E+00 + 7.150000000000000E-07 1.046380363905696E+00 2.906318716768448E-04 1.327846471718530E+00 + 7.160000000000000E-07 1.046394729835494E+00 2.908544747485000E-04 1.327822016955050E+00 + 7.170000000000000E-07 1.046408872465385E+00 2.910765471158864E-04 1.327797609196200E+00 + 7.180000000000001E-07 1.046422737012463E+00 2.922469425013393E-04 1.327773248142210E+00 + 7.190000000000000E-07 1.046436253332092E+00 2.924692924222563E-04 1.327748933495680E+00 + 7.200000000000000E-07 1.046449665635920E+00 2.932620306364937E-04 1.327724664961610E+00 + 7.210000000000000E-07 1.046462990264283E+00 2.950087154876983E-04 1.327700442247290E+00 + 7.220000000000000E-07 1.046476396774112E+00 2.975225220322731E-04 1.327676265062370E+00 + 7.230000000000000E-07 1.046489793136021E+00 2.975577990728036E-04 1.327652133118790E+00 + 7.240000000000000E-07 1.046503047673474E+00 2.991230299716263E-04 1.327628046130760E+00 + 7.250000000000000E-07 1.046516479213546E+00 3.008831249119806E-04 1.327604003814760E+00 + 7.260000000000000E-07 1.046530020906885E+00 3.030308005065120E-04 1.327580005889490E+00 + 7.270000000000000E-07 1.046543461106678E+00 3.026849499014217E-04 1.327556052075900E+00 + 7.280000000000000E-07 1.046556651803606E+00 3.046463468257975E-04 1.327532142097090E+00 + 7.290000000000000E-07 1.046569517786500E+00 3.050703016468484E-04 1.327508275678380E+00 + 7.300000000000000E-07 1.046581771727620E+00 3.056872461229752E-04 1.327484452547220E+00 + 7.310000000000000E-07 1.046593741831144E+00 3.053384705679599E-04 1.327460672433210E+00 + 7.320000000000000E-07 1.046605199691983E+00 3.073097980091759E-04 1.327436935068060E+00 + 7.330000000000000E-07 1.046616621102129E+00 3.085102646589948E-04 1.327413240185570E+00 + 7.340000000000000E-07 1.046628017565891E+00 3.110712691868297E-04 1.327389587521650E+00 + 7.350000000000000E-07 1.046639910218690E+00 3.138324950843339E-04 1.327365976814230E+00 + 7.360000000000000E-07 1.046651943151609E+00 3.173785201096158E-04 1.327342407803330E+00 + 7.370000000000000E-07 1.046664212729765E+00 3.187896703635894E-04 1.327318880230940E+00 + 7.380000000000000E-07 1.046676007554394E+00 3.196181235299881E-04 1.327295393841100E+00 + 7.390000000000000E-07 1.046687437975362E+00 3.202522588413494E-04 1.327271948379810E+00 + 7.400000000000000E-07 1.046697279881755E+00 3.202999472436175E-04 1.327248543595060E+00 + 7.410000000000000E-07 1.046706135060215E+00 3.189750555484093E-04 1.327225179236760E+00 + 7.419999999999999E-07 1.046713862580072E+00 3.211769207969832E-04 1.327201855056800E+00 + 7.430000000000000E-07 1.046721871614868E+00 3.271164739876647E-04 1.327178570808950E+00 + 7.440000000000000E-07 1.046730272175331E+00 3.346449073390914E-04 1.327155326248880E+00 + 7.450000000000000E-07 1.046740501259060E+00 3.421926333679019E-04 1.327132121134170E+00 + 7.460000000000000E-07 1.046751680426605E+00 3.497596524765032E-04 1.327108955224250E+00 + 7.470000000000000E-07 1.046763868052194E+00 3.573459650688986E-04 1.327085828280380E+00 + 7.480000000000000E-07 1.046777912415171E+00 3.649515715506583E-04 1.327062740065700E+00 + 7.490000000000000E-07 1.046792102242083E+00 3.707938097818990E-04 1.327039690345110E+00 + 7.500000000000000E-07 1.046804357080523E+00 3.631632982600370E-04 1.327016678885350E+00 + 7.510000000000000E-07 1.046815572343059E+00 3.688176560008942E-04 1.326993705454940E+00 + 7.520000000000000E-07 1.046825449654038E+00 3.744859510124802E-04 1.326970769824140E+00 + 7.530000000000000E-07 1.046836760483276E+00 3.801681835983927E-04 1.326947871765000E+00 + 7.540000000000000E-07 1.046848587411650E+00 3.858643540632899E-04 1.326925011051280E+00 + 7.550000000000000E-07 1.046860864366256E+00 3.889786095689200E-04 1.326902187458480E+00 + 7.560000000000000E-07 1.046872894123829E+00 3.920998378577586E-04 1.326879400763780E+00 + 7.570000000000000E-07 1.046884411012915E+00 3.942269549022739E-04 1.326856650746090E+00 + 7.580000000000000E-07 1.046895701267821E+00 3.965588506254991E-04 1.326833937185960E+00 + 7.590000000000000E-07 1.046906835569215E+00 3.988955762257961E-04 1.326811259865620E+00 + 7.600000000000000E-07 1.046917101288162E+00 4.010361112175749E-04 1.326788618568950E+00 + 7.610000000000000E-07 1.046926788333640E+00 3.981487270334402E-04 1.326766013081440E+00 + 7.620000000000000E-07 1.046933540127413E+00 3.952522405187550E-04 1.326743443190220E+00 + 7.630000000000000E-07 1.046938628178191E+00 3.921448270953880E-04 1.326720908684020E+00 + 7.640000000000000E-07 1.046941069178355E+00 3.910486986268732E-04 1.326698409353150E+00 + 7.650000000000000E-07 1.046941607565751E+00 3.915671726970346E-04 1.326675944989510E+00 + 7.660000000000000E-07 1.046941050594559E+00 3.955303342174100E-04 1.326653515386550E+00 + 7.670000000000000E-07 1.046939652476820E+00 4.031547449986886E-04 1.326631120339270E+00 + 7.680000000000000E-07 1.046938078300245E+00 4.116105767484131E-04 1.326608759644200E+00 + 7.690000000000000E-07 1.046935899260892E+00 4.200873169509417E-04 1.326586433099410E+00 + 7.700000000000000E-07 1.046933800541355E+00 4.330663678085595E-04 1.326564140504460E+00 + 7.710000000000000E-07 1.046931842347250E+00 4.495455754419827E-04 1.326541881660400E+00 + 7.720000000000000E-07 1.046931171415364E+00 4.695386071357624E-04 1.326519656369770E+00 + 7.730000000000001E-07 1.046931603412790E+00 4.914230986127136E-04 1.326497464436590E+00 + 7.740000000000000E-07 1.046935057936450E+00 5.133633386726560E-04 1.326475305666330E+00 + 7.750000000000000E-07 1.046939554268154E+00 5.425357270477449E-04 1.326453179865870E+00 + 7.760000000000000E-07 1.046948930406268E+00 5.699348540678753E-04 1.326431086843580E+00 + 7.770000000000000E-07 1.046959325716217E+00 5.974036690350955E-04 1.326409026409200E+00 + 7.780000000000000E-07 1.046974508683986E+00 6.249421734908905E-04 1.326386998373890E+00 + 7.790000000000000E-07 1.046991148046091E+00 6.527564809111620E-04 1.326365002550200E+00 + 7.800000000000000E-07 1.047011022384206E+00 6.794027373289624E-04 1.326343038752090E+00 + 7.810000000000000E-07 1.047030804279599E+00 6.933038356028019E-04 1.326321106794840E+00 + 7.820000000000000E-07 1.047050489959484E+00 7.070318001879980E-04 1.326299206495120E+00 + 7.830000000000000E-07 1.047069563350390E+00 7.253510589942010E-04 1.326277337670950E+00 + 7.840000000000000E-07 1.047088291011400E+00 7.435079064039188E-04 1.326255500141660E+00 + 7.850000000000000E-07 1.047106424922806E+00 7.567239792560870E-04 1.326233693727920E+00 + 7.860000000000000E-07 1.047124393007007E+00 7.695557173642177E-04 1.326211918251710E+00 + 7.870000000000000E-07 1.047141029801092E+00 7.822097693438653E-04 1.326190173536300E+00 + 7.880000000000000E-07 1.047157976621515E+00 8.038603920522829E-04 1.326168459406250E+00 + 7.890000000000000E-07 1.047179224925566E+00 8.278608148273258E-04 1.326146775687410E+00 + 7.900000000000000E-07 1.047201127980211E+00 8.517111621615959E-04 1.326125122206880E+00 + 7.910000000000000E-07 1.047225524410025E+00 8.655722047617181E-04 1.326103498793040E+00 + 7.920000000000000E-07 1.047248608976219E+00 8.725492888296510E-04 1.326081905275470E+00 + 7.930000000000000E-07 1.047269117642915E+00 8.797512649178575E-04 1.326060341485040E+00 + 7.940000000000000E-07 1.047289463317727E+00 8.924322472499099E-04 1.326038807253790E+00 + 7.950000000000000E-07 1.047309613452492E+00 9.112443464370387E-04 1.326017302415020E+00 + 7.960000000000000E-07 1.047331668415949E+00 9.296801609789082E-04 1.325995826803190E+00 + 7.970000000000000E-07 1.047355299038291E+00 9.452068444725653E-04 1.325974380253980E+00 + 7.980000000000000E-07 1.047379136045100E+00 9.552786828842519E-04 1.325952962604240E+00 + 7.990000000000000E-07 1.047403072830967E+00 9.657960036606947E-04 1.325931573691990E+00 + 8.000000000000000E-07 1.047427193866072E+00 9.790895649030144E-04 1.325910213356430E+00 + 8.010000000000000E-07 1.047451361506094E+00 9.924137067527394E-04 1.325888881437880E+00 + 8.020000000000000E-07 1.047477610041293E+00 1.005556152948590E-03 1.325867577777830E+00 + 8.030000000000000E-07 1.047504359325982E+00 1.018516099946340E-03 1.325846302218880E+00 + 8.040000000000001E-07 1.047536241473395E+00 1.031292744178663E-03 1.325825054604780E+00 + 8.050000000000000E-07 1.047568907885968E+00 1.048146908027978E-03 1.325803834780350E+00 + 8.060000000000000E-07 1.047601530876990E+00 1.040078355092751E-03 1.325782642591560E+00 + 8.070000000000000E-07 1.047634095398597E+00 1.034976873960194E-03 1.325761477885440E+00 + 8.080000000000000E-07 1.047660247782939E+00 1.027292698511066E-03 1.325740340510110E+00 + 8.090000000000000E-07 1.047686099812982E+00 1.019800174887994E-03 1.325719230314780E+00 + 8.100000000000000E-07 1.047707779443527E+00 1.012071263221199E-03 1.325698147149700E+00 + 8.110000000000000E-07 1.047729237681430E+00 9.995966685213660E-04 1.325677090866210E+00 + 8.120000000000000E-07 1.047743536155764E+00 9.827886098305182E-04 1.325656061316680E+00 + 8.130000000000000E-07 1.047757488790403E+00 9.661507499608610E-04 1.325635058354500E+00 + 8.140000000000000E-07 1.047763804444236E+00 9.490374090644322E-04 1.325614081834110E+00 + 8.150000000000000E-07 1.047769969728286E+00 9.448242560608527E-04 1.325593131610980E+00 + 8.160000000000000E-07 1.047773137187235E+00 9.401656527267664E-04 1.325572207541580E+00 + 8.170000000000000E-07 1.047776216953897E+00 9.354925875861295E-04 1.325551309483370E+00 + 8.180000000000000E-07 1.047777270288282E+00 9.308050603294224E-04 1.325530437294830E+00 + 8.190000000000000E-07 1.047778186624244E+00 9.269704113819045E-04 1.325509590835410E+00 + 8.200000000000000E-07 1.047774802403966E+00 9.233405483764064E-04 1.325488769965540E+00 + 8.210000000000000E-07 1.047771400536673E+00 9.194815325215046E-04 1.325467974546630E+00 + 8.220000000000000E-07 1.047764370361179E+00 9.177865956932137E-04 1.325447204441030E+00 + 8.230000000000000E-07 1.047757377048223E+00 9.256726507307197E-04 1.325426459512080E+00 + 8.240000000000000E-07 1.047749487584204E+00 9.337935272091834E-04 1.325405739624020E+00 + 8.250000000000000E-07 1.047741618195193E+00 9.417131416291529E-04 1.325385044642050E+00 + 8.260000000000000E-07 1.047731838654609E+00 9.498681147314757E-04 1.325364374432310E+00 + 8.270000000000000E-07 1.047721973752322E+00 9.578212902678193E-04 1.325343728861850E+00 + 8.279999999999999E-07 1.047707491516565E+00 9.668873570831467E-04 1.325323107798620E+00 + 8.290000000000000E-07 1.047692462294499E+00 9.779484068345122E-04 1.325302511111490E+00 + 8.300000000000000E-07 1.047673226736550E+00 9.929897488926683E-04 1.325281938670230E+00 + 8.310000000000000E-07 1.047653141035367E+00 1.015986829694636E-03 1.325261390345490E+00 + 8.320000000000000E-07 1.047633752137463E+00 1.044104435780425E-03 1.325240866008810E+00 + 8.330000000000000E-07 1.047614638723221E+00 1.085522902123505E-03 1.325220365532600E+00 + 8.340000000000000E-07 1.047600760948479E+00 1.132339618517545E-03 1.325199888790140E+00 + 8.350000000000001E-07 1.047589518185891E+00 1.182141510046346E-03 1.325179435655570E+00 + 8.360000000000000E-07 1.047582936855137E+00 1.232060910894234E-03 1.325159006003890E+00 + 8.370000000000000E-07 1.047579827948053E+00 1.281876161129743E-03 1.325138599710930E+00 + 8.380000000000000E-07 1.047578951020551E+00 1.332030321081121E-03 1.325118216653370E+00 + 8.390000000000000E-07 1.047580594183633E+00 1.382301999067603E-03 1.325097856708720E+00 + 8.400000000000000E-07 1.047583955198737E+00 1.432691198007155E-03 1.325077519755300E+00 + 8.410000000000000E-07 1.047590436261947E+00 1.488543580847245E-03 1.325057205672260E+00 + 8.420000000000000E-07 1.047598319842640E+00 1.544526367039972E-03 1.325036914339570E+00 + 8.430000000000000E-07 1.047610935731388E+00 1.600639559836751E-03 1.325016645637990E+00 + 8.440000000000000E-07 1.047624127180442E+00 1.656883162494959E-03 1.324996399449070E+00 + 8.450000000000000E-07 1.047643117657286E+00 1.713257178278004E-03 1.324976175655160E+00 + 8.460000000000000E-07 1.047662084005794E+00 1.769761610455194E-03 1.324955974139380E+00 + 8.470000000000000E-07 1.047681999317462E+00 1.813160079636676E-03 1.324935794785660E+00 + 8.480000000000000E-07 1.047701946094193E+00 1.859801906390260E-03 1.324915637478650E+00 + 8.490000000000000E-07 1.047715295897636E+00 1.906550012807259E-03 1.324895502103800E+00 + 8.500000000000000E-07 1.047726844952281E+00 1.945298989124846E-03 1.324875388547310E+00 + 8.510000000000000E-07 1.047738549879133E+00 1.989770398984366E-03 1.324855296696110E+00 + 8.520000000000000E-07 1.047750348897479E+00 2.095729044907462E-03 1.324835226437890E+00 + 8.530000000000000E-07 1.047767605294426E+00 2.133017683744185E-03 1.324815177661090E+00 + 8.540000000000000E-07 1.047791905510931E+00 2.218348403502016E-03 1.324795150254850E+00 + 8.550000000000000E-07 1.047818418882227E+00 2.310217348244026E-03 1.324775144109060E+00 + 8.560000000000000E-07 1.047851843468609E+00 2.395721729565107E-03 1.324755159114300E+00 + 8.570000000000000E-07 1.047885588632721E+00 2.458038095689459E-03 1.324735195161900E+00 + 8.580000000000000E-07 1.047928039204696E+00 2.551634788315200E-03 1.324715252143870E+00 + 8.589999999999999E-07 1.047970412470177E+00 2.624737227819212E-03 1.324695329952930E+00 + 8.600000000000000E-07 1.048023220548083E+00 2.723522375062108E-03 1.324675428482480E+00 + 8.610000000000000E-07 1.048078126196254E+00 2.789002004916494E-03 1.324655547626620E+00 + 8.620000000000000E-07 1.048140300421103E+00 2.899160408289475E-03 1.324635687280130E+00 + 8.630000000000000E-07 1.048207241938472E+00 2.960412745879289E-03 1.324615847338480E+00 + 8.640000000000000E-07 1.048278409160125E+00 3.054534627917653E-03 1.324596027697780E+00 + 8.650000000000000E-07 1.048357052707288E+00 3.132368488989489E-03 1.324576228254840E+00 + 8.660000000000001E-07 1.048436830487660E+00 3.213359081793805E-03 1.324556448907090E+00 + 8.670000000000000E-07 1.048526428644327E+00 3.259154960386772E-03 1.324536689552640E+00 + 8.680000000000000E-07 1.048615843313936E+00 3.357024712042136E-03 1.324516950090250E+00 + 8.690000000000000E-07 1.048715148218825E+00 3.489421142093250E-03 1.324497230419300E+00 + 8.700000000000000E-07 1.048815971181798E+00 3.531293371955145E-03 1.324477530439830E+00 + 8.710000000000000E-07 1.048926857924925E+00 3.564945421308734E-03 1.324457850052500E+00 + 8.720000000000000E-07 1.049044908795749E+00 3.679304100653368E-03 1.324438189158590E+00 + 8.730000000000000E-07 1.049166413128009E+00 3.726604082701499E-03 1.324418547660020E+00 + 8.740000000000000E-07 1.049296428879817E+00 3.770298531458578E-03 1.324398925459290E+00 + 8.750000000000000E-07 1.049426171616863E+00 3.818953357363171E-03 1.324379322459550E+00 + 8.760000000000000E-07 1.049564210263465E+00 3.868175047738017E-03 1.324359738564520E+00 + 8.770000000000000E-07 1.049701958513664E+00 3.907274996952621E-03 1.324340173678550E+00 + 8.780000000000000E-07 1.049845994882515E+00 3.934123795458572E-03 1.324320627706550E+00 + 8.790000000000000E-07 1.049992662084591E+00 3.955200886907125E-03 1.324301100554040E+00 + 8.800000000000000E-07 1.050141519573849E+00 3.967221321383652E-03 1.324281592127130E+00 + 8.810000000000000E-07 1.050294307337905E+00 3.988597805776602E-03 1.324262102332480E+00 + 8.820000000000000E-07 1.050447085007972E+00 3.978224625708395E-03 1.324242631077360E+00 + 8.830000000000000E-07 1.050607609986077E+00 3.968051170294418E-03 1.324223178269580E+00 + 8.840000000000000E-07 1.050767796843387E+00 3.955969767958174E-03 1.324203743817520E+00 + 8.850000000000000E-07 1.050923092122999E+00 3.927901360070584E-03 1.324184327630120E+00 + 8.860000000000000E-07 1.051076031953918E+00 3.841759824159816E-03 1.324164929616880E+00 + 8.870000000000000E-07 1.051227183758931E+00 3.798665924324868E-03 1.324145549687840E+00 + 8.880000000000000E-07 1.051375087152381E+00 3.745578981754750E-03 1.324126187753590E+00 + 8.890000000000000E-07 1.051522682487672E+00 3.701078944638891E-03 1.324106843725250E+00 + 8.899999999999999E-07 1.051666311265850E+00 3.661657321875855E-03 1.324087517514480E+00 + 8.910000000000000E-07 1.051809500042060E+00 3.571127079571336E-03 1.324068209033480E+00 + 8.920000000000000E-07 1.051943498022589E+00 3.351531098795739E-03 1.324048918194950E+00 + 8.930000000000000E-07 1.052070762485258E+00 3.324330688851290E-03 1.324029644912130E+00 + 8.940000000000000E-07 1.052194726230713E+00 3.208200037721505E-03 1.324010389098780E+00 + 8.950000000000000E-07 1.052304860963421E+00 3.056687924568558E-03 1.323991150669150E+00 + 8.960000000000000E-07 1.052414773434226E+00 2.942349178462081E-03 1.323971929538010E+00 + 8.970000000000001E-07 1.052502780022198E+00 2.834877091402793E-03 1.323952725620640E+00 + 8.980000000000000E-07 1.052584273861886E+00 2.666932219079857E-03 1.323933538832790E+00 + 8.990000000000000E-07 1.052659626877916E+00 2.501460898299947E-03 1.323914369090740E+00 + 9.000000000000000E-07 1.052723802933926E+00 2.379029853183195E-03 1.323895216311230E+00 + 9.010000000000000E-07 1.052787858771010E+00 2.261571369321692E-03 1.323876080411500E+00 + 9.020000000000000E-07 1.052821869541360E+00 2.080000935934572E-03 1.323856961309270E+00 + 9.030000000000000E-07 1.052853084720993E+00 1.890117567359862E-03 1.323837858922710E+00 + 9.040000000000000E-07 1.052866397922094E+00 1.724726050468831E-03 1.323818773170500E+00 + 9.050000000000000E-07 1.052859379171112E+00 1.640289877936772E-03 1.323799703971770E+00 + 9.060000000000000E-07 1.052852413253518E+00 1.487688805686849E-03 1.323780651246100E+00 + 9.070000000000000E-07 1.052848611642345E+00 1.376099936169961E-03 1.323761614913550E+00 + 9.080000000000000E-07 1.052844838670080E+00 1.333344566482168E-03 1.323742594894630E+00 + 9.090000000000000E-07 1.052835368829740E+00 1.189275458028569E-03 1.323723591110280E+00 + 9.100000000000000E-07 1.052820869248812E+00 1.054290942605935E-03 1.323704603481930E+00 + 9.110000000000000E-07 1.052805450313970E+00 9.366342882581003E-04 1.323685631931410E+00 + 9.120000000000000E-07 1.052767449737131E+00 8.436181167945962E-04 1.323666676381010E+00 + 9.130000000000000E-07 1.052729551679777E+00 7.317542083422428E-04 1.323647736753440E+00 + 9.140000000000000E-07 1.052687896630913E+00 6.763454072570361E-04 1.323628812971860E+00 + 9.150000000000000E-07 1.052643047484057E+00 5.652567173572225E-04 1.323609904959850E+00 + 9.160000000000000E-07 1.052598290617709E+00 5.129373259408466E-04 1.323591012641400E+00 + 9.170000000000000E-07 1.052552554797870E+00 4.583126029690992E-04 1.323572135940930E+00 + 9.180000000000000E-07 1.052506937663857E+00 4.038093454985267E-04 1.323553274783280E+00 + 9.190000000000000E-07 1.052462313531449E+00 3.586879638784663E-04 1.323534429093690E+00 + 9.200000000000000E-07 1.052418708673314E+00 3.120024706391773E-04 1.323515598797820E+00 + 9.209999999999999E-07 1.052375217601570E+00 2.700975173655584E-04 1.323496783821710E+00 + 9.220000000000000E-07 1.052334018043352E+00 2.391012133636555E-04 1.323477984091820E+00 + 9.230000000000000E-07 1.052293047791844E+00 2.095046342505333E-04 1.323459199535020E+00 + 9.240000000000000E-07 1.052254431963422E+00 1.862129546266315E-04 1.323440430078530E+00 + 9.250000000000000E-07 1.052219445428962E+00 1.636055562032376E-04 1.323421675650000E+00 + 9.260000000000000E-07 1.052184553785031E+00 1.480692683441966E-04 1.323402936177440E+00 + 9.270000000000000E-07 1.052153554316969E+00 1.324985958014863E-04 1.323384211589260E+00 + 9.280000000000001E-07 1.052123495606230E+00 1.168935376472338E-04 1.323365501814220E+00 + 9.290000000000000E-07 1.052094536477733E+00 1.010077326289579E-04 1.323346806781490E+00 + 9.300000000000000E-07 1.052069073221142E+00 9.248587260842015E-05 1.323328126420590E+00 + 9.310000000000000E-07 1.052043684243811E+00 8.468590033687339E-05 1.323309460661400E+00 + 9.320000000000000E-07 1.052021364441905E+00 7.686873437855818E-05 1.323290809434190E+00 + 9.330000000000000E-07 1.052000888614132E+00 6.903437426761178E-05 1.323272172669560E+00 + 9.340000000000000E-07 1.051980527802691E+00 6.440296793432885E-05 1.323253550298500E+00 + 9.350000000000000E-07 1.051964754164010E+00 6.124918707549857E-05 1.323234942252320E+00 + 9.360000000000000E-07 1.051949034063168E+00 5.808842130256970E-05 1.323216348462700E+00 + 9.370000000000000E-07 1.051934790474261E+00 5.492067042606937E-05 1.323197768861670E+00 + 9.380000000000000E-07 1.051922913853669E+00 5.199471278637858E-05 1.323179203381600E+00 + 9.390000000000000E-07 1.051911082461253E+00 4.632278740626288E-05 1.323160651955180E+00 + 9.400000000000000E-07 1.051900307289650E+00 3.839465831134583E-05 1.323142114515460E+00 + 9.410000000000000E-07 1.051889940149579E+00 3.369393273225441E-05 1.323123590995830E+00 + 9.420000000000000E-07 1.051879766217186E+00 3.148153037298949E-05 1.323105081329980E+00 + 9.430000000000000E-07 1.051871533893133E+00 3.051490221754005E-05 1.323086585451950E+00 + 9.440000000000000E-07 1.051863339025340E+00 2.929573388502167E-05 1.323068103296110E+00 + 9.450000000000000E-07 1.051856470570902E+00 2.882585761204386E-05 1.323049634797120E+00 + 9.460000000000000E-07 1.051851586647983E+00 2.835490692158482E-05 1.323031179890000E+00 + 9.470000000000000E-07 1.051846733136190E+00 2.763168465136746E-05 1.323012738510040E+00 + 9.480000000000000E-07 1.051844453224846E+00 2.715831628037134E-05 1.322994310592870E+00 + 9.490000000000000E-07 1.051843253457583E+00 2.693560805431627E-05 1.322975896074430E+00 + 9.500000000000000E-07 1.051842076403402E+00 2.671236282563884E-05 1.322957494890950E+00 diff --git a/src/test_sschiff_estimator_cylinder.c b/src/test_sschiff_estimator_cylinder.c @@ -37,568 +37,16 @@ #include <star/ssp.h> struct result { - double extinction_E; /* Expected value */ - double extinction_SE; /* Standard error */ + double absorption_E; + double absorption_SE; + double scattering_E; + double scattering_SE; + double extinction_E; + double extinction_SE; + double avg_proj_area_E; + double avg_proj_area_SE; }; -struct optical_properties { double wavelength, n_r, kappa_r, n_medium; }; -static struct optical_properties properties[] = { - { 4.00e-1, 1.035286987991962, 1.642464685700659e-03, 1.341779843605380 }, - { 4.01e-1, 1.035394484151721, 1.624259338990003e-03, 1.341687493674650 }, - { 4.02e-1, 1.035498485175403, 1.594876393429627e-03, 1.341595857289120 }, - { 4.03e-1, 1.035593745454516, 1.563339730364777e-03, 1.341504925814500 }, - { 4.04e-1, 1.035689619601582, 1.544600372647700e-03, 1.341414690748000 }, - { 4.05e-1, 1.035787642729216, 1.511779431504190e-03, 1.341325143715940 }, - { 4.06e-1, 1.035873513905482, 1.471484460290955e-03, 1.341236276471320 }, - { 4.07e-1, 1.035966600915103, 1.462298106531343e-03, 1.341148080891490 }, - { 4.08e-1, 1.036060494953379, 1.407755470265253e-03, 1.341060548975860 }, - { 4.09e-1, 1.036137924348892, 1.368311430748199e-03, 1.340973672843670 }, - { 4.10e-1, 1.036210685221302, 1.330046390562492e-03, 1.340887444731820 }, - { 4.11e-1, 1.036280450985456, 1.300974328047466e-03, 1.340801856992660 }, - { 4.12e-1, 1.036362608081751, 1.286212700588483e-03, 1.340716902091980 }, - { 4.13e-1, 1.036447639684917, 1.236376528960547e-03, 1.340632572606910 }, - { 4.14e-1, 1.036513150647102, 1.179708889220116e-03, 1.340548861223910 }, - { 4.15e-1, 1.036573630808206, 1.146669303767855e-03, 1.340465760736840 }, - { 4.16e-1, 1.036625512534661, 1.088836702310341e-03, 1.340383264045010 }, - { 4.17e-1, 1.036667015511619, 1.062695995814180e-03, 1.340301364151300 }, - { 4.18e-1, 1.036722768320762, 1.052568689767748e-03, 1.340220054160350 }, - { 4.19e-1, 1.036779280227207, 1.006697219981148e-03, 1.340139327276700 }, - { 4.20e-1, 1.036822939800069, 9.722694769944399e-04, 1.340059176803090 }, - { 4.21e-1, 1.036865833600609, 9.469232619753497e-04, 1.339979596138690 }, - { 4.22e-1, 1.036897653219950, 9.040007897511271e-04, 1.339900578777420 }, - { 4.23e-1, 1.036927468207731, 9.031823993726944e-04, 1.339822118306300 }, - { 4.24e-1, 1.036971750881576, 8.924646106519200e-04, 1.339744208403820 }, - { 4.25e-1, 1.037006936112468, 8.594846944392758e-04, 1.339666842838340 }, - { 4.26e-1, 1.037033040892262, 8.523605328290717e-04, 1.339590015466570 }, - { 4.27e-1, 1.037064992312203, 8.534479686759623e-04, 1.339513720232010 }, - { 4.28e-1, 1.037100162121758, 8.506283100332938e-04, 1.339437951163470 }, - { 4.29e-1, 1.037135678599242, 8.522590014316155e-04, 1.339362702373630 }, - { 4.30e-1, 1.037173395465157, 8.565016952270480e-04, 1.339287968057550 }, - { 4.31e-1, 1.037214737528130, 8.639171883406666e-04, 1.339213742491340 }, - { 4.32e-1, 1.037259758643682, 8.691165576390881e-04, 1.339140020030710 }, - { 4.33e-1, 1.037303754455989, 8.678624596520924e-04, 1.339066795109690 }, - { 4.34e-1, 1.037346655520205, 8.725074419904887e-04, 1.338994062239260 }, - { 4.35e-1, 1.037390939381917, 8.756139553518682e-04, 1.338921816006080 }, - { 4.36e-1, 1.037435399005391, 8.788395944917436e-04, 1.338850051071230 }, - { 4.37e-1, 1.037478748753387, 8.799399526191331e-04, 1.338778762168930 }, - { 4.38e-1, 1.037521676917342, 8.850314449839643e-04, 1.338707944105350 }, - { 4.39e-1, 1.037565219615768, 8.888250046471429e-04, 1.338637591757410 }, - { 4.40e-1, 1.037608696500256, 8.936647271660125e-04, 1.338567700071630 }, - { 4.41e-1, 1.037652289517359, 8.980316117027969e-04, 1.338498264062920 }, - { 4.42e-1, 1.037695839416368, 9.024090285591408e-04, 1.338429278813520 }, - { 4.43e-1, 1.037737726850386, 9.038232640135763e-04, 1.338360739471860 }, - { 4.44e-1, 1.037774845995166, 9.041864444660936e-04, 1.338292641251510 }, - { 4.45e-1, 1.037810077590406, 9.145542051765026e-04, 1.338224979430060 }, - { 4.46e-1, 1.037846610963693, 9.257787115899892e-04, 1.338157749348150 }, - { 4.47e-1, 1.037878259630217, 9.335983043538121e-04, 1.338090946408400 }, - { 4.48e-1, 1.037914617228649, 9.702316913975443e-04, 1.338024566074410 }, - { 4.49e-1, 1.037970208539347, 1.007116632500707e-03, 1.337958603869820 }, - { 4.50e-1, 1.038027021857111, 1.009173031428449e-03, 1.337893055377290 }, - { 4.51e-1, 1.038079358354862, 1.027814707250541e-03, 1.337827916237600 }, - { 4.52e-1, 1.038129195790007, 1.028696923944767e-03, 1.337763182148700 }, - { 4.53e-1, 1.038171798198206, 1.042122191740612e-03, 1.337698848864820 }, - { 4.54e-1, 1.038225010269563, 1.074869388950115e-03, 1.337634912195530 }, - { 4.55e-1, 1.038283232720702, 1.078816153504866e-03, 1.337571368004930 }, - { 4.56e-1, 1.038338570451783, 1.092245635394565e-03, 1.337508212210750 }, - { 4.57e-1, 1.038393288814794, 1.094812200306328e-03, 1.337445440783500 }, - { 4.58e-1, 1.038440433963812, 1.093485331340192e-03, 1.337383049745670 }, - { 4.59e-1, 1.038482662253094, 1.103439733475396e-03, 1.337321035170910 }, - { 4.60e-1, 1.038522598315466, 1.113089490368987e-03, 1.337259393183210 }, - { 4.61e-1, 1.038570039357183, 1.148516858552576e-03, 1.337198119956160 }, - { 4.62e-1, 1.038624907459639, 1.162434238059191e-03, 1.337137211712120 }, - { 4.63e-1, 1.038687104626991, 1.195966780613389e-03, 1.337076664721550 }, - { 4.64e-1, 1.038755514562735, 1.199699923934052e-03, 1.337016475302180 }, - { 4.65e-1, 1.038817135628012, 1.200459607388018e-03, 1.336956639818350 }, - { 4.66e-1, 1.038874078123107, 1.201186493947777e-03, 1.336897154680260 }, - { 4.67e-1, 1.038928224588712, 1.201880424343069e-03, 1.336838016343290 }, - { 4.68e-1, 1.038980178159473, 1.204040146415200e-03, 1.336779221307310 }, - { 4.69e-1, 1.039031437175215, 1.207552177169499e-03, 1.336720766115970 }, - { 4.70e-1, 1.039082221078632, 1.211092670367667e-03, 1.336662647356080 }, - { 4.71e-1, 1.039132490661054, 1.215106981436297e-03, 1.336604861656970 }, - { 4.72e-1, 1.039182661408432, 1.219870817381094e-03, 1.336547405689770 }, - { 4.73e-1, 1.039232954048947, 1.224643796577837e-03, 1.336490276166880 }, - { 4.74e-1, 1.039283376885394, 1.229326344827146e-03, 1.336433469841280 }, - { 4.75e-1, 1.039334064443123, 1.234117398218758e-03, 1.336376983505970 }, - { 4.76e-1, 1.039385526833898, 1.238492582229907e-03, 1.336320813993330 }, - { 4.77e-1, 1.039436291084956, 1.239066817439063e-03, 1.336264958174580 }, - { 4.78e-1, 1.039484948734477, 1.239657200984551e-03, 1.336209412959190 }, - { 4.79e-1, 1.039531959467774, 1.240540672112777e-03, 1.336154175294280 }, - { 4.80e-1, 1.039578615206627, 1.244442394725602e-03, 1.336099242164120 }, - { 4.81e-1, 1.039626229298974, 1.248349409793080e-03, 1.336044610589550 }, - { 4.82e-1, 1.039673858625137, 1.250235937286709e-03, 1.335990277627470 }, - { 4.83e-1, 1.039720595019222, 1.251865430581906e-03, 1.335936240370270 }, - { 4.84e-1, 1.039766629301611, 1.253490577728770e-03, 1.335882495945360 }, - { 4.85e-1, 1.039811982223968, 1.255111379800570e-03, 1.335829041514630 }, - { 4.86e-1, 1.039856756652771, 1.256727837858235e-03, 1.335775874273970 }, - { 4.87e-1, 1.039900966005627, 1.258467903243629e-03, 1.335722991452750 }, - { 4.88e-1, 1.039944611565546, 1.260075944188013e-03, 1.335670390313400 }, - { 4.89e-1, 1.039987632682271, 1.261679644222744e-03, 1.335618068150840 }, - { 4.90e-1, 1.040030012263802, 1.263407757979589e-03, 1.335566022292100 }, - { 4.91e-1, 1.040071688175533, 1.264745004221800e-03, 1.335514250095810 }, - { 4.92e-1, 1.040112378686561, 1.266076841520105e-03, 1.335462748951790 }, - { 4.93e-1, 1.040151955742799, 1.267403270885834e-03, 1.335411516280550 }, - { 4.94e-1, 1.040190287896555, 1.268724293318496e-03, 1.335360549532910 }, - { 4.95e-1, 1.040226903080167, 1.270352131527932e-03, 1.335309846189550 }, - { 4.96e-1, 1.040262146697474, 1.273305523147441e-03, 1.335259403760560 }, - { 4.97e-1, 1.040295747313031, 1.276077378609423e-03, 1.335209219785090 }, - { 4.98e-1, 1.040325648091718, 1.278823648922353e-03, 1.335159291830880 }, - { 4.99e-1, 1.040353759919826, 1.289072947859305e-03, 1.335109617493900 }, - { 5.00e-1, 1.040379481152520, 1.299274503512561e-03, 1.335060194397920 }, - { 5.01e-1, 1.040408126491521, 1.325811156480386e-03, 1.335011020194160 }, - { 5.02e-1, 1.040441036092354, 1.344764441067616e-03, 1.334962092560890 }, - { 5.03e-1, 1.040477445253276, 1.374045011533224e-03, 1.334913409203050 }, - { 5.04e-1, 1.040512648034075, 1.387903827917666e-03, 1.334864967851890 }, - { 5.05e-1, 1.040545900251519, 1.418801858069103e-03, 1.334816766264610 }, - { 5.06e-1, 1.040590906700738, 1.468464949020522e-03, 1.334768802224000 }, - { 5.07e-1, 1.040639195455336, 1.482563959735157e-03, 1.334721073538070 }, - { 5.08e-1, 1.040686458775797, 1.525559361876046e-03, 1.334673578039750 }, - { 5.09e-1, 1.040741671425000, 1.559105378939951e-03, 1.334626313586500 }, - { 5.10e-1, 1.040807713567804, 1.603233722101158e-03, 1.334579278060000 }, - { 5.11e-1, 1.040875584410848, 1.613985571121084e-03, 1.334532469365830 }, - { 5.12e-1, 1.040943220348469, 1.649380251826485e-03, 1.334485885433120 }, - { 5.13e-1, 1.041018230706438, 1.669091265929868e-03, 1.334439524214260 }, - { 5.14e-1, 1.041092713713321, 1.676566182400091e-03, 1.334393383684570 }, - { 5.15e-1, 1.041165314107454, 1.687389518793340e-03, 1.334347461841990 }, - { 5.16e-1, 1.041238597163079, 1.692922480695940e-03, 1.334301756706790 }, - { 5.17e-1, 1.041311956633380, 1.695363690601763e-03, 1.334256266321260 }, - { 5.18e-1, 1.041389428013160, 1.700417308885367e-03, 1.334210988749400 }, - { 5.19e-1, 1.041464755955924, 1.682652012825744e-03, 1.334165922076680 }, - { 5.20e-1, 1.041530069541219, 1.658812761280325e-03, 1.334121064409700 }, - { 5.21e-1, 1.041594181696574, 1.657124535655681e-03, 1.334076413875950 }, - { 5.22e-1, 1.041657845088931, 1.630122704418439e-03, 1.334031968623510 }, - { 5.23e-1, 1.041711910952757, 1.607543355228116e-03, 1.333987726820790 }, - { 5.24e-1, 1.041759947842067, 1.582244525828878e-03, 1.333943686656260 }, - { 5.24e-1, 1.041810066348116, 1.583879702931137e-03, 1.333899846338170 }, - { 5.26e-1, 1.041862095522510, 1.557720519804618e-03, 1.333856204094340 }, - { 5.27e-1, 1.041913531175414, 1.546395579258294e-03, 1.333812758171830 }, - { 5.28e-1, 1.041954667678475, 1.499811551277397e-03, 1.333769506836780 }, - { 5.29e-1, 1.041989797310988, 1.498636177737314e-03, 1.333726448374050 }, - { 5.30e-1, 1.042024857461786, 1.469804039124607e-03, 1.333683581087100 }, - { 5.31e-1, 1.042050856148998, 1.445880383759585e-03, 1.333640903297660 }, - { 5.32e-1, 1.042077186838255, 1.447977496199134e-03, 1.333598413345510 }, - { 5.33e-1, 1.042102162421808, 1.427209064601831e-03, 1.333556109588290 }, - { 5.34e-1, 1.042120822250241, 1.426699703280908e-03, 1.333513990401220 }, - { 5.35e-1, 1.042141900864721, 1.433471379762237e-03, 1.333472054176910 }, - { 5.36e-1, 1.042166620404672, 1.444941845941278e-03, 1.333430299325120 }, - { 5.37e-1, 1.042194724365646, 1.457717474653508e-03, 1.333388724272540 }, - { 5.38e-1, 1.042225414187512, 1.471493967735956e-03, 1.333347327462570 }, - { 5.39e-1, 1.042256744777569, 1.484092000868983e-03, 1.333306107355130 }, - { 5.40e-1, 1.042288345839686, 1.500763725452133e-03, 1.333265062426450 }, - { 5.41e-1, 1.042324503859437, 1.528281898750163e-03, 1.333224191168830 }, - { 5.42e-1, 1.042364912605382, 1.549843563702498e-03, 1.333183492090460 }, - { 5.43e-1, 1.042406562887524, 1.568902404377878e-03, 1.333142963715240 }, - { 5.44e-1, 1.042454714949787, 1.608815343234295e-03, 1.333102604582520 }, - { 5.45e-1, 1.042513875213928, 1.637618035281262e-03, 1.333062413247000 }, - { 5.46e-1, 1.042576056800957, 1.653580138686441e-03, 1.333022388278450 }, - { 5.47e-1, 1.042640409898161, 1.675120248233485e-03, 1.332982528261570 }, - { 5.48e-1, 1.042707136014031, 1.685157909842487e-03, 1.332942831795790 }, - { 5.49e-1, 1.042774026206271, 1.695250147639022e-03, 1.332903297495090 }, - { 5.50e-1, 1.042842766321411, 1.705310235117845e-03, 1.332863923987830 }, - { 5.51e-1, 1.042912300283037, 1.705530449406497e-03, 1.332824709916540 }, - { 5.52e-1, 1.042981728226964, 1.708356091134469e-03, 1.332785653937790 }, - { 5.53e-1, 1.043051509491291, 1.704016278805774e-03, 1.332746754721980 }, - { 5.54e-1, 1.043124661120986, 1.708897901535750e-03, 1.332708010953190 }, - { 5.55e-1, 1.043199194602040, 1.692798767844739e-03, 1.332669421329010 }, - { 5.55e-1, 1.043265008027065, 1.670274669430007e-03, 1.332630984560360 }, - { 5.57e-1, 1.043336588355702, 1.676935291139613e-03, 1.332592699371340 }, - { 5.58e-1, 1.043412220918312, 1.644255097722219e-03, 1.332554564499090 }, - { 5.59e-1, 1.043477770799564, 1.615743588209385e-03, 1.332516578693600 }, - { 5.60e-1, 1.043536695477078, 1.579035215844224e-03, 1.332478740717540 }, - { 5.61e-1, 1.043600994183831, 1.581811014957712e-03, 1.332441049346180 }, - { 5.62e-1, 1.043666231315476, 1.535235605622221e-03, 1.332403503367140 }, - { 5.63e-1, 1.043718864830297, 1.473563431577708e-03, 1.332366101580320 }, - { 5.64e-1, 1.043771791283996, 1.484273905192923e-03, 1.332328842797720 }, - { 5.65e-1, 1.043825841694447, 1.417653789592510e-03, 1.332291725843300 }, - { 5.66e-1, 1.043867370062960, 1.377114388401143e-03, 1.332254749552820 }, - { 5.67e-1, 1.043898295606015, 1.324202628263995e-03, 1.332217912773740 }, - { 5.68e-1, 1.043933664352148, 1.332080947209094e-03, 1.332181214365030 }, - { 5.69e-1, 1.043977438602228, 1.304393676650687e-03, 1.332144653197070 }, - { 5.70e-1, 1.044018794045142, 1.259390337125543e-03, 1.332108228151510 }, - { 5.71e-1, 1.044062030345026, 1.259437872346123e-03, 1.332071938121120 }, - { 5.72e-1, 1.044109863822610, 1.216136708576497e-03, 1.332035782009670 }, - { 5.73e-1, 1.044157258612584, 1.187015609995578e-03, 1.331999758731800 }, - { 5.74e-1, 1.044202028283407, 1.143113983172242e-03, 1.331963867212890 }, - { 5.75e-1, 1.044232131241066, 1.063840725432586e-03, 1.331928106388940 }, - { 5.76e-1, 1.044246895319073, 1.023527710863081e-03, 1.331892475206440 }, - { 5.77e-1, 1.044252610112750, 9.697779794138100e-04, 1.331856972622240 }, - { 5.78e-1, 1.044253976775146, 9.497659966550409e-04, 1.331821597603440 }, - { 5.79e-1, 1.044259957725371, 9.352618179707669e-04, 1.331786349127290 }, - { 5.80e-1, 1.044274809362838, 9.395000916553773e-04, 1.331751226181010 }, - { 5.81e-1, 1.044292271809736, 9.117482360018617e-04, 1.331716227761740 }, - { 5.82e-1, 1.044308146494882, 9.039549474109366e-04, 1.331681352876410 }, - { 5.83e-1, 1.044328059952875, 9.029813184193805e-04, 1.331646600541600 }, - { 5.84e-1, 1.044352735493407, 8.971352837393257e-04, 1.331611969783430 }, - { 5.85e-1, 1.044378685632916, 8.831537746544346e-04, 1.331577459637510 }, - { 5.86e-1, 1.044405836424728, 8.786320819678296e-04, 1.331543069148750 }, - { 5.86e-1, 1.044435280426248, 8.692005288892336e-04, 1.331508797371310 }, - { 5.88e-1, 1.044465717884262, 8.558876109011710e-04, 1.331474643368460 }, - { 5.89e-1, 1.044495153084782, 8.403500870864586e-04, 1.331440606212510 }, - { 5.90e-1, 1.044523085786671, 8.224523096408955e-04, 1.331406684984670 }, - { 5.91e-1, 1.044550545262079, 8.104699876011172e-04, 1.331372878774990 }, - { 5.92e-1, 1.044578615724940, 7.960072814810816e-04, 1.331339186682210 }, - { 5.93e-1, 1.044607054549288, 7.816152660959925e-04, 1.331305607813710 }, - { 5.94e-1, 1.044635270377074, 7.646952835940411e-04, 1.331272141285380 }, - { 5.95e-1, 1.044663935653333, 7.513518759906050e-04, 1.331238786221540 }, - { 5.96e-1, 1.044693946287332, 7.357902455879724e-04, 1.331205541754850 }, - { 5.97e-1, 1.044724372554313, 7.167716979679335e-04, 1.331172407026180 }, - { 5.98e-1, 1.044753615647646, 6.946257364021612e-04, 1.331139381184580 }, - { 5.99e-1, 1.044780691526674, 6.712950232982840e-04, 1.331106463387120 }, - { 6.00e-1, 1.044803444015409, 6.412373047096696e-04, 1.331073652798870 }, - { 6.01e-1, 1.044824539307607, 6.249866054301384e-04, 1.331040948592740 }, - { 6.02e-1, 1.044845766620673, 6.034077108891513e-04, 1.331008349949450 }, - { 6.03e-1, 1.044866538700980, 5.784136786976474e-04, 1.330975856057410 }, - { 6.04e-1, 1.044887221935235, 5.676336347762863e-04, 1.330943466112670 }, - { 6.05e-1, 1.044906979832908, 5.365222304755139e-04, 1.330911179318770 }, - { 6.06e-1, 1.044925229429125, 5.166167506612064e-04, 1.330878994886740 }, - { 6.07e-1, 1.044942099825731, 4.987864403064095e-04, 1.330846912034950 }, - { 6.08e-1, 1.044955182314457, 4.646041110856508e-04, 1.330814929989080 }, - { 6.09e-1, 1.044965016968371, 4.527564487138666e-04, 1.330783047981990 }, - { 6.10e-1, 1.044979238163159, 4.409632490707041e-04, 1.330751265253680 }, - { 6.11e-1, 1.044994223279349, 4.219740270792474e-04, 1.330719581051210 }, - { 6.12e-1, 1.045007044211732, 4.020154703506148e-04, 1.330687994628590 }, - { 6.13e-1, 1.045017963621394, 3.842194369970159e-04, 1.330656505246750 }, - { 6.14e-1, 1.045027441486893, 3.679492624268454e-04, 1.330625112173430 }, - { 6.15e-1, 1.045035617615133, 3.519800711874577e-04, 1.330593814683110 }, - { 6.16e-1, 1.045043007111847, 3.395630359374884e-04, 1.330562612056960 }, - { 6.17e-1, 1.045051084952903, 3.326688305061111e-04, 1.330531503582750 }, - { 6.17e-1, 1.045061222694419, 3.259457256833445e-04, 1.330500488554770 }, - { 6.19e-1, 1.045073433158974, 3.187415925020136e-04, 1.330469566273780 }, - { 6.20e-1, 1.045088360723040, 3.171051848434899e-04, 1.330438736046940 }, - { 6.21e-1, 1.045103671053048, 3.017345937277039e-04, 1.330407997187720 }, - { 6.22e-1, 1.045116656866957, 2.864761737697034e-04, 1.330377349015850 }, - { 6.23e-1, 1.045126365845574, 2.711663840968767e-04, 1.330346790857260 }, - { 6.24e-1, 1.045133870688024, 2.592620517622472e-04, 1.330316322043970 }, - { 6.25e-1, 1.045140983506716, 2.527585737208309e-04, 1.330285941914100 }, - { 6.26e-1, 1.045149305253753, 2.462326221730596e-04, 1.330255649811740 }, - { 6.27e-1, 1.045158155208996, 2.400150240452965e-04, 1.330225445086910 }, - { 6.28e-1, 1.045167543308349, 2.351014711807756e-04, 1.330195327095500 }, - { 6.29e-1, 1.045177106856348, 2.295070011777853e-04, 1.330165295199200 }, - { 6.30e-1, 1.045186321622641, 2.225635381454606e-04, 1.330135348765460 }, - { 6.31e-1, 1.045194792170897, 2.172613687535073e-04, 1.330105487167400 }, - { 6.32e-1, 1.045203500771309, 2.149425311305904e-04, 1.330075709783760 }, - { 6.33e-1, 1.045213256484557, 2.126151330176184e-04, 1.330046015998860 }, - { 6.34e-1, 1.045223553101800, 2.102791743765743e-04, 1.330016405202500 }, - { 6.35e-1, 1.045234013122680, 2.072644387124910e-04, 1.329986876789960 }, - { 6.36e-1, 1.045244216794131, 2.040711801010856e-04, 1.329957430161900 }, - { 6.37e-1, 1.045253740403829, 2.006985957335607e-04, 1.329928064724300 }, - { 6.38e-1, 1.045263382110148, 2.008497337040627e-04, 1.329898779888440 }, - { 6.39e-1, 1.045273947322694, 2.008317140156538e-04, 1.329869575070820 }, - { 6.40e-1, 1.045285057551241, 2.008126259646895e-04, 1.329840449693100 }, - { 6.41e-1, 1.045296551000445, 2.009616291914741e-04, 1.329811403182080 }, - { 6.42e-1, 1.045308363549875, 2.012795264863135e-04, 1.329782434969600 }, - { 6.43e-1, 1.045320492129501, 2.015974258466275e-04, 1.329753544492540 }, - { 6.44e-1, 1.045332967563536, 2.022552520949957e-04, 1.329724731192700 }, - { 6.45e-1, 1.045345759899224, 2.030843808416504e-04, 1.329695994516830 }, - { 6.46e-1, 1.045358824317712, 2.035741221492084e-04, 1.329667333916530 }, - { 6.47e-1, 1.045372099342502, 2.042351660342150e-04, 1.329638748848190 }, - { 6.48e-1, 1.045385551462448, 2.047262497537979e-04, 1.329610238772990 }, - { 6.48e-1, 1.045399124910904, 2.052178709316164e-04, 1.329581803156790 }, - { 6.50e-1, 1.045412800938988, 2.058815975994723e-04, 1.329553441470150 }, - { 6.51e-1, 1.045426582720485, 2.063745614018979e-04, 1.329525153188210 }, - { 6.52e-1, 1.045440471398717, 2.070401659792191e-04, 1.329496937790700 }, - { 6.53e-1, 1.045454454497395, 2.075344725484678e-04, 1.329468794761880 }, - { 6.54e-1, 1.045468464599443, 2.080293167479006e-04, 1.329440723590460 }, - { 6.55e-1, 1.045482498057374, 2.086976046627449e-04, 1.329412723769620 }, - { 6.56e-1, 1.045496553339788, 2.091937918392058e-04, 1.329384794796890 }, - { 6.57e-1, 1.045510635129543, 2.098639580351773e-04, 1.329356936174180 }, - { 6.58e-1, 1.045524666346868, 2.103614883294010e-04, 1.329329147407660 }, - { 6.59e-1, 1.045538559400521, 2.108595564235613e-04, 1.329301428007780 }, - { 6.60e-1, 1.045552348654994, 2.115324065074228e-04, 1.329273777489210 }, - { 6.61e-1, 1.045566110625959, 2.129043768487645e-04, 1.329246195370780 }, - { 6.62e-1, 1.045580172412483, 2.148038996275277e-04, 1.329218681175440 }, - { 6.63e-1, 1.045594198084655, 2.156579601845352e-04, 1.329191234430250 }, - { 6.64e-1, 1.045608130833515, 2.168642586307475e-04, 1.329163854666310 }, - { 6.65e-1, 1.045621924109802, 2.191267300925995e-04, 1.329136541418720 }, - { 6.66e-1, 1.045636494884311, 2.226260420764523e-04, 1.329109294226560 }, - { 6.67e-1, 1.045651816979515, 2.257827567840655e-04, 1.329082112632820 }, - { 6.68e-1, 1.045667654294284, 2.284188832691462e-04, 1.329054996184400 }, - { 6.69e-1, 1.045683780252616, 2.298254013270867e-04, 1.329027944432050 }, - { 6.70e-1, 1.045699609361488, 2.305274520667250e-04, 1.329000956930320 }, - { 6.71e-1, 1.045715078138228, 2.314077646618846e-04, 1.328974033237530 }, - { 6.72e-1, 1.045730267886692, 2.321122304430071e-04, 1.328947172915770 }, - { 6.73e-1, 1.045745266858880, 2.335286639357697e-04, 1.328920375530800 }, - { 6.74e-1, 1.045760788195037, 2.363722416512695e-04, 1.328893640652060 }, - { 6.75e-1, 1.045776625515359, 2.379755054093924e-04, 1.328866967852620 }, - { 6.76e-1, 1.045792600957640, 2.386898873686543e-04, 1.328840356709140 }, - { 6.77e-1, 1.045808247514113, 2.390477550340525e-04, 1.328813806801840 }, - { 6.78e-1, 1.045823119905555, 2.386893788440745e-04, 1.328787317714450 }, - { 6.79e-1, 1.045837555771692, 2.388668530728809e-04, 1.328760889034230 }, - { 6.79e-1, 1.045851677897483, 2.410193638974664e-04, 1.328734520351850 }, - { 6.81e-1, 1.045866374227087, 2.440765557054830e-04, 1.328708211261430 }, - { 6.82e-1, 1.045881785624778, 2.467815163039944e-04, 1.328681961360470 }, - { 6.83e-1, 1.045897803967143, 2.502150407965931e-04, 1.328655770249830 }, - { 6.84e-1, 1.045914494666557, 2.523929950830751e-04, 1.328629637533700 }, - { 6.85e-1, 1.045930993921758, 2.525860157679779e-04, 1.328603562819540 }, - { 6.86e-1, 1.045946946210705, 2.524160983200198e-04, 1.328577545718090 }, - { 6.87e-1, 1.045962222578169, 2.520631075330667e-04, 1.328551585843320 }, - { 6.88e-1, 1.045976831730619, 2.529801492320116e-04, 1.328525682812380 }, - { 6.89e-1, 1.045991376191637, 2.537167945523808e-04, 1.328499836245610 }, - { 6.90e-1, 1.046005982980502, 2.557304325072635e-04, 1.328474045766460 }, - { 6.91e-1, 1.046021138681571, 2.586616014485225e-04, 1.328448311001520 }, - { 6.92e-1, 1.046036945389053, 2.619658887977861e-04, 1.328422631580440 }, - { 6.93e-1, 1.046052941616904, 2.625325898993536e-04, 1.328397007135900 }, - { 6.94e-1, 1.046068836917074, 2.630998304190630e-04, 1.328371437303630 }, - { 6.95e-1, 1.046084299017366, 2.636676103898158e-04, 1.328345921722330 }, - { 6.96e-1, 1.046099553526171, 2.644198101017839e-04, 1.328320460033680 }, - { 6.97e-1, 1.046114617760768, 2.649889367910856e-04, 1.328295051882260 }, - { 6.98e-1, 1.046129367196580, 2.655586030323879e-04, 1.328269696915600 }, - { 6.99e-1, 1.046143900278031, 2.661288088582095e-04, 1.328244394784070 }, - { 7.00e-1, 1.046158216699328, 2.668845054481957e-04, 1.328219145140910 }, - { 7.01e-1, 1.046171853334797, 2.674560582699037e-04, 1.328193947642190 }, - { 7.02e-1, 1.046185486369156, 2.683991239950843e-04, 1.328168801946760 }, - { 7.03e-1, 1.046199164931322, 2.725016242455469e-04, 1.328143707716250 }, - { 7.04e-1, 1.046214030872363, 2.768008603975542e-04, 1.328118664615030 }, - { 7.05e-1, 1.046229607267121, 2.794347376261044e-04, 1.328093672310210 }, - { 7.06e-1, 1.046245758846082, 2.828212754461518e-04, 1.328068730471570 }, - { 7.07e-1, 1.046261938294243, 2.832271806183120e-04, 1.328043838771560 }, - { 7.08e-1, 1.046277640167732, 2.838201834351253e-04, 1.328018996885300 }, - { 7.09e-1, 1.046292931827498, 2.842263655385839e-04, 1.327994204490490 }, - { 7.10e-1, 1.046307752529465, 2.848201808467138e-04, 1.327969461267460 }, - { 7.10e-1, 1.046322261845741, 2.854145363243779e-04, 1.327944766899080 }, - { 7.12e-1, 1.046336495160005, 2.858212679038598e-04, 1.327920121070800 }, - { 7.13e-1, 1.046350900922971, 2.881123228439462e-04, 1.327895523470560 }, - { 7.14e-1, 1.046365517830787, 2.900313386182835e-04, 1.327870973788830 }, - { 7.15e-1, 1.046380363905696, 2.906318716768448e-04, 1.327846471718530 }, - { 7.16e-1, 1.046394729835494, 2.908544747485000e-04, 1.327822016955050 }, - { 7.17e-1, 1.046408872465385, 2.910765471158864e-04, 1.327797609196200 }, - { 7.18e-1, 1.046422737012463, 2.922469425013393e-04, 1.327773248142210 }, - { 7.19e-1, 1.046436253332092, 2.924692924222563e-04, 1.327748933495680 }, - { 7.20e-1, 1.046449665635920, 2.932620306364937e-04, 1.327724664961610 }, - { 7.21e-1, 1.046462990264283, 2.950087154876983e-04, 1.327700442247290 }, - { 7.22e-1, 1.046476396774112, 2.975225220322731e-04, 1.327676265062370 }, - { 7.23e-1, 1.046489793136021, 2.975577990728036e-04, 1.327652133118790 }, - { 7.24e-1, 1.046503047673474, 2.991230299716263e-04, 1.327628046130760 }, - { 7.25e-1, 1.046516479213546, 3.008831249119806e-04, 1.327604003814760 }, - { 7.26e-1, 1.046530020906885, 3.030308005065120e-04, 1.327580005889490 }, - { 7.27e-1, 1.046543461106678, 3.026849499014217e-04, 1.327556052075900 }, - { 7.28e-1, 1.046556651803606, 3.046463468257975e-04, 1.327532142097090 }, - { 7.29e-1, 1.046569517786500, 3.050703016468484e-04, 1.327508275678380 }, - { 7.30e-1, 1.046581771727620, 3.056872461229752e-04, 1.327484452547220 }, - { 7.31e-1, 1.046593741831144, 3.053384705679599e-04, 1.327460672433210 }, - { 7.32e-1, 1.046605199691983, 3.073097980091759e-04, 1.327436935068060 }, - { 7.33e-1, 1.046616621102129, 3.085102646589948e-04, 1.327413240185570 }, - { 7.34e-1, 1.046628017565891, 3.110712691868297e-04, 1.327389587521650 }, - { 7.35e-1, 1.046639910218690, 3.138324950843339e-04, 1.327365976814230 }, - { 7.36e-1, 1.046651943151609, 3.173785201096158e-04, 1.327342407803330 }, - { 7.37e-1, 1.046664212729765, 3.187896703635894e-04, 1.327318880230940 }, - { 7.38e-1, 1.046676007554394, 3.196181235299881e-04, 1.327295393841100 }, - { 7.39e-1, 1.046687437975362, 3.202522588413494e-04, 1.327271948379810 }, - { 7.40e-1, 1.046697279881755, 3.202999472436175e-04, 1.327248543595060 }, - { 7.41e-1, 1.046706135060215, 3.189750555484093e-04, 1.327225179236760 }, - { 7.41e-1, 1.046713862580072, 3.211769207969832e-04, 1.327201855056800 }, - { 7.43e-1, 1.046721871614868, 3.271164739876647e-04, 1.327178570808950 }, - { 7.44e-1, 1.046730272175331, 3.346449073390914e-04, 1.327155326248880 }, - { 7.45e-1, 1.046740501259060, 3.421926333679019e-04, 1.327132121134170 }, - { 7.46e-1, 1.046751680426605, 3.497596524765032e-04, 1.327108955224250 }, - { 7.47e-1, 1.046763868052194, 3.573459650688986e-04, 1.327085828280380 }, - { 7.48e-1, 1.046777912415171, 3.649515715506583e-04, 1.327062740065700 }, - { 7.49e-1, 1.046792102242083, 3.707938097818990e-04, 1.327039690345110 }, - { 7.50e-1, 1.046804357080523, 3.631632982600370e-04, 1.327016678885350 }, - { 7.51e-1, 1.046815572343059, 3.688176560008942e-04, 1.326993705454940 }, - { 7.52e-1, 1.046825449654038, 3.744859510124802e-04, 1.326970769824140 }, - { 7.53e-1, 1.046836760483276, 3.801681835983927e-04, 1.326947871765000 }, - { 7.54e-1, 1.046848587411650, 3.858643540632899e-04, 1.326925011051280 }, - { 7.55e-1, 1.046860864366256, 3.889786095689200e-04, 1.326902187458480 }, - { 7.56e-1, 1.046872894123829, 3.920998378577586e-04, 1.326879400763780 }, - { 7.57e-1, 1.046884411012915, 3.942269549022739e-04, 1.326856650746090 }, - { 7.58e-1, 1.046895701267821, 3.965588506254991e-04, 1.326833937185960 }, - { 7.59e-1, 1.046906835569215, 3.988955762257961e-04, 1.326811259865620 }, - { 7.60e-1, 1.046917101288162, 4.010361112175749e-04, 1.326788618568950 }, - { 7.61e-1, 1.046926788333640, 3.981487270334402e-04, 1.326766013081440 }, - { 7.62e-1, 1.046933540127413, 3.952522405187550e-04, 1.326743443190220 }, - { 7.63e-1, 1.046938628178191, 3.921448270953880e-04, 1.326720908684020 }, - { 7.64e-1, 1.046941069178355, 3.910486986268732e-04, 1.326698409353150 }, - { 7.65e-1, 1.046941607565751, 3.915671726970346e-04, 1.326675944989510 }, - { 7.66e-1, 1.046941050594559, 3.955303342174100e-04, 1.326653515386550 }, - { 7.67e-1, 1.046939652476820, 4.031547449986886e-04, 1.326631120339270 }, - { 7.68e-1, 1.046938078300245, 4.116105767484131e-04, 1.326608759644200 }, - { 7.69e-1, 1.046935899260892, 4.200873169509417e-04, 1.326586433099410 }, - { 7.70e-1, 1.046933800541355, 4.330663678085595e-04, 1.326564140504460 }, - { 7.71e-1, 1.046931842347250, 4.495455754419827e-04, 1.326541881660400 }, - { 7.72e-1, 1.046931171415364, 4.695386071357624e-04, 1.326519656369770 }, - { 7.73e-1, 1.046931603412790, 4.914230986127136e-04, 1.326497464436590 }, - { 7.74e-1, 1.046935057936450, 5.133633386726560e-04, 1.326475305666330 }, - { 7.75e-1, 1.046939554268154, 5.425357270477449e-04, 1.326453179865870 }, - { 7.76e-1, 1.046948930406268, 5.699348540678753e-04, 1.326431086843580 }, - { 7.77e-1, 1.046959325716217, 5.974036690350955e-04, 1.326409026409200 }, - { 7.78e-1, 1.046974508683986, 6.249421734908905e-04, 1.326386998373890 }, - { 7.79e-1, 1.046991148046091, 6.527564809111620e-04, 1.326365002550200 }, - { 7.80e-1, 1.047011022384206, 6.794027373289624e-04, 1.326343038752090 }, - { 7.81e-1, 1.047030804279599, 6.933038356028019e-04, 1.326321106794840 }, - { 7.82e-1, 1.047050489959484, 7.070318001879980e-04, 1.326299206495120 }, - { 7.83e-1, 1.047069563350390, 7.253510589942010e-04, 1.326277337670950 }, - { 7.84e-1, 1.047088291011400, 7.435079064039188e-04, 1.326255500141660 }, - { 7.85e-1, 1.047106424922806, 7.567239792560870e-04, 1.326233693727920 }, - { 7.86e-1, 1.047124393007007, 7.695557173642177e-04, 1.326211918251710 }, - { 7.87e-1, 1.047141029801092, 7.822097693438653e-04, 1.326190173536300 }, - { 7.88e-1, 1.047157976621515, 8.038603920522829e-04, 1.326168459406250 }, - { 7.89e-1, 1.047179224925566, 8.278608148273258e-04, 1.326146775687410 }, - { 7.90e-1, 1.047201127980211, 8.517111621615959e-04, 1.326125122206880 }, - { 7.91e-1, 1.047225524410025, 8.655722047617181e-04, 1.326103498793040 }, - { 7.92e-1, 1.047248608976219, 8.725492888296510e-04, 1.326081905275470 }, - { 7.93e-1, 1.047269117642915, 8.797512649178575e-04, 1.326060341485040 }, - { 7.94e-1, 1.047289463317727, 8.924322472499099e-04, 1.326038807253790 }, - { 7.95e-1, 1.047309613452492, 9.112443464370387e-04, 1.326017302415020 }, - { 7.96e-1, 1.047331668415949, 9.296801609789082e-04, 1.325995826803190 }, - { 7.97e-1, 1.047355299038291, 9.452068444725653e-04, 1.325974380253980 }, - { 7.98e-1, 1.047379136045100, 9.552786828842519e-04, 1.325952962604240 }, - { 7.99e-1, 1.047403072830967, 9.657960036606947e-04, 1.325931573691990 }, - { 8.00e-1, 1.047427193866072, 9.790895649030144e-04, 1.325910213356430 }, - { 8.01e-1, 1.047451361506094, 9.924137067527394e-04, 1.325888881437880 }, - { 8.02e-1, 1.047477610041293, 1.005556152948590e-03, 1.325867577777830 }, - { 8.03e-1, 1.047504359325982, 1.018516099946340e-03, 1.325846302218880 }, - { 8.04e-1, 1.047536241473395, 1.031292744178663e-03, 1.325825054604780 }, - { 8.05e-1, 1.047568907885968, 1.048146908027978e-03, 1.325803834780350 }, - { 8.06e-1, 1.047601530876990, 1.040078355092751e-03, 1.325782642591560 }, - { 8.07e-1, 1.047634095398597, 1.034976873960194e-03, 1.325761477885440 }, - { 8.08e-1, 1.047660247782939, 1.027292698511066e-03, 1.325740340510110 }, - { 8.09e-1, 1.047686099812982, 1.019800174887994e-03, 1.325719230314780 }, - { 8.10e-1, 1.047707779443527, 1.012071263221199e-03, 1.325698147149700 }, - { 8.11e-1, 1.047729237681430, 9.995966685213660e-04, 1.325677090866210 }, - { 8.12e-1, 1.047743536155764, 9.827886098305182e-04, 1.325656061316680 }, - { 8.13e-1, 1.047757488790403, 9.661507499608610e-04, 1.325635058354500 }, - { 8.14e-1, 1.047763804444236, 9.490374090644322e-04, 1.325614081834110 }, - { 8.15e-1, 1.047769969728286, 9.448242560608527e-04, 1.325593131610980 }, - { 8.16e-1, 1.047773137187235, 9.401656527267664e-04, 1.325572207541580 }, - { 8.17e-1, 1.047776216953897, 9.354925875861295e-04, 1.325551309483370 }, - { 8.18e-1, 1.047777270288282, 9.308050603294224e-04, 1.325530437294830 }, - { 8.19e-1, 1.047778186624244, 9.269704113819045e-04, 1.325509590835410 }, - { 8.20e-1, 1.047774802403966, 9.233405483764064e-04, 1.325488769965540 }, - { 8.21e-1, 1.047771400536673, 9.194815325215046e-04, 1.325467974546630 }, - { 8.22e-1, 1.047764370361179, 9.177865956932137e-04, 1.325447204441030 }, - { 8.23e-1, 1.047757377048223, 9.256726507307197e-04, 1.325426459512080 }, - { 8.24e-1, 1.047749487584204, 9.337935272091834e-04, 1.325405739624020 }, - { 8.25e-1, 1.047741618195193, 9.417131416291529e-04, 1.325385044642050 }, - { 8.26e-1, 1.047731838654609, 9.498681147314757e-04, 1.325364374432310 }, - { 8.27e-1, 1.047721973752322, 9.578212902678193e-04, 1.325343728861850 }, - { 8.27e-1, 1.047707491516565, 9.668873570831467e-04, 1.325323107798620 }, - { 8.29e-1, 1.047692462294499, 9.779484068345122e-04, 1.325302511111490 }, - { 8.30e-1, 1.047673226736550, 9.929897488926683e-04, 1.325281938670230 }, - { 8.31e-1, 1.047653141035367, 1.015986829694636e-03, 1.325261390345490 }, - { 8.32e-1, 1.047633752137463, 1.044104435780425e-03, 1.325240866008810 }, - { 8.33e-1, 1.047614638723221, 1.085522902123505e-03, 1.325220365532600 }, - { 8.34e-1, 1.047600760948479, 1.132339618517545e-03, 1.325199888790140 }, - { 8.35e-1, 1.047589518185891, 1.182141510046346e-03, 1.325179435655570 }, - { 8.36e-1, 1.047582936855137, 1.232060910894234e-03, 1.325159006003890 }, - { 8.37e-1, 1.047579827948053, 1.281876161129743e-03, 1.325138599710930 }, - { 8.38e-1, 1.047578951020551, 1.332030321081121e-03, 1.325118216653370 }, - { 8.39e-1, 1.047580594183633, 1.382301999067603e-03, 1.325097856708720 }, - { 8.40e-1, 1.047583955198737, 1.432691198007155e-03, 1.325077519755300 }, - { 8.41e-1, 1.047590436261947, 1.488543580847245e-03, 1.325057205672260 }, - { 8.42e-1, 1.047598319842640, 1.544526367039972e-03, 1.325036914339570 }, - { 8.43e-1, 1.047610935731388, 1.600639559836751e-03, 1.325016645637990 }, - { 8.44e-1, 1.047624127180442, 1.656883162494959e-03, 1.324996399449070 }, - { 8.45e-1, 1.047643117657286, 1.713257178278004e-03, 1.324976175655160 }, - { 8.46e-1, 1.047662084005794, 1.769761610455194e-03, 1.324955974139380 }, - { 8.47e-1, 1.047681999317462, 1.813160079636676e-03, 1.324935794785660 }, - { 8.48e-1, 1.047701946094193, 1.859801906390260e-03, 1.324915637478650 }, - { 8.49e-1, 1.047715295897636, 1.906550012807259e-03, 1.324895502103800 }, - { 8.50e-1, 1.047726844952281, 1.945298989124846e-03, 1.324875388547310 }, - { 8.51e-1, 1.047738549879133, 1.989770398984366e-03, 1.324855296696110 }, - { 8.52e-1, 1.047750348897479, 2.095729044907462e-03, 1.324835226437890 }, - { 8.53e-1, 1.047767605294426, 2.133017683744185e-03, 1.324815177661090 }, - { 8.54e-1, 1.047791905510931, 2.218348403502016e-03, 1.324795150254850 }, - { 8.55e-1, 1.047818418882227, 2.310217348244026e-03, 1.324775144109060 }, - { 8.56e-1, 1.047851843468609, 2.395721729565107e-03, 1.324755159114300 }, - { 8.57e-1, 1.047885588632721, 2.458038095689459e-03, 1.324735195161900 }, - { 8.58e-1, 1.047928039204696, 2.551634788315200e-03, 1.324715252143870 }, - { 8.58e-1, 1.047970412470177, 2.624737227819212e-03, 1.324695329952930 }, - { 8.60e-1, 1.048023220548083, 2.723522375062108e-03, 1.324675428482480 }, - { 8.61e-1, 1.048078126196254, 2.789002004916494e-03, 1.324655547626620 }, - { 8.62e-1, 1.048140300421103, 2.899160408289475e-03, 1.324635687280130 }, - { 8.63e-1, 1.048207241938472, 2.960412745879289e-03, 1.324615847338480 }, - { 8.64e-1, 1.048278409160125, 3.054534627917653e-03, 1.324596027697780 }, - { 8.65e-1, 1.048357052707288, 3.132368488989489e-03, 1.324576228254840 }, - { 8.66e-1, 1.048436830487660, 3.213359081793805e-03, 1.324556448907090 }, - { 8.67e-1, 1.048526428644327, 3.259154960386772e-03, 1.324536689552640 }, - { 8.68e-1, 1.048615843313936, 3.357024712042136e-03, 1.324516950090250 }, - { 8.69e-1, 1.048715148218825, 3.489421142093250e-03, 1.324497230419300 }, - { 8.70e-1, 1.048815971181798, 3.531293371955145e-03, 1.324477530439830 }, - { 8.71e-1, 1.048926857924925, 3.564945421308734e-03, 1.324457850052500 }, - { 8.72e-1, 1.049044908795749, 3.679304100653368e-03, 1.324438189158590 }, - { 8.73e-1, 1.049166413128009, 3.726604082701499e-03, 1.324418547660020 }, - { 8.74e-1, 1.049296428879817, 3.770298531458578e-03, 1.324398925459290 }, - { 8.75e-1, 1.049426171616863, 3.818953357363171e-03, 1.324379322459550 }, - { 8.76e-1, 1.049564210263465, 3.868175047738017e-03, 1.324359738564520 }, - { 8.77e-1, 1.049701958513664, 3.907274996952621e-03, 1.324340173678550 }, - { 8.78e-1, 1.049845994882515, 3.934123795458572e-03, 1.324320627706550 }, - { 8.79e-1, 1.049992662084591, 3.955200886907125e-03, 1.324301100554040 }, - { 8.80e-1, 1.050141519573849, 3.967221321383652e-03, 1.324281592127130 }, - { 8.81e-1, 1.050294307337905, 3.988597805776602e-03, 1.324262102332480 }, - { 8.82e-1, 1.050447085007972, 3.978224625708395e-03, 1.324242631077360 }, - { 8.83e-1, 1.050607609986077, 3.968051170294418e-03, 1.324223178269580 }, - { 8.84e-1, 1.050767796843387, 3.955969767958174e-03, 1.324203743817520 }, - { 8.85e-1, 1.050923092122999, 3.927901360070584e-03, 1.324184327630120 }, - { 8.86e-1, 1.051076031953918, 3.841759824159816e-03, 1.324164929616880 }, - { 8.87e-1, 1.051227183758931, 3.798665924324868e-03, 1.324145549687840 }, - { 8.88e-1, 1.051375087152381, 3.745578981754750e-03, 1.324126187753590 }, - { 8.89e-1, 1.051522682487672, 3.701078944638891e-03, 1.324106843725250 }, - { 8.89e-1, 1.051666311265850, 3.661657321875855e-03, 1.324087517514480 }, - { 8.91e-1, 1.051809500042060, 3.571127079571336e-03, 1.324068209033480 }, - { 8.92e-1, 1.051943498022589, 3.351531098795739e-03, 1.324048918194950 }, - { 8.93e-1, 1.052070762485258, 3.324330688851290e-03, 1.324029644912130 }, - { 8.94e-1, 1.052194726230713, 3.208200037721505e-03, 1.324010389098780 }, - { 8.95e-1, 1.052304860963421, 3.056687924568558e-03, 1.323991150669150 }, - { 8.96e-1, 1.052414773434226, 2.942349178462081e-03, 1.323971929538010 }, - { 8.97e-1, 1.052502780022198, 2.834877091402793e-03, 1.323952725620640 }, - { 8.98e-1, 1.052584273861886, 2.666932219079857e-03, 1.323933538832790 }, - { 8.99e-1, 1.052659626877916, 2.501460898299947e-03, 1.323914369090740 }, - { 9.00e-1, 1.052723802933926, 2.379029853183195e-03, 1.323895216311230 }, - { 9.01e-1, 1.052787858771010, 2.261571369321692e-03, 1.323876080411500 }, - { 9.02e-1, 1.052821869541360, 2.080000935934572e-03, 1.323856961309270 }, - { 9.03e-1, 1.052853084720993, 1.890117567359862e-03, 1.323837858922710 }, - { 9.04e-1, 1.052866397922094, 1.724726050468831e-03, 1.323818773170500 }, - { 9.05e-1, 1.052859379171112, 1.640289877936772e-03, 1.323799703971770 }, - { 9.06e-1, 1.052852413253518, 1.487688805686849e-03, 1.323780651246100 }, - { 9.07e-1, 1.052848611642345, 1.376099936169961e-03, 1.323761614913550 }, - { 9.08e-1, 1.052844838670080, 1.333344566482168e-03, 1.323742594894630 }, - { 9.09e-1, 1.052835368829740, 1.189275458028569e-03, 1.323723591110280 }, - { 9.10e-1, 1.052820869248812, 1.054290942605935e-03, 1.323704603481930 }, - { 9.11e-1, 1.052805450313970, 9.366342882581003e-04, 1.323685631931410 }, - { 9.12e-1, 1.052767449737131, 8.436181167945962e-04, 1.323666676381010 }, - { 9.13e-1, 1.052729551679777, 7.317542083422428e-04, 1.323647736753440 }, - { 9.14e-1, 1.052687896630913, 6.763454072570361e-04, 1.323628812971860 }, - { 9.15e-1, 1.052643047484057, 5.652567173572225e-04, 1.323609904959850 }, - { 9.16e-1, 1.052598290617709, 5.129373259408466e-04, 1.323591012641400 }, - { 9.17e-1, 1.052552554797870, 4.583126029690992e-04, 1.323572135940930 }, - { 9.18e-1, 1.052506937663857, 4.038093454985267e-04, 1.323553274783280 }, - { 9.19e-1, 1.052462313531449, 3.586879638784663e-04, 1.323534429093690 }, - { 9.20e-1, 1.052418708673314, 3.120024706391773e-04, 1.323515598797820 }, - { 9.20e-1, 1.052375217601570, 2.700975173655584e-04, 1.323496783821710 }, - { 9.22e-1, 1.052334018043352, 2.391012133636555e-04, 1.323477984091820 }, - { 9.23e-1, 1.052293047791844, 2.095046342505333e-04, 1.323459199535020 }, - { 9.24e-1, 1.052254431963422, 1.862129546266315e-04, 1.323440430078530 }, - { 9.25e-1, 1.052219445428962, 1.636055562032376e-04, 1.323421675650000 }, - { 9.26e-1, 1.052184553785031, 1.480692683441966e-04, 1.323402936177440 }, - { 9.27e-1, 1.052153554316969, 1.324985958014863e-04, 1.323384211589260 }, - { 9.28e-1, 1.052123495606230, 1.168935376472338e-04, 1.323365501814220 }, - { 9.29e-1, 1.052094536477733, 1.010077326289579e-04, 1.323346806781490 }, - { 9.30e-1, 1.052069073221142, 9.248587260842015e-05, 1.323328126420590 }, - { 9.31e-1, 1.052043684243811, 8.468590033687339e-05, 1.323309460661400 }, - { 9.32e-1, 1.052021364441905, 7.686873437855818e-05, 1.323290809434190 }, - { 9.33e-1, 1.052000888614132, 6.903437426761178e-05, 1.323272172669560 }, - { 9.34e-1, 1.051980527802691, 6.440296793432885e-05, 1.323253550298500 }, - { 9.35e-1, 1.051964754164010, 6.124918707549857e-05, 1.323234942252320 }, - { 9.36e-1, 1.051949034063168, 5.808842130256970e-05, 1.323216348462700 }, - { 9.37e-1, 1.051934790474261, 5.492067042606937e-05, 1.323197768861670 }, - { 9.38e-1, 1.051922913853669, 5.199471278637858e-05, 1.323179203381600 }, - { 9.39e-1, 1.051911082461253, 4.632278740626288e-05, 1.323160651955180 }, - { 9.40e-1, 1.051900307289650, 3.839465831134583e-05, 1.323142114515460 }, - { 9.41e-1, 1.051889940149579, 3.369393273225441e-05, 1.323123590995830 }, - { 9.42e-1, 1.051879766217186, 3.148153037298949e-05, 1.323105081329980 }, - { 9.43e-1, 1.051871533893133, 3.051490221754005e-05, 1.323086585451950 }, - { 9.44e-1, 1.051863339025340, 2.929573388502167e-05, 1.323068103296110 }, - { 9.45e-1, 1.051856470570902, 2.882585761204386e-05, 1.323049634797120 }, - { 9.46e-1, 1.051851586647983, 2.835490692158482e-05, 1.323031179890000 }, - { 9.47e-1, 1.051846733136190, 2.763168465136746e-05, 1.323012738510040 }, - { 9.48e-1, 1.051844453224846, 2.715831628037134e-05, 1.322994310592870 }, - { 9.49e-1, 1.051843253457583, 2.693560805431627e-05, 1.322975896074430 }, - { 9.50e-1, 1.051842076403402, 2.671236282563884e-05, 1.322957494890950 } -}; - -static const size_t nproperties = - sizeof(properties)/sizeof(struct optical_properties); - struct geometry { float* vertices; unsigned* indices; @@ -608,6 +56,7 @@ struct sampler_context { struct geometry geometry; double aspect_ratio; double mean_radius; + double sigma; }; struct cylinder { @@ -622,33 +71,11 @@ get_material_property const double wavelength, struct sschiff_material_properties* props) { - double n_medium, kappa_r, n_r; - size_t i; (void)mtl; - - FOR_EACH(i, 0, nproperties) { - if(properties[i].wavelength >= wavelength) - break; - } - CHECK(i < nproperties, 1); - - if(eq_eps(properties[i].wavelength, wavelength, 1.e-12)) { - n_medium = properties[i].n_medium; - kappa_r = properties[i].kappa_r; - n_r = properties[i].n_r; - } else { - const size_t i_prev = i - 1; - double d = properties[i].wavelength - properties[i_prev].wavelength; - double u = (wavelength - properties[i_prev].wavelength) / d; - n_medium = u * properties[i].n_medium + (1-u) * properties[i_prev].n_medium; - kappa_r = u * properties[i].kappa_r + (1-u) * properties[i_prev].kappa_r; - n_r = u * properties[i].n_r + (1-u) * properties[i_prev].n_r; - } - - NCHECK(props, NULL); - props->medium_refractive_index = n_medium; - props->relative_imaginary_refractive_index = kappa_r; - props->relative_real_refractive_index = n_r; + CHECK(eq_eps(wavelength, 0.6, 1.e-6), 1); + props->medium_refractive_index = 4.0/3.0; + props->relative_imaginary_refractive_index = 4.e-3; + props->relative_real_refractive_index = 1.01; } static void @@ -713,7 +140,7 @@ sample_cylinder size_t nverts, nprims; (void)rng, (void)sampler_context; - sample = ssp_ran_lognormal(rng, log(sampler_ctx->mean_radius), log(1.1374)); + sample = ssp_ran_lognormal(rng, log(sampler_ctx->mean_radius), log(sampler_ctx->sigma)); cylinder.geometry = &sampler_ctx->geometry; cylinder.radius = (float)(sample / pow(3.0 / (2.0*sampler_ctx->aspect_ratio), 1.0/3.0)); cylinder.height = (float)(2.f * cylinder.radius / sampler_ctx->aspect_ratio); @@ -799,6 +226,19 @@ cylinder_release(struct geometry* geometry) geometry->indices = NULL; } +static INLINE void +compute_intersection + (double intersection[2], double E0, double SE0, double E1, double SE1) +{ + double interval0[2], interval1[2]; + interval0[0] = E0 - 4*SE0; + interval0[1] = E0 + 4*SE0; + interval1[0] = E1 - 4*SE1; + interval1[1] = E1 + 4*SE1; + intersection[0] = MMAX(interval0[0], interval1[0]); + intersection[1] = MMIN(interval0[1], interval1[1]); +} + int main(int argc, char** argv) { @@ -810,28 +250,318 @@ main(int argc, char** argv) struct sschiff_estimator* estimator; struct ssp_rng* rng; struct time t0, t1; - struct result results[] = { - { 2.76750904574817, 2.34364058801497e-3 }, - { 2.50099694990826, 2.20010054373943e-3 }, - { 2.32782435081181, 2.05630225764867e-3 }, - { 1.96661661205101, 1.82091931178789e-3 }, - { 1.71149696366194, 1.63845385406345e-3 }, - { 1.54229472112306, 1.49426946826458e-3 }, - { 1.39877652964399, 1.36729344819050e-3 } + const struct result results[] = { + { 5.29715E-19, 4.28609E-22, 1.71856E-25, 2.19472E-28, 5.29715E-19, 4.28609E-22, 3.81205E-12, 4.22555E-16 }, + { 1.13812E-04, 9.18322E-08, 2.21384E-06, 2.82025E-09, 1.16026E-04, 9.45526E-08, 1.37239E-02, 1.52125E-06 }, + { 9.05619E-04, 7.28692E-07, 3.52032E-05, 4.47348E-08, 9.40822E-04, 7.71882E-07, 5.48945E-02, 6.08489E-06 }, + { 3.04018E-03, 2.43947E-06, 1.77116E-04, 2.24508E-07, 3.21730E-03, 2.65641E-06, 1.23512E-01, 1.36909E-05 }, + { 7.16810E-03, 5.73592E-06, 5.56304E-04, 7.03373E-07, 7.72440E-03, 6.41615E-06, 2.19576E-01, 2.43394E-05 }, + { 1.39261E-02, 1.11131E-05, 1.34971E-03, 1.70216E-06, 1.52758E-02, 1.27606E-05, 3.43087E-01, 3.80302E-05 }, + { 2.39373E-02, 1.90499E-05, 2.78123E-03, 3.49844E-06, 2.67186E-02, 2.24385E-05, 4.94045E-01, 5.47634E-05 }, + { 3.78116E-02, 3.00096E-05, 5.12017E-03, 6.42370E-06, 4.29318E-02, 3.62361E-05, 6.72450E-01, 7.45390E-05 }, + { 5.61457E-02, 4.44398E-05, 8.67959E-03, 1.08606E-05, 6.48252E-02, 5.49743E-05, 8.78301E-01, 9.73571E-05 }, + { 7.95235E-02, 6.27738E-05, 1.38147E-02, 1.72400E-05, 9.33382E-02, 7.95073E-05, 1.11160E+00, 1.23217E-04 }, + { 1.08517E-01, 8.54299E-05, 2.09215E-02, 2.60386E-05, 1.29438E-01, 1.10720E-04, 1.37234E+00, 1.52120E-04 }, + { 1.43684E-01, 1.12813E-04, 3.04347E-02, 3.77758E-05, 1.74119E-01, 1.49525E-04, 1.66054E+00, 1.84065E-04 }, + { 1.85573E-01, 1.45313E-04, 4.28269E-02, 5.30113E-05, 2.28400E-01, 1.96862E-04, 1.97617E+00, 2.19053E-04 }, + { 2.34719E-01, 1.83309E-04, 5.86064E-02, 7.23425E-05, 2.93325E-01, 2.53696E-04, 2.31926E+00, 2.57083E-04 }, + { 2.91645E-01, 2.27165E-04, 7.83160E-02, 9.64016E-05, 3.69961E-01, 3.21011E-04, 2.68979E+00, 2.98155E-04 }, + { 3.56865E-01, 2.77233E-04, 1.02531E-01, 1.25854E-04, 4.59396E-01, 3.99816E-04, 3.08777E+00, 3.42270E-04 }, + { 4.30879E-01, 3.33853E-04, 1.31860E-01, 1.61393E-04, 5.62739E-01, 4.91132E-04, 3.51320E+00, 3.89427E-04 }, + { 5.14178E-01, 3.97355E-04, 1.66938E-01, 2.03742E-04, 6.81116E-01, 5.96001E-04, 3.96607E+00, 4.39627E-04 }, + { 6.07243E-01, 4.68053E-04, 2.08431E-01, 2.53648E-04, 8.15674E-01, 7.15476E-04, 4.44639E+00, 4.92869E-04 }, + { 7.10542E-01, 5.46255E-04, 2.57033E-01, 3.11882E-04, 9.67574E-01, 8.50622E-04, 4.95415E+00, 5.49153E-04 }, + { 8.24534E-01, 6.32256E-04, 3.13460E-01, 3.79232E-04, 1.13799E+00, 1.00252E-03, 5.48937E+00, 6.08480E-04 }, + { 9.49670E-01, 7.26339E-04, 3.78455E-01, 4.56509E-04, 1.32813E+00, 1.17224E-03, 6.05203E+00, 6.70849E-04 }, + { 1.08639E+00, 8.28779E-04, 4.52783E-01, 5.44535E-04, 1.53917E+00, 1.36088E-03, 6.64213E+00, 7.36261E-04 }, + { 1.23512E+00, 9.39839E-04, 5.37230E-01, 6.44150E-04, 1.77235E+00, 1.56954E-03, 7.25969E+00, 8.04715E-04 }, + { 1.39628E+00, 1.05978E-03, 6.32603E-01, 7.56203E-04, 2.02888E+00, 1.79931E-03, 7.90469E+00, 8.76211E-04 }, + { 1.57029E+00, 1.18883E-03, 7.39725E-01, 8.81552E-04, 2.31001E+00, 2.05129E-03, 8.57713E+00, 9.50750E-04 }, + { 1.75754E+00, 1.32725E-03, 8.59440E-01, 1.02107E-03, 2.61698E+00, 2.32657E-03, 9.27703E+00, 1.02833E-03 }, + { 1.95844E+00, 1.47525E-03, 9.92605E-01, 1.17562E-03, 2.95104E+00, 2.62624E-03, 1.00044E+01, 1.10895E-03 }, + { 2.17335E+00, 1.63305E-03, 1.14009E+00, 1.34608E-03, 3.31344E+00, 2.95141E-03, 1.07592E+01, 1.19262E-03 }, + { 2.40267E+00, 1.80087E-03, 1.30279E+00, 1.53333E-03, 3.70546E+00, 3.30314E-03, 1.15414E+01, 1.27933E-03 }, + { 2.64676E+00, 1.97891E-03, 1.48159E+00, 1.73825E-03, 4.12835E+00, 3.68251E-03, 1.23511E+01, 1.36908E-03 }, + { 2.90597E+00, 2.16735E-03, 1.67741E+00, 1.96171E-03, 4.58338E+00, 4.09060E-03, 1.31882E+01, 1.46187E-03 }, + { 3.18066E+00, 2.36640E-03, 1.89117E+00, 2.20460E-03, 5.07183E+00, 4.52845E-03, 1.40528E+01, 1.55771E-03 }, + { 3.47118E+00, 2.57622E-03, 2.12379E+00, 2.46776E-03, 5.59497E+00, 4.99711E-03, 1.49448E+01, 1.65659E-03 }, + { 3.77786E+00, 2.79699E-03, 2.37621E+00, 2.75208E-03, 6.15407E+00, 5.49762E-03, 1.58643E+01, 1.75851E-03 }, + { 4.10102E+00, 3.02889E-03, 2.64937E+00, 3.05839E-03, 6.75039E+00, 6.03099E-03, 1.68112E+01, 1.86347E-03 }, + { 4.44100E+00, 3.27205E-03, 2.94421E+00, 3.38755E-03, 7.38521E+00, 6.59822E-03, 1.77855E+01, 1.97147E-03 }, + { 4.79810E+00, 3.52665E-03, 3.26169E+00, 3.74039E-03, 8.05979E+00, 7.20030E-03, 1.87873E+01, 2.08252E-03 }, + { 5.17263E+00, 3.79282E-03, 3.60277E+00, 4.11772E-03, 8.77539E+00, 7.83819E-03, 1.98166E+01, 2.19661E-03 }, + { 5.56489E+00, 4.07070E-03, 3.96838E+00, 4.52037E-03, 9.53327E+00, 8.51284E-03, 2.08733E+01, 2.31374E-03 }, + { 5.97518E+00, 4.36043E-03, 4.35950E+00, 4.94911E-03, 1.03347E+01, 9.22518E-03, 2.19575E+01, 2.43392E-03 }, + { 6.40378E+00, 4.66214E-03, 4.77708E+00, 5.40474E-03, 1.11809E+01, 9.97612E-03, 2.30690E+01, 2.55714E-03 }, + { 6.85097E+00, 4.97595E-03, 5.22207E+00, 5.88800E-03, 1.20730E+01, 1.07665E-02, 2.42081E+01, 2.68339E-03 }, + { 7.31703E+00, 5.30197E-03, 5.69542E+00, 6.39965E-03, 1.30125E+01, 1.15973E-02, 2.53746E+01, 2.81270E-03 }, + { 7.80222E+00, 5.64032E-03, 6.19810E+00, 6.94041E-03, 1.40003E+01, 1.24693E-02, 2.65685E+01, 2.94504E-03 }, + { 8.30681E+00, 5.99110E-03, 6.73103E+00, 7.51098E-03, 1.50378E+01, 1.33832E-02, 2.77899E+01, 3.08043E-03 }, + { 8.83105E+00, 6.35442E-03, 7.29517E+00, 8.11205E-03, 1.61262E+01, 1.43399E-02, 2.90387E+01, 3.21886E-03 }, + { 9.37519E+00, 6.73037E-03, 7.89146E+00, 8.74428E-03, 1.72666E+01, 1.53402E-02, 3.03150E+01, 3.36033E-03 }, + { 9.93948E+00, 7.11905E-03, 8.52081E+00, 9.40831E-03, 1.84603E+01, 1.63848E-02, 3.16187E+01, 3.50484E-03 }, + { 1.05241E+01, 7.52055E-03, 9.18416E+00, 1.01048E-02, 1.97083E+01, 1.74744E-02, 3.29499E+01, 3.65240E-03 }, + { 1.11294E+01, 7.93493E-03, 9.88241E+00, 1.08342E-02, 2.10118E+01, 1.86096E-02, 3.43085E+01, 3.80300E-03 }, + { 1.17556E+01, 8.36230E-03, 1.06165E+01, 1.15973E-02, 2.23720E+01, 1.97912E-02, 3.56946E+01, 3.95664E-03 }, + { 1.24028E+01, 8.80272E-03, 1.13873E+01, 1.23945E-02, 2.37900E+01, 2.10198E-02, 3.71081E+01, 4.11332E-03 }, + { 1.30712E+01, 9.25625E-03, 1.21956E+01, 1.32263E-02, 2.52669E+01, 2.22959E-02, 3.85490E+01, 4.27305E-03 }, + { 1.37612E+01, 9.72298E-03, 1.30425E+01, 1.40933E-02, 2.68037E+01, 2.36201E-02, 4.00174E+01, 4.43582E-03 }, + { 1.44729E+01, 1.02030E-02, 1.39287E+01, 1.49959E-02, 2.84016E+01, 2.49931E-02, 4.15133E+01, 4.60163E-03 }, + { 1.52065E+01, 1.06963E-02, 1.48550E+01, 1.59346E-02, 3.00616E+01, 2.64152E-02, 4.30366E+01, 4.77048E-03 }, + { 1.59622E+01, 1.12029E-02, 1.58225E+01, 1.69097E-02, 3.17846E+01, 2.78869E-02, 4.45873E+01, 4.94237E-03 }, + { 1.67402E+01, 1.17230E-02, 1.68317E+01, 1.79217E-02, 3.35719E+01, 2.94088E-02, 4.61655E+01, 5.11731E-03 }, + { 1.75407E+01, 1.22566E-02, 1.78836E+01, 1.89710E-02, 3.54243E+01, 3.09811E-02, 4.77712E+01, 5.29529E-03 }, + { 1.83638E+01, 1.28036E-02, 1.89790E+01, 2.00578E-02, 3.73428E+01, 3.26045E-02, 4.94043E+01, 5.47631E-03 }, + { 1.92098E+01, 1.33643E-02, 2.01186E+01, 2.11824E-02, 3.93284E+01, 3.42790E-02, 5.10648E+01, 5.66038E-03 }, + { 2.00788E+01, 1.39385E-02, 2.13033E+01, 2.23453E-02, 4.13821E+01, 3.60052E-02, 5.27528E+01, 5.84749E-03 }, + { 2.09710E+01, 1.45264E-02, 2.25337E+01, 2.35465E-02, 4.35047E+01, 3.77833E-02, 5.44682E+01, 6.03764E-03 }, + { 2.18866E+01, 1.51280E-02, 2.38106E+01, 2.47864E-02, 4.56972E+01, 3.96136E-02, 5.62111E+01, 6.23083E-03 }, + { 2.28257E+01, 1.57432E-02, 2.51347E+01, 2.60651E-02, 4.79604E+01, 4.14962E-02, 5.79814E+01, 6.42706E-03 }, + { 2.37885E+01, 1.63722E-02, 2.65067E+01, 2.73829E-02, 5.02952E+01, 4.34315E-02, 5.97791E+01, 6.62634E-03 }, + { 2.47751E+01, 1.70150E-02, 2.79274E+01, 2.87397E-02, 5.27025E+01, 4.54196E-02, 6.16044E+01, 6.82866E-03 }, + { 2.57857E+01, 1.76715E-02, 2.93973E+01, 3.01359E-02, 5.51830E+01, 4.74606E-02, 6.34570E+01, 7.03402E-03 }, + { 2.68205E+01, 1.83419E-02, 3.09171E+01, 3.15714E-02, 5.77375E+01, 4.95546E-02, 6.53371E+01, 7.24243E-03 }, + { 2.78795E+01, 1.90261E-02, 3.24874E+01, 3.30463E-02, 6.03669E+01, 5.17018E-02, 6.72447E+01, 7.45387E-03 }, + { 2.89630E+01, 1.97242E-02, 3.41089E+01, 3.45606E-02, 6.30719E+01, 5.39023E-02, 6.91797E+01, 7.66836E-03 }, + { 3.00710E+01, 2.04361E-02, 3.57822E+01, 3.61145E-02, 6.58531E+01, 5.61559E-02, 7.11421E+01, 7.88589E-03 }, + { 3.12037E+01, 2.11619E-02, 3.75077E+01, 3.77077E-02, 6.87114E+01, 5.84628E-02, 7.31320E+01, 8.10647E-03 }, + { 3.23612E+01, 2.19015E-02, 3.92861E+01, 3.93404E-02, 7.16474E+01, 6.08230E-02, 7.51494E+01, 8.33008E-03 }, + { 3.35437E+01, 2.26551E-02, 4.11179E+01, 4.10124E-02, 7.46617E+01, 6.32362E-02, 7.71941E+01, 8.55674E-03 }, + { 3.47513E+01, 2.34226E-02, 4.30036E+01, 4.27236E-02, 7.77550E+01, 6.57026E-02, 7.92664E+01, 8.78644E-03 }, + { 3.59841E+01, 2.42040E-02, 4.49438E+01, 4.44740E-02, 8.09278E+01, 6.82220E-02, 8.13661E+01, 9.01918E-03 }, + { 3.72422E+01, 2.49994E-02, 4.69387E+01, 4.62633E-02, 8.41809E+01, 7.07942E-02, 8.34932E+01, 9.25497E-03 }, + { 3.85257E+01, 2.58086E-02, 4.89890E+01, 4.80915E-02, 8.75147E+01, 7.34190E-02, 8.56478E+01, 9.49380E-03 }, + { 3.98347E+01, 2.66318E-02, 5.10951E+01, 4.99583E-02, 9.09298E+01, 7.60964E-02, 8.78298E+01, 9.73567E-03 }, + { 4.11695E+01, 2.74689E-02, 5.32573E+01, 5.18635E-02, 9.44268E+01, 7.88261E-02, 9.00392E+01, 9.98058E-03 }, + { 4.25299E+01, 2.83199E-02, 5.54761E+01, 5.38068E-02, 9.80060E+01, 8.16078E-02, 9.22762E+01, 1.02285E-02 }, + { 4.39163E+01, 2.91849E-02, 5.77517E+01, 5.57881E-02, 1.01668E+02, 8.44413E-02, 9.45405E+01, 1.04795E-02 }, + { 4.53286E+01, 3.00637E-02, 6.00847E+01, 5.78070E-02, 1.05413E+02, 8.73264E-02, 9.68323E+01, 1.07336E-02 }, + { 4.67669E+01, 3.09565E-02, 6.24752E+01, 5.98633E-02, 1.09242E+02, 9.02627E-02, 9.91516E+01, 1.09907E-02 }, + { 4.82314E+01, 3.18632E-02, 6.49236E+01, 6.19566E-02, 1.13155E+02, 9.32498E-02, 1.01498E+02, 1.12508E-02 }, + { 4.97222E+01, 3.27838E-02, 6.74302E+01, 6.40865E-02, 1.17152E+02, 9.62875E-02, 1.03872E+02, 1.15139E-02 }, + { 5.12393E+01, 3.37182E-02, 6.99952E+01, 6.62528E-02, 1.21235E+02, 9.93754E-02, 1.06274E+02, 1.17802E-02 }, + { 5.27829E+01, 3.46665E-02, 7.26189E+01, 6.84550E-02, 1.25402E+02, 1.02513E-01, 1.08703E+02, 1.20494E-02 }, + { 5.43529E+01, 3.56287E-02, 7.53015E+01, 7.06928E-02, 1.29654E+02, 1.05700E-01, 1.11160E+02, 1.23217E-02 }, + { 5.59496E+01, 3.66048E-02, 7.80432E+01, 7.29657E-02, 1.33993E+02, 1.08936E-01, 1.13643E+02, 1.25970E-02 }, + { 5.75729E+01, 3.75947E-02, 8.08442E+01, 7.52733E-02, 1.38417E+02, 1.12220E-01, 1.16155E+02, 1.28754E-02 }, + { 5.92230E+01, 3.85984E-02, 8.37046E+01, 7.76152E-02, 1.42928E+02, 1.15553E-01, 1.18694E+02, 1.31568E-02 }, + { 6.09000E+01, 3.96159E-02, 8.66246E+01, 7.99908E-02, 1.47525E+02, 1.18933E-01, 1.21260E+02, 1.34413E-02 }, + { 6.26039E+01, 4.06472E-02, 8.96044E+01, 8.23998E-02, 1.52208E+02, 1.22360E-01, 1.23854E+02, 1.37288E-02 }, + { 6.43347E+01, 4.16923E-02, 9.26439E+01, 8.48416E-02, 1.56979E+02, 1.25834E-01, 1.26475E+02, 1.40194E-02 }, + { 6.60926E+01, 4.27512E-02, 9.57432E+01, 8.73157E-02, 1.61836E+02, 1.29353E-01, 1.29123E+02, 1.43130E-02 }, + { 6.78777E+01, 4.38238E-02, 9.89026E+01, 8.98216E-02, 1.66780E+02, 1.32919E-01, 1.31800E+02, 1.46096E-02 }, + { 6.96899E+01, 4.49101E-02, 1.02122E+02, 9.23587E-02, 1.71812E+02, 1.36529E-01, 1.34503E+02, 1.49093E-02 }, + { 7.15294E+01, 4.60101E-02, 1.05401E+02, 9.49266E-02, 1.76931E+02, 1.40183E-01, 1.37234E+02, 1.52120E-02 }, + { 7.33962E+01, 4.71238E-02, 1.08740E+02, 9.75247E-02, 1.82137E+02, 1.43881E-01, 1.39992E+02, 1.55177E-02 }, + { 7.52904E+01, 4.82512E-02, 1.12140E+02, 1.00152E-01, 1.87430E+02, 1.47622E-01, 1.42778E+02, 1.58265E-02 }, + { 7.72121E+01, 4.93922E-02, 1.15599E+02, 1.02809E-01, 1.92811E+02, 1.51406E-01, 1.45592E+02, 1.61384E-02 }, + { 7.91612E+01, 5.05468E-02, 1.19118E+02, 1.05494E-01, 1.98279E+02, 1.55232E-01, 1.48432E+02, 1.64533E-02 }, + { 8.11380E+01, 5.17150E-02, 1.22697E+02, 1.08207E-01, 2.03835E+02, 1.59099E-01, 1.51300E+02, 1.67712E-02 }, + { 8.31423E+01, 5.28967E-02, 1.26335E+02, 1.10947E-01, 2.09478E+02, 1.63007E-01, 1.54196E+02, 1.70922E-02 }, + { 8.51743E+01, 5.40921E-02, 1.30033E+02, 1.13714E-01, 2.15208E+02, 1.66954E-01, 1.57119E+02, 1.74162E-02 }, + { 8.72341E+01, 5.53009E-02, 1.33791E+02, 1.16507E-01, 2.21025E+02, 1.70941E-01, 1.60070E+02, 1.77433E-02 }, + { 8.93216E+01, 5.65232E-02, 1.37608E+02, 1.19325E-01, 2.26930E+02, 1.74967E-01, 1.63048E+02, 1.80734E-02 }, + { 9.14369E+01, 5.77590E-02, 1.41484E+02, 1.22168E-01, 2.32921E+02, 1.79031E-01, 1.66053E+02, 1.84065E-02 }, + { 9.35801E+01, 5.90083E-02, 1.45419E+02, 1.25036E-01, 2.38999E+02, 1.83132E-01, 1.69086E+02, 1.87427E-02 }, + { 9.57512E+01, 6.02710E-02, 1.49413E+02, 1.27926E-01, 2.45164E+02, 1.87269E-01, 1.72146E+02, 1.90819E-02 }, + { 9.79503E+01, 6.15471E-02, 1.53465E+02, 1.30840E-01, 2.51416E+02, 1.91443E-01, 1.75234E+02, 1.94242E-02 }, + { 1.00177E+02, 6.28365E-02, 1.57576E+02, 1.33776E-01, 2.57754E+02, 1.95652E-01, 1.78349E+02, 1.97695E-02 }, + { 1.02432E+02, 6.41393E-02, 1.61745E+02, 1.36733E-01, 2.64178E+02, 1.99896E-01, 1.81492E+02, 2.01178E-02 }, + { 1.04716E+02, 6.54554E-02, 1.65972E+02, 1.39711E-01, 2.70687E+02, 2.04173E-01, 1.84662E+02, 2.04692E-02 }, + { 1.07027E+02, 6.67848E-02, 1.70256E+02, 1.42710E-01, 2.77283E+02, 2.08484E-01, 1.87860E+02, 2.08237E-02 }, + { 1.09366E+02, 6.81274E-02, 1.74597E+02, 1.45728E-01, 2.83964E+02, 2.12827E-01, 1.91085E+02, 2.11812E-02 }, + { 1.11734E+02, 6.94833E-02, 1.78996E+02, 1.48765E-01, 2.90730E+02, 2.17202E-01, 1.94337E+02, 2.15417E-02 }, + { 1.14130E+02, 7.08524E-02, 1.83451E+02, 1.51821E-01, 2.97581E+02, 2.21608E-01, 1.97617E+02, 2.19053E-02 }, + { 1.16554E+02, 7.22346E-02, 1.87962E+02, 1.54894E-01, 3.04516E+02, 2.26044E-01, 2.00924E+02, 2.22719E-02 }, + { 1.19007E+02, 7.36300E-02, 1.92529E+02, 1.57984E-01, 3.11536E+02, 2.30510E-01, 2.04259E+02, 2.26415E-02 }, + { 1.21488E+02, 7.50386E-02, 1.97152E+02, 1.61091E-01, 3.18640E+02, 2.35005E-01, 2.07621E+02, 2.30142E-02 }, + { 1.23997E+02, 7.64602E-02, 2.01830E+02, 1.64213E-01, 3.25827E+02, 2.39528E-01, 2.11011E+02, 2.33899E-02 }, + { 1.26534E+02, 7.78949E-02, 2.06563E+02, 1.67351E-01, 3.33098E+02, 2.44079E-01, 2.14428E+02, 2.37687E-02 }, + { 1.29100E+02, 7.93426E-02, 2.11351E+02, 1.70503E-01, 3.40451E+02, 2.48656E-01, 2.17873E+02, 2.41505E-02 }, + { 1.31695E+02, 8.08033E-02, 2.16192E+02, 1.73669E-01, 3.47887E+02, 2.53260E-01, 2.21345E+02, 2.45354E-02 }, + { 1.34318E+02, 8.22770E-02, 2.21087E+02, 1.76848E-01, 3.55405E+02, 2.57889E-01, 2.24844E+02, 2.49233E-02 }, + { 1.36969E+02, 8.37636E-02, 2.26035E+02, 1.80039E-01, 3.63004E+02, 2.62543E-01, 2.28371E+02, 2.53143E-02 }, + { 1.39649E+02, 8.52632E-02, 2.31036E+02, 1.83243E-01, 3.70685E+02, 2.67221E-01, 2.31925E+02, 2.57082E-02 }, + { 1.42358E+02, 8.67757E-02, 2.36089E+02, 1.86458E-01, 3.78447E+02, 2.71922E-01, 2.35507E+02, 2.61053E-02 }, + { 1.45095E+02, 8.83010E-02, 2.41195E+02, 1.89684E-01, 3.86289E+02, 2.76646E-01, 2.39117E+02, 2.65054E-02 }, + { 1.47860E+02, 8.98391E-02, 2.46351E+02, 1.92920E-01, 3.94211E+02, 2.81392E-01, 2.42753E+02, 2.69085E-02 }, + { 1.50655E+02, 9.13901E-02, 2.51558E+02, 1.96165E-01, 4.02213E+02, 2.86159E-01, 2.46417E+02, 2.73146E-02 }, + { 1.53478E+02, 9.29538E-02, 2.56816E+02, 1.99419E-01, 4.10293E+02, 2.90947E-01, 2.50109E+02, 2.77238E-02 }, + { 1.56329E+02, 9.45302E-02, 2.62123E+02, 2.02682E-01, 4.18453E+02, 2.95755E-01, 2.53828E+02, 2.81361E-02 }, + { 1.59209E+02, 9.61194E-02, 2.67480E+02, 2.05953E-01, 4.26690E+02, 3.00583E-01, 2.57574E+02, 2.85514E-02 }, + { 1.62118E+02, 9.77213E-02, 2.72886E+02, 2.09231E-01, 4.35005E+02, 3.05429E-01, 2.61348E+02, 2.89697E-02 }, + { 1.65056E+02, 9.93358E-02, 2.78340E+02, 2.12515E-01, 4.43397E+02, 3.10294E-01, 2.65150E+02, 2.93911E-02 }, + { 1.68023E+02, 1.00963E-01, 2.83842E+02, 2.15806E-01, 4.51865E+02, 3.15175E-01, 2.68979E+02, 2.98155E-02 }, + { 1.71018E+02, 1.02603E-01, 2.89392E+02, 2.19103E-01, 4.60409E+02, 3.20074E-01, 2.72835E+02, 3.02429E-02 }, + { 1.74042E+02, 1.04255E-01, 2.94988E+02, 2.22404E-01, 4.69029E+02, 3.24989E-01, 2.76719E+02, 3.06734E-02 }, + { 1.77094E+02, 1.05920E-01, 3.00630E+02, 2.25710E-01, 4.77724E+02, 3.29920E-01, 2.80630E+02, 3.11070E-02 }, + { 1.80176E+02, 1.07597E-01, 3.06318E+02, 2.29020E-01, 4.86494E+02, 3.34866E-01, 2.84568E+02, 3.15436E-02 }, + { 1.83286E+02, 1.09287E-01, 3.12051E+02, 2.32334E-01, 4.95337E+02, 3.39826E-01, 2.88534E+02, 3.19832E-02 }, + { 1.86425E+02, 1.10989E-01, 3.17828E+02, 2.35651E-01, 5.04254E+02, 3.44800E-01, 2.92528E+02, 3.24259E-02 }, + { 1.89593E+02, 1.12704E-01, 3.23650E+02, 2.38971E-01, 5.13243E+02, 3.49788E-01, 2.96549E+02, 3.28716E-02 }, + { 1.92790E+02, 1.14431E-01, 3.29515E+02, 2.42292E-01, 5.22305E+02, 3.54788E-01, 3.00597E+02, 3.33203E-02 }, + { 1.96016E+02, 1.16171E-01, 3.35423E+02, 2.45616E-01, 5.31438E+02, 3.59801E-01, 3.04673E+02, 3.37721E-02 }, + { 1.99270E+02, 1.17922E-01, 3.41373E+02, 2.48940E-01, 5.40643E+02, 3.64825E-01, 3.08776E+02, 3.42270E-02 }, + { 2.02554E+02, 1.19687E-01, 3.47364E+02, 2.52266E-01, 5.49918E+02, 3.69861E-01, 3.12907E+02, 3.46848E-02 }, + { 2.05866E+02, 1.21463E-01, 3.53397E+02, 2.55592E-01, 5.59263E+02, 3.74907E-01, 3.17065E+02, 3.51458E-02 }, + { 2.09207E+02, 1.23252E-01, 3.59470E+02, 2.58918E-01, 5.68678E+02, 3.79963E-01, 3.21251E+02, 3.56097E-02 }, + { 2.12577E+02, 1.25053E-01, 3.65584E+02, 2.62244E-01, 5.78161E+02, 3.85030E-01, 3.25464E+02, 3.60767E-02 }, + { 2.15976E+02, 1.26866E-01, 3.71736E+02, 2.65568E-01, 5.87712E+02, 3.90105E-01, 3.29705E+02, 3.65468E-02 }, + { 2.19404E+02, 1.28691E-01, 3.77928E+02, 2.68892E-01, 5.97332E+02, 3.95189E-01, 3.33973E+02, 3.70199E-02 }, + { 2.22861E+02, 1.30528E-01, 3.84157E+02, 2.72215E-01, 6.07018E+02, 4.00282E-01, 3.38268E+02, 3.74960E-02 }, + { 2.26346E+02, 1.32378E-01, 3.90424E+02, 2.75535E-01, 6.16771E+02, 4.05382E-01, 3.42591E+02, 3.79752E-02 }, + { 2.29861E+02, 1.34240E-01, 3.96728E+02, 2.78853E-01, 6.26589E+02, 4.10490E-01, 3.46941E+02, 3.84574E-02 }, + { 2.33405E+02, 1.36113E-01, 4.03069E+02, 2.82169E-01, 6.36473E+02, 4.15605E-01, 3.51319E+02, 3.89427E-02 }, + { 2.36977E+02, 1.37999E-01, 4.09445E+02, 2.85482E-01, 6.46422E+02, 4.20727E-01, 3.55724E+02, 3.94310E-02 }, + { 2.40578E+02, 1.39897E-01, 4.15856E+02, 2.88792E-01, 6.56435E+02, 4.25855E-01, 3.60157E+02, 3.99223E-02 }, + { 2.44209E+02, 1.41807E-01, 4.22302E+02, 2.92099E-01, 6.66511E+02, 4.30988E-01, 3.64617E+02, 4.04167E-02 }, + { 2.47868E+02, 1.43728E-01, 4.28783E+02, 2.95402E-01, 6.76651E+02, 4.36127E-01, 3.69105E+02, 4.09141E-02 }, + { 2.51557E+02, 1.45662E-01, 4.35296E+02, 2.98701E-01, 6.86853E+02, 4.41272E-01, 3.73620E+02, 4.14146E-02 }, + { 2.55274E+02, 1.47608E-01, 4.41842E+02, 3.01996E-01, 6.97116E+02, 4.46421E-01, 3.78162E+02, 4.19181E-02 }, + { 2.59020E+02, 1.49565E-01, 4.48421E+02, 3.05287E-01, 7.07441E+02, 4.51574E-01, 3.82732E+02, 4.24247E-02 }, + { 2.62795E+02, 1.51534E-01, 4.55031E+02, 3.08573E-01, 7.17827E+02, 4.56732E-01, 3.87329E+02, 4.29343E-02 }, + { 2.66599E+02, 1.53516E-01, 4.61673E+02, 3.11854E-01, 7.28272E+02, 4.61893E-01, 3.91954E+02, 4.34469E-02 }, + { 2.70432E+02, 1.55509E-01, 4.68345E+02, 3.15131E-01, 7.38777E+02, 4.67058E-01, 3.96606E+02, 4.39626E-02 }, + { 2.74294E+02, 1.57513E-01, 4.75047E+02, 3.18402E-01, 7.49341E+02, 4.72227E-01, 4.01286E+02, 4.44814E-02 }, + { 2.78185E+02, 1.59530E-01, 4.81778E+02, 3.21668E-01, 7.59963E+02, 4.77398E-01, 4.05993E+02, 4.50031E-02 }, + { 2.82105E+02, 1.61558E-01, 4.88538E+02, 3.24928E-01, 7.70643E+02, 4.82572E-01, 4.10728E+02, 4.55279E-02 }, + { 2.86054E+02, 1.63598E-01, 4.95326E+02, 3.28183E-01, 7.81381E+02, 4.87749E-01, 4.15490E+02, 4.60558E-02 }, + { 2.90032E+02, 1.65650E-01, 5.02142E+02, 3.31432E-01, 7.92175E+02, 4.92927E-01, 4.20279E+02, 4.65867E-02 }, + { 2.94039E+02, 1.67713E-01, 5.08986E+02, 3.34676E-01, 8.03025E+02, 4.98108E-01, 4.25096E+02, 4.71206E-02 }, + { 2.98075E+02, 1.69788E-01, 5.15855E+02, 3.37913E-01, 8.13930E+02, 5.03291E-01, 4.29940E+02, 4.76576E-02 }, + { 3.02140E+02, 1.71874E-01, 5.22751E+02, 3.41144E-01, 8.24891E+02, 5.08475E-01, 4.34812E+02, 4.81976E-02 }, + { 3.06234E+02, 1.73973E-01, 5.29672E+02, 3.44369E-01, 8.35906E+02, 5.13661E-01, 4.39711E+02, 4.87407E-02 }, + { 3.10356E+02, 1.76082E-01, 5.36619E+02, 3.47588E-01, 8.46975E+02, 5.18848E-01, 4.44638E+02, 4.92868E-02 }, + { 3.14508E+02, 1.78204E-01, 5.43589E+02, 3.50801E-01, 8.58097E+02, 5.24036E-01, 4.49592E+02, 4.98360E-02 }, + { 3.18689E+02, 1.80336E-01, 5.50584E+02, 3.54008E-01, 8.69273E+02, 5.29225E-01, 4.54574E+02, 5.03882E-02 }, + { 3.22898E+02, 1.82481E-01, 5.57602E+02, 3.57208E-01, 8.80500E+02, 5.34415E-01, 4.59583E+02, 5.09434E-02 }, + { 3.27137E+02, 1.84637E-01, 5.64643E+02, 3.60401E-01, 8.91779E+02, 5.39605E-01, 4.64619E+02, 5.15017E-02 }, + { 3.31404E+02, 1.86804E-01, 5.71706E+02, 3.63589E-01, 9.03110E+02, 5.44797E-01, 4.69683E+02, 5.20630E-02 }, + { 3.35701E+02, 1.88983E-01, 5.78791E+02, 3.66770E-01, 9.14491E+02, 5.49988E-01, 4.74775E+02, 5.26274E-02 }, + { 3.40026E+02, 1.91173E-01, 5.85897E+02, 3.69944E-01, 9.25923E+02, 5.55180E-01, 4.79894E+02, 5.31948E-02 }, + { 3.44381E+02, 1.93374E-01, 5.93024E+02, 3.73113E-01, 9.37405E+02, 5.60373E-01, 4.85040E+02, 5.37652E-02 }, + { 3.48764E+02, 1.95587E-01, 6.00172E+02, 3.76275E-01, 9.48935E+02, 5.65566E-01, 4.90214E+02, 5.43387E-02 }, + { 3.53176E+02, 1.97811E-01, 6.07339E+02, 3.79430E-01, 9.60515E+02, 5.70759E-01, 4.95415E+02, 5.49152E-02 }, + { 3.57617E+02, 2.00047E-01, 6.14526E+02, 3.82580E-01, 9.72143E+02, 5.75952E-01, 5.00643E+02, 5.54948E-02 }, + { 3.62087E+02, 2.02294E-01, 6.21732E+02, 3.85723E-01, 9.83819E+02, 5.81145E-01, 5.05899E+02, 5.60774E-02 }, + { 3.66586E+02, 2.04552E-01, 6.28956E+02, 3.88861E-01, 9.95542E+02, 5.86339E-01, 5.11183E+02, 5.66631E-02 }, + { 3.71114E+02, 2.06821E-01, 6.36198E+02, 3.91992E-01, 1.00731E+03, 5.91532E-01, 5.16494E+02, 5.72518E-02 }, + { 3.75671E+02, 2.09102E-01, 6.43458E+02, 3.95117E-01, 1.01913E+03, 5.96726E-01, 5.21832E+02, 5.78436E-02 }, + { 3.80257E+02, 2.11394E-01, 6.50736E+02, 3.98237E-01, 1.03099E+03, 6.01920E-01, 5.27198E+02, 5.84383E-02 }, + { 3.84871E+02, 2.13697E-01, 6.58029E+02, 4.01351E-01, 1.04290E+03, 6.07113E-01, 5.32591E+02, 5.90362E-02 }, + { 3.89515E+02, 2.16011E-01, 6.65340E+02, 4.04459E-01, 1.05485E+03, 6.12307E-01, 5.38012E+02, 5.96370E-02 }, + { 3.94187E+02, 2.18337E-01, 6.72666E+02, 4.07562E-01, 1.06685E+03, 6.17502E-01, 5.43460E+02, 6.02410E-02 }, + { 3.98888E+02, 2.20673E-01, 6.80007E+02, 4.10659E-01, 1.07890E+03, 6.22696E-01, 5.48936E+02, 6.08479E-02 }, + { 4.03618E+02, 2.23021E-01, 6.87364E+02, 4.13751E-01, 1.09098E+03, 6.27891E-01, 5.54439E+02, 6.14579E-02 }, + { 4.08377E+02, 2.25380E-01, 6.94736E+02, 4.16839E-01, 1.10311E+03, 6.33086E-01, 5.59970E+02, 6.20710E-02 }, + { 4.13165E+02, 2.27750E-01, 7.02121E+02, 4.19921E-01, 1.11529E+03, 6.38282E-01, 5.65528E+02, 6.26870E-02 }, + { 4.17982E+02, 2.30131E-01, 7.09521E+02, 4.22999E-01, 1.12750E+03, 6.43478E-01, 5.71113E+02, 6.33062E-02 }, + { 4.22828E+02, 2.32523E-01, 7.16935E+02, 4.26072E-01, 1.13976E+03, 6.48675E-01, 5.76726E+02, 6.39283E-02 }, + { 4.27702E+02, 2.34926E-01, 7.24361E+02, 4.29141E-01, 1.15206E+03, 6.53873E-01, 5.82366E+02, 6.45536E-02 }, + { 4.32606E+02, 2.37341E-01, 7.31801E+02, 4.32205E-01, 1.16441E+03, 6.59071E-01, 5.88034E+02, 6.51818E-02 }, + { 4.37538E+02, 2.39766E-01, 7.39253E+02, 4.35266E-01, 1.17679E+03, 6.64270E-01, 5.93729E+02, 6.58131E-02 }, + { 4.42499E+02, 2.42202E-01, 7.46717E+02, 4.38323E-01, 1.18922E+03, 6.69471E-01, 5.99452E+02, 6.64474E-02 }, + { 4.47488E+02, 2.44649E-01, 7.54193E+02, 4.41376E-01, 1.20168E+03, 6.74672E-01, 6.05202E+02, 6.70848E-02 }, + { 4.52507E+02, 2.47107E-01, 7.61681E+02, 4.44426E-01, 1.21419E+03, 6.79876E-01, 6.10979E+02, 6.77253E-02 }, + { 4.57555E+02, 2.49577E-01, 7.69180E+02, 4.47473E-01, 1.22673E+03, 6.85080E-01, 6.16784E+02, 6.83687E-02 }, + { 4.62631E+02, 2.52057E-01, 7.76690E+02, 4.50516E-01, 1.23932E+03, 6.90286E-01, 6.22617E+02, 6.90152E-02 }, + { 4.67736E+02, 2.54547E-01, 7.84211E+02, 4.53557E-01, 1.25195E+03, 6.95495E-01, 6.28477E+02, 6.96648E-02 }, + { 4.72870E+02, 2.57049E-01, 7.91743E+02, 4.56596E-01, 1.26461E+03, 7.00705E-01, 6.34364E+02, 7.03174E-02 }, + { 4.78032E+02, 2.59562E-01, 7.99284E+02, 4.59632E-01, 1.27732E+03, 7.05917E-01, 6.40279E+02, 7.09730E-02 }, + { 4.83224E+02, 2.62086E-01, 8.06835E+02, 4.62667E-01, 1.29006E+03, 7.11132E-01, 6.46221E+02, 7.16317E-02 }, + { 4.88444E+02, 2.64620E-01, 8.14396E+02, 4.65699E-01, 1.30284E+03, 7.16349E-01, 6.52191E+02, 7.22934E-02 }, + { 4.93693E+02, 2.67165E-01, 8.21967E+02, 4.68730E-01, 1.31566E+03, 7.21569E-01, 6.58188E+02, 7.29582E-02 }, + { 4.98971E+02, 2.69721E-01, 8.29546E+02, 4.71760E-01, 1.32852E+03, 7.26792E-01, 6.64212E+02, 7.36260E-02 }, + { 5.04277E+02, 2.72288E-01, 8.37135E+02, 4.74789E-01, 1.34141E+03, 7.32019E-01, 6.70265E+02, 7.42968E-02 }, + { 5.09612E+02, 2.74866E-01, 8.44732E+02, 4.77817E-01, 1.35434E+03, 7.37248E-01, 6.76344E+02, 7.49707E-02 }, + { 5.14976E+02, 2.77454E-01, 8.52338E+02, 4.80845E-01, 1.36731E+03, 7.42481E-01, 6.82451E+02, 7.56477E-02 }, + { 5.20369E+02, 2.80053E-01, 8.59952E+02, 4.83872E-01, 1.38032E+03, 7.47719E-01, 6.88585E+02, 7.63276E-02 }, + { 5.25790E+02, 2.82663E-01, 8.67574E+02, 4.86899E-01, 1.39336E+03, 7.52960E-01, 6.94747E+02, 7.70106E-02 }, + { 5.31240E+02, 2.85284E-01, 8.75204E+02, 4.89927E-01, 1.40644E+03, 7.58205E-01, 7.00936E+02, 7.76967E-02 }, + { 5.36719E+02, 2.87915E-01, 8.82842E+02, 4.92955E-01, 1.41956E+03, 7.63456E-01, 7.07153E+02, 7.83858E-02 }, + { 5.42227E+02, 2.90557E-01, 8.90488E+02, 4.95985E-01, 1.43271E+03, 7.68710E-01, 7.13397E+02, 7.90780E-02 }, + { 5.47763E+02, 2.93210E-01, 8.98141E+02, 4.99015E-01, 1.44590E+03, 7.73970E-01, 7.19669E+02, 7.97731E-02 }, + { 5.53328E+02, 2.95874E-01, 9.05801E+02, 5.02047E-01, 1.45913E+03, 7.79236E-01, 7.25968E+02, 8.04714E-02 }, + { 5.58921E+02, 2.98548E-01, 9.13468E+02, 5.05080E-01, 1.47239E+03, 7.84507E-01, 7.32294E+02, 8.11726E-02 }, + { 5.64544E+02, 3.01232E-01, 9.21143E+02, 5.08115E-01, 1.48569E+03, 7.89783E-01, 7.38648E+02, 8.18770E-02 }, + { 5.70195E+02, 3.03928E-01, 9.28824E+02, 5.11153E-01, 1.49902E+03, 7.95066E-01, 7.45030E+02, 8.25843E-02 }, + { 5.75874E+02, 3.06634E-01, 9.36512E+02, 5.14193E-01, 1.51239E+03, 8.00355E-01, 7.51438E+02, 8.32947E-02 }, + { 5.81582E+02, 3.09350E-01, 9.44207E+02, 5.17236E-01, 1.52579E+03, 8.05651E-01, 7.57875E+02, 8.40082E-02 }, + { 5.87319E+02, 3.12078E-01, 9.51908E+02, 5.20282E-01, 1.53923E+03, 8.10954E-01, 7.64338E+02, 8.47246E-02 }, + { 5.93085E+02, 3.14816E-01, 9.59616E+02, 5.23332E-01, 1.55270E+03, 8.16264E-01, 7.70830E+02, 8.54442E-02 }, + { 5.98879E+02, 3.17564E-01, 9.67330E+02, 5.26385E-01, 1.56621E+03, 8.21581E-01, 7.77348E+02, 8.61667E-02 }, + { 6.04701E+02, 3.20323E-01, 9.75051E+02, 5.29442E-01, 1.57975E+03, 8.26906E-01, 7.83894E+02, 8.68923E-02 }, + { 6.10553E+02, 3.23092E-01, 9.82777E+02, 5.32503E-01, 1.59333E+03, 8.32240E-01, 7.90468E+02, 8.76210E-02 }, + { 6.16433E+02, 3.25872E-01, 9.90510E+02, 5.35568E-01, 1.60694E+03, 8.37581E-01, 7.97069E+02, 8.83527E-02 }, + { 6.22341E+02, 3.28663E-01, 9.98250E+02, 5.38639E-01, 1.62059E+03, 8.42932E-01, 8.03697E+02, 8.90874E-02 } }; - double wavelengths[] = { - 4.00e-1, - 4.50e-1, - 5.00e-1, - 6.00e-1, - 6.50e-1, - 7.00e-1, - 7.50e-1 + const double x[] = { + 1.396263402e-5, 8.377720036e-1, 1.675530045e+0, 2.513288086e+0, + 3.351046126e+0, 4.188804167e+0, 5.026562208e+0, 5.864320249e+0, + 6.702078290e+0, 7.539836331e+0, 8.377594372e+0, 9.215352413e+0, + 1.005311045e+1, 1.089086850e+1, 1.172862654e+1, 1.256638458e+1, + 1.340414262e+1, 1.424190066e+1, 1.507965870e+1, 1.591741674e+1, + 1.675517478e+1, 1.759293282e+1, 1.843069086e+1, 1.926844890e+1, + 2.010620695e+1, 2.094396499e+1, 2.178172303e+1, 2.261948107e+1, + 2.345723911e+1, 2.429499715e+1, 2.513275519e+1, 2.597051323e+1, + 2.680827127e+1, 2.764602931e+1, 2.848378736e+1, 2.932154540e+1, + 3.015930344e+1, 3.099706148e+1, 3.183481952e+1, 3.267257756e+1, + 3.351033560e+1, 3.434809364e+1, 3.518585168e+1, 3.602360972e+1, + 3.686136776e+1, 3.769912581e+1, 3.853688385e+1, 3.937464189e+1, + 4.021239993e+1, 4.105015797e+1, 4.188791601e+1, 4.272567405e+1, + 4.356343209e+1, 4.440119013e+1, 4.523894817e+1, 4.607670622e+1, + 4.691446426e+1, 4.775222230e+1, 4.858998034e+1, 4.942773838e+1, + 5.026549642e+1, 5.110325446e+1, 5.194101250e+1, 5.277877054e+1, + 5.361652858e+1, 5.445428662e+1, 5.529204467e+1, 5.612980271e+1, + 5.696756075e+1, 5.780531879e+1, 5.864307683e+1, 5.948083487e+1, + 6.031859291e+1, 6.115635095e+1, 6.199410899e+1, 6.283186703e+1, + 6.366962508e+1, 6.450738312e+1, 6.534514116e+1, 6.618289920e+1, + 6.702065724e+1, 6.785841528e+1, 6.869617332e+1, 6.953393136e+1, + 7.037168940e+1, 7.120944744e+1, 7.204720548e+1, 7.288496353e+1, + 7.372272157e+1, 7.456047961e+1, 7.539823765e+1, 7.623599569e+1, + 7.707375373e+1, 7.791151177e+1, 7.874926981e+1, 7.958702785e+1, + 8.042478589e+1, 8.126254394e+1, 8.210030198e+1, 8.293806002e+1, + 8.377581806e+1, 8.461357610e+1, 8.545133414e+1, 8.628909218e+1, + 8.712685022e+1, 8.796460826e+1, 8.880236630e+1, 8.964012435e+1, + 9.047788239e+1, 9.131564043e+1, 9.215339847e+1, 9.299115651e+1, + 9.382891455e+1, 9.466667259e+1, 9.550443063e+1, 9.634218867e+1, + 9.717994671e+1, 9.801770475e+1, 9.885546280e+1, 9.969322084e+1, + 1.005309789e+2, 1.013687369e+2, 1.022064950e+2, 1.030442530e+2, + 1.038820110e+2, 1.047197691e+2, 1.055575271e+2, 1.063952852e+2, + 1.072330432e+2, 1.080708012e+2, 1.089085593e+2, 1.097463173e+2, + 1.105840754e+2, 1.114218334e+2, 1.122595915e+2, 1.130973495e+2, + 1.139351075e+2, 1.147728656e+2, 1.156106236e+2, 1.164483817e+2, + 1.172861397e+2, 1.181238977e+2, 1.189616558e+2, 1.197994138e+2, + 1.206371719e+2, 1.214749299e+2, 1.223126879e+2, 1.231504460e+2, + 1.239882040e+2, 1.248259621e+2, 1.256637201e+2, 1.265014781e+2, + 1.273392362e+2, 1.281769942e+2, 1.290147523e+2, 1.298525103e+2, + 1.306902684e+2, 1.315280264e+2, 1.323657844e+2, 1.332035425e+2, + 1.340413005e+2, 1.348790586e+2, 1.357168166e+2, 1.365545746e+2, + 1.373923327e+2, 1.382300907e+2, 1.390678488e+2, 1.399056068e+2, + 1.407433648e+2, 1.415811229e+2, 1.424188809e+2, 1.432566390e+2, + 1.440943970e+2, 1.449321550e+2, 1.457699131e+2, 1.466076711e+2, + 1.474454292e+2, 1.482831872e+2, 1.491209453e+2, 1.499587033e+2, + 1.507964613e+2, 1.516342194e+2, 1.524719774e+2, 1.533097355e+2, + 1.541474935e+2, 1.549852515e+2, 1.558230096e+2, 1.566607676e+2, + 1.574985257e+2, 1.583362837e+2, 1.591740417e+2, 1.600117998e+2, + 1.608495578e+2, 1.616873159e+2, 1.625250739e+2, 1.633628319e+2, + 1.642005900e+2, 1.650383480e+2, 1.658761061e+2, 1.667138641e+2, + 1.675516222e+2, 1.683893802e+2, 1.692271382e+2, 1.700648963e+2, + 1.709026543e+2, 1.717404124e+2, 1.725781704e+2, 1.734159284e+2, + 1.742536865e+2, 1.750914445e+2, 1.759292026e+2, 1.767669606e+2, + 1.776047186e+2, 1.784424767e+2, 1.792802347e+2, 1.801179928e+2, + 1.809557508e+2, 1.817935089e+2, 1.826312669e+2, 1.834690249e+2, + 1.843067830e+2, 1.851445410e+2, 1.859822991e+2, 1.868200571e+2, + 1.876578151e+2, 1.884955732e+2, 1.893333312e+2, 1.901710893e+2, + 1.910088473e+2, 1.918466053e+2, 1.926843634e+2, 1.935221214e+2, + 1.943598795e+2, 1.951976375e+2, 1.960353955e+2, 1.968731536e+2, + 1.977109116e+2, 1.985486697e+2, 1.993864277e+2, 2.002241858e+2, + 2.010619438e+2, 2.018997018e+2, 2.027374599e+2 }; - const size_t nwavelengths = sizeof(wavelengths) / sizeof(double); + const size_t nx = sizeof(x)/sizeof(double); const size_t nscatt_angles = 1; - const size_t ngeoms = 1000; - const size_t ndirs = 100; + const size_t ngeoms = 100; + const size_t ndirs = 10; size_t i; (void)argc, (void)argv; @@ -842,51 +572,68 @@ main(int argc, char** argv) cylinder_init(&sampler_ctx.geometry, 64); - sampler_ctx.aspect_ratio = 0.263; - sampler_ctx.mean_radius = 0.983; + FOR_EACH(i, 0, nx) { + const double wavelength = 0.6; /* In micron */ + struct sschiff_estimator_state state; + double interval[2]; + struct sschiff_estimator_value* val; - distrib.sample = sample_cylinder; - distrib.context = &sampler_ctx; + sampler_ctx.aspect_ratio = 0.837; + sampler_ctx.mean_radius = (x[i] * 0.450) / (2*PI); + sampler_ctx.sigma = 1.18; - time_current(&t0); - CHECK(sschiff_integrate(dev, rng, &distrib, wavelengths, nwavelengths, - nscatt_angles, ngeoms, ndirs, &estimator), RES_OK); - time_current(&t1); - time_sub(&t0, &t1, &t0); - time_dump(&t0, TIME_MIN|TIME_SEC|TIME_MSEC, NULL, buf, sizeof(buf)); + distrib.sample = sample_cylinder; + distrib.context = &sampler_ctx; - FOR_EACH(i, 0, nwavelengths) { - struct sschiff_estimator_state state; - double interval0[2], interval1[2], interval2[2]; - struct sschiff_estimator_value* val; + time_current(&t0); + CHECK(sschiff_integrate(dev, rng, &distrib, &wavelength, 1, nscatt_angles, + ngeoms, ndirs, &estimator), RES_OK); + time_current(&t1); + time_sub(&t0, &t1, &t0); + time_dump(&t0, TIME_MIN|TIME_SEC|TIME_MSEC, NULL, buf, sizeof(buf)); CHECK(sschiff_estimator_get_wavelength_state - (estimator, wavelengths[i], &state), RES_OK); + (estimator, wavelength, &state), RES_OK); + CHECK(eq_eps(state.wavelength, wavelength, 1.e-6), 1); - CHECK(eq_eps(state.wavelength, wavelengths[i], 1.e-12), 1); - val = state.values + SSCHIFF_EXTINCTION_CROSS_SECTION; + printf("%u - x = %g - Wavelength = %g micron - %s\n", + (unsigned)i, x[i], wavelength, buf); - interval0[0] = results[i].extinction_E - 3*results[i].extinction_SE; - interval0[1] = results[i].extinction_E + 3*results[i].extinction_SE; - interval1[0] = val->E - 3*val->SE; - interval1[1] = val->E + 3*val->SE; - interval2[0] = MMAX(interval0[0], interval1[0]); - interval2[1] = MMIN(interval0[1], interval1[1]); - - printf( -"%2.1d - Wavelength: %7.2g - " -"Extinction cross section ~ %12.7g +/- %12.7g ~ %12.7g +/- %12.7g (%6.2g)\n", - (int)i, - wavelengths[i], - results[i].extinction_E, - results[i].extinction_SE, - val->E, val->SE, interval2[1] - interval2[0]); - - CHECK(interval2[0] <= interval2[1], 1); + val = state.values + SSCHIFF_EXTINCTION_CROSS_SECTION; + compute_intersection(interval, results[i].extinction_E, + results[i].extinction_SE, val->E, val->SE); + printf(" Extinction ~ %9.3g +/- %9.3g ~ %9.3g +/- %9.3g (%9.3g)\n", + results[i].extinction_E, results[i].extinction_SE, + val->E, val->SE, interval[1] - interval[0]); + CHECK(interval[0] <= interval[1], 1); + + val = state.values + SSCHIFF_ABSORPTION_CROSS_SECTION; + compute_intersection(interval, results[i].absorption_E, + results[i].absorption_SE, val->E, val->SE); + printf(" Absorption ~ %9.3g +/- %9.3g ~ %9.3g +/- %9.3g (%9.3g)\n", + results[i].absorption_E, results[i].absorption_SE, + val->E, val->SE, interval[1] - interval[0]); + CHECK(interval[0] <= interval[1], 1); + + val = state.values + SSCHIFF_SCATTERING_CROSS_SECTION; + compute_intersection(interval, results[i].scattering_E, + results[i].scattering_SE, val->E, val->SE); + printf(" Scattering ~ %9.3g +/- %9.3g ~ %9.3g +/- %9.3g (%9.3g)\n", + results[i].scattering_E, results[i].scattering_SE, + val->E, val->SE, interval[1] - interval[0]); + CHECK(interval[0] <= interval[1], 1); + + val = state.values + SSCHIFF_AVERAGE_PROJECTED_AREA; + compute_intersection(interval, results[i].avg_proj_area_E, + results[i].avg_proj_area_SE, val->E, val->SE); + printf(" Proj Area ~ %9.3g +/- %9.3g ~ %9.3g +/- %9.3g (%9.3g)\n", + results[i].avg_proj_area_E, results[i].avg_proj_area_SE, + val->E, val->SE, interval[1] - interval[0]); + CHECK(interval[0] <= interval[1], 1); + + CHECK(sschiff_estimator_ref_put(estimator), RES_OK); } - printf("Elapsed time: %s\n", buf); - CHECK(sschiff_estimator_ref_put(estimator), RES_OK); CHECK(sschiff_device_ref_put(dev), RES_OK); CHECK(ssp_rng_ref_put(rng), RES_OK);