o Nrf#@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 einsumchunkr_/var/www/html/software/conda/envs/catlas/lib/python3.10/site-packages/dask/array/einsumfuncs.py chunk_einsum s   rcCst|dkr tdt|dtr9|ddd}dd|ddD}|D]}|d vr-q&|tvr7td |q&nt|}g}g}tt|d D]}|| d|| dqIt|rd|d nd}d d|D}d}t|d}t |D])\} } | D]}|t ur|d7}qt|t r|t |7}qtd| |kr|d7}qy|dur|d7}|D]}|t ur|d7}qt|t r|t |7}qtdd|vsd|vr|ddkp|ddk} | s|ddkrtdd|vr|dddddd} ttt| } d| }d}d|vr!|d\}}|d}d}n|d}d}t |D]g\} } d| vr| ddksE| ddkrItd|| jdkrTd}nt|| jd}|t| d8}||krk|}|dkrttd|dkr| dd|| <q,|| d}| d||| <q,d|}|dkrd}n|| d}|r|d|d|7}n>d}|dd}tt|D]}|tvrtd |||dkr||7}qdtt|t|}|d||7}d|vr|d\}}n+|}|dd}d}tt|D]}|tvrtd |||dkr*||7}q|D]}||vr;td|q.t|dt|krMtd|||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 cSg|]}t|qSrr.0orrr 9z&parse_einsum_input..r Nz.,->z#Character %s is not a valid symbol.cSrrr)rvrrrr Kr!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                             rIF)r optimize split_everyc Os |}t|\}}}d||f} |durtjdd|D}|dur7dd|D} tj| g| Rd|i\}} dd|d D}d d |D} | t|} t| }tt t |t | gd d t ||DRdd| D|| |||d|}|dkrt|}|j t t||||dS|S)aDask added an additional keyword-only argument ``split_every``. split_every: int >= 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 rrrrr zeinsum..FcSs$g|]}tj|jd|jdqS)r)r)npZ broadcast_tor rrrrrrr s$rJcSrr)tuple)rirrrr r!r%cSsh|] }|D]}|qqSrr)rrOarrr szeinsum..css|] }|D]}|VqqdS)Nr)rZaprPrrr szeinsum..cSsi|]}|dqSr r)rindrrr rLzeinsum..)Z adjust_chunksr rr r rJr)ZaxisrK)rIr;rMZ result_typeZ einsum_pathr<r:r+rrrNzipsumr1r2)r rJrKrrZ einsum_dtypeZinputsZoutputsopsrZfake_opsrAZall_indsZ contract_indsr resultsizerrrrsB    r) __future__rnumpyrMZdask.array.corerrrZ dask.utilsrr7r:r0rrIrrrrrs  +