ported to numpy with numba

[mikofski]

As indicated in #15 and #16, the calculation is too trivial to take advantage of parallelization, therefore it is a better candidate for SIMD vectorization.

• use numpy
• use numpy with Numba
• use numpy with Cython?
• use Intel MKL VML (vector math library)

easiest are first 2 options