Rewrote tesseroid code for gzz in python using numpy fancy indexing. Now only 2x slower, not 20x