blob: 415ac4a22d0894f5c1a61568a727eb6ae99c9546 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
# cython: profile=True
cimport cython
@cython.profile(False)
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
|
