import ctypes
intersect_dso = ctypes.CDLL('./intersect_dso.so')
intersect_c = intersect_dso.intersect
intersect_c.argtypes = ([ctypes.c_double] * 6) + ([ctypes.POINTER(ctypes.c_double)] * 2)
def intersect(lat1, lon1, az1, lat2, lon2, az2):
dX = ctypes.c_double()
dY = ctypes.c_double()
intersect_c(
lat1, lon1, az1, lat2, lon2, az2,
ctypes.byref(dX), ctypes.byref(dY),
)
return dX.value, dY.value