pub fn morton_encode_u64(x: u32, y: u32, z: u32) -> u64
Encode x,y,z position into a u64 morton value. Input should be 0..=2u32.pow(21) (or 1u32 << 21)