[ViewVC] Diff of: cvs/deliantra/server/include/noise.h

Comparing deliantra/server/include/noise.h (file contents):
Revision 1.6 by root, Sat Apr 23 04:56:49 2011 UTC vs.
Revision 1.22 by root, Sat Nov 17 23:40:01 2018 UTC

1	/*	1	/*
2	* This file is part of Deliantra, the Roguelike Realtime MMORPG.	2	* This file is part of Deliantra, the Roguelike Realtime MMORPG.
3	*	3	*
		4	* Copyright (©) 2017,2018 Marc Alexander Lehmann / the Deliantra team
4	* Copyright (©) 2010,2011 Marc Alexander Lehmann / Robin Redeker / the Deliantra team	5	* Copyright (©) 2010,2011,2012,2013,2014,2015,2016 Marc Alexander Lehmann / Robin Redeker / the Deliantra team
5	*	6	*
6	* Deliantra is free software: you can redistribute it and/or modify it under	7	* Deliantra is free software: you can redistribute it and/or modify it under
7	* the terms of the Affero GNU General Public License as published by the	8	* the terms of the Affero GNU General Public License as published by the
8	* Free Software Foundation, either version 3 of the License, or (at your	9	* Free Software Foundation, either version 3 of the License, or (at your
9	* option) any later version.	10	* option) any later version.
10	*	11	*
11	* This program is distributed in the hope that it will be useful,	12	* This program is distributed in the hope that it will be useful,
12	* but WITHOUT ANY WARRANTY; without even the implied warranty of	13	* but WITHOUT ANY WARRANTY; without even the implied warranty of
13	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the	14	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14	* GNU General Public License for more details.	15	* GNU General Public License for more details.
15	*	16	*
16	* You should have received a copy of the Affero GNU General Public License	17	* You should have received a copy of the Affero GNU General Public License
17	* and the GNU General Public License along with this program. If not, see	18	* and the GNU General Public License along with this program. If not, see
18	* <http://www.gnu.org/licenses/>.	19	* <http://www.gnu.org/licenses/>.
19	*	20	*
20	* The authors can be reached via e-mail to <support@deliantra.net>	21	* The authors can be reached via e-mail to <support@deliantra.net>
21	*/	22	*/
22		23
23	#ifndef NOISE_H_	24	#ifndef NOISE_H_
24	#define NOISE_H_	25	#define NOISE_H_
25		26
		27	#if 0 // blitz++0.09
		28	#include <string.h> // workaround for tinyvec using memcpy without including string.h
		29	#include <cstdlib> // workaround for tinyvec using labs without including cstdlib
26	#include <blitz/tinyvec.h>	30	#include <blitz/tinyvec.h>
27	#include <blitz/tinyvec-et.h>	31	#include <blitz/tinyvec-et.h>
28
29	#include "global.h"
30
31	typedef blitz::TinyVector<float,2> vec2d;	32	typedef blitz::TinyVector<float,2> vec2d;
32	typedef blitz::TinyVector<float,3> vec3d;	33	typedef blitz::TinyVector<float,3> vec3d;
		34	#else
33		35
34	/////////////////////////////////////////////////////////////////////////////	36	// blitz++ 0.10 - not working
		37	#include <blitz/array.h>
		38	#include <blitz/tinyvec2.h>
		39	#include <blitz/tinyvec2.cc>
35		40
		41	template<typename T, int length>
		42	struct vecnd
		43	: blitz::TinyVector<T, length>
		44	{
		45	static const int numElements = length;
		46	using blitz::TinyVector<T, length>::TinyVector;
		47	};
		48
		49	typedef vecnd<float,2> vec2d;
		50	typedef vecnd<float,3> vec3d;
		51
		52	#endif
		53
		54	#include "global.h"
		55
		56	vec2d
		57	inline floor (vec2d v)
		58	{
		59	return vec2d (fastfloor (v[0]), fastfloor (v[1]));
		60	}
		61
		62	vec3d
		63	inline floor (vec3d v)
		64	{
		65	return vec3d (fastfloor (v[0]), fastfloor (v[1]), fastfloor (v[2]));
		66	}
		67
		68	vec2d
		69	inline pow (vec2d v, vec2d p)
		70	{
		71	return vec2d (pow (v[0], p[0]), pow (v[1], p[1]));
		72	}
		73
		74	/////////////////////////////////////////////////////////////////////////////
		75
		76	// various s-shaped curves, smooth to, first, or second derivative
		77	// used for smooth interpolation from 0..1
		78
		79	// linear
		80	template<typename T>
		81	inline T
		82	sigmoid0 (T x)
		83	{
		84	return x;
		85	}
		86
		87	// 3x²-2x³
		88	template<typename T>
		89	inline T
		90	sigmoid1 (T x)
		91	{
		92	return (3 - 2 * x) * x * x;
		93	}
		94
		95	// 6x⁵ - 15x⁴ + 10x³
		96	template<typename T>
		97	inline T
		98	sigmoid2 (T x)
		99	{
		100	#ifdef MORE_PARALLELITY
		101	float x2 = x * x;
		102	float x4 = x2 * x2;
		103
		104	return (6 * x4 + 10 * x2) * x - 15 * x4;
		105	#endif
		106
		107	// simple horner
		108	return ((6 * x - 15) * x + 10) * x * x * x;
		109	}
		110
		111	// blend between a and b
		112	// c is the control function - if lower than ca
		113	// then return a, if higher than cb, return b
		114	template<typename T, typename U>
		115	inline T
		116	blend (T a, T b, U c, U ca, U cb, U weight (U) = sigmoid1)
		117	{
		118	if (c <= ca) return a;
		119	if (c >= cb) return b;
		120
		121	U w = weight ((c - ca) / (cb - ca));
		122	return (U(1) - w) * a + w * b;
		123	}
		124
		125	// blend between a and b
		126	// c is the control function - if lower than -c_w
		127	// then return a, if higher than +c_w then b.
		128	template<typename T, typename U>
		129	inline T
		130	blend0 (T a, T b, U c, U c_w, U weight (U) = sigmoid1)
		131	{
		132	return blend<T,U> (a, b, c, -c_w, c_w, weight);
		133	}
		134
		135	template<class vec_t, int a, int b>
		136	struct rotate_nn
		137	{
		138	typename vec_t::T_numtype s, c;
		139
		140	void set (typename vec_t::T_numtype angle)
		141	{
		142	s = sin (angle);
		143	c = cos (angle);
		144	}
		145
		146	void operator ()(vec_t &P) const
		147	{
		148	vec_t o = P;
		149
		150	P[a] = o[a] * c - o[b] * s;
		151	P[b] = o[a] * s + o[b] * c;
		152	}
		153	};
		154
		155	template<class vec_t>
		156	struct rotate_xy : rotate_nn<vec_t, 0, 1>
		157	{
		158	};
		159
		160	template<class vec_t>
		161	struct rotate_xz : rotate_nn<vec_t, 0, 2>
		162	{
		163	};
		164
		165	template<class vec_t>
		166	struct rotate_yz : rotate_nn<vec_t, 1, 2>
		167	{
		168	};
		169
		170	/////////////////////////////////////////////////////////////////////////////
		171
36	template<int N, typename T = uint8_t>	172	template< int N, typename T>
37	struct permutation	173	struct permutation
38	{	174	{
39	T pmap[N];	175	T pmap[N];
40		176
41	template<class generator>	177	void seed (seedable_rand_gen &rng);
42	void seed (generator &rng)
43	{
44	for (int i = 0; i < N; ++i)
45	pmap[i] = i;
46		178
47	// fisher-yates to randomly perturb
48	for (int i = N; --i; )
49	::swap (pmap[i], pmap[rng (i + 1)]);
50	}
51
52	T operator ()(T v)	179	ecb_pure T operator ()(T v)
53	{	180	{
54	return pmap[v & (N - 1)];	181	return pmap[v & T(N - 1)];
55	}	182	}
		183	};
		184
		185	/////////////////////////////////////////////////////////////////////////////
		186
		187	template<class vec_t>
		188	struct noise_gen_base
		189	{
		190	permutation<256, uint8_t> rvmap[vec_t::numElements + 1];
		191
		192	typedef typename vec_t::T_numtype value_t;
		193
		194	void seed (seedable_rand_gen &rng);
		195	void seed (seed_t seed);
		196
		197	ecb_pure value_t operator ()(vec_t P, uint32_t seed = 0);
56	};	198	};
57		199
58	template<class vec_t>	200	template<class vec_t>
59	struct noise_gen;	201	struct noise_gen;
60		202
61	// modelled after 2d/3d kensler noise without projection	203	// modelled after 2d/3d kensler noise without projection
62	template<>	204	template<>
63	struct noise_gen<vec2d>	205	struct noise_gen<vec2d>
		206	: noise_gen_base<vec2d>
64	{	207	{
65	permutation<256, uint8_t> rvmap[2];	208	static value_t abs_avg() { return 0.2231; } // avg(abs(noise))
66	vec2d rvec[256]; // random unit vectors
67
68	noise_gen<vec2d> (uint32_t seed);
69	vec2d::T_numtype operator() (vec2d P);
70	};	209	};
71		210
72	template<>	211	template<>
73	struct noise_gen<vec3d>	212	struct noise_gen<vec3d>
		213	: noise_gen_base<vec3d>
74	{	214	{
75	permutation<256, uint8_t> rvmap[3];	215	static vec3d::T_numtype abs_avg() { return 0.415; } // avg(abs(noise))
76		216
77	noise_gen<vec3d> (uint32_t seed);	217	using noise_gen_base<vec3d>::operator ();
78	vec2d::T_numtype operator() (vec3d P);
79		218
80	// noise projected on a surface with normal n	219	// noise projected on a surface with normal n
81	vec2d::T_numtype operator() (vec3d P, vec3d N);	220	ecb_pure vec2d::T_numtype operator() (vec3d P, vec3d N, uint32_t seed = 0);
82	};	221	};
83		222
84	template<class vec_t, int a, int b>	223	typedef noise_gen<vec2d> noise2d;
85	struct rotate_nn	224	typedef noise_gen<vec3d> noise3d;
86	{
87	typename vec_t::T_numtype s, c;
88		225
89	void set (typename vec_t::T_numtype angle)	226	/////////////////////////////////////////////////////////////////////////////
90	{
91	s = sin (angle);
92	c = cos (angle);
93	}
94
95	void operator ()(vec_t &P) const
96	{
97	vec_t o = P;
98
99	P[a] = o[a] * c - o[b] * s;
100	P[b] = o[a] * s + o[b] * c;
101	}
102	};
103
104	template<class vec_t>
105	struct rotate_xy : rotate_nn<vec_t, 0, 1>
106	{
107	};
108
109	template<class vec_t>
110	struct rotate_xz : rotate_nn<vec_t, 0, 2>
111	{
112	};
113
114	template<class vec_t>
115	struct rotate_yz : rotate_nn<vec_t, 1, 2>
116	{
117	};
118		227
119	template<class vec_t>	228	template<class vec_t>
120	struct frac_gen	229	struct frac_gen
121	: noise_gen<vec_t>	230	: noise_gen<vec_t>
122	{	231	{
123	enum { MAX_OCTAVES = 64 };	232	enum { MAX_OCTAVES = 32 };
124		233
125	typedef typename vec_t::T_numtype value_t;	234	typedef typename vec_t::T_numtype value_t;
126		235
		236	int octaves;
127	value_t h, lac, ex[MAX_OCTAVES];	237	value_t h, lac, ex[MAX_OCTAVES];
128	value_t fbm_mul[MAX_OCTAVES];	238	value_t fbm_mul[MAX_OCTAVES];
129	rotate_xy<vec_t> rot[MAX_OCTAVES];	239	rotate_xy<vec_t> rot[MAX_OCTAVES];
130		240
131	value_t noise (vec_t P)	241	frac_gen (int octaves = 3, value_t lacunarity = 2, value_t hurst_expo = .5, uint32_t seed = 0);
		242
		243	ecb_pure value_t noise (vec_t P, uint32_t seed = 0)
132	{	244	{
133	return operator() (P);	245	return this->operator() (P, seed);
134	}	246	}
135		247
136	frac_gen (value_t hurst_expo = .5f, value_t lacunarity = 2.f);	248	ecb_pure value_t fBm (vec_t P);
137		249	ecb_pure value_t turbulence (vec_t P);
138	value_t fBm (vec_t P, int octaves);
139	value_t turbulence (vec_t P, int octaves);
140	value_t multifractal (vec_t P, int octaves, value_t offset = 1);	250	ecb_pure value_t multifractal (vec_t P, value_t offset = 1);
141	value_t heterofractal (vec_t P, int octaves, value_t offset = 1);	251	ecb_pure value_t heterofractal (vec_t P, value_t offset = 1);
142	value_t hybridfractal (vec_t P, int octaves, value_t offset = 1, value_t gain = 1);	252	ecb_pure value_t hybridfractal (vec_t P, value_t offset = 1, value_t gain = 1);
143	value_t terrain (vec_t P, int octaves);
144	value_t terrain2 (vec_t P, int octaves);
145	value_t ridgedmultifractal (vec_t P, int octaves, value_t offset = 1, value_t gain = 8);	253	ecb_pure value_t ridgedmultifractal (vec_t P, value_t offset = 1, value_t gain = 8);
		254	ecb_pure value_t billowfractal (vec_t P, value_t offset = 1, value_t gain = 2);
		255	ecb_pure value_t terrain (vec_t P);
		256	ecb_pure value_t terrain2 (vec_t P);
146	};	257	};
		258
		259	typedef frac_gen<vec2d> frac2d;
		260	typedef frac_gen<vec3d> frac3d;
		261
		262	/////////////////////////////////////////////////////////////////////////////
		263
		264	template<typename T, typename U>
		265	inline T
		266	border_blend (T a, T b, vec2d P, U N, U W)
		267	{
		268	U border = W; // within n places of the border
		269
		270	min_it (border, P [0]);
		271	min_it (border, N - P [0]);
		272	min_it (border, P [1]);
		273	min_it (border, N - P [1]);
		274
		275	return blend (a, b, border, U(0), W);
		276	}
147		277
148	#endif	278	#endif
149		279

Diff Legend

-–
+Removed lines
-+
+Added lines
-<
+Changed lines
->
+Changed lines

Comparing deliantra/server/include/noise.h (file contents): Revision 1.6 by root, Sat Apr 23 04:56:49 2011 UTC vs. Revision 1.22 by root, Sat Nov 17 23:40:01 2018 UTC

Diff Legend

Comparing deliantra/server/include/noise.h (file contents):
Revision 1.6 by root, Sat Apr 23 04:56:49 2011 UTC vs.
Revision 1.22 by root, Sat Nov 17 23:40:01 2018 UTC