a /'a@sddlZddlmZddlZddlZddlZddlmZmZddl m Z m Z m Z m Z mZmZmZddl mZddl mZdZeeddZed d d gZed gd ZedgdZejereeeddZne Zded<e!dge!dge!dgd\Z"Z#de_$gdZ%gdZ&gdZ'dde(e%e&e'DZ)d\Z*Z+dde,d DZ-d!d"Z.d#d$Z/d%d&Z0d'd(Z1d)d*Z2d+d,Z3d-d.Z4d>d/d0Z5d1d2Z6d3d4Z7d5d6Z8d7d8Z9d9d:Z:d;d<Z;eAs r0)& .>r1cCsg|]}t|ddttqS)r%r$)rmom_atolmom_rtol)r-ir*r*r/r0KscCs@ttdd }tjt|ddWdn1s20YdS)z3Writes the current res_dict to the target JSON filewr)indentN)openresults_filepathjsondumpres_dict)fr*r*r/ write_dataQsr?cCs t|S)zConverts a namedtuple to a dict)dict_asdict)Z named_tupler*r*r/to_dictWsrBcCst|jt|jdS)z;Formats an MPResult namedtuple into a dict for JSON dumping)rr)rBrfloatr)rr*r*r/mp_res_to_dict\srDcshtttfddfddfdd}t|dtgt tgdd d }|S) z7Straightforward implementation of studentized range CDFcs(t|t||t|dSNrphiPhiszrrr*r/innerjszcdf_mp..innercs&|dtt|||SrErGrrIrMnur*r/outermszcdf_mp..outercsBtdtdtdddd||SNr7rrrrrIrrPrQr*r/wholeps zcdf_mp..wholergauss-legendrer$methodZ maxdegreer rr rrrPrUresr*rMrrPrQrr/cdf_mpfsr]cshtttfddfddfdd}t|dtgt tgdd d }|S) z7Straightforward implementation of studentized range PDFcs8t||t|t||t|dSNr7rFrIrLr*r/rM}szpdf_mp..innercs"|tt|||SNrNrIrOr*r/rQszpdf_mp..outercsJtdtddtdddd||SrRrSrIrTr*r/rUs "zpdf_mp..wholerrVr$rWrYrZr*r\r/pdf_mpysr`csztttfddfddfddfdd}t|d tgd tgt tgd d d }|S) z.Implementation of the studentized range momentcs8t|||t|t|||t|dSr^rFrrJrK)rr*r/rMszmoment_mp..innercs$|tt||||Sr_rNrarOr*r/rQszmoment_mp..outercsLtdtddtdddd|||SrRrSrarTr*r/pdfs " zmoment_mp..pdfcs||||Sr_r*ra)rrbr*r/rUszmoment_mp..wholerrVr$rWrY)rrrPrUr[r*)rMrrrPrQrbr/ moment_mps rccsD|tvrt|dt|ttfddt|}t|dkS)znSearches the results dict for a result in the set that matches a case. Returns True if such a case exists.z not present in data structure!cs |dkS)Nrr*)r[Z case_dictr*r/zresult_exists..r)r= ValueErrorrBlistfilterlen)set_keyr.Z existing_resr*rdr/ result_existss rlcCsLt}||}td|dd|dt|d||tt||fS)zARuns the single passed case, returning an mp dictionary and indexz Finished r/z in batch. (Took zs))time perf_counterprintrDr)r. run_lambdarkindex total_casest_startr[r*r*r/runs  rucCs$|\}}}t|||tdS)zSA callback for completed jobs. Inserts and writes a calculated result to file.N)r=insertr?)r[rrrkZ result_dictr*r*r/ write_results rwcs~fddtD}ttt|dtdg}tt}|D]}||jt|tdqN| | dS)z)Runs an array of cases and writes to filecs,g|]$\}}t|s||tfqSr*)rlrj)r-rrr.casesrqrkr*r/r0s zrun_cases..rmz? cases won't be calculated because their results already exist.)argscallbackN) enumeraterprjr num_poolsappendZ apply_asyncrurwclosejoin)ryrqrkZjob_argZjobsZpoolr.r*rxr/ run_casess$rcCst|j|j|jSr_)r`rrrr.r*r*r/run_pdfsrcCst|j|j|jSr_)r]rrrrr*r*r/run_cdfsrcCst|j|j|jSr_)rcrrrrr*r*r/ run_momentsrcCst}dtttt}td|dtdttdtttdtdttdtttd td ttd ttt d td t|ddS)Nr7z Processing z test caseszRunning 1st batch (z' PDF cases). These take about 30s each.rzRunning 2nd batch (z' CDF cases). These take about 30s each.rzRunning 3rd batch (zH moment cases). These take about anywhere from a few hours to days each.r zTest data generated in rJ) rnrorj cdf_pdf_cases moment_casesrprrrr)rtrsr*r*r/mains   r__main__)rr)= itertools collectionsrr;rnosZmultiprocessingrrZmpmathrrrrr r r r rGr rHr:maxr}rrrpathisfileloadr9r=r@ setdefaultr+r,ZdpsZcp_qZcp_kZcp_nuproductrr2r3rangerr?rBrDr]r`rcrlrurwrrrrr__name__r*r*r*r/sh $