Как реализовать сортировку координат по оси z для массовой вставки Quad-Tree — C++

Это непроверенный код, вы должны дважды проверить, работает ли он.

Кроме того, этот код почти наверняка не является переносимым и может даже иметь неопределенное поведение. Это, безусловно, поведение, определяемое реализацией, по крайней мере, но оно, вероятно, не указано... Мне нужно более внимательно прочитать правила, касающиеся reinterpret_cast в и из char*, чтобы знать наверняка.

#include <cstdint>
#include <vector>

uint64_t reinterpretDoubleAsUInt(double d) {
  int const doubleSize = sizeof(double);
  char* array = reinterpret_cast<char*>(&d);
  uint64_t result = 0;
  for (auto i = 0; i < doubleSize; ++i) {
    result += (uint64_t)array[i] << (8*i);
  }

  return result;
}

bool lessThanZOrderDouble(std::vector<double> const& a, std::vector<double> const& b) {
  uint64_t j = 0;
  uint64_t x = 0;
  if (a.size() != b.size() || a.size() == 0) {
    throw std::exception();
  }

  int dimensions = a.size();

  for (auto i = 0; i < dimensions; ++i) {
    auto y = reinterpretDoubleAsUInt(a[i]) ^ reinterpretDoubleAsUInt(b[i]);
    if (x < y && x < (x ^ y)) {
      j = i;
      x = y;
    }
  }
  return (a[j] - b[j]) > 0;
}

int main() {
  // blank for compilation's sake
}