import numpy as np
nn = 15
print (" nn,", nn)
xvec = np.zeros(nn)
for i in range(nn): xvec[i] = i+1
print (" Vector:", *xvec)
import modVecsum
xsumf = modVecsum.vec_sum(xvec)
print (" sumF,", xsumf)
import ctypes
my_code = ctypes.CDLL("./vec_sum.so")
my_code.vec_sum.restype = ctypes.c_double
my_code.vec_sum.argtypes = \
[ctypes.POINTER(ctypes.c_double), ctypes.c_int]
xsumc = my_code.vec_sum \
(xvec.ctypes.data_as\
(ctypes.POINTER(ctypes.c_double)), nn)
print (" sumC,", xsumc)