Add pure Python implementation

Author: gjbex

source-code/cython/Numpy/.gitignore

@@ -2,3 +2,5 @@ build
 *.so
 array_init.c
 array_sum.c
+array_init_pure.c
+array_sum_pure.c

source-code/cython/Numpy/Makefile

@@ -1,11 +1,14 @@
 VERSION = cpython-311-x86_64-linux-gnu
 ARRAY_SUM_LIB = array_sum.$(VERSION).so
 ARRAY_INIT_LIB = array_init.$(VERSION).so
+ARRAY_SUM_PURE_LIB = array_sum_pure.$(VERSION).so
+ARRAY_INIT_PURE_LIB = array_init_pure.$(VERSION).so
 PROFILE = array_sums_prof.txt
 override DIM = 1000
 override ITER = 10
 
-all: $(ARRAY_SUM_LIB) $(ARRAY_INIT_LIB)
+all: $(ARRAY_SUM_LIB) $(ARRAY_SUM_PURE_LIB) \
+	 $(ARRAY_INIT_LIB) $(ARRAY_INIT_PURE_LIB)
 
 $(ARRAY_SUM_LIB): array_sum.pyx
 	python setup.py build_ext --inplace
@@ -19,5 +22,4 @@ profile: $(ARRAY_SUM_LIB) compute_sums.py
 
 clean:
 	python setup.py clean
-	rm -f array_sum.c $(ARRAY_SUM_LIB) $(PROFILE) \
-	rm -f array_init.c $(ARRAY_INIT_LIB)
+	$(RM) $(wildcard *.c) $(wildcard *.so)

source-code/cython/Numpy/array_init_pure.py

@@ -0,0 +1,9 @@
+import cython
+import numpy as np
+
+def cons(dims: tuple[int]) -> np.array:
+    array = np.empty(np.prod(dims))
+    mv: cython.double[:] = array
+    for i in range(mv.nbytes//mv.itemsize):
+        array[i] = np.float64(i)
+    return array.reshape(dims)

source-code/cython/Numpy/array_sum_pure.py

@@ -0,0 +1,12 @@
+import cython
+import numpy as np
+
+def array_sum(array: np.array) -> np.float64:
+    shape = array.shape
+    array.reshape((-1, ))
+    mv: cython.double[:] = array.reshape((-1, ))
+    total: cython.double = 0.0
+    with cython.boundscheck(False):
+        for i in range(mv.shape[0]):
+            total += mv[i]
+    return total

source-code/cython/Numpy/setup.py

@@ -6,3 +6,8 @@
 setup(
     ext_modules=cythonize('*.pyx')
 )
+setup(
+    ext_modules=cythonize(['array_init.pyx', 'array_init_pure.py',
+                           'array_sum.pyx', 'array_sum_pure.py'],
+                          language_level='3str')
+)