NWc@@sddlmZddlmZddlmZddlmZddlmZddl mZm Z ddl m Z dd l m Z dd l mZdd l mZdd l mZdd ZdefdYZdDZidd6dd6ed6ed6ZdZdZdZdZdZdZdZdZdZdCedZ!de"fd YZ#d!Z$d"Z%d#Z&d$Z'd%Z(d&Z)d'Z*d(Z+d)Z,d*Z-d+Z.e.e j/Z0d,Z1d-Z2d.Z3d/Z4d0Z5d1e"fd2YZ6d3Z7dCd4Z8dCd5Z9e:d6e;e j/Z<e9d7d8e<Z=e9d9Z>e9d:Z?e9d;Z@e9d<d=d>ZAe9d?d8e<d@eAgZBe9dAd8id@eAgZCe9dBd8id@eAgZDdCS(Ei(tabsolute_importi(tCythonTransform(t ModuleNode(t CompileError(tCythonUtilityCode(t UtilityCodetTempitaUtilityCode(tOptions(t Interpreter(t PyrexTypes(tNaming(tSymtabcC@shddlm}||}|dkrdd|}djg|jdD]}||^qH}n|S(Ni(tdedentt s (ttextwrapR tjointsplit(ttexttreindentR tindenttx((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyR s    2tIntroduceBufferAuxiliaryVarscB@s8eZeZeZdZdZdZdZRS(cC@spt|tstd|_tt|j|}|jrlt|j t |j |j j t n|S(Ni( t isinstanceRtAssertionErrortmax_ndimtsuperRt__call__tbuffers_existstuse_bufstruct_declare_codetscopetuse_py2_buffer_functionstuse_utility_codetempty_bufstruct_utility(tselftnodetresult((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyR!s    c @sZjj}g|D]\jjr^q}t|dkrk|jddt|_ng|D]\jjrr^qr}t|dkrt|_nx?|D]7\dkrt j t rt|_ PqqW~t t r/t|dkr/tjdnx|D] jjjr`tjdnjj}|jtjkrtjdtjn|j|jkr|j|_nfd}tjtjftjtjff}g|D]\} } || | ^q\} } tj| | _q6W|_ |_!dS( NitkeycS@s|jS(N(tname(tentry((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyt7st memoryviews'Buffer vars not allowed in module scopes-Buffers with pointer types not yet supported.s1Buffer ndims exceeds Options.buffer_max_dims = %dc @sRj|}jddd|d|dj}jrNt|_n|S(NR%tcnamettypetpos(tmanglet declare_vartNoneR+tis_argtTruetused(R*tprefixR)taux_var(R&R%R"R(s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytdecvarUs   ("tentriestitemsR*t is_buffertlentsortR0Rtis_memoryviewsliceRtutility_code_definitionRtusing_memoryviewRRR+tdtypetis_ptrR%tndimRtbuffer_max_dimsRR tc_pyx_buffer_nd_typeR tpybuffernd_prefixtc_pyx_buffer_typetpybufferstruct_prefixR t BufferAuxt buffer_auxtbuffer_entriesR( R!R"Rt scope_itemstbufvarstmemviewslicevarstbuftypeR4tauxvarsR*R2t pybufferndtrcbuffer((R&R%R"Rs9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyt handle_scope0s@+ +  !     . cC@s$|j||j|j||S(N(RORt visitchildren(R!R"((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytvisit_ModuleNodegs cC@s$|j||j|j||S(N(ROt local_scopeRP(R!R"((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytvisit_FuncDefNodels ( t__name__t __module__tFalseRR<RRORQRS(((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyRs   7 R=R?tmodetnegative_indicestcasttfulls"%s" is not a buffer optionsToo many buffer optionss#"%s" buffer option already supplieds "%s" missings[Only allowed buffer modes are: "c", "fortran", "full", "strided" (as a compile-time string)s#ndim must be a non-negative integers0dtype must be "object", numeric type or a structs"%s" must be a booleanc @s`|dkrt}ntj||d|dd\}}t|tkret|ddtnixL|jD]>\}\}}|t krt|t |n||pybuffer.buf( R&R*RFtbuflocal_nd_varR)tbuf_ptrtbuffer_ptr_typet buf_ptr_typetinit_attributes(R!R&((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyt__init__s   cC@s1|j|_|j|_|j|_dS(N(tget_buf_shapevarstshapetget_buf_stridevarststridestget_buf_suboffsetvarst suboffsets(R!((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyRscC@s |jdS(Ns%s.diminfo[%d].suboffsets(t _for_all_ndim(R!((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyRscC@s |jdS(Ns%s.diminfo[%d].strides(R(R!((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyRscC@s |jdS(Ns%s.diminfo[%d].shape(R(R!((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyRscC@s0gt|jjD]}||j|f^qS(N(trangeR*R?R)(R!tsti((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyRscC@sg}|jj}|jj}|dkrxVt||j|jD]6\}}}|j||j||j|qIWd|} t} n|dkrd|} t} nJ|dkrd|} t } n+|dkrd|} t } n t st x=t||jD]&\}}|j||j|qW| |j jkr|j jj| |j d } |j d } | | | d | d |n|jj} d | | |jdj|f}|S(NRZs__Pyx_BufPtrFull%ddR]s__Pyx_BufPtrStrided%ddR^s__Pyx_BufPtrCContig%ddR_s__Pyx_BufPtrFortranContig%ddtutility_code_prototutility_code_defR%tnds%s(%s, %s, %s)s, (R*R?RWRlRRtappendtbuf_lookup_full_codetbuf_lookup_strided_codetbuf_lookup_c_codetbuf_lookup_fortran_codeRVRt globalstatet utility_codestaddRtempty_declaration_codeR~R(R!tcodet index_cnamestparamsRRWRRtotfuncnametfuncgent protocodetdefcodetbuf_ptr_type_codetptrcode((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytgenerate_buffer_lookup_codesD                   "   ( RTRURRRRRRR(((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyR|s      cC@sd}|j}|dkr(|d7}nW|dkrA|d7}n>|dkrZ|d7}n%|dkrs|d 7}n tst|jr|d 7}n|S( Nt PyBUF_FORMATRZs| PyBUF_INDIRECTR]s| PyBUF_STRIDESR^s| PyBUF_C_CONTIGUOUSR_s| PyBUF_F_CONTIGUOUSs| PyBUF_WRITABLE(RWRVRtwritable_needed(RFt buffer_typetflagsRW((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyt get_flags s           cC@s%|j}t|j_t|j_dS(N(RFR0R}R1t rcbuf_var(R&RF((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytused_buffer_aux_varss  c C@s|j|jj}}|jj}ddg}|dkrJ|jdng}xNt|jjD]:}x1|D])}|jd||||||fqpWqcW|jdj |dS(NRRRZRs1%s.diminfo[%d].%s = %s.rcbuffer->pybuffer.%s[%d];R ( RFR*RWR}R)RRR?tputlnR( t buf_entryRRFRWtpybuffernd_structtfldnamestlnRtfldname((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyt put_unpack_buffer_aux_into_scope"s      cC@so|j}|jj}|jj}|jd||jd||jd||jd||fdS(Ns%s.pybuffer.buf = NULL;s%s.refcount = 0;s%s.data = NULL;s%s.rcbuffer = &%s;(RFR}R)RR(R&RtbufauxRtpybuffer_struct((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyt put_init_vars4s   cC@s|jjt|j}t||j||j}|jd|jd|jjj |j|j d|||jdt ||dS(Nt{s'__Pyx_BufFmt_StackElem __pyx_stack[%d];s%s == -1t}( RRtacquire_utility_codeRFtget_getbuffer_callR)R*RR=tstruct_nesting_deptht error_goto_ifR(R&RR+RFt getbuffer((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytput_acquire_arg_bufferAs   cC@s.|jjt|jd|jjjdS(Ns0__Pyx_SafeReleaseBuffer(&%s.rcbuffer->pybuffer);(RRRRRFR}R)(RR&((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytput_release_buffer_codeOsc C@sP|j}t|j}t||}|jj}t||j}dtS(Ns__Pyx_GetBufferAndValidate(&%(pybuffernd_struct)s.rcbuffer->pybuffer, (PyObject*)%(obj_cname)s, &%(dtype_typeinfo)s, %(flags)s, %(ndim)d, %(cast)d, __pyx_stack)( R?RrRYRR}R)tget_type_information_cnameR=tlocals( Rt obj_cnameRFRR?RYRRtdtype_typeinfo((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyRSs  cC@s|j|j}}|jjt|jj}t||} |jd|jd|j j t |d||} |r9|jd||j j tjdt} |jd| | |f|jd|jd| gtd D]!} |j j tjdt^q\} }}|jd | ||f|jd|jd | ||jd | ||f|jjt|jd |jd|jd| ||fx'| ||fD]}|j j|qW|jd|jdt|||j|j| ||j j| n|jd|jd | ||jd|tj|tjd|f|j|j||jdt|||jd|jddS(sJ Generate code for reassigning a buffer variables. This only deals with getting the buffer auxiliary structure and variables set up correctly, the assignment itself and refcounting is the responsibility of the caller. However, the assignment operation may throw an exception so that the reassignment never happens. Depending on the circumstances there are two possible outcomes: - Old buffer released, new acquired, rhs assigned to lhs - Old buffer released, new acquired which fails, reaqcuire old lhs buffer (which may or may not succeed). Rs'__Pyx_BufFmt_StackElem __pyx_stack[%d];s%ss0__Pyx_SafeReleaseBuffer(&%s.rcbuffer->pybuffer);t manage_refs%s = %s;s if (%s) {s%s < 0isPyErr_Fetch(&%s, &%s, &%s);s%s == -1s/Py_XDECREF(%s); Py_XDECREF(%s); Py_XDECREF(%s);s!__Pyx_RaiseBufferFallbackError();s} else {sPyErr_Restore(%s, %s, %s);RsA%s = %s; __Pyx_INCREF(Py_None); %s.rcbuffer->pybuffer.buf = NULL;tPy_NoneN(RFR*RRRR}R)RRR=RRt funcstatet allocate_tempR t c_int_typeRVtunlikelyRtpy_object_typetraise_buffer_fallback_codet release_tempRterror_goto_if_negttypecastt error_gototput(t lhs_cnamet rhs_cnameRtis_initializedR+RRFRRRRt retcode_cnameRR*Rzttbtt((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytput_assign_to_buffer_sL  :"     "     cC@sb|do |}|dr|jjtjdt}|jd|x tt|||jD]\} \} } } | dkr|jd| |r|jd| | f|jd|j d | || fn|jd || f|j d n| dkrd } nd } |jd|j d| | | f|| fqeW|rs|j j t d}n|j j td}|jd|j d||jd||f|j|j||jd|jj|n^|rRxUt|||jD]8\} } } | dkr|jd| | | fqqWn|j||S(se Generates code to process indices and calculate an offset into a buffer. Returns a C string which gives a pointer which can be read from or written to at will (it is an expression so caller should store it in a temporary if it is used more than once). As the bounds checking can have any number of combinations of unsigned arguments, smart optimizations etc. we insert it directly in the function body. The lookup however is delegated to a inline function that is instantiated once per ndim (lookup with suboffsets tend to get quite complicated). entry is a BufferEntry t wraparoundt boundscheckRs%s = -1;is if (%s < 0) {s %s += %s;sif (%s) %s = %d;s%s < 0s%s = %d;s} else ts(size_t)s %s >= %s%st __Pyx_RaiseBufferIndexErrorNogilt__Pyx_RaiseBufferIndexErrors if (%s) {s%s != -1s%s(%s);Rsif (%s < 0) %s += %s;(RRR RRVRt enumerateRlRRRRRtraise_indexerror_nogiltraise_indexerror_codeRRR(R&t index_signedsRt directivesR+RRXtin_nogil_contexttfailed_dim_temptdimtsignedR)RRYtfunc((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytput_buffer_lookup_codesF 4       ( $cC@s|jtdS(N(Rtbuffer_struct_declare_code(Ru((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyRscC@sCtddjdg|djdg|f}td|S(Nsi static Py_ssize_t __Pyx_zeros[] = {%s}; static Py_ssize_t __Pyx_minusones[] = {%s}; s, t0s-1tproto(R RR(RR((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytget_empty_bufstruct_codes3c C@sdjgt|D]}d|||f^q}|jd||||fdjgt|D]}d|||f^qe}|jd||f|jtd||fdjgt|D]"}td||||f^qd d S( s Generates a buffer lookup function for the right number of dimensions. The function gives back a void* at the right location. s, s i%d, s%d, o%ds1#define %s(type, buf, %s) (type)(%s_imp(buf, %s))s.Py_ssize_t i%d, Py_ssize_t s%d, Py_ssize_t o%ds1static CYTHON_INLINE void* %s_imp(void* buf, %s);sf static CYTHON_INLINE void* %s_imp(void* buf, %s) { char* ptr = (char*)buf; RsX ptr += s%d * i%d; if (o%d >= 0) ptr = *((char**)ptr) + o%d; s return ptr; }N(RRRR (RtdefinR%RRt macroargstfuncargs((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyRs 55cC@sdjgt|D]}d||f^q}djgt|D]}d||f^qE}|jd|||fdS(s Generates a buffer lookup function for the right number of dimensions. The function gives back a void* at the right location. s, si%d, s%ds + s i%d * s%ds1#define %s(type, buf, %s) (type)((char*)buf + %s)N(RRR(RRR%RRtargstoffset((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyR s22cC@s|dkr |jd|ndjgt|D]}d||f^q3}djgt|dD]}d||f^qi}|jd||||dfdS( s Similar to strided lookup, but can assume that the last dimension doesn't need a multiplication as long as. Still we keep the same signature for now. is.#define %s(type, buf, i0, s0) ((type)buf + i0)s, si%d, s%ds + s i%d * s%ds9#define %s(type, buf, %s) ((type)((char*)buf + %s) + i%d)N(RRR(RRR%RRRR((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyRs  26cC@s|dkr |jd|ndjgt|D]}d||f^q3}djgtd|D]}d||f^qh}|jd|||dfd S( sB Like C lookup, but the first index is optimized instead. is.#define %s(type, buf, i0, s0) ((type)buf + i0)s, si%d, s%ds + s i%d * s%ds9#define %s(type, buf, %s) ((type)((char*)buf + %s) + i%d)iN(RRR(RRR%RRRR((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyR!s  25cC@s|jtdS(N(RtGetAndReleaseBufferUtilityCode(Ru((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyR-sRcB@sAeZdZeZdZdZdZdZ dZ RS(cC@sdS(N((R!((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyR8scC@s t|tS(N(RR(R!tother((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyt__eq__;scC@sdS(NiFos((R!((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyt__hash__>scK@sdS(N((R!tkwargs((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytget_treeAsc@s|d}|d}|jj}|jjgtfd|tjddddtd}|j|j }|j|j |j |d }|j ||j |dS( NRRc@s|krdSj|x|jD]}|q'Wx|jD]}t|jtrcqEn|j}|jrE|kr|j rqEnd}}xJ|j j D]<}|j dkr|j }q|j dkr|j }qqW|r j|j||fq qEqEWdS(Nu __getbuffer__u__releasebuffer__(Rtcimported_modulest type_entriesRR;RR*RpR1R.Rtpyfunc_entriesR%t func_cnameRt typeptr_cname(RtmteRtreleaseR`R(t cython_scopetfind_buffer_typesttypestvisited_scopes(s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyRLs(      tGetAndReleaseBuffert from_filesBuffer.ctcontextRi( t module_nodeRRRtsetRtloadtdictt format_codeRtinject_string_constantstimplR(R!toutputRt proto_codeRut util_codeRR ((RRRRs9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytput_codeCs         N( RTRUR.trequiresRVtis_cython_utilityRRRRR(((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyR0s    cC@sM|jr dS|jrdS|js,|jr5d}nd}||jSdS(Ntobjecttptrtnn_R(t is_pyobjectR>t is_typedeftis_struct_or_uniontspecialization_name(R=R2((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytmangle_dtype_namens   c C@s9t|}d|}d|}|jr-dS|dkrH|j}n|dkrctsctn||jjkr5|jjj||jd}g}|j rx)|j r|j |j |j }qWn|j o|j}|j} |jrd}n|jr|jj} t| dks3tg| D]} t|| j|d^q:} |jd|d txLt| | D];\} } |jd | | j|j| jfd tqW|jd d t|jd d tn tstt|}d }d }|tjkr4d| }d}n|jrTd| }d|}nr|sc|j rld}nZ|j!r~d}nH|jrd}|j"rd}qn$|j#rd}ntst|d} |||| dj$g|D]}t|^qpd t||||f }|j| |d tn|S(sr Output the run-time type information (__Pyx_TypeInfo) for given dtype, and return the name of the type info struct. Structs with two floats of the same size are encoded as complex numbers. One can seperate between complex numbers declared as struct or with native encoding by inspecting to see if the fields field of the type is filled in. s__Pyx_TypeInfo_%ss__Pyx_StructFields_%ssittypeinfotNULLis!static __Pyx_StructField %s[] = {tsafes {&%s, "%s", offsetof(%s, %s)},s {NULL, NULL, 0}s};RsIS_UNSIGNED(%s)s'H's%s ? 'U' : 'I's'C's'R's'S't__PYX_BUF_FLAGS_PACKED_STRUCTs'O'sLstatic __Pyx_TypeInfo %s = { "%s", %s, sizeof(%s), { %s }, %s, %s, %s, %s };s, N(%Rtis_errorR.RRVRRRRtis_arrayRtsizet base_typeRtcan_be_complexRtis_simple_buffer_dtypet is_structRt var_entriesR8RR*RR0RlR%R)tstrR t c_char_typetis_intt is_complextis_floattpackedRR(RR=tmaxdeptht namesuffixR%tstructinfo_namettypecodet arraysizestcomplex_possibletdeclcodetfieldstfRRtrepRt is_unsignedt typegroupRttup((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pyR|sv              ) *                4cK@s<|dkrtj|d|Stj|dd||SdS(NsBuffer.cR(R.RRR(tutil_code_nameRR((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytload_buffer_utilitys tmax_dimstBufferStructDeclareRtBufferIndexErrortBufferIndexErrorNogiltBufferFallbackErrortBufferFormatStructst proto_blocktutility_code_proto_before_typestBufferFormatCheckRtTypeInfoToFormattTypeInfoCompareN(sdtypesndimsmodesnegative_indicesscast(Et __future__RtVisitorRRtErrorsRRRtCodeRRRRR R R R RRjR0RVRfRhRkRiRmRoRtRsRqRbR.R{RR|RRRRRRRRRRRR@R RRRRRRRRR9RR%RRRRRtbuffer_structs_codeRt_typeinfo_to_format_codettypeinfo_compare_code(((s9/bar/jli/Chip-seq/script/cython/Cython/Compiler/Buffer.pytsz \"EG      H B    >  U