ó î&]\c@sòddlZddlZyddlmZWnDek rrddlZddlZddlZded„ZnXdfd„ƒYZ de fd„ƒYZ d e fd „ƒYZ d e fd „ƒYZ d fd„ƒYZdfd„ƒYZdS(iÿÿÿÿN(t lru_cacheiÿcs)dtfd„ƒY‰‡‡‡fd†S(s Thanks to ilialuk @ https://stackoverflow.com/users/2121105/ilialuk for this code snippet. Modifications by S. Endres t LruCacheClasscBsAeZd„Zdd„Zd„Zd„Zee_d„ZRS(cSs(||_||_||_i|_dS(N(t _input_funct _max_sizet_timeoutt _caches_dict(tselft input_functmax_sizettimeout((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyt__init__s   cSs?||jkr;|j|=tjƒtjƒg|j|;scSstt|ƒƒt|ƒS(N(RR(R((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR >sii(tjointmaptsortedRR R R RRRtlenRtpopitemtFalsetKeyErrorR(RRRRt kwargs_keytkeytcur_caller_cache_dict((RsElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR7s0       N(t__name__t __module__R RRRRR(((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyRs    cstj|ƒˆ|ˆˆƒƒS(N(RR(R(RtmaxsizeR (sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR as (tobject(R-R ((RR-R sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR sNtComplexcBsÝeZdedddd„Zd„Zeed„Zd„Zd„Zd„Z d„Z d„Z d„Z d „Z ed „Zd „Zed dƒd „ƒZed dƒd„ƒZed dƒd„ƒZd„ZRS(cCsù||_||_||_d|_d|_g|_t|||||ƒ|_|j|d|ƒ|rvd|_ n |j ƒ|jj gƒ|jdj |j ƒ|j j ƒ|_d|_|jƒg|_|jj dƒ|jj dƒdS(Nitsymmetryi(tdimtboundsR0tgent perm_cycletHt VertexCachetVtn_cubetgeneration_cyclet add_centroidtappendtC0thomology_group_rankthgrthgrdt graph_mapt performance(RR1tfunct func_argsR0R2tg_constg_args((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR fs&           cCs|jS(N(R5(R((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyRsc Csfttj|dtƒƒ}||_ttj|dtƒƒ}||_t|ƒ}t|ƒ}|g}|rÔtddd|j ƒ|_ |j j |j |ƒd} |j | ||ƒ|j j |j |ƒn_tdd||ƒ|_ |j j |j |ƒ|j j |j |ƒg} |j| ||ƒ|rbdGHx!|j ƒD]} | jƒqKWndS(sw Generate the simplicial triangulation of the n dimensional hypercube containing 2**n vertices tdtypeisInitial hyper cube:N(tlisttnptzerostinttorigintonestsupremumttupletSimplexR1R<t add_vertexR7t perm_symmetrytCelltpermt print_out( RR1R0tprintoutRKRMt origintuplet supremumtuplet x_parentsti_st i_parentstv((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR8“s,     c Cst|ƒ}gt|jƒD]}||kr|^q}xÜ|D]Ô}tj|ƒ}|j|ƒtj|ƒ} d| |RUtC_newtV_newRbR[tt1tt2tvecR|R}((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyRvVs0&#   +c Csy|j|Wn!tk r2|jjgƒnX|ƒ}|dj}|dj}|jttj|ƒtj|ƒdƒ}|j|j|j|ƒx&|D]}|j |j|jƒq¦Wt ||j |j |j ƒ}|j|dƒ|j|ƒx"|dd!D]}|j|ƒqWt ||j |j |j ƒ} | j|| j dƒx"|dd!D]}| j|ƒqlWxMt|dd!ƒD]8\} }| |j krÂ| j|ƒq—| j|ƒq—W| j|dƒ|j|j|ƒ|j|j| ƒdS(sI Split a hypersimplex S into two sub simplices by building a hyperplane which connects to a new vertex on an edge (the longest edge in dim = {2, 3}) and every other vertex in the simplex that is not connected to the edge being split. This function utilizes the knowledge that the problem is specified with symmetric constraints The longest edge is tracked by an ordering of the vertices in every simplices, the edge between first and second vertex is the longest edge to be split in the next iteration. iiÿÿÿÿg@iN(R5RtR;RR7RNRHRhRjR^RORwR9R1RPRu( RtSR3tstfirstxtlastxR‡R[tS_new_ltS_new_utk((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyRŒs:    -   R-cCs8|j||ƒ}|j||ƒ}||}t|ƒS(s- Use the origin and supremum vectors to find a new cell in that subspace direction NOTE: NOT CURRENTLY IN USE! Parameters ---------- origin : tuple vector (hashable) supremum : tuple vector (hashable) Returns ------- (RƒR„RN(RRKRMtv_x_tRˆR‰RŠ((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pytgenerate_sub_cell_2Ìs cCs$tj|ƒ}||tj|ƒS(N(RHRh(RRKtv_xtv_o((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyRƒâscCs tj|ƒ}|tj|ƒS(N(RHRh(RRMR”tv_s((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR„èsc CsJddlm}|jdkrØ|jƒx|jD]ú}xñ|D]é}xà|ƒD]Õ}|jdkr€tj|j dt ƒ}nstj|j dt ƒ}xXt t |jƒƒD]A}|||j|d|j|d|j|d||_sN(RwRtsumR¥(R((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR=ZscCs,|jdkr%|j|j|_n|jS(su Returns the difference between the current homology group of the cell and it's parent group N(R³RRwR²thgd(R((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pythomology_group_differentialcscCsdS(s– Returns the number of stationary points theoretically contained in the cell based information currently known about the cell N((R((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pytpolytopial_sperner_lemmamscCs"x|ƒD]}|jƒq WdS(s3 Print the current cell to console N(RT(RR[((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyRTts( R+R,R RRPR=R¸R¹RT(((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR°Gs   RRcBseZdZd„ZRS(sR Contains a cell that is symmetric to the initial hypercube triangulation cCs8tt|ƒj||ƒ||_||_d|_dS(N(tsuperRRR RKRMRRi(RR±R²RKRM((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR s  (R+R,t__doc__R (((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyRR|sROcBseZdZd„ZRS(sq Contains a simplex that is symmetric to the initial symmetry constrained hypersimplex triangulation cCs2tt|ƒj||ƒ|d|d|_dS(Ni(RºROR R9(RR±R²R9R1((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR s(R+R,R»R (((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyROŠstVertexcBsSeZdddddddd„Zd„Zd„Zd„Zd„Zd„ZRS( c Csj||_t|ƒ|_tj|dtƒ} |dk ryx=t|ƒD],\} \} } | | | | | | | ˆjr7t‡fd†ˆjDƒƒˆ_tˆ_nˆjS(sBCheck whether this vertex is strictly less than all its neighboursc3s!|]}ˆj|jkVqdS(N(RÁ(RµR[(R(sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pys Ös(RÄtallRoRÌR&(R((RsElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR´Ós " cCs[dj|jƒGHd}x'|jD]}|dj|jƒ7}q!W|GHdj|jƒGHdS(Ns Vertex: {}s Connections: s{} s Order = {}(R…RRoR½(Rtconstrtvc((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyRTÛs N((( R+R,RR RÊR^RjR´RT(((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR¼–s $   R6cBs,eZdddded„Zed„ZRS(cCs^i|_||_||_||_||_||_d|_d|_|rZd|_ndS(Niiÿÿÿÿ( tcacheRBRDRÅRCR2tnfevtsizeRp(RRBRCR2RDRÅtindexed((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR æs        cCsJy|j|SWn4tk rE|ry|jd7_t|d|jd|jd|jd|jd|jd|jƒ}n9t|d|jd|jd|jd|jd|jƒ}||j|<|jdk r:|jdk r|j r|j d7_ |j d7_ q7|j d7_ q:|j d7_ |j d7_ n|j|SXdS(NiR2RBRCRDRÅRp( RÑR'RpR¼R2RBRCRDRÅRR¾RÒRÓ(RRRÔtxval((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyt __getitem__õs0       (N((R+R,RRkR RÖ(((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyR6ås (tnumpyRHR]RRt ImportErrorR R RR/R.R°RRROR¼R6(((sElib/python2.7/site-packages/scipy/optimize/_shgo_lib/triangulation.pyts      Xÿã5 O