a tDf#@snddlmZddlZddlmZmZmZddlm Z dZ e e Z ddZ dd Ze edd dd d d ZdS)) annotationsN)asarray blockwise einsum_lookup) derived_from4abcdefghijklmnopqrstuvwxyzABCDEFGHIJKLMNOPQRSTUVWXYZcOs^|d}|d}|d}tt|d}||g|Rd|i|}||jd|S)N subscriptsncontract_inds kernel_dtyperdtype)poprdispatchtypeZreshapeshape)operandskwargsrr r einsumchunkrc/nfs/NAS7/SABIOD/METHODE/ermites/ermites_venv/lib/python3.9/site-packages/dask/array/einsumfuncs.py chunk_einsum s    rcCst|dkrtdt|dtrt|ddd}dd|ddD}|D]"}|d vrZqL|tvrLtd |qLn0t|}g}g}tt|d D]$}|| d|| dqt|r|d nd}d d|D}d}t|d}t |D]Z\} } | D]<}|t ur|d7}n"t|t r2|t |7}ntdq| |kr|d7}q|dur|d7}|D]<}|t ur~|d7}n"t|t r|t |7}ntdqfd|vsd|vr|ddkp|ddk} | s|ddkrtdd|vr|dddddd} ttt| } d| }d}d|vr`|d\}}|d}d}n|d}d}t |D]\} } d| vrv| ddks| ddkrtd|| jdkrd}n t|| jd}|t| d8}||kr|}|dkrtdn:|dkr$| dd|| <n|| d}| d||| <qvd|}|dkr`d}n|| d}|r|d|d|7}n|d}|dd}tt|D]4}|tvrtd |||dkr||7}qdtt|t|}|d||7}d|vr |d\}}nV|}|dd}d}tt|D]4}|tvrZtd |||dkr@||7}q@|D]}||vrztd|qzt|dt|krtd|||fS)a A reproduction of numpy's _parse_einsum_input() which in itself is a reproduction of c side einsum parsing in python. Returns ------- input_strings : str Parsed input strings output_string : str Parsed output string operands : list of array_like The operands to use in the numpy contraction Examples -------- The operand list is simplified to reduce printing: >> a = np.random.rand(4, 4) >> b = np.random.rand(4, 4, 4) >> __parse_einsum_input(('...a,...a->...', a, b)) ('za,xza', 'xz', [a, b]) >> __parse_einsum_input((a, [Ellipsis, 0], b, [Ellipsis, 0])) ('za,xza', 'xz', [a, b]) rzNo input operands cSsg|] }t|qSrr.0orrr 9z&parse_einsum_input..r Nz.,->z#Character %s is not a valid symbol.cSsg|] }t|qSrr)rvrrrrKr z...z=For this input type lists must contain either int or Ellipsis,->->z%Subscripts can only contain one '->'..TFzInvalid Ellipses.rzEllipses lengths do not match.z/Output character %s did not appear in the inputzDNumber of einsum subscripts must be equal to the number of operands.)len ValueError isinstancestrreplaceeinsum_symbols_setlistrangeappendr enumerateEllipsisinteinsum_symbols TypeErrorcountsetjoinsplitrmaxndimsorted)rrsZ tmp_operandsZ operand_listZsubscript_list_Z output_listlastnumsubinvalidusedZunusedZ ellipse_indslongestZ input_tmpZ output_subZsplit_subscriptsZout_subZ ellipse_countZrep_indsZ out_ellipseZoutput_subscriptZtmp_subscriptsZ normal_indsZinput_subscriptscharrrrparse_einsum_inputs                                 rHF)r optimize split_everyc Os"|}t|\}}}d||f} |dur= 2 or dict(axis: int), optional Determines the depth of the recursive aggregation. Defaults to ``None`` which would let dask heuristically decide a good default. r%NcSsg|] }|jqSr)r rrrrrr zeinsum..FcSs$g|]}tj|jd|jdqS)r)r)npZ broadcast_tor rrrrrrrr rIcSsg|] }t|qSr)tuple)rirrrrr r$cSsh|]}|D]}|q qSrr)rrMarrr r zeinsum..css|]}|D] }|Vq qdS)Nr)rZaprNrrr r zeinsum..cSsi|] }|dqSr r)rindrrr r zeinsum..)Z adjust_chunksr rr r rIr)ZaxisrJ)rHr:rKZ result_typeZ einsum_pathr;r9r*rrrLzipsumr0r1)r rIrJrrZ einsum_dtypeinputsoutputsopsrZfake_opsr@Zall_indsZ contract_indsr resultsizerrrrsB     r) __future__rnumpyrKZdask.array.corerrrZ dask.utilsrr6r9r/rrHrrrrrs  +