Constrained version suffers from integer overflow

[KilianB]

Please see pull request: #2

awarp sometimes returns np.nan for calculations since values in the distance matrix [-1,-1] are negative and the square root fails. This is related to an uncaught overflow in numpy or numba