Accelerate the rowwise AVL implementation of distance covariance using the GPU

[vnmabus]

Numba supports GPU programming, but most numpy functions are unsupported, which makes almost impossible sharing code for the GPU and CPU implementations.

If, however, any charitable soul wants to try to implement an alternative version of rowwise using the AVL implementation, accelerated using the GPU via Numba, it would be very helpful.