From 081e3308453610e174223d6cd46b2f11ff4d634a Mon Sep 17 00:00:00 2001
From: Eric Fischer <enf@pobox.com>
Date: Wed, 12 Oct 2016 09:49:25 -0700
Subject: [PATCH] Reset counter with each tile, to match old low-zoom feature
 density

---
 main.cpp | 25 +++++++++++++++++++++++++
 tile.cpp |  1 -
 2 files changed, 25 insertions(+), 1 deletion(-)

diff --git a/main.cpp b/main.cpp
index 0814d1f..853fd54 100644
--- a/main.cpp
+++ b/main.cpp
@@ -217,15 +217,38 @@ struct drop_state {
 	double scale;
 	double seq;
 	long long included;
+	unsigned x;
+	unsigned y;
 };
 
 int calc_feature_minzoom(struct index *ix, struct drop_state *ds, int maxzoom, int basezoom, double droprate, double gamma) {
 	int feature_minzoom = 0;
+	unsigned xx, yy;
+	decode(ix->index, &xx, &yy);
 
 	if (gamma >= 0 && (ix->t == VT_POINT ||
 			   (additional[A_LINE_DROP] && ix->t == VT_LINE) ||
 			   (additional[A_POLYGON_DROP] && ix->t == VT_POLYGON))) {
 		for (ssize_t i = maxzoom; i >= 0; i--) {
+			// XXX This resets the feature counter at the start of each tile,
+			// which makes the feature count come out close to what it is if
+			// feature dropping happens during tiling. It means that the low
+			// zooms are heavier than they legitimately should be though.
+			{
+				unsigned xxx = 0, yyy = 0;
+				if (i != 0) {
+					xxx = xx >> (32 - i);
+					yyy = yy >> (32 - i);
+				}
+				if (ds[i].x != xxx || ds[i].y != yyy) {
+					ds[i].seq = 0;
+					ds[i].gap = 0;
+					ds[i].previndex = 0;
+				}
+				ds[i].x = xxx;
+				ds[i].y = yyy;
+			}
+
 			ds[i].seq++;
 		}
 		for (ssize_t i = maxzoom; i >= 0; i--) {
@@ -836,6 +859,8 @@ void prep_drop_states(struct drop_state *ds, int maxzoom, int basezoom, double d
 		ds[i].scale = (double) (1LL << (64 - 2 * (i + 8)));
 		ds[i].seq = 0;
 		ds[i].included = 0;
+		ds[i].x = 0;
+		ds[i].y = 0;
 	}
 }
 
diff --git a/tile.cpp b/tile.cpp
index f705533..599a190 100644
--- a/tile.cpp
+++ b/tile.cpp
@@ -596,7 +596,6 @@ long long write_tile(FILE *geoms, long long *geompos_in, char *metabase, char *s
 		double accum_area = 0;
 
 		double interval = 0;
-		double seq = 0;
 		if (z < basezoom) {
 			interval = std::exp(std::log(droprate) * (basezoom - z));
 		}