diff --git a/.gitignore b/.gitignore index bd1d077..b29b1d4 100644 --- a/.gitignore +++ b/.gitignore @@ -69,6 +69,11 @@ relposdensity_TestData* *~ *.pylint* +# Not test data # +################# +!tests/data/ +!tests/data/** + # Notebooks with output # ######################### *.ipynb diff --git a/tests/data/Z1Z2_all_z_disks_7.0precisionfilter_1numberoflocalisations_PERPL-relpos_150.0filter.csv b/tests/data/Z1Z2_all_z_disks_7.0precisionfilter_1numberoflocalisations_PERPL-relpos_150.0filter.csv new file mode 100644 index 0000000..a2d3ff7 --- /dev/null +++ b/tests/data/Z1Z2_all_z_disks_7.0precisionfilter_1numberoflocalisations_PERPL-relpos_150.0filter.csv @@ -0,0 +1,270 @@ +xx_separation,yy_separation,zz_separation,xy_separation,xz_separation,yz_separation,xyz_separation +-7.151650434930343181e+01,6.660300718149301247e+00,6.016999999999995907e+00,7.182597023361478250e+01,7.176917641539392889e+01,8.975739226168508722e+00,7.207755745722879226e+01 +6.355380152069847099e+01,2.228663977067299129e+01,4.496600000000000819e+01,6.734819967898235404e+01,7.785259689780640713e+01,5.018600869034858647e+01,8.097975769289557491e+01 +1.657676939183147624e+00,6.906013399324729107e+01,6.066340000000000288e+01,6.908002605674067809e+01,6.068604446159512378e+01,9.192034707683207273e+01,9.193529300306803975e+01 +2.019242826622939901e+01,-2.146871772400481859e+01,-1.306738000000000284e+02,2.947269923165875483e+01,1.322247184369353477e+02,1.324256313828746840e+02,1.339562690076124341e+02 +-1.350703058700019028e+02,-1.562633905252369004e+01,-3.894900000000001228e+01,1.359712101880405157e+02,1.405738671617732507e+02,4.196673770004559145e+01,1.414397207328984223e+02 +1.367213010107516311e+02,1.433764101684573689e+01,4.870519999999999072e+01,1.374710227648012051e+02,1.451375577068615144e+02,5.077168952248818101e+01,1.458440214991361756e+02 +-2.805145105601695832e+00,-1.985235738270148431e+00,1.215600000000000591e+01,3.436568055485904605e+00,1.247546291980707522e+01,1.231704091640947496e+01,1.263243191155156353e+01 +-5.996352202346315607e+00,3.677466555305386464e+00,1.255299999999999727e+01,7.034202158025682827e+00,1.391165154590149555e+01,1.308057985203215701e+01,1.438950343132010623e+01 +-1.353587268048067926e+01,-1.014988427411560679e+01,-2.240000000000009095e+00,1.691862878604890241e+01,1.371996535061890121e+01,1.039412097187344308e+01,1.706627082874646462e+01 +-1.634101778608237510e+01,-1.213512001238575522e+01,9.915999999999996817e+00,2.035411506305458929e+01,1.911428571213322769e+01,1.567125373781578368e+01,2.264104803228122975e+01 +-1.436368583147850586e+01,-3.332993443339819351e+01,6.028999999999996362e+00,3.629325006113449348e+01,1.557768633865171637e+01,3.387083362326092129e+01,3.679060805423089420e+01 +-2.036003803382482147e+01,-2.965246787809280704e+01,1.858199999999999363e+01,3.596943146617859810e+01,2.756486663742077781e+01,3.499370765239552838e+01,4.048568542090054478e+01 +4.927954652530752355e+01,-1.587029597265257053e+01,3.909700000000000841e+01,5.177200015451924742e+01,6.290508019818391006e+01,4.219528058041080243e+01,6.487615439435002429e+01 +-5.331082733183939126e+01,-1.119534229912994761e+01,4.498900000000000432e+01,5.447366336129673670e+01,6.975711025985232538e+01,4.636103763069468187e+01,7.064977084888444381e+01 +6.562056431138989865e+01,-3.735175960266815309e+00,2.918100000000001160e+01,6.572678297314887175e+01,7.181635761820045616e+01,2.941908055079484896e+01,7.191342545727754043e+01 +6.364323235678602941e+01,1.745963846074562298e+01,3.306800000000001205e+01,6.599469675663192447e+01,7.172136117517470666e+01,3.739427227771585649e+01,7.381593746610420226e+01 +6.963958455913234502e+01,1.378217190544023651e+01,2.051500000000001478e+01,7.099028102493784331e+01,7.259846391328500204e+01,2.471464115521620997e+01,7.389509608221408143e+01 +6.281541920578820282e+01,-5.720411698536963740e+00,4.133700000000001751e+01,6.307535176278975086e+01,7.519657212266332635e+01,4.173093192106737348e+01,7.541384202518570135e+01 +7.998425014286840451e+01,2.959475847313137820e+01,2.315200000000001523e+01,8.528382027090441397e+01,8.326761300119599696e+01,3.757478986079605932e+01,8.837049905935762695e+01 +7.717910503726670868e+01,2.760952273486122976e+01,3.530800000000002115e+01,8.196889654009174819e+01,8.487207502090100775e+01,4.482120714178524423e+01,8.924995722127975739e+01 +8.317545723961302429e+01,2.393205617955584330e+01,2.275500000000002387e+01,8.654998555748066735e+01,8.623193556924658765e+01,3.302322422146902170e+01,8.949128463152214863e+01 +8.598060234521472012e+01,2.591729191782599173e+01,1.059900000000001796e+01,8.980183739768219198e+01,8.663141913097085478e+01,2.800080036987887411e+01,9.042515579748675236e+01 +-9.136597650400290149e+01,2.782280247332550971e+01,6.159999999999996589e-01,9.550837659598083462e+01,9.136805305209261974e+01,2.782962079277560008e+01,9.551036308170805000e+01 +-8.934521593947283691e+01,3.555562386921803864e+01,6.851800000000000068e+01,9.616012687179826912e+01,1.125934453477245256e+02,7.719403288291982790e+01,1.180740713450685604e+02 +1.063974556169705465e+02,6.835308203349541145e+00,6.450899999999998613e+01,1.066167904225220724e+02,1.244260006701381371e+02,6.487012038862557972e+01,1.246136071261882421e+02 +-8.504841211299572024e+01,-5.847510236894493119e+01,-8.875600000000000023e+01,1.032112881423371533e+02,1.229262377970706410e+02,1.062871823554401516e+02,1.361256681746706647e+02 +3.173758478115632897e+01,4.727976006981498358e+01,1.337450000000000045e+02,5.694427100244902817e+01,1.374590823217625086e+02,1.418559154115868068e+02,1.453629080095756763e+02 +-2.392274805290799122e+01,-2.704352280301463907e+01,1.857199999999997431e+01,3.610609366852202129e+01,3.028559159737369555e+01,3.280657418258071800e+01,4.060257607591032070e+01 +-2.308351217383460607e+01,2.931981353488572495e+01,2.863499999999999091e+01,3.731621631409110051e+01,3.678058943627154775e+01,4.098310250237855712e+01,4.703682839010340899e+01 +1.433813713006748003e+01,-1.915509915504026139e+01,1.091899999999998272e+01,2.392697222801171364e+01,1.802239543902583918e+01,2.204863679775744600e+01,2.630065704502538892e+01 +-4.642081764560862212e+01,-2.169072818310542061e+01,-1.509010000000000673e+01,5.123846211587748911e+01,4.881191892250138409e+01,2.642345183967007571e+01,5.341433438703708703e+01 +1.842371682521115872e+01,-7.212639363190191943e+01,4.352199999999999136e+01,7.444225950358861610e+01,4.726095455717720029e+01,8.424002102530637615e+01,8.623116886601752640e+01 +-1.216452657805348281e+00,6.385158019292248355e+00,-9.795490000000003761e+01,6.500000000001061373e+00,9.796245296581078321e+01,9.816278661968262043e+01,9.817032359124635832e+01 +-5.968007803614818840e+01,-6.088463094738381187e+01,5.687999999999999545e+01,8.525637806052907308e+01,8.244420000461364850e+01,8.332030176133021371e+01,1.024889476968119908e+02 +8.571554184072010685e+01,4.135197561123459309e+01,3.835200000000000387e+01,9.516900756023665053e+01,9.390436633643953712e+01,5.639912934569239411e+01,1.026061202073267111e+02 +8.007056944696887513e+01,6.232610936388027767e+01,-3.083600000000001273e+01,1.014684187321366977e+02,8.580300104053392829e+01,6.953706065428970362e+01,1.060504544827613245e+02 +8.816597044536683825e+01,-1.757047681272433692e+01,-6.320830000000002258e+01,8.989972191280561731e+01,1.084828444200432074e+02,6.560496051608053847e+01,1.098965385664616150e+02 +-6.573243231690139510e+01,-8.148120851892053906e+01,4.175499999999999545e+01,1.046897320657562460e+02,7.787318333865658815e+01,9.155690780440222909e+01,1.127094495816558890e+02 +-1.400362548214980052e+01,1.111225830934281475e+02,-1.718400000000002592e+01,1.120014732045974455e+02,2.216734947268879452e+01,1.124434005593740693e+02,1.133120463851931277e+02 +-1.725554349859157810e+01,-7.817043058963099611e+01,-8.066000000000002501e+01,8.005229540744039696e+01,8.248508581211387991e+01,1.123238702082879570e+02,1.136415663390825870e+02 +-1.239693795620987657e+02,3.521779520093332394e+00,-2.943500000000000227e+01,1.240193936447032428e+02,1.274159577682941205e+02,2.964493474420459762e+01,1.274646195028246183e+02 +1.640122038670961047e+01,-1.307608502948278328e+02,-3.263299999999998136e+01,1.317854316682985711e+02,3.652276987268925268e+01,1.347713347074457602e+02,1.357656535689339705e+02 +1.021167622274297173e+02,-8.940887468359323975e+01,5.719000000000022510e+00,1.357268580642751772e+02,1.022767817679722384e+02,8.959159465701277725e+01,1.358472927996722888e+02 +-7.608129842285779887e+01,6.987621934744402097e+01,8.951299999999997681e+01,1.033007744404660997e+02,1.174774069287705203e+02,1.135573123990353395e+02,1.366880652032212140e+02 +-1.345867880909754604e+02,-4.120251370381083689e+00,4.811820000000001585e+01,1.346498421833470331e+02,1.429299293356193346e+02,4.829428167594097943e+01,1.429893043945604347e+02 +1.383075166921662458e+02,-2.267687867513359379e+01,4.035399999999998499e+01,1.401542364682586879e+02,1.440743366792082156e+02,4.628915793624601349e+01,1.458480555784017554e+02 +-7.321879195266956231e+01,-2.866737701625061163e+01,-1.160175900000000304e+01,7.863084636451628739e+01,7.413226225404422109e+01,3.092602976274088178e+01,7.948214146520112422e+01 +-8.525306828610337107e+01,5.878447369676359813e+01,6.909200000000004138e+01,1.035552992366904164e+02,1.097350906328281610e+02,9.071559299153314271e+01,1.244885716200507346e+02 +-8.183884806976493564e+01,-9.031131128830838861e+01,3.305400000000003047e+01,1.218758794840053952e+02,8.826190553906067748e+01,9.617016097841232636e+01,1.262786479021683590e+02 +2.767760709706635680e+01,1.643776339349415139e+01,-2.444299999999998363e+01,3.219083720564089646e+01,3.692573876877180794e+01,2.945607432059634689e+01,4.041918169631930624e+01 +-6.945947218197034090e+01,4.200205257204743248e+00,8.920999999999992269e+00,6.958634923604313371e+01,7.003001154360829617e+01,9.860322773755950010e+00,7.015585678331183317e+01 +6.574859475281482446e+01,-3.568233180426886975e+00,-2.629799999999997340e+01,6.584534911441937766e+01,7.081286970579476758e+01,2.653897307790747817e+01,7.090271365751642918e+01 +-6.216487691273141536e+01,-3.506861620536710689e+00,-4.718700300000000425e+01,6.226371334895342358e+01,7.804540456488810207e+01,4.731713569678116471e+01,7.812415280899146808e+01 +6.821080798048569704e+01,-3.892268843038982595e+01,1.580200000000002092e+01,7.853464203776505315e+01,7.001726593741497595e+01,4.200808111124813848e+01,8.010863376690363680e+01 +2.743080318425927544e+01,9.090963115460363042e+01,3.548699999999999477e+01,9.495793805680305866e+01,4.485282747312111695e+01,9.759041041857585697e+01,1.013722702172524492e+02 +-5.028714114400645485e+01,-9.976614373404390790e+01,-2.081599999999997408e+01,1.117232294556469441e+02,5.442520023332230039e+01,1.019146176540044166e+02,1.136458791861769839e+02 +-1.079192256406240631e+02,1.331543229217959379e+01,3.946970999999999918e+01,1.087375740027316908e+02,1.149104750244990925e+02,4.165523670094456321e+01,1.156793758951164506e+02 +1.253460281918705732e+02,-1.164015534790632955e+01,-8.160000000000025011e+00,1.258853446593391396e+02,1.256113545165293317e+02,1.421544288875280593e+01,1.261495366618546399e+02 +-2.867946738574391929e+01,-1.256321143277887131e+02,-1.076399999999998158e+01,1.288640368760831620e+02,3.063291604679431757e+01,1.260923940865212245e+02,1.293128133480998940e+02 +9.689027536622961634e+01,8.670942589739888717e+01,2.656600000000000250e+01,1.300240362394596332e+02,1.004663018954306040e+02,9.068779904406385128e+01,1.327102194859171789e+02 +4.738127995079776156e+01,-1.278440624754443888e+02,-4.533000000000015461e+00,1.363418131022218631e+02,4.759762366732051930e+01,1.279244011095042310e+02,1.364171473422648546e+02 +7.505888704786411836e+01,-1.114062990819502375e+02,-2.897599999999999909e+01,1.343324234874104093e+02,8.045772244392718164e+01,1.151128665751007532e+02,1.374220090669648471e+02 +4.078000479622642160e+01,-1.298323195849934564e+02,-1.968499999999997385e+01,1.360861491849928484e+02,4.528253544337208325e+01,1.313161468891768777e+02,1.375025062498866930e+02 +7.796474824107281165e+01,1.162039071275380593e+02,-3.627000000000009550e+00,1.399351635579826620e+02,7.804906852290977781e+01,1.162604969914781208e+02,1.399821600383397424e+02 +-1.238992705003329320e+01,-1.871522662686857075e+01,1.414575999999999851e+02,2.244482122896497955e+01,1.419991651034087567e+02,1.426902670312696841e+02,1.432271712970698445e+02 +-1.127308171569202386e+02,-9.467229195034951772e+01,-1.496999999999985675e+00,1.472110050233989682e+02,1.127407563655086591e+02,9.468412682246277257e+01,1.472186163805352521e+02 +4.089250384919068892e+01,-1.828422528798000712e+00,5.324000000000012278e+00,4.093336047773125586e+01,4.123762659339253389e+01,5.629218857338574189e+00,4.127814162483446125e+01 +2.361894162649332429e+01,-5.992725253542448627e+01,-2.602500000000000568e+01,6.441374077012015675e+01,3.514477242145266445e+01,6.533434182299949100e+01,6.947251704811219497e+01 +6.882032895231895964e+01,-4.134407240336122413e+01,2.237000000000023192e+00,8.028430730846325503e+01,6.885667611862622550e+01,4.140454675388169647e+01,8.031546656155192920e+01 +-4.369166050781859667e+01,3.899459965263758932e+01,5.857500000000001705e+01,5.856227454599043369e+01,7.307524767614870598e+01,7.036767316793617510e+01,8.282856164995237691e+01 +8.102272922051633941e+01,2.815665728845897320e+01,2.781500000000002615e+01,8.577575414999766679e+01,8.566421000243400385e+01,3.957867575677249761e+01,9.017291292289965554e+01 +1.046416708470096637e+02,-3.177059524696551307e+01,1.790000000000020464e+00,1.093583558764506165e+02,1.046569795935938174e+02,3.182098085142109056e+01,1.093730044389401428e+02 +-1.247143897283349361e+02,1.083794236417861612e+01,3.075999999999999091e+01,1.251844239512259520e+02,1.284517676223687204e+02,3.261347259476146121e+01,1.289081750704751528e+02 +-9.884907505150476936e+01,2.567996030880567560e+01,9.237999999999999545e+01,1.021303089195360201e+02,1.352967258973328626e+02,9.588287001055941516e+01,1.377122521782279705e+02 +-2.025743376255195471e+01,-5.496032874333650398e+00,-7.001000000000033197e+00,2.098975940786211325e+01,2.143309645487978088e+01,8.900583034597040921e+00,2.212654516638185953e+01 +-9.132544436713942559e+01,-1.213759096692228923e+00,1.572000000000002728e+00,9.133350973219128832e+01,9.133897291329412838e+01,1.986050136528064503e+00,9.134703708386102505e+01 +-1.721486911250758567e+01,-9.353126900368215502e+01,3.865999999999985448e+00,9.510231332622734612e+01,1.764362985784863014e+01,9.361113308490156726e+01,9.518085918922942312e+01 +1.127304211767295783e+02,-2.690468623339393162e+01,-1.085900000000003729e+01,1.158965486975354935e+02,1.132522217869602770e+02,2.901344554370228224e+01,1.164041574901868330e+02 +7.411057525463183993e+01,-9.231750990698992609e+01,2.293999999999982720e+00,1.183845429099577728e+02,7.414607070217847706e+01,9.234600733885132229e+01,1.184067668505463189e+02 +1.329878549392815330e+02,-2.140865335906028122e+01,-3.858000000000004093e+00,1.347000371195189530e+02,1.330438037841348944e+02,2.175349633158778317e+01,1.347552750878413121e+02 +9.598493343102745712e+00,-2.889450733687226602e+00,-3.926099999999999568e+01,1.002397126891300339e+01,4.041728832390401038e+01,3.936718235462636528e+01,4.052044077993220839e+01 +-2.824457196938601555e+01,2.860566647128143813e+01,2.633500000000000796e+01,4.020000000000026574e+01,3.861719915703140060e+01,3.888208301089584751e+01,4.805800895792523164e+01 +5.457297344588732813e+01,8.492383015080122277e+00,-1.604400000000003956e+01,5.522979268474891512e+01,5.688250492660748137e+01,1.815297510808690262e+01,5.751295450592279934e+01 +-5.669823801369420835e+01,1.315901995370859368e+01,1.701499999999998636e+01,5.820524031390661435e+01,5.919628720500567454e+01,2.150976594810135367e+01,6.064124194803087420e+01 +1.290178862918000959e+02,-4.991494445215721498e+00,-1.758999999999986130e+00,1.291144066322596018e+02,1.290298766340720817e+02,5.292362213286176953e+00,1.291263880119028045e+02 +1.243107414870246430e+02,-6.697611179328850994e+01,2.872899999999998499e+01,1.412053823336807170e+02,1.275872873371554590e+02,7.287767142237105134e+01,1.440982839627209273e+02 +2.365081087111320812e+00,-6.143556905517471023e+00,6.599800000000001887e+00,6.583076788245985078e+00,7.010775177440201844e+00,9.016687390130114110e+00,9.321709072908417326e+00 +1.394790977223237860e+01,-3.547775780083611608e+00,4.299199999999999022e+00,1.439204294046134969e+01,1.459545503416510392e+01,5.574032079720018018e+00,1.502045340993684341e+01 +3.750824548224045429e+00,1.990570057064905996e+01,6.749699999999998568e+01,2.025600157977626026e+01,6.760113677884092453e+01,7.037103043304351502e+01,7.047092030759849024e+01 +-1.166569134282144660e+02,-5.632137293605046580e+01,1.784000000000000341e+01,1.295412386076326072e+02,1.180131393133744382e+02,5.907929120598586081e+01,1.307639021289881782e+02 +1.266309021252818638e+02,3.326010563629824901e+01,4.349999999999994316e+00,1.309260096390340209e+02,1.267055952713324984e+02,3.354336189080812147e+01,1.309982538051573897e+02 +7.205892179119473440e+01,4.196765766981116030e+01,-1.081559999999999775e+02,8.338928288455285553e+01,1.299623351040966668e+02,1.160129502525061014e+02,1.365704683890334081e+02 +1.174443140617986501e+02,4.434527589894423727e+01,5.427400000000000091e+01,1.255375258637825766e+02,1.293786457706463864e+02,7.008687873313724026e+01,1.367674576644600677e+02 +-6.632768207694971352e+01,-8.815389401443098905e+00,1.279376000000000033e+02,6.691092960047701865e+01,1.441089549731761679e+02,1.282409473766436179e+02,1.443783293772302727e+02 +-1.129060888799904205e+02,-3.641567236540140584e+01,8.533699999999998909e+01,1.186334105553766705e+02,1.415280483691352913e+02,9.278202823189552362e+01,1.461379131813526158e+02 +-6.396260098983839271e+01,-1.495894630696056993e+01,1.345373999999999910e+02,6.568854085759406303e+01,1.489682057492311742e+02,1.353664732249995666e+02,1.497173884315372447e+02 +1.982782070284883957e+01,-5.608197148973931689e+01,2.319700000000000273e+01,5.948386335805917469e+01,3.051627898064115385e+01,6.069010080050892242e+01,6.384693265146145080e+01 +-1.368842256047610135e+02,-1.674630647590265653e+01,-4.120400000000000773e+01,1.379047859938152101e+02,1.429512533537748027e+02,4.447705472021343098e+01,1.439288005091405296e+02 +3.375707901839632541e+00,3.131229177413842990e+00,3.236400000000000432e+01,4.604345773291826482e+00,3.253957436474151166e+01,3.251512097719287198e+01,3.268988369511323810e+01 +2.205239535231521586e+01,-2.833150824126641965e+01,-1.049999999999897682e-01,3.590240242657531411e+01,2.205264532374322073e+01,2.833170281195514306e+01,3.590255596750404976e+01 +2.542810325415484840e+01,-2.520027906385257666e+01,3.225900000000001455e+01,3.580003491618418110e+01,4.107592380098056140e+01,4.093527996601522290e+01,4.819009837093101112e+01 +1.071760102077314514e+01,4.356516301311967254e+01,-3.849899999999999523e+01,4.486413266742340511e+01,3.996298252934178663e+01,5.813859672678460555e+01,5.911821547543673461e+01 +3.152577849038061686e+01,-5.317330994564690627e+01,-3.262770000000000437e+01,6.181646706178198514e+01,4.537005087846547013e+01,6.238563695167207612e+01,6.989880118636068573e+01 +5.152175447471381631e+01,4.617230322070099646e+00,5.069499999999999318e+01,5.172823213681032684e+01,7.228052441116273030e+01,5.090483121322613869e+01,7.242784702722934753e+01 +-5.755035457751364447e+01,-1.953763516913932108e+01,4.981400000000002137e+01,6.077633174188727594e+01,7.611489938243069275e+01,5.350844591279391693e+01,7.858242230931760730e+01 +-6.265581091228159494e+01,4.642768310958672373e+01,1.128600000000000136e+01,7.798256471801751388e+01,6.366415347018748605e+01,4.777974000477821903e+01,7.879501377625230418e+01 +7.341893118933512596e+00,4.043393383570582955e+01,-7.086299999999999955e+01,4.109508973101385010e+01,7.124232003219604792e+01,8.158717898316031381e+01,8.191685521917989377e+01 +2.456729653622096521e+01,-6.149470579571789131e+01,5.594619999999999749e+01,6.622047191013149359e+01,6.110261331185937195e+01,8.313588957449063344e+01,8.668983904957093500e+01 +-7.933884616664363421e+01,2.399848930548978387e+01,-2.851500000000001478e+01,8.288896187068607446e+01,8.430751885836959048e+01,3.726972919066770373e+01,8.765663252144729256e+01 +7.337341193305474008e+01,-2.862520096467051189e+00,-4.978499999999999659e+01,7.342922851290498443e+01,8.866906903592560241e+01,4.986722617413843750e+01,8.871526263840070214e+01 +-6.958481954159651650e+00,-8.321395850070985034e+00,8.857390000000000896e+01,1.084740061028190006e+01,8.884681329297279717e+01,8.896393308584990223e+01,8.923565353158986113e+01 +-5.743397783390901168e+01,-5.726175154650991317e+01,-4.663399999999998613e+01,8.110221945175226210e+01,7.398237469712600500e+01,7.384875182543188998e+01,9.355372764353209902e+01 +-4.794530867889989167e+01,-2.233775896738552547e+01,8.204400000000001114e+01,5.289355442773727844e+01,9.502615724270341957e+01,8.503054399264428298e+01,9.761632054118823021e+01 +-9.446232929057441652e+01,2.249507379402348306e+01,1.487299999999999045e+01,9.710386192113821835e+01,9.562603089118000810e+01,2.696728525451833036e+01,9.823627705180747682e+01 +-7.362062769183103228e+01,1.616532789706980111e+00,6.623799999999999955e+01,7.363837314878293228e+01,9.903266867927571582e+01,6.625772273675120516e+01,9.904586131686369299e+01 +1.471050223338170326e+01,-6.876544207697224920e+01,7.075800000000000978e+01,7.032129762738848910e+01,7.227097231917063880e+01,9.866803224976884223e+01,9.975860596459713747e+01 +2.901928726723781438e+01,9.608474679417849984e+01,1.610900000000000887e+01,1.003712986864266270e+02,3.319064498467107427e+01,9.742575864473108993e+01,1.016557793782521060e+02 +-6.999770403121510753e+01,5.993749273880894179e+00,8.214900000000000091e+01,7.025385113998900977e+01,1.079265341361501669e+02,8.236736751504291476e+01,1.080928388053516471e+02 +9.763095582071400713e+01,-4.660329243233718444e+01,3.687000000000000455e+01,1.081835033634987440e+02,1.043609143044761254e+02,5.942443744398374150e+01,1.142937762960002317e+02 +-3.610425006214063615e+01,-2.899263919429540692e+01,1.141864999999999952e+02,4.630431945294279217e+01,1.197583970116483414e+02,1.178097190799663707e+02,1.232178833702729008e+02 +-6.961429286644124659e+01,3.810628725951573870e+01,9.985990000000001032e+01,7.936144466930642238e+01,1.217298213228979336e+02,1.068835289308634913e+02,1.275548451765333340e+02 +-1.146266288783117488e+02,-7.536922406688972842e-01,-6.448300000000000409e+01,1.146291066876116673e+02,1.315192812366544786e+02,6.448740451432081500e+01,1.315214407957875835e+02 +-8.310085038793113199e+01,-5.392700218631580356e+01,-9.711070000000000846e+01,9.906499331247250950e+01,1.278132989547148668e+02,1.110792942869773157e+02,1.387240460572363361e+02 +-7.947108716928050853e+01,3.083555097826138081e+01,1.146717000000000155e+02,8.524367953109683071e+01,1.395179288720893567e+02,1.187452314201414936e+02,1.428848616225327817e+02 +-1.227386545295048563e+02,-6.287950925609584374e+01,3.878200000000001069e+01,1.379079765640863400e+02,1.287199317887993004e+02,7.387744045571315610e+01,1.432572983271728049e+02 +-1.310546055257400440e+02,-5.564521875680293306e+01,1.960400000000001342e+01,1.423787203201369778e+02,1.325127406912532990e+02,5.899751847741100619e+01,1.437220122876095445e+02 +-1.234816165578122309e+02,-7.358967300015501678e+01,-1.236000000000018417e+00,1.437468246605827176e+02,1.234878023277220933e+02,7.360005209420536687e+01,1.437521384049653648e+02 +-9.418158940266221180e+01,9.960099305523363000e+01,4.391370000000000573e+01,1.370785526623342605e+02,1.039162394917255199e+02,1.088520595362287935e+02,1.439407608972884134e+02 +1.048991904234353569e+02,-5.603583004211395746e+01,-8.241270000000000095e+01,1.189279378447336626e+02,1.334004995222362311e+02,9.965875460690185150e+01,1.446917672892650444e+02 +8.656964184475145885e+01,1.156223819633178209e+02,-3.370500000000001251e+01,1.444397386455677008e+02,9.289956896632267558e+01,1.204348879514211319e+02,1.483201440297301872e+02 +6.352040690951980650e+00,1.694573415490594925e+00,-2.263000000000005230e+00,6.574191965557226069e+00,6.743114261193394299e+00,2.827144860188010256e+00,6.952781385891347199e+00 +-6.723306220832455438e+00,2.228711165862478083e+00,-3.527999999999991587e+00,7.083078426780581971e+00,7.592735379234833104e+00,4.173001013759765065e+00,7.913076772022651362e+00 +-2.950766233738249866e+00,9.768468592050339794e+00,-4.591999999999984539e+00,1.020441081101911429e+01,5.458340898676966368e+00,1.079395398516566829e+01,1.119001626451202469e+01 +-1.531642873297823826e+01,-1.156291860921555781e+00,1.143600000000000705e+01,1.536001302082894782e+01,1.911478708048827713e+01,1.149430758539345909e+01,1.914972835316561373e+01 +4.283528320106779574e+01,4.003871268015609530e+01,-1.930799999999999272e+01,5.863411975974205603e+01,4.698574625262092042e+01,4.445106722098010721e+01,6.173134425881048060e+01 +4.955858942190025118e+01,3.781000151429361722e+01,-1.578000000000000114e+01,6.233498215287599464e+01,5.201021231920211818e+01,4.097077756780905133e+01,6.430130947344206049e+01 +-7.305041833349241642e+01,-3.203211484279199794e+01,1.673199999999999932e+01,7.976477919483046719e+01,7.494213396146552952e+01,3.613884620877948350e+01,8.150079646236630992e+01 +-7.117615445162300603e+01,-6.702935951864128583e+01,5.278999999999996362e+00,9.777003630970257575e+01,7.137165266071183112e+01,6.723691603932520877e+01,9.791244987743161232e+01 +-5.970734023952900316e+01,5.898443456981840427e+01,5.312799999999998590e+01,8.392931549822188231e+01,7.992215501648387033e+01,7.938354933814176206e+01,9.933133636471455929e+01 +-6.605938093048098381e+01,5.728986115432780934e+01,5.539099999999999113e+01,8.744123741119261695e+01,8.620907544985267634e+01,7.968871358029414864e+01,1.035090956438155274e+02 +1.358041727549061761e+01,4.861792124950989091e+01,-9.373699999999999477e+01,5.047900553695555459e+01,9.471564233206910899e+01,1.055951108509458720e+02,1.064648071852853661e+02 +-5.759573717613238841e+01,-1.841143826913139492e+01,-8.845799999999999841e+01,6.046693311224039746e+01,1.055560832205426181e+02,9.035374271793081391e+01,1.071497445820574654e+02 +1.545468115736002801e+01,1.362067657366060303e+01,-1.051899999999999977e+02,2.060024271701491116e+01,1.063192516418153559e+02,1.060681805742149351e+02,1.071881807850097061e+02 +-8.649258318460124428e+01,-6.818565137956284161e+01,1.671500000000000341e+01,1.101374141697563545e+02,8.809289512182674287e+01,7.020451751173339972e+01,1.113985692233093516e+02 +2.889684600846885587e+01,4.977421311043144669e+01,-1.051730000000000018e+02,5.755432216610502394e+01,1.090705626566451798e+02,1.163564876565232424e+02,1.198910335638149718e+02 +-8.633367249649381847e+01,-6.897606101589735772e+01,-5.252100000000001501e+01,1.105042985589277862e+02,1.010542351795898526e+02,8.669574634472441232e+01,1.223505432803657413e+02 +9.192585346796840895e+01,2.089802536569004587e+01,7.984799999999998477e+01,9.427136362649808632e+01,1.217623325984452265e+02,8.253744948921683999e+01,1.235426772576967522e+02 +-9.268571318744579912e+01,-7.067063443138795265e+01,-5.025800000000000978e+01,1.165546223879617003e+02,1.054348518899963381e+02,8.671911631776973195e+01,1.269285096579974095e+02 +9.487661970170665882e+01,1.112955677363970608e+01,8.443999999999996930e+01,9.552716891021077572e+01,1.270105765911731766e+02,8.517030370955397700e+01,1.274972689903589043e+02 +-7.087756956621888094e+01,-6.159237073198164580e+01,-1.050807000000000073e+02,9.390021299230404850e+01,1.267500823672481829e+02,1.218013696346468464e+02,1.409226862946135554e+02 +-8.151548386075592134e+01,4.990617087041209743e+01,1.079506999999999834e+02,9.557928645893905184e+01,1.352705723339084329e+02,1.189284638824397007e+02,1.441830559756933781e+02 +1.606464688440973987e+01,-3.163822246080235345e+01,1.437754999999999939e+02,3.548309456628599179e+01,1.446702017686115482e+02,1.472153915890903590e+02,1.480893122418022188e+02 +6.890679743988584960e+00,3.956138605481100967e+00,3.070599999999998886e+01,7.945596264597927672e+00,3.146966640011002525e+01,3.095980407989975092e+01,3.171736016758002208e+01 +4.501738904155354248e+01,-6.326221927924052579e+00,4.702299999999999613e+01,4.545972283241432876e+01,6.509783287574622079e+01,4.744663963740052992e+01,6.540450236031104225e+01 +-5.453677876037181704e+01,1.687554035999710322e+01,5.184800000000001319e+01,5.708803815160960227e+01,7.524942087190932227e+01,5.452521404306353503e+01,7.711847511458998383e+01 +3.864463427785813110e+01,-1.967127706913197471e+01,9.089650000000000318e+01,4.336319752970213415e+01,9.877034712260257265e+01,9.300071426489452620e+01,1.007101812740398401e+02 +4.259475376662521739e+01,-7.597352796573795786e+01,-5.590500000000002956e+01,8.709931113390163659e+01,7.028287183545815253e+01,9.432574397565475977e+01,1.034971450089334866e+02 +-1.187702890961991216e+02,2.532090890953986673e+01,1.429599999999999227e+01,1.214394087600887957e+02,1.196275770380505321e+02,2.907789614131694833e+01,1.222779850013890695e+02 +-1.118796093522105366e+02,2.927704751502096769e+01,4.500199999999998113e+01,1.156468438825722131e+02,1.205911563623271405e+02,5.368729379654737244e+01,1.240942081807206137e+02 +-1.241522331289015710e+02,-5.549153909136293805e+00,-3.402400000000000091e+01,1.242761847660299139e+02,1.287299870538839457e+02,3.447355051495686240e+01,1.288495350243857445e+02 +1.146764396608741663e+02,6.231110805230491678e+01,-1.566000000000002501e+00,1.305119151648619891e+02,1.146871316639060154e+02,6.233078326722694129e+01,1.305213099689094065e+02 +-5.050056208414025605e+01,9.110694446191155293e+01,-8.771600000000000819e+01,1.041670874124833119e+02,1.012146403778332200e+02,1.264696484741925673e+02,1.361795827427883694e+02 +-1.099175587182653544e+02,-2.653480517013304052e+01,-2.214740000000000109e+01,1.130750441078854038e+02,1.121266116555000139e+02,3.456303245053765494e+01,1.152235779984297750e+02 +7.600589674652292160e+01,-7.717326000472894520e+01,-6.102600000000001046e+01,1.083171657679455677e+02,9.747342723144141985e+01,9.838640523851604769e+01,1.243253034422220935e+02 +2.041981696859147632e+01,-1.108443390388492844e+01,2.514999999999999147e+01,2.323431944345196598e+01,3.239585505941734311e+01,2.748430779498718479e+01,3.423968603828559765e+01 +-2.021868852919942583e+01,3.340797411037419806e+01,-1.398300000000000409e+01,3.904981562056616440e+01,2.458291388018882984e+01,3.621625357708100523e+01,4.147785419956309028e+01 +3.258849689599446720e+01,3.010325514058786212e+01,4.361000000000018417e+00,4.436458159387475320e+01,3.287899710970882694e+01,3.041749975029718911e+01,4.457840756464469223e+01 +-5.018309526492521400e+01,-3.386067556372290710e+01,-2.478100000000000591e+01,6.053832174747051198e+01,5.596821429497781963e+01,4.196002038407156931e+01,6.541396151434540229e+01 +-1.823588114970334573e+01,-2.449393065014010062e+01,-6.147309999999998809e+01,3.053686296920458787e+01,6.412089663219084912e+01,6.617321710710362481e+01,6.863994481065654440e+01 +-2.584713758285579388e+01,-6.878000493438230478e+01,-3.521500000000000341e+01,7.347627916546080939e+01,4.368261377512884991e+01,7.727085675578894097e+01,8.147919872581425693e+01 +6.795323959824236226e+01,-3.326294226467507542e+01,-4.157000000000000739e+01,7.565755811549327348e+01,7.965995023784620344e+01,5.323991198436664973e+01,8.632572617707447193e+01 +2.654477488956763409e+01,-3.594046919645916205e+01,-7.661000000000001364e+01,4.468044762532544212e+01,8.107846307089091908e+01,8.462156596318477852e+01,8.868728488345695382e+01 +-1.226787035677261883e+01,-2.015775178211697494e+01,-8.927379999999999427e+01,2.359736425959518868e+01,9.011277384217278552e+01,9.152128890782428527e+01,9.233984495568540751e+01 +-4.576559484432073077e+01,5.558589055278571323e+01,-6.300660000000000593e+01,7.200195066802557164e+01,7.787375241385582569e+01,8.402156194755173146e+01,9.567712654318580690e+01 +6.941749190347400145e+01,-2.395896111752517754e+01,-6.168009999999999593e+01,7.343582232126199472e+01,9.286131012525549977e+01,6.616998227173014868e+01,9.590231872071687746e+01 +-5.628449053656368051e+00,-1.021879790447565028e+02,-2.123199999999999932e+01,1.023428673626120258e+02,2.196536507207665778e+01,1.043704023430569805e+02,1.045220566387841359e+02 +3.782267651017536991e+01,-9.715165228449404822e+01,2.244499999999999318e+01,1.042544886323874778e+02,4.398105141300479914e+01,9.971068932971648735e+01,1.066432202486431322e+02 +8.277159216091968119e+01,6.396393070431076922e+01,2.914200000000002433e+01,1.046065050558496097e+02,8.775188107871888121e+01,7.028969053243777410e+01,1.085899491849935998e+02 +-2.784035911797400331e+01,3.225115198224739288e+01,-1.059170000000000016e+02,4.260542688437445946e+01,1.095148231282768307e+02,1.107183259139245735e+02,1.141649389655150628e+02 +-3.159481539329863153e+01,-7.319269116696887068e+01,8.412509999999998911e+01,7.972077771823292380e+01,8.986247720682207785e+01,1.115087552180243904e+02,1.158983815676901514e+02 +-3.268768732536409516e+01,-9.073317878232046496e+00,-1.144237999999999857e+02,3.392359061184729541e+01,1.190012223009527048e+02,1.147829739280152950e+02,1.193466212610990027e+02 +-6.841897641462855972e+01,-5.835460621386300772e+01,-8.625409999999999400e+01,8.992450389076671513e+01,1.100950775486147677e+02,1.041394729830385586e+02,1.246041177763019050e+02 +-5.008562527978210710e+01,-5.189454923528319341e+01,-1.031051999999999964e+02,7.212221641075639411e+01,1.146265768777326315e+02,1.154284475654646798e+02,1.258264533674854135e+02 +-1.195031171832561085e+02,-2.793558811775801587e+01,-4.142510000000000048e+01,1.227248634140615593e+02,1.264793814284567190e+02,4.996434722374568338e+01,1.295277229399562202e+02 +-1.000732624686716008e+02,-1.286987718962336658e+01,-8.284199999999998454e+01,1.008974310872166029e+02,1.299132588503715056e+02,8.383573643068918102e+01,1.305491806331990006e+02 +-2.827736806309985695e+00,-1.016718289623713645e+02,-9.717899999999997362e+01,1.017111444238017413e+02,9.722013236179920170e+01,1.406446545217901019e+02,1.406730782381598885e+02 +-1.137188344425630930e+02,8.884883281746078865e+01,-2.143659999999999854e+01,1.443124679298367994e+02,1.157216536631544841e+02,9.139826536967261461e+01,1.458959088513456948e+02 +-1.123062768384625087e+01,-1.477633860857167747e+00,1.799000000000006594e+00,1.132741806414530927e+01,1.137380319739944490e+01,2.328046998398375145e+00,1.146938538021656484e+01 +4.757541141575711663e+01,3.192635789783889777e+01,6.661780000000000257e+01,5.729495702066153484e+01,8.186178014323016328e+01,7.387302352997946286e+01,8.786719169769483528e+01 +-6.231246802394889528e+01,-5.870295332233217778e+01,2.558199999999999363e+01,8.560888038048133808e+01,6.735935269311643481e+01,6.403495492903786612e+01,8.934942150903700053e+01 +-2.751885023101203842e+01,-1.640111221726328949e+01,8.856600000000003092e+01,3.203566137915891687e+01,9.274278125027778685e+01,9.007182044326221160e+01,9.418184515075152774e+01 +-7.354309570779514615e+01,-6.018058718318934552e+01,2.738100000000000023e+01,9.502783802654552403e+01,7.847487551621803448e+01,6.611673188167617354e+01,9.889392883791892075e+01 +9.385453126516222255e+01,-7.967474230266452651e+01,6.562999999999988177e+00,1.231127028376816668e+02,9.408371808130944203e+01,7.994459037355811404e+01,1.232875118128325056e+02 +-1.112226880569287459e+02,6.369114978857396636e+01,5.135799999999999699e+00,1.281680494507092192e+02,1.113411998329858221e+02,6.389787948774642246e+01,1.282709060607254230e+02 +1.050851589490084734e+02,-7.819710844180735876e+01,4.763999999999981583e+00,1.309873215238717421e+02,1.051930906825081138e+02,7.834209254710891912e+01,1.310739260722671418e+02 +-5.603449523317976855e+00,4.018488949786842568e+00,1.324259000000000128e+02,6.895426020191314009e+00,1.325443987400839774e+02,1.324868568736143004e+02,1.326053011414322498e+02 +-4.807791488412476610e+00,-2.817241808940161718e+01,1.597110000000000696e+01,2.857971308463566373e+01,1.667905555497827663e+01,3.238458238443169535e+01,3.273951794407020799e+01 +-2.606647449750016676e+01,-3.846165502511576051e+01,3.047399999999996112e+01,4.646245796339060519e+01,4.010144347437912415e+01,4.907100552537119142e+01,5.556459912570085891e+01 +4.781018683944421355e+01,1.713230966266928590e+01,-3.355900000000002592e+01,5.078710466250420552e+01,5.841250248553443924e+01,3.767920534429558899e+01,6.087311788466365670e+01 +7.084589041357685346e+01,-5.473198161502477888e+01,8.693999999999988404e+00,8.952502443451189151e+01,7.137734811894135589e+01,5.541818697420017514e+01,8.994618188672575343e+01 +-6.930855019379669102e+01,2.189166211218071112e+01,-7.319529999999998893e+01,7.268369830986810598e+01,1.008029120217071579e+02,7.639893200905288495e+01,1.031526632816134423e+02 +-9.356496466668249923e+01,3.074357472582778428e+01,4.324299999999999500e+01,9.848639499951300991e+01,1.030745344984761118e+02,5.305774623862718187e+01,1.075617359891523108e+02 +-7.387666133694438031e+01,-5.559396468778504641e+01,6.403299999999998704e+01,9.245782822454928862e+01,9.776495373749006035e+01,8.479925706459160040e+01,1.124663286899695294e+02 +-1.106042329591946327e+01,1.875625326431145368e+01,2.262070000000000647e+01,2.177452640127706829e+01,2.517993312093815561e+01,2.938525318939747066e+01,3.139786726021231544e+01 +-3.271265133234555833e+01,-5.008237495194407529e+00,-1.069499999999999318e+01,3.309380606699508576e+01,3.441657423671928484e+01,1.180954985629304232e+01,3.477906015118696104e+01 +-4.141827944840406417e+01,-2.223989103811618406e+00,3.387699999999998113e+01,4.147794594721348460e+01,5.350817695891058179e+01,3.394992277655240542e+01,5.355437544216123058e+01 +4.268213477456447436e+01,-4.562998324660293292e+01,-3.536000000000001364e+00,6.248087707450456918e+01,4.282835421673457432e+01,4.576678563199806149e+01,6.258085406895107639e+01 +5.931320495613908861e+01,-4.403389283076103311e+01,3.456100000000000705e+01,7.387178080972402938e+01,6.864778949950947151e+01,5.597719570352688123e+01,8.155674540465621192e+01 +4.060557062391308136e+01,-3.974968722274570609e+01,6.403000000000000114e+01,5.682297070727388899e+01,7.581987381744706056e+01,7.536496888015089723e+01,8.560777359562452205e+01 +1.805370266365207499e+01,5.465705645323578210e+01,-6.431999999999999318e+01,5.756153229370949731e+01,6.680567775172673350e+01,8.440649394526586491e+01,8.631565559039542279e+01 +4.433627796834844048e+01,1.714830766907834914e+01,-7.648990000000003420e+01,4.753703819128859465e+01,8.841046513901318349e+01,7.838857861910354075e+01,9.005817454295882385e+01 +-5.626549560700368602e+00,-9.682061732937472698e+01,-3.130500000000000682e+01,9.698396774725313207e+01,3.180662014359617373e+01,1.017557613358635251e+02,1.019112016659613147e+02 +5.860389333148850710e+01,-8.875152779751078924e+01,-2.080000000000040927e+00,1.063543605123965392e+02,5.864079393739893220e+01,8.877589811650643981e+01,1.063746981194344556e+02 +-9.842800072096360964e+01,-2.491721106902332394e+00,-4.059300000000001774e+01,9.845953483538603734e+01,1.064700097441810840e+02,4.066940278728695546e+01,1.064991626680726426e+02 +1.649335030491783982e+01,1.062207578381858184e+02,2.505999999999971806e+00,1.074936277181082431e+02,1.668264488265374723e+01,1.062503149723261515e+02,1.075228349514615758e+02 +5.709892092883092118e+01,6.647107061544011231e+01,6.653599999999997294e+01,8.762813475134404939e+01,8.767739769881906398e+01,9.405021278425063258e+01,1.100260391725509095e+02 +7.258636386481521185e+01,-5.448871241720553371e+01,-6.938400000000004297e+01,9.076232698647760344e+01,1.004137424614547314e+02,8.822221509849394749e+01,1.142450850408904302e+02 +3.464337842907843878e+01,-9.671418252748480882e+00,-1.088306000000000182e+02,3.596804137008204805e+01,1.142114843846289460e+02,1.092594885004483274e+02,1.146202403433178034e+02 +-1.907644091331894742e+01,-1.141717101653539430e+02,4.175540000000000873e+01,1.157544383598354045e+02,4.590668826085529730e+01,1.215676471403538557e+02,1.230552860675271205e+02 +-5.987125491380356834e+01,1.085546536775182176e+02,-7.096999999999979991e+00,1.239704803572187757e+02,6.029041859162737182e+01,1.087863973300243146e+02,1.241734569422932992e+02 +9.851720037723498535e+01,6.869505971925173071e+01,3.265899999999999181e+01,1.201026644167390316e+02,1.037894457600013851e+02,7.606327307466830234e+01,1.244638914745952718e+02 +1.228862099236030190e+02,-5.074387708045833278e+00,-1.993699999999998340e+01,1.229909346252779443e+02,1.244929899969786504e+02,2.057263667135464758e+01,1.245963641885243476e+02 +5.791162975332190399e+01,1.084447469419974368e+02,-3.137100000000000932e+01,1.229390906099426957e+02,6.586270949244222095e+01,1.128911102802778714e+02,1.268785231668452980e+02 +-1.497692698779064813e+01,6.118220049983938225e+01,-1.110509000000000412e+02,6.298864977121003506e+01,1.120562837720742095e+02,1.267894477029242211e+02,1.276709535909016608e+02 +5.818958955188281834e+01,9.333108628845729982e+01,6.920210000000001571e+01,1.099850898985860397e+02,9.041547974006769550e+01,1.161878750653159358e+02,1.299447984507654326e+02 +-5.910828942869557068e+01,4.062683990680488932e+01,1.120381000000000284e+02,7.172398483073565956e+01,1.266740925793291410e+02,1.191766586728423079e+02,1.330295675840886247e+02 +1.172978789805783890e+02,5.270424638165241049e+01,-4.283600000000001273e+01,1.285944400042254188e+02,1.248747985357430537e+02,6.791657001540878014e+01,1.355413327955732257e+02 +1.000221727798925713e+02,-8.652753869369917084e+01,-3.595700000000002206e+01,1.322552456426598155e+02,1.062889500212072136e+02,9.370123158950265463e+01,1.370560317862745023e+02 +3.870972840764807188e+01,-7.967230966029637784e+01,-1.077949000000000410e+02,8.857832692030275723e+01,1.145346390372968557e+02,1.340425954412110059e+02,1.395201077479874812e+02 +5.015125721365984646e+01,-8.407360703508311417e+01,1.056049600000000339e+02,9.789545443992915352e+01,1.169083238127733750e+02,1.349843656742927465e+02,1.439997485296479738e+02 +3.908625965839746641e+00,-5.197368859256584983e+00,2.457000000000022055e+00,6.503076195156165085e+00,4.616731088209145462e+00,5.748868763432525064e+00,6.951751505916109686e+00 +-7.116404000764305238e+00,-2.518108008203580539e+01,2.332699999999999818e+01,2.616734606336673608e+01,2.438836064400586423e+01,3.432543842542874302e+01,3.505536947458968910e+01 +-2.629930529215926072e+01,-4.433296175175200915e-01,-3.308399999999997476e+01,2.630304164920787713e+01,4.226351280774226638e+01,3.308697020202615846e+01,4.226583793088648378e+01 +-3.020793125799900736e+01,4.754039241739064892e+00,-3.554099999999999682e+01,3.057973184970705915e+01,4.664420426899781091e+01,3.585754551153766556e+01,4.688584734224164663e+01 +4.970742368248465937e+00,5.152583546443975138e+01,-3.061899999999999977e+01,5.176504612187665799e+01,3.101985558785695574e+01,5.993692418958881518e+01,6.014269000468816273e+01 +4.757239816735091154e+01,-4.017557631953968666e+01,-4.043500000000000227e+01,6.226724660686491575e+01,6.243494448137975894e+01,5.700058032693314658e+01,7.424418647274768546e+01 +2.549998086550112930e+01,-7.195950928028105409e+01,-5.846000000000003638e+00,7.634408948962440888e+01,2.616151257364382232e+01,7.219658365780790632e+01,7.656758920065185237e+01 +5.148102413319065818e+01,-4.537294517879627165e+01,-3.797799999999998022e+01,6.862215385719258620e+01,6.397362214070856368e+01,5.916952457302704005e+01,7.843040535404755076e+01 +-7.993244266581314150e-01,-7.240283889779857418e+01,-3.892999999999997840e+01,7.240725101811024444e+01,3.893820514018398882e+01,8.220532817561506533e+01,8.220921420375957211e+01 +2.940860683134087594e+01,-7.715687813953763907e+01,-3.388999999999981583e+00,8.257148418188842243e+01,2.960323422804320970e+01,7.723127064369367645e+01,8.264100266211592327e+01 +7.778032942534991889e+01,-4.492961556127875156e+01,-4.894000000000005457e+00,8.982455120956772987e+01,7.793414451648233410e+01,4.519537133915708438e+01,8.995777473904217914e+01 +3.953905227060022298e+01,9.985886713528861947e+01,-2.353000000000000114e+01,1.074017225187775182e+02,4.601084170559430930e+01,1.025936364768460010e+02,1.099490377402207315e+02 +7.698100499869178748e+01,-1.173324544590773257e+02,-4.382399999999998386e+01,1.403316785333948360e+02,8.858113854883899307e+01,1.252495422961355729e+02,1.470153834671734785e+02 +-1.222840332536579808e+02,5.680902402978881582e+01,-6.291199999999997772e+01,1.348356406889471089e+02,1.375183788909021700e+02,8.476547030021788487e+01,1.487903550099899803e+02 +4.848472252429928631e+00,7.311109137295716209e+00,1.219527999999999857e+01,8.772684879786339707e+00,1.312373946179148021e+01,1.421890189486687817e+01,1.502281113102444010e+01 +-7.087166985522344476e+00,1.679900032542718691e+00,-1.479599999999999227e+01,7.283543093853080741e+00,1.640577800290731858e+01,1.489106040949860699e+01,1.649156196362293514e+01 +-1.007586459898811881e+01,-1.533384989435944590e+01,1.361799999999999500e+01,1.834802441681090457e+01,1.694027660391476076e+01,2.050797104988075148e+01,2.284950599028103824e+01 +-1.008103075273538707e+00,2.560593150404201879e+01,2.249199999999996180e+01,2.562576828116705485e+01,2.251458051597618137e+01,3.408157555321775334e+01,3.409648169533134876e+01 +-1.658646057582882349e+01,-3.444908308745652903e+01,6.635999999999995680e+00,3.823414704161722710e+01,1.786469060559191036e+01,3.508241470546865770e+01,3.880575338786772477e+01 +-3.434883847532910295e+01,9.675603102423337987e+00,1.533199999999999363e+01,3.568557131390603132e+01,3.761532837294191012e+01,1.812974129422766723e+01,3.883980205922621565e+01 +-9.067761523714580107e+00,-4.093978139840146468e+01,-8.873999999999966803e+00,4.193196871123340230e+01,1.268748103647678249e+01,4.189049506688716917e+01,4.286067983595053477e+01 +2.129211357074382249e+01,-6.660469878087133111e+01,1.757999999999981355e+00,6.992524579863882650e+01,2.136456562416968907e+01,6.662789553700913814e+01,6.994734136477279662e+01 +-8.431016299629845889e+01,-3.446050515499518951e+01,-9.747000000000014097e+00,9.108089810712158396e+01,8.487171256350619331e+01,3.581243393763470806e+01,9.160094982586078061e+01 +-9.139732998182080337e+01,-3.278060512245247082e+01,-2.454300000000000637e+01,9.709809472899041793e+01,9.463525123761145608e+01,4.095029818199321880e+01,1.001518789089850117e+02 +-9.691478054719118518e+01,-3.696721941788200638e+01,-1.178240000000000265e+01,1.037258405605857661e+02,9.762837516967199747e+01,3.879948789932455355e+01,1.043928874481401579e+02 +1.017632527996211138e+02,4.427832855517772259e+01,2.397767999999999944e+01,1.109789619702755203e+02,1.045499342828200753e+02,5.035374383124522524e+01,1.135396808969545361e+02 +9.512397995032370090e+01,2.555481243152826210e+01,-7.325700000000000500e+01,9.849680197854230812e+01,1.200631484286064676e+02,7.758631636706689960e+01,1.227526294993316185e+02 +-4.316132828115951270e+01,-1.046261905117717106e+02,-6.191399999999998727e+01,1.131792383787749685e+02,7.547346325029756997e+01,1.215729539700567869e+02,1.290072997779566322e+02 +-1.223700456531623786e+02,3.174180093886570830e+01,3.114040000000000319e+01,1.264198164845987264e+02,1.262701571445804660e+02,4.446646420621473084e+01,1.301986732350204079e+02 +-1.175215734007324500e+02,3.905291007616142451e+01,4.333568000000000353e+01,1.238404215109125630e+02,1.252569414269969457e+02,5.833618899516106637e+01,1.312037772362629084e+02 +-1.225129042414046125e+02,-2.035800320110865869e+01,-4.423699999999999477e+01,1.241928339317527872e+02,1.302548420392254513e+02,4.869661654916437499e+01,1.318361565315066741e+02 +-1.325887688403927314e+02,-3.569185309546810458e+01,-3.061899999999999977e+01,1.373087397072727356e+02,1.360783038680711172e+02,4.702586031524012355e+01,1.406812466571134053e+02 +3.472558804828804568e+01,4.505389588815341995e+01,1.288085999999999842e+02,5.688338949113328624e+01,1.334073532428381554e+02,1.364606498909504921e+02,1.408097135639440580e+02 +3.585177655033476185e+00,7.654835150535745925e+00,3.978999999999984993e+00,8.452810183603377325e+00,5.355925673322139424e+00,8.627220999944158919e+00,9.342507211665871836e+00 +-3.844098935697184061e+01,-2.386755825922114127e+01,1.465899999999999181e+01,4.524787287818221415e+01,4.114117090874811566e+01,2.800975933951084329e+01,4.756318198985573531e+01 +6.019351231143082259e+01,4.324790255738867017e+01,3.660999999999998522e+01,7.411909335656878284e+01,7.045247351503262223e+01,5.666280239816403963e+01,8.266760006193332799e+01 +-9.863450166840266320e+01,-6.711546081660981145e+01,-2.195099999999999341e+01,1.193031852047547545e+02,1.010475695916241534e+02,7.061396095409098450e+01,1.213057970626301625e+02 +-1.327972344454428821e+02,2.845021131101589162e+01,9.379999999999881766e-01,1.358106034151950041e+02,1.328005471387747036e+02,2.846566998406073878e+01,1.358138426081795842e+02 +-1.363824121004763583e+02,2.079537616048014570e+01,-3.040999999999996817e+00,1.379587257117143224e+02,1.364163113830019824e+02,2.101654944694456617e+01,1.379922377563318889e+02 diff --git a/tests/data/Z1Z2_perpl_models_config.yaml b/tests/data/Z1Z2_perpl_models_config.yaml new file mode 100644 index 0000000..8b07561 --- /dev/null +++ b/tests/data/Z1Z2_perpl_models_config.yaml @@ -0,0 +1,226 @@ +############# VALUES FOR MODEL CONFIGURATION ######################################## +# Read comments around the values and set as required +# Lists of values where indicated generate multiple models for fitting and comparison +# Optionally generate few models or one model first for testing rough parameter +# values and ranges +##################################################################################### + + +# Model direction +## xx, yy or zz (distances in 1D); +## xy, xz or yz (distances across 2D plane); +### or xyz (3D Euclidean distances) +model_direction: xz + + +# Dimensionality of the space being modelled +## 1 for distances along a single direction +## 2 for distances across a plane +## Not yet implemented for KDE: 3 for distances in a 3D volume +# +## Used to calculate different KDE kernel shapes +## appropriate to different dimensionality +#dimension: [2] + + +# Limits for distance components in x, y, and z +## +## One use is to restrict modelling to the use of pairs of points +## with relative postitions close to the modelling direction. +## +## Use "null" to apply no limits +limits: + x_limit: 200. # e.g. 10 or null + y_limit: 200. + z_limit: 200. + + +#### FROM HERE, INCLUDE ALL DESIRED VALUES IN A LIST #### +#### FOR A SWEEP THROUGH MULTIPLE MODEL FORMS ########### + + +## EXPERIMENTAL DATA PREPARATION PARAMETERS ############### +## DO NOT COMPARE INFORMATION CRITERIA BETWEEN DIFFERENT ## +## bin_size, kde_kernel_size or fitlength values ########## + +# Bin size for histogram analysis (not used for KDE) +bin_sizes: [5] # e.g. [1, 2, 5] + + +# kernel size for building experimental KDE +## sqrt(2) * average or max estimated localisation precision +## are reasonable to try and several values can be tried +## +## Smaller values keep more detail, larger values smoothen noise +kde_kernel_size: [6.693] # e.g. [5, 6, 7] + + +# Maximum distance for the fit (less than relpos_filter if using a KDE) +fitlength: [60] # Can try multiple fit lengths, e.g. [90, 100, 110] + + +## MODEL PARAMETERS ##################################### +## Define the form of the model ######################### +## Include all desired values each parameter in a list ## + +# Background model options +## From: +## null, flat, linear, linear_flat, trans_bg_0, trans_bg_1, trans_bg_2 +# +# TODO: Change model names and add more models +background: [null, flat, linear, trans_bg_0, trans_bg_1, trans_bg_2] # e.g. [flat, linear] + + +# Number of characteristic distance peaks in the model +n_peaks: [0, 2, 3] # Can try models with different numbers of peaks + + +# Relative amplitudes of peaks +# +## Indepedent amplitude ("indep_amp") for characteristic distance peaks +## or amplitudes decrease linearly ("lin_decrease") (1, 1 - 1/n, 1 - 2/n, ...), +## as used for 1D analysis through a finite space. +## +## So "indep_amp", "lin_decrease", or both in the list +peak_amps: [indep_amp] # e.g. [indep_amp, lin_decrease] + + +# Relationship between the characteristic distances: +# +## Distance ratios: +## Multipliers for fixed relationships between distances +## +### "int_ratios": Distances are at integer multiples of the first distance +## +### "custom_ratios": Distances are at other multiples of the first distance +#### For custom_ratios, also include these distance ratios +#### in `custom_ratios_list` below +## +### "indep_dists": Distances are at independent (multiple) distances +dist_ratios: [int_ratios, indep_dists, custom_ratios] # int_ratios, custom_ratios, indep_dists +# +# +## Ratios between characteristic distances +### null if dist_ratios is not custom_ratios +### or must be list with length corresponding to max of n_peaks +custom_ratios_list: [1, 1.4142, 2] # null or e.g. [1.0, 0.7, 0.4] for max n_peaks: 3 + + +# Model repeated localisations of single underlying objects +## These localisations have a spread determined by the experimental precision +repeats: [True, False] # True, False or both + + +# Include offset - Not implemented yet +offset: [False] # True, False (or both) + + +# Normalise the model RPD +## This can account for finding more pairs and larger distances in 2D and 3D +## e.g. divides counts by distance for 2D distances +## Use this if you normalised (sometimes called standardised) +## the experimental RPD +normalise: [True] # True, False + + +#### PARAMETER VALUE INITIALISATION AND ALLOWED RANGE FOR FIT #### + +# Initial guesses +params_initial: + # amplitude of peak 1 + amp_peak_1: 10 + # amplitude of peak 2 + amp_peak_2: 10 + # amplitude of peak 3 + amp_peak_3: 10 + # etc... + amp_peak_4: 10 + amp_peak_5: 10 + amp_peak_6: 10 + + # First characteristic distance + ## Can use a single value or a list corresponding to + ## a list of different charac_dist relationships used. + ## e.g. [20, 10, 24] would be the distances used for models with + ## dist_ratios of [int_ratios, custom_ratios, indep_dists], + ## if that was the list specified (the lists must have the same length). + characteristic_distance_1: [20, 15, 20] + # Second characteristic distance + characteristic_distance_2: 40 # Can be list as for distance_1 + characteristic_distance_3: 55 + # etc. + + # Sigma for characteristic distances + characteristic_distance_broadening: 5 + + # Standard deviation of localisation precision + ## (Spread of repeated localisations) + loc_prec_sd: 5 + + # amplitude of peak due to repeated localisations + loc_prec_amp: 10 + + # Background slope + bg_slope: 0.1 + + # Background offset + bg_offset: 0 + + # Quadratic and cubic background model coefficients + trans_bg_quad_coeff: -0.1 + trans_bg_cubic_coeff: 0.1 + + +# Parameter lower bounds +params_lower: + amp_peak_1: 0 + amp_peak_2: 0 + amp_peak_3: 0 + amp_peak_4: 0 + amp_peak_5: 0 + amp_peak_6: 0 + # etc... + + # Can use lists as for initial guesses + characteristic_distance_1: [5, 5, 5] + characteristic_distance_2: 20 + characteristic_distance_3: 50 + # etc. + + characteristic_distance_broadening: 0 + + loc_prec_sd: 0 + loc_prec_amp: 0 + + bg_slope: 0 + bg_offset: 0 + + trans_bg_quad_coeff: -10 + trans_bg_cubic_coeff: 0 + +# upper bounds for parameters +params_upper: + amp_peak_1: 100 + amp_peak_2: 100 + amp_peak_3: 100 + amp_peak_4: 100 + amp_peak_5: 100 + amp_peak_6: 100 + # etc... + + # Can use lists as for initial guesses + characteristic_distance_1: [40, 20, 40] + characteristic_distance_2: 50 + characteristic_distance_3: 60 + # etc. + + characteristic_distance_broadening: 15 + + loc_prec_sd: 10 + loc_prec_amp: 100 + + bg_slope: 10 + bg_offset: 10 + + trans_bg_quad_coeff: 0 + trans_bg_cubic_coeff: 10 diff --git a/tests/test_utils.py b/tests/test_utils.py index d8d06d8..eb16a4e 100644 --- a/tests/test_utils.py +++ b/tests/test_utils.py @@ -21,6 +21,7 @@ """ import unittest + import perpl.io.utils as ut @@ -33,7 +34,6 @@ class TestFindExponent(unittest.TestCase): of significant figures needed. """ - def test_positve_exponent_float(self): """ Test the exponent for a large number in float format. The exponent will @@ -80,7 +80,7 @@ def test_positive_exponent_complex(self): """ print("Start TestFindExponent test_positve_exponent_complex", flush=True) - value = 2.088696E+4 + value = 2.088696e4 result = ut.find_exponent(str(value)) @@ -93,7 +93,7 @@ def test_negative_exponent_complex(self): """ print("Start TestFindExponent test_negative_exponent_complex", flush=True) - value = 6.0E-2 + value = 6.0e-2 result = ut.find_exponent(str(value)) @@ -106,7 +106,7 @@ def test_find_exponent_size_zero_complex(self): """ print("Start TestFindExponent test_zero_exponent_complex", flush=True) - value = 1.06E0 + value = 1.06e0 result = ut.find_exponent(str(value)) @@ -136,8 +136,13 @@ def test_larger_pos_expo_value_smaller_neg_expo_uncertainty_floats(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_larger_pos_expo_value_smaller_neg_" - "expo_uncertainty_floats"), flush=True) + print( + ( + "Start TestPlusAndMinus test_larger_pos_expo_value_smaller_neg_" + "expo_uncertainty_floats" + ), + flush=True, + ) value = 86.96 @@ -145,11 +150,11 @@ def test_larger_pos_expo_value_smaller_neg_expo_uncertainty_floats(self): result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "86.960".strip()) self.assertIn(str(result2), "0.01") @@ -168,8 +173,13 @@ def test_larger_pos_expo_value_smaller_pos_expo_uncertainty_floats(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_larger_pos_expo_value_smaller_pos_" - "expo_uncertainty_floats"), flush=True) + print( + ( + "Start TestPlusAndMinus test_larger_pos_expo_value_smaller_pos_" + "expo_uncertainty_floats" + ), + flush=True, + ) value = 86.96 @@ -177,11 +187,11 @@ def test_larger_pos_expo_value_smaller_pos_expo_uncertainty_floats(self): result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "87.00".strip()) self.assertIn(str(result2), "10.00") @@ -199,8 +209,13 @@ def test_smaller_pos_expo_value_larger_pos_expo_uncertainty_floats(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_smaller_pos_expo_value_larger_pos_" - "expo_uncertainty_floats"), flush=True) + print( + ( + "Start TestPlusAndMinus test_smaller_pos_expo_value_larger_pos_" + "expo_uncertainty_floats" + ), + flush=True, + ) value = 86.96 @@ -208,11 +223,11 @@ def test_smaller_pos_expo_value_larger_pos_expo_uncertainty_floats(self): result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "86.96".strip()) self.assertIn(str(result2), "300.0") @@ -231,8 +246,13 @@ def test_larger_neg_expo_value_smaller_neg_expo_uncertainty_floats(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_larger_neg_expo_value_smaller_neg_" - "expo_uncertainty_floats"), flush=True) + print( + ( + "Start TestPlusAndMinus test_larger_neg_expo_value_smaller_neg_" + "expo_uncertainty_floats" + ), + flush=True, + ) value = 0.8696 @@ -240,11 +260,11 @@ def test_larger_neg_expo_value_smaller_neg_expo_uncertainty_floats(self): result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "0.870".strip()) self.assertIn(str(result2), "0.011") @@ -263,8 +283,13 @@ def test_smaller_neg_expo_value_larger_neg_expo_uncertainty_floats(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_smaller_neg_expo_value_larger_neg_" - "expo_uncertainty_floats"), flush=True) + print( + ( + "Start TestPlusAndMinus test_smaller_neg_expo_value_larger_neg_" + "expo_uncertainty_floats" + ), + flush=True, + ) value = 0.08696 @@ -272,11 +297,11 @@ def test_smaller_neg_expo_value_larger_neg_expo_uncertainty_floats(self): result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "0.09".strip()) self.assertIn(str(result2), "0.101") @@ -295,8 +320,13 @@ def test_smaller_neg_expo_value_larger_neg_expo_uncertainty_floats1(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_smaller_neg_expo_value_larger_neg_" - "expo_uncertainty_floats1"), flush=True) + print( + ( + "Start TestPlusAndMinus test_smaller_neg_expo_value_larger_neg_" + "expo_uncertainty_floats1" + ), + flush=True, + ) value = 0.08696 @@ -304,16 +334,15 @@ def test_smaller_neg_expo_value_larger_neg_expo_uncertainty_floats1(self): result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "0.09".strip()) self.assertIn(str(result2), "0.330") - def test_smaller_neg_expo_value_larger_pos_expo_uncertainty_floats(self): """ Tests that when there is a fraction float value (smaller than 1.0 so with @@ -328,8 +357,13 @@ def test_smaller_neg_expo_value_larger_pos_expo_uncertainty_floats(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_smaller_neg_expo_value_larger_pos_" - "expo_uncertainty_floats"), flush=True) + print( + ( + "Start TestPlusAndMinus test_smaller_neg_expo_value_larger_pos_" + "expo_uncertainty_floats" + ), + flush=True, + ) value = 0.08696 @@ -337,16 +371,15 @@ def test_smaller_neg_expo_value_larger_pos_expo_uncertainty_floats(self): result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "0.09".strip()) self.assertIn(str(result2), "33.000") - def test_larger_pos_expo_value_smaller_zero_expo_uncertainty_floats(self): """ Tests that when there is a large float (bigger than 10 and so with a @@ -361,8 +394,13 @@ def test_larger_pos_expo_value_smaller_zero_expo_uncertainty_floats(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_larger_pos_expo_value_smaller_zero_" - "expo_uncertainty_floats"), flush=True) + print( + ( + "Start TestPlusAndMinus test_larger_pos_expo_value_smaller_zero_" + "expo_uncertainty_floats" + ), + flush=True, + ) value = 86.96 @@ -370,11 +408,11 @@ def test_larger_pos_expo_value_smaller_zero_expo_uncertainty_floats(self): result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "87.0".strip()) self.assertIn(str(result2), "1.10") @@ -392,8 +430,13 @@ def test_larger_zero_expo_value_smaller_pos_expo_uncertainty_floats(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_larger_zero_expo_value_smaller_pos_" - "expo_uncertainty_floats"), flush=True) + print( + ( + "Start TestPlusAndMinus test_larger_zero_expo_value_smaller_pos_" + "expo_uncertainty_floats" + ), + flush=True, + ) value = 8.696 @@ -401,17 +444,15 @@ def test_larger_zero_expo_value_smaller_pos_expo_uncertainty_floats(self): result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "8.70".strip()) self.assertIn(str(result2), "1.10") - - def test_larger_pos_expo_value_smaller_neg_expo_uncertainty_complex(self): """ Tests that when there is a large complex (bigger than 10 and so with a @@ -426,20 +467,25 @@ def test_larger_pos_expo_value_smaller_neg_expo_uncertainty_complex(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_larger_pos_expo_value_smaller_neg_" - "expo_uncertainty_complex"), flush=True) + print( + ( + "Start TestPlusAndMinus test_larger_pos_expo_value_smaller_neg_" + "expo_uncertainty_complex" + ), + flush=True, + ) - value = 8.696E+1 + value = 8.696e1 - uncertainty = 1.1E-2 + uncertainty = 1.1e-2 result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "86.960".strip()) self.assertIn(str(result2), "0.01") @@ -458,20 +504,25 @@ def test_larger_pos_expo_value_smaller_pos_expo_uncertainty_complex(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_larger_pos_expo_value_smaller_pos_" - "expo_uncertainty_complex"), flush=True) + print( + ( + "Start TestPlusAndMinus test_larger_pos_expo_value_smaller_pos_" + "expo_uncertainty_complex" + ), + flush=True, + ) - value = 8.696E+1 + value = 8.696e1 - uncertainty = 1.00E+1 + uncertainty = 1.00e1 result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "87.00".strip()) self.assertIn(str(result2), "10.00") @@ -489,20 +540,25 @@ def test_smaller_pos_expo_value_larger_pos_expo_uncertainty_complex(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_smaller_pos_expo_value_larger_pos_" - "expo_uncertainty_complex"), flush=True) + print( + ( + "Start TestPlusAndMinus test_smaller_pos_expo_value_larger_pos_" + "expo_uncertainty_complex" + ), + flush=True, + ) - value = 8.696E+1 + value = 8.696e1 - uncertainty = 3.0001E+2 + uncertainty = 3.0001e2 result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "86.96".strip()) self.assertIn(str(result2), "300.0") @@ -521,20 +577,25 @@ def test_larger_neg_expo_value_smaller_neg_expo_uncertainty_complex(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_larger_neg_expo_value_smaller_neg_" - "expo_uncertainty_complex"), flush=True) + print( + ( + "Start TestPlusAndMinus test_larger_neg_expo_value_smaller_neg_" + "expo_uncertainty_complex" + ), + flush=True, + ) - value = 8.696E-1 + value = 8.696e-1 - uncertainty = 1.1E-2 + uncertainty = 1.1e-2 result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "0.870".strip()) self.assertIn(str(result2), "0.011") @@ -553,20 +614,25 @@ def test_smaller_neg_expo_value_larger_neg_expo_uncertainty_complex(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_smaller_neg_expo_value_larger_neg_" - "expo_uncertainty_complex"), flush=True) + print( + ( + "Start TestPlusAndMinus test_smaller_neg_expo_value_larger_neg_" + "expo_uncertainty_complex" + ), + flush=True, + ) - value = 8.696E-3 + value = 8.696e-3 - uncertainty = 1.01E-1 + uncertainty = 1.01e-1 result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "0.01".strip()) self.assertIn(str(result2), "0.1010") @@ -585,25 +651,29 @@ def test_smaller_neg_expo_value_larger_neg_expo_uncertainty_complex1(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_smaller_neg_expo_value_larger_neg_" - "expo_uncertainty_complex1"), flush=True) + print( + ( + "Start TestPlusAndMinus test_smaller_neg_expo_value_larger_neg_" + "expo_uncertainty_complex1" + ), + flush=True, + ) - value = 8.696E-2 + value = 8.696e-2 - uncertainty = 3.30001E-1 + uncertainty = 3.30001e-1 result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "0.09".strip()) self.assertIn(str(result2), "0.330") - def test_smaller_neg_expo_value_larger_pos_expo_uncertainty_complex(self): """ Tests that when there is a fraction float value (smaller than 1.0 so with @@ -618,25 +688,29 @@ def test_smaller_neg_expo_value_larger_pos_expo_uncertainty_complex(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_smaller_neg_expo_value_larger_pos_" - "expo_uncertainty_complex"), flush=True) + print( + ( + "Start TestPlusAndMinus test_smaller_neg_expo_value_larger_pos_" + "expo_uncertainty_complex" + ), + flush=True, + ) - value = 8.696E-2 + value = 8.696e-2 - uncertainty = 3.30001E+1 + uncertainty = 3.30001e1 result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "0.09".strip()) self.assertIn(str(result2), "33.000") - def test_larger_pos_expo_value_smaller_zero_expo_uncertainty_complex(self): """ Tests that when there is a large float (bigger than 10 and so with a @@ -651,20 +725,25 @@ def test_larger_pos_expo_value_smaller_zero_expo_uncertainty_complex(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_larger_pos_expo_value_smaller_zero_" - "expo_uncertainty_complex"), flush=True) + print( + ( + "Start TestPlusAndMinus test_larger_pos_expo_value_smaller_zero_" + "expo_uncertainty_complex" + ), + flush=True, + ) - value = 8.696E+1 + value = 8.696e1 - uncertainty = 1.1E-0 + uncertainty = 1.1e-0 result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "87.0".strip()) self.assertIn(str(result2), "1.10") @@ -682,20 +761,25 @@ def test_larger_zero_expo_value_smaller_pos_expo_uncertainty_complex(self): so we just check that the returned strings are in rather than equal to each other. """ - print(("Start TestPlusAndMinus test_larger_zero_expo_value_smaller_pos_" - "expo_uncertainty_complex"), flush=True) + print( + ( + "Start TestPlusAndMinus test_larger_zero_expo_value_smaller_pos_" + "expo_uncertainty_complex" + ), + flush=True, + ) - value = 8.696E+0 + value = 8.696e0 - uncertainty = 1.1E-0 + uncertainty = 1.1e-0 result1, result2 = ut.plus_and_minus(value, uncertainty) - #print(result1) - #print(value) + # print(result1) + # print(value) - #print(result2) - #print(uncertainty) + # print(result2) + # print(uncertainty) self.assertIn(str(result1).strip(), "8.70".strip()) self.assertIn(str(result2), "1.10") @@ -714,41 +798,45 @@ def test_linux_filenames(self): """ print("Start TestPrimaryFilenameAndPathSetup test_linux_filename", flush=True) - prog = 'relative_positions' - prog_short_name = 'rp' - description = 'Calculating the relative positions of points as vectors.' - info = {'prog':prog, 'prog_short_name':prog_short_name, 'description':description} - info['start'] = '2019-11-07_15-49-55' - info['dims'] = 2 - info['filter_dist'] = 200 - info['zoom'] = 10 - info['verbose'] = True - info['in_file_and_path'] = ("/localhome/joanna/PERPL_data/Nup107_SNAP_3D" - "_GRROUPED_10nmZprec.csv") - info['colours_analysed'] = None - info['nns'] = 0 - info['bin_size'] = 1 + prog = "relative_positions" + prog_short_name = "rp" + description = "Calculating the relative positions of points as vectors." + info = { + "prog": prog, + "prog_short_name": prog_short_name, + "description": description, + } + info["start"] = "2019-11-07_15-49-55" + info["dims"] = 2 + info["filter_dist"] = 200 + info["zoom"] = 10 + info["verbose"] = True + info["in_file_and_path"] = ( + "/localhome/joanna/PERPL_data/Nup107_SNAP_3D" "_GRROUPED_10nmZprec.csv" + ) + info["colours_analysed"] = None + info["nns"] = 0 + info["bin_size"] = 1 ut.primary_filename_and_path_setup(info) - result1 = info['results_dir'] - result2 = info['in_file_no_extension'] - result3 = info['in_file_no_path'] - result4 = info['short_results_dir'] - result5 = info['short_filename_without_extension'] + result1 = info["results_dir"] + result2 = info["in_file_no_extension"] + result3 = info["in_file_no_path"] + result4 = info["short_results_dir"] + result5 = info["short_filename_without_extension"] self.assertIn( result1, - '/localhome/joanna/PERPL_data/PERPL/' - 'rel_posns_filter200_2D_bin1_2019-11-07_15-49-55' + "/localhome/joanna/PERPL_data/PERPL/" + "rel_posns_filter200_2D_bin1_2019-11-07_15-49-55", ) - self.assertIn(result2, 'Nup107_SNAP_3D_GRROUPED_10nmZprec') - self.assertIn(result3, 'Nup107_SNAP_3D_GRROUPED_10nmZprec.csv') + self.assertIn(result2, "Nup107_SNAP_3D_GRROUPED_10nmZprec") + self.assertIn(result3, "Nup107_SNAP_3D_GRROUPED_10nmZprec.csv") self.assertIn( - result4, - '/localhome/joanna/PERPL_data/PERPL/rel_posns_f200_2D_b1' + result4, "/localhome/joanna/PERPL_data/PERPL/rel_posns_f200_2D_b1" ) - self.assertIn(result5, 'Nup107') + self.assertIn(result5, "Nup107") class TestSecondaryFilenameAndPathSetup(unittest.TestCase): @@ -758,35 +846,40 @@ class TestSecondaryFilenameAndPathSetup(unittest.TestCase): rot_2d_symm script that analyses experimental data. """ - def test_linux_filenames(self): """ Test the input expected from a Linux system. """ print("Start TestSecondaryFilenameAndPathSetup test_linux_filename", flush=True) - prog = 'relative_positions' - prog_short_name = 'rp' - description = 'Calculating the relative positions of points as vectors.' - info = {'prog':prog, 'prog_short_name':prog_short_name, 'description':description} - info['start'] = '2019-11-07_16-49-55' - info['dims'] = 2 - info['filter_dist'] = 200 - info['zoom'] = 10 - info['verbose'] = True - info['in_file_and_path'] = ('/localhome/joanna/PERPL_data/PERPL_relative' - '_positions/Nup107_SNAP_3D_GRROUPED_10nmZprec/' - 'filter_200_2019-11-07_15-49-55' - '/Nup107_SNAP_3D_GRROUPED_10nmZprec_PERPL' - '-relpos_200.0filter.csv') + prog = "relative_positions" + prog_short_name = "rp" + description = "Calculating the relative positions of points as vectors." + info = { + "prog": prog, + "prog_short_name": prog_short_name, + "description": description, + } + info["start"] = "2019-11-07_16-49-55" + info["dims"] = 2 + info["filter_dist"] = 200 + info["zoom"] = 10 + info["verbose"] = True + info["in_file_and_path"] = ( + "/localhome/joanna/PERPL_data/PERPL_relative" + "_positions/Nup107_SNAP_3D_GRROUPED_10nmZprec/" + "filter_200_2019-11-07_15-49-55" + "/Nup107_SNAP_3D_GRROUPED_10nmZprec_PERPL" + "-relpos_200.0filter.csv" + ) ut.secondary_filename_and_path_setup(info) - result1 = info['results_dir'] - result2 = info['in_file_no_extension'] - result3 = info['in_file_no_path'] - result4 = info['short_results_dir'] - result5 = info['short_filename_without_extension'] + result1 = info["results_dir"] + result2 = info["in_file_no_extension"] + result3 = info["in_file_no_path"] + result4 = info["short_results_dir"] + result5 = info["short_filename_without_extension"] print("result1: ", result1) print("result2: ", result2) @@ -794,20 +887,26 @@ def test_linux_filenames(self): print("result4: ", result4) print("result5: ", result5) - self.assertIn(result1, '/localhome/joanna/PERPL_data/PERPL_relative_' - 'positions/Nup107_SNAP_3D_GRROUPED_10nmZprec/filter_200' - '_2019-11-07_15-49-55/relative_positions/filter_200_' - '2019-11-07_16-49-55') - self.assertIn(result2, 'Nup107_SNAP_3D_GRROUPED_10nmZprec_PERPL-relpos_200') - self.assertIn(result3, 'Nup107_SNAP_3D_GRROUPED_10nmZprec_PERPL-relpos_200.0filter.csv') - self.assertIn(result4, '/localhome/joanna/PERPL_data/PERPL_relative' - '_positions/Nup107_SNAP_3D_GRROUPED_10nmZprec/' - 'filter_200_2019-11-07_15-49-55' - '/rp/Nup10-s-s_200/f_200_2019-11-07_16-49-55') - self.assertIn(result5, 'Nup10-s-s_200') - + self.assertIn( + result1, + "/localhome/joanna/PERPL_data/PERPL_relative_" + "positions/Nup107_SNAP_3D_GRROUPED_10nmZprec/filter_200" + "_2019-11-07_15-49-55/relative_positions/filter_200_" + "2019-11-07_16-49-55", + ) + self.assertIn(result2, "Nup107_SNAP_3D_GRROUPED_10nmZprec_PERPL-relpos_200") + self.assertIn( + result3, "Nup107_SNAP_3D_GRROUPED_10nmZprec_PERPL-relpos_200.0filter.csv" + ) + self.assertIn( + result4, + "/localhome/joanna/PERPL_data/PERPL_relative" + "_positions/Nup107_SNAP_3D_GRROUPED_10nmZprec/" + "filter_200_2019-11-07_15-49-55" + "/rp/Nup107/f_200_2019-11-07_16-49-55", + ) + self.assertIn(result5, "Nup107") -if __name__ == '__main__': +if __name__ == "__main__": unittest.main() - \ No newline at end of file