1662 GistEntryVector *entryvec = (GistEntryVector *) PG_GETARG_POINTER(0);
1663 GIST_SPLITVEC *v = (GIST_SPLITVEC *) PG_GETARG_POINTER(1);
1677 POSTGIS_DEBUG(3,
"[GIST] 'picksplit' entered");
1681 maxoff = entryvec->n - 1;
1682 nentries = context.
entriesCount = maxoff - FirstOffsetNumber + 1;
1691 for (i = FirstOffsetNumber; i <= maxoff; i = OffsetNumberNext(i))
1693 box = (BOX2DF *) DatumGetPointer(entryvec->vector[i].key);
1694 if (i == FirstOffsetNumber)
1706 context.
first =
true;
1707 for (dim = 0; dim < 2; dim++)
1715 for (i = FirstOffsetNumber; i <= maxoff; i = OffsetNumberNext(i))
1717 box = (BOX2DF *) DatumGetPointer(entryvec->vector[i].key);
1720 intervalsLower[i - FirstOffsetNumber].
lower = box->xmin;
1721 intervalsLower[i - FirstOffsetNumber].
upper = box->xmax;
1725 intervalsLower[i - FirstOffsetNumber].
lower = box->ymin;
1726 intervalsLower[i - FirstOffsetNumber].
upper = box->ymax;
1734 memcpy(intervalsUpper, intervalsLower,
1777 rightLower = intervalsLower[i1].
lower;
1778 leftUpper = intervalsUpper[i2].
lower;
1784 while (i1 < nentries && (rightLower == intervalsLower[i1].lower ||
1785 isnan(intervalsLower[i1].lower)))
1787 leftUpper = Max(leftUpper, intervalsLower[i1].upper);
1792 rightLower = intervalsLower[i1].
lower;
1798 while (i2 < nentries && intervalsUpper[i2].upper <= leftUpper)
1813 rightLower = intervalsLower[i1].
upper;
1814 leftUpper = intervalsUpper[i2].
upper;
1820 while (i2 >= 0 && (leftUpper == intervalsUpper[i2].upper ||
1821 isnan(intervalsUpper[i2].upper)))
1823 rightLower = Min(rightLower, intervalsUpper[i2].lower);
1828 leftUpper = intervalsUpper[i2].
upper;
1834 while (i1 >= 0 && intervalsLower[i1].lower >= rightLower)
1841 rightLower, i1 + 1, leftUpper, i2 + 1);
1850 POSTGIS_DEBUG(4,
"no acceptable splits, trivial split");
1852 PG_RETURN_POINTER(v);
1863 POSTGIS_DEBUGF(4,
"split direction: %d", context.
dim);
1866 v->spl_left = (OffsetNumber *) palloc(nentries *
sizeof(OffsetNumber));
1867 v->spl_right = (OffsetNumber *) palloc(nentries *
sizeof(OffsetNumber));
1872 leftBox = palloc0(
sizeof(BOX2DF));
1873 rightBox = palloc0(
sizeof(BOX2DF));
1879 commonEntriesCount = 0;
1883 #define PLACE_LEFT(box, off) \
1885 if (v->spl_nleft > 0) \
1886 adjustBox(leftBox, box); \
1888 *leftBox = *(box); \
1889 v->spl_left[v->spl_nleft++] = off; \
1892 #define PLACE_RIGHT(box, off) \
1894 if (v->spl_nright > 0) \
1895 adjustBox(rightBox, box); \
1897 *rightBox = *(box); \
1898 v->spl_right[v->spl_nright++] = off; \
1905 for (i = FirstOffsetNumber; i <= maxoff; i = OffsetNumberNext(i))
1913 box = (BOX2DF *) DatumGetPointer(entryvec->vector[i].key);
1914 if (context.
dim == 0)
1925 if (upper <= context.
leftUpper || isnan(upper))
1928 if (lower >= context.
rightLower || isnan(lower))
1931 commonEntries[commonEntriesCount++].
index = i;
1959 if (commonEntriesCount > 0)
1971 for (i = 0; i < commonEntriesCount; i++)
1973 box = (BOX2DF *) DatumGetPointer(entryvec->vector[
1974 commonEntries[i].index].key);
1987 for (i = 0; i < commonEntriesCount; i++)
1989 box = (BOX2DF *) DatumGetPointer(entryvec->vector[
1990 commonEntries[i].index].key);
1996 if (v->spl_nleft + (commonEntriesCount - i) <= m)
1998 else if (v->spl_nright + (commonEntriesCount - i) <= m)
2010 v->spl_ldatum = PointerGetDatum(leftBox);
2011 v->spl_rdatum = PointerGetDatum(rightBox);
2013 POSTGIS_DEBUG(4,
"[GIST] 'picksplit' completed");
2015 PG_RETURN_POINTER(v);
static char * box2df_to_string(const BOX2DF *a)
static void fallbackSplit(GistEntryVector *entryvec, GIST_SPLITVEC *v)
static int interval_cmp_upper(const void *i1, const void *i2)
#define PLACE_RIGHT(box, off)
static void g_box_consider_split(ConsiderSplitContext *context, int dimNum, float rightLower, int minLeftCount, float leftUpper, int maxLeftCount)
static int interval_cmp_lower(const void *i1, const void *i2)
static int common_entry_cmp(const void *i1, const void *i2)
#define PLACE_LEFT(box, off)
static float box2df_penalty(const BOX2DF *b1, const BOX2DF *b2)
static void adjustBox(BOX2DF *b, BOX2DF *addon)