cython

1
import random
2
from libc.stdlib cimport malloc, free
3

4
def random_noise(int number=1):
5
    cdef int i
6
    # allocate number * sizeof(double) bytes of memory
7
    cdef double *my_array = <double *> malloc(
8
        number * sizeof(double))
9
    if not my_array:
10
        raise MemoryError()
11

12
    try:
13
        ran = random.normalvariate
14
        for i in range(number):
15
            my_array[i] = ran(0, 1)
16

17
        # ... let's just assume we do some more heavy C calculations here to make up
18
        # for the work that it takes to pack the C double values into Python float
19
        # objects below, right after throwing away the existing objects above.
20

21
        return [x for x in my_array[:number]]
22
    finally:
23
        # return the previously allocated memory to the system
24
        free(my_array)
25