server/include/noise.h

/*
 * This file is part of Deliantra, the Roguelike Realtime MMORPG.
 *
 * Copyright (©) 2010,2011,2012 Marc Alexander Lehmann / Robin Redeker / the Deliantra team
 *
 * Deliantra is free software: you can redistribute it and/or modify it under
 * the terms of the Affero GNU General Public License as published by the
 * Free Software Foundation, either version 3 of the License, or (at your
 * option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the Affero GNU General Public License
 * and the GNU General Public License along with this program. If not, see
 * <http://www.gnu.org/licenses/>.
 *
 * The authors can be reached via e-mail to <support@deliantra.net>
 */

#ifndef NOISE_H_
#define NOISE_H_

#if 0 // blitz++0.09
#include <string.h> // workaround for tinyvec using memcpy without including string.h
#include <cstdlib> // workaround for tinyvec using labs without including cstdlib
#include <blitz/tinyvec.h>
#include <blitz/tinyvec-et.h>
typedef blitz::TinyVector<float,2> vec2d;
typedef blitz::TinyVector<float,3> vec3d;
#else

// blitz++ 0.10 - not working
#include <blitz/array.h>
#include <blitz/tinyvec2.h>
#include <blitz/tinyvec2.cc>

template<typename T, int length>
struct vecnd
: blitz::TinyVector<T, length>
{
  static const int numElements = length;
  using blitz::TinyVector<T, length>::TinyVector;
};

typedef vecnd<float,2> vec2d;
typedef vecnd<float,3> vec3d;

#endif

#include "global.h"

vec2d
inline floor (vec2d v)
{
  return vec2d (fastfloor (v[0]), fastfloor (v[1]));
}

vec3d
inline floor (vec3d v)
{
  return vec3d (fastfloor (v[0]), fastfloor (v[1]), fastfloor (v[2]));
}

vec2d
inline pow (vec2d v, vec2d p)
{
  return vec2d (pow (v[0], p[0]), pow (v[1], p[1]));
}

/////////////////////////////////////////////////////////////////////////////

// various s-shaped curves, smooth to, first, or second derivative
// used for smooth interpolation from 0..1

// linear
template<typename T>
inline T
sigmoid0 (T x)
{
  return x;
}

// 3x²-2x³
template<typename T>
inline T
sigmoid1 (T x)
{
  return (3 - 2 * x) * x * x;
}

// 6x⁵ - 15x⁴ + 10x³
template<typename T>
inline T
sigmoid2 (T x)
{
#ifdef MORE_PARALLELITY
  float x2 = x  * x;
  float x4 = x2 * x2;

  return (6 * x4 + 10 * x2) * x - 15 * x4;
#endif

  // simple horner
  return ((6 * x - 15) * x + 10) * x * x * x;
}

// blend between a and b
// c is the control function - if lower than ca
// then return a, if higher than cb, return b
template<typename T, typename U>
inline T
blend (T a, T b, U c, U ca, U cb, U weight (U) = sigmoid1)
{
  if (c <= ca) return a;
  if (c >= cb) return b;

  U w = weight ((c - ca) / (cb - ca));
  return (U(1) - w) * a + w * b;
}

// blend between a and b
// c is the control function - if lower than -c_w
// then return a, if higher than +c_w then b.
template<typename T, typename U>
inline T
blend0 (T a, T b, U c, U c_w, U weight (U) = sigmoid1)
{
  return blend<T,U> (a, b, c, -c_w, c_w, weight);
}

template<class vec_t, int a, int b>
struct rotate_nn
{
  typename vec_t::T_numtype s, c;

  void set (typename vec_t::T_numtype angle)
  {
    s = sin (angle);
    c = cos (angle);
  }

  void operator ()(vec_t &P) const
  {
    vec_t o = P;

    P[a] = o[a] * c - o[b] * s;
    P[b] = o[a] * s + o[b] * c;
  }
};

template<class vec_t>
struct rotate_xy : rotate_nn<vec_t, 0, 1>
{
};

template<class vec_t>
struct rotate_xz : rotate_nn<vec_t, 0, 2>
{
};

template<class vec_t>
struct rotate_yz : rotate_nn<vec_t, 1, 2>
{
};

/////////////////////////////////////////////////////////////////////////////

template< int N, typename T>
struct permutation
{
  T pmap[N];

  void seed (seedable_rand_gen &rng);

  T operator ()(T v) func_pure
  {
    return pmap[v & T(N - 1)];
  }
};

/////////////////////////////////////////////////////////////////////////////

template<class vec_t>
struct noise_gen_base
{
  permutation<256, uint8_t> rvmap[vec_t::numElements + 1];

  typedef typename vec_t::T_numtype value_t;

  void seed (seedable_rand_gen &rng);
  void seed (seed_t seed);

  value_t operator ()(vec_t P, uint32_t seed = 0) func_pure;
};

template<class vec_t>
struct noise_gen;

// modelled after 2d/3d kensler noise without projection
template<>
struct noise_gen<vec2d>
: noise_gen_base<vec2d>
{
  static value_t abs_avg() { return 0.2231; } // avg(abs(noise))
};

template<>
struct noise_gen<vec3d>
: noise_gen_base<vec3d>
{
  static vec3d::T_numtype abs_avg() { return 0.415; } // avg(abs(noise))

  using noise_gen_base<vec3d>::operator ();

  // noise projected on a surface with normal n
  vec2d::T_numtype operator() (vec3d P, vec3d N, uint32_t seed = 0) func_pure;
};

typedef noise_gen<vec2d> noise2d;
typedef noise_gen<vec3d> noise3d;

/////////////////////////////////////////////////////////////////////////////

template<class vec_t>
struct frac_gen
: noise_gen<vec_t>
{
  enum { MAX_OCTAVES = 32 };

  typedef typename vec_t::T_numtype value_t;

  int octaves;
  value_t h, lac, ex[MAX_OCTAVES];
  value_t fbm_mul[MAX_OCTAVES];
  rotate_xy<vec_t> rot[MAX_OCTAVES];

  frac_gen (int octaves = 3, value_t lacunarity = 2, value_t hurst_expo = .5, uint32_t seed = 0);

  value_t noise (vec_t P, uint32_t seed = 0) func_pure
  {
    return this->operator() (P, seed);
  }

  value_t fBm (vec_t P) func_pure;
  value_t turbulence (vec_t P) func_pure;
  value_t multifractal (vec_t P, value_t offset = 1) func_pure;
  value_t heterofractal (vec_t P, value_t offset = 1) func_pure;
  value_t hybridfractal (vec_t P, value_t offset = 1, value_t gain = 1) func_pure;
  value_t ridgedmultifractal (vec_t P, value_t offset = 1, value_t gain = 8) func_pure;
  value_t billowfractal (vec_t P, value_t offset = 1, value_t gain = 2) func_pure;
  value_t terrain (vec_t P) func_pure;
  value_t terrain2 (vec_t P) func_pure;
};

typedef frac_gen<vec2d> frac2d;
typedef frac_gen<vec3d> frac3d;

/////////////////////////////////////////////////////////////////////////////

template<typename T, typename U>
inline T
border_blend (T a, T b, vec2d P, U N, U W)
{
  U border = W; // within n places of the border

  min_it (border,     P [0]);
  min_it (border, N - P [0]);
  min_it (border,     P [1]);
  min_it (border, N - P [1]);

  return blend (a, b, border, U(0), W);
}

#endif

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