Next step of index widening; tests still pass

This commit is contained in:
Erica Fischer
2024-04-05 12:44:29 -07:00
parent 47124087ca
commit 77f5e9677f
5 changed files with 26 additions and 26 deletions
+17 -17
View File
@@ -9,8 +9,8 @@
#define UINT_BITS 32
index_t (*encode_index)(unsigned int wx, unsigned int wy) = NULL;
void (*decode_index)(index_t index, unsigned *wx, unsigned *wy) = NULL;
index_t (*encode_index)(unsigned long long wx, unsigned long long wy) = NULL;
void (*decode_index)(index_t index, unsigned long long *wx, unsigned long long *wy) = NULL;
struct projection projections[] = {
{"EPSG:4326", lonlat2tile, tile2lonlat, "urn:ogc:def:crs:OGC:1.3:CRS84"},
@@ -107,7 +107,7 @@ void tiletoepsg3857(long long ix, long long iy, int zoom, double *ox, double *oy
// https://en.wikipedia.org/wiki/Hilbert_curve
void hilbert_rot(unsigned long long n, unsigned *x, unsigned *y, unsigned long long rx, unsigned long long ry) {
void hilbert_rot(unsigned long long n, unsigned long long *x, unsigned long long *y, unsigned long long rx, unsigned long long ry) {
if (ry == 0) {
if (rx == 1) {
*x = n - 1 - *x;
@@ -120,11 +120,11 @@ void hilbert_rot(unsigned long long n, unsigned *x, unsigned *y, unsigned long l
}
}
unsigned long long hilbert_xy2d(unsigned long long n, unsigned x, unsigned y) {
unsigned long long d = 0;
index_t hilbert_xy2d(unsigned long long n, unsigned long long x, unsigned long long y) {
index_t d = 0;
unsigned long long rx, ry;
for (unsigned long long s = n / 2; s > 0; s /= 2) {
for (index_t s = n / 2; s > 0; s /= 2) {
rx = (x & s) != 0;
ry = (y & s) != 0;
@@ -135,12 +135,12 @@ unsigned long long hilbert_xy2d(unsigned long long n, unsigned x, unsigned y) {
return d;
}
void hilbert_d2xy(unsigned long long n, unsigned long long d, unsigned *x, unsigned *y) {
void hilbert_d2xy(index_t n, index_t d, unsigned long long *x, unsigned long long *y) {
unsigned long long rx, ry;
unsigned long long t = d;
index_t t = d;
*x = *y = 0;
for (unsigned long long s = 1; s < n; s *= 2) {
for (index_t s = 1; s < n; s *= 2) {
rx = 1 & (t / 2);
ry = 1 & (t ^ rx);
hilbert_rot(s, x, y, rx, ry);
@@ -150,20 +150,20 @@ void hilbert_d2xy(unsigned long long n, unsigned long long d, unsigned *x, unsig
}
}
index_t encode_hilbert(unsigned int wx, unsigned int wy) {
index_t encode_hilbert(unsigned long long wx, unsigned long long wy) {
return hilbert_xy2d(1LL << UINT_BITS, wx, wy);
}
void decode_hilbert(index_t index, unsigned *wx, unsigned *wy) {
void decode_hilbert(index_t index, unsigned long long *wx, unsigned long long *wy) {
hilbert_d2xy(1LL << UINT_BITS, index, wx, wy);
}
index_t encode_quadkey(unsigned int wx, unsigned int wy) {
unsigned long long out = 0;
index_t encode_quadkey(unsigned long long wx, unsigned long long wy) {
index_t out = 0;
int i;
for (i = 0; i < UINT_BITS; i++) {
unsigned long long v = ((wx >> (UINT_BITS - (i + 1))) & 1) << 1;
index_t v = ((wx >> (UINT_BITS - (i + 1))) & 1) << 1;
v |= (wy >> (UINT_BITS - (i + 1))) & 1;
v = v << (64 - 2 * (i + 1));
@@ -176,7 +176,7 @@ index_t encode_quadkey(unsigned int wx, unsigned int wy) {
static std::atomic<unsigned char> decodex[256];
static std::atomic<unsigned char> decodey[256];
void decode_quadkey(index_t index, unsigned *wx, unsigned *wy) {
void decode_quadkey(index_t index, unsigned long long *wx, unsigned long long *wy) {
static std::atomic<int> initialized(0);
if (!initialized) {
for (size_t ix = 0; ix < 256; ix++) {
@@ -197,8 +197,8 @@ void decode_quadkey(index_t index, unsigned *wx, unsigned *wy) {
*wx = *wy = 0;
for (size_t i = 0; i < 8; i++) {
*wx |= ((unsigned) decodex[(index >> (8 * i)) & 0xFF]) << (4 * i);
*wy |= ((unsigned) decodey[(index >> (8 * i)) & 0xFF]) << (4 * i);
*wx |= ((unsigned long long) decodex[(index >> (8 * i)) & 0xFF]) << (4 * i);
*wy |= ((unsigned long long) decodey[(index >> (8 * i)) & 0xFF]) << (4 * i);
}
}