Consider speeding it up with raw numpy and not function invocations
Consider speeding it up with raw numpy and not function invocations