U |evc@s dZddlZddlZddlZddlmZddlmZddlZ ddl Z ddl m Z m Z mZddlmZmZmZddlmZmZmZmZmZddlmZmZmZmZddlmZdd l m!Z!m"Z"dd lm#Z#dd l$m%Z%ed Z&Gd dde'Z(Gddde'Z)GdddeddZ*GdddeddZ+d9ddZ,e-ej-e.ej-ej-e.ej-Z/ddZ0ddZ1ddZ2d d!Z3d"d#Z4d$d%Z5Gd&d'd'e'Z6d(d)Z7d*d+Z8d:d-d.Z9d;d/d0Z:e&;e jd3d4Z?e?e @e&j;e"d5d6ZAe!e jBd7d8ZCdS))arrayr*r)r-r r.r/r0r4widthlenr@zipr1 icmp_unsignedif_thenrB) rrrr8r:ONEr@srcdstdimcondrrrr^s z#_ArrayIndexingHelper.update_indicescs|jjfdd|jDS)z The indexing helper is built using alloca for each value, so it actually contains pointers to the actual indices to load. Note that update_indices assumes the same. This method returns the indices as values csg|]}|qSrrD.0indexr8rr tsz2_ArrayIndexingHelper.as_values..)rJr*r@r rrYrr!lsz_ArrayIndexingHelper.as_valuesNr"rrrrrH\srH)rJr@c@s0eZdZdZddZddZddZdd Zd S) _ArrayHelperzHelper class to handle array arguments/result. It provides methods to generate code loading/storing specific items as well as support code for handling indices. cCsn|jtj}tt|jd}g}t|j D]2}t |j t|j}|j ||||q0t||SNr)r)r-r r.r/r0r4rKrangendimr r5r*rBappendrH)rr:ZEROr@ixrrrr=s z _ArrayHelper.create_iter_indicesc Cs$tj|j|j|j|j|j|j|dS)N)datar1strideslayoutinds)r get_item_pointer2r)r*rcr1rdrer?rrr_load_effective_addresssz$_ArrayHelper._load_effective_addresscCs&|jj|j}||}||j|Sr)r)data_model_managerr,rhload_from_data_pointerr*)rr@modelptrrrrrAs z_ArrayHelper.load_datacCsH|j}|j}|||j|}||j|jks2t||||dSr) r)r*get_value_as_datar,r3typeAssertionErrorrBrh)rr@valuectxr8Z store_valuerrrrCs z_ArrayHelper.store_dataN)r#r$r%rFr=rhrArCrrrrr[ws   r[) r)r*r1rdrcrer,r^rE input operandc Cst|tjr&|}|j}|||||}t|tjr|||||}t||j |j }t||j |j }t |||||j |j|j|j | St|tjtjhBkst|tjjrt||||Std|t|dS)zreturns an instance of the appropriate Helper (either _ScalarHelper or _ArrayHelper) class to handle the argument. using the polymorphic interface of the Helper classes, scalar and array cases can be handled with the same codezunsupported type for {0}: {1}N) isinstancer OptionalrncastArrayCompatible make_arrayr unpack_tupler1r^rdr[rcredtype unliteral number_domainr2scalars_NPDatetimeBaser&NotImplementedErrorformatstr) r7r8inpZtyinpwhereotyaryr1rdrrr_prepare_arguments*   rcCsz||kr dSd}||}||krv||}||}|dkrT||krd|dkrd|d Sn|dkrd|||<|d7}|d7}q|S)aLow-level utility function used in calculating a shape for an implicit output array. This function assumes that the destination shape is an LLVM pointer to a C-style array that was already initialized to a size of one along all axes. Returns an integer value: >= 1 : Succeeded. Return value should equal the number of dimensions in the destination shape. 0 : Failed to broadcast because source shape is larger than the destination shape (this case should be weeded out at type checking). < 0 : Failed to broadcast onto destination axis, at axis number == -(return_value + 1). rr(r)src_ndim src_shape dest_ndim dest_shape src_index dest_indexZ src_dim_sizeZ dest_dim_sizerrr_broadcast_ontos  rc sRdd|D}tj}fdd}|d}|d}t||jd} ||j} t||jdtfd d t|jD} | D]} || qt |D]\} }t |d sq||j}t|jD] }|j |t | |qȈ tt|| | g}td ||Nd | f}tjdd}|dk rP|d|j|jf7}jt|fW5QRXq|j}tfdd | D}t||}t|}||}zd|||}Wn0tk r|j tj!j kr|"}YnX||j#|"f}||}|j}t$|j |}t$|j%|}t&|||j'|j(|j)|| S)zUtility function to handle allocation of an implicit output array given the target context, builder, output array type, and a list of _ArrayHelper instances. cSs"g|]}t|tjr|jn|qSr)rsr rtrnrWrbrrrrZsz _build_array..cstj|Sr) get_constantr r.r>)r)rrmake_intp_constsz%_build_array..make_intp_constrr(rrc3s|]}t|VqdSr)r gep_inboundsrV)r*rrr sz_build_array..r^.z+Not enough inputs to {}, expected {} got {}cSsg|] }|jqSrr,rWarrrrZvscSsg|] }|jqSrrrrrrrZwscSsg|] }|qSr)r=rWrrrrrZ|srcSsg|]}t|tr|jqSr)rsr[rerrrrrZs cSsg|]}|dkr|qS)CrrrrrrZscSsg|]}|dkr|qS)FrrrrrrZsrr)r.ordercSsg|] }|jqSr)rE)rWorrrrZs)%rMargsrLnin RuntimeErrorrr#rrrsr rvrrr/r0r-r_ enable_nrtnrtincrefrrorr.r1r loop_nestrrAr!rgeneraterCrrr)r)r*rrrZ kernel_class argumentsout_iret_tyoutputroutputs outer_sigkernelr:r@Z loopshapeZ input_layoutsZ num_c_layoutZ num_f_layoutrrZvals_inrarXrZvals_outZval_outoutrrrnumpy_ufunc_kernelOsf           rc@seZdZddZddZdS)_KernelcCs||_||_||_dSr)r)r*r)rr)r*rrrrr<sz_Kernel.__init__cCsXt|tjrDt|tjsD|j}|j|d}||j|j||d}|}|j|j|||S)zNumpy uses cast semantics that are different from standard Python (for example, it does allow casting from complex to float). This method acts as a patch to context.cast so that it allows complex to real/int casts. real)rsr Complexunderlying_floatr) get_getattrr*ru)rr+fromtytotynewtyattrrrrrus  z _Kernel.castN)r#r$r%r<rurrrrrsrcsGfdddtS)aUse the ufunc loop type information to select the code generation function from the table provided by the dict_of_kernels. The dict of kernels maps the loop identifier to a function with the following signature: (context, builder, signature, args). The loop type information has the form 'AB->C'. The letters to the left of '->' are the input types (specified as NumPy letter types). The letters to the right of '->' are the output types. There must be 'ufunc.nin' letters to the left of '->', and 'ufunc.nout' letters to the right. For example, a binary float loop resulting in a float, will have the following signature: 'ff->f'. A given ufunc implements many loops. The list of loops implemented for a given ufunc can be accessed using the 'types' attribute in the ufunc object. The NumPy machinery selects the first loop that fits a given calling signature (in our case, what we call the outer_sig). This logic is mimicked by 'ufunc_find_matching_loop'. cs(eZdZfddZddZZS)z'_ufunc_db_function.._KernelImplcsvt||||t|jtt|}||j|_ t |j |j |_ |j dkrrd}t|j|dS)Nz2Don't know how to lower ufunc '{0}' for loop '{1}')superr<rrrrget_ufunc_infor ufunc_sigfnrrr inner_sigr~rr#)rr)r*rloopr) _KernelImpl __class__rrrr<s z0_ufunc_db_function.._KernelImpl.__init__c sj}j}fddt||j|jD}tjdjj||}W5QRXjj}||j  j|} ||j |j S)Ncs g|]\}}}|||qSr)ru)rWr+intyouttyr rrrZszD_ufunc_db_function.._KernelImpl.generate..numpy) rrrMrrr)rr*rir from_returnru)rrisigosig cast_argsresdmmrr rrs z0_ufunc_db_function.._KernelImpl.generate)r#r$r%r<r __classcell__rrr)rrrs r)r)rrrr_ufunc_db_functionsrcsXfdd}tj}|fj}tjdD]$}tjf|}|f|||q.S)Ncst||||Srrr)r*rrrrrrdo_ufuncsz'register_ufunc_kernel..do_ufuncr()r Anyrr]rr)rrlowerr_anyin_argsZn_explicit_outout_argsrrrregister_ufunc_kernels  rFcs.|rtfdd}tj}||||dS)Ncst||||Srrrrrrlower_unary_operatorsz.lower_unary_operator)ror r)operatorrrrinplacer _arr_kindrrrregister_unary_operator_kernelsrc svfdd}fdd}tj}tj}||f||f||fg} | D].} |s^||f| |qB||f| |qBdS)Ncst||||Srrrrrrlower_binary_operator sz>register_binary_operator_kernel..lower_binary_operatorcsDt||df}tj|jf|j|jdf}t||||Sr\)rr signaturerrrrrrrlower_inplace_operator s z?register_binary_operator_kernel..lower_inplace_operator)r rr) oprrrrrrrrZ formal_sigsrrrrregister_binary_operator_kernelsrcCs$Gdddt}t||||tj|S)zLowering function for +(array) expressions. Defined here (numba.targets.npyimpl) since the remaining array-operator lowering functions are also registered in this module. c@seZdZddZdS)z1array_positive_impl.._UnaryPositiveKernelcWs |\}|Srr)rrr+rrrr(sz:array_positive_impl.._UnaryPositiveKernel.generateN)r#r$r%rrrrr_UnaryPositiveKernel'sr)rrnppositive)r)r*rrrrrrarray_positive_impl!src Csi}|D]}t|}t|||||<qtjjtjjfD]d}|D]V\}}tt|}||}|j dkrvt ||||qB|j dkrt ||||qBt dqBq6tj jfD]n}|D]`\}}tt|}||}|j dkrt ||||ddq|j dkr t ||||ddqt dqqdS)Nr(z4There shouldn't be any non-unary or binary operatorsT)r)rrrNumpyRulesUnaryArrayOperator_op_mapNumpyRulesArrayOperatoritemsgetattrrrrrrNumpyRulesInplaceArrayOperator) ZufuncsrZkernelsrZdb_funcr r ufunc_namerrrrregister_ufuncs0s8         rcsbfdd}ttjr0ttj}||Sttjjr^t j}tt|}||SdS)zPGiven a string or NumberClass description *desc*, returns the dtype object. cs"t|}|}dd}||fS)NcSs|Sr)get_dummy_value)r)r*rrrrrcodegen]sz9_make_dtype_object..from_nb_type..codegen)r DType)nb_typerrrdescrr from_nb_typeYs z(_make_dtype_object..from_nb_typeN) rsr Literalrrry literal_value functions NumberClassr) typingctxrrrthestrrrr_make_dtype_objectUs  rcCs4t|tjtjjfr dd}|Std|dS)zKProvide an implementation so that numpy.dtype function can be lowered. cSst|Sr)rrrrrimpssznumpy_dtype..impzunknown dtype descriptor: {}N)rsr rrrrNumbaTypeErrorr)rr rrr numpy_dtypensr")rr)F)F)DrFmathsys itertools collectionsr llvmlite.irr/rrrnumba.nprrrnumba.core.imputilsrrr numba.corer r r r r numba.np.numpy_supportrrrrnumba.core.typingrnumba.core.extendingrrr numba.cpythonrregistryobjectrr&rHr[rr.CPointerrrrrrrrrrrrrrposrrr get_ufuncsrryr"rrrrs^     ))  (UQ5  "