U
    |e                 	   @   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 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.m/Z/m0Z0m1Z1m2Z2m3Z3 ddl4m5Z5m6Z6 ddl7m8Z8 ddl9m:Z:m;Z; ddl.m<Z< ddl=m>Z?m@ZAmBZBmCZCmDZD dd ZEdd ZFdd ZGdd ZHdd ZIdd ZJdd ZKd d! ZLd"d# ZMd4d$d%ZNd&d' ZOd(d) ZPe0ePd*d+ ZQed,ejRd-d. ZSd/d0 ZTed1ejUe%e)jVd2d3 ZWd5d4d5ZXd6d7 ZYd8d9 ZZeej[ejRej\eej[ejRej]d:d; Z^eej[ejRej_d<d= Z`eejaejRejbejbd>d? ZceedejRd@dA ZeedBejfdCdD ZgedEejfejbdFdG ZhG dHdI dIeiZjG dJdK dKejZkG dLdM dMejZlG dNdO dOejZmG dPdQ dQejZnG dRdS dSejZoG dTdU dUeiZpdVdW Zqeej[ejRejfdXdY ZrdZd[ Zsd\d] Ztd^d_ Zud`da Zvdbdc Zwddde Zxdfdg Zydhdi Zze2djdk Z{e2dldm Z|e/dndo Z}dpdq Z~e0e	jdrds Ze/dtdu Ze0e	jdvdw Ze0e	jdxdy Zdzd{ Zd|d} Zed~ejfdd Zdd Zed~ejfej_dd Zed~ejfeejbdd Ze0e	jd6ddZe ejfddd Ze0e	jd7ddZe0e	jd8ddZdd Zdd Zedejfej_dd Zedejfeejbdd Ze0e	jdd Ze0e	jd9ddZedejfdd Zee	jejfdd Zedejfdd Ze/dd Ze/dd Ze0e	jd:ddZe1ejfdd;ddZdd Ze0e	jdd Ze0e	jdd Ze0e	jdd Ze/dd Ze1ejfddd Ze2dd Zdd Ze0eddddĄ ZedejfejddǄ Ze ejfdȃddʄ Ze ejfd˃e ejd˃dd̈́ Ze ejfd΃e ejd΃ddЄ Ze ejfdуe ejdуddӄ Ze ejfdԃddք Ze ejfd׃e ejd׃ddل Ze ejfdڃe ejdڃdd܄ Ze ejd݃dd߄ Ze ejddd Ze ejddd Ze ejddd Ze ejfddd Ze ejddd Ze#ejeje#ejejÃdd Zdd Ze ejfddd Ze ejd݃e ejddd Ze ejddd Ze ejfddd Ze ejfddd Zːd d Ze1ejfde1ejfddd Z͐dd Ze ejϐde΃ e ejϐd	e΃ ed
ejejbdd Ze!ejfdd Zed
ejfe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jbdd Zedejejejbdd Ze$ejfdd Ze$ejԃd d! Ze$ejރd"d# Zeejejfejfd$d% Ze3ejfd&d'd( Zd)d* Zd+d, Zd<d-d.Zd=d/d0Zd1d2 Zd3d4 Zd5d6 Ze ejfd7d8d9 Zed1eje%e)jVd:d; Zeej[ejej\d<d= Zeejaejej\ejbd>d? Zeedejd@dA Zee	jejfdBdC Zed1eje%e)jVdDdE Zeeeej\ee	jeej\dFdG Zeeej_ee	jej_dHdI Zed1eje%e)jVdJdK Zee	jejbdLdM Zed1eje%e)jVdNdO ZdPdQ Ze<ejfdRdSdT ZdUdV Ze2dWdX Z dYdZ Zd[d\ Zd]d^ Zd_d` Ze2dadb Ze0e	jefdcddZe2dedf Z	e0e	j
d>dgdhZe2didj Ze1ejfdkdldm Ze0e	jefdndoZe0e	jd?dpdqZe0e	jd@drdsZe0e	jdAdtduZe0e	jdBdvdwZe0e	jdCdxdyZe0e	jdDdzd{Zd|d} Ze0ed~d Ze0e	jddefddZe0e	j dEddZ!e0e	j"e1ejfddd Z#dd Z$e0e	j%dFddZ&e0e	j'dGddZ(dd Z)e2dd Z*edejfdd Z+e0e	j,dd Z-dd Z.e2dd Z/e0e	j0dd Z1e0e	j2dd Z3edejfejedejfej-dd Z4e2dd Z5e0e	j6efddZ7e0edHddZ8e0edIddZ9dd Z:dd Z;dd Z<dd Z=dd Z>dd Z?dd Z@dd ZAe2dd ZBe0e	jCdJddZDdd ZEdd ZFdd ZGdÐdĄ ZHe2dŐdƄ ZIe0e	jHdǐdȄ ZJdɐdʄ ZKdːd̄ ZLe0e	jMd͐d΄ ZNe0e	jOdϐdЄ ZPe0e	jQdѐd҄ ZRdӐdԄ ZSdՐdք ZTdאd؄ ZUdِdڄ ZVe2dېd܄ ZWe0e	jXdKdݐdބZYdߐd ZZe2dd Z[e0e	j\dd Z]dd Z^e2dd Z_e0e	j`dLddZadd Zbe2dd Zce0e	jddd Zee2dd Zfe0e	jgdd Zhe2dd Zie0e	jjdd Zke1ejfddd Zle1ejfddd Zme0e	jndd  Zoe0e	jpdd Zqe2dd Zre0e	jsdd Zte0e	judMddZve0e	jwdNd	d
Zxi Zydd ZzdOddZ{edejfdd Z|e0e	j}dd Z~edejfej-ee	jejfej-dd Ze#ejfejfdd Ze#ejfeje#ejfeje#ejfej\e#ejfeje#ejfeje#ejfeje#ejfejdd Ze#ejfejdd Zdd Zeed d! Zeeejfej_ej_d"d# Ze0e	jjjdPd$d%Ze0e	jjjdQd&d'Ze0ed(d) Ze0e	jd*d+ Ze/d,d- Ze0e	jd.d/ Ze0e	jdRd2d3ZdS (S  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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moduleadd_metadataset_metadata)builderloadlower_boundupper_boundZrange_operandsmd r<   N/var/www/website-v5/atlas_env/lib/python3.8/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)r3   widthr>   )r7   r8   r:   r<   r<   r=   mark_positive6   s    rA   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k
r@   t | | Y S X | jj }| |}|d krttj	| j
|dd}n|}|| j
|f}||fS )N	__array__T)zfill)r    _contextget_functionNotImplementedErrorsuper
_make_refsdata_model_manager_get_be_typer	   alloca_once_builder)selfrefsigZ
array_impl	datamodelbe_type	outer_ref)ArrayStruct	__class__
array_typereal_array_typer<   r=   rH   L   s    


z*make_array.<locals>.ArrayStruct._make_refsc                    sp   | j }dkr | dS | d}g }tD ]2}t||d|}||}|| t|| q0t	||S )z\
            Override .shape to inform LLVM that its elements are all positive.
            r   shape)
rL   __getattr___get_ptr_by_namer2   r	   gep_inboundsr8   appendrA   
pack_array)rM   r7   ptrdimsiZdimptrr8   )basendimr<   r=   rW   `   s    


z%make_array.<locals>.ArrayStruct.shape)__name__
__module____qualname__rH   propertyrW   __classcell__r<   rS   rU   r`   ra   rV   )rT   r=   rS   J   s   rS   )as_arrayr	   create_struct_proxyra   )rU   r<   rg   r=   
make_array>   s
    
-rj   c                 C   s   |  |j}| |S )zB
    Return the item size for the given array or buffer type.
    )get_data_typedtypeget_abi_sizeof)contextrU   lltyr<   r<   r=   get_itemsizez   s    rp   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_valuerl   )rn   r7   arraytyr]   rr   r<   r<   r=   	load_item   s    rv   c                 C   s$   |j r
dnd}| j||j|||dS )z4
    Store the item at the given array pointer.
    Nr?   rq   )rs   
pack_valuerl   )rn   r7   ru   valr]   rr   r<   r<   r=   
store_item   s    ry   c                 C   s>   |j r(| |||tj}t|||}n| |||tj}|S )zM
    Fix the integer index' type and value for the given dimension size.
    )signedcastr   intpr)   Z	fix_indexuintp)rn   r7   idxtyidxsizeindr<   r<   r=   fix_integer_index   s
    r   c                 C   s^   t |tjrR|jdkrRt |jtjs(tt|| ||}t| |||j	}|j|fS ||fS dS )zS
    Normalize the index type and value.  0-d arrays are converted to scalars.
    r   N)

isinstancer   Arrayra   rl   IntegerAssertionErrorrj   rv   data)rn   r7   r~   r   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   ).0r~   r   r7   rn   r<   r=   
<listcomp>   s   z%normalize_indices.<locals>.<listcomp>)lenzip)rn   r7   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r|| }d| j	 dt
|dkrndnd d| d}t||d	krt||	d
d	}|tj}t|ttfrt|||}t|ttfrt|||}t|tr||}t|||||d}|d	kr2t||	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)rW   stridesr   itemsizer   parentnswflagsnitemszmissing {0}) rD   rL   
_datamodelr   r   float64rI   set_fields_fe_typer   
ValueErrorr   get_value_typeget_typer|   r   tuplelistr	   r\   intdictget_constantunpack_tupler3   countmulkeysformatitemssetattr)arrayr   rW   r   r   r   r   rn   r7   rP   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   s`    


*



r   c                 C   sh   |j }|j}|tjd}t||j| j}|D ]}|j	||dgd}q0||_
|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)rD   rL   r   r   r|   r	   r   rW   ra   r   r   rp   r   )arytyr   rn   r7   r   r   r   r<   r<   r=   update_array_info   s    r   c                 C   s
   t  dS )z/Constrain axis values to valid positive values.NrF   	func_namearg_namera   axisr<   r<   r=   normalize_axis  s    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   ra   r   r   r<   r   r=   normalize_axis_overloads  s    



r   getiterc                 C   st   |j \}|\}| ||j}| tjd}t||}||_||_	| j
rX| j||| | }	t| ||j|	}
|
S r   )argsmake_helperreturn_typer   r   r|   r	   alloca_once_valueindexr   
enable_nrtnrtincref	_getvaluer   )rn   r7   rO   r   ru   r   iterobjzero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
rt| |||||	|
|}| S |
rtt| |||	S dS )z1 Evaluate `ary[idx]`, where idx is a single int. r   r   r?   N)r	   r   rW   ra   r   r   pointer_addr   r   r   Buffer	make_viewr   r   rv   )rn   r7   r   r   aryr   shapesr   offsetdataptrview_shapesview_stridesretaryr<   r<   r=   _getitem_array_single_intE  s      r   iternextc              	   C   s   |j \}|\}|j}| j|||d}t|| ||jd}	tj||	jdd\}
||j	}|
d||
}|| ||> t| ||j||	|}|| t||}|||j	 W 5 Q R X d S )Nvaluer?   r   <)r   rU   r   rj   r   r	   r   rW   r8   r   icmp_signed	set_validif_thenr   
yield_typeyield_increment_indexstore)rn   r7   rO   r   resultitertyiterru   r   r   r   r   is_validr   nindexr<   r<   r=   iternext_arrayZ  s*    
     
r   c              	   C   sD  |  tjd}|  tjd}t||j|j}	t||j|j}
g }g }g }tdd |D }d}t	||D ]n\}}|tj
kr|jt| d | }t|D ]2}|| ||	|  ||
|  |d7 }qqlt|tjrT| j|||d}t| ||| t|||	|  ||j t||}t|||
| }|| || nt|tjrt| ||||	| }|rt| |||	| | || n6t|r|| || |d8 }ntd|f |d7 }ql||jkst||jk r"||	|  ||
|  |d7 }qt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   s   g | ]}t |r|qS r<   r   r   r   r<   r<   r=   r     s      z"basic_indexing.<locals>.<listcomp>r   unexpected index type: %sF
wraparoundboundscheck)r   r   r|   r	   r   rW   ra   r   r   r   ellipsisr2   r[   r   	SliceTyper   r)   guard_invalid_slice	fix_slicestartget_slice_lengthZ
fix_strider   r   do_boundscheckr   rF   r   get_item_pointer)rn   r7   r   r   r   r   r   r   oner   r   Zoutput_indicesZoutput_shapesZoutput_stridesnum_newaxesaxindexvalr~   	n_missingr_   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   rW   r   r   r   r   )rj   r   r   r   r   )	rn   r7   r   r   r   r   r   r   r   r<   r<   r=   r     s    r   c              	   C   sb   t | |||||| jd\}}}	t|tjrHt| |||||||	}
|
 S |rPtt| |||S dS )zp
    Return the result of indexing *ary* with the given *indices*,
    returning either a scalar or a view.
    r   N)	r  enable_boundscheckr   r   r   r   r   r   rv   )rn   r7   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 tt|| ||}t| ||j|||f|f}t| ||j|S )z4
    Basic indexing with an integer or a slice.
    r?   )r   ra   r   rj   r  r   r   	rn   r7   rO   r   r   r~   r   r   r   r<   r<   r=   getitem_arraynd_intp  s    

   r  c              	   C   s   |j \}}|\}}t|| ||}|j}tj||t|d}	t| |||	\}}	tdd |D rvt| |||||||	S t	| ||j
||||	}
t| ||j
|
S )z2
    Basic or advanced indexing with a tuple.
    r   c                 s   s   | ]}t |tjV  qd S Nr   r   r   )r   tyr<   r<   r=   	<genexpr>  s     z&getitem_array_tuple.<locals>.<genexpr>)r   rj   r   r	   r   r   r   anyfancy_getitemr  r   r   )rn   r7   rO   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k
r   d}Y n
X t|}|rt| ||||
|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   rj   r   r  r  rF   boolfancy_setslicer{   rl   ry   )rn   r7   rO   r   r   r~   valtyr   r   rx   r   r   r   r   r   Zuse_fancy_indexingr<   r<   r=   setitem_array	  s6    
 

 r#  c           	      C   sD   |j \}|\}t|}|| ||}|j}||d}t| ||j|S r   )r   rj   rW   extract_valuer   r   )	rn   r7   rO   r   r   r   ZarysttyZshapearyr   r<   r<   r=   	array_len2  s    r%  z
array.itemc              	   C   st   |j \}|\}t|| ||}|j}|j|d||ddd d}| j|t|f W 5 Q R X t	| |||j
S )N!=r?   Flikelyz>item(): can only convert an array of size 1 to a Python scalar)r   rj   r   r   r   r3   	call_convreturn_user_excr   rv   r   )rn   r7   rO   r   r   r   r   r   r<   r<   r=   
array_item=  s    r+  zarray.itemsetc           
   	   C   s   |j \}}|\}}||jks tt|| ||}|j}|j|d||ddd d}	| j	|t
|	f W 5 Q R X t| ||||j |  S )Nr&  r?   Fr'  z/itemset(): can only write to an array of size 1)r   rl   r   rj   r   r   r   r3   r)  r*  r   ry   r   get_dummy_value)
rn   r7   rO   r   r   r"  r   rx   r   r   r<   r<   r=   array_itemsetL  s    
r-  c                   @   s@   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   s   t dS )ze
        Prepare the indexer by initializing any required variables, basic
        blocks...
        Nr   rM   r<   r<   r=   prepareg  s    zIndexer.preparec                 C   s   t dS )z=
        Return this dimension's size as an integer.
        Nr   r/  r<   r<   r=   get_sizen  s    zIndexer.get_sizec                 C   s   t dS )z;
        Return this dimension's shape as a tuple.
        Nr   r/  r<   r<   r=   	get_shapet  s    zIndexer.get_shapec                 C   s   t dS )z}
        Return a half-open [lower, upper) range of indices this dimension
        is guaranteed not to step out of.
        Nr   r/  r<   r<   r=   get_index_boundsz  s    zIndexer.get_index_boundsc                 C   s   t dS )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.
        Nr   r/  r<   r<   r=   	loop_head  s    	zIndexer.loop_headc                 C   s   t dS )z)
        Finish indexation loop.
        Nr   r/  r<   r<   r=   	loop_tail  s    zIndexer.loop_tailN)
rb   rc   rd   __doc__r0  r1  r2  r3  r4  r5  r<   r<   r<   r=   r.  a  s   r.  c                   @   sH   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  )	rn   r7   r   r   dimr   r   r|   ll_intp)rM   rn   r7   r   r   r8  r<   r<   r=   __init__  s    zEntireIndexer.__init__c                 C   sB   | j }|| jj| j| _t|| j| _	|
 | _|
 | _d S r  )r7   r$  r   rW   r8  r   r	   rK   r9  r   append_basic_blockbb_startbb_endrM   r7   r<   r<   r=   r0    s
    
zEntireIndexer.preparec                 C   s   | j S r  r   r/  r<   r<   r=   r1    s    zEntireIndexer.get_sizec                 C   s   | j fS r  r?  r/  r<   r<   r=   r2    s    zEntireIndexer.get_shapec                 C   s   |  d| jfS r   r9  r   r/  r<   r<   r=   r3    s    zEntireIndexer.get_index_boundsc              	   C   s|   | j }| j t| jd| j || j || j || j}|j	|
d|| jdd || j W 5 Q R X ||fS Nr   >=Fr'  )r7   r   r   r9  r   branchr<  position_at_endr8   r   r   r   r=  )rM   r7   	cur_indexr<   r<   r=   r4    s    zEntireIndexer.loop_headc                 C   sD   | j }t||| j}||| j || j || j	 d S r  )
r7   r	   r   r8   r   r   rC  r<  rD  r=  rM   r7   
next_indexr<   r<   r=   r5    s
    zEntireIndexer.loop_tailNrb   rc   rd   r6  r:  r0  r1  r2  r3  r4  r5  r<   r<   r<   r=   r7    s   r7  c                   @   sH   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 )IntegerIndexerz0
    Compute indices from a single integer.
    c                 C   s&   || _ || _|| _| j tj| _d S r  )rn   r7   r   r   r   r|   r9  )rM   rn   r7   r   r<   r<   r=   r:    s    zIntegerIndexer.__init__c                 C   s   d S r  r<   r/  r<   r<   r=   r0    s    zIntegerIndexer.preparec                 C   s   t | jdS Nr?   )r   r9  r/  r<   r<   r=   r1    s    zIntegerIndexer.get_sizec                 C   s   dS Nr<   r<   r/  r<   r<   r=   r2    s    zIntegerIndexer.get_shapec                 C   s   | j | j| j |  fS r  )r   r7   addr1  r/  r<   r<   r=   r3    s    zIntegerIndexer.get_index_boundsc                 C   s
   | j d fS r  )r   r/  r<   r<   r=   r4    s    zIntegerIndexer.loop_headc                 C   s   d S r  r<   r/  r<   r<   r=   r5    s    zIntegerIndexer.loop_tailNrH  r<   r<   r<   r=   rI    s   rI  c                   @   sH   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 )IntegerArrayIndexerz;
    Compute indices from an array of integer indices.
    c                 C   s@   || _ || _|| _|| _|| _|jdks,t| j tj	| _
d S rJ  )rn   r7   r~   r   r   ra   r   r   r   r|   r9  )rM   rn   r7   r~   r   r   r<   r<   r=   r:    s    zIntegerArrayIndexer.__init__c                 C   sD   | j }t|| jjd | _t|| j| _|	 | _
|	 | _d S r   )r7   r	   r   r   rW   idx_sizerK   r9  	idx_indexr;  r<  r=  r>  r<   r<   r=   r0    s
    
zIntegerArrayIndexer.preparec                 C   s   | j S r  rN  r/  r<   r<   r=   r1    s    zIntegerArrayIndexer.get_sizec                 C   s   | j fS r  rP  r/  r<   r<   r=   r2    s    zIntegerArrayIndexer.get_shapec                 C   s   |  d| jfS r   r@  r/  r<   r<   r=   r3    s    z$IntegerArrayIndexer.get_index_boundsc              	   C   s   | j }| j t| jd| j || j || j || j}|j	|
d|| jdd || j W 5 Q R X t| j|| jj| j| j|}t| j|| jj|| j}||fS rA  )r7   r   r   r9  rO  rC  r<  rD  r8   r   r   rN  r=  r   rn   r~   rl   r   r   r   )rM   r7   rE  r   r<   r<   r=   r4    s0          zIntegerArrayIndexer.loop_headc                 C   sD   | j }t||| j}||| j || j || j	 d S r  )
r7   r	   r   r8   rO  r   rC  r<  rD  r=  rF  r<   r<   r=   r5    s    
zIntegerArrayIndexer.loop_tailNrH  r<   r<   r<   r=   rM    s   	rM  c                   @   sH   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 )BooleanArrayIndexerz>
    Compute indices from an array of boolean predicates.
    c                 C   sH   || _ || _|| _|| _|jdks&t| j tj| _	t
| j	d| _d S Nr?   r   )rn   r7   r~   r   ra   r   r   r   r|   r9  r   r   )rM   rn   r7   r~   r   r<   r<   r=   r:  &  s    zBooleanArrayIndexer.__init__c                 C   s^   | j }t|| jjd | _t|| j| _t|| j| _	|
 | _|
 | _|
 | _d S r   )r7   r	   r   r   rW   r   rK   r9  rO  r   r;  r<  bb_tailr=  r>  r<   r<   r=   r0  /  s    

zBooleanArrayIndexer.preparec              	   C   s   | j }t|| j}t|| jP}||}t| j|| j	j
| j	| j|j}|||||j}||| W 5 Q R X ||S r  )r7   r	   r   r   	for_ranger   r8   r   rn   r~   rl   r   r   rL  zextr3   r   )rM   r7   r   loopcpredr<   r<   r=   r1  8  s    
    zBooleanArrayIndexer.get_sizec                 C   s
   |   fS r  r1  r/  r<   r<   r=   r2  G  s    zBooleanArrayIndexer.get_shapec                 C   s   |  d| jfS r   r@  r/  r<   r<   r=   r3  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 5 Q R X t| j|| jj| j| j|}|	|| || j W 5 Q R X t||}||| j ||fS )NrB  Fr'  )r7   r   r   rO  r   rC  r<  rD  r8   r   r   r   r=  r   rn   r~   rl   r   not_rS  r	   r   )rM   r7   rE  	cur_countrX  
next_countr<   r<   r=   r4  O  s0        zBooleanArrayIndexer.loop_headc                 C   s\   | j }|| j || j t||| j}||| j || j	 || j
 d S r  )r7   rC  rS  rD  r	   r   r8   rO  r   r<  r=  rF  r<   r<   r=   r5  g  s    
zBooleanArrayIndexer.loop_tailNrH  r<   r<   r<   r=   rQ  !  s   		rQ  c                   @   sH   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 )SliceIndexerz(
    Compute indices along a slice.
    c                 C   sL   || _ || _|| _|| _|| _|| _|| _| j tj	| _
t| j
d| _d S r   )rn   r7   r   r   r8  r~   r  r   r   r|   r9  r   r   )rM   rn   r7   r   r   r8  r~   r  r<   r<   r=   r:  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  )r7   r$  r   rW   r8  dim_sizer)   r  rn   r~   r  r  r	   
is_neg_intstepis_step_negativerK   r9  r   r   r;  r<  r=  r>  r<   r<   r=   r0    s    
zSliceIndexer.preparec                 C   s   t | j| jS r  )r)   r  r7   r  r/  r<   r<   r=   r1    s    zSliceIndexer.get_sizec                 C   s
   |   fS r  rY  r/  r<   r<   r=   r2    s    zSliceIndexer.get_shapec                 C   s   t | j| j\}}||fS r  )r)   Zget_slice_boundsr7   r  )rM   lowerupperr<   r<   r=   r3    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 5 Q R X ||fS )Nz<=rB  Fr'  )r7   r   r  r  r   r   r   rC  r<  rD  r8   selectra  r   stopr   r=  )rM   r7   rE  r[  Zis_finishedr<   r<   r=   r4    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   )r7   rL  r8   r   r  r`  r   r	   r   r   rC  r<  rD  r=  )rM   r7   rG  r\  r<   r<   r=   r5    s    zSliceIndexer.loop_tailNrH  r<   r<   r<   r=   r]  r  s   r]  c                   @   s@   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 )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kr|jt| d | }t|D ].}t|||||	}|| |	d7 }	|
d7 }
qqxt|t
jr||||}t|||||	||}|| nt|t
jrTt||||| j|	 }t|||}|| nt|t
jrt||||}t|jt
jrt||||| j|	 }n*t|jt
jrt ||||}n
dst!|| n.t"|r| j|
 |	d8 }	nt!d|f |	d7 }	|
d7 }
qx|	|jks*t!|	|jf|	|jk r\t|||||	}|| |	d7 }	q*t||jks~t!t||jf|| _#d S )Nc                 S   s   g | ]}t |r|qS r<   r   r   r<   r<   r=   r     s      z)FancyIndexer.__init__.<locals>.<listcomp>r   r?   r   )$rn   r7   r   r	   r   rW   ra   r   r   r   r   r|   r9  newaxesr   r   r   r2   r7  r[   r   r   r   r]  r   r   rI  r   rj   rl   rM  BooleanrQ  r   r   indexers)rM   rn   r7   r   r   r   r   ri  r  r	  Znew_axr
  r~   r  r_   indexerr  r   r   r<   r<   r=   r:    sz    


 
  




"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   s   g | ]}|  qS r<   )r2  r   r_   r<   r<   r=   r     s     z(FancyIndexer.prepare.<locals>.<listcomp>r<   )
ri  r0  rn   r   r   r|   rg  insertsumindexers_shape)rM   r_   r  	res_shaper<   r<   r=   r0    s    


zFancyIndexer.preparec                 C   s   | j S )z?
        Get the resulting data shape as Python tuple.
        )rn  r/  r<   r<   r=   r2    s    zFancyIndexer.get_shapec              	   C   s  t || jjkst| j}tj}| d}| d}|}|}t| j	| j
|D ]\}	}
}|||d|
|}|	 \}}|||}|||||}|d||}||||||}||||||}qL|||}||||}||||}||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   ra   r   r7   r	   	false_bitr9  r   ri  rn  or_icmp_unsignedr3  r   subr   rL  rd  )rM   r   r   r7   is_emptyr   r  rb  rc  rj  rW   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  )r4  rk  r<   r<   r=   r  >  s     z+FancyIndexer.begin_loops.<locals>.<genexpr>)r   ri  )rM   r   countsr<   r<   r=   begin_loops=  s    zFancyIndexer.begin_loopsc                 C   s   t | jD ]}|  q
d S r  )reversedri  r5  )rM   r_   r<   r<   r=   	end_loopsA  s    zFancyIndexer.end_loopsN)
rb   rc   rd   r6  r:  r0  r2  rw  ry  r{  r<   r<   r<   r=   rf    s   B&rf  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   rW   r   r   rf  r0  r   r2  _empty_nd_implr   r   r   r|   ry  get_item_pointer2layoutr  rv   r8   gepry   r   r   r{  r   r   )rn   r7   rO   r   r   r   r   r   r   r   r   rj  Zout_tyZ
out_shapesr   Zout_dataZout_idx_r]   rx   curZnext_idxr<   r<   r=   r  F  s<    
   
r  c           	   	   C   s   |j \}}|\}}t|| ||}|jdkrht| |||\}}t| ||j|||f|f}t| ||j|S t| ||||||f|fS dS )z3
    Advanced or basic indexing with an array.
    r   N)r   rj   ra   r   r  r   r   r  r  r<   r<   r=   fancy_getitem_arrayq  s$    


   
   r  c              	   C   s   |j }|d}|d}|jdkr8|}	|||j}
n|}	|}
t|jD ]Z}||| ||| |}|d||}|	||
|
||
}
|	||	|
|	|}	qJ|
|
|}
|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?   CFrB  rp  )r   r3   r~  r   r   r2   ra   rt  r   rd  rL  )rn   r7   arrtyarrr   r   r   r   r  rb  rc  r_   Zmax_axis_offsetZ
is_upwardsru  r<   r<   r=   offset_bounds_from_strides  s4    




  
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.
    )ptrtointr3   rL  )rn   r7   rb  rc  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  )	rn   r7   r  r  r   r   r   rb  rc  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_rs  )rn   r7   a_starta_endZb_startZb_endZmay_overlapr<   r<   r=   extents_may_overlap  s
    r  c                    s  j dt d  jdd  |j|j}j |}	|	j  	|	}	 
|	 tj}
t |
F}t 	
j|}t |	|} 
 || W 5 Q R X W 5 Q R X  	
fdd} fdd}||fS )Nr   Fr'  c                    s   t  } jddv\}}|.  t j  | dd| W 5 Q R X |*  t j 	
j| dd| W 5 Q R X W 5 Q R X t 
 |S )NFr'  r   )r	   rK   if_elser   r}  r8   r~  rv   )source_indicessrc_ptrZif_copy	otherwise)r7   rn   	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               	      sB    j dd*  }  | } j |  W 5 Q R X d S )NFr'  )r   r8   bitcastr   freer   )r7   rn   r  r  	voidptrtyr<   r=   src_cleanup  s    
z&maybe_copy_source.<locals>.src_cleanup)r3   r	   r   r   r   r   r   r   allocater  r   r   r   r|   	loop_nestr}  r~  r8   )rn   r7   r  r  srcr  r  r  	allocsizer   r   r   r  dest_ptrr  r  r<   )r7   rn   r  r  r  r  r  r  r  r  r  r  r  r=   maybe_copy_source  s:    
  
  &"r  c           
   	      s   |  tjd}|  tjdt|t|krZt|t| }g| | }|g| | }nt|t|k rt|t| } fdd|d| D }t j|tj} j	 
|dd d}	| j t|	f W 5 Q R X ||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                    s   g | ]}  d |qS rp  )rs  r   r  r7   r  r<   r=   r   %  s   z(_bc_adjust_dimension.<locals>.<listcomp>NFr'  z,cannot broadcast source array for assignment)r   r   r}   r   	functoolsreducer  r	   true_bitr   rZ  r)  r*  r   )
rn   r7   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r7   r<   r=   r   @  s   z,_bc_adjust_shape_strides.<locals>.<listcomp>c                    s   g | ]}  d |qS r  r  )r   r  r  r<   r=   r   B  s     c                    s   g | ]\}}  ||qS r<   )r  )r   xyr  r<   r=   r   C  s     c                    s    g | ]\}}}  |||qS r<   rd  )r   pr  r  r  r<   r=   r   D  s   c                    s   g | ]\}}  ||qS r<   r  )r   r  r  )r7   r   r<   r=   r   F  s   )r   r   r}   r   )
rn   r7   r   r   r  Z	bc_shapesZ
bc_stridesmismatchZ
src_is_onepredsr<   )r7   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ra   r~  r  )r	   r   rW   r   r  r  copyr   rj   r   r   r\   r   r   r   )	rn   r7   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)ra   r~  readonlyc           
   	   S   sz   |\}}|j d }t|| ||}t||}t| ||||\}}t|dt| |j	dd  |
 }	t| ||j|	S )Nr   r   )r   rj   r	   r   r  r   r   r   r   r   r   r   r   )
rn   r7   rO   r   r  Zshape_r  r  destr   r<   r<   r=   codegenj  s    
z$_numpy_broadcast_to.<locals>.codegen)r  r   )	typingctxr   rW   retrO   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   rj   r	   get_null_valuer   r3   r   r   r   )rn   r7   rO   r   r  r  r  r   r<   r<   r=   r    s    
z#get_readonly_array.<locals>.codegen)r  )r  r  r  rO   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 r|| }|| }	||	kst|dkr|d7 }|d7 }qLtdqLd 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)rW   r   r   )
r   Z
dest_shapeZ	src_shapeZsrc_ndimZ	dest_ndimr   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   rW   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rVt|jtjsRd}t|tS t|tj	r|j
dkrtdd |D sd| d	}t|tS t|tj	r|j
dkrt| tjo| jdk}t| s|rd
d }|S d}t|nd| }t|d S )Nz-The first argument "array" must be array-likec                 S   s   t | |fS r  )r  broadcast_tor  r<   r<   r=   r     s    z numpy_broadcast_to.<locals>.implz7The second argument "shape" must be a tuple of integersr   c                 S   s   g | ]}t |tjqS r<   )r   r   IntegerLiteral)r   typr<   r<   r=   r     s     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   UniTuplerl   r  Tupler   allr   ra   r   )r   rW   r   r   Zis_scalar_arrayr<   r<   r=   numpy_broadcast_to  s4    


r  c                 C   sr   t t|D ]`}|t| | }|| }|dk r8td|dkrBq| | dkrX|| |< q| | |krtdqd S )Nr   #negative dimensions are not allowedr?   z=shape mismatch: objects cannot be broadcast to a single shape)r2   r   r   )rmrW   r_   r   tmpr<   r<   r=   numpy_broadcast_shapes_list  s    
r  c                     s   t | D ]n\}}t|tj}t|tjo4t|jtj}t|tjoNt|jdk}|s|s|sd| d| }t	|qd| D ]6}t|tjrt
dqt|tjrt
t|qdkrdd S d ddlm   fd	d
}|S d S )Nr   z	Argument z* must be either an int or tuple[int]. Got r?   c                  W   s   dS rK  r<   r   r<   r<   r=   <lambda>      z+ol_numpy_broadcast_shapes.<locals>.<lambda>r?   )_isinstance_no_warnc                     sz   dg }}t | D ]@} |tr>t|dkr>t|| q |trt||f qt|D ]\}}t|||}q`|S rR  )r   r   r   r  r   	enumerater*   )r   r  r  argr   elemr  r  tup_initr<   r=   r     s    

z'ol_numpy_broadcast_shapes.<locals>.impl)r  r   r   r   r  rl   r  r   r   r   maxr  numba.cpython.builtinsr  )r   r   r  Zis_intZis_int_tupleZis_empty_tupler   r   r<   r  r=   ol_numpy_broadcast_shapes  s*    r  c                     s   t | D ]"\}}t|std| dqd }d }| D ]H}t|tjtjfrV|j}n|}|d krh|}q8||kr8td|  q8d t | D ]T\}}t|tj	rt
 |j qt|tjtjtjfrt
 d qtd| qd   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|||| }q8g }t | D ]}|tt|| qZ|S rJ  )	r   r  r  r  rW   r2   r*   r[   r  )r   rW   r   r  r_   outsr  r  r<   r=   r   >  s    
z$numpy_broadcast_arrays.<locals>.impl)r  r   r   r   r   r   r   r  rl   ArrayCompatibler  ra   Numberrh  )r   r   r  Zunified_dtypedtr   r<   r  r=   numpy_broadcast_arrays  s.    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stt||D ]\}} | d||}q j|dd d}j t|f W 5 Q R X t |||\}}||
|j\}}t |||\}}t ||||}t ||||\}}nt	t
jrRj| }t|dkst tt!t
j"}| f} d|d |} j|dd d}j t|f W 5 Q R X  fdd}d	d
 }nfdd}dd
 }#t
j$d} |% \}!}"t&|"}"|j'D ]}#|"(|#|  qd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 sizer?   r   c                    s,   | \} tjttj}| |fS r  )rE   operatorgetitemr    r   r|   )r  r   getitem_implr7   rn   r  Z	src_dtyper  r<   r=   r    s    z#fancy_setslice.<locals>.src_getitemc                   S   s   d S r  r<   r<   r<   r<   r=   r    s    z#fancy_setslice.<locals>.src_cleanupc                    s    S r  r<   )r  )r  r<   r=   r    s    c                   S   s   d S r  r<   r<   r<   r<   r=   r    s    c                 S   s   g | ]}|d k	r|qS r  r<   )r   rW  r<   r<   r=   r     s      z"fancy_setslice.<locals>.<listcomp>r   )0r   rj   r	   r   rW   r   r   rf  r0  r   r   r   rl   r2  r  rq  r   r   r   rr  r   r   r)  r*  r   r  rw  r   r  r  r  SequencerE   r    r|   r   r}   ry  r   rg  rl  r{   r}  r~  r  ry   r{  r,  )'rn   r7   rO   r   r   r   r   r  r   Zdest_shapesdest_strides	dest_datarj  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_indicesrx  r_   r  rx   r  r<   r  r=   r!  U  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                    s"   g | ]\}}  ||tjqS r<   r{   r   r|   )r   r  rx   r   r<   r=   r     s   z#vararg_to_tuple.<locals>.<listcomp>rl   r   )r   r   r	   r\   r3   r   r  r|   r   r   r    r   )rn   r7   rO   r   r   Zdimtysr   r^   rW   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   rn   r7   rO   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   rW   r   r8  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 ]\}} || qtj|	dddtj|
} fdd|D }t||D ]6\}}t| ||
 tj|g|g|d d qވ tttjd	d |D  t|j }t||j |d  |d
 |j|j|jd | }t |j|S )Nr   r?   c                    s   g | ]}t  qS r<   r	   rK   r   r  )r7   ll_ary_sizer<   r=   r     s     z)array_transpose_tuple.<locals>.<listcomp>r   rl   ra   r~  c                    s   g | ]}t  qS r<   rj   r  )r7   rn   	np_ary_tyr<   r=   r      s     r   rW   r   r   r   c                 S   s   g | ]}|  qS r<   r   r   ar<   r<   r=   r   0  s        r  ) r   rj   r   rl   r   r   r|   r   	ArrayTyperW   r   r   r   r   r   rm   r   r  
as_pointercompile_internalr  r   r    voidr   r   r8   r   r   r   r   r   )rn   r7   rO   r   r   r   Zaxistyr   Znum_axisrl   r9  arysZll_arysr  dstZnp_itemsizeZnp_arysZnp_aryZll_aryr  r   r<   )r7   rn   r  r
  r=   array_transpose_tuple  sT    
	  r  c                 C   s    t | |||\}}t| |||S r  )r  r  rn   r7   rO   r   r  r  r<   r<   r=   array_transpose_vararg?  s    r  c                 C   s8   t | tjrtd|d kr*ddd}n
ddd}|S )Nz#np.transpose does not accept tuplesc                 S   s   |   S r  	transposer  axesr<   r<   r=   np_transpose_implK  s    z*numpy_transpose.<locals>.np_transpose_implc                 S   s
   |  |S r  r  r  r<   r<   r=   r  N  s    )N)N)r   r   r  r   UnsupportedError)r  r  r  r<   r<   r=   numpy_transposeE  s    

r  Tc           	      C   s   |j dkr|}nt|| ||}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  )ra   rj   r	   r   rW   r   r   r   r\   r   r   r   r   r   )	rn   r7   r  r   r   r   r  r   r   r<   r<   r=   r   T  s     
r   2   c                 C   sT   t | tjstdt |tjs,tdt |ttjfsFt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  re  numr  r<   r<   r=   r   q  s    znumpy_logspace.<locals>.impl)r"  )r   r   r  r   r   r   r   )r  re  r%  r   r<   r<   r=   numpy_logspaceh  s    



r&  r?   c                 C   sR   t |ttjfstdt | tjs0td| jdk rDtdddd}|S )	Nz*The second argument "k" must be an integer)The first argument "arr" must be an arrayr  Input must be >= 2-d.r?   c                 S   st   |d }|dkr| d d  S |dkr8t t | ddS |dkrPt t | S |dkrlt t | ddS td S )N   r   r?   r     )r  swapaxesfliplrflipudr   )r  r   r<   r<   r=   r     s    znumpy_rot90.<locals>.impl)r?   )	r   r   r   r   r   r   r   ra   NumbaValueError)r  r   r   r<   r<   r=   numpy_rot90x  s    




r/  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_reshaperW   r   r   )r   r   r|   r  intcr   FunctionTyper	   get_or_insert_functionr4   ra   rZ   rY   callr   )rn   r7   r   r   newndnewshape
newstridesr9  Zll_intp_starZll_intcfntyfnndrW   r   Z
is_f_orderr   r<   r<   r=   _attempt_nocopy_reshape  sD    	     

    r:  c           	      C   s   d}d}t |D ]&\}}|dk r.|d7 }|}q||9 }q|dkrR| |krtdnP|dkr|dkrpd}| dk}n| | }| | dk}|std|||< ntdd S )Nr   r?   z)total size of new array must be unchangedzmultiple negative shape values)r  r   )	ZorigsizerW   Znum_neg_valueZ
known_sizer	  r   Zneg_axinferredokr<   r<   r=   normalize_reshape_value  s(    



r=  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 5 Q R X t|| |}t||
j"|#||#||
j$|
j%|
j&d | }t'| ||j|S )	Nr   r?   r   r  r  rp  zincompatible shape for arrayr  )(r   r   r   r   r|   r   r  r   rj   r	   rK   r   r   rl   r   rm   r   r  r  r   r  r=  r   r    r  r}   r   r:  rs  r3   r   r)  r*  rF   r   r8   r   r   r   r   )rn   r7   rO   r   r   rettyr  rW   r9  Zll_shaper   r5  Zshape_ary_ty	shape_aryZshape_itemsizer   r4  r6  r<  failr   r  r   r<   r<   r=   array_reshape  sd    

 
 rA  c                 C   s    t | |||\}}t| |||S r  )r  rA  r  r<   r<   r=   array_reshape_vararg  s    rB  c                 C   s   dd }|S )Nc                 S   s
   |  |S r  )reshape)r  rW   r<   r<   r=   np_reshape_impl  s    z#np_reshape.<locals>.np_reshape_implr<   )r  rW   rD  r<   r<   r=   
np_reshape  s    rE  c                 C   s\   t | stdt |s$tdt|r8ddd}n t|tjsNt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  ravelr  concatenate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  rG  rH  r<   r<   r=   r   8  s    )N)N)r   r   r   r   r   r   r   )r  rI  r   r   r<   r<   r=   	np_append#  s    



rM  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)rC  r   r   r<   r<   r=   
imp_nocopy@  s    zarray_ravel.<locals>.imp_nocopyc                 S   s   |   S )zCopy version)flattenrN  r<   r<   r=   imp_copyD  s    zarray_ravel.<locals>.imp_copyr   r   )r   r~  r  r   r   )rn   r7   rO   r   rO  rQ  impr   r<   r<   r=   array_ravel=  s    rS  c                 C   s   dd }|  ||||S )Nc                 S   s   |   S r  )rF  r  r<   r<   r=   np_ravel_implV  s    znp_ravel.<locals>.np_ravel_impl)r  )rn   r7   rO   r   rU  r<   r<   r=   np_ravelT  s    rV  zarray.flattenc                 C   s,   dd }|  ||||}t| ||j|}|S )Nc                 S   s   |   | jS r  )r  rC  r   rN  r<   r<   r=   rR  _  s    zarray_flatten.<locals>.imp)r  r   r   )rn   r7   rO   r   rR  r   r<   r<   r=   array_flatten\  s    rW  c                 C   sl   |d krt | n|}t | ||\}}}t |jD ]0}|| }	|| }
|| }tt|	|
|||< q6|S r  )r  
empty_likebroadcast_arraysndindexrW   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_implg  s    rb  c                 C   sH   t | jD ]6}| | }|| }|r4t||||< qt||||< q|S r  )r  rZ  rW   r[  r  )r  bZuse_minr   r   ra  Zval_br<   r<   r=   _np_clip_impl_noneu  s    rd  c                 C   sv  t | stdt|tjs0t |s0tdt|tjsNt |sNtdt|tjspt|spd}t||d kpt|tj}|d kpt|tj}|r|rddd}|S t|tj}t|tj}	|r|	rddd}
|
S |r|	s| rdd	d
}|S ddd}|S nf|s<|	r<|r,ddd}|S ddd}|S n6|rPddd}|S |rdddd}|S ddd}|S d S ) N#The argument "a" must be array-likez6The 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d S )Nz&array_clip: must set either max or minr   r  r\  r]  r   r<   r<   r=   
np_clip_nn  s    znp_clip.<locals>.np_clip_nnc                 S   sH   |d krt | n|}t | jD ] }| | }tt|||||< q"|S r  )r  rX  rZ  rW   r[  r  r  r\  r]  r   r  r   ra  r<   r<   r=   
np_clip_ss  s
    znp_clip.<locals>.np_clip_ssc                 S   sB   |d krt | n|}t | jD ]}| | }t||||< q"|S r  )r  rX  rZ  rW   r  rh  r<   r<   r=   
np_clip_sn  s
    znp_clip.<locals>.np_clip_snc                 S   s   t | |}t| |||S r  r  	full_likerb  )r  r\  r]  r   Z
a_min_fullr<   r<   r=   
np_clip_sa  s    znp_clip.<locals>.np_clip_sac                 S   sB   |d krt | n|}t | jD ]}| | }t||||< q"|S r  )r  rX  rZ  rW   r[  rh  r<   r<   r=   
np_clip_ns  s
    znp_clip.<locals>.np_clip_nsc                 S   s   t | |}t| |||S r  rk  )r  r\  r]  r   Z
a_max_fullr<   r<   r=   
np_clip_as  s    znp_clip.<locals>.np_clip_asc                 S   s4   |d krt | n|}t | |\}}t||d|S )NTr  rX  rY  rd  )r  r\  r]  r   r  r^  r`  r<   r<   r=   
np_clip_na  s    znp_clip.<locals>.np_clip_nac                 S   s4   |d krt | n|}t | |\}}t||d|S )NFrp  )r  r\  r]  r   r  r^  r_  r<   r<   r=   
np_clip_an  s    znp_clip.<locals>.np_clip_anc                 S   s   t | |||S r  )rb  rf  r<   r<   r=   
np_clip_aa  s    znp_clip.<locals>.np_clip_aa)N)N)N)N)N)N)N)N)N)	r   r   r   r   r   NoneTyper   r   r  )r  r\  r]  r   r   Za_min_is_noneZa_max_is_nonerg  Za_min_is_scalarZa_max_is_scalarri  rj  rm  rn  ro  rq  rr  rs  r<   r<   r=   np_clip  sV    












ru  clipc                 C   s   ddd}|S )Nc                 S   s   t | |||S r  )r  rv  rf  r<   r<   r=   r     s    zarray_clip.<locals>.impl)NNNr<   )r  r\  r]  r   r   r<   r<   r=   
array_clip  s    
rw  c              	      s  |j |j kst|j|jks t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r2|d |kr}n|d |kr2}||krN| ksJ| dkrNdS |kr`| d }nd}||k r|| dkr|dS || }||  |9  < |||< n:||kr|| | }|| dkrdS || ||< |||< n dS )Nr!  r   Fr?   Tr<   )	r9  r^   r   old_itemsizenew_itemsizer~  r_   newdimZ
bytelengthZ
any_layoutZc_layoutZf_layoutint8r<   r=   rR  	  s.    



z_change_dtype.<locals>.imprW   r   r   r?   r  Nr  )ra   r   r~  ordr   r}  r   r|   rp   r	   rZ   rY   r   rj   r\   rm   rk   r   r   r    booleanr  r   r   r   )rn   r7   oldtynewtyr   Z
new_layoutrR  ry  rz  r9  Z
shape_dataZstrides_dataZshape_strides_array_typearyclsZshape_constantsizeof_intpZstrides_constantr?  Zstrides_aryrW   r   r   rO   r   r<   r|  r=   _change_dtype	  s    	
0   



	
r  c                 C   s   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  rW   rT  r<   r<   r=   r   	  s    znp_shape.<locals>.implr   r   r   r  r   r<   r<   r=   np_shape	  s    
r  c                 C   s   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   r_   r  rc  r<   r=   r   	  s      z5np_unique.<locals>.np_unique_impl.<locals>.<listcomp>)r  sortrF  r   r  r   )r  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rdt |tjr6tS t |tjtjfrZt |j	tjrZ|S t
dd S )Nc           	      S   s   t j|t jd}|jd dkr,t| |d S t |dk rBtdt | }| }|jd }|j|jkrrtdt jt 	||j
d}d}t|D ](}|| |||||  < ||| 7 }q|S )Nrl   r   r?   r  z(operands could not be broadcast together)r  r  int64rW   np_repeat_impl_repeats_scalerr  r   rF  emptyrm  rl   r2   )	r  repeatsZrepeats_arrayasaaraveln	to_returnposr_   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  rh  r   r  rl   r   r   )r  r  r  r<   r<   r=   	np_repeat	  s    r  c                 C   s   |dk rt dt| }| }|jd }|dkrDtjd|jdS |dkrVt|S tj|| |jd}t|D ] }|| ||| |d | < qr|S d S )Nr   r  r  r?   )	r   r  r  rF  rW   r  rl   r  r2   )r  r  r  r  r  r  r_   r<   r<   r=   r  	  s    


r  repeatc                 C   s   dd }|S )Nc                 S   s   t | |S r  )r  r  )r  r  r<   r<   r=   array_repeat_impl	  s    z'array_repeat.<locals>.array_repeat_implr<   )r  r  r  r<   r<   r=   array_repeat	  s    r  c                 C   s   t |}dd }||fS )z"Computes the itemsize of the dtypec                 S   s,   |  |jd j}|  |j}|| |S r   )rk   r   rl   r   rm   )cgctxr7   rO   llargsro   llintpr<   r<   r=   r  	  s    z%_intrin_get_itemsize.<locals>.codegen)r   r|   )tyctxrl   rO   r  r<   r<   r=   _intrin_get_itemsize	  s    
r  c                 C   s   d S r  r<   )r  rl   r<   r<   r=   _compatible_view	  s    r  generic)targetc                 C   s   dd }|S )zDDetermines if the array and dtype are compatible for forming a view.c                 S   s   t |}|| jkr| jdkr*d}t|n | jd }| j| dk}| jdk}| j| | jk}|rt|rt|rtd}t||| jk r|dks| j| dkrd}	t|	n(| j| | j }
|
| dkrd}t|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   ra   r   rW   r   r   )r  rl   Z
dtype_sizeZmsg1r   p1p2p3Zmsg2Zmsg3r{  Zmsg4r<   r<   r=   r   	  s(    






z ol_compatible_view.<locals>.implr<   r  rl   r   r<   r<   r=   ol_compatible_view	  s    #r  z
array.viewc              	   C   s,  |j d }|j}t|| ||d }t|| |}t|jj}t|D ]:}	t||	}
|	dkrt|jj	}|
|
||_qFt||	|
 qFtdkr| j}|t}|||j i }| ||}||| t| ||||}|d|t|j	d}|| d}| j|t|f W 5 Q R X | }t| ||j|S )Nr   r   )r?      rp  z"new type not compatible with array)r   r   rj   r   r   r   sortedgetattrr   r3   r  r   r   typing_contextresolve_value_typer  get_call_typerE   r  rs  r   r   r)  r*  r   r   r   )rn   r7   rO   r   r   r>  r   r  fieldsr   rx   r  r  r7  Z_compatible_view_sigr   r<  r@  r   r   r<   r<   r=   
array_view
  s0    



r  rl   c                 C   s   |   }t| |||S r  )r,  r   rn   r7   r  r   r   r<   r<   r=   array_dtypeD
  s    r  rW   c                 C   s(   t |}|| ||}|j}t| |||S r  )rj   rW   r   rn   r7   r  r   ru   r   r   r<   r<   r=   array_shapeJ
  s    r  r   c                 C   s(   t |}|| ||}|j}t| |||S r  )rj   r   r   r  r<   r<   r=   array_stridesS
  s    r  ra   c                 C   s   |  tj|j}t| |||S r  )r   r   r|   ra   r   r  r<   r<   r=   
array_ndim\
  s    r  r   c                 C   s(   t |}|| ||}|j}t| |||S r  )rj   r   r   r  r<   r<   r=   
array_sizec
  s    r  r   c                 C   s(   t |}|| ||}|j}t| |||S r  )rj   r   r   r  r<   r<   r=   array_itemsizek
  s    r  nbytesc                 C   s2   t |}|| ||}||j|j}t| |||S )z"
    nbytes = size * itemsize
    )rj   r   r   r   r   r  r<   r<   r=   array_nbytest
  s    r  
contiguousc                 C   s   |  tj|j}t| |||S r  )r   r   r  	is_contigr   r  r<   r<   r=   array_contiguous
  s    r  c_contiguousc                 C   s   |  tj|j}t| |||S r  )r   r   r  is_c_contigr   r  r<   r<   r=   array_c_contiguous
  s    r  f_contiguousc                 C   s   |  tj|j}t| |||S r  )r   r   r  is_f_contigr   r  r<   r<   r=   array_f_contiguous
  s    r  r  c                 C   s    |  tj|j }t| |||S r  )r   r   r  mutabler   r  r<   r<   r=   array_readonly
  s    r  ctypesc           	      C   sP   t |}|| ||}t|}| ||}|j|_|j|_| }t| |||S r  )rj   r   ArrayCTypesr   r   r   r   r   )	rn   r7   r  r   ru   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   r|   r   )rn   r7   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   )rn   r7   fromtytotyrx   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   )rj   r   r  r|   ra   rm   r   rl   r    bool_rW   r   r   r  )
checkerrn   r7   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   )rn   r7   r  r   flagsobjr   r<   r<   r=   array_flags
  s
    r  c                 C   sZ   |j jdkr2| j|||d}tt| ||j |j}n|j jdk}| tj|}t	| |||S )Nr   r   )
rU   r~  r   r  r   r   r   r   r  r   )rn   r7   r  r   r  r   rx   r<   r<   r=   array_flags_c_contiguous
  s     r  c                 C   sr   |j jdkr2| j|||d}tt| ||j |j}n2|j j}|j jdkrN|dkn|dk}| tj	|}t
| |||S )Nrx  r   r?   r  )rU   r~  r   r  r   r   ra   r   r   r  r   )rn   r7   r  r   r  r   r~  rx   r<   r<   r=   array_flags_f_contiguous
  s     r  realc                 C   sL   |j tjkrt| |||ddS |j tjkr8t| |||S tdtj d S )Nr  attrzunsupported .real for {})	rl   r   complex_domainarray_complex_attrnumber_domainr   rF   r   r3   )rn   r7   r  r   r<   r<   r=   array_real_part
  s
    r  imagc                 C   s   |j tjkrt| |||ddS |j tjkrt|jdd|}t| |||g\}}t| |||}t	
||j||j|jd t| ||j| S tdtj d S )Nr  r  Tr  r   zunsupported .imag for {})rl   r   r  r  r  r    r  _parse_empty_like_argsr|  r	   memsetr   r   r   r   r   r   r   rF   r   r3   )rn   r7   r  r   rO   r  r   r   r<   r<   r=   array_imag_part  s    
r  c                 C   s   |dks|j tjkr"td|t|}|| ||}|j j}| | |}|j	
|}	| | }
||j|
}|dkr||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  rl   r~  )r   r   )rl   r   r  rF   r   rj   underlying_floatrm   rk   r   r3   r   r  r  r   r  r   IntTyper  r   r	   copy_structr   r   )rn   r7   r  r   r  ru   r   ZfltyZsizeof_fltyr   Z
llfltptrtyr   Zresulttyr   replr<   r<   r=   r    s     r  conj	conjugatec                 C   s   dd }|S )Nc                 S   s
   t | S r  )r  r  r  r<   r<   r=   r   Q  s    zarray_conj.<locals>.implr<   r  r   r<   r<   r=   
array_conjN  s    r  c                 C   s   |   S r  )r,  )rn   r7   ZdtypetyZdtypevalr<   r<   r=   
dtype_typeY  s    r  r3   kindstatic_getitemc                 C   s@   |j }t|tjr.| |d}t| |||S d}t|dS )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   LoweringError)rn   r7   rO   r   r>  r   r   r<   r<   r=   static_getitem_number_clazzf  s    r  c              	      s  t |}| ||}|j}t|tjs6td||f ||}||}	t|tjrr|j	|j|j
|j
 dd}
n|j	|dd}
t |
}| |} tj|	}tj||j||jjd}t|tjr.t||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                    s   g | ]}  tj|qS r<   r   r   r|   rk  rn   r<   r=   r     s     z(array_record_getattr.<locals>.<listcomp>c                    s   g | ]}  tj|qS r<   r  rk  r  r<   r=   r     s     r  )rj   rl   r   r   RecordrF   typeofr   NestedArrayr  ra   r   r|   r	   r   r   r3   r   rW   r   rm   rk   r   r   r   r   r   )rn   r7   r  r   r  ru   r   Zrectyperl   r   restyZrarytyZraryZconstoffsetZ
newdataptrrW   r   Zdatasizer   r<   r  r=   array_record_getattr~  sX    

 
 
   r  c                 C   s0   |d }t |tstt| ||jd |d |S rR  )r   strrF   r  r   )rn   r7   rO   r   r   r<   r<   r=   array_record_getitem  s    
r  c              	      s    || ||}||}t|tjrt|}| |}|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rdnd} ||||}t |||S dS )zT
    Generic getattr() implementation for records: get the given record member.
    c                    s   g | ]}  tj|qS r<   r  r   r   r  r<   r=   r     s     z"record_getattr.<locals>.<listcomp>c                    s   g | ]}  tj|qS r<   r  r  r  r<   r=   r     s     Nr  r?   )sentry_record_alignmentr   r  r   r   r  rj   rl   rW   r   r	   get_record_memberrk   r   r\   r   r|   r   r   r   rs   rt   )rn   r7   r  r   r  r   elemtyr   r   rl   r5  r6  newdatar   dptrrr   r<   r  r=   record_getattr  sB    







	
r  c              	   C   s   |j \}}|\}}| || ||}	||}
t|
tjrt|| ||d d}|j	}t
|||	|jj}t|||| tj|
j nHt
|||	| |
}| ||||
}|jrdnd}| j||
|||d dS )zT
    Generic setattr() implementation for records: set the given record member.
    r?   r   Nrq   )r   r  r   r  r   r   r  r	   ri   r   r  r3   pointeememcpyr   r|   r   rk   r{   rs   rw   )rn   r7   rO   r   r  r  r"  r  rx   r   r   Z
val_structr  r  r  rr   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   )rn   r7   rO   r   r   r<   r<   r=   record_static_getitem_str  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  insert_const_stringr4   r  )rn   r7   rO   r   r   r  Zll_fieldr   r<   r<   r=   record_static_getitem_int  s
    r  static_setitemc                 C   sJ   |j \}}}|\}}}	t|j||}
| ||
}|dk	s<t||||	fS )3
    Record.__setitem__ redirects to setattr()
    N)r   r    r   get_setattrr   )rn   r7   rO   r   rectyr  r"  recr   rx   getattr_sigr   r<   r<   r=   record_static_setitem_str  s    
r  c                 C   s^   |j \}}}|\}}}	t|j||}
t|j d j}| || |
}|dk	sPt||||	fS )r  r   N)r   r    r   r   r  r  r   )rn   r7   rO   r   r  r  r"  r  r   rx   r  r  r   r<   r<   r=   record_static_setitem_int'  s    
r  c                 C   s   |  |||S )zB
    Create a constant array (mechanism is target-dependent).
    )make_constant_array)rn   r7   r  pyvalr<   r<   r=   constant_array9  s    r  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  	bytearraytostringr	   r   )rn   r7   r  r  ltyrx   r<   r<   r=   constant_recordA  s    r  c                 C   s"   t jt|t jd}| |||S )zM
    Create a constant array from bytes (mechanism is target-dependent).
    r  )r  r   r  uint8r  )rn   r7   r  r  bufr<   r<   r=   constant_bytesK  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  )rW   r   r  r   r  rc  r<   r<   r=   array_is_impl]  s
    
zarray_is.<locals>.array_is_impl)r   r	   rq  r  )rn   r7   rO   r   atybtyr   r<   r<   r=   array_isW  s
    
r#  __hash__c                 C   s   dd S )Nc                 S   s   d S r  r<   r  r<   r<   r=   r  j  r  zol_array_hash.<locals>.<lambda>r<   r  r<   r<   r=   ol_array_hashh  s    r%  c                 C   s
   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_clsp  s    r*  c                 C   s
   t | dS )zv
    Return the Structure representation of the given *nditerty* (an
    instance of types.NumpyNdEnumerateType).
    ndenumerater'  nditertyr<   r<   r=   make_array_ndenumerate_clsx  s    r.  c              
   C   s   |  tjd}|d}	|d k	r.|tj| tt|D ]}
t	|||
}t
|||}||
 }|d||}t||, ||| |d k	r||
 ||	 W 5 Q R X ||| |d k	r:||
 q:|d k	r|tj| ||	 ||	 d S )Nr   Zend_incrementr   )r   r   r|   r;  r   r	   
false_byterz  r2   rZ   r   r8   r   	if_likelyrC  	true_byterD  )rn   r7   ra   rW   r   end_flagloop_continue
loop_breakr   bbendr8  idxptrr   r   	in_boundsr<   r<   r=   _increment_indices  s*    


r8  c                 C   s*   t ||j|j}t| ||j||| d S r  )r	   r   rW   ra   r8  )rn   r7   r  r  r   r2  rW   r<   r<   r=   _increment_indices_array  s    r9  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_dimra   )rM   r:  r;  r<  r=  r<   r<   r=   r:    s
    z-make_nditer_cls.<locals>.BaseSubIter.__init__c                 S   s   t | j| j| d S r  )r   r:  r;  )rM   r]   r<   r<   r=   set_member_ptr  s    z3make_nditer_cls.<locals>.BaseSubIter.set_member_ptrc                 S   s   t | j| jS r  )r  r:  r;  r/  r<   r<   r=   
member_ptr  s    z/make_nditer_cls.<locals>.BaseSubIter.member_ptrc                 S   s   d S r  r<   rM   rn   r7   r<   r<   r=   init_specific  s    z2make_nditer_cls.<locals>.BaseSubIter.init_specificc                 S   s   d S r  r<   rM   rn   r7   logical_dimr<   r<   r=   r3    s    z2make_nditer_cls.<locals>.BaseSubIter.loop_continuec                 S   s   d S r  r<   rB  r<   r<   r=   r4    s    z/make_nditer_cls.<locals>.BaseSubIter.loop_breakN)rb   rc   rd   r6  r:  r>  r  cached_propertyr?  rA  r3  r4  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   r|   r>  r	   r   )rM   rn   r7   r   r<   r<   r=   rA    s    z2make_nditer_cls.<locals>.FlatSubIter.init_specificc                 S   s   | | j}||j|gS r  )r8   r?  r  r   )rM   rn   r7   r   r  r  r   r<   r<   r=   compute_pointer  s    z4make_nditer_cls.<locals>.FlatSubIter.compute_pointerc                 S   s8   || j d kr4|| j}t||}||| j d S rJ  )ra   r8   r?  r	   r   r   )rM   rn   r7   rC  r   r<   r<   r=   r3    s    z2make_nditer_cls.<locals>.FlatSubIter.loop_continuec                 S   s^   |dkr&| tjd}||| j n4|| jd krZ|| j}t||}||| j d S Nr   r?   )	r   r   r|   r   r?  ra   r8   r	   r   )rM   rn   r7   rC  r   r   r<   r<   r=   r4    s    z/make_nditer_cls.<locals>.FlatSubIter.loop_breakN)rb   rc   rd   r6  rA  rF  r3  r4  r<   r<   r<   r=   FlatSubIter  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
td S r  )need_shaped_indexingr   r@  r,  r<   r=   rA    s    z9make_nditer_cls.<locals>.TrivialFlatSubIter.init_specificc                 S   s&   t |dkstt |||j|S rJ  )r   r   r  r   rM   rn   r7   r   r  r  r<   r<   r=   rF    s    z;make_nditer_cls.<locals>.TrivialFlatSubIter.compute_pointerN)rb   rc   rd   r6  rA  rF  r<   r,  r<   r=   TrivialFlatSubIter  s   rK  c                   @   s   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ttj|||||ddS NFr  )r   ra   r   r	   r  rJ  r<   r<   r=   rF    s
     z7make_nditer_cls.<locals>.IndexedSubIter.compute_pointerNrb   rc   rd   r6  rF  r<   r<   r<   r=   IndexedSubIter  s   rN  c                   @   s   e Zd ZdZdd ZdS )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  	  s    z7make_nditer_cls.<locals>.ZeroDimSubIter.compute_pointerNrM  r<   r<   r<   r=   ZeroDimSubIter  s   rO  c                   @   s   e Zd ZdZdd ZdS )z&make_nditer_cls.<locals>.ScalarSubIterz8
        Sub-iterator "walking" a scalar value.
        c                 S   s   |S r  r<   rJ  r<   r<   r=   rF    s    z6make_nditer_cls.<locals>.ScalarSubIter.compute_pointerNrM  r<   r<   r<   r=   ScalarSubIter  s   rP  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 ]8\}}|\}}}}d| }	|| }
||
| |	|| q&|S )N)r&  indexed0dscalarzindex%d)rI  r  ri  r[   )rM   lZ	factoriesr_   rt  r  r<  r=  r  r;  factory)rH  rN  rP  rK  rO  r-  r<   r=   subiters  s    
z(make_nditer_cls.<locals>.NdIter.subitersc                    s@  | tjd}||t||| _t|D ]8\}}t|tjs,d| }t	
||| }	t| ||	 q,| ||||}ttj}
d}d}t|D ]6\}}t|tjr|jkr|| j}|| j} qqdkst|||
d}| tjd}dd }t||D ]T\}}t|tjr|jdkrttjttj|j|
}|||||j|f qt	||} dkr|ddd	 }|d
||}||t	jt	j}js|f}t|kst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   r2   r   )rW   
main_shaper  r_   r<   r<   r=   check_shapeO  s    zBmake_nditer_cls.<locals>.NdIter.init_specific.<locals>.check_shaperx  r!  rp  r?  )'r   r   r|   
make_tupler  arraysr  r   r   r	   r   r   _arrays_or_scalarsr  ra   rW   r   r   r   r    noner  r   r   rd  r1  r/  rI  r   rK   r3   r2   rZ   r   r   r\   	exhaustedrV  rA  )rM   rn   r7   arrtysr[  r   r_   r  r;  slotZmain_shape_tyrX  Zmain_nitemsr  rY  r  rO   r   Zshape_is_emptyr^  r   r8  r6  subiter)r~  ra   r-  nshapesr<   r=   rA  ,  sf    

 


z-make_nditer_cls.<locals>.NdIter.init_specificc                    s*  | d}t||| j}t|| |d || W 5 Q R X  j}t	|| j}| 
||||}| j}|d | |||||}	dd |	D }	t|	dkr||	d  n||| j|	 t	|| j}
t||t|
|
|| jt| j||t| j|| || || dS )	zJ
            Compute next iteration of the nditer() instance.
            r  FTc                 S   s   g | ]}|  qS r<   r  r   r   r<   r<   r=   r     s     zEmake_nditer_cls.<locals>.NdIter.iternext_specific.<locals>.<listcomp>r?   r   N)r;  r	   as_bool_bitr8   r^  if_unlikelyr   rC  r[  r   r\  r   _make_viewsr   r   rZ  r   rW   r8  r  partial_loop_continue_loop_breakrD  )rM   rn   r7   r   r5  r^  r_  r[  r   viewsrW   r,  r<   r=   iternext_specific{  sD    



 

z1make_nditer_cls.<locals>.NdIter.iternext_specificc                 S   s@   | j D ]4}|j|  kr"|jk rn q|||||j  qd S r  )rV  r<  r=  r3  rM   rn   r7   r8  rt  r<   r<   r=   rh    s    
z.make_nditer_cls.<locals>.NdIter._loop_continuec                 S   s@   | j D ]4}|j|  kr"|jk rn q|||||j  qd S r  )rV  r<  r=  r4  rl  r<   r<   r=   ri    s    
z+make_nditer_cls.<locals>.NdIter._loop_breakc                    s   dg }j }| j}j}	t|	tjr2t|	}	n|	g}	 fddtD t||D ]z\}
}|
\}}}}|j	|j
 }dkr|ddd }|D ]:}|| dkst| | ||	| || || |||< qqZtdd |D st|S )z:
            Compute the views to be yielded.
            Nc              	      s    g | ]}  t |qS r<   r8   r	   rZ   rk  r7   r   r<   r=   r     s   z?make_nditer_cls.<locals>.NdIter._make_views.<locals>.<listcomp>rx  r!  c                 s   s   | ]
}|V  qd S r  r<   rc  r<   r<   r=   r    s     z>make_nditer_cls.<locals>.NdIter._make_views.<locals>.<genexpr>)ri  rV  r   r   r   r  r   r2   r   r<  r=  r   
_make_viewr  )rM   rn   r7   r   r_  r[  rj  ri  rV  rettysrt  ra  r  array_indicesZsub_indicesr_   )r~  narraysr-  rb  rn  r=   rf    s2    


  z+make_nditer_cls.<locals>.NdIter._make_viewsc                 S   s   t |tjr|jdkst||||||}||||}	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   ra   r   rF  rj   rp   rZ  r  r|   r   )rM   rn   r7   r   r>  r  r  ra  r]   viewr   rW   r   r<   r<   r=   ro    s    
z*make_nditer_cls.<locals>.NdIter._make_viewc           	      S   s`   g }t t||D ]H\}\}}t|tjrF||||||d q|t| d|  q|S )Nr   rW  )r  r   r   r   r   r[   rj   r  )	rM   rn   r7   r_  r[  rT  r_   r  r  r<   r<   r=   r\    s    
z2make_nditer_cls.<locals>.NdIter._arrays_or_scalarsN)rb   rc   rd   r6  r  rD  rV  rA  rk  rh  ri  rf  ro  r\  r<   
rH  rN  rP  rK  rO  r~  rr  ra   r-  rb  r<   r=   NdIter  s   O*ru  )ra   r~  r   r[  rI  objectr	   ri   )r-  rE  ru  r<   rt  r=   make_nditer_cls  s    
 
, Yrw  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 ]Z}t	|||}|
|| || }	|d|	|}
t||
 |
tj| W 5 Q R X q@|| _|| _t|||j| _d S Nr   r?  rp  )r   r   r|   r	   rK   r3   r   r/  r2   rZ   r   rs  re  r1  r   r^  r\   rW   )rM   rn   r7   r   r   r   r^  r8  r6  r^  dim_is_emptyra   r<   r=   rA    s"    
z3make_ndindex_cls.<locals>.NdIndexIter.init_specificc           
   	      s   | tjd} d}t  j}t | |	d  
| W 5 Q R X  fddtD }|D ]}t | qt|t ||j |	d t j}	t| |	jj  
|  | d S )Nr   r  Fc              	      s"   g | ]}  t j|qS r<   )r8   r	   rZ   r   r   r8  r7   rM   r<   r=   r     s
   
zKmake_ndindex_cls.<locals>.NdIndexIter.iternext_specific.<locals>.<listcomp>T)r   r   r|   r;  r	   rd  r8   r^  re  r   rC  r2   rA   r   r\   r3   r   rW   r8  r   rD  )
rM   rn   r7   r   r   r5  r^  r   r8   rW   rz  r|  r=   rk    s,    




 
z7make_ndindex_cls.<locals>.NdIndexIter.iternext_specificN)rb   rc   rd   r6  rA  rk  r<   rz  r<   r=   NdIndexIter  s   r}  )ra   r	   ri   )r-  r}  r<   rz  r=   make_ndindex_cls  s    7r~  c                    s\    dkst | j}|jdkr:G  fdddt| }|S G  fdddt| }|S d S )N)r&  r+  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rztj||j	| tj|j
d}t|j
D ]}t|||}||| qT|| _d S )Nr   r+  r?  )r   r   r|   r	   r   r   r   rv  rK   r3   ra   r2   rZ   r   r   )	rM   rn   r7   r  r  r   r   r8  r6  r  r<   r=   rA  ;  s     zD_make_flattening_iter_cls.<locals>.CContiguousFlatIter.init_specificc              	      s   |j }|j} | j} d||}	||	 t |	  |j	|g}
t
| ||
}dkrp|| nP| j fddt|D }t |}|t ||g t| || t |} || j W 5 Q R X d S )Nr   r&  c              	      s    g | ]}  t |qS r<   rm  r{  rn  r<   r=   r   d  s
   
z\_make_flattening_iter_cls.<locals>.CContiguousFlatIter.iternext_specific.<locals>.<listcomp>)ra   r   r8   r   r   r   r	   r0  r  r   rv   r   r   r2   r\   make_anonymous_structr9  r   r   )rM   rn   r7   r  r  r   ra   r   r   r   r]   r   idxvalsidxtupler  rn  r=   rk  T  s4    
 zH_make_flattening_iter_cls.<locals>.CContiguousFlatIter.iternext_specificc                 S   s   | |j|g}t||||S r  )r  r   rv   rM   rn   r7   r  r  r   r]   r<   r<   r=   r  r  s    z>_make_flattening_iter_cls.<locals>.CContiguousFlatIter.getitemc                 S   s$   | |j|g}t||||| d S r  )r  r   ry   rM   rn   r7   r  r  r   r   r]   r<   r<   r=   setitemv  s    z>_make_flattening_iter_cls.<locals>.CContiguousFlatIter.setitemN)rb   rc   rd   r6  rA  rk  r  r  r<   r  r<   r=   CContiguousFlatIter6  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 ]t}t||	|}t||
|}||| ||| || }|d||}t|| |tj| W 5 Q R X q||	| _|
| _|| _d S rx  )r   r   r|   r   ra   r	   r   rW   rK   r3   r   r/  r2   rZ   r   rs  re  r1  r   pointersr^  )rM   rn   r7   r  r  r   r   ra   r   r   r  r^  r8  r6  ptrptrr^  ry  r<   r<   r=   rA    sB        z9_make_flattening_iter_cls.<locals>.FlatIter.init_specificc              
      s"  |j }t |j|}t |j|}| j| j}	|tj	d}
 
d}t  | j}t | |d  | W 5 Q R X |d t |	|d } |}t| ||}dkr|| n: fddt|D }t |}|t ||g tt|D ]}t |}t  |}|| }|| } d	||}t |~  || t |	|} |}t ||} || t|d |D ] }t |	|} || q | W 5 Q R X  |
| q tj| j  |  | d S )
Nr   r  FTr?   r&  c              	      s    g | ]}  t |qS r<   rm  r{  rn  r<   r=   r     s
   
zQ_make_flattening_iter_cls.<locals>.FlatIter.iternext_specific.<locals>.<listcomp>r   )ra   r	   r   rW   r   r   r  r   r   r|   r;  rd  r8   r^  re  r   rC  rZ   rv   r   r2   r\   r  rz  r   r   r0  r   r   r1  rD  )rM   rn   r7   r  r  r   ra   r   r   r  r   r5  Zis_exhaustedZlast_ptrr]   r   r  r  r8  r6  r   r   rv  r7  r  Z	inner_dimr  rn  r=   rk    sl    
 





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   )ra   r	   r   rW   r   rz  r2   r[   uremudivreverser}  r   r~  )rM   rn   r7   r  r  r   ra   r   r   r   r8  r]   r<   r<   r=   _ptr_for_index  s      z:_make_flattening_iter_cls.<locals>.FlatIter._ptr_for_indexc                 S   s    |  |||||}t||||S r  )r  rv   r  r<   r<   r=   r    s    z3_make_flattening_iter_cls.<locals>.FlatIter.getitemc                 S   s&   |  |||||}t||||| d S r  )r  ry   r  r<   r<   r=   r     s    z3_make_flattening_iter_cls.<locals>.FlatIter.setitemN)	rb   rc   rd   r6  rA  rk  r  r  r  r<   r  r<   r=   FlatIter}  s   !Er  )r   rU   r~  r	   ri   )r)  r  rU   r  r  r<   r  r=   r(  0  s    
D r(  r&  c                 C   sh   t t|}|| |}||_| |}|| ||dd}|| ||| | }t| |t||S Nr   )rN   )	r*  r   NumpyFlatTyper   rj   rY   rA  r   r   )rn   r7   r  r  flatiterclsflatiterarrclsr   r<   r<   r=   make_array_flatiter  s    

r  c                 C   sZ   |j \}|\}t|}|| ||d}|j}| |}	|	| ||jd}
|| |||
| d S r  )r   r*  rU   rj   r   rk  )rn   r7   rO   r   r   r)  r  r  r  r  r  r<   r<   r=   iternext_numpy_flatiter  s    
r  c                 C   sj   |j d }|\}}t|}|| ||d}|j}| |}	|	| ||jd}
|| |||
|}t| ||j|S Nr   r   )r   r*  rU   rj   r   r  r   r   )rn   r7   rO   r   r)  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   r*  rU   rj   r   r  r,  )rn   r7   rO   r   r)  r  r   r   r  r  r  r  r<   r<   r=   iternext_numpy_getitem_any7  s    


r  c           	      C   sF   |j d }t|}|| ||d d}| |j}|| ||jd}|jS r  )r   r*  rj   rU   r   r   )	rn   r7   rO   r   r)  r  r  r  r  r<   r<   r=   iternext_numpy_getitem_flatH  s    
r  c           
      C   sr   |j \}|\}tt|}|| |}||_| |}|| ||dd}|| ||| | }	t	| ||j
|	S r  )r   r.  r   NumpyNdEnumerateTyper   rj   rY   rA  r   r   r   )
rn   r7   rO   r   r  r  	nditerclsr:  r  r   r<   r<   r=   make_array_ndenumerateS  s    

r  c                 C   sZ   |j \}|\}t|}|| ||d}|j}| |}	|	| ||jd}
|| |||
| d S r  )r   r.  rU   rj   r   rk  )rn   r7   rO   r   r   r-  r:  r  r  r  r  r<   r<   r=   iternext_numpy_nditere  s    
r  c                    s^    fddt |j|D }ttt|}| }| | | }t |j	|S )zndindex(*shape)c                    s"   g | ]\}}  ||tjqS r<   r  )r   argtyr  r   r<   r=   r   y  s   z&make_array_ndindex.<locals>.<listcomp>)
r   r   r~  r   NumpyNdIndexTyper   rA  r   r   r   )rn   r7   rO   r   rW   r  r:  r   r<   r   r=   make_array_ndindexu  s    

r  c           
         s   |j j}|dkrJ|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   r7   rn   r~   r<   r=   r     s   z,make_array_ndindex_tuple.<locals>.<listcomp>)r   ra   r   rl   r	   r   r~  r   r  r   rA  r   r   )
rn   r7   rO   r   ra   r  rW   r  r:  r   r<   r  r=   make_array_ndindex_tuple  s    
r  c                 C   s6   |j \}|\}t|}|| ||d}|| || d S r  )r   r~  rk  )rn   r7   rO   r   r   r-  r:  r  r<   r<   r=   iternext_numpy_ndindex  s
    r  c           	      C   sn   |j }|j}t|jd tjr0t||d }n
|d g}t|| |}|	| ||| |
 }t| |||S )z
    nditer(...)
    r   )r   r[  r   r   r   r  r	   r   rw  rA  r   r   )	rn   r7   rO   r   r-  r_  r[  r:  r   r<   r<   r=   make_array_nditer  s    
r  c                 C   s2   |j \}|\}t|| ||d}|| || d S r  )r   rw  rk  )rn   r7   rO   r   r   r-  r:  r<   r<   r=   iternext_numpy_nditer2  s    r  c              	   C   sP  t |}|| |}| |j}| tjt| |}| tjd}tt	dd}	|D ]0}
|
||
}||d}||	||d}	qT|jdkrd}n|jdkr|g}t|dd D ]}|||d | qtt|}nT|jdkr$|g}|dd D ]}|||d | qt|}ntd|j|
||}||d}||	||d}	|j|	d	d
 | j|td W 5 Q R X |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!  rx  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  )+rj   rk   rl   r   r   r|   rp   r   r   r  smul_with_overflowr$  rr  ra   r~  rz  r[   r   r   rF   r   r   r)  r*  r   get_preferred_array_alignmentuint32r,  MemInfoPointervoidptrTypeRefr3   r    r  _call_allocatorr   meminfo_datar   r	   r\   r   r  r  )rn   r7   r  r   r  r   datatyper   Zarrlenoverflowr   Zarrlen_multr   Zdimension_sizeZallocsize_multr  rl   	align_valrr   r   mipZarytypeclassargtypesr   r   r   Zshape_arrayZstrides_arrayr<   r<   r=   r|    sv    
 



 
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  rr   r<   r<   r=   r     s    z _ol_array_allocate.<locals>.implr<   )r  r  rr   r   r<   r<   r=   _ol_array_allocate  s    r  c                 C   s   |  ||S )z9Trampoline to call the intrinsic used for allocation
    )r  )r  r   rr   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)rn   r7   r    r   r  rr   r   r<   r<   r=   r  +  s    zintrin_alloc.<locals>.codegen)r   r  r  r    )r  r  rr   r  r  rO   r<   r<   r=   r  '  s    r  c              
   C   s   dd }t |tjr.d}| |||tjg}n&t |tjs>t|j}tj	|||d}g }|D ]}|
|| ||j| q\| |tjd}	t|D ]>}
|d||
 |	}t|| | j|td W 5 Q R X q|S )z;
    Parse the shape argument to an array constructor.
    c           
   	   S   s   |  tj}|j}t|}t|d|d > d }|j|k rJ|||}n^|j|kr|d||}	|j	|	dd | j
|td W 5 Q R X |j|kr|||}n|}|S )z0Cast src to intp only if value can be maintainedr?   >Fr'  )z#Cannot safely convert value to intp)r   r   r|   r@   r   r  r   sextr   r   r)  r*  r   trunc)
rn   r7   src_tr  r   Z
intp_widthZintp_irmaxvalr   Z	is_largerr<   r<   r=   safecast_intp9  s$    


 
z#_parse_shape.<locals>.safecast_intpr?   r   r   r   )znegative dimensions not allowed)r   r   r   r{   r|   r  r   r   r	   r   r[   r3   get_constant_genericr2   r   re  r)  r*  r   )rn   r7   r  rx   r  ra   Zpassed_shapesr   r   r   r8  Zis_negr<   r<   r=   _parse_shape5  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  )rn   r7   rO   r   ZarrshapetypeZarrshaper  r<   r<   r=   _parse_empty_argse  s    
r  c                 C   sX   |j d }t|tjrJt|| ||d d}tj||j|jd}|j	|fS |j	dfS dS )z_
    Parse the arguments of a np.empty_like(), np.zeros_like() or
    np.ones_like() call.
    r   r   r   r<   N)
r   r   r   r   rj   r	   r   rW   ra   r   )rn   r7   rO   r   arytyper   r   r<   r<   r=   r  o  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   UnicodeTyper   r   )fnamerl   r   r<   r<   r=   _check_const_str_dtype}  s    r  c                 C   s"   |j }||||}dd }||fS )Nc                 S   s(   t | |||\}}t| |||}| S r  )r  r|  r   r  r7   rO   r  r  r   r   r<   r<   r=   r    s    znumpy_empty_nd.<locals>.codegeninstance_type)r  Zty_shapety_dtypety_retty_refty_rettyrO   r  r<   r<   r=   numpy_empty_nd  s    r  c                    s   t d| |tks0t|tjr(|jtks0t|r8tj}nt|}t	| }|d k	r||d k	r|tj
||dd tf fdd	}|S d|  d| d}t|d S )	Nr  r   r  c                    s   t | | S r  )r  )rW   rl   r>  r<   r=   r     s    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   )rW   rl   nb_dtypera   r   r   r<   r  r=   ol_np_empty  s"    

r  c                 C   s"   |j }||||}dd }||fS )Nc                 S   s(   t | |||\}}t| |||}| S r  )r  r|  r   r  r<   r<   r=   r    s    z$numpy_empty_like_nd.<locals>.codegenr  )r  Zty_prototyper  r  r  rO   r  r<   r<   r=   numpy_empty_like_nd  s    r  c                    s   t d| t|st|}nt| tjr0| j}n| }|d k	r~t| tjrn| jdkrX| jnd}| j||dd qt|dd nd|  d| d	}t	
|d fd
d	}|S )NrX  r  r   F)rl   r~  r  r   z3Cannot parse input types to function np.empty_like(r  r  c                    s   t | | S r  )r  )r  rl   r  r<   r=   r     s    zol_np_empty_like.<locals>.impl)N)r  r   r  r   r   r   rl   r~  r  r   r   )r  rl   r  r~  r   r   r<   r  r=   ol_np_empty_like  s    


r  c                 C   s   t |}dd }||fS )Nc                 S   s<   t |jd | ||d }t||j||j|jd d S r   )rj   r   r	   r  r   r   r   r   )r  r7   rO   r  r   r<   r<   r=   r    s    z(_zero_fill_array_method.<locals>.codegen)r   r]  )r  rM   rO   r  r<   r<   r=   _zero_fill_array_method  s    
r  
_zero_fillc                 C   s   dd }|S )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     s    z ol_array_zero_fill.<locals>.implr<   )rM   r   r<   r<   r=   ol_array_zero_fill  s    r  c                 C   s   t d| tfdd}|S )Nzerosc                 S   s   t j| |d}|  |S Nr  )r  r  r  )rW   rl   r  r<   r<   r=   r     s    zol_np_zeros.<locals>.impl)r  r  rW   rl   r   r<   r<   r=   ol_np_zeros  s    
r  c                 C   s   t d| ddd}|S )N
zeros_likec                 S   s   t j| |d}|  |S r  )r  rX  r  )r  rl   r  r<   r<   r=   r     s    zol_np_zeros_like.<locals>.impl)Nr  r  r<   r<   r=   ol_np_zeros_like  s    

r  c                 C   s   t d| ddd}|S )N	ones_likec                 S   s2   t j| |d}|j}tt|D ]}d||< q |S Nr  r?   r  rX  r&  r2   r   )r  rl   r  arr_flatr   r<   r<   r=   r     s
    
zol_np_ones_like.<locals>.impl)Nr  r  r<   r<   r=   ol_np_ones_like  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  r&  r2   r   )rW   
fill_valuerl   r  r  r   r  r<   r=   r  
  s
    
zimpl_np_full.<locals>.full)N)r  r   r  )rW   r  rl   r  r<   r  r=   impl_np_full  s    

r  c                 C   s   t d| ddd}|S )Nrl  c                 S   s0   t | |}|j}tt|D ]}|||< q|S r  r  )r  r  rl   r  r  r   r<   r<   r=   rl    s
    
z$impl_np_full_like.<locals>.full_like)Nr  )r  r  rl   rl  r<   r<   r=   impl_np_full_like  s    

r  c                 C   s   t d| ddd}|S )Nonesc                 S   s2   t j| |d}|j}tt|D ]}d||< q |S r  r  )rW   rl   r  r  r   r<   r<   r=   r   '  s
    
zol_np_ones.<locals>.impl)Nr  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 rJ  )r  r  r2   )r  rl   r  r_   r  r<   r=   r  8  s    z"impl_np_identity.<locals>.identity)N)r  r   r  r   r  )r  rl   r  r<   r  r=   impl_np_identity0  s    

r  c                 C   s   d S r  r<   NMr<   r<   r=   _eye_none_handler@  s    r  c                 C   s"   t |tjrdd }ndd }|S )Nc                 S   s   | S r  r<   r  r<   r<   r=   r   G  s    z$_eye_none_handler_impl.<locals>.implc                 S   s   |S r  r<   r  r<   r<   r=   r   J  s    r   r   rt  )r  r  r   r<   r<   r=   _eye_none_handler_implD  s    
r  c                    sf   |d kst |tjr tt n.t |tjtjfrDtt	|d| n
t| d dtf fdd	}|S )Nrl   r   c                    s   t | |}t| |f }|dkrPt| || }t|D ]}d|||| f< q8n,t| | |}t|D ]}d||| |f< qf|S rG  )r  r  r  r[  r2   )r  r  r   rl   _Mr  dr_   r  r<   r=   r   Z  s    
znumpy_eye.<locals>.impl)
r   r   rt  r  rl   r  	DTypeSpecr  r   r  )r  r  r   rl   r   r<   r  r=   	numpy_eyeO  s    
r  c                 C   sD   t | stdt| tjr@| jdkr2tdddd}|S d S )Nz#The argument "v" must be array-like)r?   r  zInput must be 1- or 2-d.r   c                 S   s.  | j dkr| j}|d t| }t||f| j}|dkrbt|| D ]}| | |||| f< qFn&t|| D ]}| | ||| |f< qn|S | j\}}|dk r|| }|dkr|| }tt||d}t	|| j}|dkrt|D ]}| ||| f ||< qn$t|D ]}| || |f ||< q
|S d S rR  )
ra   rW   r  r  r  rl   r2   r  r[  r  )r   r   r   r  r  r_   rowscolsr<   r<   r=   	diag_implr  s.    


zimpl_np_diag.<locals>.diag_impl)r   )r   r   r   r   r   r   ra   NumbaTypeError)r   r   r  r<   r<   r=   impl_np_diagi  s    



r  takec                    s   t | tjr$t |tjr$dd }|S tdd | |fD rT|jdk  fdd}|S t | tjr~t |tjtjfr~dd }|S d S )Nc                 S   s.   || j d ks|| j  k r"td|  | S )Nr?   Index out of bounds)r   
IndexErrorrF  )r  r   r<   r<   r=   	take_impl  s    znumpy_take.<locals>.take_implc                 s   s   | ]}t |tjV  qd S r  r  r   r  r<   r<   r=   r    s     znumpy_take.<locals>.<genexpr>rx  c                    s   t j|j| jd} r | }n|}t |}d}|  }|D ]:}|| jd ks\|| j k rdtd|| ||< |d }q>||j	S Nr  r   r?   r  )
r  r  r   rl   r  r:  rF  r  rC  rW   )r  r   r  Zwalkeritr_   r&  r  F_orderr<   r=   r    s    


c                 S   s~   t |}t j|j| jd}t |}d}|  }|D ]:}|| jd ksT|| j k r\td|| ||< |d }q6||j	S r  )
r  r   r  r   rl   r:  rF  r  rC  rW   )r  r   convertr  r  r_   r&  r  r<   r<   r=   r    s    


)r   r   r   r   r  r~  r  r  )r  r   r  r<   r  r=   
numpy_take  s    
r
  c                  G   s|   dd | D }t dd |D r(tj}nPt dd |D rBtj}n6t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     s      z!_arange_dtype.<locals>.<listcomp>c                 s   s   | ]}t |tjV  qd S r  r   r   Complexr  r<   r<   r=   r    s     z _arange_dtype.<locals>.<genexpr>c                 s   s   | ]}t |tjV  qd S r  )r   r   Floatr  r<   r<   r=   r    s     zint%sr  c                 S   s   g | ]}t |qS r<   )r   	unliteralr   r  r<   r<   r=   r     s     )
r  r   
complex128r   r  r  rl   r   r   r  )r   boundsrl   ZNPY_TYZunliteral_boundsr<   r<   r=   _arange_dtype  s    r  c                    s  t |tjr|j}t |tjr$|j}t |tjr6|j}|d krDtj}|d krRtj}|d kr`tj}t | tjrt |tjtjfrt |tjtjfrt |tjtjfsd S t |tjrt| ||n|j	t
dd | ||fD t| dd  t|dd t|dd d fdd	}|S )Nc                 S   s   g | ]}t |tjqS r<   r  r  r<   r<   r=   r     s   znp_arange.<locals>.<listcomp>r   c                    s    d k	r n| }d k	rn|}d k	r,n|}|d k	r<|nd}|d krTd| }}	n
|| }}	|dkrnt d|	| | }
tt|
j}dkrtt|
j}tt||d}n
t|d}t	|}|}t
|D ]}|||  ||< q|S )Nr?   r   zMaximum allowed size exceededT)r   r   mathceilr  r  r  r[  r  r  r2   )r  re  r`  rl   Z	lit_startZlit_stopZlit_stepZ_step_start_stopnitems_cnitems_rnitems_ir   r  rx   r_   Zstart_valueZ
step_valueZ
stop_valueZ
true_dtypeZuse_complexr<   r=   r     s(    

znp_arange.<locals>.impl)NNN)r   r   Optionalr3   r]  r  rt  r  r  rl   r  r  )r  re  r`  rl   r   r<   r  r=   	np_arange  s>    r  c                    sn   t dd | |fD sd S t|ttjfs8d}t|tdd | |fD rVtj ntj	 d fdd	}|S )	Nc                 s   s   | ]}t |tjV  qd S r  )r   r   r  r  r<   r<   r=   r     s     z!numpy_linspace.<locals>.<genexpr>z%The argument "num" must be an integerc                 s   s   | ]}t |tjV  qd S r  r  r  r<   r<   r=   r  '  s     r"  c                    s   t | }| d } |d }|dkr(|S |d }|dkrn||  }t ||}td|D ]}| ||  ||< qVn| |d< |dkr||d< |S )Ng      ?r   r?   r!  )r  r  divider2   )r  re  r%  r  divdeltar`  r_   r  r<   r=   r#  -  s    z numpy_linspace.<locals>.linspace)r"  )
r  r   r   r   r   r   r   r  r  r   )r  re  r%  r   r#  r<   r  r=   numpy_linspace  s    
r   c              
   C   s  |j d }t|| ||d d}t||j}|j}t| |||}|j}	|j}
|jdks\t	|j|jkrtj
||
|	|j|jdd nt||j}t||j}| tj}t|||H}t| ||	|||j|}t| ||
|||j|}|||| W 5 Q R X t| ||j| S )z
    Array copy.
    r   r   r  r?   rq   )r   rj   r	   r   rW   r   r|  r   r~  r   
raw_memcpyr   r   r   r   r   r|   r  r}  r   r8   r   r   )rn   r7   rO   r   r  r   r   rettyper  r  r  r  r  r   r   r  r  r<   r<   r=   _array_copyF  s>    
 
  
  r#  c                 C   s@   t |tjst|jdkrdnd}|j|dd}||}|tfS )Nrx  r   F)r~  r  )r   r   r   r   r~  r  r#  )r  r  r~  r  rO   r<   r<   r=   _array_copy_intrinsick  s
    r$  z
array.copyc                 C   s   t | |||S r  )r#  r  r<   r<   r=   
array_copyt  s    r%  c                 C   s   t | tjrdd }|S )Nc                 S   s   t | S r  )r$  rT  r<   r<   r=   
numpy_copy|  s    z#impl_numpy_copy.<locals>.numpy_copyr  )r  r&  r<   r<   r=   impl_numpy_copyy  s    r'  c              
   C   s  |j }|jd }|j|ks"td|jdkr|jdks:t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s|jdkr|jdkrt| |||d S |jdkr|dkst|dkrtnt}t|| |||d }||V\}}| t| |||d }|j}W 5 Q R X | t| |||}|j}W 5 Q R X W 5 Q R X ||j}||| ||| |S t| |||S d	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   ra   rj   r  r   r  r|   rZ  r   r   r   r   r   r   r   r   r   r  r  blockr#  phir3   add_incoming)rn   r7   rO   r   output_layoutr>  r   r   r  rW   r   
check_funcr  thenorelseZout_thenZthen_blkZ
out_orelseZ
orelse_blkZret_phir<   r<   r=   _as_layout_array  sd    

   r/  c                    sH   t  tjst |j jt|jdd}|| }| fddfS )Nr?   )r~  ra   c                    s   t | ||| jdS )Nr+  )r/  r   )rW  rc  r   r  r0  r<   r=   r    s       z,_as_layout_array_intrinsic.<locals>.<lambda>)	r   r   r!   r   RequireLiteralValuer  r   r  ra   )r  r  r+  r  rO   r<   r0  r=   _as_layout_array_intrinsic  s
    

r2  c                 C   sF   t | stdt| tjtjfr.dd }nt| tjrBdd }|S )Nre  c                 S   s   t t | S r  )r  ascontiguousarrayr   rT  r<   r<   r=   r     s    z%array_ascontiguousarray.<locals>.implc                 S   s
   t | dS Nr   r2  rT  r<   r<   r=   r     s    r   r   r   r   r   r  rh  r   r  r<   r<   r=   array_ascontiguousarray  s    

r7  c                 C   sL   t | stdt| tjtjfr0dd }|S t| tjrHdd }|S d S )Nre  c                 S   s   t t | S r  )r  asfortranarrayr   rT  r<   r<   r=   r     s    z"array_asfortranarray.<locals>.implc                 S   s
   t | dS )Nrx  r5  rT  r<   r<   r=   r     s    r6  r  r<   r<   r=   array_asfortranarray  s    
r9  zarray.astypec              
   C   s  |j d }t|| ||d d}t||j}|j}t| |||}|j}	|j}
t||j}t||j}| 	t
j}t|||h}t| ||	|||j|}t| ||
|||j|}t| |||}| |||j|j}t| |||| W 5 Q R X t| ||j| S r  )r   rj   r	   r   rW   r   r|  r   r   r   r   r|   r  r}  r~  rv   r{   rl   ry   r   r   )rn   r7   rO   r   r  r   r   r"  r  r  r  r  r  r   r   r  r  itemr<   r<   r=   array_astype  s6    

  
  r;  c                 C   s"   |j }||||}dd }||fS )Nc              	   S   s  |j d }|j}t|| ||d d}t|}|| |}|j}	t| |}
t|jj|
}||j	|j}|
||}t||}|j|dd d}| j|t|f W 5 Q R X 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   rj   r   rp   r   r   r3   r   r   sremr	   is_not_nullr   r)  r*  r   r\   sdivr  r   r   r   r   r   r   r   r   )rn   r7   rO   r   Zbuftyr   r  Z
out_ary_tyout_aryZout_datamodelr   ll_itemsizer  remZis_incompatibler   rW   r   r   r   r<   r<   r=   r    s<    


 znp_frombuffer.<locals>.codegenr  )r  bufferrl   r>  r  rO   r  r<   r<   r=   np_frombuffer
  s    $rC  c                    s   t d| t| tjr | jdkr4d|  }t||tksZt|tjrR|j	tksZt
|rbtj}nt|}|d k	rt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?   )rl   ra   r~  r  z3Cannot parse input types to function np.frombuffer(r  r  c                    s   t | | S r  )rC  )rB  rl   r  r<   r=   r   M  s    z impl_np_frombuffer.<locals>.impl)r  r   r   r   r~  r   r   r  r  r  r   r  r  r   r  )rB  rl   r   r  r   r<   r  r=   impl_np_frombuffer6  s*    





rE  c                    sP   t |r$tdd  d fdd	}|S t|tjrLtd| d fdd	}|S d S )Nr   c                    s
    | |S r  r<   r]   rW   rl   intrinsic_cfarrayr<   r=   r   W  s    zimpl_carray.<locals>.implc                    s
    | |S r  r<   rF  rG  r<   r=   r   ]  s    )N)Nr   get_cfarray_intrinsicr   r   r  r]   rW   rl   r   r<   rG  r=   impl_carrayR  s    

rL  c                    sP   t |r$tdd  d fdd	}|S t|tjrLtd| d fdd	}|S d S )Nrx  c                    s
    | |S r  r<   rF  rG  r<   r=   r   g  s    zimpl_farray.<locals>.implc                    s
    | |S r  r<   rF  rG  r<   r=   r   m  s    )N)NrI  rK  r<   rG  r=   impl_farrayb  s    

rM  c                    s   t  fdd}|S )Nc           	         s   |t jkrd }n*t|t jr$|j}nd| d}t| d kr^|d krXd}t||}nVt t jr j}|d k	r||krd| d| d}t|nd  d}t|t|}|d krd| 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   CPointerrl   r   r  r  r  r   r    
np_cfarray)	r  r]   rW   	ptr_dtyper   rl   ra   r>  rO   dtype_r~  r<   r=   rH  s  s2    




z0get_cfarray_intrinsic.<locals>.intrinsic_cfarrayr%   )r~  rS  rH  r<   rR  r=   rJ  r  s    rJ  c                    s<  |j dd \}}|dd \}}|j}|jdks6tt| }	t|}
t|
}t|t	j
rrt |}n|f}|f} fddt||D }|}g }|jdkr|D ]}||  ||}qn,t|D ]}||  ||}q|   ||j }t|	||||dd |	 }t |j|S )zR
    numba.numpy_support.carray(...) and
    numba.numpy_support.farray(...).
    Nr  r  c                    s"   g | ]\}}  ||tjqS r<   r  )r   r  r   r   r<   r=   r     s   znp_cfarray.<locals>.<listcomp>rx  r  )r   r   r~  r   rj   rp   r	   r   r   r   r  r   r   r[   r   rz  r  r  rk   rl   r  r   r   r   )rn   r7   rO   r   r  r  r]   rW   r   r?  r   r@  r   offr   r   r   r   r<   r   r=   rP    sH    




	rP  c                 C   sV   t |tjr| tjt|S t |tjrJ| tttj|}|||fS dsRt	d S r   )
r   r   r  r   r|   r   r  rE   r    r   )rn   r7   seqtyseqr  r<   r<   r=   _get_seq_size  s    rX  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   decref)r7   r   r  rn   r  r>  r<   r=   wrap  s    
z$_get_borrowing_getitem.<locals>.wrap)rl   rE   r  r  r    r   r|   )rn   rV  r[  r<   rZ  r=   _get_borrowing_getitem  s    r\  c                    st    tj}t|d fdd}g }|| }}	t|D ].}
|
dkrV|||	\}}	|t ||	 q<t|S )zF
    Compute the likely shape of a nested sequence (possibly 0d).
    r   c                    sT   t | tjr2t| dkrdS | d  |dfS nt| }| j| |ffS d S )Nr   )NN)r   r   r  r   r$  r\  rl   )rV  rW  r  r7   rn   r   r<   r=   get_first_item  s    
z.compute_sequence_shape.<locals>.get_first_item)r   r   r|   r   r2   r[   rX  r   )rn   r7   ra   rV  rW  r   r^  r   innertyinnerr_   r<   r]  r=   compute_sequence_shape  s    

ra  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 5 Q R X t |dkrdd S t| tjrt| }t	|0}| j
}|||jf}	||	|dd   W 5 Q R X nXt| tjrtt | D ],}
| |
 }||
}	||	|dd   qndst| d S )Nr   r&  Fr'  r?   )r   rX  r   r   r   r   r  r\  r	   rT  rl   r   r  r2   r$  r   )rV  rW  r   r   expectedr  r  rV  r_  r`  r_   rb  r7   check_seq_sizern   r<   r=   re    s*    
 z,check_sequence_shape.<locals>.check_seq_sizeNr<   )rn   r7   rV  rW  r   r<   rd  r=   check_sequence_shape  s    rf  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 rL  )r	   r}  r~  r{   rl   ry   )r   r"  rx   r]   )r  r7   rn   r   r   r   r<   r=   assign_item3  s      z-assign_sequence_to_array.<locals>.assign_itemc              	      s  t |dkr2t| tjtjfr"t|| | d S |d }t| tjrt| }t|:}| j	}|||j
f} |||dd  ||j
f  W 5 Q R X nlt| tjrtt | D ]B}	| |	 }||	}tj|	}
 |||dd  ||
f  qnds
t| d S rG  )r   r   r   r  r  r   r\  r	   rT  rl   r   r2   r$  r   r|   )rV  rW  r   r   r   r  rV  r_  r`  r_   r   )assignrg  r7   rn   r<   r=   rh  9  s$    
* z(assign_sequence_to_array.<locals>.assignr<   Nr<   )rn   r7   r   r   r   r  rV  rW  r<   )r  rh  rg  r7   rn   r   r   r   r=   assign_sequence_to_array,  s    ri  c                 C   s>   t | |\}}t|r|}nt|}|d kr0d S t||dS r4  )r/   r   r  r   r   )r  rv  rl   ra   Z	seq_dtyper<   r<   r=   np_array_typerU  s    rj  c                 C   s0   t d| t| ||}|||}dd }||fS )Nr   c           
   	   S   s   |j }|j}|jd }|d }t| ||||}t||ks>tt| |||| t| |||}	t| ||	j	||	j
||| t| ||j |	 S r   )r   ra   r   ra  r   r   rf  r|  ri  r   r   r   r   )
rn   r7   rO   r   r  ra   rV  rW  r   r  r<   r<   r=   r  f  s"    
   
znp_array.<locals>.codegen)r  rj  )r  objrl   r  rO   r  r<   r<   r=   np_array`  s
    

rl  c                 C   sL   t d| t| stdt|s>t|d k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   s
   t | |S r  )rl  )rv  rl   r<   r<   r=   r     s    zimpl_np_array.<locals>.impl)N)r  r   r   r   r   r  )rv  rl   r   r   r<   r<   r=   impl_np_arrayz  s    



rm  c           
   	   C   s   | d}| |}|d||}||||||}||d|||d||}|j|dd  d| }	| j|t|	f W 5 Q R X |S )Nr   r   rB  Fr'  z%s(): axis out of bounds)	r3   r   rd  rL  rr  r   r)  r*  r  )
rn   r7   r   ra   r   r   Zll_ndimZis_neg_axisZaxis_out_of_boundsr   r<   r<   r=   _normalize_axis  s    

rn  c              
   C   s   t ||d kstttj|}t||}td}t|d D ]R}t|}	|d|	|}
|| }|	|
|
|	||	}||t||d| qD||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?   rB  r   r   r   r   r  r	   r   rK   r2   r   rd  rL  r   rZ   r   r8   )rn   r7   
orig_shapera   r   ll_shtyr   r  r8  ll_dim
after_axisr  r   r<   r<   r=   _insert_axis_in_shape  s    


rt  c              
   C   s   t ||d kstttj|}t||}td}td}t|d D ]N}	t|	}
|d|
|}|	||
|
||
}|||	 t||d| qN||t||d| t|||S )zD
    Same as _insert_axis_in_shape(), but with a strides array.
    r?   r   rB  ro  )rn   r7   orig_stridesra   r   rq  r   r  r   r8  rr  rs  r   r<   r<   r=   _insert_axis_in_strides  s"    




rv  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   ra   r   rj   r	   r   rW   r   rt  rv  r   r   r   r   r   r   )rn   r7   rO   r   r   r>  ra   r  r  r  r   r   
new_shapesZnew_stridesr<   r<   r=   expand_dims  s$    
rx  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)	r{   r   r   r|   rn  r   ra   rx  r   )rn   r7   rO   r   r   r  r<   r<   r=   r    s     znp_expand_dims.<locals>.codegen)ra   r~  r  )r  r  r   r~  r  rO   r  r<   r<   r=   np_expand_dims  s
    
ry  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   s
   t | |S r  )ry  )r  r   r<   r<   r=   r     s    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    



rz  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 )rz  )r  r  ra   r  )minimumr<   r=   r     s     z-_atleast_nd.<locals>.impl.<locals>.<listcomp>c                    sn   t t |d } fddt|D }t|dkrV |j|}n|d }t |j|S )Nr   c                    s"   g | ]\}}} |||qS r<   r<   )r   r  r  r>  r7   rn   	transformr<   r=   r     s   z>_atleast_nd.<locals>.impl.<locals>.codegen.<locals>.<listcomp>r?   )_atleast_nd_transformr	   r   r   r   rZ  r   r   )rn   r7   rO   r   arrsretsr  )r_  r  r{  rp  r|  r=   r    s    

z*_atleast_nd.<locals>.impl.<locals>.codegenr?   r   )r    r   r   r  StarArgTuple
from_types)typingcontextr   r  )r  r{  )r_  rp  r=   r     s    
z_atleast_nd.<locals>.implc                     s    |  S r  r<   r  )r   r<   r=   r  (  r  z_atleast_nd.<locals>.<lambda>rT  )r{  r  r<   )r  r   r{  r=   _atleast_nd  s    r  c                    s"   t  kst fdd}|S )z`
    Return a callback successively inserting 1-sized dimensions at the
    following axes.
    c           	         sb   t D ]T}|d }|j|k rt | }|j|jd d}t| |t|||f|}|}q|S )Nr?   rz  )r2   ra   r	   r   r  rx  r   r    )	rn   r7   r  r  r>  r_   ra   r   Znewarrtyr  min_ndimr<   r=   r}  2  s    

 z(_atleast_nd_transform.<locals>.transform)r   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   s   | ]}t |tjV  qd S r  r  r  r<   r<   r=   r  D  s     z np_atleast_1d.<locals>.<genexpr>r?   r   r  r  r  r<   r<   r=   np_atleast_1dB  s    r  c                  G   s$   t dd | D r tdddgS d S )Nc                 s   s   | ]}t |tjV  qd S r  r  r  r<   r<   r=   r  J  s     z np_atleast_2d.<locals>.<genexpr>r  r   r  r  r<   r<   r=   np_atleast_2dH  s    r  c                  G   s&   t dd | D r"tddddgS d S )Nc                 s   s   | ]}t |tjV  qd S r  r  r  r<   r<   r=   r  P  s     z np_atleast_3d.<locals>.<genexpr>r*  r   r  r  r  r<   r<   r=   np_atleast_3dN  s    r  c	                 C   s  t |t |  kr.t |  kr.t |ks4n ttd}	t| |||}
t||
j}g }|D ]d}|	}tt||D ]B\}\}}|	d|
||}|||}||||||}qx|| qb|
j}t|||||D ]\}}}}}|j}tj||tj|jd}|h}t| |||||j|}t| |||}| |||j|j}t| |||||j|}t| |||| W 5 Q R X t|||}q|
S )z2
    Concatenate arrays along the given axis.
    r   rp  )order)r   r   r	   r   r|  r   r   r  r   r   r3   r   rd  rL  r[   r   r  r~  r}  rv   r{   rl   ry   r   )rn   r7   r   r_  r  
arr_shapesarr_stridesr>  
ret_shapesr   r  Zret_stridesZcopy_offsetsZarr_shr   r8  r   rv  is_axisaddendZret_datar  r  Zarr_stZarr_datar  r   r  rx   r  r<   r<   r=   _do_concatenateT  sT    4

 
  
  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 } |
\}}|& t	 j
||g } || W 5 Q R X |d tj}|D ]} | d||}q j |dd j td f W 5 Q R X W 5 Q R X W 5 Q R X qr fdd|	D }	t |||||||		}t || S )Nc                    s"   g | ]\}}t | |d qS r   r	  r   r!  r  r   r<   r=   r     s   z#_np_concatenate.<locals>.<listcomp>np.concatenatec                    s   g | ]}t  |jqS r<   )r	   r   rW   r   r  r  r<   r=   r     s     c                    s   g | ]}t  |jqS r<   r	   r   r   r  r  r<   r=   r     s     c                    s   g | ]}t  |qS r<   )r	   r   r  r  r<   r=   r     s   r   rp  c                    s   g | ]}|  qS r<   r<   r  )r8  r<   r=   r     s     r?   Fr'  z<np.concatenate(): input sizes over dimension %d do not matchc                    s   g | ]}  |qS r<   )r8   r  r  r<   r=   r     s     )ra   r   rn  r2   r   r3   r8   r  r  r  rL  r   r	   r  r  r   rZ  r)  r*  r   r  r   r   )rn   r7   r_  r  r>  r   ra   r  r  r  r  Zret_shape_ptrZret_shZother_shapesZon_axisZon_other_dimr  is_okr  r<   )r7   rn   r8  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 ]n}tj}tt |j|
D ]N\}} 	| 
d||} j |dd j td	 W 5 Q R X qqp fd
d|D }ttj|t }t }t|d D ]l}t|} 
d||}|
| } | |||} |t |d|  |t |d| q$ |t |d|  |	t |d| t  |}|gt| }t  |} fddtt|D }t|d D ]l}t|} 
d||} | |||}tt|D ]*} || | t || d| qZqtt|D ]"} |t || d| q fdd|D }t |||||||	}t || S )Nr   r?   c                    s"   g | ]\}}t | |d qS r  r	  r  r   r<   r=   r     s   z_np_stack.<locals>.<listcomp>np.stackrp  Fr'  )z5np.stack(): all input arrays must have the same shapec                    s   g | ]}t  |jqS r<   r  r  r  r<   r=   r     s     rB  c                    s   g | ]}t  qS r<   r  rk  )r7   rq  r<   r=   r     s   c                    s   g | ]}t   |qS r<   )r	   r   r8   )r   r  r  r<   r=   r     s   )ra   r	   r   r   r   rn  r   rW   r  r  r   r   rZ  r)  r*  r   r   r  rK   r2   rd  rL  r   rZ   r8   r  r   r   )rn   r7   r_  r  r>  r   ra   r   r  Z
ll_narraysrp  r  r  r  Zorig_shru  input_shapesr  r8  rr  rs  r   Zinput_stridesr_   r  r<   )r7   rn   rq  r=   	_np_stack  s    

 






    r  c                 C   sN   |d k	rt |tjsd S t| d|\}}|dkr8tdt|}t|||S )Nr  r   z.zero-dimensional arrays cannot be concatenated)r   r   r   r0   	TypeErrorr1   r   )r  r[  r   rl   ra   r~  r<   r<   r=   np_concatenate_typer  s     r  c                 C   s6   t | ||}t|tjst|||}dd }||fS )Nc              	   S   sF   |  ||d |jd tj}t| |t|jd t||d |j|S rR  )	r{   r   r   r|   r  r   r	   r   r   rn   r7   rO   r   r   r<   r<   r=   r  3  s    znp_concatenate.<locals>.codegen)r  r   r   r   r   )r  r[  r   r  rO   r  r<   r<   r=   np_concatenate-  s
    
r  c                 C   s   t | tjrddd}|S d S )Nr   c                 S   s
   t | |S r  )r  r[  r   r<   r<   r=   r   A  s    z!impl_np_concatenate.<locals>.impl)r   r   r   r  r[  r   r   r<   r<   r=   impl_np_concatenate>  s    
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)ra   r  )rn   r   r[  r  r<   r<   r=   _column_stack_dimsF  s    
r  c                 C   sB   t | d|td\}}t|}t|||}||}dd }||fS )Nznp.column_stack)dim_chooserc                 S   s   t |jd }t||d }g }g }| tjd}t||D ]t\}	}
|	jdkrf|	|	 |	|
 q>|	jdkstt
|	jdd}t||	}t| |||
f|}|	| |	| q>t| ||||j|S )Nr   r?   r  rz  )r   r   r	   r   r   r   r|   r   ra   r[   r   r  r   r    rx  r  r   )rn   r7   rO   r   Zorig_arrtysZ	orig_arrsr_  r  r   r  r  r  
expand_sigZnewarrr<   r<   r=   r  X  s&    



 z np_column_stack.<locals>.codegen)r0   r  r1   r   r   )r  r  rl   ra   r~  r  rO   r  r<   r<   r=   np_column_stackO  s     
r  c                 C   s   t | tjrdd }|S d S )Nc                 S   s   t | S r  )r  r  r<   r<   r=   r   x  s    zimpl_column_stack.<locals>.implr  r  r   r<   r<   r=   impl_column_stacku  s    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 )rx  N)r~  r  r<   r<   r=   r    s     z"np_stack_common.<locals>.<genexpr>rx  r   r?   c                 S   s,   |  ||d |jd tj}t| ||||S rJ  )r{   r   r   r|   r  r  r<   r<   r=   r    s    z np_stack_common.<locals>.codegen)r0   r  r   r   )	r  r[  r   rl   ra   r~  r  rO   r  r<   r<   r=   np_stack_common  s     
r  c                 C   s   t | tjrddd}|S d S )Nr   c                 S   s
   t | |S r  )r  r  r<   r<   r=   r     s    zimpl_np_stack.<locals>.impl)r   r  r  r<   r<   r=   impl_np_stack  s    
r  c                 C   s4   t | ||\}}t||}t|}t|||}|S r  )r0   r  r1   r   r   )r  r   r[  ndim_minrl   ra   r~  r  r<   r<   r=   NdStack_typer  s
    
r  c                 C   s&   t | d|d}||}dd }||fS )Nz	np.hstackr?   c                    sj   |j d }|d j}|dkr:| tjd t| ||| S |dkrFdnd  fdd}| ||||S d S )Nr   r?   c                    s   t j|  dS rJ  rL  r[  rK  r<   r=   np_hstack_impl  s    z3_np_hstack.<locals>.codegen.<locals>.np_hstack_implr   ra   r   r   r|   r  r  )rn   r7   rO   r   r  ra   r  r<   rK  r=   r    s    

z_np_hstack.<locals>.codegenr  r  r  r  rO   r  r<   r<   r=   
_np_hstack  s    r  c                 C   s   t | tjrdd }|S d S )Nc                 S   s   t | S r  )r  r  r<   r<   r=   r     s    zimpl_np_hstack.<locals>.implr  r  r<   r<   r=   impl_np_hstack  s    r  c                 C   s&   t | d|d}||}dd }||fS )Nz	np.vstackr  c                 S   sd   |j d }|d j}|dkr&dd }n.|dkrL| tjd}t| ||||S dd }| ||||S )Nr   c                 S   s   t t | dS rJ  )r  rx  hstackr  r<   r<   r=   np_vstack_impl  s    z3_np_vstack.<locals>.codegen.<locals>.np_vstack_implr?   c                 S   s   t j| ddS )Nr   rK  rL  r  r<   r<   r=   r    s    r  )rn   r7   rO   r   r  ra   r  r   r<   r<   r=   r    s    


z_np_vstack.<locals>.codegenr  r  r<   r<   r=   
_np_vstack  s    r  c                 C   s   t | tjrdd }|S d S )Nc                 S   s   t | S r  )r  r  r<   r<   r=   r     s    zimpl_np_vstack.<locals>.implr  r  r<   r<   r=   impl_np_vstack  s    r  c                 C   s&   t | d|d}||}dd }||fS )Nz	np.dstackr*  c                 S   s   |j d }|j}|d j}|dkr:dd }| ||||S |dkr| tjd}|j|jd d}	tj	|	f|j  }
t
| ||
||}| tjd}t	||	}t| |||f|S |dkr| tjd}t
| ||||S dd }| ||||S d S )Nr   c                 S   s   t | dddS )Nr?   r!  )r  r  rC  r  r<   r<   r=   r    s    z3_np_dstack.<locals>.codegen.<locals>.np_vstack_implr?   rz  r  c                 S   s   t j| ddS )Nr  rK  rL  r  r<   r<   r=   r    s    )r   r   ra   r  r   r   r|   r  r   r    r  rx  )rn   r7   rO   r   r  r>  ra   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   s   t | tjrdd }|S d S )Nc                 S   s   t | S r  )r  r  r<   r<   r=   r     s    zimpl_np_dstack.<locals>.implr  r  r<   r<   r=   impl_np_dstack  s    r  fillc                 C   s   dd }|S )Nc                 S   s   || d d < d S r  r<   )r  rx   r<   r<   r=   	fill_impl$  s    zarr_fill.<locals>.fill_implr<   )r  rx   r  r<   r<   r=   arr_fill!  s    r  dotc                 C   s   dd }|S )Nc                 S   s   t | |S r  )r  r  )r  otherr<   r<   r=   dot_impl-  s    zarray_dot.<locals>.dot_implr<   )r  r  r  r<   r<   r=   	array_dot+  s    r  c                 C   s"   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  r(  r!  .r  r  ra   r   r  r  r<   r<   r=   r   9  s    

znp_flip_lr.<locals>.implr  r  r<   r<   r=   
np_flip_lr3  s    r  c                 C   s"   t | std|  dd }|S )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   J  s    

znp_flip_ud.<locals>.implr  r  r<   r<   r=   
np_flip_udD  s    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                    sX   dd }t jg}tjf| }| t j}| }| |g}	| ||||	}
|
S )Nc                 S   s*   |}t | D ]}t||td d d}q|S )Nr!  )r2   r*   r  )lengthempty_tupler   r_   r<   r<   r=   r   `  s    z6_build_flip_slice_tuple.<locals>.codegen.<locals>.impl)r   r|   r   r    r   get_constant_undefr  )rn   r7   r    r   r   inner_argtypes	inner_sigll_idx_typer  
inner_argsr   r   
tuple_typer<   r=   r  _  s    

z(_build_flip_slice_tuple.<locals>.codegen)	r   r   r  r   r1  r   r   r  slice3_type)r  szrO   r  r<   r  r=   _build_flip_slice_tupleU  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  ra   )r  slr<   r<   r=   r   z  s    
znp_flip.<locals>.impl)r   r   r   r   r   r  r<   r<   r=   np_flips  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rht |jjtjrhddd}|S t |tj	rt
dd |jD rddd}|S d S )Nr   c                 S   s   t jt | ||dS rJ  )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?   rK  )divmodrW   r  cumsumr   r  )r   r  r   rT  rA  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+   ra   r   r*   r  r[   r   r  r   Z	slice_tupr   prevr  r   r<   r<   r=   r     s    
c                 s   s   | ]}t |tjV  qd S r  r   r   r   )r   tr<   r<   r=   r    s     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+   ra   r   r   r*   r  r[   r  r<   r<   r=   r     s    
)r   )r   )r   )r   )r   r   r  ListTyper  r   IterableTypeiterator_typer   r  r  r   r  r   r   r<   r<   r=   np_array_split  s*    





r  c                 C   sP   t | tjtjtjfr$ddd}|S t |tjr>ddd}|S t| ||dS d S )Nr   c                 S   s   t jt | ||dS rJ  )r  splitr  r  r<   r<   r=   r     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 divisionrK  )r  rW   r   r  r  )r   r  r   r  rA  r<   r<   r=   r     s      rK  )r   )r   )r   r   r  r  r  r   r  r  r<   r<   r=   np_split  s    


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k
r   | dkrLtj|r8tnd|dd}|j}n$| dkrptj|r`tnd|d}|j}|t |< | Y S X dS )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  s    z#array_sort.<locals>.array_sort_impl)r   r  r   rl   r   r  r  )rn   r7   rO   r   r  r  r<   r  r=   
array_sort   s    
r  c                 C   s   t | stddd }|S )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"  s    z)array_argsort.<locals>.array_argsort_implr?   r  )	r   r  r   r   rl   r   r  replacer  )	rn   r7   rO   r   r  r  r  Zinnersig	innerargsr<   r  r=   array_argsort  s    
 r  c                 C   s   |j |j ks|jdkst|S )Nr  )r  r~  r   )rn   r7   r  r  rx   r<   r<   r=   array_to_array.  s    r  c                 C   s4   dd }t ||}| ||||g}t| ||j|S )Nc                 S   s
   |  dS r   )r   rT  r<   r<   r=   r   =  s    zarray0d_to_scalar.<locals>.impl)r    r  r   r   rn   r7   r  r  rx   r   rO   r   r<   r<   r=   array0d_to_scalar5  s    
r  c                 C   s4   dd }t ||}| ||||g}t| ||j|S )Nc                 S   s   t | d S rK  )r  rT  r<   r<   r=   r   H  s    z array_to_unichrseq.<locals>.impl)r    r  r   r   r  r<   r<   r=   array_to_unichrseqF  s    
r  c                 C   s   t dS )zR
    An intrinsic returning a derived array with the given shape and strides.
    Nr   r  rW   r   r<   r<   r=   reshape_uncheckedS  s    r  c                    s   dd   fdd}|S )Nc                 S   s   t | tjotdd | D S )Nc                 s   s   | ]}t |tjV  qd S r  r  rc  r<   r<   r=   r  ^  s     z>type_reshape_unchecked.<locals>.check_shape.<locals>.<genexpr>)r   r   r  r  rW   r<   r<   r=   rY  \  s    z+type_reshape_unchecked.<locals>.check_shapec                    sJ   t | tjsd S  |r  |s$d S t|t|kr8d S | jt|ddS )Nr  r  )r   r   r   r   r  r  rY  r<   r=   typer`  s    z%type_reshape_unchecked.<locals>.typerr<   )rn   r  r<   r  r=   type_reshape_uncheckedZ  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   rj   r	   r   r   r   r   r   r   r   )rn   r7   rO   r   r   r>  r   r   rW   r   r   r<   r<   r=   impl_shape_uncheckedl  s    
r  c                    s`   |d t jfkrtdd  ntdd  |d t jfkr@tdntdd d	 fdd	}|S )
Nc                 S   s   | j S r  r   r  rW   r<   r<   r=   r2    s    zas_strided.<locals>.get_shapec                 S   s   |S r  r<   r  r<   r<   r=   r2    s    z*as_strided() strides argument is mandatoryc                 S   s   |S r  r<   )r  r   r<   r<   r=   get_strides  s    zas_strided.<locals>.get_stridesc                    s   t |  | || |} | S r  )r  )r  rW   r   r2  r  r<   r=   as_strided_impl  s    z#as_strided.<locals>.as_strided_impl)NN)r   r]  r"   rF   )r  rW   r   r  r<   r  r=   
as_strided  s    


r	  c                    s  t |tjr>tt| jd tt| jd tdd n^t |tjrt |jtjrtt| jt	| tt| jt	| tdd n
t
dt|rtdd  nLt |tjrtdd  n2t |tjrt |jtjrtd	d  n
t
d
d fdd	}|S )Nr?   c                 S   s   | fS r  r<   window_shaper<   r<   r=   get_window_shape  s    z-sliding_window_view.<locals>.get_window_shapec                 S   s   | S r  r<   r
  r<   r<   r=   r    s    z4window_shape must be an integer or tuple of integersc                 S   s   t t|S r  )r   r2   r  r   ra   r<   r<   r=   get_axis  s    z%sliding_window_view.<locals>.get_axisc                 S   s   t dd||gS )Nsliding_window_viewr   r   r  r<   r<   r=   r    s    c                    s    fdd|D S )Nc                    s   g | ]}t d d |qS )r  r   r  r  rz  r<   r=   r     s   z9sliding_window_view.<locals>.get_axis.<locals>.<listcomp>r<   r  r<   rz  r=   r    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 ]n\}}|dk rtd|| |k rtd|| | d }t|||}t|||}t||| j| }|d7 }qzt| ||}	|	S )Nz2Must provide matching length window_shape and axisr   z-`window_shape` cannot contain negative valuesz4window_shape cannot be larger than input array shaper?   )	ra   r   r   r2   r*   rW   r   r   r  )
r  r  r   	out_shapeout_stridesr_   r	  r8  trimmedrs  r  r  Zshape_bufferZstride_bufferr<   r=   sliding_window_view_impl  s8    
z5sliding_window_view.<locals>.sliding_window_view_impl)N)r   r   r   r   r2   ra   r"   r  rl   r   r   r   r   )r  r  r   r  r<   r  r=   r    s>    'r  c                 C   s   t | tjrdd }|S d S )Nc                 S   s6   | j dkrdS | j dkr&t| dS d}t|d S )Nr   Fr?   z[The truth value of an array with more than one element is ambiguous. Use a.any() or a.all())r   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s4tdt | tjsJtd| jtt	  fdd}|S )Nz.The second argument "axis1" must be an integerz-The third argument "axis2" must be an integerr'  c                    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	axes_listra   r<   r=   r     s    znumpy_swapaxes.<locals>.impl)
r   r   r   r   r   r   r   ra   r   r2   )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 ]H\}\}}	|dkrL|	}
n"|	dkrZ|}
n||	krjtd|}
t|||
}q2t||| j| }t	| |} t	||}|}t
|dkrtt
|D ]}t||| j| }q|}t
|dkrtt
|D ]}t||| j|d |  }q|j| }t||f | | j}td d d f}t|D ]h}t|D ]V}| || |  }||| |  }||| |  }t|D ]}|||  ||< qq`qR|S )Nznp.take_along_axisr   r?   z*`arr` and `indices` dimensions don't matchr   )r   ra   r   rW   r  r   r   r*   r  r  r   r2   r  rl   r  rZ  )r  r   r   ZNi_origZNk_origindices_broadcast_shape	arr_shaper_   d1d2new_valZarr_broadcast_shapeNiNkJr   Znp_s_iikkZa_1dZ
indices_1dZout_1djr<   r<   r=   _take_along_axis_impl-  sX    
    
r'  c                    s  t | tjstdt |tjs,tdt |jtjsDtdt|rRd}n| j}||jkrltdt	t
|jt|rfdd}nt|d t |tjstd	|j}|d
k r| j| }|d
k s|| jkrtdt	t
| t	t
|d | j fdd}|S )Nr'  z.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<   )r'  rP  r  r   r   )r  r<   r=   take_along_axis_impl~  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  )r'  r(  r!  r"  r  r<   r=   r)    s    )r   r   r   r   r   rl   r   r   ra   r   r2   r   r  r.  r   )r  r   r   arr_ndimr)  r<   r*  r=   arr_take_along_axisi  s<    






r,  T        c                 C   s   t | tjr`t | tjr$ddd}qt | tjr<ddd}qt | tjrTddd}qtdnlt| rt | j	tjrddd}qt | j	tjrdd	d}qt | j	tjrdd
d}qtdn
td|S )NTr-  c                 S   s   | S r  r<   r  r  nanr<   r<   r=   r     s    znan_to_num_impl.<locals>.implc                 S   sF   t | r|S t | r(t t| jS t | rBt t| jS | S r  )r  r  isneginffinfor3   r[  isposinfr  r.  r<   r<   r=   r     s    


c                 S   s*   t j| j|d}t j| j|d}t||S )N)r/  )r  
nan_to_numr  r  complex)r  r  r/  r  rW  r<   r<   r=   r     s    z4Only Integer, Float, and Complex values are acceptedc                 S   s   | S r  r<   r.  r<   r<   r=   r     s    c           	      S   s   t | jj}t | jj}t | }|r4t |n|}|j}t|j	D ]J}t 
|| rd|||< qHt || r||||< qHt || rH|||< qH|S r  )r  r1  rl   r[  r  r  r  r&  r2   r   r  r0  r2  )	r  r  r/  Zmin_infZmax_infx_outputZoutput_flatr_   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  r/  )r  r  r  r3  r  r  )r  r  r/  r5  r6  r<   r<   r=   r     s
    
z4The first argument must be a scalar or an array-like)Tr-  )Tr-  )Tr-  )Tr-  )Tr-  )Tr-  )
r   r   r  r   r  r  r   r   r   rl   )r  r  r/  r   r<   r<   r=   nan_to_num_impl  s,    
r7  )N)N)N)r"  )r?   )N)N)NNN)NNN)N)N)N)N)N)N)N)N)r   )NNN)r"  )N)N)N)r   )r   )r   )r   )F)NN)N)Tr-  (  r6  r  r  r  llvmliter   llvmlite.irr   numpyr  numbar   r   
numba.corer   r   r   r	   r
   numba.np.numpy_supportr   r   r   r   r   r   r   r   r   r   numba.core.imputilsr   r   r   r   r   r   r   r   r   r   r   numba.core.typingr    numba.core.typesr!   numba.core.extendingr"   r#   r$   r%   r&   
numba.miscr'   r(   numba.cpythonr)   numba.cpython.unsafe.tupler*   r+   r,   Znumba.core.typing.npydeclr-   r  r.   r  r/   r0   r1   r>   rA   rj   rp   rv   ry   r   r   r   r   r   r   r   r   r   r   ArrayIteratorBORROWEDr   r  r   r  r  r   r   r  r  r  r  Anyr#  r   r%  r   r+  r-  rv  r.  r7  rI  rM  rQ  r]  rf  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  broadcast_shapesr  rY  r  r!  r  r  r  r  VarArgr  r  r  r   logspacer&  rot90r/  r:  r=  rA  rB  rC  rE  r[   rM  rS  rF  rV  rW  rb  rd  rv  ru  rw  r  rW   r  uniquer  r  r  r  r  r  r  r  r  r  r  
MemoryViewr  r  r  r  r  r  r  r  r  r  r  r  r  rO  r  r  r  r  r  r  r  r  r  r  r  r  DTypeNumberClassr  r  r  r  r  r  r	  r  r  r  r  r  r  Bytesr  is_r#  r%  r*  r.  r8  r9  rw  r~  r(  r  r  r  r  r  r  r+  r  r  r  rZ  r  r  r  r  r:  r  NumpyNdIterTyper  r|  r  r  r  r  r  r  r  r  r  r  r  r  rX  r  r  r  r  r  r  r  r  r  r  r  rl  r  r  r  r  r  r  r  eyer  diagr  r   r
  r  aranger  r#  r   r#  r$  r%  r  r'  r/  r2  r3  r7  r8  r9  r;  rC  rD  rE  rL  rM  rJ  rP  rX  r\  ra  rf  ri  rj  rl  r   rm  rn  rt  rv  rx  ry  rz  r  r~  
atleast_1dr  
atleast_2dr  
atleast_3dr  r  r  r  r  r  rG  r  r  r  column_stackr  r  r  stackr  r  r  r  r  r  vstackr  r  dstackr  r  r  r,  r  r-  r  r  flipr  r  r  r  r  r  r  r  r  r  r  argsortr  r  UnicodeCharSeqr  r  rh  NPTimedelta
NPDatetimer  r  r  type_callabler  r  libstride_tricksr	  r  r   r  r+  r  r'  take_along_axisr,  r3  r7  r<   r<   r<   r=   <module>   s  $4<	
E




D

(



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

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






















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

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