gG! ZdZddlmZddlmZddlmZddlmZm Z m Z m Z m Z ddl ZddlmZGddeZeGd d Zd ed efd Zded efdZded e efdZded e eeffdZddde eefdej0ej2ded e eej0effdZy)z Classes and functions for reshaping Stan output. Especially with the addition of tuples, Stan writes flat arrays of data with a rich internal structure. ) dataclass)Enum)prod)AnyDictIterableListTupleNceZdZdZdZdZy) VariableTypeN)__name__ __module__ __qualname__SCALARCOMPLEXTUPLEK/var/www/dash_apps/app1/venv/lib/python3.12/site-packages/stanio/reshape.pyr r s FG Err c0eZdZUdZeed<eed<eed<eedfed<eed<e ded<dd e d e jfd Z d e efd Zd efd Zd efdZdde j"ded e j"fdZdde j"de d ej(efdZy)Variablea This class represents a single output variable of a Stan model. It contains information about the name, dimensions, and type of the variable, as well as the indices of where that variable is located in the flattened output array Stan models write. Generally, this class should not be instantiated directly, but rather created by the :func:`parse_header()` function. name start_idxend_idx. dimensionstypecontentstopreturnc|jtjk(r[t|jDcgc]%\}}t |dz|j df'}}}tj |}n[|jtjk(rtj}n-|jtjk(rtj}|rStj |jfScc}}w)Nr F)r!) rr r enumerater strdtypenprfloat64r complex128r)selfr!iparameltsr&s rr&zVariable.dtype2s 99 ** *!*$-- 8AuQUU[[U[34DHHTNE YY,-- -JJE YY,.. .MME L88UDOO45 5s*C9cBt|j|jSN)rangerrr*s rcolumnszVariable.columnsCsT^^T\\22rc,t|jSr/)rrr1s rnum_eltszVariable.num_eltsFsDOO$$rc4|j|jz Sr/)rrr1s relt_sizezVariable.elt_sizeIs||dnn,,rsrcoffsetc J |j|z}|j|z}|jtjk(r'|d||fj dg|j ddiS|jtjk(rr|d||fj ddg|j ddi}|dddddfd|dddddfzz}|jj dg|j ddiS|jtjk(rtjt|jddt|j ft}t|j!D]}||j#z|j!z}|j$D cgc]} | j'|||z } } t| d jd D] t) fd | D| |f<|j dg|j ddiSycc} w) N.orderFry?r )r&)r8rc3(K|] }| ywr/r).0eltr+s r z+Variable._extract_helper..bs'?3A'?s)rrrr rreshaperrsqueezerr'emptyrshapeobjectr0r4r6r _extract_helpertuple) r*r7r8startendretoutidxoffr,r-r+s @rrFzVariable._extract_helperNs 'llV# 99 ++ +.3sE#I~&..rODOOO3O O YY,.. .-#c59n%--b!QdooQSQCa1f+SADqD\ 11C(3;;=((IdooISI I YY,,, , hhciin%tDOO'<=VCT]]_- @DMMO+t}}>"&))#eck)BtAw}}Q/0@A"''?$'?"?C3K@  @3;;r?DOO?3? ?- s%H rEc |j|}|s|j|j}|jdkDr0|jg|j dd|j ddi}|S|jd}|S)a Given an array where the final dimension is the flattened output of a Stan model, (e.g. one row of a Stan CSV file), extract the variable and reshape it to the correct type and dimensions. This will most likely result in copies of the data being made if the variable is not a scalar. Parameters ---------- src : np.ndarray The array to extract from. Indicies besides the final dimension are preserved in the output. object : bool If True, the output of tuple types will be an object array, otherwise it will use custom dtypes to represent tuples. Returns ------- npt.NDArray[Any] The extracted variable, reshaped to the correct dimensions. If the variable is a tuple, this will be an object array, otherwise it will have a dtype of either float64 or complex128. r Nr:r;r<r)axis)rFastyper&ndimrArDrrB)r*r7rErKs rextract_reshapezVariable.extract_reshapees8""3'**TZZ\*C 88a<#++Ksyy"~KKsKC ++1+%C rN)T)r)rrr__doc__r%__annotations__intr r r boolr'r&rr2r4r6ndarrayrFnptNDArrayrrRrrrrrs  IN Lc3h :666"3#3%#%-#- @2::@s@2::@.$2::$t$s{{SVGW$rrtupr"c2d|jdddzS)Ndummy_:r split)rZs r_munge_first_tupler`s ciiQ'* **rr,cN|jddjddS)N.rr]r^)r,s r_get_base_namercs& ;;s A  $ $S )! ,,rheadercV|jdz}|jd}g}d}t|d}tdt |dz D]V}||}t||dz}||k7s d|vr}|jddd}d|vsd|vrt j } |dd }nt j} |jt|||dzttt|| g n|jddjddd}djtjtt |||dz} |jt|||dzttt|t j"t%|  |dz}|}Y|S) Nz,__dummy,rr r]rbz.realz.imagr:)rrrrrr )stripr_rcr0lenr rrappendrrGmaprUjoindictfromkeysr`r _from_header) rdentriesparamsrrr+entry next_namedimsr munged_headers rrnrns \\^j (Fll3G FI '!* %D 1c'lQ& '( "71q5>2  %{{3'+e#w%'7'//D9D'..D !"+ !A#(S$#8!!#  {{3'*005ab9 #MM#&8')aRSe:T"UV!  !"+ !A#(S$#8)//!-m!<  AIDQ(T MrcTt|Dcic]}|j|c}Scc}w)a Given a comma-separated list of names of Stan outputs, like that from the header row of a CSV file, parse it into a dictionary of :class:`Variable` objects. Parameters ---------- header : str Comma separated list of Stan variables, including index information. For example, an ``array[2] real foo` would be represented as ``foo.1,foo.2``. Returns ------- Dict[str, Variable] A dictionary mapping the base name of each variable to a :class:`Variable`. )rnr)rdr,s r parse_headerrvs&$,8+? @%EJJ  @@ @s%TrE parameterssourcerEc|jDcic] }|j|j||"c}Scc}w)a Given a dictionary of :class:`Variable` objects and a source array, extract the variables from the source array and reshape them to the correct dimensions. Parameters ---------- parameters : Dict[str, Variable] A dictionary of :class:`Variable` objects, like that returned by :func:`parse_header()`. source : npt.NDArray[np.float64] The array to extract from. object : bool If True, the output of tuple types will be an object array, otherwise it will use custom dtypes to represent tuples. Returns ------- Dict[str, npt.NDArray[Any]] A dictionary mapping the base name of each variable to the extracted and reshaped data. rw)valuesrrR)rxryrEr,s rstan_variablesr|sF< &&(    E))&)@@  s%;)rS dataclassesrenumrmathrtypingrrrr r numpyr' numpy.typingrXr rr%r`rcrnrvrYr(rVr|rrrrs "334  rr rj+C+C+-#-#-11h1hAAc8m!4A2 S(]# KK #    #s{{3  r