l jddlmZddlmZddlZerddlmZmZGddZ GddZ Gd d Z Gd d Z Gd dZ dS)) annotations) TYPE_CHECKINGN)PointND PointND_Arrayc&eZdZd dZddZdd Zd S)QuickHullPoint coordinatesrreturnNonec||_dSN)r selfr s L/home/agentuser/manim-venv/lib/python3.11/site-packages/manim/utils/qhull.py__init__zQuickHullPoint.__init__ s&intcNt|jSr )hashr tobytesrs r__hash__zQuickHullPoint.__hash__sD$,,..///rotherobjectboolc|t|tsttj|j|j}|Sr ) isinstancer ValueErrornp array_equalr )rrare_coordinates_equals r__eq__zQuickHullPoint.__eq__s?%00  &(n  e/' ' %$rNr rr r r rrrr r__name__ __module__ __qualname__rrr"rrrr sP''''0000%%%%%%rrc&eZdZd dZddZdd Zd S)SubFacetr rr r cP||_td|D|_dS)Nc34K|]}t|VdSr )r).0cs r z$SubFacet.__init__..s*GGaq 1 1GGGGGGr)r frozensetpointsrs rrzSubFacet.__init__s,&GG;GGGGG rrc*t|jSr )rr3rs rrzSubFacet.__hash__!sDK   rrrrcZt|tst|j|jkSr )rr,rr3rrs rr"zSubFacet.__eq__$s)%**  {el**rNr#r$r%r&r*rrr,r,sTHHHH!!!!++++++rr,c.eZdZddZddZdd ZddZdS)Facetr rinternalrr r c|_tj|d_|_t fdtjjdD_ dS)Nraxisc3jK|]-}ttjj|dV.dS)rr;N)r,rdeleter )r/irs rr1z!Facet.__init__../sS# #  RYt/;;; < <# # # # # # r) r rmeancentercompute_normalnormalr2rangeshape subfacets)rr r9s` rrzFacet.__init__+s&!w{;;; ))(33 "# # # # 4+1!455# # #   rc|j|jz }tj|\}}}|dddf}|tj|z}tj||j|z dkr|dz}|S)Nr)r rArlinalgsvdnormdot)rr9centered_vhrCs rrBzFacet.compute_normal4s#dk19==**1bRU)")..((( 6&$+0 1 1A 5 5 bLF rrc*t|jSr )rrFrs rrzFacet.__hash__@sDN###rrrrcZt|tst|j|jkSr )rr8rrFr6s rr"z Facet.__eq__Cs)%''  ~00rN)r rr9rr r )r9rr rr$r%)r'r(r)rrBrr"r*rrr8r8*sd        $$$$111111rr8ceZdZddZdS)Horizonr r c:t|_g|_dSr )setfacetsboundaryrs rrzHorizon.__init__Js"%%% (* rN)r r )r'r(r)rr*rrrSrSIs(++++++rrScDeZdZdZdddZdd Zdd ZddZddZddZ dS) QuickHulla QuickHull algorithm for constructing a convex hull from a set of points. Parameters ---------- tolerance A tolerance threshold for determining when points lie on the convex hull (default is 1e-5). Attributes ---------- facets List of facets considered. removed Set of internal facets that have been removed from the hull during the construction process. outside Dictionary mapping each facet to its outside points and eye point. neighbors Mapping of subfacets to their neighboring facets. Each subfacet links precisely two neighbors. unclaimed Points that have not yet been classified as inside or outside the current hull. internal An internal point (i.e., the center of the initial simplex) used as a reference during hull construction. tolerance The tolerance used to determine if points are considered outside the current hull. h㈵> tolerancefloatr r cg|_t|_i|_i|_d|_d|_||_dSr )rVrUremovedoutside neighbors unclaimedr9r[)rr[s rrzQuickHull.__init__js;#% #&55 QS 57/3(, "rr3rcxtj}|||jd|jddzd}||_tj|d}||_t|jdD]W}ttj ||d|}| ||j |X|j D]F}|jD]<}|j|t#|=GdS)NrF)replacer;r9)rrandom default_rngchoicerErar@r9rDr8r>classifyrVappendrFr` setdefaultrUadd) rr3rngsimplex new_internalr0facetfsfs r initializezQuickHull.initializess7i##%% JJv|A Q!(;UJ K K   "a 8 8 8 $ w}Q'(( & &A")GQQ777,OOOE MM% K  u % % % % <  5 5!!! 2 2Hx0E7:??AAHw~--d6M6MXw77- ''111trc |j\}}|dks ||dzkrtd|dkrtd|||jJ|jJ d}|jD]}||jvr |j|\}}|nd}|||}|jD]} |j| d} | Jtj |j| f|_|j | | j D]G} |j | | |j | tkr|j | =H|jD]} t#tj | j|f|j} || |j| | j D]<} |j | t | =|sdS)Nrrcz0Not enough points supplied to build Convex Hull!z*The Convex Hull of 1D data is its min-max!TFre)rErrsrar9rVr^r_rrvstackrlrFr`discardrUrWr8r rirjrk)rr3numdimupdatedrpr_rzr~rqpoints_to_appendrrnfnsfs rbuildzQuickHull.buildsI<S 1HH#a--OPP P !88IJJ J ~)))}((( G J JDL((#|E2 ?"G"223>>G$^ 7 7+/<?1+=(/;;;)+DNDT3U)V)V ((+++"#+77B N2.66q999#~b1SUU::$(N2$67&.JJ"Ir~s&;<rs1"""""" 433333333 % % % % % % % % + + + + + + + +11111111>++++++++ LLLLLLLLLLr