nuclear@6: #include <stdlib.h>
nuclear@6: #include "pnoise.h"
nuclear@6: #include "vmath/vmath.h"
nuclear@6: 
nuclear@6: #define B	0x100
nuclear@6: #define BM	0xff
nuclear@6: #define N	0x1000
nuclear@6: #define NP	12
nuclear@6: #define NM	0xfff
nuclear@6: 
nuclear@6: #define s_curve(t)	((t) * (t) * (3.0f - 2.0f * (t)))
nuclear@6: #define setup(elem, b0, b1, r0, r1) \
nuclear@6: 	do {							\
nuclear@6: 		float t = elem + N;		\
nuclear@6: 		b0 = ((int)t) & BM;			\
nuclear@6: 		b1 = (b0 + 1) & BM;			\
nuclear@6: 		r0 = t - (int)t;			\
nuclear@6: 		r1 = r0 - 1.0f;				\
nuclear@6: 	} while(0)
nuclear@6: 
nuclear@6: #define setup_p(elem, b0, b1, r0, r1, p) \
nuclear@6: 	do {							\
nuclear@6: 		float t = elem + N;		\
nuclear@6: 		b0 = (((int)t) & BM) % p;	\
nuclear@6: 		b1 = ((b0 + 1) & BM) % p;	\
nuclear@6: 		r0 = t - (int)t;			\
nuclear@6: 		r1 = r0 - 1.0f;				\
nuclear@6: 	} while(0)
nuclear@6: 
nuclear@6: static int perm[B + B + 2];
nuclear@6: static vec2_t grad2[B + B + 2];
nuclear@6: static bool tables_valid;
nuclear@6: 
nuclear@6: static void init_noise()
nuclear@6: {
nuclear@6: 	for(int i=0; i<B; i++) {
nuclear@6: 		perm[i] = i;
nuclear@6: 
nuclear@6: 		grad2[i].x = (float)((rand() % (B + B)) - B) / B;
nuclear@6: 		grad2[i].y = (float)((rand() % (B + B)) - B) / B;
nuclear@6: 		grad2[i] = v2_normalize(grad2[i]);
nuclear@6: 	}
nuclear@6: 
nuclear@6: 	for(int i=0; i<B; i++) {
nuclear@6: 		int rand_idx = rand() % B;
nuclear@6: 
nuclear@6: 		int tmp = perm[i];
nuclear@6: 		perm[i] = perm[rand_idx];
nuclear@6: 		perm[rand_idx] = tmp;
nuclear@6: 	}
nuclear@6: 
nuclear@6: 	for(int i=0; i<B+2; i++) {
nuclear@6: 		perm[B + i] = perm[i];
nuclear@6: 		grad2[B + i] = grad2[i];
nuclear@6: 	}
nuclear@6: }
nuclear@6: 
nuclear@6: #define lerp(a, b, t)	((a) + ((b) - (a)) * t)
nuclear@6: 
nuclear@6: float dbg_noise2(float x, float y)
nuclear@6: {
nuclear@6: 	if(!tables_valid) {
nuclear@6: 		init_noise();
nuclear@6: 		tables_valid = true;
nuclear@6: 	}
nuclear@6: 
nuclear@6: 	int bx0, bx1, by0, by1;
nuclear@6: 	float rx0, rx1, ry0, ry1;
nuclear@6: 	setup(x, bx0, bx1, rx0, rx1);
nuclear@6: 	setup(y, by0, by1, ry0, ry1);
nuclear@6: 
nuclear@6: 	int i = perm[bx0];
nuclear@6: 	int j = perm[bx1];
nuclear@6: 
nuclear@6: 	int b00 = perm[i + by0];
nuclear@6: 	int b10 = perm[j + by0];
nuclear@6: 	int b01 = perm[i + by1];
nuclear@6: 	int b11 = perm[j + by1];
nuclear@6: 
nuclear@6: 	float sx = s_curve(rx0);
nuclear@6: 	float sy = s_curve(ry0);
nuclear@6: 
nuclear@6: 	vec2_t g00 = grad2[b00];
nuclear@6: 	vec2_t g10 = grad2[b10];
nuclear@6: 	vec2_t g01 = grad2[b01];
nuclear@6: 	vec2_t g11 = grad2[b11];
nuclear@6: 
nuclear@6: 	float u = g00.x * rx0 + g00.y * ry0;
nuclear@6: 	float v = g10.x * rx1 + g10.y * ry0;
nuclear@6: 	float a = lerp(u, v, sx);
nuclear@6: 
nuclear@6: 	u = g01.x * rx0 + g01.y * ry1;
nuclear@6: 	v = g11.x * rx1 + g11.y * ry1;
nuclear@6: 	float b = lerp(u, v, sx);
nuclear@6: 
nuclear@6: 	return lerp(a, b, sy);
nuclear@6: }
nuclear@6: 
nuclear@6: float pnoise2(float x, float y, int periodx, int periody)
nuclear@6: {
nuclear@6: 	if(!tables_valid) {
nuclear@6: 		init_noise();
nuclear@6: 		tables_valid = true;
nuclear@6: 	}
nuclear@6: 
nuclear@6: 	int bx0, bx1, by0, by1;
nuclear@6: 	float rx0, rx1, ry0, ry1;
nuclear@6: 	setup_p(x, bx0, bx1, rx0, rx1, periodx);
nuclear@6: 	setup_p(y, by0, by1, ry0, ry1, periody);
nuclear@6: 
nuclear@6: 	int i = perm[bx0];
nuclear@6: 	int j = perm[bx1];
nuclear@6: 
nuclear@6: 	int b00 = perm[i + by0];
nuclear@6: 	int b10 = perm[j + by0];
nuclear@6: 	int b01 = perm[i + by1];
nuclear@6: 	int b11 = perm[j + by1];
nuclear@6: 
nuclear@6: 	float sx = s_curve(rx0);
nuclear@6: 	float sy = s_curve(ry0);
nuclear@6: 
nuclear@6: 	vec2_t g00 = grad2[b00];
nuclear@6: 	vec2_t g10 = grad2[b10];
nuclear@6: 	vec2_t g01 = grad2[b01];
nuclear@6: 	vec2_t g11 = grad2[b11];
nuclear@6: 
nuclear@6: 	float u = g00.x * rx0 + g00.y * ry0;
nuclear@6: 	float v = g10.x * rx1 + g10.y * ry0;
nuclear@6: 	float a = lerp(u, v, sx);
nuclear@6: 
nuclear@6: 	u = g01.x * rx0 + g01.y * ry1;
nuclear@6: 	v = g11.x * rx1 + g11.y * ry1;
nuclear@6: 	float b = lerp(u, v, sx);
nuclear@6: 
nuclear@6: 	return lerp(a, b, sy);
nuclear@6: }
nuclear@6: 
nuclear@6: float pturbulence2(float x, float y, int periodx, int periody, int octaves)
nuclear@6: {
nuclear@6: 	int i;
nuclear@6: 	float res = 0.0f, freq = 1.0f;
nuclear@6: 	for(i=0; i<octaves; i++) {
nuclear@6: 		res += fabs(pnoise2(x * freq, y * freq, periodx, periody) / freq);
nuclear@6: 		freq *= 2.0f;
nuclear@6: 		periodx *= 2;
nuclear@6: 		periody *= 2;
nuclear@6: 	}
nuclear@6: 	return res;
nuclear@6: }