U
    |eCj                     @   s(  d dl mZ d dlZd dlmZmZ d dlmZ d dl	m
Z
mZ d dlmZ d dlmZmZ edd	d
dgZdd Zdd Zdd Zdd Zdd Zdd ZG dd deZedZG dd deZG dd deZdd  Zd!d" ZG d#d$ d$eZG d%d& d&eZ G d'd( d(eZ!G d)d* d*e!Z"dS )+    )
namedtupleN)Constant	IRBuilder)ir)typescgutils)global_compiler_lock)make_library_cache	NullCache_wrapper_infolibraryenvnamec              
   C   s   |  }|j |||	j|	j|\}}|j|jddT\}}| || W 5 Q R X |( |
 }|j ||
| |
| W 5 Q R X W 5 Q R X t	||D ]$\}}|
||||j| q|
||||j| |jS NTlikely)	call_convcall_functionreturn_typeargsif_elseis_ok
gil_ensureraise_errorgil_releasezipstoreaddloadstepcode)r   r   contextfuncbuilderarraysoutoffsetsstore_offset	signaturepyapir   elemsstatusretvalif_okif_errorgiloffary r2   T/var/www/website-v5/atlas_env/lib/python3.8/site-packages/numba/np/ufunc/wrappers.py_build_ufunc_loop_body   s&    
 r4   c              	   C   s   |  }t jgt|	j }|jdd4 |j||t j||\}}|D ]}|| qDW 5 Q R X || t||D ]$\}}|	|
|||j| qp|	|
|||j| |jS )NT)keep_new)r   pyobjectlenr   err_pushr   r   decrefr   r   r   r   r   r    )r   r   r!   r"   r#   r$   r%   r&   r'   r(   r   r)   r*   Z_objargsr+   r,   elemr0   r1   r2   r2   r3   _build_ufunc_loop_body_objmode,   s"    
 r;   c
                    s@    fdd}
fdd}t |
|| | |||	dS )Nc                     s   fddt  D } | S )Nc                    s    g | ]\}}|  |qS r2   load_directr   .0r0   r1   r#   r2   r3   
<listcomp>O   s   z6build_slow_loop_body.<locals>.load.<locals>.<listcomp>)r   r*   )r$   r#   r&   r2   r3   r   N   s    
z"build_slow_loop_body.<locals>.loadc                    s    |   d S N)store_directr   r,   )r#   r%   r'   r2   r3   r   S   s    z#build_slow_loop_body.<locals>.storer   r4   )r!   r"   r#   r$   r%   r&   r'   r(   r)   r   r   r   r2   )r$   r#   r&   r%   r'   r3   build_slow_loop_bodyL   s        rH   c                    sb   |  |	}|
||	 fdd}fdd}t||| | |	S )Nc                     s8   fddt  D } fddt | jD } | S )Nc                    s    g | ]\}}|  |qS r2   r<   r>   r@   r2   r3   rA   b   s   z5build_obj_loop_body.<locals>.load.<locals>.<listcomp>c                    s   g | ]\}} || qS r2   )from_native_value)r?   vt)env_managerr)   r2   r3   rA   e   s   )r   r   rB   )r$   r#   rL   r&   r)   r(   r2   r3   r   `   s    

z!build_obj_loop_body.<locals>.loadc              	      sd   t  | } j|dd@ j| }|jd ks8t|j 	 
|  W 5 Q R X d S r   )r   is_not_nullif_thento_native_valuer   cleanupAssertionErrorrD   valuer   r9   )r,   r   native)r#   r%   r)   r(   r'   r2   r3   r   i   s    z"build_obj_loop_body.<locals>.store)get_env_bodyget_env_managerr;   )r!   r"   r#   r$   r%   r&   r'   r(   r)   envptrr   env_bodyr   r   r2   )r$   r#   rL   r&   r%   r)   r(   r'   r3   build_obj_loop_body[   s    	     rX   c                    s<    fdd}fdd}t ||| || ||||	|
dS )Nc                     s   fdd D } | S )Nc                    s   g | ]}|  qS r2   )load_aligned)r?   r1   )indr2   r3   rA   ~   s   z6build_fast_loop_body.<locals>.load.<locals>.<listcomp>r2   rB   )r$   rZ   r2   r3   r   }   s    
z"build_fast_loop_body.<locals>.loadc                    s    |   d S rC   )store_alignedrE   )rZ   r%   r2   r3   r      s    z#build_fast_loop_body.<locals>.storerF   rG   )r!   r"   r#   r$   r%   r&   r'   r(   rZ   r)   r   r   r   r2   )r$   rZ   r%   r3   build_fast_loop_body{   s        r\   c           *      C   sB  t |tsttd}t|}t|}|tj}	t|	}
t	t
 ||
|
|g}| d}|d}|r|jtjtjgt|j }n|j|j|j}tj|||d}|jd t||d|j }|j\}}}}d|_d|_d	|_d
|_t|d}||j}|j}|||j|}|j|dd}g }t |jD ]"\}}|!t"|||||| qNt"||||t||j}g }|#tjd}|D ](} t$%||	}!|!|! |&||! qt$%||	}"|&||" t$j'}#|D ]}$|(|#|$j)}#q|*|}%|rf|%+ }&t$j,|||	d" t-|||||||"||%|| W 5 Q R X |%.|& |/  n|0|#\}'}(|'B t$j,|||	d&})t1|||||||"||)j2|%|d W 5 Q R X W 5 Q R X |(> t$j,|||	d" t3|||||||"||%|d
 W 5 Q R X W 5 Q R X W 5 Q R X |/  ~|4| |5|  t6|||jdS )z
    Wrap the scalar function with a loop that iterates over the arguments

    Returns
    -------
    (library, env, name)
       Zufunc_wrapper r   alwaysinlinez
__ufunc__.r   dimsstepsdataentry	loopcountr   intprF   r   r   r   )7
isinstancestrrQ   r   IntTypePointerTypeget_value_typer   rg   FunctionTypeVoidTypecodegencreate_librarycreate_ir_moduler   get_function_typer6   r7   r   r   Function
attributesr   r   r   append_basic_blockget_env_namefndescenvironmentr   declare_env_globalmodule	enumerateappend	UArrayArgget_constantr   alloca_oncer   true_bitand_is_unit_stridedget_python_apir   	for_rangerX   r   ret_voidr   r\   indexrH   add_ir_moduleadd_linking_libraryr   )*r   r!   fnamer(   objmodecresbyte_t
byte_ptr_tbyte_ptr_ptr_tintp_t
intp_ptr_tfnty
wrapperlibwrapper_module	func_typer"   wrapperarg_argsarg_dims	arg_stepsarg_datar#   envnamer   rV   re   r$   itypr%   r&   zero_pr'   Zunit_stridedr1   r)   r/   r   Z
is_stridedloopr2   r2   r3   build_ufunc_wrapper   s    



 
  

         

               $

r   c                   @   s4   e Zd Zdd Zdd Zdd Zdd Zd	d
 ZdS )r~   c                 C   s   || _ || _|| _| j tj|}| j|||g}||}	| j	||	
 | _| j |	}
| j tj|
| _| j||g}| j|| _|d| j| j| _|| _d S )Nz==)r!   r#   fe_typer   r   rg   r   gepget_data_typebitcast
as_pointerdataptrget_abi_sizeofZabisizer   icmp_unsignedr   )selfr!   r#   r   rb   r   r   offsetZoffseted_args	data_typesizeofZoffseted_stepr2   r2   r3   __init__   s$    
 zUArrayArg.__init__c                 C   s&   t | j| j|}| j| j| j|S )zm
        Generic load from the given *byteoffset*.  load_aligned() is
        preferred if possible.
        )r   pointer_addr#   r   r!   unpack_valuer   )r   
byteoffsetptrr2   r2   r3   r=     s    zUArrayArg.load_directc                 C   s&   | j | j|g}| j| j | j|S rC   )r#   r   r   r!   r   r   )r   rZ   r   r2   r2   r3   rY     s    zUArrayArg.load_alignedc                 C   s,   t | j| j|}| j| j| j|| d S rC   )r   r   r#   r   r!   
pack_valuer   )r   rR   r   r   r2   r2   r3   rD     s    zUArrayArg.store_directc                 C   s,   | j | j|g}| j| j | j|| d S rC   )r#   r   r   r!   r   r   )r   rR   rZ   r   r2   r2   r3   r[     s    zUArrayArg.store_alignedN)__name__
__module____qualname__r   r=   rY   rD   r[   r2   r2   r2   r3   r~      s
   r~   Zgufc                   @   s   e Zd Zdd Zedd Zedd Zedd Zed	d
 Zedd Z	edd Z
dd Zdd Zdd Zedd Zdd Zdd Zdd ZdS )_GufuncWrapperc                 C   sN   || _ || _|| _|| _| jjtjk| _|r8t	| j dnt
 | _t|| _dS )a,  
        The *is_parfors* argument is a boolean that indicates if the GUfunc
        being built is to be used as a ParFors kernel. If True, it disables
        the caching on the wrapper as a separate unit because it will be linked
        into the caller function and cached along with it.
        )py_funcN)r   r   sinsoutr(   r   r   r6   Zis_objectmodeGufWrapperCacher
   cachebool
is_parfors)r   r   r   r   r   r   r   r2   r2   r3   r   &  s    z_GufuncWrapper.__init__c                 C   s   | j jS rC   )r   r   r   r2   r2   r3   r   6  s    z_GufuncWrapper.libraryc                 C   s   | j jS rC   )r   target_contextr   r2   r2   r3   r!   :  s    z_GufuncWrapper.contextc                 C   s   | j jS rC   )r!   r   r   r2   r2   r3   r   >  s    z_GufuncWrapper.call_convc                 C   s   | j jS rC   )r   r(   r   r2   r2   r3   r(   B  s    z_GufuncWrapper.signaturec                 C   s   | j jS rC   )r   rx   r   r2   r2   r3   rx   F  s    z_GufuncWrapper.fndescc                 C   s   | j jS rC   )r   ry   r   r2   r2   r3   r   J  s    z_GufuncWrapper.envc                 C   sR   t d}t |}t |}| jtj}t |}t t  ||||g}|S )Nr]   )	r   rk   rl   r!   rm   r   rg   rn   ro   )r   r   r   r   r   r   r   r2   r2   r3   _wrapper_function_typeN  s    



 z%_GufuncWrapper._wrapper_function_typec           !         sr  | j tj}|  }|d}| j| jj	| jj
}| jj}tj|||d}|jd t|||}	d|	_|	j\}
}}}d|
_d|_d|_d|_t|	d	}|j|d
d}| j |}t }| j| jfD ]}|D ]}|t|O }qqi }| jD ]"}|D ]}||krt|||< qqi }| D ]2\}}|||| j tj|d g||< q$g }t| jt| j }tt| j j| j| j D ]>\}\}}t!| j ||
||||||	}|t|7 }|"| q|d}| #|| t$j%|||d:  fdd|D }| &||||\}} t$'|| | W 5 Q R X |(| |)| | *|| |+  |,| |-| j. dS )z
        The LLVM IRBuilder code to create the gufunc wrapper.
        The *library* arg is the CodeLibrary to which the wrapper should
        be added.  The *name* arg is the name of the wrapper function being
        created.
        Z_gufunc_wrapperr_   r`   weak_odrr   ra   rb   rc   rd   re      z.returnrf   c                    s   g | ]}|  jqS r2   )get_array_at_offsetr   )r?   ar   r2   r3   rA     s     z1_GufuncWrapper._build_wrapper.<locals>.<listcomp>N)/r!   rm   r   rg   r   rr   r   rs   rx   restypeargtypesllvm_func_namer   rt   ru   r   linkager   r   r   rv   r   r   setr   r   r7   itemsr   r   r|   r   r(   
GUArrayArgr}   gen_prologuer   r   gen_loop_bodycbranch_or_continuebranchposition_at_endgen_epiloguer   r   r   r   )!r   r   r   r   r   r   r   r   r"   r   r   r   r   r   r#   re   r)   Zunique_symsgrpsymsZsym_mapssym_dimr   r$   step_offsetr   symr1   Zbbreturnr   	innercallerrorr2   r   r3   _build_wrapperY  s~    





     



z_GufuncWrapper._build_wrapperc                 C   s   | j r(| j t| }| || nX| j| jj	| jj
}|d kr| j t| }|  | || | j| jj	| |S rC   )r   r!   rp   rq   rj   r   r   load_overloadr   r(   r   enable_object_cachingsave_overloadr   wrapper_namer   r2   r2   r3   _compile_wrapper  s     z_GufuncWrapper._compile_wrapperc                 C   s&   d| j j }| |}t|| j|dS )Nz__gufunc__.rh   )rx   mangled_namer   r   r   r   r2   r2   r3   build  s    
  z_GufuncWrapper.buildc              	   C   sl   | j ||| jj| jj|\}}|j|jdd* | }| jj 	||| |
| W 5 Q R X |j|jfS )NFr   )r   r   r(   r   r   rN   is_errorr   r!   r   r   r    )r   r#   r)   r"   r   r+   r,   r/   r2   r2   r3   r     s       z_GufuncWrapper.gen_loop_bodyc                 C   s   d S rC   r2   r   r#   r)   r2   r2   r3   r     s    z_GufuncWrapper.gen_prologuec                 C   s   d S rC   r2   r   r2   r2   r3   r     s    z_GufuncWrapper.gen_epilogueN)r   r   r   r   propertyr   r!   r   r(   rx   r   r   r   r   r   r   r   r   r   r2   r2   r2   r3   r   %  s*   





T
r   c                   @   s$   e Zd Zdd Zdd Zdd ZdS )_GufuncObjectWrapperc                 C   s"   t | j|||| j|\}}||fS rC   )_prepare_call_to_object_moder!   r(   )r   r#   r)   r"   r   r   r   r2   r2   r3   r     s      z"_GufuncObjectWrapper.gen_loop_bodyc                 C   s   |  | _d S rC   )r   r/   r   r2   r2   r3   r     s    z!_GufuncObjectWrapper.gen_prologuec                 C   s   | | j d S rC   )r   r/   r   r2   r2   r3   r     s    z!_GufuncObjectWrapper.gen_epilogueN)r   r   r   r   r   r   r2   r2   r2   r3   r     s   r   c                 C   s2   |j }|jtjkrtnt}|| |||||d S )N)r   )r(   r   r   r6   r   r   r   )r   r   r   r   r   r   r(   Zwrapclsr2   r2   r3   build_gufunc_wrapper  s    
     r   c           &   
   C   s\  |j }|d}| tj}| tj}	t|	}
| tj}| tj	}t
|||
|
|||g}t||d}tj|tddd}|tj| g }g }tt||jD ]&\}\}}tj||dd}|| t|tjr| |}|| ||d}t|d	}t||j}||d
||g}||d||g}||j|}tt |j}t||j!}t||j"}|#|||||||g} n|$||} || | ||  t%|| }!||!| t&||!| qtj	gt'| }"| j()||tj	|"|\}#}$||#j*| |+|$ |,| |-| |D ]}|+|.| q0|#j/}%|%|.|fS )Nzufunc.core.returnZnumba_ndarray_newr   r   r_   T)zfill)rR   r   shapestrides)0r{   rv   rm   r   int32rg   r   rl   voidptrr6   rn   r   get_or_insert_functionr   rk   r   r   r|   r   r   r}   ri   Array
make_arrayr   ndimr   _get_ptr_by_namer   rc   npdtyperj   numitemsizecallrI   is_nullr   r7   r   r   r   r9   r   r   r   r    )&r!   r#   r)   r"   r(   r   modZbb_core_returnZll_intZll_intpZll_intp_ptrZ
ll_voidptrZll_pyobjr   Zfn_array_newZerror_pointerZobject_argsZobject_pointersr   argargtyZobjptraryclsarrayr   ndra   r   rc   r   Ztype_numr   objobj_is_nullZ
object_sigr+   r,   r   r2   r2   r3   r     s|    



  


 
   


r   c                   @   s   e Zd Zdd Zdd ZdS )r   c
                    sV  || _ || _|tj|}
|j|j||
gdddd}|| _|j||
gdd}||}t|tj	r*| }t
||jkrt
|dkr|jdkrntd|d |j} fdd	|D }g }t|D ]8}|j||tj|| gd
d}||}|| q|rtnt}||j|||||d| _n(|rDtd||d t||d| _d S )Nzdata.ptrr_   rc   zcore.step.ptrr   r   z.type and shape signature mismatch for arg #{0}c                    s   g | ]} | qS r2   r2   )r?   r   r   r2   r3   rA   {  s     z'GUArrayArg.__init__.<locals>.<listcomp>zstep.ptrr   r   	core_step	as_scalarr   r   z2scalar type {0} given for non scalar argument #{1}r   stride)r!   r#   r   r   rg   r   r   rc   ri   r   r7   r   	TypeErrorformatranger}   _ArrayAsScalarArgLoader_ArrayArgLoaderr   _loader_ScalarArgLoader)r   r!   r#   r   rb   r   r   r   r   r   r   rc   Zcore_step_ptrr  r  r   r   r   jZstepptrr   Zldclsr2   r  r3   r   \  s\    


 zGUArrayArg.__init__c                 C   s   | j j| j| j| j|dS )N)r!   r#   rc   rZ   )r  r   r!   r#   rc   )r   rZ   r2   r2   r3   r     s     zGUArrayArg.get_array_at_offsetN)r   r   r   r   r   r2   r2   r2   r3   r   [  s   ;r   c                   @   s    e Zd ZdZdd Zdd ZdS )r  z
    Handle GFunc argument loading where a scalar type is used in the core
    function.
    Note: It still has a stride because the input to the gufunc can be an array
          for this argument.
    c                 C   s   || _ || _d S rC   r  )r   r   r	  r2   r2   r3   r     s    z_ScalarArgLoader.__init__c                 C   s:   | |||| jg}|||| j }||S rC   )r   mulr	  r   r   r   r   r   )r   r!   r#   rc   rZ   Zdptrr2   r2   r3   r     s
    z_ScalarArgLoader.loadN)r   r   r   __doc__r   r   r2   r2   r2   r3   r    s   r  c                   @   s(   e Zd ZdZdd Zdd Zdd ZdS )	r  zD
    Handle GUFunc argument loading where an array is expected.
    c                 C   s(   || _ || _|| _|| _|| _|| _d S rC   r  )r   r   r   r  r  r   r   r2   r2   r3   r     s    z_ArrayArgLoader.__init__c              	   C   s   t j| j| jdd}||}|||}t|||| j|}| 	||\}	}
|
|| j}|j||||jj|	|
|t j|d d | S )NA)r   r   layout)rc   r   r   r   meminfo)r   r   r   r   r   r   r   r  r  _shape_and_stridesr   r   populate_arrayr   rc   typer   rg   	_getvalue)r   r!   r#   rc   rZ   Zarytypr   r   Zoffseted_datar   r   r   r2   r2   r3   r     s.    

	z_ArrayArgLoader.loadc                 C   s$   t || j}t || j}||fS rC   )r   
pack_arrayr   r   )r   r!   r#   r   r   r2   r2   r3   r    s    z"_ArrayArgLoader._shape_and_stridesN)r   r   r   r  r   r   r  r2   r2   r2   r3   r    s   r  c                   @   s   e Zd ZdZdd ZdS )r  z
    Handle GUFunc argument loading where the shape signature specifies
    a scalar "()" but a 1D array is used for the type of the core function.
    c                 C   s@   | tjd}| tjd}t||g}t||g}||fS )Nr   r   )r   r   rg   r   r  )r   r!   r#   oner   r   r   r2   r2   r3   r    s
    z*_ArrayAsScalarArgLoader._shape_and_stridesN)r   r   r   r  r  r2   r2   r2   r3   r    s   r  )#collectionsr   numpyr   llvmlite.irr   r   llvmliter   
numba.corer   r   numba.core.compiler_lockr   numba.core.cachingr	   r
   r   r4   r;   rH   rX   r\   r   objectr~   r   r   r   r   r   r   r  r  r  r2   r2   r2   r3   <module>   s0     o) <
`A+