summaryrefslogtreecommitdiff
path: root/libphysics/blockmap.c
blob: 7f02b1a4ae0697a4b46c403a5fd42694e7128731 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
#include <vecmath.h>
#include <physics.h>

#include <stdlib.h>
#include <string.h>

#include "dat.h"

#define BUCKETS      0x10000 /* must be a power of 2 */     

#define BUCKETMASK  (BUCKETS - 1)

#define PHI   1.618033988749895 
#define HMUL  (int)(BUCKETS / PHI)

static int hash(int a);
static void makeblockmap(Body *b);

static BlockmapNode *allocnode(void);

static BlockmapNode nodes[(1 << 20)];
static int nodesi;
static BlockmapNode *lists[BUCKETS];

void
phxmakeblkmap()
{
	nodesi = 0;

	for(int i = 0; i < BUCKETS; i++)
		lists[i] = NULL;

	for(int i = 0; i < phxbodypoolsize; i++) {
		Body *b = phxbodypool + i;
		if(!b->active)
			continue;
		makeblockmap(b);
	}
}

void
makeblockmap(Body *b)
{
	BlockmapNode *n;
	if((n = allocnode()) == NULL)
		return;

	n->id = b - phxbodypool;
	n->next = NULL;

	int mapx = (int)floorf(b->pos[0] / BLOCK_SIZE);
	int mapy = (int)floorf(b->pos[1] / BLOCK_SIZE);
	int h = hash(mapx + hash(mapy));

	n->next = lists[h];
	lists[h] = n;
}

int
hash(int a)
{
	return (a * HMUL) & BUCKETMASK;
}

BlockmapNode *
allocnode(void)
{
	if(nodesi == (sizeof(nodes) / sizeof(nodes[0])))
		return NULL;
	return nodes + nodesi++;
}

BlockmapNode *
phxnodelist(float x, float y)
{
	int mapx = (int)floorf(x / BLOCK_SIZE);
	int mapy = (int)floorf(y / BLOCK_SIZE);
	int h = hash(mapx + hash(mapy));
	return lists[h];
}