B T\@sddlmZmZmZddlZddlZddlZddlmZddl m Z m Z ddl m Z mZddlZddlmZmZmZddlmZdd lmZdd lmZdd lmZdd lmZd dlm Z d dl!m"Z"m#Z#m$Z$m%Z%d dlm&Z&m'Z'd dl(m)Z)d dl*m+Z+m,Z,d dlm-Z-m.Z.m/Z/m0Z0m1Z1m2Z2m3Z3m4Z4m5Z5m6Z6m7Z7m8Z8m9Z9d dl:m;Z;e ej<dddZddZ>e ej?ddZ?e ej@ddZ@e ejAdd d!ZAe ejBdd"d#ZBe ejCdd$d%ZCe ejDd&d'ZDe ejEdd(d)ZEd*d+ZFe ejGd,d-ZGe ejHd.d/ZHd0ZIeIJZKd1d2ZLe ejMdd3d4ZMe ejNd5d6ZNe ejOd7d8ZOe&ejPd9d:ZPe ejQd;d<ZQd=d>ZRe ejSd?d@ZSe ejTdAdBZTe ejUddCdDZUe ejVddFdGZVe ejWddHdIZWdJdKZXe ejYdLdMZYe ejZddNdOZZe ej[ddPdQZ[ddRdSZ\e ej]ddTdUZ]e ej^ddVdWZ^e ej_ddXdYZ_ddZd[Z`e ejadd\d]Zadd^d_Zbe&eced`dddad`Zde ejeddbdcZee ejfdddeZfe ejgddfdgZge ejhddhdiZhe ejidjdkZie ejjddldmZjdndoZke ejlddpdqZldrdsZmdtduZndvdwZoe ejpddzd{Zpe ejqdd|d}Zqd~dZre ejsddZsddZtejuetevgdZtddZwe ejxddZxe ejydddZye ejzdddZze ej{ddZ{e ej|ddZ|ddZ}ddZ~e ejdddZe ejddZe e jdddZdddZe ejddZdddZe ejdddZdS))divisionprint_functionabsolute_importN) LooseVersion)wrapspartial)RealIntegral)concatsliding_window interleave)Iterable)flatten)tokenize)HighLevelGraph)funcname)chunk)arangediagemptyindices) safe_wraps validate_axis)ones)multiplysqrt) Array map_blockselemwise from_arrayasarray asanyarray concatenatestack blockwisebroadcast_shapesis_scalar_for_elemwise broadcast_totensordot_lookup)einsumcCsHx&|dk r&|j|kr&|dddf}qW|dk rD|j|krD||}|S)N)ndimdtypeastype)xr-ndminr12lib/python3.7/site-packages/dask/array/routines.pyarrays  r3cGsdd|D}tj|S)NcSsg|]}t|r|n|jqSr1)r(r-).0ar1r1r2 *szresult_type..)np result_type)argsr1r1r2r8(sr8cGsg}xp|D]h}t|}|jdkr*|d}n>|jdkrH|ddddf}n |jdkrh|dddddf}||q Wt|dkr|dS|SdS)Nr)NNNrr )r#r,appendlen)arysnew_arysr/r1r1r2 atleast_3d.s      r>cGsng}xL|D]D}t|}|jdkr*|d}n|jdkrD|dddf}||q Wt|dkrf|dS|SdS)Nr)NNr)r#r,r:r;)r<r=r/r1r1r2 atleast_2dBs     r?cGsRg}x0|D](}t|}|jdkr(|d}||q Wt|dkrJ|dS|SdS)Nrr)r#r,r:r;)r<r=r/r1r1r2 atleast_1dTs   r@FcCs tdd|D}t|d|dS)Ncss|]}t|VqdS)N)r?)r4r/r1r1r2 fszvstack..r)axisallow_unknown_chunksizes)tupler$)tuprCr1r1r2vstackdsrFcCs2tdd|Dr t|d|dSt|d|dSdS)Ncss|]}|jdkVqdS)rN)r,)r4r/r1r1r2rAlszhstack..r)rBrCr)allr$)rErCr1r1r2hstackjsrHcCs tdd|D}t|d|dS)Ncss|]}t|VqdS)N)r>)r4r/r1r1r2rAtszdstack..r )rBrC)rDr$)rErCr1r1r2dstackrsrIc Csr||kr |S|dkr||j}|dkr0||j}tt|j}t|}||||<||<ttj||||||jdS)Nr)axis1axis2r-)r,listranger&r7swapaxesr-)r5rJrKindoutr1r1r2rNxs  rNcsl|rt|jkr4tdnttjddd}tfdd|D}ttj|ttjj|dS)Nzaxes don't match arrayc3s$|]}|dkr|jn|VqdS)rN)r,)r4d)r5r1r2rAsztranspose..)r-axes) r;r, ValueErrorrDrMr&r7 transposer-)r5rSr1)r5r2rUs rUcCslt|}|jtdg}ytddd||<Wn.tk rZtdt|t|jfYnXt|}||S)z Reverse element order along axis. Parameters ---------- axis : int Axis to reverse element order of. Returns ------- reversed array : ndarray NrQz#`axis` of %s invalid for %s-D array)r#r,slice IndexErrorrTstrrD)mrBslr1r1r2flipsr[cCs t|dS)Nr)r[)rYr1r1r2flipudsr\cCs t|dS)Nr)r[)rYr1r1r2fliplrsr]Zabcdefghijklmnopqrstuvwxyzcs tgddd}tt|}tfdddD}tfdddD}t|dkr||kr|dkrt t fd dt j Dfd dt j D}n|d }t ddg|j}x tdD]|dqW|t |}|S) NcSs|jS)N)Z__array_priority__)r/r1r1r2sz_tensordot..)keycsg|]}j|qSr1)shape)r4i)r5r1r2r6sz_tensordot..rcsg|]}j|qSr1)r`)r4ra)br1r2r6srcs g|]\}}|dkr|qS)rr1)r4ras)rSr1r2r6scs g|]\}}|dkr|qS)rr1)r4rarc)rSr1r2r6s)rS)maxr*dispatchtyper7r3r;rGminzerosrD enumerater`rVr,sortedinsert)r5rbrSr/ tensordotZa_dimsZb_dimsrOr1)r5rSrbr2 _tensordots$  rmc Cs"t|tr|\}}n.tt|jd|j|dd}ttd|}t|trR|f}t|trb|f}t|trtt|}t|trt|}t|j |j }tt|j}tt|j|j|j}||}x0t ||D]"\} } | || || || <qWt t ||||||||fd} | j|d} | S)NrrQr)r-rS)rB) isinstancerrDrMr,r rLr7Z promote_typesr-zipremover&rmsum) ZlhsZrhsrSZ left_axesZ right_axesdtZ left_indexZ right_indexZ out_indexlrZ intermediateresultr1r1r2rls2         rlcCs"t|||jdf|jdffdS)Nrr )rS)rlr,)r5rbr1r1r2dotsrvcCst||S)N)rvconjravel)r5rbr1r1r2vdotsryc CsPt|}t|}|jdks$|jdkr,tdd}|jdkrPd}|tjddf}d}|jdkrtd}|ddtjf}|j|jkr||j|jtjf}n$|j|jkr||j|jtjf}ttjttd|jd|ttd|jd|jddf|ttd|jdd|jft ||dd}|r>|ddddf}|rL|d}|S) Nrz"`matmul` does not support scalars.FrT)r-r$.).r) r#r,rTr7newaxisr&matmulrDrMr8)r5rbZa_is_1dZb_is_1drPr1r1r2r{s6    "r{c CsB|}|}t|j|jj}ttjd|d|d|dS)NZijraj)r-)rr7outerr-rfr&)r5rbr-r1r1r2r}&sr}cCstj|||f||S)N)r7apply_along_axis)arrfunc1d func1d_axis func1d_args func1d_kwargsr1r1r2_inner_apply_along_axis0src Ost|}|j|t|jd|}tjd|jd}t||f||}ttjtdkrt|jdkj t ddkrt d| |j d||j||df|j |dd}|jtt|d|j|j d||j|j |dd|tt|||jd||||d }|S)N)r)r-z1.13.0rzdNo more than one non-trivial dimension allowed in result. Need NumPy 1.13.0+ for this functionality.z -along-axis) namer-chunksZ drop_axisnew_axisrrrr)r"r`r;r7rr-r3r __version__rqintrTrechunkrrrrrLrMr,)rrBrr9kwargsZ test_dataZ test_resultrur1r1r2r~:s. 6 $r~cCst|}y t|}Wntk r.|f}YnX|jtdf}|}xT|D]L}t|||d}|j|jdkr||d|d}qJ|j|jkrJtdqJW|S)Nrr)NzJfunc must either preserve dimensionality of the input or reduce it by one.)r"rD TypeErrorr,rVr~rT)funcr5rSrZrurar1r1r2apply_over_axesds     rcCs|j|d|j|dS)N)rB)rdrg)r5rBr1r1r2ptpsrrQcCst|}t|}t|}|jtdg}|jtdg}tdd||<tdd||<t|}t|}|}x t|D]}||||}qrW|S)NrrQ)r"rr,rVrDrM)r5nrBZsl_1Zsl_2rtrar1r1r2diffsrcCsnt|}|}|dd|dd}|g}|dk rHt|g|}|dk rb|t|g}t|}|S)NrrQ)r"rr$)ZaryZto_endZto_beginZaryfrtr1r1r2ediff1dsrcCsH||}|dk r,||d||d|}tj||fd|i|}|S)aU x: nd-array array of one block coord: 1d-array or scalar coordinate along which the gradient is computed. axis: int axis along which the gradient is computed array_locs: actual location along axis. None if coordinate is scalar grad_kwargs: keyword to be passed to np.gradient NrrrB)r7gradient)r/Zblock_idcoordrB array_locs grad_kwargsZ block_locZgradr1r1r2_gradient_kernels rc sZtt|dd|d<|ddkr2tdd}|dd}|dkrZttj}nt |t rnd}|f}t |j}t |t t |krtdtfd d |D}|d krd }t |dkrt ||}t |t |krtd tjjtjt fr tn(tjjtr2jjdkr2tg}xt|D]\}xBjD]4}t||ddkrXtdt|qXWt||rd}nbt ||trtdtj} t | d} | | d} | dd8<d| d<| | f}|!j"t#jfddtjDd||||dqBW|rV|d}|S)NZ edge_orderrr z+edge_order must be less than or equal to 2.FrBTzduplicate axes not allowedc3s|]}|jVqdS)N)r,)r4ax)fr1r2rAszgradient..r1)rzGSpacing must either be a single scalar, or a scalar / 1d-array per axiszcChunk size must be larger than edge_order + 1. Minimum chunk for aixs {} is {}. Rechunk to proceed.z(dask array coordinated is not supported.rQrcsi|]}|krdnd|qS)rrr1)r4r|)rr1r2 szgradient..Znone)r-ZdepthboundaryrrBrr)$r"mathZceilgetrTpoprDrMr,rnr rr;setr issubclassr-rfr7Zbool8r.floatritemsizerirrgformatisscalarrNotImplementedErrorr3cumsumr:Z map_overlapr) rZvarargsrZdrop_result_listrBZresultsracrrZarray_loc_stopZarray_loc_startr1)rrr2rsp         rcs dkrtdjdkstdk r6jjks6tt}d|dk rfddtD}tjdgdgdj }n*fddtD}tgj }d|tj t |d f|d f<ff}t j |dkrgngd }t|||S) Nz.Must specify minlength argument in da.bincountrz bincount-cs0i|](\}}tjj|fj|ff|fqSr1)r7bincountr)r4ra_) minlengthrweightsr/r1r2rszbincount..)rcs*i|]"\}}tjj|fdf|fqS)N)r7rr)r4rar)rrr/r1r2rsz bincount-sum-r) dependencies)rr,AssertionErrorrrri __dask_keys__r7rr-rqrLrfrom_collectionsr)r/rrtokendskr-rgraphr1)rrrr/r2rs&  "rcCs2t|}tjdg|ddj}|jtj|||dS)NrF)right)r-binsr)r7r"digitizer-r)r5rrr-r1r1r2r-s rcsts"dksdkr"td|dk r>|j|jkr>tdts}\}}||krl|d8}|d7}tj||dddn}t|||||} ttt| } d| tdff} d | dfd d |dkrfd d t t| D} t gdj } n>t| }t| }fdd t t ||D} |j } tj| |dkrn|gn||gd}t|| | d}|jdd}|dk r|rttt|jd}|||fS|fSn<|rttt|jd}|||fS|fSdS)aq Blocked variant of :func:`numpy.histogram`. Follows the signature of :func:`numpy.histogram` exactly with the following exceptions: - Either an iterable specifying the ``bins`` or the number of ``bins`` and a ``range`` argument is required as computing ``min`` and ``max`` over blocked arrays is an expensive operation that must be performed explicitly. - ``weights`` must be a dask.array.Array with the same block structure as ``a``. Examples -------- Using number of bins and range: >>> import dask.array as da >>> import numpy as np >>> x = da.from_array(np.arange(10000), chunks=10) >>> h, bins = da.histogram(x, bins=10, range=[0, 10000]) >>> bins array([ 0., 1000., 2000., 3000., 4000., 5000., 6000., 7000., 8000., 9000., 10000.]) >>> h.compute() array([1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000, 1000]) Explicitly specifying the bins: >>> h, bins = da.histogram(x, bins=np.array([0, 5000, 10000])) >>> bins array([ 0, 5000, 10000]) >>> h.compute() array([5000, 5000]) Nzudask.array.histogram requires either specifying bins as an iterable or specifying both a range and the number of binsz.block_histcs"i|]\}}|f|dfqS)rr1)r4rak)rrrMr1r2rxszhistogram..rcs(i|] \}\}}||f|dfqS)rr1)r4rarw)rrrMr1r2r~s)r)r-)rB)r)NN)r7iterablerTrZlinspacerr;rLrrrirr-rorrrrqr!rr.r)r5rrMZnormedrZdensityZ bin_tokenZmnZmxrZnchunksrrr-Za_keysZw_keysrmappedrZdbr1)rrrrMr2r4sL%    $   rc CsH|dk r|t|krtdt|}|dkr    rcCsPt|||}|jdkr||St|}||jddf}t|}|||jS)Nr1rr)rr`rreshaperr)r/rrrrRZsqr_dr1r1r2corrcoefs  rcCs|jtj||jdS)N)decimalsr-)rr7roundr-)r5rr1r1r2rsrc Cs|dk }|dk }t|}d|jfg}|r:|dtjf|rN|dtjf|rb|dtjftj|j|d}||d<|rtjt|tjd|d<|s|rxlt |dD]\\} } || k} |r|| j d|d| | dd |r|| j d|d| | dd qW|S) a Helper/wrapper function for :func:`numpy.unique`. Uses :func:`numpy.unique` to find the unique values for the array chunk. Given this chunk may not represent the whole array, also take the ``indices`` and ``counts`` that are in 1-to-1 correspondence to ``ar`` and reduce them in the same fashion as ``ar`` is reduced. Namely sum any counts that correspond to the same value and take the smallest index that corresponds to the same value. To handle the inverse mapping from the unique values to the original array, simply return a NumPy array created with ``arange`` with enough values to correspond 1-to-1 to the unique values. While there is more work needed to be done to create the full inverse mapping for the original array, this provides enough information to generate the inverse mapping in Dask. Given Dask likes to have one array returned from functions like ``blockwise``, some formatting is done to stuff all of the resulting arrays into one big NumPy structured array. Dask is then able to handle this object and can split it apart into the separate results on the Dask side, which then can be passed back to this function in concatenated chunks for further reduction or can be return to the user to perform other forms of analysis. By handling the problem in this way, it does not matter where a chunk is in a larger array or how big it is. The chunk can still be computed on the same way. Also it does not matter if the chunk is the result of other chunks being run through this function multiple times. The end result will still be just as accurate using this strategy. Nvaluesrinversecounts)r-Tr)rrP) r7uniquer-r:intprr`rr;rirgrq) arrrreturn_inverse return_index return_countsurrrtravrYr1r1r2_unique_internals,!   "&rcCs~|}|dg}d|jfg}|rX|t|jdtj|jdddg|dtjfn|ddg|r|t |jdftj|jdddg|dtjfn|ddgt t df||dd}t d d |jD|_ |dg}|r||dn |d|r||dn |dd |j}|dft ft d d |D|fi} d|jfg}|rz|dtjf|r|d tjf|r|dtjfdd|D} tj|| | d} tjff} t| || |}|dg} |r| |d|rF|dddf|ddddfktj}| ||d jdd|rZ| |dt| dkrr| d} nt | } | S)Nrarr)r-rrrF)r-rcss|]}tjft|VqdS)N)r7nanr;)r4rr1r1r2rA5szunique..zunique-aggregate-css*|]"}t|drtj|fn|VqdS)rN)hasattrr7r$r)r4or1r1r2rAOsrcSsg|]}t|dr|qS)r)r)r4rr1r1r2r6^szunique..)rr)rB)rxr-extendrr`r7rrr:rr&rrDZ_chunksrrrrrr.rqr;)rrrrr9Z out_dtyperPZ out_partsrrrrrruZmtchesr1r1r2rsn          , rcCs*tj|||d}||jd|jS)N) assume_unique)r)r7Zin1drxrr`r,)element test_elementsrrr1r1r2 _isin_kernel|s risinc s|tt|}ttj}tfddt|jD}tt|||||dd|Dt|d }|j|d}|rx|}|S)Nc3s|]}|jVqdS)N)r,)r4ra)rr1r2rAszisin..cSsi|]}dd|qS)cSsdS)Nrr1)rr1r1r2r^sz!isin...r1)r4rBr1r1r2rszisin..)Z adjust_chunksr-r)rB)r"rDrMr,r&rboolany)rrrinvertZ element_axesZ test_axesrrur1)rr2rs$  cCs:|}|dkr2t|}t|ts&td|f}d}nPy t|Wntk rX|f}YnXy t|Wntk r|f}YnXt|t|krtdxt||D]\}}| }||j|;}|jt dg}|jt dg}t |d||<t d|||<t |}t |}t ||||g|d}qW| |j}|S)NzAExpect `shift` to be an instance of Integral when `axis` is None.)rz,Must have the same number of shifts as axes.)rB) rxrnr rr;rTror`r,rVrDr$r)r3shiftrBrurarcZsl1Zsl2r1r1r2rolls<      rcCs |dS)N)rQ)r)r3r1r1r2rxsrxcsdkr"tddtjDntts2ftfddDrPtdtjtfddtjD}|S)Ncss|]\}}|dkr|VqdS)rNr1)r4rarRr1r1r2rAszsqueeze..c3s|]}j|dkVqdS)rN)r`)r4ra)r5r1r2rAsz,cannot squeeze axis with size other than onec3s&|]\}}|krdntdVqdS)rN)rV)r4rarc)rBr1r2rAs)rDrir`rnrrTrr,)r5rBrZr1)r5rBr2rs  rcsttt|}jdkr(tddkr<|}dt|j|tfddt|jD}|tfddt|jD}|S)Nrz!Condition must be one dimensionalrc3s,|]$}|krtdtntdVqdS)N)rVr;)r4ra)rB conditionr1r2rAszcompress..c3s"|]}|krntdVqdS)N)rV)r4ra)rBrr1r2rAs) r"r.rr,rTrxrrDrM)rr5rBr1)rBrr2compresss  rcCs(t|t}t|}t||S)N)r"r.rrrx)rrr1r1r2extractsrcCsJt||j}t|tjr.t|tr.t|||S|tdf||fSdS)N)rr,rnr7ndarrayr_take_dask_array_from_numpyrV)r5rrBr1r1r2takes  rcs<ttjstt|tst|jfdd|jjdS)Ncst|S)N)r7r)block)r5rBr1r2r^sz-_take_dask_array_from_numpy..)rr-)rnr7rrrrrr-)r5rrBr1)r5rBr2r s rcCstttj|d||jdS)N)r)r-)rrr7aroundr-)r/rr1r1r2rsrcCsddl}t||S)Nr)pandasr7r"isnull)rpdr1r1r2_asarray_isnullsrcCsddl}tt|ddS)z pandas.isnull for dask arrays rNr)r-)rr r)rrr1r1r2rsrcCs t|S)z pandas.notnull for dask arrays )r)rr1r1r2notnull%srh㈵>:0yE>cCs"ttj|||d}t|||ddS)N)rtolatol equal_nanr)r-)rr7iscloser )arr1arr2rrrrr1r1r2r*srcCst|||||dS)N)rrr)rrG)rrrrrr1r1r2allclose0srcGs t||S)N)r7choose)r5choicesr1r1r2variadic_choose5srcCstt|f|S)N)r r)r5rr1r1r2r9srcCstt|S)N)rr7 count_nonzero)rr1r1r2_isnonzero_vec>sr)ZotypescCsJytjt|jdtWntk r:|jttdSX|tSdS)N)r-) r7rhrDr-r.rrTrr)r5r1r1r2 isnonzeroEs rcsft|}t|}t|jtj|jdjdkrTt fddt t Dddt |ddS)N)r-rrcsg|]}|qSr1)rx)r4ra)rOr1r2r6\szargwhere..)rBr) r"rrrr`r7rrr,r%rMr;r)r5Znzr1)rOr2argwhereTs  "rcCs|dk|dkkrtd|dkr0|dkr0t|St|r~t||}t|}t|}t|j|j}|rj|n|}t|| |St tj |||SdS)Nz1either both or neither of x and y should be given) rTnonzeror7rr8r"r'r`r)r.r where)rr/rr-r`rPr1r1r2rbs   rcCstt|tjj|dS)N)rB)rr"r.r7rrq)r5rBr1r1r2rvsrcCstt|dddfS)Nr)rr"rx)r5r1r1r2 flatnonzero{srcs<t|jdkr2tfddtjdDSfSdS)Nrc3s|]}dd|fVqdS)Nr1)r4ra)rOr1r2rAsznonzero..)rr,rDrMr`)r5r1)rOr2rs  rcOs&tj|t||df|d|dS)Nfunclist func_argsfunc_kw)r7 piecewiserL)r/condlistrr1r1r2_int_piecewisesr cCsttj|f|S)N)r7r% unravel_index)r func_kwargsr1r1r2_unravel_index_kernelsrCc CsP|r:|jr:t|jttjt|ff|jd||dd}ntdd|D}|S)Nr)dimsorder)r-rrrcss|]}tdtjddVqdS))rr)r-rN)rr7r)r4rar1r1r2rAsz unravel_index..)sizerDrrr7rr;r)rrrZunraveled_indicesr1r1r2r s r cOs tt|f||jd|||dS)Nr )r-rrrr )rr r-)r/r rr9kwr1r1r2r s r c ss*tfddDs*d}t|dtkrDttjdtfddt  D}t fddt j D}tjd jjd j}tj|gd }t|||d S) Nc3s.|]&\}}j|D]}||dkVqqdS)rN)r)r4raZdivbd)r/r1r2rAszcoarsen..z6Coarsening factor does not align with block dimensionsZdaskzcoarsen-cs,i|]$}tj|ff|ddqS)rN)rcoarsen)r4r_)rSr reduction trim_excessr1r2rszcoarsen..c3s*|]"\}tfdd|DVqdS)c3s"|]}t|dVqdS)rN)rr)r4r)rSrar1r2rAsz$coarsen...N)rD)r4Zbds)rS)rar2rAs)r)r-)r)rGitemsrTinspectZgetfilegetattrr7__name__rrrrDrirrr,r-rrr) rr/rSrmsgrrrrrr1)rSrrrr/r2rs  rcsLtdg|dgg}ddtd|D}tdf|fdd|D}|S)z Split an array into a list of arrays (using slices) at the given breaks >>> split_at_breaks(np.arange(6), [3, 5]) [array([0, 1, 2]), array([3, 4]), array([5])] NcSsg|]\}}t||qSr1)rV)r4rar|r1r1r2r6sz#split_at_breaks..r csg|]}|fqSr1r1)r4rc)r3preslicer1r2r6s)r r rV)r3ZbreaksrBZ padded_breaksZslicesZ split_arrayr1)r3rr2split_at_breakss rcst|jttr,tj|jtjdk}|rNt t dk|jt dk rt dt|t}t|dddkrdt|}t|dt|}t|f|i|d|d}tfdd t|jD}t||}n|r |tdfd }tfd d tt|j|jD} || }tdd } t| | dk} t|| } tt|| g} d d| D} t| dS)Nrz9da.insert only implemented for monotonic ``obj`` argumentr,zvalues-r-r1)rr-c3s&|]\}}|krtn|VqdS)N)r;)r4rrc)rBobjr1r2rAszinsert..)Nc3s&|]\}\}}|kr|n|VqdS)Nr1)r4rZarr_bdZ values_bd)rBr1r2rAsrQcSsg|]}|jr|qSr1)nbytes)r4rar1r1r2r6szinsert..)rB)rr,rnrVr7rrr`r"r@rrrrrrrrrfrrDrir)rorrrrrLr r$)rr rrBZ scalar_objZ split_arrrr-Z values_shapeZ values_chunksrZ values_breaksZ split_valuesZ interleavedr1)rBr r2rks>          rkc CsVt|}|dkr0||}|j|j|j}nt|}t|jjtjtjfr`t |j|jd}nt |j|j}|j |j kr|dkrt d|j dkrt d|j d|j |krt dt||j dd|j }|d|}|rdd lm} || |}|j||d }t|||d ||}|rN|j |j krFt||j }||fS|SdS) NZf8z;Axis must be specified when shapes of a and weights differ.rz81D weights expected when shapes of a and weights differ.rz5Length of weights not compatible with specified axis.)rrQ) getmaskarray)rBr-)r-)r#rr-rfrrr7ZintegerZbool_r8r`rr,rTr)rNZmar"rqrcopy) r5rBrreturned is_maskedZavgZsclZwgtZ result_dtyper"r1r1r2_averages<     r&cCst||||ddS)NF)r%)r&)r5rBrr$r1r1r2average3sr')NN)F)F)F)N)r )N)rrQ)NN)NN)F)NNFNN)NrrN)Nr)r)F)FFF)F)FF)N)N)N)r)r)rrF)rrF)NN)N)r)F)r)NNFF)NNF)Z __future__rrrrrrZdistutils.versionr functoolsrrZnumbersrr Znumpyr7Ztoolzr r r Z compatibilityrZcorerbaserZhighlevelgraphrZutilsrrZcreationrrrrrrZwraprZufuncrrrrr r!r"r#r$r%r&r'r(r)r*Z einsumfuncsr+r3r8r>r?r@rFrHrIrNrUr[r\r]ZalphabetupperZALPHABETrmrlrvryr{r}rr~rrrrrrrrrrrrrrrrrrrxrrrrrrrrrrrrrrZ vectorizerrrrrrrr rr r rrrkr&r'r1r1r1r2s        <       #'  *   N   c 0   = f   /              / 0