Cython:对于我来自1 <= i <N

Question

我正在学习Cython并遇到了这些代码:

import numpy as np
cimport numpy as np

def mean(np.ndarray[np.double_t] input):

    cdef np.double_t cur
    # Py_ssize_t is numpy's index type
    cdef Py_ssize_t i
    cdef Py_ssize_t N = len(input)

    for i from 0 <= i < N:
        cur += input[i]

    return cur / N

a=np.array([1,2,3,4], dtype=np.double)

显然,这将返回的平均一个是2.5.我的问题是:

for循环是Python循环,Cython还是C？

Answer 1

编译它并看到:Cython生成的C代码很好地注释.

  /* "cyexample.pyx":11
 *     cdef Py_ssize_t N = len(input)
 * 
 *     for i from 0 <= i < N:             # <<<<<<<<<<<<<<
 *         cur += input[i]
 * 
 */
  __pyx_t_1 = __pyx_v_N;
  for (__pyx_v_i = 0; __pyx_v_i < __pyx_t_1; __pyx_v_i++) {
    /* "cyexample.pyx":12
 * 
 *     for i from 0 <= i < N:
 *         cur += input[i]             # <<<<<<<<<<<<<<
 * 
 *     return cur / N
 */
    __pyx_t_2 = __pyx_v_i;
    __pyx_t_3 = -1;
    if (__pyx_t_2 < 0) {
      __pyx_t_2 += __pyx_bshape_0_input;
      if (unlikely(__pyx_t_2 < 0)) __pyx_t_3 = 0;
    } else if (unlikely(__pyx_t_2 >= __pyx_bshape_0_input)) __pyx_t_3 = 0;
    if (unlikely(__pyx_t_3 != -1)) {
      __Pyx_RaiseBufferIndexError(__pyx_t_3);
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    }
    __pyx_v_cur = (__pyx_v_cur + (*__Pyx_BufPtrStrided1d(__pyx_t_5numpy_double_t *, __pyx_bstruct_input.buf, __pyx_t_2, __pyx_bstride_0_input)));
  }

因此循环本身成功转换为C.注意,现在Cython可以自然地处理范围,因此不需要较旧的"从0 <= i <N"样式.引入(非Python)"for/from"语法的重点是表示哪些循环应该是C-ified.

Answer 2

for..from似乎是一个Pyrex/Cython循环:http://docs.cython.org/src/userguide/language_basics.html#integer-for-loops