El código de morton 2D codifica / decodifica 64 bits

¿Cómo codificar / decodificar códigos de mortón (orden z) dados [x, y] como enteros sin signo de 32 bits que producen código de mortón de 64 bits y viceversa? Tengo xy2d y d2xy, pero solo para coordenadas de 16 bits de ancho que producen un número de morton de 32 bits. Busqué mucho en la red, pero no pude encontrar. Por favor ayuda.