开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 游戏 > 其他游戏 > 查看源码
Noise.cpp
资源名称：MMOPRG_server_dev.rar　[点击查看]
上传用户：ghyvgy
上传日期：2009-05-26
资源大小：547k
文件大小：11k
源码类别：

其他游戏

开发平台：

Python

Noise.cpp：源码内容
  1. /*
  2. s_p_oneil@hotmail.com
  3. Copyright (c) 2000, Sean O'Neil
  4. All rights reserved.
  6. Redistribution and use in source and binary forms, with or without
  7. modification, are permitted provided that the following conditions are met:
  9. * Redistributions of source code must retain the above copyright notice,
  10.   this list of conditions and the following disclaimer.
  11. * Redistributions in binary form must reproduce the above copyright notice,
  12.   this list of conditions and the following disclaimer in the documentation
  13.   and/or other materials provided with the distribution.
  14. * Neither the name of this project nor the names of its contributors
  15.   may be used to endorse or promote products derived from this software
  16.   without specific prior written permission.
  18. THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
  19. AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  20. IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  21. ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
  22. LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
  23. CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
  24. SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
  25. INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
  26. CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
  27. ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  28. POSSIBILITY OF SUCH DAMAGE.
  29. */
  31. #include "Master.h"
  32. #include "Noise.h"
  34. void CNoise::Init(int nDimensions, unsigned int nSeed)
  35. {
  36. m_nDimensions = MIN(nDimensions, MAX_DIMENSIONS);
  37. CRandom r(nSeed);
  39. int i, j, k;
  40. for(i=0; i<256; i++)
  41. {
  42. m_nMap[i] = i;
  43. for(j=0; j<m_nDimensions; j++)
  44. m_nBuffer[i][j] = (float)r.RandomD(-0.5, 0.5);
  45. Normalize(m_nBuffer[i], m_nDimensions);
  46. }
  48. while(--i)
  49. {
  50. j = r.RandomI(0, 255);
  51. SWAP(m_nMap[i], m_nMap[j], k);
  52. }
  53. //_fpreset(); // Bug in CRandom! Causes messed up floating point operations!
  54. }
  56. float CNoise::Noise(float *f)
  57. {
  58. int n[MAX_DIMENSIONS]; // Indexes to pass to lattice function
  59. float r[MAX_DIMENSIONS]; // Remainders to pass to lattice function
  60. float w[MAX_DIMENSIONS]; // Cubic values to pass to interpolation function
  62. for(int i=0; i<m_nDimensions; i++)
  63. {
  64. n[i] = Floor(f[i]);
  65. r[i] = f[i] - n[i];
  66. w[i] = Cubic(r[i]);
  67. }
  69. float fValue;
  70. switch(m_nDimensions)
  71. {
  72. case 1:
  73. fValue = Lerp(Lattice(n[0], r[0]),
  74.   Lattice(n[0]+1, r[0]-1),
  75.   w[0]);
  76. break;
  77. case 2:
  78. fValue = Lerp(Lerp(Lattice(n[0], r[0], n[1], r[1]),
  79.    Lattice(n[0]+1, r[0]-1, n[1], r[1]),
  80.    w[0]),
  81.   Lerp(Lattice(n[0], r[0], n[1]+1, r[1]-1),
  82.    Lattice(n[0]+1, r[0]-1, n[1]+1, r[1]-1),
  83.    w[0]),
  84.   w[1]);
  85. break;
  86. case 3:
  87. fValue = Lerp(Lerp(Lerp(Lattice(n[0], r[0], n[1], r[1], n[2], r[2]),
  88. Lattice(n[0]+1, r[0]-1, n[1], r[1], n[2], r[2]),
  89. w[0]),
  90.    Lerp(Lattice(n[0], r[0], n[1]+1, r[1]-1, n[2], r[2]),
  91. Lattice(n[0]+1, r[0]-1, n[1]+1, r[1]-1, n[2], r[2]),
  92. w[0]),
  93.    w[1]),
  94.   Lerp(Lerp(Lattice(n[0], r[0], n[1], r[1], n[2]+1, r[2]-1),
  95. Lattice(n[0]+1, r[0]-1, n[1], r[1], n[2]+1, r[2]-1),
  96. w[0]),
  97.    Lerp(Lattice(n[0], r[0], n[1]+1, r[1]-1, n[2]+1, r[2]-1),
  98. Lattice(n[0]+1, r[0]-1, n[1]+1, r[1]-1, n[2]+1, r[2]-1),
  99. w[0]),
  100.    w[1]),
  101.   w[2]);
  102. break;
  103. case 4:
  104. fValue = Lerp(Lerp(Lerp(Lerp(Lattice(n[0], r[0], n[1], r[1], n[2], r[2], n[3], r[3]),
  105.  Lattice(n[0]+1, r[0]-1, n[1], r[1], n[2], r[2], n[3], r[3]),
  106.  w[0]),
  107. Lerp(Lattice(n[0], r[0], n[1]+1, r[1]-1, n[2], r[2], n[3], r[3]),
  108.  Lattice(n[0]+1, r[0]-1, n[1]+1, r[1]-1, n[2], r[2], n[3], r[3]),
  109.  w[0]),
  110. w[1]),
  111. Lerp(Lerp(Lattice(n[0], r[0], n[1], r[1], n[2]+1, r[2]-1, n[3], r[3]),
  112.  Lattice(n[0]+1, r[0]-1, n[1], r[1], n[2]+1, r[2]-1, n[3], r[3]),
  113.  w[0]),
  114. Lerp(Lattice(n[0], r[0], n[1]+1, r[1]-1, n[2]+1, r[2]-1),
  115.  Lattice(n[0]+1, r[0]-1, n[1]+1, r[1]-1, n[2]+1, r[2]-1, n[3], r[3]),
  116.  w[0]),
  117. w[1]),
  118.    w[2]),
  119.   Lerp(Lerp(Lerp(Lattice(n[0], r[0], n[1], r[1], n[2], r[2], n[3]+1, r[3]-1),
  120.  Lattice(n[0]+1, r[0]-1, n[1], r[1], n[2], r[2], n[3]+1, r[3]-1),
  121.  w[0]),
  122. Lerp(Lattice(n[0], r[0], n[1]+1, r[1]-1, n[2], r[2], n[3]+1, r[3]-1),
  123.  Lattice(n[0]+1, r[0]-1, n[1]+1, r[1]-1, n[2], r[2], n[3]+1, r[3]-1),
  124.  w[0]),
  125. w[1]),
  126. Lerp(Lerp(Lattice(n[0], r[0], n[1], r[1], n[2]+1, r[2]-1, n[3]+1, r[3]-1),
  127.  Lattice(n[0]+1, r[0]-1, n[1], r[1], n[2]+1, r[2]-1, n[3]+1, r[3]-1),
  128.  w[0]),
  129. Lerp(Lattice(n[0], r[0], n[1]+1, r[1]-1, n[2]+1, r[2]-1),
  130.  Lattice(n[0]+1, r[0]-1, n[1]+1, r[1]-1, n[2]+1, r[2]-1, n[3]+1, r[3]-1),
  131.  w[0]),
  132. w[1]),
  133.    w[2]),
  134.   w[3]);
  135. break;
  136. }
  137. return CLAMP(-0.99999f, 0.99999f, fValue);
  138. }
  140. void CSeededNoise::Init(unsigned int nSeed)
  141. {
  142. /*
  143. for(int y=0; y<64; y++)
  144. {
  145. float fDy = (float)y/m_nHeight - 0.5f;
  146. for(int x=0; x<64; x++)
  147. {
  148. float fDx = (float)x/m_nWidth - 0.5f;
  149. float fDist = sqrtf(fDx*fDx + fDy*fDy);
  150. float fIntensity = expf(-Max(fDist-fSizeDisc,0)*fExpose);
  151. m_nBuffer[x][y] = 0;
  152. }
  153. }
  154. */
  155. }
  157. float CSeededNoise::Noise(float *f)
  158. {
  159. int n[2]; // Indexes to pass to lattice function
  160. float r[2]; // Remainders to pass to lattice function
  161. float w[2]; // Cubic values to pass to interpolation function
  163. for(int i=0; i<2; i++)
  164. {
  165. n[i] = Floor(f[i]);
  166. r[i] = f[i] - n[i];
  167. w[i] = Cubic(r[i]);
  168. }
  169. float fValue = Lerp(Lerp(m_nBuffer[n[0]][n[1]],
  170.  m_nBuffer[n[0]+1][n[1]],
  171.  w[0]),
  172. Lerp(m_nBuffer[n[0]][n[1]+1],
  173.  m_nBuffer[n[0]+1][n[1]+1],
  174.  w[0]),
  175. w[1]);
  176. return CLAMP(-0.99999f, 0.99999f, fValue);
  177. }
  179. float CFractal::fBm(float *f, float fOctaves)
  180. {
  181. // Initialize locals
  182. float fValue = 0;
  183. float fTemp[MAX_DIMENSIONS];
  184. for(int i=0; i<m_nDimensions; i++)
  185. fTemp[i] = f[i];
  187. // Inner loop of spectral construction, where the fractal is built
  188. for(i=0; i<fOctaves; i++)
  189. {
  190. fValue += Noise(fTemp) * m_fExponent[i];
  191. for(int j=0; j<m_nDimensions; j++)
  192. fTemp[j] *= m_fLacunarity;
  193. }
  195. // Take care of remainder in fOctaves
  196. fOctaves -= (int)fOctaves;
  197. if(fOctaves > DELTA)
  198. fValue += fOctaves * Noise(fTemp) * m_fExponent[i];
  199. return CLAMP(-0.99999f, 0.99999f, fValue);
  200. }
  202. float CFractal::fBmTest(float *f, int nStart, int nEnd, float fInitial)
  203. {
  204. float fTemp[MAX_DIMENSIONS];
  205. float fValue = 0, fExp = 2;
  206. int i;
  208. // Initialize locals
  209. for(i=0; i<nStart; i++)
  210. fExp *= m_fLacunarity;
  211. for(i=0; i<m_nDimensions; i++)
  212. fTemp[i] = f[i] * fExp;
  214. // Inner loop of spectral construction, where the fractal is built
  215. for(i=nStart; i<nEnd; i++)
  216. {
  217. fValue += Noise(fTemp) * m_fExponent[i];
  218. for(int j=0; j<m_nDimensions; j++)
  219. fTemp[j] *= m_fLacunarity;
  220. }
  222. if(nStart)
  223. fValue += fInitial;
  224. if(fValue <= 0.0f)
  225. fValue = (float)-pow(-fValue, 0.7f);
  226. else
  227. fValue = (float)pow(fValue, 1 + Noise(fTemp) * fValue);
  228. return fValue * 1.33333f;
  229. }
  231. float CFractal::fBmTest(float *f, float fOctaves)
  232. {
  233. float fTemp[MAX_DIMENSIONS];
  234. float fValue = 0;
  235. int i;
  237. // Initialize locals
  238. for(i=0; i<m_nDimensions; i++)
  239. fTemp[i] = f[i] * 2;
  241. // Inner loop of spectral construction, where the fractal is built
  242. for(i=0; i<fOctaves; i++)
  243. {
  244. fValue += Noise(fTemp) * m_fExponent[i];
  245. for(int j=0; j<m_nDimensions; j++)
  246. fTemp[j] *= m_fLacunarity;
  247. }
  249. // Take care of remainder in fOctaves
  250. fOctaves -= (int)fOctaves;
  251. if(fOctaves > DELTA)
  252. fValue += fOctaves * Noise(fTemp) * m_fExponent[i];
  254. if(fValue <= 0.0f)
  255. fValue = (float)-pow(-fValue, 0.7f);
  256. else
  257. fValue = (float)pow(fValue, 1 + Noise(fTemp) * fValue);
  258. return fValue * 1.33333f;
  259. }
  261. float CFractal::Turbulence(float *f, float fOctaves)
  262. {
  263. // Initialize locals
  264. float fValue = 0;
  265. float fTemp[MAX_DIMENSIONS];
  266. for(int i=0; i<m_nDimensions; i++)
  267. fTemp[i] = f[i];
  269. // Inner loop of spectral construction, where the fractal is built
  270. for(i=0; i<fOctaves; i++)
  271. {
  272. fValue += Abs(Noise(fTemp)) * m_fExponent[i];
  273. for(int j=0; j<m_nDimensions; j++)
  274. fTemp[j] *= m_fLacunarity;
  275. }
  277. // Take care of remainder in fOctaves
  278. fOctaves -= (int)fOctaves;
  279. if(fOctaves > DELTA)
  280. fValue += fOctaves * Abs(Noise(fTemp) * m_fExponent[i]);
  281. return CLAMP(-0.99999f, 0.99999f, fValue);
  282. }
  284. float CFractal::Multifractal(float *f, float fOctaves, float fOffset)
  285. {
  286. // Initialize locals
  287. float fValue = 1;
  288. float fTemp[MAX_DIMENSIONS];
  289. for(int i=0; i<m_nDimensions; i++)
  290. fTemp[i] = f[i];
  292. // Inner loop of spectral construction, where the fractal is built
  293. for(i=0; i<fOctaves; i++)
  294. {
  295. fValue *= Noise(fTemp) * m_fExponent[i] + fOffset;
  296. for(int j=0; j<m_nDimensions; j++)
  297. fTemp[j] *= m_fLacunarity;
  298. }
  300. // Take care of remainder in fOctaves (shouldn't that be a multiply?)
  301. fOctaves -= (int)fOctaves;
  302. if(fOctaves > DELTA)
  303. fValue *= fOctaves * (Noise(fTemp) * m_fExponent[i] + fOffset);
  304. return CLAMP(-0.99999f, 0.99999f, fValue);
  305. }
  307. float CFractal::Heterofractal(float *f, float fOctaves, float fOffset)
  308. {
  309. // Initialize locals
  310. float fValue = Noise(f) + fOffset;
  311. float fTemp[MAX_DIMENSIONS];
  312. for(int i=0; i<m_nDimensions; i++)
  313. fTemp[i] = f[i] * m_fLacunarity;
  315. // Inner loop of spectral construction, where the fractal is built
  316. for(i=1; i<fOctaves; i++)
  317. {
  318. fValue += (Noise(fTemp) + fOffset) * m_fExponent[i] * fValue;
  319. for(int j=0; j<m_nDimensions; j++)
  320. fTemp[j] *= m_fLacunarity;
  321. }
  323. // Take care of remainder in fOctaves
  324. fOctaves -= (int)fOctaves;
  325. if(fOctaves > DELTA)
  326. fValue += fOctaves * (Noise(fTemp) + fOffset) * m_fExponent[i] * fValue;
  327. return CLAMP(-0.99999f, 0.99999f, fValue);
  328. }
  330. float CFractal::HybridMultifractal(float *f, float fOctaves, float fOffset, float fGain)
  331. {
  332. // Initialize locals
  333. float fValue = (Noise(f) + fOffset) * m_fExponent[0];
  334. float fWeight = fValue;
  335. float fTemp[MAX_DIMENSIONS];
  336. for(int i=0; i<m_nDimensions; i++)
  337. fTemp[i] = f[i] * m_fLacunarity;
  339. // Inner loop of spectral construction, where the fractal is built
  340. for(i=1; i<fOctaves; i++)
  341. {
  342. if(fWeight > 1)
  343. fWeight = 1;
  344. float fSignal = (Noise(fTemp) + fOffset) * m_fExponent[i];
  345. fValue += fWeight * fSignal;
  346. fWeight *= fGain * fSignal;
  347. for(int j=0; j<m_nDimensions; j++)
  348. fTemp[j] *= m_fLacunarity;
  349. }
  351. // Take care of remainder in fOctaves
  352. fOctaves -= (int)fOctaves;
  353. if(fOctaves > DELTA)
  354. {
  355. if(fWeight > 1)
  356. fWeight = 1;
  357. float fSignal = (Noise(fTemp) + fOffset) * m_fExponent[i];
  358. fValue += fOctaves * fWeight * fSignal;
  359. }
  360. return CLAMP(-0.99999f, 0.99999f, fValue);
  361. }
  363. float CFractal::RidgedMultifractal(float *f, float fOctaves, float fOffset, float fGain)
  364. {
  365. // Initialize locals
  366. float fSignal = fOffset - Abs(Noise(f));
  367. fSignal *= fSignal;
  368. float fValue = fSignal;
  369. float fTemp[MAX_DIMENSIONS];
  370. for(int i=0; i<m_nDimensions; i++)
  371. fTemp[i] = f[i];
  373. // Inner loop of spectral construction, where the fractal is built
  374. for(i=1; i<fOctaves; i++)
  375. {
  376. for(int j=0; j<m_nDimensions; j++)
  377. fTemp[j] *= m_fLacunarity;
  378. float fWeight = Clamp(0, 1, fSignal * fGain);
  379. fSignal = fOffset - Abs(Noise(fTemp));
  380. fSignal *= fSignal;
  381. fSignal *= fWeight;
  382. fValue += fSignal * m_fExponent[i];
  383. }
  384. return CLAMP(-0.99999f, 0.99999f, fValue);
  385. }