2012-03-21 14:53:57 +08:00
|
|
|
/* Copyright (c) 2010 Scott Lembcke
|
|
|
|
*
|
|
|
|
* Permission is hereby granted, free of charge, to any person obtaining a copy
|
|
|
|
* of this software and associated documentation files (the "Software"), to deal
|
|
|
|
* in the Software without restriction, including without limitation the rights
|
|
|
|
* to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
|
|
|
|
* copies of the Software, and to permit persons to whom the Software is
|
|
|
|
* furnished to do so, subject to the following conditions:
|
|
|
|
*
|
|
|
|
* The above copyright notice and this permission notice shall be included in
|
|
|
|
* all copies or substantial portions of the Software.
|
|
|
|
*
|
|
|
|
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
|
|
|
|
* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
|
|
|
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
|
|
|
|
* AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
|
|
|
|
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
|
|
|
|
* OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
|
|
|
|
* SOFTWARE.
|
|
|
|
*/
|
|
|
|
|
|
|
|
#include "chipmunk_private.h"
|
|
|
|
|
|
|
|
static inline cpSpatialIndexClass *Klass();
|
|
|
|
|
|
|
|
//MARK: Basic Structures
|
|
|
|
|
|
|
|
typedef struct Bounds {
|
2012-06-12 09:38:38 +08:00
|
|
|
cpFloat min, max;
|
2012-03-21 14:53:57 +08:00
|
|
|
} Bounds;
|
|
|
|
|
|
|
|
typedef struct TableCell {
|
2012-06-12 09:38:38 +08:00
|
|
|
void *obj;
|
|
|
|
Bounds bounds;
|
2012-03-21 14:53:57 +08:00
|
|
|
} TableCell;
|
|
|
|
|
|
|
|
struct cpSweep1D
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
cpSpatialIndex spatialIndex;
|
|
|
|
|
|
|
|
int num;
|
|
|
|
int max;
|
|
|
|
TableCell *table;
|
2012-03-21 14:53:57 +08:00
|
|
|
};
|
|
|
|
|
|
|
|
static inline cpBool
|
|
|
|
BoundsOverlap(Bounds a, Bounds b)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
return (a.min <= b.max && b.min <= a.max);
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static inline Bounds
|
|
|
|
BBToBounds(cpSweep1D *sweep, cpBB bb)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
Bounds bounds = {bb.l, bb.r};
|
|
|
|
return bounds;
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static inline TableCell
|
|
|
|
MakeTableCell(cpSweep1D *sweep, void *obj)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
TableCell cell = {obj, BBToBounds(sweep, sweep->spatialIndex.bbfunc(obj))};
|
|
|
|
return cell;
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
//MARK: Memory Management Functions
|
|
|
|
|
|
|
|
cpSweep1D *
|
|
|
|
cpSweep1DAlloc(void)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
return (cpSweep1D *)cpcalloc(1, sizeof(cpSweep1D));
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static void
|
|
|
|
ResizeTable(cpSweep1D *sweep, int size)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
sweep->max = size;
|
|
|
|
sweep->table = (TableCell *)cprealloc(sweep->table, size*sizeof(TableCell));
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
cpSpatialIndex *
|
|
|
|
cpSweep1DInit(cpSweep1D *sweep, cpSpatialIndexBBFunc bbfunc, cpSpatialIndex *staticIndex)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
cpSpatialIndexInit((cpSpatialIndex *)sweep, Klass(), bbfunc, staticIndex);
|
|
|
|
|
|
|
|
sweep->num = 0;
|
|
|
|
ResizeTable(sweep, 32);
|
|
|
|
|
|
|
|
return (cpSpatialIndex *)sweep;
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
cpSpatialIndex *
|
|
|
|
cpSweep1DNew(cpSpatialIndexBBFunc bbfunc, cpSpatialIndex *staticIndex)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
return cpSweep1DInit(cpSweep1DAlloc(), bbfunc, staticIndex);
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static void
|
|
|
|
cpSweep1DDestroy(cpSweep1D *sweep)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
cpfree(sweep->table);
|
|
|
|
sweep->table = NULL;
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
//MARK: Misc
|
|
|
|
|
|
|
|
static int
|
|
|
|
cpSweep1DCount(cpSweep1D *sweep)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
return sweep->num;
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static void
|
|
|
|
cpSweep1DEach(cpSweep1D *sweep, cpSpatialIndexIteratorFunc func, void *data)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
TableCell *table = sweep->table;
|
|
|
|
for(int i=0, count=sweep->num; i<count; i++) func(table[i].obj, data);
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static int
|
|
|
|
cpSweep1DContains(cpSweep1D *sweep, void *obj, cpHashValue hashid)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
TableCell *table = sweep->table;
|
|
|
|
for(int i=0, count=sweep->num; i<count; i++){
|
|
|
|
if(table[i].obj == obj) return cpTrue;
|
|
|
|
}
|
|
|
|
|
|
|
|
return cpFalse;
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
//MARK: Basic Operations
|
|
|
|
|
|
|
|
static void
|
|
|
|
cpSweep1DInsert(cpSweep1D *sweep, void *obj, cpHashValue hashid)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
if(sweep->num == sweep->max) ResizeTable(sweep, sweep->max*2);
|
|
|
|
|
|
|
|
sweep->table[sweep->num] = MakeTableCell(sweep, obj);
|
|
|
|
sweep->num++;
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static void
|
|
|
|
cpSweep1DRemove(cpSweep1D *sweep, void *obj, cpHashValue hashid)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
TableCell *table = sweep->table;
|
|
|
|
for(int i=0, count=sweep->num; i<count; i++){
|
|
|
|
if(table[i].obj == obj){
|
|
|
|
int num = --sweep->num;
|
|
|
|
|
|
|
|
table[i] = table[num];
|
|
|
|
table[num].obj = NULL;
|
|
|
|
|
|
|
|
return;
|
|
|
|
}
|
|
|
|
}
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
//MARK: Reindexing Functions
|
|
|
|
|
|
|
|
static void
|
|
|
|
cpSweep1DReindexObject(cpSweep1D *sweep, void *obj, cpHashValue hashid)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
// Nothing to do here
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static void
|
|
|
|
cpSweep1DReindex(cpSweep1D *sweep)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
// Nothing to do here
|
|
|
|
// Could perform a sort, but queries are not accelerated anyway.
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
//MARK: Query Functions
|
|
|
|
|
|
|
|
static void
|
|
|
|
cpSweep1DQuery(cpSweep1D *sweep, void *obj, cpBB bb, cpSpatialIndexQueryFunc func, void *data)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
// Implementing binary search here would allow you to find an upper limit
|
|
|
|
// but not a lower limit. Probably not worth the hassle.
|
|
|
|
|
|
|
|
Bounds bounds = BBToBounds(sweep, bb);
|
|
|
|
|
|
|
|
TableCell *table = sweep->table;
|
|
|
|
for(int i=0, count=sweep->num; i<count; i++){
|
|
|
|
TableCell cell = table[i];
|
2014-03-26 01:37:35 +08:00
|
|
|
if(BoundsOverlap(bounds, cell.bounds) && obj != cell.obj) func(obj, cell.obj, 0, data);
|
2012-06-12 09:38:38 +08:00
|
|
|
}
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static void
|
|
|
|
cpSweep1DSegmentQuery(cpSweep1D *sweep, void *obj, cpVect a, cpVect b, cpFloat t_exit, cpSpatialIndexSegmentQueryFunc func, void *data)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
cpBB bb = cpBBExpand(cpBBNew(a.x, a.y, a.x, a.y), b);
|
|
|
|
Bounds bounds = BBToBounds(sweep, bb);
|
|
|
|
|
|
|
|
TableCell *table = sweep->table;
|
|
|
|
for(int i=0, count=sweep->num; i<count; i++){
|
|
|
|
TableCell cell = table[i];
|
|
|
|
if(BoundsOverlap(bounds, cell.bounds)) func(obj, cell.obj, data);
|
|
|
|
}
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
//MARK: Reindex/Query
|
|
|
|
|
|
|
|
static int
|
|
|
|
TableSort(TableCell *a, TableCell *b)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
return (a->bounds.min < b->bounds.min ? -1 : (a->bounds.min > b->bounds.min ? 1 : 0));
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static void
|
|
|
|
cpSweep1DReindexQuery(cpSweep1D *sweep, cpSpatialIndexQueryFunc func, void *data)
|
|
|
|
{
|
2012-06-12 09:38:38 +08:00
|
|
|
TableCell *table = sweep->table;
|
|
|
|
int count = sweep->num;
|
|
|
|
|
|
|
|
// Update bounds and sort
|
|
|
|
for(int i=0; i<count; i++) table[i] = MakeTableCell(sweep, table[i].obj);
|
|
|
|
qsort(table, count, sizeof(TableCell), (int (*)(const void *, const void *))TableSort); // TODO use insertion sort instead
|
|
|
|
|
|
|
|
for(int i=0; i<count; i++){
|
|
|
|
TableCell cell = table[i];
|
|
|
|
cpFloat max = cell.bounds.max;
|
|
|
|
|
|
|
|
for(int j=i+1; table[j].bounds.min < max && j<count; j++){
|
2014-03-26 01:37:35 +08:00
|
|
|
func(cell.obj, table[j].obj, 0, data);
|
2012-06-12 09:38:38 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
// Reindex query is also responsible for colliding against the static index.
|
|
|
|
// Fortunately there is a helper function for that.
|
|
|
|
cpSpatialIndexCollideStatic((cpSpatialIndex *)sweep, sweep->spatialIndex.staticIndex, func, data);
|
2012-03-21 14:53:57 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
static cpSpatialIndexClass klass = {
|
2012-06-12 09:38:38 +08:00
|
|
|
(cpSpatialIndexDestroyImpl)cpSweep1DDestroy,
|
|
|
|
|
|
|
|
(cpSpatialIndexCountImpl)cpSweep1DCount,
|
|
|
|
(cpSpatialIndexEachImpl)cpSweep1DEach,
|
|
|
|
(cpSpatialIndexContainsImpl)cpSweep1DContains,
|
|
|
|
|
|
|
|
(cpSpatialIndexInsertImpl)cpSweep1DInsert,
|
|
|
|
(cpSpatialIndexRemoveImpl)cpSweep1DRemove,
|
|
|
|
|
|
|
|
(cpSpatialIndexReindexImpl)cpSweep1DReindex,
|
|
|
|
(cpSpatialIndexReindexObjectImpl)cpSweep1DReindexObject,
|
|
|
|
(cpSpatialIndexReindexQueryImpl)cpSweep1DReindexQuery,
|
|
|
|
|
|
|
|
(cpSpatialIndexQueryImpl)cpSweep1DQuery,
|
|
|
|
(cpSpatialIndexSegmentQueryImpl)cpSweep1DSegmentQuery,
|
2012-03-21 14:53:57 +08:00
|
|
|
};
|
|
|
|
|
|
|
|
static inline cpSpatialIndexClass *Klass(){return &klass;}
|
|
|
|
|