docs/examples/tutorial/profiling_tutorial/calc_pi_3.pyx


1
2
3
4
5
6
7
8
9
10
11
12
13
14

# cython: profile=True


cdef inline double recip_square(long long i) except -1.0:
    return 1. / (i * i)

def approx_pi(int n=10000000):
    cdef double val = 0.
    cdef int k
    for k in range(1, n + 1):
        val += recip_square(k)
    return (6 * val) ** .5