o
    i                 	   @   s  d Z ddlZddlZddlZddlmZ ddlmZ ddlZ	ddl
mZmZ ddlmZmZmZmZmZ ddlmZmZmZmZmZmZmZmZ ddlmZmZmZ dd	lm Z m!Z!m"Z"m#Z#m$Z$m%Z%m&Z&m'Z'm(Z(m)Z)m*Z* dd
l+m,Z, ddl-m.Z. ddl/m0Z0m1Z1m2Z2m3Z3m4Z4 ddl5m6Z6m7Z7 ddl8m9Z9 ddl:m;Z;m<Z< ddl/m=Z= ddl>m?Z@mAZBmCZCmDZDmEZE dd ZFdd ZGdd ZHdd ZIdd ZJdd ZKdd ZLd d! ZMd"d# ZN	d@d$d%ZOd&d' ZPd(d) ZQe1eQd*d+ ZRe d,ejSd-d. ZTd/d0 ZUe d1ejVe&e*jWd2d3 ZX	d@d4d5ZYd6d7 ZZd8d9 Z[e ej\ejSej]e ej\ejSej^d:d; Z_e ej\ejSej`d<d= Zae ejbejSejcejcd>d? Zde eeejSd@dA Zfe dBejgdCdD Zhe dEejgejcdFdG ZiG dHdI dIejZkG dJdK dKekZlG dLdM dMekZmG dNdO dOekZnG dPdQ dQekZoG dRdS dSekZpG dTdU dUejZqdVdW Zre ej\ejSejgdXdY ZsdZd[ Ztd\d] Zud^d_ Zvd`da Zwdbdc Zxddde Zydfdg Zzdhdi Z{e3djdk Z|e3dldm Z}e0dndo Z~dpdq Ze1e	jdrds Ze0dtdu Ze1e	jdvdw Ze1e	jdxdy Zdzd{ Zd|d} Ze d~ejgdd Zdd Ze d~ejgej`dd Ze d~ejgeejcdd Ze1e	jd@ddZe!ejgddd Ze1e	jdAddZe1e	jdAddZe1e	jdBddZdd Zdd Ze dejgej`dd Ze dejgeejcdd Ze1e	jdd Ze1e	jdd Ze1e	jd@ddZe dejgdd Ze e	jejgdd Ze dejgdd Ze0dd Ze0dd Ze1e	jd@ddZe2ejgddCddZdd Ze1e	jdd Ze1e	jdd Ze1e	jdd Ze0dd Ze2ejgddd Ze3dd ZddĄ Ze1eddƍddȄ Ze dejgejdd˄ Ze!ejgd̃dd΄ Ze!ejgdσe!ejdσddф Ze!ejgd҃e!ejd҃ddԄ Ze!ejgdՃe!ejdՃddׄ Ze!ejgd؃ddڄ Ze!ejgdۃe!ejdۃdd݄ Ze!ejgdރe!ejdރdd Ze!ejddd Ze!ejddd Ze!ejddd Ze!ejddd Ze!ejgddd Ze!ejddd Ze$ejejǃe$ejejȃdd Zdd Ze!ejgddd Ze!ejde!ejddd Ze!ejddd Ze!ejgddd  Ze!ejgddd ZАdd Ze2ejgde2ejgddd	 ZҐd
d Ze!ejԐdeӃ e!ejԐdeӃ e dejejcdd Ze"ejgdd Ze dejgej.dd Ze"ejكdd Ze#ejكdd Ze dejej.dd Ze dejej݃dd Ze dejej.ejcdd Ze dejejejcd d! Ze%ejgd"d# Ze%ejكd$d% Ze%ejd&d' Ze ejejgejgd(d) Ze4ejgd*d+d, Zd-d. Zd/d0 Z		dCd1d2Z	d@d3d4Zd5d6 Zd7d8 Zd9d: Ze!ejgd;d<d= Ze d1eje&e*jWd>d? Ze ej\ejej]d@dA Ze ejbejej]ejcdBdC Ze eeejdDdE Ze e	jejgdFdG Ze d1eje&e*jWdHdI Ze eeej]e e	jeej]dJdK Ze eej`e e	jej`dLdM Ze d1eje&e*jWdNdO Ze e	jejcdPdQ Ze d1ej e&e*jWdRdS ZdTdU Ze=ejgdVdWdX ZdYdZ Ze3d[d\ Zd]d^ Zd_d` Zdadb Zdcdd Z	e3dedf Z
e1e	jefdgdhZe3didj Ze1e	jd@dkdlZe3dmdn Ze2ejgdodpdq Ze1e	jefdrdsZe1e	jd@dtduZe1e	jd@dvdwZe1e	jd@dxdyZe1e	jd@dzd{Ze1e	jd@d|d}Ze1e	jd@d~dZ dd Z!e1e!dd Z"e1e	j#ddefddZ$e1e	j%dDddZ&e1e	j'dDddZ(e1e	j)e2ejgddd Z*dd Z+e1e	j,dCddZ-e1e	j.dAddZ/dd Z0e3dd Z1e dejgdd Z2e1e	j3dd Z4dd Z5e3dd Z6e1e	j7dd Z8e1e	j9dd Z:e dejgeje dejgej.dd Z;e3dd Z<e1e	j=efddZ>e1ed@ddZ?e1ed@ddZ@dd ZAdd ZBdd ZCdd ZDdd ZEdd ZFdd ZGdd ZHe3dd ZIe1e	jJd@ddZKdÐdĄ ZLdŐdƄ ZMdǐdȄ ZNdɐdʄ ZOe3dːd̄ ZPe1e	jOd͐d΄ ZQdϐdЄ ZRdѐd҄ ZSe1e	jTdӐdԄ ZUe1e	jVdՐdք ZWe1e	jXdאd؄ ZYdِdڄ ZZdېd܄ Z[dݐdބ Z\dߐd Z]e3dd Z^e1e	j_dDddZ`dd Zae3dd Zbe1e	jcdd Zddd Zee3dd Zfe1e	jgdDddZhdd Zie3dd Zje1e	jkdd Zle3dd Zme1e	jndd Zoe3dd Zpe1e	jqdd Zre2ejgdd d Zse2ejgddd Zte1e	judd Zve1e	jwdd Zxe3d	d
 Zye1e	jzdd Z{e1e	j|dDddZ}e1e	j~dDddZe1e	jdd Ze1e	jdd Ze1e	jdd Zi Zdd ZdEddZe dejgdd Ze1e	jdd  Ze d!ejgej.e e	jejgej.d"d# Ze$ejgejgd$d% Ze$ejgeje$ejgeje$ejgej]e$ejgeje$ejgeje$ejgeje$ejgejd&d' Ze$ejgejd(d) Zd*d+ Zeed,d- Ze eejgej`ej`d.d/ Ze1e	jjjdFd0d1Ze1e	jjjd@d2d3Ze1ed4d5 Ze1e	jd6d7 Ze0d8d9 Ze1e	jd:d; Ze1e	jdGd>d?ZdS (H  z[
Implementation of operations on Array objects and objects supporting
the buffer protocol.
    N)ir)Constant)pndindexliteral_unroll)typestypingerrorscgutils	extending)as_dtype
from_dtypecarrayfarrayis_contiguous
is_fortrancheck_is_integertype_is_scalar)type_can_asarrayis_nonelikenumpy_version)lower_builtinlower_getattrlower_getattr_genericlower_setattr_generic
lower_castlower_constantiternext_implimpl_ret_borrowedimpl_ret_new_refimpl_ret_untrackedRefType)	signature)StringLiteral)register_jitableoverloadoverload_method	intrinsicoverload_attribute)	quicksort	mergesort)slicing)tuple_setitembuild_full_slice_tuple)overload_classmethod)parse_dtypeparse_shape_parse_nested_sequence_sequence_of_arrays_choose_concatenation_layoutc                 C   s4   t |j|t |j|g}| j|}|d| dS )zz
    Set the "range" metadata on a load instruction.
    Note the interval is in the form [lower_bound, upper_bound).
    rangeN)r   typemoduleZadd_metadataZset_metadata)builderloadlower_boundupper_boundZrange_operandsmd r;   `/var/www/html/eduruby.in/lip-sync/lip-sync-env/lib/python3.10/site-packages/numba/np/arrayobj.pyset_range_metadata+   s
   

r=   c                 C   s&   d|j jd > d }t| |d| dS )zF
    Mark the result of a load instruction as positive (or zero).
       r   N)r4   widthr=   )r6   r7   r9   r;   r;   r<   mark_positive6   s   r@   c                    s6   j tjG  fddd  S )z
    Return the Structure representation of the given *array_type*
    (an instance of types.ArrayCompatible).

    Note this does not call __array_wrap__ in case a new array structure
    is being created (rather than populated).
    c                       s4   e Zd Z fddZefddZ  ZS )zmake_array.<locals>.ArrayStructc                    s   t }z	| jd|}W n ty   t | | Y S w | jj }| |}|d u r9tj	| j
|dd}n|}|| j
|f}||fS )NZ	__array__T)zfill)r!   _contextget_functionNotImplementedErrorsuper
_make_refsdata_model_managerZ_get_be_typer	   alloca_once_builder)selfrefsigZ
array_impl	datamodelZbe_typeZ	outer_ref)ArrayStruct	__class__
array_typereal_array_typer;   r<   rF   L   s   


z*make_array.<locals>.ArrayStruct._make_refsc                    sp   | j }dkr | dS | d}g }tD ]}t||d|}||}|| t|| qt	||S )z\
            Override .shape to inform LLVM that its elements are all positive.
            r   shape)
rI   __getattr___get_ptr_by_namer3   r	   gep_inboundsr7   appendr@   
pack_array)rJ   r6   ptrdimsiZdimptrr7   )basendimr;   r<   rR   `   s   


z%make_array.<locals>.ArrayStruct.shape)__name__
__module____qualname__rF   propertyrR   __classcell__r;   rN   rP   r[   r\   rQ   )rO   r<   rN   J   s    rN   )Zas_arrayr	   create_struct_proxyr\   )rP   r;   rb   r<   
make_array>   s
   
-rd   c                 C   s   |  |j}| |S )zB
    Return the item size for the given array or buffer type.
    )get_data_typedtypeget_abi_sizeof)contextrP   lltyr;   r;   r<   get_itemsizez   s   
rj   c                 C   s"   |j rdnd}| j||j||dS )z3
    Load the item at the given array pointer.
    Nr>   align)alignedunpack_valuerf   )rh   r6   arraytyrX   rl   r;   r;   r<   	load_item   s   rp   c                 C   s$   |j rdnd}| j||j|||dS )z4
    Store the item at the given array pointer.
    Nr>   rk   )rm   
pack_valuerf   )rh   r6   ro   valrX   rl   r;   r;   r<   
store_item   s   rs   c                 C   s@   |j r| |||tj}t|||}|S | |||tj}|S )zM
    Fix the integer index' type and value for the given dimension size.
    )signedcastr   intpr*   Z	fix_indexuintp)rh   r6   idxtyidxsizeindr;   r;   r<   fix_integer_index   s   r|   c                 C   sZ   t |tjr)|jdkr)t |jtjsJ t|| ||}t| |||j}|j|fS ||fS )zS
    Normalize the index type and value.  0-d arrays are converted to scalars.
    r   )	
isinstancer   Arrayr\   rf   Integerrd   rp   data)rh   r6   rx   ry   idxaryZidxvalr;   r;   r<   normalize_index   s   
r   c                    s2   t |rt fddt||D  \}}||fS )z^
    Same as normalize_index(), but operating on sequences of
    index types and values.
    c                    s   g | ]\}}t  ||qS r;   )r   ).0rx   ry   r6   rh   r;   r<   
<listcomp>   s
    z%normalize_indices.<locals>.<listcomp>)lenzip)rh   r6   index_typesindicesr;   r   r<   normalize_indices   s
   r   c                 C   s  | j }| j}| j}	ttjdd}
|j|
 }t|j}t|	j}||@ |krC|| }d| j	 dt
|dkr7dnd d| d}t||d	u rRt||	d
d	}|tj}t|ttfrft|||}t|ttfrtt|||}t|tr}||}t|||||d}|d	u rt||	dd	|d< n||d< |tjd}t|||jj}|D ]}|j||dgd}q||d< t| }||krtd|| | D ]
\}}t| || q| S )z
    Helper function for populating array structures.
    This avoids forgetting to set fields.

    *shape* and *strides* can be Python tuples or LLVM arrays.
    r>   CzThe datamodel for type z is missing fields  .Nmeminfo)rR   stridesr   itemsizer   parentnswflagsnitemszmissing {0}) rB   rI   
_datamodelr   r~   float64rG   set_fieldsZ_fe_typer   
ValueErrorr   get_value_typeget_typerv   r}   tuplelistr	   rW   intdictget_constantunpack_tupler4   countmulkeysformatitemssetattr)arrayr   rR   r   r   r   r   rh   r6   rM   Zstandard_arrayZstandard_array_type_datamodelZrequired_fieldsZdatamodel_fieldsmissingmsgintp_tattrsr   unpacked_shapeaxlenZ
got_fieldskvr;   r;   r<   populate_array   sh   





r   c                 C   sh   |j }|j}|tjd}t||j| j}|D ]}|j	||dgd}q||_
|tjt|| |_dS )z
    Update some auxiliary information in *array* after some of its fields
    were changed.  `itemsize` and `nitems` are updated.
    r>   r   r   N)rB   rI   r   r   rv   r	   r   rR   r\   r   r   rj   r   )arytyr   rh   r6   r   r   r   r;   r;   r<   update_array_info   s   
r   c                 C   s   t  )z/Constrain axis values to valid positive values.rD   	func_namearg_namer\   axisr;   r;   r<   normalize_axis     r   c                    sL   t | ts
tdt |tstd| j d|j d  fdd}|S )Nz func_name must be a str literal.zarg_name must be a str literal.z: Argument z* out of bounds for dimensions of the arrayc                    s,   |dk r||7 }|dk s||krt  |S Nr   r   r   r   r;   r<   impl"  s
   z&normalize_axis_overloads.<locals>.impl)r}   r"   r   TypingErrorliteral_value)r   r   r\   r   r   r;   r   r<   normalize_axis_overloads  s   



r   Zgetiterc                 C   st   |j \}|\}| ||j}| tjd}t||}||_||_	| j
r,| j||| | }	t| ||j|	}
|
S r   )argsmake_helperreturn_typer   r   rv   r	   alloca_once_valueindexr   
enable_nrtnrtincref	_getvaluer   )rh   r6   rL   r   ro   r   iterobjzeroZindexptrresoutr;   r;   r<   getiter_array-  s   r   c              	   C   s   t j||j|jd}t j||j|jd}||d |}t ||j|}	|dd }
|dd }t|t	j
rEt| |||||	|
|}| S |
rIJ t| |||	S )z1 Evaluate `ary[idx]`, where idx is a single int. r   r   r>   N)r	   r   rR   r\   r   r   pointer_addr   r}   r   Buffer	make_viewr   rp   )rh   r6   r   r   aryry   shapesr   offsetdataptrview_shapesview_stridesretaryr;   r;   r<   _getitem_array_single_intE  s   r   Ziternextc                 C   s   |j \}|\}|j}| j|||d}t|| ||jd}	tj||	jdd\}
||j	}|
d||
}|| ||% t| ||j||	|}|| t||}|||j	 W d    d S 1 sew   Y  d S )Nvaluer>   r   <)r   rP   r   rd   r   r	   r   rR   r7   r   icmp_signed	set_validif_thenr   
yield_typeyield_increment_indexstore)rh   r6   rL   r   resultZitertyiterro   r   r   r   r   is_validr   Znindexr;   r;   r<   iternext_arrayZ  s"   

"r   c              	   C   s>  |  tjd}|  tjd}t||j|j}	t||j|j}
g }g }g }tdd |D }d}t	||D ]\}}|tj
u ri|jt| d | }t|D ]}|| ||	|  ||
|  |d7 }qNq6t|tjr| j|||d}t| ||| t|||	|  ||j t||}t|||
| }|| || n=t|tjrt| ||||	| }|rt| |||	| | || nt|r|| || |d8 }ntd|f |d7 }q6||jksJ ||jk r||	|  ||
|  |d7 }||jk stj| ||||ddd}|||fS )	z
    Perform basic indexing on the given array.
    A (data pointer, shapes, strides) tuple is returned describing
    the corresponding view.
    r   r>   c                 S      g | ]}t |r|qS r;   r   r   ry   r;   r;   r<   r         z"basic_indexing.<locals>.<listcomp>r   unexpected index type: %sF
wraparoundboundscheck)r   r   rv   r	   r   rR   r\   r   r   r   ellipsisr3   rV   r}   	SliceTyper   r*   guard_invalid_slice	fix_slicestartget_slice_lengthZ
fix_strider   r|   Zdo_boundscheckr   rD   get_item_pointer)rh   r6   r   r   r   r   r   r   oner   r   Zoutput_indicesZoutput_shapesZoutput_stridesnum_newaxesaxindexvalrx   	n_missingrZ   sliceshstr{   r   r;   r;   r<   basic_indexingv  sf   










r  c           	   	   C   s.   t || |}t|||||j|j|jd |S )zF
    Build a view over the given array with the given parameters.
    r   rR   r   r   r   r   )rd   r   r   r   r   )	rh   r6   r   r   r   r   r   r   r   r;   r;   r<   r     s   r   c              	   C   s^   t | |||||| jd\}}}	t|tjr$t| |||||||	}
|
 S |r(J t| |||S )zp
    Return the result of indexing *ary* with the given *indices*,
    returning either a scalar or a view.
    r   )r  enable_boundscheckr}   r   r   r   r   rp   )rh   r6   r   r   r   r   r   r   r   r   r   r;   r;   r<   _getitem_array_generic  s   r  c           	      C   sZ   |j \}}|\}}|jdksJ t|| ||}t| ||j|||f|f}t| ||j|S )z4
    Basic indexing with an integer or a slice.
    r>   )r   r\   rd   r  r   r   	rh   r6   rL   r   r   rx   r   ry   r   r;   r;   r<   getitem_arraynd_intp  s   

r
  c              	   C   s   |j \}}|\}}t|| ||}|j}tj||t|d}	t| |||	\}}	tdd |D r;t| |||||||	S t	| ||j
||||	}
t| ||j
|
S )z2
    Basic or advanced indexing with a tuple.
    r   c                 s       | ]	}t |tjV  qd S Nr}   r   r~   )r   tyr;   r;   r<   	<genexpr>      z&getitem_array_tuple.<locals>.<genexpr>)r   rd   r   r	   r   r   r   anyfancy_getitemr  r   r   )rh   r6   rL   r   r   tuptyr   tupr   r   r   r;   r;   r<   getitem_array_tuple  s    


r  c              	   C   s   |j \}}}|\}}}	t|tjr|j}
tj||t|d}n|f}
|f}t|| ||}t| ||
|\}
}zt	| ||||
|| j
d\}}}W n tyR   d}Y nw t|}|rbt| ||||
|S | ||	||j}	t| |||	| dS )zH
    array[a] = scalar_or_array
    array[a,..,b] = scalar_or_array
    r   r  TN)r   r}   r   	BaseTupler	   r   r   rd   r   r  r  rD   boolfancy_setsliceru   rf   rs   )rh   r6   rL   r   r   rx   valtyr   ry   rr   r   r   r   r   r   Zuse_fancy_indexingr;   r;   r<   setitem_array	  s4   

r  c           	      C   sD   |j \}|\}t|}|| ||}|j}||d}t| ||j|S r   )r   rd   rR   extract_valuer   r   )	rh   r6   rL   r   r   r   ZarysttyZshapearyr   r;   r;   r<   	array_len2  s   r  z
array.itemc                 C   s   |j \}|\}t|| ||}|j}|j|d||ddd d}| j|t|f W d    n1 s7w   Y  t	| |||j
S )N!=r>   FZlikelyz>item(): can only convert an array of size 1 to a Python scalar)r   rd   r   r   r   r4   	call_convreturn_user_excr   rp   r   )rh   r6   rL   r   r   r   r   r   r;   r;   r<   
array_item=  s   r!  zarray.itemsetc           
      C   s   |j \}}|\}}||jksJ t|| ||}|j}|j|d||ddd d}	| j|t	|	f W d    n1 s@w   Y  t
| ||||j |  S )Nr  r>   Fr  z/itemset(): can only write to an array of size 1)r   rf   rd   r   r   r   r4   r  r   r   rs   r   get_dummy_value)
rh   r6   rL   r   r   r  r   rr   r   r   r;   r;   r<   array_itemsetL  s   
r#  c                   @   @   e Zd ZdZdd Zdd Zdd Zdd	 Zd
d Zdd Z	dS )Indexerzq
    Generic indexer interface, for generating indices over a fancy indexed
    array on a single dimension.
    c                 C      t )ze
        Prepare the indexer by initializing any required variables, basic
        blocks...
        r   rJ   r;   r;   r<   prepareg     zIndexer.preparec                 C   r&  )z=
        Return this dimension's size as an integer.
        r   r'  r;   r;   r<   get_sizen     zIndexer.get_sizec                 C   r&  )z;
        Return this dimension's shape as a tuple.
        r   r'  r;   r;   r<   	get_shapet  r+  zIndexer.get_shapec                 C   r&  )z}
        Return a half-open [lower, upper) range of indices this dimension
        is guaranteed not to step out of.
        r   r'  r;   r;   r<   get_index_boundsz  r)  zIndexer.get_index_boundsc                 C   r&  )aW  
        Start indexation loop.  Return a (index, count) tuple.
        *index* is an integer LLVM value representing the index over this
        dimension.
        *count* is either an integer LLVM value representing the current
        iteration count, or None if this dimension should be omitted from
        the indexation result.
        r   r'  r;   r;   r<   	loop_head  s   	zIndexer.loop_headc                 C   r&  )z)
        Finish indexation loop.
        r   r'  r;   r;   r<   	loop_tail  r+  zIndexer.loop_tailN)
r]   r^   r_   __doc__r(  r*  r,  r-  r.  r/  r;   r;   r;   r<   r%  a  s    r%  c                   @   H   e Zd ZdZdd Zdd Zdd Zdd	 Zd
d Zdd Z	dd Z
dS )EntireIndexerz:
    Compute indices along an entire array dimension.
    c                 C   s2   || _ || _|| _|| _|| _| j tj| _d S r  )	rh   r6   r   r   dimr   r   rv   ll_intp)rJ   rh   r6   r   r   r3  r;   r;   r<   __init__  s   zEntireIndexer.__init__c                 C   sB   | j }|| jj| j| _t|| j| _	|
 | _|
 | _d S r  )r6   r  r   rR   r3  rz   r	   rH   r4  r   append_basic_blockbb_startbb_endrJ   r6   r;   r;   r<   r(    s
   
zEntireIndexer.preparec                 C      | j S r  rz   r'  r;   r;   r<   r*       zEntireIndexer.get_sizec                 C      | j fS r  r;  r'  r;   r;   r<   r,       zEntireIndexer.get_shapec                 C      |  d| jfS r   r4  rz   r'  r;   r;   r<   r-       zEntireIndexer.get_index_boundsc                 C   s   | j }| j t| jd| j || j || j || j}|j	|
d|| jdd || j W d    ||fS 1 sBw   Y  ||fS Nr   >=Fr  )r6   r   r   r4  r   branchr7  position_at_endr7   r   r   rz   r8  )rJ   r6   	cur_indexr;   r;   r<   r.    s   
zEntireIndexer.loop_headc                 C   D   | j }t||| j}||| j || j || j	 d S r  )
r6   r	   r   r7   r   r   rD  r7  rE  r8  rJ   r6   
next_indexr;   r;   r<   r/    s
   zEntireIndexer.loop_tailNr]   r^   r_   r0  r5  r(  r*  r,  r-  r.  r/  r;   r;   r;   r<   r2    s    r2  c                   @   r1  )IntegerIndexerz0
    Compute indices from a single integer.
    c                 C   s&   || _ || _|| _| j tj| _d S r  )rh   r6   ry   r   r   rv   r4  )rJ   rh   r6   ry   r;   r;   r<   r5    s   zIntegerIndexer.__init__c                 C      d S r  r;   r'  r;   r;   r<   r(       zIntegerIndexer.preparec                 C   s   t | jdS Nr>   )r   r4  r'  r;   r;   r<   r*       zIntegerIndexer.get_sizec                 C      dS Nr;   r;   r'  r;   r;   r<   r,    rM  zIntegerIndexer.get_shapec                 C   s   | j | j| j |  fS r  )ry   r6   addr*  r'  r;   r;   r<   r-    s   zIntegerIndexer.get_index_boundsc                 C   s
   | j d fS r  )ry   r'  r;   r;   r<   r.       
zIntegerIndexer.loop_headc                 C   rL  r  r;   r'  r;   r;   r<   r/    rM  zIntegerIndexer.loop_tailNrJ  r;   r;   r;   r<   rK    s    rK  c                   @   r1  )IntegerArrayIndexerz;
    Compute indices from an array of integer indices.
    c                 C   s@   || _ || _|| _|| _|| _|jdksJ | j tj| _	d S rN  )
rh   r6   rx   r   rz   r\   r   r   rv   r4  )rJ   rh   r6   rx   r   rz   r;   r;   r<   r5    s   zIntegerArrayIndexer.__init__c                 C   sD   | j }t|| jjd | _t|| j| _|	 | _
|	 | _d S r   )r6   r	   r   r   rR   idx_sizerH   r4  	idx_indexr6  r7  r8  r9  r;   r;   r<   r(    s
   
zIntegerArrayIndexer.preparec                 C   r:  r  rU  r'  r;   r;   r<   r*    r<  zIntegerArrayIndexer.get_sizec                 C   r=  r  rW  r'  r;   r;   r<   r,    r>  zIntegerArrayIndexer.get_shapec                 C   r?  r   r@  r'  r;   r;   r<   r-    rA  z$IntegerArrayIndexer.get_index_boundsc                 C   s   | j }| j t| jd| j || j || j || j}|j	|
d|| jdd || j W d    n1 s?w   Y  t| j|| jj| j| j|}t| j|| jj|| j}||fS rB  )r6   r   r   r4  rV  rD  r7  rE  r7   r   r   rU  r8  r   rh   rx   rf   r   r|   rz   )rJ   r6   rF  r   r;   r;   r<   r.    s&   zIntegerArrayIndexer.loop_headc                 C   rG  r  )
r6   r	   r   r7   rV  r   rD  r7  rE  r8  rH  r;   r;   r<   r/    s   
zIntegerArrayIndexer.loop_tailNrJ  r;   r;   r;   r<   rT    s    	rT  c                   @   r1  )BooleanArrayIndexerz>
    Compute indices from an array of boolean predicates.
    c                 C   sH   || _ || _|| _|| _|jdksJ | j tj| _t	| jd| _
d S Nr>   r   )rh   r6   rx   r   r\   r   r   rv   r4  r   r   )rJ   rh   r6   rx   r   r;   r;   r<   r5  &  s   zBooleanArrayIndexer.__init__c                 C   s^   | j }t|| jjd | _t|| j| _t|| j| _	|
 | _|
 | _|
 | _d S r   )r6   r	   r   r   rR   rz   rH   r4  rV  r   r6  r7  bb_tailr8  r9  r;   r;   r<   r(  /  s   

zBooleanArrayIndexer.preparec                 C   s   | j }t|| j}t|| j-}||}t| j|| j	j
| j	| j|j}|||||j}||| W d    n1 sAw   Y  ||S r  )r6   r	   r   r   	for_rangerz   r7   r   rh   rx   rf   r   r   rR  Zzextr4   r   )rJ   r6   r   loopcpredr;   r;   r<   r*  8  s   

	zBooleanArrayIndexer.get_sizec                 C   
   |   fS r  r*  r'  r;   r;   r<   r,  G  rS  zBooleanArrayIndexer.get_shapec                 C   r?  r   r@  r'  r;   r;   r<   r-  J  s   z$BooleanArrayIndexer.get_index_boundsc                 C   s  | j }| j | j| j | j | j| j || j || j || j}|| j}|j	|
d|| jdd || j W d    n1 sKw   Y  t| j|| jj| j| j|}|	|| || j W d    n1 sww   Y  t||}||| j ||fS )NrC  Fr  )r6   r   r   rV  r   rD  r7  rE  r7   r   r   rz   r8  r   rh   rx   rf   r   not_rZ  r	   r   )rJ   r6   rF  	cur_countr^  
next_countr;   r;   r<   r.  O  s,   zBooleanArrayIndexer.loop_headc                 C   s\   | j }|| j || j t||| j}||| j || j	 || j
 d S r  )r6   rD  rZ  rE  r	   r   r7   rV  r   r7  r8  rH  r;   r;   r<   r/  g  s   
zBooleanArrayIndexer.loop_tailNrJ  r;   r;   r;   r<   rX  !  s    		rX  c                   @   r1  )SliceIndexerz(
    Compute indices along a slice.
    c                 C   sL   || _ || _|| _|| _|| _|| _|| _| j tj	| _
t| j
d| _d S r   )rh   r6   r   r   r3  rx   r  r   r   rv   r4  r   r   )rJ   rh   r6   r   r   r3  rx   r  r;   r;   r<   r5  w  s   zSliceIndexer.__init__c                 C   s   | j }|| jj| j| _t| j|| j	| j
 t|| j
| j t|| j
j| _t|| j| _t|| j| _| | _| | _d S r  )r6   r  r   rR   r3  dim_sizer*   r   rh   rx   r  r   r	   Z
is_neg_intstepis_step_negativerH   r4  r   r   r6  r7  r8  r9  r;   r;   r<   r(    s   
zSliceIndexer.preparec                 C   s   t | j| jS r  )r*   r   r6   r  r'  r;   r;   r<   r*       zSliceIndexer.get_sizec                 C   r_  r  r`  r'  r;   r;   r<   r,    rS  zSliceIndexer.get_shapec                 C   s   t | j| j\}}||fS r  )r*   Zget_slice_boundsr6   r  )rJ   lowerupperr;   r;   r<   r-    s   zSliceIndexer.get_index_boundsc              	   C   s   | j }| j | jj| j | j | j| j || j |	| j |
| j}|
| j}|| j|d|| jj|d|| jj}|j|dd || j W d    ||fS 1 s_w   Y  ||fS )Nz<=rC  Fr  )r6   r   r  r   r   r   r   rD  r7  rE  r7   selectrg  r   stopr   r8  )rJ   r6   rF  rb  is_finishedr;   r;   r<   r.    s*   
zSliceIndexer.loop_headc                 C   sp   | j }|j|| j| jjdgd}||| j t||| j	}||| j	 |
| j || j d S )Nr   r   )r6   rR  r7   r   r  rf  r   r	   r   r   rD  r7  rE  r8  )rJ   r6   rI  rc  r;   r;   r<   r/    s   zSliceIndexer.loop_tailNrJ  r;   r;   r;   r<   rd  r  s    rd  c                   @   r$  )FancyIndexerz4
    Perform fancy indexing on the given array.
    c              	   C   s|  || _ || _|| _t||j|j| _t||j|j| _| j 	t
j| _g | _g }tdd |D }d}	d}
t||D ]\}}|t
ju rm|jt| d | }t|D ]}t|||||	}|| |	d7 }	|
d7 }
qTq<t|t
jr||||}t|||||	||}|| nkt|t
jrt||||| j|	 }t|||}|| nNt|t
jrt||||}t|jt
jrt||||| j|	 }nt|jt
jrt ||||}nJ || nt!|r| j|
 |	d8 }	nt"d|f |	d7 }	|
d7 }
q<|	|jksJ |	|jf|	|jk r(t|||||	}|| |	d7 }	|	|jk st||jks9J t||jf|| _#d S )Nc                 S   r   r;   r   r   r;   r;   r<   r     r   z)FancyIndexer.__init__.<locals>.<listcomp>r   r>   r   )$rh   r6   r   r	   r   rR   r\   r   r   r   r   rv   r4  newaxesr   r   r   r3   r2  rV   r}   r   r   rd  r   r|   rK  r~   rd   rf   rT  BooleanrX  r   AssertionErrorindexers)rJ   rh   r6   r   r   r   r   rr  r   r   Znew_axr   rx   r   rZ   indexerr  r{   r   r;   r;   r<   r5    sv   






"
zFancyIndexer.__init__c                 C   s^   | j D ]}|  q| jtjd}dd | j D }| jD ]	}|||f qt|d| _	d S )Nr>   c                 S      g | ]}|  qS r;   )r,  r   rZ   r;   r;   r<   r         z(FancyIndexer.prepare.<locals>.<listcomp>r;   )
rr  r(  rh   r   r   rv   ro  insertsumindexers_shape)rJ   rZ   r   Z	res_shaper;   r;   r<   r(    s   


zFancyIndexer.preparec                 C   r:  )z?
        Get the resulting data shape as Python tuple.
        )ry  r'  r;   r;   r<   r,    s   zFancyIndexer.get_shapec              	   C   s  t || jjks
J | j}tj}| d}| d}|}|}t| j| j	|D ]C\}	}
}|
||d|
|}|	 \}}|||}|||||}|d||}||||||}||||||}q&|||}||||}||||}||fS )z
        Get a half-open [lower, upper) range of byte offsets spanned by
        the indexer with the given strides and itemsize.  The indexer is
        guaranteed to not go past those bounds.
        r   r>   ==r   )r   r   r\   r6   r	   	false_bitr4  r   rr  ry  or_icmp_unsignedr-  r   subr   rR  rk  )rJ   r   r   r6   is_emptyr   r   ri  rj  rs  rR   strideZlower_indexZupper_indexZlower_offsetZupper_offsetZis_downwardsr;   r;   r<   get_offset_bounds  sB   


zFancyIndexer.get_offset_boundsc                 C   s    t dd | jD  \}}||fS )Nc                 s   s    | ]}|  V  qd S r  )r.  ru  r;   r;   r<   r  >  s    z+FancyIndexer.begin_loops.<locals>.<genexpr>)r   rr  )rJ   r   countsr;   r;   r<   begin_loops=  s   zFancyIndexer.begin_loopsc                 C   s   t | jD ]}|  qd S r  )reversedrr  r/  )rJ   rZ   r;   r;   r<   	end_loopsA  s   
zFancyIndexer.end_loopsN)
r]   r^   r_   r0  r5  r(  r,  r  r  r  r;   r;   r;   r<   rn    s    B&rn  c                 C   s  t ||j}t ||j}	|j}
t| |||||}|  |j}| }t	| |||}|j}t 
|| tjd}| \}}t j| ||
||	|j|d| jd	}t| |||}||}|||g}t| |||| t ||}||| |  t| ||| S )Nr   Fr   )r	   r   rR   r   r   rn  r(  r   r,  _empty_nd_implr   r   r   rv   r  get_item_pointer2layoutr  rp   r7   geprs   r   r   r  r   r   )rh   r6   rL   r   r   r   r   r   r   r   r   rs  Zout_tyZ
out_shapesr   Zout_dataZout_idx_rX   rr   curZnext_idxr;   r;   r<   r  F  s6   

r  c           	   	   C   s   |j \}}|\}}t|| ||}|jdkr4t| |||\}}t| ||j|||f|f}t| ||j|S t| ||||||f|fS )z3
    Advanced or basic indexing with an array.
    r   )r   rd   r\   r   r  r   r   r  r	  r;   r;   r<   fancy_getitem_arrayq  s   



r  c              	   C   s   |j }|d}|d}|jdv r|}	|||j}
|	|
fS |}	|}
t|jD ]-}||| ||| |}|d||}|	||
|
||
}
|	||	|
|	|}	q(|
|
|}
|d|j|}|	|||
}
|	|||	}	|	|
fS )z
    Compute a half-open range [lower, upper) of byte offsets from the
    array's data pointer, that bound the in-memory extent of the array.

    This mimics offset_bounds_from_strides() from
    numpy/core/src/private/mem_overlap.c
    r   r>   CFrC  rz  )r   r4   r  r   r   r3   r\   r~  r   rk  rR  )rh   r6   arrtyarrr   r   r   r   r   ri  rj  rZ   Zmax_axis_offsetZ
is_upwardsr  r;   r;   r<   offset_bounds_from_strides  s2   



r  c                 C   s.   | ||j}|||}|||}||fS )z
    Given [lower, upper) byte offsets and a base data pointer,
    compute the memory pointer bounds as pointer-sized integers.
    )ptrtointr4   rR  )rh   r6   ri  rj  r   Zdata_ptr_as_intr   endr;   r;   r<   compute_memory_extents  s   r  c           	      C   s&   t | |||||\}}t| ||||S )zr
    Compute a half-open range [start, end) of pointer-sized integers
    which fully contain the array data.
    )r  r  )	rh   r6   r  r  r   r   r   ri  rj  r;   r;   r<   get_array_memory_extents  s   
r  c                 C   s$   | |d|||d||}|S )z[
    Whether two memory extents [a_start, a_end) and [b_start, b_end)
    may overlap.
    r   )and_r}  )rh   r6   Za_startZa_endZb_startZb_endZmay_overlapr;   r;   r<   extents_may_overlap  s
   r  c                    s@  j dt d  jddc  |j|j}j |}	|	j  	|	}	 
|	 tj}
t |
(}t 	
j|}t |	|} 
 || W d    n1 smw   Y  W d    n1 s|w   Y   	
fdd} fdd}||fS )Nr   Fr  c                    s   t  } jddT\}}|  t j  | dd| W d    n1 s1w   Y  |  t j 	
j| dd| W d    n1 sUw   Y  W d    n1 sdw   Y  t 
 |S )NFr  r   )r	   rH   if_elser   r  r7   r  rp   )source_indicessrc_ptrZif_copyZ	otherwise)r6   rh   	copy_datacopy_layoutcopy_shapescopy_stridesptrtysrc_data
src_shapessrc_stridessrctyuse_copyr;   r<   src_getitem  s4   
z&maybe_copy_source.<locals>.src_getitemc                     sX    j dd  }  | } j |  W d    d S 1 s%w   Y  d S )NFr  )r   r7   bitcastr   freer   )r6   rh   r  r  	voidptrtyr;   r<   src_cleanup  s
   
"z&maybe_copy_source.<locals>.src_cleanup)r4   r	   r   r   r   r   r   r   allocater  r   r   r   rv   	loop_nestr  r  r7   )rh   r6   r  r  srcr  r  r  	allocsizer   r   r   r  dest_ptrr  r  r;   )r6   rh   r  r  r  r  r  r  r  r  r  r  r  r<   maybe_copy_source  s8   

"r  c           
         s  |  tjd}|  tjdt|t|kr0t|t| }g| | }|g| | }||fS t|t|k rt|t| } fdd|d| D }t j|tj} j	 
|dd d}	| j t|	f W d   n1 sww   Y  ||d }||d }||fS )	ae  
    Preprocess dimension for broadcasting.
    Returns (shapes, strides) such that the ndim match *target_shape*.
    When expanding to higher ndim, the returning shapes and strides are
    prepended with ones and zeros, respectively.
    When truncating to lower ndim, the shapes are checked (in runtime).
    All extra dimension must have size of 1.
    r   r>   c                       g | ]	}  d |qS rz  )r}  r   r  r6   r   r;   r<   r   %      z(_bc_adjust_dimension.<locals>.<listcomp>NFr  z,cannot broadcast source array for assignment)r   r   rw   r   	functoolsreducer  r	   true_bitr   ra  r  r   r   )
rh   r6   r   r   target_shaper   Znd_diffZ
dim_is_oneacceptedr   r;   r  r<   _bc_adjust_dimension  s,   	

r  c           
         s   g }g }|  tjd|  tjd fddt||D } fdd|D } fddt||D }	 fddt|	||D } fddt|	|D }||fS )	z
    Broadcast shapes and strides to target_shape given that their ndim already
    matches.  For each location where the shape is 1 and does not match the
    dim for target, it is set to the value at the target and the stride is
    set to zero.
    r   r>   c                    s   g | ]\}}  d ||qS )r  r   )r   taroldr6   r;   r<   r   @      z,_bc_adjust_shape_strides.<locals>.<listcomp>c                    r  r  r  )r   r  r  r;   r<   r   B      c                    s   g | ]
\}}  ||qS r;   )r  )r   xyr  r;   r<   r   C      c                    s    g | ]\}}}  |||qS r;   rk  )r   pr  r  r  r;   r<   r   D  s    c                    s   g | ]\}}  ||qS r;   r  )r   r  r  )r6   r   r;   r<   r   F  r  )r   r   rw   r   )
rh   r6   r   r   r  Z	bc_shapesZ
bc_stridesmismatchZ
src_is_onepredsr;   )r6   r   r   r<   _bc_adjust_shape_strides4  s    


r  c           	   	   C   s   t ||j}t ||j}t| ||||\}}t| ||||\}}|jt|dd}t|| |}t	||j
t ||t |||j|j|jd ||fS )zX
    Broadcast the given array to the target_shape.
    Returns (array_type, array)
    Ar\   r  r  )r	   r   rR   r   r  r  copyr   rd   r   r   rW   r   r   r   )	rh   r6   arrtyper  r  r   r   Znew_arrtypeZnew_arrr;   r;   r<   _broadcast_to_shapeK  s&   


r  c                 C   s,   |j |jddd}|||}dd }||fS )Nr  T)r\   r  readonlyc           
   	   S   sz   |\}}|j d }t|| ||}t||}t| ||||\}}t|dt| |j	dd  |
 }	t| ||j|	S )Nr   r   )r   rd   r	   r   r  r   r   r   r   r   r   r   r   )
rh   r6   rL   r   r  Zshape_r  r  destr   r;   r;   r<   codegenj  s   
z$_numpy_broadcast_to.<locals>.codegen)r  r   )	typingctxr   rR   retrL   r  r;   r;   r<   _numpy_broadcast_toe  s   
r  c                 C   s$   |j dd}||}dd }||fS )NTr  c                 S   sH   |\}|j d }t|| ||}t|jj|_| }t| ||j|S r   )	r   rd   r	   Zget_null_valuer   r4   r   r   r   )rh   r6   rL   r   r  r  r  r   r;   r;   r<   r    s   
z#get_readonly_array.<locals>.codegen)r  )r  r  r  rL   r  r;   r;   r<   get_readonly_array|  s   	r  c           
      C   s   | j }t|}t|}||krtd|D ]
}|dk rtdqd}|| }||k rM|| }|| }	||	ks:|dkrC|d7 }|d7 }ntd||k s*d S d S )NzDinput operand has more dimensions than allowed by the axis remappingr   z4all elements of broadcast shape must be non-negativer>   z=operands could not be broadcast together with remapped shapes)rR   r   r   )
r   Z
dest_shapeZ	src_shapeZsrc_ndimZ	dest_ndimrz   Z	src_indexZ
dest_indexZsrc_dimZdest_dimr;   r;   r<   _can_broadcast  s&   
r  c                 C   s   t | } t| | t| |S r  )npasarrayr  r  r   rR   r;   r;   r<   _default_broadcast_to_impl  s   


r  c                 C   s   t | s	tdt|tjrdd }|S t|tjr+t|jtjs)d}t|tS t|tj	rL|j
dkrLtdd |D sJd| d	}t|tS t|tj	ru|j
dkrut| tjoa| jdk}t| sh|rnd
d }|S d}t|d| }t|)Nz-The first argument "array" must be array-likec                 S   s   t | |fS r  )r  broadcast_tor  r;   r;   r<   r        z numpy_broadcast_to.<locals>.implz7The second argument "shape" must be a tuple of integersr   c                 S      g | ]}t |tjqS r;   )r}   r   IntegerLiteral)r   typr;   r;   r<   r     r   z&numpy_broadcast_to.<locals>.<listcomp>"z," object cannot be interpreted as an integerc                 S   s   t | } t| S r  )r  r  r  r  r;   r;   r<   r     s   
z/Cannot broadcast a non-scalar to a scalar arrayz:The argument "shape" must be a tuple or an integer. Got %s)r   r   r   r}   r   r   UniTuplerf   r  Tupler   allr~   r\   r   )r   rR   r   r   Zis_scalar_arrayr;   r;   r<   numpy_broadcast_to  s4   




r  c                 C   sr   t t|D ]0}|t| | }|| }|dk rtd|dkr!q| | dkr,|| |< q| | |kr6tdqd S )Nr   #negative dimensions are not allowedr>   z=shape mismatch: objects cannot be broadcast to a single shape)r3   r   r   )rmrR   rZ   r   tmpr;   r;   r<   numpy_broadcast_shapes_list  s   
r  c                     s   t | D ]7\}}t|tj}t|tjot|jtj}t|tjo't|jdk}|s;|s;|s;d| d| }t	|qd | D ]}t|tjrNt
 d q@t|tjr[t
 t| q@ dkrddd S d   fdd	}|S )
Nr   z	Argument z* must be either an int or tuple[int]. Got r>   c                  W   rP  rQ  r;   r   r;   r;   r<   <lambda>      z+ol_numpy_broadcast_shapes.<locals>.<lambda>r>   c                     sz   dg  }}t | D ] }t|trt|dkrt| | qt|tr+t| |f qt|D ]
\}}t|||}q0|S rY  )r   r}   r   r   r  r   	enumerater+   )r   r  r  argry   elemr  Ztup_initr;   r<   r   
  s   

z'ol_numpy_broadcast_shapes.<locals>.impl)r  r}   r   r   r  rf   r  r   r   r   maxr  )r   ry   r  Zis_intZis_int_tupleZis_empty_tupler   r   r;   r  r<   ol_numpy_broadcast_shapes  s0   
r  c                     s   t | D ]\}}t|std| dqd }d }| D ]$}t|tjtjfr+|j}n|}|d u r4|}q||kr@td|  qd t | D ])\}}t|tj	rXt
 |j qGt|tjtjtjfrit
 d qGtd| d   fdd	}|S )
Nz
Argument "z" must be array-likezTMismatch of argument types. Numba cannot broadcast arrays with different types. Got r   r>   zUnhandled type r   c                     s|   dg  }t | D ]}t| t|j q	}t D ]
}t|||| }qg }t | D ]}|tt|| q-|S rN  )	r   r  r  r  rR   r3   r+   rV   r  )r   rR   r   r  rZ   Zoutsr  r;   r<   r   <  s   
z$numpy_broadcast_arrays.<locals>.impl)r  r   r   r   r}   r   r~   r  rf   ZArrayCompatibler  r\   Numberrp  )r   ry   r  Zunified_dtypedtr   r;   r  r<   numpy_broadcast_arrays  s6   r  c           '         s:  |j \}}|\}}t| |}t |j}	t |j}
|j}t ||||}|  t	t
jrۈj| }t t |\t j}t j}j}tj}t|t|ksnJ t||D ]\}} | d||}qs j|dd d}j t|f W d   n1 sw   Y  t |||\}}||
|j\}}t |||\}}t ||||}t ||||\}}nmt	t
jr<j| }t|dksJ tt t
j!}| f} d|d |} j|dd d}j t|f W d   n	1 s(w   Y   fdd	}d
d }nfdd	}dd }"t
j#d} |$ \}!}"t%|"}"|j&D ]	}#|"'|#|  q\dd |"D }$||$}%( |%|j}%tj) ||	|
|j*|!dj+d	}&t, ||%|& |-  |  . S )z
    Implement slice assignment for arrays.  This implementation works for
    basic as well as fancy indexing, since there's no functional difference
    between the two for indexed assignment.
    r  Fr  z0cannot assign slice from input of different sizeNr>   r   c                    s,   | \} tjttj}| |fS r  )rC   operatorgetitemr!   r   rv   )r  ry   getitem_implr6   rh   r  Z	src_dtyper  r;   r<   r    s   z#fancy_setslice.<locals>.src_getitemc                   S   rL  r  r;   r;   r;   r;   r<   r    rM  z#fancy_setslice.<locals>.src_cleanupc                    s    S r  r;   )r  )r  r;   r<   r    rM  c                   S   rL  r  r;   r;   r;   r;   r<   r    rM  c                 S   s   g | ]}|d ur|qS r  r;   )r   r]  r;   r;   r<   r     r   z"fancy_setslice.<locals>.<listcomp>r   )/r   rd   r	   r   rR   r   r   rn  r(  r}   r   r   rf   r,  r  r{  r   r   r|  r   r   r  r   r   r  r  r   r  r  r  SequencerC   r!   rv   r   rw   r  r   ro  rw  ru   r  r  r  rs   r  r"  )'rh   r6   rL   r   r   r   r   r  r   Zdest_shapesdest_strides	dest_datars  Zindex_shaper  r  r  Zshape_errorur   r   Z	src_startZsrc_endZ
dest_lowerZ
dest_upperZ
dest_startZdest_endr  r  r  len_implZseq_lenr   Zdest_indicesr  rZ   r  rr   r  r;   r  r<   r  S  s   







r  c                    s   |j d }|j dd  }|d }|dd  } fddt||D }t ||d j}tjtjt|d}	t	
|j||	}
||f}|
|fS )Nr   r>   c                    "   g | ]\}}  ||tjqS r;   ru   r   rv   )r   r  rr   r   r;   r<   r         z#vararg_to_tuple.<locals>.<listcomp>rf   r   )r   r   r	   rW   r4   r   r  rv   r   r   r!   r   )rh   r6   rL   r   r   Zdimtysr   rY   rR   shapetynew_signew_argsr;   r   r<   vararg_to_tuple  s   
r	  zarray.transposec                 C   s   t | ||jd |d S r   )array_Tr   rh   r6   rL   r   r;   r;   r<   array_transpose  s   r  c                 C   sn   t | t t| krtdt |}| D ]}||ks t||kr$tdq||  |d d < ||  |d d < d S )Nzrepeated axis in transposez2axis is out of bounds for array of given dimension)r   r   r   abs)r   rR   r   r3  r  r;   r;   r<   permute_arrays  s   r  c              
      s  |j d }t| |d }|j d |d }}|j|j}}	tj}
t|
|||j	|j
g} fdd|D }t||D ]
\}} || qBtj|	dddtj|
} fdd|D }t||D ]\}}t| ||
 tj|g|g|d d qo tttjd	d |D  t|j }t||j |d  |d
 |j|j|jd | }t |j|S )Nr   r>   c                       g | ]}t  qS r;   r	   rH   r   r  )r6   ll_ary_sizer;   r<   r     r   z)array_transpose_tuple.<locals>.<listcomp>r   rf   r\   r  c                    s   g | ]	}t  qS r;   rd   r  )r6   rh   	np_ary_tyr;   r<   r     r  r   rR   r   r   r   c                 S   rt  r;   r   r   ar;   r;   r<   r   .  rv     r  ) r   rd   r   rf   r   r   rv   r   	ArrayTyperR   r   r   r   r~   r   rg   r   r  
as_pointercompile_internalr  r   r!   voidr   r   r7   r   r   r   r   r   )rh   r6   rL   r   r   r   Zaxistyr   Znum_axisrf   r4  ZarysZll_arysr  dstZnp_itemsizeZnp_arysZnp_aryZll_aryr  r   r;   )r6   rh   r  r  r<   array_transpose_tuple  sP   
	r   c                 C       t | |||\}}t| |||S r  )r	  r   rh   r6   rL   r   r  r  r;   r;   r<   array_transpose_vararg=     r#  c                 C   s:   t | tjrtd|d u rddd}|S ddd}|S )Nz#np.transpose does not accept tuplesc                 S      |   S r  	transposer  axesr;   r;   r<   np_transpose_implI  r>  z*numpy_transpose.<locals>.np_transpose_implc                 S   
   |  |S r  r&  r(  r;   r;   r<   r*  L  rS  r  )r}   r   r  r   ZUnsupportedError)r  r)  r*  r;   r;   r<   numpy_transposeC  s   


r,  Tc           	      C   s   |j dkr|}nFt|| ||}t|| |}t||j|j }t||j|j }t||jt||d d d t||d d d |j	|j
|jd | }t| |||S )Nr>   r  )r\   rd   r	   r   rR   r   r   r   rW   r   r   r   r   r   )	rh   r6   r  r   r   r   r  r   r   r;   r;   r<   r
  R  s    
r
  2   c                 C   sT   t | tjstdt |tjstdt |ttjfs#tdddd}|S )Nz+The first argument "start" must be a numberz+The second argument "stop" must be a numberz+The third argument "num" must be an integerr/  c                 S   s   t | ||}t d|S )Ng      $@)r  linspacepower)r   rl  numr  r;   r;   r<   r   o  s   znumpy_logspace.<locals>.implr/  )r}   r   r  r   r   r   r   )r   rl  r2  r   r;   r;   r<   numpy_logspacef  s   



r4  c                    s   t | tjsd}t|t |tjsd}t|t |ttjfs)d}t|tdd | |fD rJtt	
t| t|d  d
 fdd	}|S d
d	d}|S )Nz%The argument "start" must be a numberz$The argument "stop" must be a number%The argument "num" must be an integerc                 s   r  r  r}   r   Complexr   r  r;   r;   r<   r    r  z"numpy_geomspace.<locals>.<genexpr>r/  c           	         s   | dks|dkrt d | }  |}| jdk|jdk@ }t| dkt|dk@ }d}|r8| j} |j}d}|rC|  } | }| }t| }t|}t|||}|dkrd| |d< |dkrd||d< || S )Nr   &Geometric sequence cannot include zeror.  r>   y              ?)r   realr  signimaglog10logspace)	r   rl  r2  Zboth_imaginaryboth_negativeout_signlogstartlogstopr   Zresult_dtyper;   r<   r     s.   

znumpy_geomspace.<locals>.implc                 S   s   | dks|dkrt dt| dkt|dk@ }d}|r'|  } | }| }t| }t|}t|||}|dkrH| |d< |dkrH||d< || S )Nr   r9  r.  r>   )r   r  r;  r=  r>  )r   rl  r2  r?  r@  rA  rB  r   r;   r;   r<   r     s    

r3  )r}   r   r  r   r   r   r   r  r   r  Zresult_typer   )r   rl  r2  r   r   r;   rC  r<   numpy_geomspacev  s"   


2
rD  r>   c                 C   sR   t |ttjfstdt | tjstd| jdk r"tdddd}|S )	Nz*The second argument "k" must be an integerz'The first argument "m" must be an arrayr  Input must be >= 2-d.r>   c                 S   sp   |d }|dkr| d d  S |dkrt t | ddS |dkr(t t | S |dkr6t t | ddS t)N   r   r>   r     )r  swapaxesfliplrflipudrq  )r  r   r;   r;   r<   r     s   znumpy_rot90.<locals>.implr  )	r}   r   r   r   r   r   r~   r\   NumbaValueError)r  r   r   r;   r;   r<   numpy_rot90  s   




rL  c                 C   s   |  tj}| }|  tj}	t|	||||||||	g}
t|j	|
d}||j
}t||ddd}t||ddd}||}t||dd}t||dd}|	d}|||||||||j|g}|S )a  
    Call into Numba_attempt_nocopy_reshape() for the given array type
    and instance, and the specified new shape.

    Return value is non-zero if successful, and the array pointed to
    by *newstrides* will be filled up with the computed results.
    Znumba_attempt_nocopy_reshaperR   r   r   )r   r   rv   r  Zintcr   FunctionTyper	   Zget_or_insert_functionr5   r\   rU   rT   callr   )rh   r6   r   r   newndnewshape
newstridesr4  Zll_intp_starZll_intcfntyfnndrR   r   Z
is_f_orderr   r;   r;   r<   _attempt_nocopy_reshape  s2   	

rU  c           	      C   s   d}d}t |D ]\}}|dk r|d7 }|}q||9 }q|dkr*| |kr(tdd S |dkrO|dkr9d}| dk}n
| | }| | dk}|sItd|||< d S td)Nr   r>   z)total size of new array must be unchangedzmultiple negative shape values)r  r   )	origsizerR   Znum_neg_valueZ
known_sizer   r   Zneg_axinferredokr;   r;   r<   normalize_reshape_value  s*   

rY  zarray.reshapec              	   C   s  |j d }|j}|j d }|d }| tj}t||j}	t|| ||d }
t	
||	}||| tj|jddd}t|| |}| tj| |}t|||| | tj|jg|g|d d |
j}| |tttjtj||| g |j}t	
||	}t| |||
|||}|d||d}|| d}| j |t!|f W d    n1 sw   Y  t|| |}t||
j"|#||#||
j$|
j%|
j&d | }t'| ||j|S )	Nr   r>   r   r  r  rz  zincompatible shape for arrayr  )(r   r   r   r   rv   r   r  r   rd   r	   rH   r   r~   rf   r   rg   r   r  r  r   r  rY  r   r!   r  rw   r   rU  r}  r4   r   r  r   rD   r   r7   r   r   r   r   )rh   r6   rL   r   r   rettyr  rR   r4  Zll_shaper   rP  Zshape_ary_ty	shape_aryZshape_itemsizerz   rO  rQ  rX  failr   r  r   r;   r;   r<   array_reshape  sb   


r]  c                 C   r!  r  )r	  r]  r"  r;   r;   r<   array_reshape_vararg[  r$  r^  c                 C      dd }|S )Nc                 S   r+  r  )reshape)r  rP  r;   r;   r<   np_reshape_implc  rS  z#np_reshape.<locals>.np_reshape_implr;   )r  rP  ra  r;   r;   r<   
np_reshapea  s   rb  c                 C   sV   t | sd}t|t|tjrt|jtjs%t|tjs%d}t|dd }|S )N#The argument "a" must be array-likezBThe argument "new_shape" must be an integer or a tuple of integersc           	      S   s   t | } t | } t|tr(d}t |D ]}||9 }|dk r&d}t|qn|dk r2d}t||}| jdkrBt || j	S | | j  }| }t
|d D ]	}t || f}qQ|d | }t ||S )Nr>   r   z0All elements of `new_shape` must be non-negative)r  r  ravelr}   r   r   rz   zerosZastyperf   r3   concatenater`  )	r  	new_shapenew_sizeZ
dim_lengthr   msg2repeatsr   rZ   r;   r;   r<   r   x  s.   



znumpy_resize.<locals>.impl)r   r   r   r}   r   r  rf   r   )r  rg  r   r   r;   r;   r<   numpy_resizeh  s   


rk  c                 C   s^   t | s	tdt |stdt|rddd}|S t|tjs(tdddd}|S )Nz+The first argument "arr" must be array-likez/The second argument "values" must be array-likec                 S   s.   t t | } t t |}t | |fS r  )r  rd  r  rf  r  valuesr   r;   r;   r<   r     s   znp_append.<locals>.implz,The third argument "axis" must be an integerc                 S   s   t j| |f|dS Nr   r  rf  rl  r;   r;   r<   r        r  )r   r   r   r   r}   r   r   )r  rm  r   r   r;   r;   r<   	np_append  s   




rr  zarray.ravelc                 C   sN   dd }dd }|j d jdkr|}n|}| ||||}t| ||j|}|S )Nc                 S   s   |  | jS )zNo copy version)r`  rz   r   r;   r;   r<   
imp_nocopy     zarray_ravel.<locals>.imp_nocopyc                 S   r%  )zCopy version)flattenrs  r;   r;   r<   imp_copy     zarray_ravel.<locals>.imp_copyr   r   )r   r  r  r   r   )rh   r6   rL   r   rt  rw  impr   r;   r;   r<   array_ravel  s   rz  c                 C   s   dd }|  ||||S )Nc                 S   r%  r  )rd  r  r;   r;   r<   np_ravel_impl  r>  znp_ravel.<locals>.np_ravel_impl)r  )rh   r6   rL   r   r|  r;   r;   r<   np_ravel  s   r}  zarray.flattenc                 C   s,   dd }|  ||||}t| ||j|}|S )Nc                 S   s   |   | jS r  )r  r`  rz   rs  r;   r;   r<   ry    rh  zarray_flatten.<locals>.imp)r  r   r   )rh   r6   rL   r   ry  r   r;   r;   r<   array_flatten  s   r~  c                 C   sl   |d u r	t | n|}t | ||\}}}t |jD ]}|| }	|| }
|| }tt|	|
|||< q|S r  )r  
empty_likebroadcast_arraysndindexrR   minr  )r  a_mina_maxr   r  a_ba_min_ba_max_br   val_aZ	val_a_minZ	val_a_maxr;   r;   r<   _np_clip_impl  s   r  c                 C   sH   t | jD ]}| | }|| }|rt||||< qt||||< q|S r  )r  r  rR   r  r  )r  bZuse_minr   r   r  Zval_br;   r;   r<   _np_clip_impl_none  s   r  c                 C   s^  t | s	tdt|tjst |stdt|tjs't |s'tdt|tjs8t|s8d}t||d u pAt|tj}|d u pKt|tj}|rW|rWddd}|S t|tj}t|tj}	|rn|	rnddd}
|
S |r|	s|r{dd	d
}|S ddd}|S |s|	r|rddd}|S ddd}|S |rddd}|S |rddd}|S ddd}|S )Nrc  z6The argument "a_min" must be a number or an array-likez6The argument "a_max" must be a number or an array-likez5The argument "out" must be an array if it is providedc                 S   s   t d)Nz&array_clip: must set either max or minr   r  r  r  r   r;   r;   r<   
np_clip_nn	  r>  znp_clip.<locals>.np_clip_nnc                 S   sH   |d u r	t | n|}t | jD ]}| | }tt|||||< q|S r  )r  r  r  rR   r  r  r  r  r  r   r  r   r  r;   r;   r<   
np_clip_ss	  s
   znp_clip.<locals>.np_clip_ssc                 S   sB   |d u r	t | n|}t | jD ]}| | }t||||< q|S r  )r  r  r  rR   r  r  r;   r;   r<   
np_clip_sn%	  
   znp_clip.<locals>.np_clip_snc                 S   s   t | |}t| |||S r  r  	full_liker  )r  r  r  r   Z
a_min_fullr;   r;   r<   
np_clip_sa2	     znp_clip.<locals>.np_clip_sac                 S   sB   |d u r	t | n|}t | jD ]}| | }t||||< q|S r  )r  r  r  rR   r  r  r;   r;   r<   
np_clip_ns=	  r  znp_clip.<locals>.np_clip_nsc                 S   s   t | |}t| |||S r  r  )r  r  r  r   Z
a_max_fullr;   r;   r<   
np_clip_asJ	  r  znp_clip.<locals>.np_clip_asc                 S   s4   |d u r	t | n|}t | |\}}t||d|S )NTr  r  r  r  )r  r  r  r   r  r  r  r;   r;   r<   
np_clip_naV	     znp_clip.<locals>.np_clip_nac                 S   s4   |d u r	t | n|}t | |\}}t||d|S )NFr  )r  r  r  r   r  r  r  r;   r;   r<   
np_clip_an^	  r  znp_clip.<locals>.np_clip_anc                 S      t | |||S r  )r  r  r;   r;   r<   
np_clip_aaf	  s   znp_clip.<locals>.np_clip_aar  )	r   r   r   r}   r   NoneTyper~   r   r  )r  r  r  r   r   Za_min_is_noneZa_max_is_noner  Za_min_is_scalarZa_max_is_scalarr  r  r  r  r  r  r  r  r;   r;   r<   np_clip  sV   












r  clipc                 C   s   ddd}|S )Nc                 S   s   t | |||S r  )r  r  r  r;   r;   r<   r   q	  rh  zarray_clip.<locals>.implNNNr;   )r  r  r  r   r   r;   r;   r<   
array_clipo	  s   
r  c              	      s  |j |j ksJ |j|jksJ t|j}td tdtdtj fdd}| tjt| |}| tjt| |}| tj|j }	t	||
ddd}
t	||
ddd}tjtjd	dd
}| |}t|| tj|j g}| | tj}| tj|}t||g}|| |}t||
|||dd || |}t|||||dd | }| }|	||||| tj|g}ttjtj||tjtjtj}| ||||}t|| t| ||j|}|S )a  
    Attempt to fix up *ary* for switching from *oldty* to *newty*.

    See Numpy's array_descr_set()
    (np/core/src/multiarray/getset.c).
    Attempt to fix the array's shape and strides for a new dtype.
    False is returned on failure, True on success.
    r  r   Fc           	         s   | kr|d |kr}n
|d |kr}||kr'| ks%| dkr'dS |kr0| d }nd}||k rP|| dkr>dS || }||  |9  < |||< dS ||krn|| | }|| dkrbdS || ||< |||< dS 	 dS )Nr.  r   Fr>   Tr;   )	rT  rY   r   old_itemsizenew_itemsizer  rZ   newdimZ
bytelengthZ
any_layoutZc_layoutZf_layoutint8r;   r<   ry  	  s4   

z_change_dtype.<locals>.imprR   r   r   r>   r  Nr  )r\   r  ordr   r  r   rv   rj   r	   rU   rT   r~   rd   rW   rg   re   r   r   r!   booleanr  r   r   r   )rh   r6   Zoldtynewtyr   Z
new_layoutry  r  r  rT  Z
shape_dataZstrides_dataZshape_strides_array_typearyclsZshape_constantZsizeof_intpZstrides_constantr[  Zstrides_aryrR   r   r   rL   r   r;   r  r<   _change_dtypev	  s   	
0



	
r  c                 C      t | s	tddd }|S )Nz+The argument to np.shape must be array-likec                 S   s   t | jS r  )r  r  rR   r{  r;   r;   r<   r   	  rO  znp_shape.<locals>.implr   r   r   r  r   r;   r;   r<   np_shape	  s   
r  c                 C      dd }|S )Nc                    sJ   t |   t d d } fddt dd  D }t || S )Nr>   c                    s    g | ]\}} | |kr|qS r;   r;   )r   rZ   r  r  r;   r<   r   
       z5np_unique.<locals>.np_unique_impl.<locals>.<listcomp>)r  sortrd  r   r  r   )arheadtailr;   r  r<   np_unique_impl
  s   z!np_unique.<locals>.np_unique_implr;   )r  r  r;   r;   r<   	np_unique 
  s   r  c                 C   sh   dd }t | tjtjtjtjtjfr2t |tjrtS t |tjtjfr-t |j	tjr-|S t
dd S )Nc           	      S   s   t j|t jd}|jd dkrt| |d S t |dk r!tdt | }| }|jd }|j|jkr9tdt jt 	||j
d}d}t|D ]}|| |||||  < ||| 7 }qJ|S )Nrf   r   r>   r  z(operands could not be broadcast together)r  r  Zint64rR   np_repeat_impl_repeats_scalerr  r   rd  emptyrx  rf   r3   )	r  rj  Zrepeats_arrayasaaraveln	to_returnposrZ   r;   r;   r<   !np_repeat_impl_repeats_array_like
  s$   

z4np_repeat.<locals>.np_repeat_impl_repeats_array_likezIThe repeats argument must be an integer or an array-like of integer dtype)r}   r   r~   Listr  r  rp  r   r  rf   r   r   )r  rj  r  r;   r;   r<   	np_repeat

  s    r  c                 C   s   |dk rt dt| }| }|jd }|dkr"tjd|jdS |dkr+t|S tj|| |jd}t|D ]}|| ||| |d | < q9|S )Nr   r  r  r>   )	r   r  r  rd  rR   r  rf   r  r3   )r  rj  r  r  r  r  rZ   r;   r;   r<   r  7
  s   


r  repeatc                 C   r_  )Nc                 S      t | |S r  )r  r  )r  rj  r;   r;   r<   array_repeat_implK
  rO  z'array_repeat.<locals>.array_repeat_implr;   )r  rj  r  r;   r;   r<   array_repeatI
     r  c                 C      t |}dd }||fS )z"Computes the itemsize of the dtypec                 S   s,   |  |jd j}|  |j}|| |S r   )re   r   rf   r   rg   )cgctxr6   rL   llargsri   Zllintpr;   r;   r<   r  V
  s   z%_intrin_get_itemsize.<locals>.codegen)r   rv   )tyctxrf   rL   r  r;   r;   r<   _intrin_get_itemsizeQ
  s   
r  c                 C   rL  r  r;   )r  rf   r;   r;   r<   _compatible_view]
  rM  r  Zgeneric)targetc                 C   r_  )zDDetermines if the array and dtype are compatible for forming a view.c                 S   s   t |}|| jkrf| jdkrd}t|	 | jd }| j| dk}| jdk}| j| | jk}|r:|r:|r:d}t||| jk rR|dksJ| j| dkrPd}	t|	d S | j| | j }
|
| dkrhd}t|d S d S )Nr   zOChanging the dtype of a 0d array is only supported if the itemsize is unchangedr>   zJTo change to a dtype of a different size, the last axis must be contiguouszZWhen changing to a smaller dtype, its size must be a divisor of the size of original dtypezuWhen changing to a larger dtype, its size must be a divisor of the total size in bytes of the last axis of the array.)r  r   r\   r   rR   rz   r   )r  rf   Z
dtype_sizeZmsg1r   p1p2Zp3ri  Zmsg3r  Zmsg4r;   r;   r<   r   g
  s0   




z ol_compatible_view.<locals>.implr;   r  rf   r   r;   r;   r<   ol_compatible_viewa
  s   #r  z
array.viewc                 C   sD  |j d }|j}t|| ||d }t|| |}t|jj}t|D ]}	t||	}
|	dkr:|jj	}|
|
||_q#t||	|
 q#tdkrc| j}|t}||g |j R i }| ||}||| t| ||||}|d|t|j	d}|| d}| j|t|f W d    n1 sw   Y  | }t| ||j|S )Nr   r   )r>      rz  z"new type not compatible with array)r   r   rd   r   r   r   sortedgetattrr   r4   r  r   r   Ztyping_contextZresolve_value_typer  Zget_call_typerC   r  r}  r   r   r  r   r   r   r   )rh   r6   rL   r   r   rZ  r   r  fieldsr   rr   r  r  rR  Z_compatible_view_sigr   rX  r\  r   r   r;   r;   r<   
array_view
  s2   



r  rf   c                 C   s   |   }t| |||S r  )r"  r   rh   r6   r  r   r   r;   r;   r<   array_dtype
  s   r  rR   c                 C   (   t |}|| ||}|j}t| |||S r  )rd   rR   r   rh   r6   r  r   ro   r   r   r;   r;   r<   array_shape
     r  r   c                 C   r  r  )rd   r   r   r  r;   r;   r<   array_strides
  r  r  r\   c                 C      |  tj|j}t| |||S r  )r   r   rv   r\   r   r  r;   r;   r<   
array_ndim
  s   r  rz   c                 C   r  r  )rd   r   r   r  r;   r;   r<   
array_size
  s   r  r   c                 C   r  r  )rd   r   r   r  r;   r;   r<   array_itemsize
  r  r  nbytesc                 C   s2   t |}|| ||}||j|j}t| |||S )z"
    nbytes = size * itemsize
    )rd   r   r   r   r   r  r;   r;   r<   array_nbytes
  s   r  
contiguousc                 C   r  r  )r   r   r  	is_contigr   r  r;   r;   r<   array_contiguous
     r  c_contiguousc                 C   r  r  )r   r   r  Zis_c_contigr   r  r;   r;   r<   array_c_contiguous
  r  r  f_contiguousc                 C   r  r  )r   r   r  Zis_f_contigr   r  r;   r;   r<   array_f_contiguous   r  r  r  c                 C   s    |  tj|j }t| |||S r  )r   r   r  mutabler   r  r;   r;   r<   array_readonly  r$  r  ctypesc           	      C   sP   t |}|| ||}t|}| ||}|j|_|j|_| }t| |||S r  )rd   r   ArrayCTypesr   r   r   r   r   )	rh   r6   r  r   ro   r   Zactctinfor   r;   r;   r<   array_ctypes  s   
r  r   c                 C   s8   | j |||d}|j}||| tj}t| |||S Nr   )r   r   r  r   r   rv   r   )rh   r6   r  r   r  r   r;   r;   r<   array_ctypes_data  s   r  c                 C   s6   | j |||d}|j}||| |}t| |||S r  )r   r   r  r   r   )rh   r6   fromtytotyrr   r  r   r;   r;   r<   array_ctypes_to_pointer$  s   r  c           
      C   sr   t ||||d}ttj|j}|||j}ttj	||tj}|j
|j|tj|g}||| ||}	|	S )a)  Helper to invoke the contiguous checker function on an array

    Args
    ----
    checker :
        ``numba.numpy_supports.is_contiguous``, or
        ``numba.numpy_supports.is_fortran``.
    context : target context
    builder : llvm ir builder
    aryty : numba type
    ary : llvm value
    r   )rd   r   r  rv   r\   rg   r   rf   r!   Zbool_rR   r   r   r  )
checkerrh   r6   r   r   Ztup_intpr   Z	check_sigZ
check_argsr  r;   r;   r<   _call_contiguous_check-  s   
r  r   c                 C   s>   |  |t|}||_| }| j||| t| |||S r  )r   r   
ArrayFlagsr   r   r   r   r   )rh   r6   r  r   flagsobjr   r;   r;   r<   array_flagsG  s
   r  c                 C   sZ   |j jdkr| j|||d}tt| ||j |j}n|j jdk}| tj|}t	| |||S )Nr   r   )
rP   r  r   r  r   r   r   r   r  r   )rh   r6   r  r   r  r   rr   r;   r;   r<   array_flags_c_contiguousP  s   r   c                 C   sr   |j jdkr| j|||d}tt| ||j |j}n|j j}|j jdkr'|dkn|dv }| tj	|}t
| |||S )Nr  r   r>   r  )rP   r  r   r  r   r   r\   r   r   r  r   )rh   r6   r  r   r  r   r  rr   r;   r;   r<   array_flags_f_contiguous^  s   r  r:  c                 C   sH   |j tjv rt| |||ddS |j tjv rt| |||S tdtj )Nr:  attrzunsupported .real for {})	rf   r   complex_domainarray_complex_attrnumber_domainr   rD   r   r4   )rh   r6   r  r   r;   r;   r<   array_real_parto  s
   r  r<  c                 C   s   |j tjv rt| |||ddS |j tjv rHt|jdd|}t| |||g\}}t| |||}t	
||j||j|jd t| ||j| S tdtj )Nr<  r  Tr  r   zunsupported .imag for {})rf   r   r  r  r  r!   r  _parse_empty_like_argsr  r	   memsetr   r   r   r   r   r   r   rD   r   r4   )rh   r6   r  r   rL   r  r   r   r;   r;   r<   array_imag_partz  s   
r
  c                 C   s   |dvs
|j tjvrtd|t|}|| ||}|j j}| | |}|j	
|}	| | }
||j|
}|dkrK||tddg}|j|dd}t|| |}t||	d}t||| t| ||| S )	a  
    Given a complex array, it's memory layout is:

        R C R C R C
        ^   ^   ^

    (`R` indicates a float for the real part;
     `C` indicates a float for the imaginary part;
     the `^` indicates the start of each element)

    To get the real part, we can simply change the dtype and itemsize to that
    of the underlying float type.  The new layout is:

        R x R x R x
        ^   ^   ^

    (`x` indicates unused)

    A load operation will use the dtype to determine the number of bytes to
    load.

    To get the imaginary part, we shift the pointer by 1 float offset and
    change the dtype and itemsize.  The new layout is:

        x C x C x C
          ^   ^   ^
    )r:  r<  zcannot get attribute `{}`r<      r>   r  rf   r  )r   r   )rf   r   r  rD   r   rd   Zunderlying_floatrg   re   r   r4   r   r  r  r   r  r   IntTyper  r   r	   Zcopy_structr   r   )rh   r6   r  r   r  ro   r   ZfltyZsizeof_fltyr   Z
llfltptrtyr   Zresulttyr   replr;   r;   r<   r    s    r  conj	conjugatec                 C   r  )Nc                 S   s
   t | S r  )r  r  r  r;   r;   r<   r     rS  zarray_conj.<locals>.implr;   r  r   r;   r;   r<   
array_conj     r  c                 C   r%  r  )r"  )rh   r6   ZdtypetyZdtypevalr;   r;   r<   
dtype_type  rx  r  r4   kindZstatic_getitemc                 C   s<   |j }t|tjr| |d}t| |||S d}t|)zThis handles the "static_getitem" when a Numba type is subscripted e.g:
    var = typed.List.empty_list(float64[::1, :])
    It only allows this on simple numerical types. Compound types, like
    records, are not supported.
    NzlUnreachable; the definition of __getitem__ on the numba.types.abstract.Type metaclass should prevent access.)r   r}   r   r~   r   r   r   ZLoweringError)rh   r6   rL   r   rZ  r   r   r;   r;   r<   static_getitem_number_clazz  s   
r  c              	      s  t |}| ||}|j}t|tjstd||f ||}||}	t|tjr9|j	|j|j
|j
 dd}
n|j	|dd}
t |
}| |} tj|	}tj||j||jjd}t|tjrt||j|j
}| fdd|jD 7 }t||j|j
}| fdd|jD 7 }  |j}n|j}|j}  |}t|||| tj||j|jd	 | }t ||
|S )
zq
    Generic getattr() implementation for record arrays: fetch the given
    record member, i.e. a subarray.
    zattribute %r of %s not definedr  r  r  )r   c                       g | ]	}  tj|qS r;   r   r   rv   ru  rh   r;   r<   r     r  z(array_record_getattr.<locals>.<listcomp>c                    r  r;   r  ru  r  r;   r<   r     r  r  )rd   rf   r}   r   RecordrD   typeofr   NestedArrayr  r\   r   rv   r	   r   r   r4   r   rR   r   rg   re   r   r   r   r   r   )rh   r6   r  r   r  ro   r   Zrectyperf   r   ZrestyZrarytyZraryZconstoffsetZ
newdataptrrR   r   Zdatasizer   r;   r  r<   array_record_getattr  sN   


r  c                 C   s0   |d }t |tstt| ||jd |d |S rY  )r}   strrD   r  r   )rh   r6   rL   r   r   r;   r;   r<   array_record_getitem*  s   
r   c              	      s    || ||}||}t|tjret|}| |}|j}	 fdd|jD }
 fdd|j	D }t
||| |	}t||t
||
t
|| tj|jddd | }t |||S t
||| |}|jrudnd} ||||}t |||S )zT
    Generic getattr() implementation for records: get the given record member.
    c                    r  r;   r  r   r   r  r;   r<   r   B  r  z"record_getattr.<locals>.<listcomp>c                    r  r;   r  r!  r  r;   r<   r   D  r  Nr  r>   )sentry_record_alignmentr   r  r}   r   r  rd   rf   rR   r   r	   get_record_memberre   r   rW   r   rv   rz   r   r   rm   rn   )rh   r6   r  r   r  r   elemtyr   r   rf   rP  rQ  newdatar   dptrrl   r;   r  r<   record_getattr3  sB   







	
r'  c              	   C   s   |j \}}|\}}| || ||}	||}
t|
tjrHt|| ||d d}|j	}t
|||	|jj}t|||| tj|
j dS t
|||	| |
}| ||||
}|jr`dnd}| j||
|||d dS )zT
    Generic setattr() implementation for records: set the given record member.
    r>   r   Nrk   )r   r"  r   r  r}   r   r  r	   rc   r   r#  r4   ZpointeeZmemcpyr   rv   r   re   ru   rm   rq   )rh   r6   rL   r   r  r  r  r  rr   r   r$  Z
val_structr  r  r&  rl   r;   r;   r<   record_setattr[  s,   





r(  c                 C   s4   |  |jd |d }|| ||jd |d |d S )3
    Record.__getitem__ redirects to getattr()
    r   r>   )get_getattrr   )rh   r6   rL   r   r   r;   r;   r<   record_static_getitem_stry  s   r+  c                 C   s^   |j d j}t|j d j}| |j|| }| |j d |}|| ||j d |d || S )r)  r>   r   )r   r   r   r  Zinsert_const_stringr5   r*  )rh   r6   rL   r   ry   r  Zll_fieldr   r;   r;   r<   record_static_getitem_int  s
   r,  Zstatic_setitemc                 C   sJ   |j \}}}|\}}}	t|j||}
| ||
}|dusJ ||||	fS )3
    Record.__setitem__ redirects to setattr()
    N)r   r!   r   get_setattr)rh   r6   rL   r   rectyr  r  recry   rr   getattr_sigr   r;   r;   r<   record_static_setitem_str  s   
r2  c                 C   s^   |j \}}}|\}}}	t|j||}
t|j d j}| || |
}|dus(J ||||	fS )r-  r   N)r   r!   r   r   r  r.  )rh   r6   rL   r   r/  r  r  r0  ry   rr   r1  r  r   r;   r;   r<   record_static_setitem_int  s   
r3  c                 C   s   |  |||S )zB
    Create a constant array (mechanism is target-dependent).
    )make_constant_array)rh   r6   r  pyvalr;   r;   r<   constant_array  s   r6  c                 C   s0   t t d|j}|t| }t||S )zG
    Create a record constant as a stack-allocated array of bytes.
       )r   r  r  r  	bytearrayZtostringr	   r   )rh   r6   r  r5  Zltyrr   r;   r;   r<   constant_record  s   r9  c                 C   s"   t jt|t jd}| |||S )zM
    Create a constant array from bytes (mechanism is target-dependent).
    r  )r  r   r8  Zuint8r4  )rh   r6   r  r5  bufr;   r;   r<   constant_bytes  s   r;  c                 C   s0   |j \}}||krtjS dd }| ||||S )Nc                 S   s(   | j |j ko| j|jko| jj|jjkS r  )rR   r   r  r   r  r  r;   r;   r<   array_is_impl  s
   
zarray_is.<locals>.array_is_impl)r   r	   r{  r  )rh   r6   rL   r   atyZbtyr=  r;   r;   r<   array_is  s
   
r?  __hash__c                 C   s   dd S )Nc                 S   rL  r  r;   r  r;   r;   r<   r    r  zol_array_hash.<locals>.<lambda>r;   r  r;   r;   r<   ol_array_hash  rx  rA  c                 C   
   t | dS )zq
    Return the Structure representation of the given *flatiterty* (an
    instance of types.NumpyFlatType).
    flat_make_flattening_iter_cls)
flatitertyr;   r;   r<   make_array_flat_cls     
rG  c                 C   rB  )zv
    Return the Structure representation of the given *nditerty* (an
    instance of types.NumpyNdEnumerateType).
    ndenumeraterD  nditertyr;   r;   r<   make_array_ndenumerate_cls  rH  rL  c              	   C   s  |  tjd}|d}	|d ur|tj| tt|D ]T}
t	|||
}t
|||}||
 }|d||}t|| ||| |d urO||
 ||	 W d    n1 s^w   Y  ||| |d urq||
 q|d ur}|tj| ||	 ||	 d S )Nr   Zend_incrementr   )r   r   rv   r6  r   r	   
false_byter  r3   rU   r   r7   r   	if_likelyrD  	true_byterE  )rh   r6   r\   rR   r   end_flagloop_continue
loop_breakr   bbendr3  idxptrry   r   	in_boundsr;   r;   r<   _increment_indices  s.   

rV  c                 C   s*   t ||j|j}t| ||j||| d S r  )r	   r   rR   r\   rV  )rh   r6   r  r  r   rP  rR   r;   r;   r<   _increment_indices_array  s   rW  c                    s   j jtjjrnd	G dd dt}G dd d| G fddd|G dd	 d	|G d
d d|G dd d|G  	f
dddt}|S )zq
    Return the Structure representation of the given *nditerty* (an
    instance of types.NumpyNdIterType).
    r>   c                   @   sF   e Zd ZdZdd Zdd Zejdd Zdd	 Z	d
d Z
dd ZdS )z$make_nditer_cls.<locals>.BaseSubIterzF
        Base class for sub-iterators of a nditer() instance.
        c                 S   s&   || _ || _|| _|| _|| | _d S r  )nditermember_name	start_dimend_dimr\   )rJ   rX  rY  rZ  r[  r;   r;   r<   r5  *  s
   z-make_nditer_cls.<locals>.BaseSubIter.__init__c                 S   s   t | j| j| d S r  )r   rX  rY  )rJ   rX   r;   r;   r<   set_member_ptr1     z3make_nditer_cls.<locals>.BaseSubIter.set_member_ptrc                 S   s   t | j| jS r  )r  rX  rY  r'  r;   r;   r<   
member_ptr4     z/make_nditer_cls.<locals>.BaseSubIter.member_ptrc                 S   rL  r  r;   rJ   rh   r6   r;   r;   r<   init_specific8  rM  z2make_nditer_cls.<locals>.BaseSubIter.init_specificc                 S   rL  r  r;   rJ   rh   r6   logical_dimr;   r;   r<   rQ  ;  rM  z2make_nditer_cls.<locals>.BaseSubIter.loop_continuec                 S   rL  r  r;   rb  r;   r;   r<   rR  >  rM  z/make_nditer_cls.<locals>.BaseSubIter.loop_breakN)r]   r^   r_   r0  r5  r\  r  cached_propertyr^  ra  rQ  rR  r;   r;   r;   r<   BaseSubIter%  s    
re  c                   @   s0   e Zd ZdZdd Zdd Zdd Zdd	 Zd
S )z$make_nditer_cls.<locals>.FlatSubIterz
        Sub-iterator walking a contiguous array in physical order, with
        support for broadcasting (the index is reset on the outer dimension).
        c                 S   s$   | tjd}| t|| d S r   )r   r   rv   r\  r	   r   )rJ   rh   r6   r   r;   r;   r<   ra  G  s   z2make_nditer_cls.<locals>.FlatSubIter.init_specificc                 S   s   | | j}||j|gS r  )r7   r^  r  r   )rJ   rh   r6   r   r  r  r   r;   r;   r<   compute_pointerK  s   z4make_nditer_cls.<locals>.FlatSubIter.compute_pointerc                 S   s<   || j d kr|| j}t||}||| j d S d S rN  )r\   r7   r^  r	   r   r   )rJ   rh   r6   rc  r   r;   r;   r<   rQ  O  s
   z2make_nditer_cls.<locals>.FlatSubIter.loop_continuec                 S   sd   |dkr| tjd}||| j d S || jd kr0|| j}t||}||| j d S d S Nr   r>   )	r   r   rv   r   r^  r\   r7   r	   r   )rJ   rh   r6   rc  r   r   r;   r;   r<   rR  V  s   z/make_nditer_cls.<locals>.FlatSubIter.loop_breakN)r]   r^   r_   r0  ra  rf  rQ  rR  r;   r;   r;   r<   FlatSubIterA  s    rh  c                       s$   e Zd ZdZ fddZdd ZdS )z+make_nditer_cls.<locals>.TrivialFlatSubIterzx
        Sub-iterator walking a contiguous array in physical order,
        *without* support for broadcasting.
        c                    s    j rJ d S r  )need_shaped_indexingr`  rJ  r;   r<   ra  g  r  z9make_nditer_cls.<locals>.TrivialFlatSubIter.init_specificc                 S   s&   t |dksJ t |||j|S rN  )r   r  r   rJ   rh   r6   r   r  r  r;   r;   r<   rf  j  s   z;make_nditer_cls.<locals>.TrivialFlatSubIter.compute_pointerN)r]   r^   r_   r0  ra  rf  r;   rJ  r;   r<   TrivialFlatSubItera  s    rk  c                   @      e Zd ZdZdd ZdS )z'make_nditer_cls.<locals>.IndexedSubIterzA
        Sub-iterator walking an array in logical order.
        c                 S   s(   t || jks	J tj|||||ddS NFr  )r   r\   r	   r   rj  r;   r;   r<   rf  s  s   z7make_nditer_cls.<locals>.IndexedSubIter.compute_pointerNr]   r^   r_   r0  rf  r;   r;   r;   r<   IndexedSubItern      ro  c                   @   rl  )z'make_nditer_cls.<locals>.ZeroDimSubIterz5
        Sub-iterator "walking" a 0-d array.
        c                 S   s   |j S r  r  rj  r;   r;   r<   rf  }  r<  z7make_nditer_cls.<locals>.ZeroDimSubIter.compute_pointerNrn  r;   r;   r;   r<   ZeroDimSubIterx  rp  rq  c                   @   rl  )z&make_nditer_cls.<locals>.ScalarSubIterz8
        Sub-iterator "walking" a scalar value.
        c                 S   s   |S r  r;   rj  r;   r;   r<   rf    rM  z6make_nditer_cls.<locals>.ScalarSubIter.compute_pointerNrn  r;   r;   r;   r<   ScalarSubIter  rp  rr  c                       s|   e Zd ZdZej fddZ	fddZfddZdd	 Z	d
d Z
	fddZdd Zdd ZdS )zmake_nditer_cls.<locals>.NdIterz
        .nditer() implementation.

        Note: 'F' layout means the shape is iterated in reverse logical order,
        so indices and shapes arrays have to be reversed as well.
        c                    sd   g }j r nd}tjD ]\}}|\}}}}d| }	|| }
||
| |	|| q|S )N)rC  ZindexedZ0dZscalarzindex%d)ri  r  rr  rV   )rJ   lZ	factoriesrZ   r~  r  rZ  r[  r  rY  factory)rh  ro  rr  rk  rq  rK  r;   r<   subiters  s   
z(make_nditer_cls.<locals>.NdIter.subitersc                    s2  | tjd}||t||| _t|D ]\}}t|tjs2d| }t	
||| }	t| ||	 q| ||||}ttj}
d}d}t|D ]\}}t|tjre|jkre|| j}|| j} nqJdkslJ |||
d}| tjd}dd }t||D ]'\}}t|tjr|jdkrttjttj|j|
}|||||j|f qt	||} dkr|ddd	 }|d
||}||t	jt	j}js|f}t|ksJ t	j||jd}tD ]}t	|||}| || q|| _!t	"|||j| _t	
||| _#| j$D ]	}|%|| qdS )zY
            Initialize the nditer() instance for the specific array inputs.
            r   scalar%dNr;   r>   c                 S   s>   t | }t|D ]}| | |t || |  krtdqd S )Nz2nditer(): operands could not be broadcast together)r   r3   r   )rR   
main_shaper  rZ   r;   r;   r<   check_shape  s   zBmake_nditer_cls.<locals>.NdIter.init_specific.<locals>.check_shaper  r.  rz  r;  )&r   r   rv   
make_tupler  arraysr  r}   r~   r	   r   r   _arrays_or_scalarsr  r\   rR   r   r   r!   noner  r   r   rk  rO  rM  ri  r   rH   r4   r3   rU   r   r   rW   	exhaustedru  ra  )rJ   rh   r6   arrtysrz  r   rZ   r  rY  ZslotZmain_shape_tyrw  Zmain_nitemsr  rx  r  rL   r   Zshape_is_emptyr}  r   r3  rT  subiter)r  r\   rK  nshapesr;   r<   ra    sl   




z-make_nditer_cls.<locals>.NdIter.init_specificc                    s>  | d}t||| j}t|| |d || W d   n1 s*w   Y   j}t	|| j}| 
||||}| j}|d | |||||}	dd |	D }	t|	dkrg||	d  n||| j|	 t	|| j}
t||t|
|
|| jt| j||t| j|| || || dS )	zJ
            Compute next iteration of the nditer() instance.
            r  FNTc                 S   rt  r;   r  r   r   r;   r;   r<   r     rv  zEmake_nditer_cls.<locals>.NdIter.iternext_specific.<locals>.<listcomp>r>   r   )r6  r	   as_bool_bitr7   r}  if_unlikelyr   rD  rz  r   r{  r   _make_viewsr   r   ry  r   rR   rV  r  partial_loop_continue_loop_breakrE  )rJ   rh   r6   r   rS  r}  r~  rz  r   viewsrR   rJ  r;   r<   iternext_specific  sD   





z1make_nditer_cls.<locals>.NdIter.iternext_specificc                 S   @   | j D ]}|j|  kr|jk rn q|||||j  qd S r  )ru  rZ  r[  rQ  rJ   rh   r6   r3  r~  r;   r;   r<   r    
   
z.make_nditer_cls.<locals>.NdIter._loop_continuec                 S   r  r  )ru  rZ  r[  rR  r  r;   r;   r<   r    r  z+make_nditer_cls.<locals>.NdIter._loop_breakc                    s   dg }j }| j}j}	t|	tjrt|	}	n|	g}	 fddtD t||D ]=\}
}|
\}}}}|j	|j
 }dkrJ|ddd }|D ]}|| du sVJ | | ||	| || || |||< qLq-tdd |D svJ |S )z:
            Compute the views to be yielded.
            Nc              	          g | ]}  t |qS r;   r7   r	   rU   ru  r6   r   r;   r<   r   /  s    z?make_nditer_cls.<locals>.NdIter._make_views.<locals>.<listcomp>r  r.  c                 s   s    | ]}|V  qd S r  r;   r  r;   r;   r<   r  <  s    z>make_nditer_cls.<locals>.NdIter._make_views.<locals>.<genexpr>)rr  ru  r   r}   r   r  r   r3   r   rZ  r[  
_make_viewr  )rJ   rh   r6   r   r~  rz  r  rr  ru  rettysr~  r  r  Zarray_indicesZsub_indicesrZ   )r  narraysrK  r  r  r<   r  #  s0   



z+make_nditer_cls.<locals>.NdIter._make_viewsc                 S   s   t |tjr|jdksJ ||||||}||||}	t||}
||ttj	dd}||ttj	dd}t
|	||||
dd |	S )zD
            Compute a 0d view for a given input array.
            r   r;   N)r   )r}   r   r~   r\   rf  rd   rj   ry  r  rv   r   )rJ   rh   r6   r   rZ  r  r  r  rX   viewr   rR   r   r;   r;   r<   r  ?  s   
z*make_nditer_cls.<locals>.NdIter._make_viewc           	      S   s`   g }t t||D ]$\}\}}t|tjr#||||||d q	|t| d|  q	|S )Nr   rv  )r  r   r}   r   r~   rV   rd   r  )	rJ   rh   r6   r~  rz  rs  rZ   r  r  r;   r;   r<   r{  S  s   
z2make_nditer_cls.<locals>.NdIter._arrays_or_scalarsN)r]   r^   r_   r0  r  rd  ru  ra  r  r  r  r  r  r{  r;   
rh  ro  rr  rk  rq  r  r  r\   rK  r  r;   r<   NdIter  s    O*r  )r\   r  r   rz  ri  objectr	   rc   )rK  re  r  r;   r  r<   make_nditer_cls  s   
 
, Yr  c                    s$   | j  G  fdddt| }|S )zr
    Return the Structure representation of the given *nditerty* (an
    instance of types.NumpyNdIndexType).
    c                       s(   e Zd ZdZ fddZ fddZdS )z%make_ndindex_cls.<locals>.NdIndexIterz,
        .ndindex() implementation.
        c              	      s   | tjd}tj||j| tj d}t|tj}t D ]7}t	|||}|
|| || }	|d|	|}
t||
 |
tj| W d    n1 sRw   Y  q || _|| _t|||j| _d S Nr   r;  rz  )r   r   rv   r	   rH   r4   r   rM  r3   rU   r   r}  r  rO  r   r}  rW   rR   )rJ   rh   r6   r   r   r   r}  r3  rT  re  dim_is_emptyr\   r;   r<   ra  o  s&   
z3make_ndindex_cls.<locals>.NdIndexIter.init_specificc           
         s   | tjd} d}t  j}t | |	d  
| W d    n1 s1w   Y   fddtD }|D ]}t | qD|t ||j |	d t j}	t| |	jj  
|  | d S )Nr   r  Fc              	      s"   g | ]}  t j|qS r;   )r7   r	   rU   r   r   r3  r6   rJ   r;   r<   r     s    

zKmake_ndindex_cls.<locals>.NdIndexIter.iternext_specific.<locals>.<listcomp>T)r   r   rv   r6  r	   r  r7   r}  r  r   rD  r3   r@   r   rW   r4   r   rR   rV  r   rE  )
rJ   rh   r6   r   r   rS  r}  r   r7   rR   r  r  r<   r    s,   





z7make_ndindex_cls.<locals>.NdIndexIter.iternext_specificN)r]   r^   r_   r0  ra  r  r;   r  r;   r<   NdIndexIterj  s    r  )r\   r	   rc   )rK  r  r;   r  r<   make_ndindex_clsc  s   7r  c                    sX    dv sJ | j }|jdkrG  fdddt| }|S G  fdddt| }|S )N)rC  rI  r   c                       s8   e Zd ZdZ fddZ fddZdd Zdd	 Zd
S )z6_make_flattening_iter_cls.<locals>.CContiguousFlatIterzZ
            .flat() / .ndenumerate() implementation for C-contiguous arrays.
            c           	         s   | tjd}t||| _|j| _ dkr?tj||j	| tj|j
d}t|j
D ]}t|||}||| q*|| _d S d S )Nr   rI  r;  )r   r   rv   r	   r   r   r   r  rH   r4   r\   r3   rU   r   r   )	rJ   rh   r6   r  r  r   r   r3  rT  r  r;   r<   ra    s   
zD_make_flattening_iter_cls.<locals>.CContiguousFlatIter.init_specificc                    s   |j }|j} | j} d||}	||	 t |	W  |j	|g}
t
| ||
}dkr8|| n(| j fddt|D }t |}|t ||g t| || t |} || j W d    d S 1 sxw   Y  d S )Nr   rC  c              	      r  r;   r  r  r  r;   r<   r         

z\_make_flattening_iter_cls.<locals>.CContiguousFlatIter.iternext_specific.<locals>.<listcomp>)r\   r   r7   r   r   r   r	   rN  r  r   rp   r   r   r3   rW   make_anonymous_structrW  r   r   )rJ   rh   r6   r  r  r   r\   r   r   r   rX   r   idxvalsidxtupler  r  r<   r    s4   
"zH_make_flattening_iter_cls.<locals>.CContiguousFlatIter.iternext_specificc                 S   s   | |j|g}t||||S r  )r  r   rp   rJ   rh   r6   r  r  r   rX   r;   r;   r<   r    s   z>_make_flattening_iter_cls.<locals>.CContiguousFlatIter.getitemc                 S   s$   | |j|g}t||||| d S r  )r  r   rs   rJ   rh   r6   r  r  r   r   rX   r;   r;   r<   setitem  s   z>_make_flattening_iter_cls.<locals>.CContiguousFlatIter.setitemN)r]   r^   r_   r0  ra  r  r  r  r;   r  r;   r<   CContiguousFlatIter  s    r  c                       s<   e Zd ZdZdd Z fddZdd Zdd	 Zd
d ZdS )z+_make_flattening_iter_cls.<locals>.FlatIterz
            Generic .flat() / .ndenumerate() implementation for
            non-contiguous arrays.
            It keeps track of pointers along each dimension in order to
            minimize computations.
            c              	   S   s  | tjd}|j}|j}t||j|}tj||j	| tj|jd}	tj||j	| tj|jd}
t
|tj}t|D ]D}t||	|}t||
|}||| ||| || }|d||}t|| |tj| W d    n1 s}w   Y  q>|	| _|
| _|| _d S r  )r   r   rv   r   r\   r	   r   rR   rH   r4   r   rM  r3   rU   r   r}  r  rO  r   pointersr}  )rJ   rh   r6   r  r  r   r   r\   r   r   r  r}  r3  rT  ptrptrre  r  r;   r;   r<   ra    s>   
z9_make_flattening_iter_cls.<locals>.FlatIter.init_specificc              	      sH  |j }t |j|}t |j|}| j| j}	|tj	d}
 
d}t  | j}t | |d  | W d    n1 sJw   Y  |d t |	|d } |}t| ||}dkrs|| n fddt|D }t |}|t ||g tt|D ]y}t |}t  |}|| }|| } d	||}t |C  || t |	|} |}t ||} || t|d |D ]}t |	|} || q | W d    n	1 sw   Y   |
| q tj| j  |  | d S )
Nr   r  FTr>   rC  c              	      r  r;   r  r  r  r;   r<   r   5  r  zQ_make_flattening_iter_cls.<locals>.FlatIter.iternext_specific.<locals>.<listcomp>r   )r\   r	   r   rR   r   r   r  r   r   rv   r6  r  r7   r}  r  r   rD  rU   rp   r   r3   rW   r  r  r   r   rN  r   r   rO  rE  )rJ   rh   r6   r  r  r   r\   r   r   r  r   rS  is_exhaustedZlast_ptrrX   r   r  r  r3  rT  ry   r   r  rU  r  Z	inner_dimr  r  r<   r    sn   





z=_make_flattening_iter_cls.<locals>.FlatIter.iternext_specificc              	   S   s   |j }tj||j|d}tj||j|d}g }	tt|D ]}
|	||||
  |	|||
 }q|	
  t|||j|||j|	}|S )Nr   )r\   r	   r   rR   r   r  r3   rV   ZuremZudivreverser  r   r  )rJ   rh   r6   r  r  r   r\   r   r   r   r3  rX   r;   r;   r<   _ptr_for_index_  s   z:_make_flattening_iter_cls.<locals>.FlatIter._ptr_for_indexc                 S   s    |  |||||}t||||S r  )r  rp   r  r;   r;   r<   r  p  s   z3_make_flattening_iter_cls.<locals>.FlatIter.getitemc                 S   s&   |  |||||}t||||| d S r  )r  rs   r  r;   r;   r<   r  t  s   z3_make_flattening_iter_cls.<locals>.FlatIter.setitemN)	r]   r^   r_   r0  ra  r  r  r  r  r;   r  r;   r<   FlatIter  s    !Er  )rP   r  r	   rc   )rF  r  rP   r  r  r;   r  r<   rE    s   
D rE  rC  c                 C   sh   t t|}|| |}||_| |}|| ||dd}|| ||| | }t| |t||S Nr   )rK   )	rG  r   NumpyFlatTyper   rd   rT   ra  r   r   )rh   r6   r  r  flatiterclsflatiterarrclsr   r;   r;   r<   make_array_flatiter{  s   

r  c                 C   Z   |j \}|\}t|}|| ||d}|j}| |}	|	| ||jd}
|| |||
| d S r  )r   rG  rP   rd   r   r  )rh   r6   rL   r   r   rF  r  r  r  r  r  r;   r;   r<   iternext_numpy_flatiter     
r  c                 C   sj   |j d }|\}}t|}|| ||d}|j}| |}	|	| ||jd}
|| |||
|}t| ||j|S Nr   r   )r   rG  rP   rd   r   r  r   r   )rh   r6   rL   r   rF  r  r   r  r  r  r  r   r;   r;   r<   iternext_numpy_getitem  s   

r  c                 C   sf   |j d }|\}}}t|}|| ||d}|j}	| |	}
|
| ||jd}|| ||	||| |  S r  )r   rG  rP   rd   r   r  r"  )rh   r6   rL   r   rF  r  r   r   r  r  r  r  r;   r;   r<   iternext_numpy_getitem_any  s   


r  c           	      C   sF   |j d }t|}|| ||d d}| |j}|| ||jd}|jS r  )r   rG  rd   rP   r   r   )	rh   r6   rL   r   rF  r  r  r  r  r;   r;   r<   iternext_numpy_getitem_flat  s   
r  c           
      C   sr   |j \}|\}tt|}|| |}||_| |}|| ||dd}|| ||| | }	t	| ||j
|	S r  )r   rL  r   NumpyNdEnumerateTyper   rd   rT   ra  r   r   r   )
rh   r6   rL   r   r  r  	nditerclsrX  r  r   r;   r;   r<   make_array_ndenumerate  s   

r  c                 C   r  r  )r   rL  rP   rd   r   r  )rh   r6   rL   r   r   rK  rX  r  r  r  r  r;   r;   r<   iternext_numpy_nditer  r  r  c                    s^    fddt |j|D }ttt|}| }| | | }t |j	|S )zndindex(*shape)c                    r  r;   r  )r   Zargtyr  r   r;   r<   r     r  z&make_array_ndindex.<locals>.<listcomp>)
r   r   r  r   NumpyNdIndexTyper   ra  r   r   r   )rh   r6   rL   r   rR   r  rX  r   r;   r   r<   make_array_ndindex  s   

r  c           
         s   |j j}|dkr%|jd j|d }t ||} fdd|D }ng }ttt	|}| }|
 | | }	t |j |	S )zndindex(shape)r   c                    s   g | ]}  |tjqS r;   r  r   r6   rh   rx   r;   r<   r         z,make_array_ndindex_tuple.<locals>.<listcomp>)r   r\   r   rf   r	   r   r  r   r  r   ra  r   r   )
rh   r6   rL   r   r\   r  rR   r  rX  r   r;   r  r<   make_array_ndindex_tuple  s   
r  c                 C   s6   |j \}|\}t|}|| ||d}|| || d S r  )r   r  r  )rh   r6   rL   r   r   rK  rX  r  r;   r;   r<   iternext_numpy_ndindex  s
   r  c           	      C   sn   |j }|j}t|jd tjrt||d }n|d g}t|| |}|	| ||| |
 }t| |||S )z
    nditer(...)
    r   )r   rz  r}   r   r   r  r	   r   r  ra  r   r   )	rh   r6   rL   r   rK  r~  rz  rX  r   r;   r;   r<   make_array_nditer  s   
r  c                 C   s2   |j \}|\}t|| ||d}|| || d S r  )r   r  r  )rh   r6   rL   r   r   rK  rX  r;   r;   r<   iternext_numpy_nditer2/  s   r  c                 C   sb  t |}|| |}| |j}| tjt| |}| tjd}tt	dd}	|D ]}
|
||
}||d}||	||d}	q*|jdkrKd}nN|jdkrp|g}t|dd D ]}|||d | q[tt|}n)|jdkr|g}|dd D ]}|||d | q~t|}ntd|j|
||}||d}||	||d}	|j|	d	d
 | j|td W d   n1 sw   Y  |j}| |}| tj|}|  ||f}ttj}tt|}t ||tjtj}| !|t"||}| j#$||}| %tj}t&j'|||d}t&j'|||d}t(||)||* ||||d |S )a	  Utility function used for allocating a new array during LLVM code
    generation (lowering).  Given a target context, builder, array
    type, and a tuple or list of lowered dimension sizes, returns a
    LLVM value pointing at a Numba runtime allocated array.
    r>   r   r;   r   Nr.  r  z3Don't know how to allocate array with layout '{0}'.Fr  )z[array is too big; `arr.size * arr.dtype.itemsize` is larger than the maximum possible size.)r  r  )+rd   re   rf   r   r   rv   rj   r   r   r  Zsmul_with_overflowr  r|  r\   r  r  rV   r   r   rD   r   r   r  r   r   Zget_preferred_array_alignmentZuint32r"  MemInfoPointervoidptrZTypeRefr4   r!   r  _call_allocatorr   Zmeminfo_datar   r	   rW   r   r  r  )rh   r6   r  r   r  r   datatyper   Zarrlenoverflowr   Zarrlen_multr   Zdimension_sizeZallocsize_multr  rf   Z	align_valrl   r   mipZarytypeclassargtypesr   r   r   Zshape_arrayZstrides_arrayr;   r;   r<   r  <  st   




	
r  	_allocatec                 C   s   dd }|S )zTImplements a Numba-only default target (cpu) classmethod on the array
    type.
    c                 S   s
   t ||S r  )intrin_alloc)clsr  rl   r;   r;   r<   r     rS  z _ol_array_allocate.<locals>.implr;   )r  r  rl   r   r;   r;   r<   _ol_array_allocate  s   r  c                 C   s   |  ||S )z9Trampoline to call the intrinsic used for allocation
    )r  )r  rz   rl   r;   r;   r<   r    s   r  c                 C   s(   dd }t t j}t|||}||fS )z3Intrinsic to call into the allocator for Array
    c                 S   s   |\}}| j |||}|S r  )r   Zmeminfo_alloc_aligned)rh   r6   r!   r   r  rl   r   r;   r;   r<   r    s   zintrin_alloc.<locals>.codegen)r   r  r  r!   )r  r  rl   r  r  rL   r;   r;   r<   r    s   r  c              	   C   s   dd }t |tjrd}| |||tjg}nt |tjsJ |j}tj|||d}g }|D ]}|	|| ||j
| q.| |tjd}	t|D ])}
|d||
 |	}t|| | j|td W d   n1 slw   Y  qH|S )	z;
    Parse the shape argument to an array constructor.
    c           
      S   s   |  tj}|j}t|}t|d|d > d }|j|k r&|||}|S |j|kr`|d||}	|j	|	dd | j
|td W d   n1 sLw   Y  |j|kr^|||}|S |}|S )z0Cast src to intp only if value can be maintainedr>   >Fr  )z#Cannot safely convert value to intpN)r   r   rv   r?   r   r  r   Zsextr   r   r  r   r   trunc)
rh   r6   Zsrc_tr  r   Z
intp_widthZintp_irmaxvalr   Z	is_largerr;   r;   r<   safecast_intp  s(   



z#_parse_shape.<locals>.safecast_intpr>   r   r   r   )znegative dimensions not allowedN)r}   r   r   ru   rv   r  r   r	   r   rV   r4   get_constant_genericr3   r   r  r  r   r   )rh   r6   r  rr   r  r\   Zpassed_shapesr   r   r   r3  Zis_negr;   r;   r<   _parse_shape  s(   r  c                 C   s*   |j d }|d }|j}|t| |||fS )zL
    Parse the arguments of a np.empty(), np.zeros() or np.ones() call.
    r   )r   r   r  )rh   r6   rL   r   ZarrshapetypeZarrshaper  r;   r;   r<   _parse_empty_args  s   
r  c                 C   sT   |j d }t|tjr%t|| ||d d}tj||j|jd}|j	|fS |j	dfS )z_
    Parse the arguments of a np.empty_like(), np.zeros_like() or
    np.ones_like() call.
    r   r   r   r;   )
r   r}   r   r~   rd   r	   r   rR   r\   r   )rh   r6   rL   r   arytyper   r   r;   r;   r<   r    s   


r  c                 C   s&   t |tjrd|  d}t|d S )NzIf np.z0 dtype is a string it must be a string constant.)r}   r   ZUnicodeTyper   r   )fnamerf   r   r;   r;   r<   _check_const_str_dtype  s   
r  c                 C   "   |j }||||}dd }||fS )Nc                 S   (   t | |||\}}t| |||}| S r  )r  r  r   r  r6   rL   r  r  r   r   r;   r;   r<   r       znumpy_empty_nd.<locals>.codegenZinstance_type)r  Zty_shapety_dtypety_retty_refty_rettyrL   r  r;   r;   r<   numpy_empty_nd     r  c                    s   t d| |tu st|tjr|jtu st|rtj}nt|}t	| }|d ur>|d ur>tj
||dd tf fdd	}|S d|  d| d}t|)	Nr  r   r  c                       t | | S r  )r  )rR   rf   rZ  r;   r<   r     rO  zol_np_empty.<locals>.implz.Cannot parse input types to function np.empty(, ))r  floatr}   r   Function
typing_keyr   doublety_parse_dtypety_parse_shaper~   r   r   )rR   rf   nb_dtyper\   r   r   r;   r  r<   ol_np_empty  s    



r  c                 C   r  )Nc                 S   r  r  )r  r  r   r  r;   r;   r<   r    r  z$numpy_empty_like_nd.<locals>.codegenr  )r  Zty_prototyper  r  r  rL   r  r;   r;   r<   numpy_empty_like_nd  r  r  c                    s   t d| t|st|}nt| tjr| j}n| }|d ur?t| tjr7| jdkr,| jnd}| j||dd nt|dd nd|  d| d	}t	
|d fd
d	}|S )Nr  r  r   F)rf   r  r  r   z3Cannot parse input types to function np.empty_like(r  r  c                    r  r  )r  )r  rf   r  r;   r<   r   9  rO  zol_np_empty_like.<locals>.implr  )r  r   r  r}   r   r~   rf   r  r  r   r   )r  rf   r  r  r   r   r;   r  r<   ol_np_empty_like%  s&   


r  c                 C   r  )Nc                 S   s<   t |jd | ||d }t||j||j|jd d S r   )rd   r   r	   r	  r   r   r   r   )r  r6   rL   r  r   r;   r;   r<   r  B  s   z(_zero_fill_array_method.<locals>.codegen)r   r|  )r  rJ   rL   r  r;   r;   r<   _zero_fill_array_method>  s   
r  
_zero_fillc                 C   r  )z?Adds a `._zero_fill` method to zero fill an array using memset.c                 S   s   t |  d S r  )r  r'  r;   r;   r<   r   L  rO  z ol_array_zero_fill.<locals>.implr;   )rJ   r   r;   r;   r<   ol_array_zero_fillI  r  r  c                 C   s   t d| tfdd}|S )Nre  c                 S      t j| |d}|  |S Nr  )r  r  r  )rR   rf   r  r;   r;   r<   r   U     zol_np_zeros.<locals>.impl)r  r  rR   rf   r   r;   r;   r<   ol_np_zerosQ  s   
r  c                 C      t d| ddd}|S )N
zeros_likec                 S   r  r  )r  r  r  )r  rf   r  r;   r;   r<   r   a  r  zol_np_zeros_like.<locals>.implr  r  r  r;   r;   r<   ol_np_zeros_like\  s   

r  c                 C   r  )N	ones_likec                 S   2   t j| |d}|j}tt|D ]}d||< q|S Nr  r>   r  r  rC  r3   r   )r  rf   r  arr_flatry   r;   r;   r<   r   m  
   
zol_np_ones_like.<locals>.implr  r  r  r;   r;   r<   ol_np_ones_likeh  s   

r  c                    s2   t d| t|st| n| d fdd	}|S )Nfullc                    s0   t |  }|j}tt|D ]}|||< q|S r  r  r  rC  r3   r   )rR   
fill_valuerf   r  r  ry   r  r;   r<   r  ~  
   
zimpl_np_full.<locals>.fullr  )r  r   r  )rR   r	  rf   r  r;   r
  r<   impl_np_fullv  s   

r  c                 C   s   t d| ddd}|S )Nr  c                 S   s0   t | |}|j}tt|D ]}|||< q|S r  r  )r  r	  rf   r  r  ry   r;   r;   r<   r    r  z$impl_np_full_like.<locals>.full_liker  r  )r  r	  rf   r  r;   r;   r<   impl_np_full_like  s   

r  c                 C   r  )Nonesc                 S   r  r  r  )rR   rf   r  r  ry   r;   r;   r<   r     r  zol_np_ones.<locals>.implr  r  r  r;   r;   r<   
ol_np_ones  s   

r  c                    s4   t d| t|st| ntj d fdd	}|S )Nidentityc                    s.   t | | f }t| D ]}d|||f< q|S rN  )r  re  r3   )r  rf   r  rZ   r
  r;   r<   r    s   z"impl_np_identity.<locals>.identityr  )r  r   r  r   r  )r  rf   r  r;   r
  r<   impl_np_identity  s   

r  c                 C   rL  r  r;   NMr;   r;   r<   _eye_none_handler  rM  r  c                 C   s$   t |tjrdd }|S dd }|S )Nc                 S      | S r  r;   r  r;   r;   r<   r     rM  z$_eye_none_handler_impl.<locals>.implc                 S      |S r  r;   r  r;   r;   r<   r     rM  r}   r   r  )r  r  r   r;   r;   r<   _eye_none_handler_impl  s
   r  c                    sf   |d u s
t |tjrtt nt |tjtjfr"tt	|d| nt| d dtf fdd	}|S )Nrf   r   c                    s   t | |}t| |f }|dkr)t| || }t|D ]
}d|||| f< q|S t| | |}t|D ]
}d||| |f< q4|S rg  )r  r  re  r  r3   )r  r  r   rf   Z_Mr  drZ   r  r;   r<   r     s   
znumpy_eye.<locals>.impl)
r}   r   r  r  rf   r  	DTypeSpecr  r   r  )r  r  r   rf   r   r;   r  r<   	numpy_eye  s   
r  c                 C   sD   t | s	tdt| tjr | jdvrtdddd}|S d S )N#The argument "v" must be array-like)r>   r  zInput must be 1- or 2-d.r   c                 S   s*  | j dkrG| j}|d t| }t||f| j}|dkr2t|| D ]}| | |||| f< q#|S t|| D ]}| | ||| |f< q8|S | j\}}|dk rT|| }|dkr\|| }tt||d}t	|| j}|dkrt|D ]}| ||| f ||< qs|S t|D ]}| || |f ||< q|S rY  )
r\   rR   r  r  re  rf   r3   r  r  r  )r   r   r   r  r  rZ   rowscolsr;   r;   r<   	diag_impl  s2   

zimpl_np_diag.<locals>.diag_implr  )r   r   r   r}   r   r~   r\   NumbaTypeError)r   r   r!  r;   r;   r<   impl_np_diag  s   



r#  c                 C   sB   t | sd}t|t|ttjfsd}t|ddd}|S )Nr  z#The argument "k" must be an integerr   c           	      S   s   t | } |  } t| }t|}|| }t ||f| j}t d| }t d|}t|D ]}| | ||| || f< q/|S r   )	r  r  rd  r   r  re  rf   maximumr3   )	r   r   r   Zabs_kr  r   rZ   jtr;   r;   r<   r     s   
znumpy_diagflat.<locals>.implr  )r   r   r   r}   r   r   r   )r   r   r   r   r;   r;   r<   numpy_diagflat  s   


r'  takec                    s   t | tjrt |tjrdd }|S tdd | |fD r*|jdk  fdd}|S t | tjr?t |tjtjfrAdd }|S d S d S )Nc                 S   s.   || j d ks|| j  k rtd|  | S )Nr>   Index out of bounds)rz   
IndexErrorrd  )r  r   r;   r;   r<   	take_impl#  s   znumpy_take.<locals>.take_implc                 s   r  r  r  r8  r;   r;   r<   r  )  r  znumpy_take.<locals>.<genexpr>r  c                    s   t j|j| jd} r| }n|}t |}d}|  }|D ]}|| jd ks.|| j k r2td|| ||< |d }q||j	S Nr  r   r>   r)  )
r  r  rz   rf   r  rX  rd  r*  r`  rR   )r  r   r  ZwalkeritrZ   rC  r  ZF_orderr;   r<   r+  ,  s   


c                 S   s~   t |}t j|j| jd}t |}d}|  }|D ]}|| jd ks*|| j k r.td|| ||< |d }q||j	S r,  )
r  r   r  rz   rf   rX  rd  r*  r`  rR   )r  r   convertr  r-  rZ   rC  r  r;   r;   r<   r+  ?  s   


)r}   r   r~   r   r  r  r  r  )r  r   r+  r;   r.  r<   
numpy_take  s   
r0  c                  G   s   dd | D }t dd |D rtj}|S t dd |D r#tj}|S ttddttj  }dd |D }t	||g }|S )	Nc                 S   s   g | ]
}t |tjs|qS r;   r  r  r;   r;   r<   r   O  r  z!_arange_dtype.<locals>.<listcomp>c                 s   r  r  r6  r  r;   r;   r<   r  Q  r  z _arange_dtype.<locals>.<genexpr>c                 s   r  r  )r}   r   Floatr  r;   r;   r<   r  S  r  zint%sr7  c                 S   s   g | ]}t |qS r;   )r   Z	unliteralr   r  r;   r;   r<   r   b      )
r  r   
complex128r   r  r  rf   r   r   r  )r   Zboundsrf   ZNPY_TYZunliteral_boundsr;   r;   r<   _arange_dtypeN  s   r5  c                    s  t |tjr	|j}t |tjr|j}t |tjr|j}|d u r"tj}|d u r)tj}|d u r0tj}t | tjrQt |tjtjfrQt |tjtjfrQt |tjtjfsSd S t |tjr`t| ||n|j	t
dd | ||fD t| dd  t|dd t|dd d fdd	}|S )Nc                 S   r  r;   r6  r2  r;   r;   r<   r         znp_arange.<locals>.<listcomp>r   c                    s    d ur n| }d urn|}d urn|}|d ur|nd}|d u r*d|}}	n||}}	|dkr7t d|	| | }
tt|
j}du rZtt|
j}tt||d}nt|d}t	|}|}t
|D ]
}|||  ||< qk|S )Nr>   r   zMaximum allowed size exceededT)r   r   mathceilr:  r<  r  r  r  r  r3   )r   rl  rf  rf   Z	lit_startZlit_stopZlit_stepZ_step_start_stopZnitems_cZnitems_rZnitems_ir   r  rr   rZ   Zstart_valueZ
step_valueZ
stop_valueZ
true_dtypeZuse_complexr;   r<   r     s(   

znp_arange.<locals>.implr  )r}   r   Optionalr4   r|  r  r  r  r5  rf   r  r  )r   rl  rf  rf   r   r;   r;  r<   	np_arangeh  s>   r=  c                    sn   t dd | |fD sd S t|ttjfsd}t|tdd | |fD r+tj ntj	 d fdd	}|S )	Nc                 s   r  r  )r}   r   r  r8  r;   r;   r<   r    r  z!numpy_linspace.<locals>.<genexpr>r5  c                 s   r  r  r6  r8  r;   r;   r<   r    r  r/  c                    s   t | }| d } |d }|dkr|S |d }|dkr7||  }t ||}td|D ]
}| ||  ||< q+n| |d< |dkrC||d< |S )Ng      ?r   r>   r.  )r  r  divider3   )r   rl  r2  r  divdeltarf  rZ   r  r;   r<   r0    s    z numpy_linspace.<locals>.linspacer3  )
r  r}   r   r   r   r   r   r  r4  r   )r   rl  r2  r   r0  r;   r  r<   numpy_linspace  s   
rA  c              
   C   s0  |j d }t|| ||d d}t||j}|j}t| |||}|j}	|j}
|jdv s.J |j|jkrBtj	||
|	|j
|jdd nLt||j}t||j}| tj}t|||)}t| ||	|||j|}t| ||
|||j|}|||| W d   n1 sw   Y  t| ||j| S )z
    Array copy.
    r   r   r  r>   rk   N)r   rd   r	   r   rR   r   r  r   r  Z
raw_memcpyr   r   r   r   r   rv   r  r  r   r7   r   r   )rh   r6   rL   r   r  r   r   rettyper  r  r  r  r  r   r   r  r  r;   r;   r<   _array_copy  s6   


	rC  c                 C   s@   t |tjsJ |jdkrdnd}|j|dd}||}|tfS )Nr  r   F)r  r  )r}   r   r~   r  r  rC  )r  r  r  r  rL   r;   r;   r<   _array_copy_intrinsic  s
   rD  z
array.copyc                 C   r  r  )rC  r  r;   r;   r<   
array_copy  r_  rE  c                 C   s   t | tjr
dd }|S )Nc                 S      t | S r  )rD  r{  r;   r;   r<   
numpy_copy
  r>  z#impl_numpy_copy.<locals>.numpy_copyr  )r  rG  r;   r;   r<   impl_numpy_copy  s   rH  c              	   C   s  |j }|jd }|j|ksJ d|jdkr`|jdksJ t|| ||d d}t|| |}| |ttjdd}	| 	|ttjd|j
f}
t||j|	|
|j
|j|j t| ||| S |j|jksp|jdkry|jdv ryt| |||d S |jdkr|dv sJ |dkrtnt}t|| |||d }||D\}}| t| |||d }|j}W d	   n1 sw   Y  | t| |||}|j}W d	   n1 sw   Y  W d	   n1 sw   Y  ||j}||| ||| |S t| |||S )
zd
    Common logic for layout conversion function;
    e.g. ascontiguousarray and asfortranarray
    r   z return-type has incorrect layoutr>   r   r  r  r  r   N)r   r   r  r\   rd   r  r   r  rv   ry  r   r   r   r   r   r   r   r   r   r  r  blockrC  phir4   Zadd_incoming)rh   r6   rL   r   output_layoutrZ  r   r   r  rR   r   Z
check_funcr  ZthenorelseZout_thenZthen_blkZ
out_orelseZ
orelse_blkZret_phir;   r;   r<   _as_layout_array  s`   


rM  c                    sH   t  tjst |j jt|jdd}|| }| fddfS )Nr>   )r  r\   c                    s   t | ||| jdS )NrK  )rM  r   )r]  r  r   r  rN  r;   r<   r  X  s    z,_as_layout_array_intrinsic.<locals>.<lambda>)	r}   r   r"   r   RequireLiteralValuer  r   r  r\   )r  r  rK  r  rL   r;   rN  r<   _as_layout_array_intrinsicP  s
   

rP  c                 C   sH   t | s	tdt| tjtjfrdd }|S t| tjr"dd }|S )Nrc  c                 S      t t | S r  )r  ascontiguousarrayr   r{  r;   r;   r<   r   b  rh  z%array_ascontiguousarray.<locals>.implc                 S   rB  Nr   rP  r{  r;   r;   r<   r   e  rS  r   r   r   r}   r   r  rp  r~   r  r;   r;   r<   array_ascontiguousarray\  s   
rV  c                 C   sL   t | s	tdt| tjtjfrdd }|S t| tjr$dd }|S d S )Nrc  c                 S   rQ  r  )r  asfortranarrayr   r{  r;   r;   r<   r   p  rh  z"array_asfortranarray.<locals>.implc                 S   rB  )Nr  rT  r{  r;   r;   r<   r   t  rS  rU  r  r;   r;   r<   array_asfortranarrayj  s   
rX  zarray.astypec              
   C   s  |j d }t|| ||d d}t||j}|j}t| |||}|j}	|j}
t||j}t||j}| 	t
j}t|||9}t| ||	|||j|}t| ||
|||j|}t| |||}| |||j|j}t| |||| W d    n1 s~w   Y  t| ||j| S r  )r   rd   r	   r   rR   r   r  r   r   r   r   rv   r  r  r  rp   ru   rf   rs   r   r   )rh   r6   rL   r   r  r   r   rB  r  r  r  r  r  r   r   r  r  itemr;   r;   r<   array_astypey  s0   


rZ  c                 C   r  )Nc              	   S   s2  |j d }|j}t|| ||d d}t|}|| |}|j}	t| |}
t|jj|
}||j	|j}|
||}t||}|j|dd d}| j|t|f W d    n1 s\w   Y  t||||g}t||g}||j| |	d}t||||||j|jd | }t| ||j|S )Nr   r   Fr  z.buffer size must be a multiple of element sizer   r  )r   r   rd   r   rj   r   r   r4   r   r   Zsremr	   Zis_not_nullr   r  r   r   rW   Zsdivr  r   r   r   r   r   r   r   r   )rh   r6   rL   r   Zbuftyr   r:  Z
out_ary_tyout_aryZout_datamodelr   ll_itemsizer  remZis_incompatibler   rR   r   r   r   r;   r;   r<   r    s<   


znp_frombuffer.<locals>.codegenr  )r  bufferrf   rZ  r  rL   r  r;   r;   r<   np_frombuffer  s   $r_  c                    s   t d| t| tjr| jdkrd|  }t||tu s-t|tjr)|j	tu s-t
|r1tj}nt|}|d urEtj|dd| j d nd|  d| d}t|tf fd	d
	}|S )N
frombufferr   z+Argument "buffer" must be buffer-like. Got r>   )rf   r\   r  r  z3Cannot parse input types to function np.frombuffer(r  r  c                    r  r  )r_  )r^  rf   r  r;   r<   r     rO  z impl_np_frombuffer.<locals>.impl)r  r}   r   r   r  r   r   r  r  r  r   r  r  r~   r  )r^  rf   r   r  r   r;   r  r<   impl_np_frombuffer  s0   






ra  c                    P   t |rtdd  d fdd	}|S t|tjr&td| d fdd	}|S d S )Nr   c                    
    | |S r  r;   rX   rR   rf   intrinsic_cfarrayr;   r<   r     rS  zimpl_carray.<locals>.implc                    rc  r  r;   rd  re  r;   r<   r     rS  r  r   get_cfarray_intrinsicr}   r   r  rX   rR   rf   r   r;   re  r<   impl_carray     

rj  c                    rb  )Nr  c                    rc  r  r;   rd  re  r;   r<   r     rS  zimpl_farray.<locals>.implc                    rc  r  r;   rd  re  r;   r<   r     rS  r  rg  ri  r;   re  r<   impl_farray  rk  rl  c                    s   t  fdd}|S )Nc           	         s   |t ju rd }nt|t jr|j}nd| d}t| d u r/|d u r,d}t||}n+t t jrO j}|d urN||krNd| d| d}t|nd  d}t|t|}|d u rmd| d}t|t 	||}t
|||}|tfS )Nz pointer argument expected, got ''z*explicit dtype required for void* argumentzmismatching dtype 'z' for pointer type 'zinvalid dtype spec 'zinvalid shape ')r   r  r}   CPointerrf   r   r"  r  r  r~   r!   
np_cfarray)	r  rX   rR   Z	ptr_dtyper   rf   r\   rZ  rL   dtype_r  r;   r<   rf    s4   





z0get_cfarray_intrinsic.<locals>.intrinsic_cfarrayr&   )r  rq  rf  r;   rp  r<   rh     s   rh  c                    s<  |j dd \}}|dd \}}|j}|jdv sJ t| }	t|}
t|
}t|tj	r9t
 |}n|f}|f} fddt||D }|}g }|jdkrf|D ]}||  ||}qWnt|D ]}||  ||}qj|   ||j }t|	||||dd |	 }t |j|S )zR
    numba.numpy_support.carray(...) and
    numba.numpy_support.farray(...).
    Nr  r  c                    r  r;   r  )r   r  r   r   r;   r<   r   9  r  znp_cfarray.<locals>.<listcomp>r  r  )r   r   r  rd   rj   r	   r   r}   r   r  r   r   rV   r   r  r  r  re   rf   r  r   r   r   )rh   r6   rL   r   r  r  rX   rR   r   r[  r   r\  r   offr   r   r   r   r;   r   r<   ro  $  sJ   




	ro  c                 C   sN   t |tjr| tjt|S t |tjr%| tttj|}|||fS J r  )	r}   r   r  r   rv   r   r  rC   r!   )rh   r6   seqtyseqr  r;   r;   r<   _get_seq_sizeX  s   rv  c                    s2   |j  tjt|tj fdd}|S )zK
    Return a getitem() implementation that doesn't incref its result.
    c                    s$   | |} j r j| | |S r  )r   r   Zdecref)r6   r   r  rh   r  rZ  r;   r<   wrapj  s   
z$_get_borrowing_getitem.<locals>.wrap)rf   rC   r  r  r!   r   rv   )rh   rt  rx  r;   rw  r<   _get_borrowing_getitemb  s   ry  c                    st    tj}t|d fdd}g }||}}	t|D ]}
|
dkr+|||	\}}	|t ||	 qt|S )zF
    Compute the likely shape of a nested sequence (possibly 0d).
    r   c                    sN   t | tjrt| dkrdS | d  |dfS t| }| j| |ffS )Nr   NN)r}   r   r  r   r  ry  rf   )rt  ru  r  r6   rh   r   r;   r<   get_first_itemz  s   
z.compute_sequence_shape.<locals>.get_first_item)r   r   rv   r   r3   rV   rv  r   )rh   r6   r\   rt  ru  r   r|  r   innertyinnerrZ   r;   r{  r<   compute_sequence_shapes  s   

r  c                    s0   fdd  fdd||| dS )z?
    Check the nested sequence matches the given *shapes*.
    c                      s   j  td d S )N)zincompatible sequence shape)r  r   r   r;   r   r;   r<   _fail  s   
z#check_sequence_shape.<locals>._failc                    s<  t |dkrd S t| |}|d }d||}j|dd    W d    n1 s/w   Y  t |dkr<d S t| tjrut| }t	|}| j
}|||jf}	||	|dd   W d    d S 1 snw   Y  d S t| tjrtt | D ]}
| |
 }||
}	||	|dd   qd S J | )Nr   r  Fr  r>   )r   rv  r   r   r}   r   r  ry  r	   r[  rf   r   r  r3   r  )rt  ru  r   rz   expectedr  r  r\  r}  r~  rZ   r  r6   check_seq_sizerh   r;   r<   r    s0   
"z,check_sequence_shape.<locals>.check_seq_sizeNr;   )rh   r6   rt  ru  r   r;   r  r<   check_sequence_shape  s   r  c                    s:    fddfdd||d dS )zl
    Assign a nested sequence contents to an array.  The shape must match
    the sequence's structure.
    c              
      sB   t j j| dd}|| j}t || d S rm  )r	   r  r  ru   rf   rs   )r   r  rr   rX   )r  r6   rh   r   r   r   r;   r<   assign_item  s
   z-assign_sequence_to_array.<locals>.assign_itemc                    s  t |dkrt| tjtjfrJ || | d S |d }t| tjr[t| }t|#}| j}|||j	f} |||dd  ||j	f  W d    d S 1 sTw   Y  d S t| tjrt
t | D ]!}	| |	 }||	}tj|	}
 |||dd  ||
f  qgd S J | rg  )r   r}   r   r  r  ry  r	   r[  rf   r   r3   r  r   rv   )rt  ru  r   r   rz   r  r\  r}  r~  rZ   r   )assignr  r6   rh   r;   r<   r    s(   
 "z(assign_sequence_to_array.<locals>.assignr;   Nr;   )rh   r6   r   r   r   r  rt  ru  r;   )r  r  r  r6   rh   r   r   r   r<   assign_sequence_to_array  s   r  c                 C   s>   t | |\}}t|r|}n
t|}|d u rd S t||dS rS  )r0   r   r  r   r~   )r  r  rf   r\   Z	seq_dtyper;   r;   r<   np_array_typer  s   r  c                 C   s0   t d| t| ||}|||}dd }||fS )Nr   c           
   	   S   s   |j }|j}|jd }|d }t| ||||}t||ksJ t| |||| t| |||}	t| ||	j||	j	||| t
| ||j |	 S r   )r   r\   r   r  r   r  r  r  r   r   r   r   )
rh   r6   rL   r   r  r\   rt  ru  r   r  r;   r;   r<   r    s   


znp_array.<locals>.codegen)r  r  )r  objrf   r  rL   r  r;   r;   r<   np_array  s
   

r  c                 C   sL   t d| t| stdt|st|d u rd}t|ddd}|S )Nr   z(The argument "object" must be array-likez:The argument "dtype" must be a data-type if it is providedc                 S   
   t | |S r  )r  )r  rf   r;   r;   r<   r     rS  zimpl_np_array.<locals>.implr  )r  r   r   r   r   r  )r  rf   r   r   r;   r;   r<   impl_np_array  s   



r  c           
      C   s   | d}| |}|d||}||||||}||d|||d||}|j|dd d| }	| j|t|	f W d    |S 1 sLw   Y  |S )Nr   r   rC  Fr  z%s(): axis out of bounds)	r4   r   rk  rR  r|  r   r  r   r*  )
rh   r6   r   r\   r   r   Zll_ndimZis_neg_axisZaxis_out_of_boundsr   r;   r;   r<   _normalize_axis  s   


r  c              
   C   s   t ||d ks
J ttj|}t||}td}t|d D ])}t|}	|d|	|}
|| }||
|	|	||	}|
|t||d| q"|
|t||d| t|||S )z
    Compute shape with the new axis inserted
    e.g. given original shape (2, 3, 4) and axis=2,
    the returned new shape is (2, 3, 1, 4).
    r>   rC  r   r   r   r  r	   r   rH   r3   r   rk  rR  r   rU   r   r7   )rh   r6   
orig_shaper\   r   ll_shtyr   r   r3  ll_dim
after_axisr  ry   r;   r;   r<   _insert_axis_in_shape*  s   


r  c              
   C   s   t ||d ks
J ttj|}t||}td}td}t|d D ]'}	t|	}
|d|
|}|||	|
||
}|
||	 t||d| q'|
|t||d| t|||S )zD
    Same as _insert_axis_in_shape(), but with a strides array.
    r>   r   rC  r  )rh   r6   orig_stridesr\   r   r  r   r   r   r3  r  r  ry   r;   r;   r<   _insert_axis_in_stridesG  s"   




r  c              	   C   s   |j }|j}|jd }t|| ||d d}t|| |}	t||j}
t||j}t| ||
||}t	| ||||}t
|	|j|||j|j|jd |	 S )z/
    np.expand_dims() with the given axis.
    r   r   r  )r   r\   r   rd   r	   r   rR   r   r  r  r   r   r   r   r   r   )rh   r6   rL   r   r   rZ  r\   r  r  r  r   r   Z
new_shapesZnew_stridesr;   r;   r<   expand_dimsd  s$   
r  c                 C   sB   |j dkr|jnd}|j|j d |d}|||}dd }||fS )Nr>   r  r  c                 S   sP   |  ||d |jd tj}t| |d|jj|}t| ||||}t| ||j|S )Nr>   znp.expand_dims)	ru   r   r   rv   r  r   r\   r  r   )rh   r6   rL   r   r   r  r;   r;   r<   r    s   znp_expand_dims.<locals>.codegen)r\   r  r  )r  r  r   r  r  rL   r  r;   r;   r<   np_expand_dims  s
   
r  c                 C   sL   t | tjsd|  }t|t |tjs d| }t|dd }|S )Nz)First argument "a" must be an array. Got z(Argument "axis" must be an integer. Got c                 S   r  r  )r  )r  r   r;   r;   r<   r     rS  z!impl_np_expand_dims.<locals>.impl)r}   r   r~   r   r   r   )r  r   r   r   r;   r;   r<   impl_np_expand_dims  s   



r  c                    s   t  fddfddS )Nc                    sX   | fdd|D  fdd}t tdkr tnd tj||fS )Nc                    s    g | ]}|j t|j d qS )r  )r  r  r\   r8  )minimumr;   r<   r     r  z-_atleast_nd.<locals>.impl.<locals>.<listcomp>c                    sn   t t |d } fddt|D }t|dkr+ |j|}n|d }t |j|S )Nr   c                    s"   g | ]\}}} |||qS r;   r;   )r   r  r  rZ  r6   rh   	transformr;   r<   r     s    z>_atleast_nd.<locals>.impl.<locals>.codegen.<locals>.<listcomp>r>   )_atleast_nd_transformr	   r   r   r   ry  r   r   )rh   r6   rL   r   arrsZretsr  )r~  r)  r  r  r  r<   r    s   

z*_atleast_nd.<locals>.impl.<locals>.codegenr>   r   )r!   r   r   r  ZStarArgTupleZ
from_types)Ztypingcontextr   r  )r)  r  )r~  r  r<   r     s   
z_atleast_nd.<locals>.implc                     s    |  S r  r;   r  )r   r;   r<   r    s    z_atleast_nd.<locals>.<lambda>rr  )r  r)  r;   )r)  r   r  r<   _atleast_nd  s   r  c                    s"   t  ksJ  fdd}|S )z`
    Return a callback successively inserting 1-sized dimensions at the
    following axes.
    c           	         sb   t D ]*}|d }|j|k r.t | }|j|jd d}t| |t|||f|}|}q|S )Nr>   r  )r3   r\   r	   r   r  r  r   r!   )	rh   r6   r  r  rZ  rZ   r\   r   Znewarrtyr)  min_ndimr;   r<   r    s   
z(_atleast_nd_transform.<locals>.transform)r   )r  r)  r  r;   r  r<   r    s   r  c                  G   s"   t dd | D rtddgS d S )Nc                 s   r  r  r  r8  r;   r;   r<   r    r  z np_atleast_1d.<locals>.<genexpr>r>   r   r  r  r  r;   r;   r<   np_atleast_1d  s   r  c                  G   s$   t dd | D rtdddgS d S )Nc                 s   r  r  r  r8  r;   r;   r<   r    r  z np_atleast_2d.<locals>.<genexpr>r  r   r  r  r;   r;   r<   np_atleast_2d     r  c                  G   s$   t dd | D rtdg dS d S )Nc                 s   r  r  r  r8  r;   r;   r<   r    r  z np_atleast_3d.<locals>.<genexpr>rG  )r   r   r  r  r  r;   r;   r<   np_atleast_3d  r  r  c	                 C   s  t |t |  krt |  krt |ksJ  J td}	t| |||}
t||
j}g }|D ]2}|	}tt||D ]!\}\}}|d|	||}|
||}||||||}q=|| q2|
j}t|||||D ]_\}}}}}|j}tj||tj|jd}|9}t| |||||j|}t| |||}| |||j|j}t| |||||j|}t| |||| W d   n1 sw   Y  t|||}qp|
S )z2
    Concatenate arrays along the given axis.
    r   rz  )orderN)r   r	   r   r  r   r   r  r   r   r4   r   rk  rR  rV   r   r  r  r  rp   ru   rf   rs   r   )rh   r6   r   r~  r  
arr_shapesarr_stridesrZ  
ret_shapesr   r  Zret_stridesZcopy_offsetsZarr_shr   r3  rz   r  is_axisaddendZret_datar  r  Zarr_stZarr_datar  r   r  rr   r  r;   r;   r<   _do_concatenate  sL   6



r  c                    s  |j } fddt||D }t d||} fdd|D } fdd|D } fdd|d D }	t|D ] d||}
|	 } |}fd	d|d
d  D } |
v\}}| t	 j
||g } || W d    n1 sw   Y  |@ tj}|D ]} | d||}q j |dd j td f W d    n1 sw   Y  W d    n1 sw   Y  W d    n1 sw   Y  q9 fdd|	D }	t |||||||		}t || S )Nc                    "   g | ]\}}t | |d qS r   r  r   r>  r  r   r;   r<   r   "  r  z#_np_concatenate.<locals>.<listcomp>np.concatenatec                       g | ]	}t  |jqS r;   )r	   r   rR   r   r  r  r;   r<   r   (  r  c                    r  r;   r	   r   r   r  r  r;   r<   r   )  r  c                    s   g | ]}t  |qS r;   )r	   r   r  r  r;   r<   r   .  r6  r   rz  c                    s   g | ]}|  qS r;   r;   r  )r3  r;   r<   r   5  rv  r>   Fr  z<np.concatenate(): input sizes over dimension %d do not matchc                    s   g | ]}  |qS r;   )r7   r  r  r;   r<   r   I  r3  )r\   r   r  r3   r   r4   r7   r  r  r  rR  r   r	   r  r  r   ra  r  r   r   r  r   r   )rh   r6   r~  r  rZ  r   r\   r  r  r  r  Zret_shape_ptrZret_shZother_shapesZon_axisZon_other_dimr  is_okr  r;   )r6   rh   r3  r<   _np_concatenate  s`   

r  c                    s  |j }td}td}tt|}	 fddt||D }t d||}t |d j}
|dd  D ]A}tj}tt |j|
D ]1\}} 	| 
d||} j |dd j td	 W d    n1 ssw   Y  qGq8 fd
d|D }ttj|t }t }t|d D ]5}t|} 
d||}|
| } | |||} |t |d|  |t |d| q |t |d|  |	t |d| t  |}|gt| }t  |} fddtt|D }t|d D ]6}t|} 
d||} | |||}tt|D ]} || | t || d| q6qtt|D ]} |t || d| qT fdd|D }t |||||||	}t || S )Nr   r>   c                    r  r  r  r  r   r;   r<   r   X  r  z_np_stack.<locals>.<listcomp>np.stackrz  Fr  )z5np.stack(): all input arrays must have the same shapec                    r  r;   r  r  r  r;   r<   r   j  r  rC  c                    r  r;   r  ru  )r6   r  r;   r<   r     r6  c                    s   g | ]}t   |qS r;   )r	   r   r7   )r   r  r  r;   r<   r     r  )r\   r	   r   r   r   r  r   rR   r  r  r   r   ra  r  r   r   r   r  rH   r3   rk  rR  r   rU   r7   r  r   r   )rh   r6   r~  r  rZ  r   r\   r   r   Z
ll_narraysr  r  r  r  Zorig_shr  Zinput_shapesr  r3  r  r  ry   Zinput_stridesrZ   r  r;   )r6   rh   r  r<   	_np_stackQ  s   








r  c                 C   sN   |d urt |tjsd S t| d|\}}|dkrtdt|}t|||S )Nr  r   z.zero-dimensional arrays cannot be concatenated)r}   r   r   r1   	TypeErrorr2   r~   )r  rz  r   rf   r\   r  r;   r;   r<   np_concatenate_typer  s   r  c                 C   s6   t | ||}t|tjsJ |||}dd }||fS )Nc              	   S   sF   |  ||d |jd tj}t| |t|jd t||d |j|S rY  )	ru   r   r   rv   r  r   r	   r   r   rh   r6   rL   r   r   r;   r;   r<   r    s   znp_concatenate.<locals>.codegen)r  r}   r   r~   )r  rz  r   r  rL   r  r;   r;   r<   np_concatenate  s
   
r  c                 C      t | tjrddd}|S d S )Nr   c                 S   r  r  )r  rz  r   r;   r;   r<   r     rS  z!impl_np_concatenate.<locals>.implr  r}   r   r  rz  r   r   r;   r;   r<   impl_np_concatenate     
r  c                 C   s*   |D ]}|j dk s|j dkrtdqdS )Nr>   r  z7np.column_stack() is only defined on 1-d and 2-d arrays)r\   r  )rh   r   rz  r  r;   r;   r<   _column_stack_dims  s
   r  c                 C   sB   t | d|td\}}t|}t|||}||}dd }||fS )Nznp.column_stack)Zdim_chooserc                 S   s   t |jd }t||d }g }g }| tjd}t||D ]:\}	}
|	jdkr3|	|	 |	|
 q|	jdks:J |	j
dd}t||	}t| |||
f|}|	| |	| qt| ||||j|S )Nr   r>   r  r  )r   r   r	   r   r   r   rv   r   r\   rV   r  r   r!   r  r  r   )rh   r6   rL   r   Zorig_arrtysZ	orig_arrsr~  r  r   r  r  r  
expand_sigZnewarrr;   r;   r<   r    s$   



z np_column_stack.<locals>.codegen)r1   r  r2   r   r~   )r  r  rf   r\   r  r  rL   r  r;   r;   r<   np_column_stack  s   
r  c                 C      t | tjrdd }|S d S )Nc                 S   rF  r  )r  r  r;   r;   r<   r     r>  zimpl_column_stack.<locals>.implr  r  r   r;   r;   r<   impl_column_stack     r  c              	   C   s*   t | |t|jd t||d |j|S )z/
    np.stack() with the given axis value.
    r   )r  r   r   r	   r   r   r  r;   r;   r<   _np_stack_common  s   r  c           	      C   sV   t | d|\}}tdd |D rdnd}t||d |}|||}dd }||fS )	Nr  c                 s   s    | ]}|j d kV  qdS )r  N)r  r  r;   r;   r<   r    s    z"np_stack_common.<locals>.<genexpr>r  r   r>   c                 S   s,   |  ||d |jd tj}t| ||||S rN  )ru   r   r   rv   r  r  r;   r;   r<   r    s   z np_stack_common.<locals>.codegen)r1   r  r   r~   )	r  rz  r   rf   r\   r  r  rL   r  r;   r;   r<   np_stack_common  s   
r  c                 C   r  )Nr   c                 S   r  r  )r  r  r;   r;   r<   r   )  rS  zimpl_np_stack.<locals>.implr  r  r  r;   r;   r<   impl_np_stack&  r  r  c                 C   s4   t | ||\}}t||}t|}t|||}|S r  )r1   r  r2   r   r~   )r  r   rz  Zndim_minrf   r\   r  r  r;   r;   r<   NdStack_typer.  s
   
r  c                 C   &   t | d|d}||}dd }||fS )Nz	np.hstackr>   c                    sf   |j d }|d j}|dkr| tjd t| ||| S |dkr#dnd  fdd}| ||||S )Nr   r>   c                    s   t j|  dS rn  rp  rz  ro  r;   r<   np_hstack_implJ  r  z3_np_hstack.<locals>.codegen.<locals>.np_hstack_implr   r\   r   r   rv   r  r  )rh   r6   rL   r   r  r\   r  r;   ro  r<   r  <  s   

z_np_hstack.<locals>.codegenr  r  r  r  rL   r  r;   r;   r<   
_np_hstack7     r  c                 C   r  )Nc                 S   rF  r  )r  r  r;   r;   r<   r   U  r>  zimpl_np_hstack.<locals>.implr  r  r;   r;   r<   impl_np_hstackR  r  r  c                 C   r  )Nz	np.vstackr  c                 S   sd   |j d }|d j}|dkrdd }n|dkr&| tjd}t| ||||S dd }| ||||S )Nr   c                 S   s   t t | dS rN  )r  r  hstackr  r;   r;   r<   np_vstack_impld  rq  z3_np_vstack.<locals>.codegen.<locals>.np_vstack_implr>   c                 S      t j| ddS )Nr   ro  rp  r  r;   r;   r<   r  m  r  r  )rh   r6   rL   r   r  r\   r  r   r;   r;   r<   r  _  s   


z_np_vstack.<locals>.codegenr  r  r;   r;   r<   
_np_vstackZ  r  r  c                 C   r  )Nc                 S   rF  r  )r  r  r;   r;   r<   r   x  r>  zimpl_np_vstack.<locals>.implr  r  r;   r;   r<   impl_np_vstacku  r  r  c                 C   r  )Nz	np.dstackrG  c                 S   s   |j d }|j}|d j}|dkrdd }| ||||S |dkrY| tjd}|j|jd d}	tj	|	g|j R  }
t
| ||
||}| tjd}t	||	}t| |||f|S |dkrl| tjd}t
| ||||S dd }| ||||S )Nr   c                 S   s   t | dddS )Nr>   r.  )r  r  r`  r  r;   r;   r<   r    r]  z3_np_dstack.<locals>.codegen.<locals>.np_vstack_implr>   r  r  c                 S   r  )Nr  ro  rp  r  r;   r;   r<   r    r  )r   r   r\   r  r   r   rv   r  r   r!   r  r  )rh   r6   rL   r   r  rZ  r\   r  r   Zstack_rettyZ	stack_sigZ	stack_retr  r;   r;   r<   r    s*   


z_np_dstack.<locals>.codegenr  r  r;   r;   r<   
_np_dstack}  s   "r  c                 C   r  )Nc                 S   rF  r  )r  r  r;   r;   r<   r     r>  zimpl_np_dstack.<locals>.implr  r  r;   r;   r<   impl_np_dstack  r  r  fillc                 C   r_  )Nc                 S   s   || d d < d S r  r;   )r  rr   r;   r;   r<   	fill_impl  s   zarr_fill.<locals>.fill_implr;   )r  rr   r  r;   r;   r<   arr_fill  s   r  dotc                 C   r_  )Nc                 S   r  r  )r  r  )r  otherr;   r;   r<   dot_impl  rO  zarray_dot.<locals>.dot_implr;   )r  r  r  r;   r;   r<   	array_dot  r  r  c                 C   "   t | std|  dd }|S )NzCannot np.fliplr on %s typec                 S   s4   t | }|jdk rtd|d d d d ddf S )Nr  rE  r.  .r  r  r\   r   r  r  r;   r;   r<   r     s   

znp_flip_lr.<locals>.implr  r  r   r;   r;   r<   
np_flip_lr     r  c                 C   r  )NzCannot np.flipud on %s typec                 S   s.   t | }|jdk rtd|d d ddf S )Nr>   zInput must be >= 1-d.r.  .r  r  r;   r;   r<   r     s   

znp_flip_ud.<locals>.implr  r  r;   r;   r<   
np_flip_ud  r  r  c                    sN   t |tjst|t|j tjtj d|} fdd}||fS )zY Creates a tuple of slices for np.flip indexing like
    `(slice(None, None, -1),) * sz` r  c                    sZ   dd }t jg}tjg|R  }| t j}| }| |g}	| ||||	}
|
S )Nc                 S   s*   |}t | D ]}t||td d d}q|S )Nr.  )r3   r+   r  )lengthempty_tupler   rZ   r;   r;   r<   r     s   z6_build_flip_slice_tuple.<locals>.codegen.<locals>.impl)r   rv   r   r!   r   Zget_constant_undefr  )rh   r6   r!   r   r   Zinner_argtypesZ	inner_sigZll_idx_typer  Z
inner_argsr   rz   Z
tuple_typer;   r<   r    s   

z(_build_flip_slice_tuple.<locals>.codegen)	r}   r   r  r   rO  r   r   r  Zslice3_type)r  szrL   r  r;   r  r<   _build_flip_slice_tuple  s   

r  c                 C   s&   t | tjstd|  dd }|S )NzCannot np.flip on %s typec                 S   s   t | j}| | S r  )r  r\   )r  slr;   r;   r<   r     s   
znp_flip.<locals>.impl)r}   r   r~   r   r   r  r;   r;   r<   np_flip  s   r  c                 C   s   t | tjtjtjfrd	dd}|S t |tjrd	dd}|S t |tjr4t |jjtjr4d	dd}|S t |tj	rKt
dd |jD rMd	dd}|S d S d S )
Nr   c                 S      t jt | ||dS rn  )r  array_splitr  r   indices_or_sectionsr   r;   r;   r<   r     s
   znp_array_split.<locals>.implc                 S   sP   t | j| |\}}tt|d g| |g|| d   }tj| ||dS )Nr>   ro  )divmodrR   r  Zcumsumr   r  )r   r  r   rs  r]  r   r;   r;   r<   r     s   c              	   S   sp   t | j}tdd| j|}g }d}|D ]}t||t||}|| |  |}q|| t||t|d   |S Nznp.splitr   r   )r,   r\   r   r+   r  rV   r   r  r   Z	slice_tupr   prevr  ry   r;   r;   r<   r   -  s   
c                 s   r  r  r}   r   r   )r   r&  r;   r;   r<   r  =  r  z!np_array_split.<locals>.<genexpr>c              	   S   st   t | j}tdd| j|}g }d}t|D ]}t||t||}|| |  |}q|| t||t|d   |S r  )r,   r\   r   r   r+   r  rV   r  r;   r;   r<   r   ?  s   
r  )r}   r   r  ListTyper  r   ZIterableTypeZiterator_typer   r  r  r   r  r   r   r;   r;   r<   np_array_split  s,   





r  c                 C   sL   t | tjtjtjfrddd}|S t |tjrddd}|S t| ||dS )Nr   c                 S   r  rn  )r  splitr  r  r;   r;   r<   r   S  s   znp_split.<locals>.implc                 S   s4   t | j| |\}}|dkrtdtj| ||dS )Nr   z0array split does not result in an equal divisionro  )r  rR   r   r  r  )r   r  r   r  r]  r;   r;   r<   r   Y  s   ro  r  )r}   r   r  r  r  r   r  r  r;   r;   r<   np_splitN  s   


r   c                 C   N   t | tjsd}t|t |tjtjtjtjfs!d}t|dd }|S )N#The argument "ary" must be an array:The argument "indices_or_sections" must be int or 1d-arrayc                 S   "   | j dk r	tdtj| |ddS )Nr  z3vsplit only works on arrays of 2 or more dimensionsr   ro  r\   r   r  r  r   r  r;   r;   r<   r   t     
znumpy_vsplit.<locals>.implr}   r   r~   r   r   r   r  r  r   r  r   r   r;   r;   r<   numpy_vspliti     

r
  c                 C   r  )Nr  r  c                 S   s<   | j dkr	td| j dkrtj| |ddS tj| |ddS )Nr   z3hsplit only works on arrays of 1 or more dimensionsr>   ro  r  r  r;   r;   r<   r     s
   

znumpy_hsplit.<locals>.implr  r	  r;   r;   r<   numpy_hsplit}  s   

r  c                 C   r  )Nr  r  c                 S   r  )NrG  z3dsplit only works on arrays of 3 or more dimensionsr  ro  r  r  r;   r;   r<   r     r  znumpy_dsplit.<locals>.implr  r	  r;   r;   r<   numpy_dsplit  r  r  c                 C   s   | |k pt |ot |  S r  )r  isnanr<  r;   r;   r<   	lt_floats  s   r  Fc                 C   s   | ||f}zt | W S  ty?   | dkr%tj|rtnd|dd}|j}n| dkr7tj|r/tnd|d}|j}|t |< | Y S w )z6
    Get a sort implementation of the given kind.
    r(   NT)lt
is_argsortZis_np_arrayr)   )r  r  )	_sortsKeyErrorr(   Zmake_jit_quicksortr  Zrun_quicksortr)   Zmake_jit_mergesortZrun_mergesort)r  is_floatr  keyr  funcr;   r;   r<   get_sort_func  s(   



r  z
array.sortc                    s<   |j d }tdt|jtjd  fdd}| ||||S )Nr   r(   )r  r  c                    s    |  d S r  r;   r  Z	sort_funcr;   r<   array_sort_impl  ru  z#array_sort.<locals>.array_sort_impl)r   r  r}   rf   r   r1  r  )rh   r6   rL   r   r  r  r;   r  r<   
array_sort  s   
r  c                 C   r  )NzArgument "a" must be array-likec                 S   s   |   }|  |S r  )r  r  )r  r   r;   r;   r<   np_sort_impl  s   z"impl_np_sort.<locals>.np_sort_implr  )r  r  r;   r;   r<   impl_np_sort  s   
r  zarray.argsortc           	         sb   |j \}}t|jt|jtjdd  fdd}|j|j d d d}|d d }| ||||S )NT)r  r  r  c                    s    | S r  r;   r  r  r;   r<   array_argsort_impl  r>  z)array_argsort.<locals>.array_argsort_implr>   r  )	r   r  r   r}   rf   r   r1  replacer  )	rh   r6   rL   r   r  r  r  ZinnersigZ	innerargsr;   r  r<   array_argsort  s   
r  c                 C   s   |j |j ks|jdksJ |S )Nr  )r  r  )rh   r6   r  r  rr   r;   r;   r<   array_to_array  s   r   c                 C   4   dd }t ||}| ||||g}t| ||j|S )Nc                 S   s
   |  dS r   )r(  r{  r;   r;   r<   r   	  s   
zarray0d_to_scalar.<locals>.impl)r!   r  r   r   rh   r6   r  r  rr   r   rL   r   r;   r;   r<   array0d_to_scalar  s   
r#  c                 C   r!  )Nc                 S   s   t | d S rQ  )r  r{  r;   r;   r<   r     rO  z array_to_unichrseq.<locals>.impl)r!   r  r   r   r"  r;   r;   r<   array_to_unichrseq  s   
r$  c                 C   r&  )zR
    An intrinsic returning a derived array with the given shape and strides.
    r   r  rR   r   r;   r;   r<   reshape_unchecked  r+  r&  c                    s   dd   fdd}|S )Nc                 S   s   t | tjotdd | D S )Nc                 s   r  r  r  r  r;   r;   r<   r  *  r  z>type_reshape_unchecked.<locals>.check_shape.<locals>.<genexpr>)r}   r   r  r  rR   r;   r;   r<   rx  (  s   z+type_reshape_unchecked.<locals>.check_shapec                    sJ   t | tjsd S  |r |sd S t|t|krd S | jt|ddS )Nr  r  )r}   r   r~   r   r  r%  rx  r;   r<   typer,  s   z%type_reshape_unchecked.<locals>.typerr;   )rh   r)  r;   r(  r<   type_reshape_unchecked&  s   	r*  c                 C   s   |j d }|j}t|| ||d }t|| |}t||d }t||d }	t||j||	|j|jd |	 }
t
| |||
S )Nr   r>   r  r  )r   r   rd   r	   r   r   r   r   r   r   r   )rh   r6   rL   r   r   rZ  r   r   rR   r   r   r;   r;   r<   impl_shape_unchecked8  s   
r+  c                    s^   |d t jfv rtdd  ntdd  |d t jfv rtdtdd d	 fdd	}|S )
Nc                 S   r:  r  r'  r  rR   r;   r;   r<   r,  Q  r   zas_strided.<locals>.get_shapec                 S   r  r  r;   r,  r;   r;   r<   r,  U     z*as_strided() strides argument is mandatoryc                 S   r  r  r;   )r  r   r;   r;   r<   get_strides`  r-  zas_strided.<locals>.get_stridesc                    s   t |  | || |} | S r  )r&  )r  rR   r   r,  r.  r;   r<   as_strided_impld  s   z#as_strided.<locals>.as_strided_implrz  )r   r|  r#   rD   )r  rR   r   r0  r;   r/  r<   
as_stridedN  s   

r1  c                    s  t |tjrtt| jd tt| jd tdd n/t |tjrIt |jtjrItt| jt	| tt| jt	| tdd nt
dt|rYtdd  n&t |tjrftdd  nt |tjrzt |jtjrztd	d  nt
d
d fdd	}|S )Nr>   c                 S   s   | fS r  r;   window_shaper;   r;   r<   get_window_shapeu  r   z-sliding_window_view.<locals>.get_window_shapec                 S   r  r  r;   r2  r;   r;   r<   r4  ~  r-  z4window_shape must be an integer or tuple of integersc                 S   s   t t|S r  )r   r3   r3  r   r\   r;   r;   r<   get_axis  ru  z%sliding_window_view.<locals>.get_axisc                 S   s   t dd||gS )Nsliding_window_viewr   r   r5  r;   r;   r<   r6    s   c                    s    fdd|D S )Nc                    s   g | ]	}t d d |qS )r7  r   r8  r  r  r;   r<   r     r  z9sliding_window_view.<locals>.get_axis.<locals>.<listcomp>r;   r5  r;   r  r<   r6    s   
z2axis must be None, an integer or tuple of integersc           
         s   |} ||| j }t|t|krtd}}t| j D ]}t||| j| }t||| j| }q | j }t||D ]7\}}|dk rItd|| |k rStd|| | d }t|||}t|||}t||| j| }|d7 }q=t| ||}	|	S )Nz2Must provide matching length window_shape and axisr   z-`window_shape` cannot contain negative valuesz4window_shape cannot be larger than input array shaper>   )	r\   r   r   r3   r+   rR   r   r   r&  )
r  r3  r   Z	out_shapeZout_stridesrZ   r   r3  Ztrimmedr  r6  r4  Zshape_bufferZstride_bufferr;   r<   sliding_window_view_impl  s8   
z5sliding_window_view.<locals>.sliding_window_view_implr  )r}   r   r   r   r3   r\   r#   r  rf   r   r   r   r   )r  r3  r   r:  r;   r9  r<   r7  k  s>   'r7  c                 C   r  )Nc                 S   s2   | j dkrdS | j dkrt| dS d}t|)Nr   Fr>   z[The truth value of an array with more than one element is ambiguous. Use a.any() or a.all())rz   r  r(  r   )r  r   r;   r;   r<   r     s   

zol_bool.<locals>.implr  r  r;   r;   r<   ol_bool  s   	r;  c                    sn   t |ttjfstdt |ttjfstdt | tjs%td| jtt	  fdd}|S )Nz.The second argument "axis1" must be an integerz-The third argument "axis2" must be an integerz'The first argument "a" must be an arrayc                    s`   t dd|}t dd|}|dk r|7 }|dk r|7 }t ||}t|||}t| |S )Nznp.swapaxesaxis1axis2r   )r   r+   r  r'  )r  r<  r=  Z
axes_tupleZ	axes_listr\   r;   r<   r     s   znumpy_swapaxes.<locals>.impl)
r}   r   r   r   r   r   r~   r\   r   r3   )r  r<  r=  r   r;   r>  r<   numpy_swapaxes  s   


r?  c                 C   s  t dd| j|}t| j}d||< tt||jD ]$\}\}}	|dkr&|	}
n|	dkr-|}
n
||	kr5td|}
t|||
}qt||| j| }t	| |} t	||}|}t
|dkrmtt
|D ]}t||| j| }qa|}t
|dkrtt
|D ]}t||| j|d |  }q{|j| }t||f | | j}td d d f}t|D ]1}t|D ])}| || |  }||| |  }||| |  }t|D ]
}|||  ||< qqq|S )Nznp.take_along_axisr   r>   z*`arr` and `indices` dimensions don't matchr   )r   r\   r   rR   r  r   r   r+   r  r  r   r3   r  rf   r  r  )r  r   r   ZNi_origZNk_origindices_broadcast_shapeZ	arr_shaperZ   Zd1Zd2new_valZarr_broadcast_shapeNiNkJr   Znp_s_iikkZa_1dZ
indices_1dZout_1dr%  r;   r;   r<   _take_along_axis_impl  sT   

rG  c                    s  t | tjstdt |tjstdt |jtjs"tdt|r)d}n| j}||jkr6tdt	t
|jt|rIfdd}|S t|d t |tjsYtd	|j}|d
k re| j| }|d
k sn|| jkrstdt	t
| t	t
|d | j fdd}|S )Nz)The first argument "arr" must be an arrayz.The second argument "indices" must be an arrayz'The indices array must contain integersr>   z;`indices` and `arr` must have the same number of dimensionsc                    s   t |  |ddd S )Nr   r;   )rG  rv  r  r   r   )r@  r;   r<   take_along_axis_implJ  s   z1arr_take_along_axis.<locals>.take_along_axis_implr   zaxis must be a literal valuer   zaxis is out of boundsc                    s   t | || S r  )rG  rH  rB  rC  r@  r;   r<   rI  [  s   )r}   r   r~   r   r   rf   r   r   r\   r   r3   r   r  rK  r   )r  r   r   Zarr_ndimrI  r;   rJ  r<   arr_take_along_axis5  s>   






rK  T        c                 C   s   t | tjr2t | tjrddd}|S t | tjr ddd}|S t | tjr-ddd}|S tdt| ret | j	tjrDddd}|S t | j	tjrRdd	d}|S t | j	tjr`dd
d}|S tdtd)NTrL  c                 S   r  r  r;   r  r  nanr;   r;   r<   r   f  rM  znan_to_num_impl.<locals>.implc                 S   sF   t | r|S t | rt t| jS t | r!t t| jS | S r  )r  r  isneginffinfor4   r  isposinfr  rM  r;   r;   r<   r   j  s   


c                 S   s*   t j| j|d}t j| j|d}t||S )N)rN  )r  
nan_to_numr:  r<  complex)r  r  rN  r  r]  r;   r;   r<   r   s  s   
z4Only Integer, Float, and Complex values are acceptedc                 S   r  r  r;   rM  r;   r;   r<   r     rM  c           	      S   s   t | jj}t | jj}t | }|rt |n|}|j}t|j	D ]%}t 
|| r2|||< q$t || r>|||< q$t || rI|||< q$|S r  )r  rP  rf   r  r  r  r  rC  r3   rz   r  rO  rQ  )	r  r  rN  Zmin_infZmax_infx_outputZoutput_flatrZ   r;   r;   r<   r     s   


c                 S   sD   t | }|rt |n|}t j|jd|d t j|jd|d |S )NF)r  rN  )r  r  r  rR  r:  r<  )r  r  rN  rT  rU  r;   r;   r<   r     s
   
z4The first argument must be a scalar or an array-likeTrL  )
r}   r   r  r   r1  r7  r   r   r   rf   )r  r  rN  r   r;   r;   r<   nan_to_num_impla  s6   
;
7
.
"


rW  r  r3  r  r  r  )Frz  rV  (  r0  r  r7  r  Zllvmliter   Zllvmlite.irr   numpyr  Znumbar   r   Z
numba.corer   r   r   r	   r
   Znumba.np.numpy_supportr   r   r   r   r   r   r   r   r   r   r   Znumba.core.imputilsr   r   r   r   r   r   r   r   r   r   r    Znumba.core.typingr!   Znumba.core.typesr"   Znumba.core.extendingr#   r$   r%   r&   r'   Z
numba.miscr(   r)   Znumba.cpythonr*   Znumba.cpython.unsafe.tupler+   r,   r-   Znumba.core.typing.npydeclr.   r  r/   r  r0   r1   r2   r=   r@   rd   rj   rp   rs   r|   r   r   r   r   r   r   r   r   r   ZArrayIteratorZBORROWEDr   r  r   r  r  r   r   r
  r  r  r  Anyr  r   r  r~   r!  r#  r  r%  r2  rK  rT  rX  rd  rn  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  Zbroadcast_shapesr  r  r  r  r	  r  r  r   ZVarArgr#  r'  r,  r
  r>  r4  Z	geomspacerD  Zrot90rL  rU  rY  r]  r^  r`  rb  resizerk  rV   rr  rz  rd  r}  r~  r  r  r  r  r  r  rR   r  uniquer  r  r  r  r  r  r  r  r  r  r  Z
MemoryViewr  r  r  r  r  r  r  r  r  r  r  r  r  rn  r  r  r  r  r  r   r  r  r
  r  r  r  ZDTypeZNumberClassr  r  r   r  r'  r(  r+  r  r,  r2  r3  r6  r9  Bytesr;  is_r?  rA  rG  rL  rV  rW  r  r  rE  r  r  r  r  r  r  rI  r  r  r  r  r  r  r  r  rX  r  ZNumpyNdIterTyper  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  re  r  r  r  r   r  r  r  r  r  r  r  r  r  r  r  eyer  Zdiagr#  Zdiagflatr'  r(  r0  r5  Zaranger=  r0  rA  rC  rD  rE  r  rH  rM  rP  rR  rV  rW  rX  rZ  r_  r`  ra  rj  rl  rh  ro  rv  ry  r  r  r  r  r  r   r  r  r  r  r  r  r  r  r  Z
atleast_1dr  Z
atleast_2dr  Z
atleast_3dr  r  r  r  r  r  rf  r  r  r  Zcolumn_stackr  r  r  stackr  r  r  r  r  r  Zvstackr  r  Zdstackr  r  r  rI  r  rJ  r  r  flipr  r  r  r  r   Zvsplitr
  Zhsplitr  Zdsplitr  r  r  r  r  r  r  Zargsortr  r   ZUnicodeCharSeqr1  r7  rp  ZNPTimedeltaZ
NPDatetimer#  r$  r&  Ztype_callabler*  r+  libZstride_tricksr1  r7  r  r;  rH  r?  rG  Ztake_along_axisrK  rR  rW  r;   r;   r;   r<   <module>   s   (4<	E



D


(




22=QE +
*=%



*

(
9 


6


F$
@


.








y~

	
,




+

*

































7	7'	"  JA X






O	0















&
.

C'%
A


+$4
))




=2Y
	
%
	


)	




>




^

;
+