llwindgen.h
上传用户:king477883
上传日期:2021-03-01
资源大小:9553k
文件大小:4k
源码类别:

游戏引擎

开发平台:

C++ Builder

  1. /** 
  2.  * @file windgen.h
  3.  * @brief Templated wind noise generation
  4.  *
  5.  * $LicenseInfo:firstyear=2002&license=viewergpl$
  6.  * 
  7.  * Copyright (c) 2002-2010, Linden Research, Inc.
  8.  * 
  9.  * Second Life Viewer Source Code
  10.  * The source code in this file ("Source Code") is provided by Linden Lab
  11.  * to you under the terms of the GNU General Public License, version 2.0
  12.  * ("GPL"), unless you have obtained a separate licensing agreement
  13.  * ("Other License"), formally executed by you and Linden Lab.  Terms of
  14.  * the GPL can be found in doc/GPL-license.txt in this distribution, or
  15.  * online at http://secondlifegrid.net/programs/open_source/licensing/gplv2
  16.  * 
  17.  * There are special exceptions to the terms and conditions of the GPL as
  18.  * it is applied to this Source Code. View the full text of the exception
  19.  * in the file doc/FLOSS-exception.txt in this software distribution, or
  20.  * online at
  21.  * http://secondlifegrid.net/programs/open_source/licensing/flossexception
  22.  * 
  23.  * By copying, modifying or distributing this software, you acknowledge
  24.  * that you have read and understood your obligations described above,
  25.  * and agree to abide by those obligations.
  26.  * 
  27.  * ALL LINDEN LAB SOURCE CODE IS PROVIDED "AS IS." LINDEN LAB MAKES NO
  28.  * WARRANTIES, EXPRESS, IMPLIED OR OTHERWISE, REGARDING ITS ACCURACY,
  29.  * COMPLETENESS OR PERFORMANCE.
  30.  * $/LicenseInfo$
  31.  */
  32. #ifndef WINDGEN_H
  33. #define WINDGEN_H
  34. #include "llcommon.h"
  35. #include "llrand.h"
  36. template <class MIXBUFFERFORMAT_T>
  37. class LLWindGen
  38. {
  39. public:
  40. LLWindGen() :
  41. mTargetGain(0.f),
  42. mTargetFreq(100.f),
  43. mTargetPanGainR(0.5f),
  44. mbuf0(0.0),
  45. mbuf1(0.0),
  46. mbuf2(0.0),
  47. mbuf3(0.0),
  48. mbuf4(0.0),
  49. mbuf5(0.0),
  50. mY0(0.0),
  51. mY1(0.0),
  52. mCurrentGain(0.f),
  53. mCurrentFreq(100.f),
  54. mCurrentPanGainR(0.5f) {};
  55. static const U32 getInputSamplingRate() {return mInputSamplingRate;}
  56. // newbuffer = the buffer passed from the previous DSP unit.
  57. // numsamples = length in samples-per-channel at this mix time.
  58. // stride = number of bytes between start of each sample.
  59. // NOTE: generates L/R interleaved stereo
  60. MIXBUFFERFORMAT_T* windGenerate(MIXBUFFERFORMAT_T *newbuffer, int numsamples, int stride)
  61. {
  62. U8 *cursamplep = (U8*)newbuffer;
  63. double bandwidth = 50.0F;
  64. double a0,b1,b2;
  65. // calculate resonant filter coeffs
  66. b2 = exp(-(F_TWO_PI) * (bandwidth / mInputSamplingRate));
  67. while (numsamples--)
  68. {
  69. mCurrentFreq = (float)((0.999 * mCurrentFreq) + (0.001 * mTargetFreq));
  70. mCurrentGain = (float)((0.999 * mCurrentGain) + (0.001 * mTargetGain));
  71. mCurrentPanGainR = (float)((0.999 * mCurrentPanGainR) + (0.001 * mTargetPanGainR));
  72. b1 = (-4.0 * b2) / (1.0 + b2) * cos(F_TWO_PI * (mCurrentFreq / mInputSamplingRate));
  73. a0 = (1.0 - b2) * sqrt(1.0 - (b1 * b1) / (4.0 * b2));
  74. double nextSample;
  75. // start with white noise
  76. nextSample = ll_frand(2.0f) - 1.0f;
  77. // apply pinking filter
  78. mbuf0 = 0.997f * mbuf0 + 0.0126502f * nextSample; 
  79. mbuf1 = 0.985f * mbuf1 + 0.0139083f * nextSample;
  80. mbuf2 = 0.950f * mbuf2 + 0.0205439f * nextSample;
  81. mbuf3 = 0.850f * mbuf3 + 0.0387225f * nextSample;
  82. mbuf4 = 0.620f * mbuf4 + 0.0465932f * nextSample;
  83. mbuf5 = 0.250f * mbuf5 + 0.1093477f * nextSample;
  84. nextSample = mbuf0 + mbuf1 + mbuf2 + mbuf3 + mbuf4 + mbuf5;
  85. // do a resonant filter on the noise
  86. nextSample = (double)( a0 * nextSample - b1 * mY0 - b2 * mY1 );
  87. mY1 = mY0;
  88. mY0 = nextSample;
  89. nextSample *= mCurrentGain;
  90. MIXBUFFERFORMAT_T sample;
  91. sample = llfloor(((F32)nextSample*32768.f*(1.0f - mCurrentPanGainR))+0.5f);
  92. *(MIXBUFFERFORMAT_T*)cursamplep = llclamp(sample, (MIXBUFFERFORMAT_T)-32768, (MIXBUFFERFORMAT_T)32767);
  93. cursamplep += stride;
  94. sample = llfloor(((F32)nextSample*32768.f*mCurrentPanGainR)+0.5f);
  95. *(MIXBUFFERFORMAT_T*)cursamplep = llclamp(sample, (MIXBUFFERFORMAT_T)-32768, (MIXBUFFERFORMAT_T)32767);
  96. cursamplep += stride;
  97. }
  98. return newbuffer;
  99. }
  100. F32 mTargetGain;
  101. F32 mTargetFreq;
  102. F32 mTargetPanGainR;
  103. private:
  104. static const U32 mInputSamplingRate = 44100;
  105. F64 mbuf0;
  106. F64 mbuf1;
  107. F64 mbuf2;
  108. F64 mbuf3;
  109. F64 mbuf4;
  110. F64 mbuf5;
  111. F64 mY0;
  112. F64 mY1;
  113. F32 mCurrentGain;
  114. F32 mCurrentFreq;
  115. F32 mCurrentPanGainR;
  116. };
  117. #endif