Add code to use GSL for geodesic generation. Much more general and roughly similar speed to prev. fn.