diff options
author | Sergey Sharybin <sergey.vfx@gmail.com> | 2012-12-10 21:05:33 +0400 |
---|---|---|
committer | Sergey Sharybin <sergey.vfx@gmail.com> | 2012-12-10 21:05:33 +0400 |
commit | decffdc556794541ba0a724c9d4eff052e6a53b5 (patch) | |
tree | fb7c6b8eaac7916b387e929f972da150f8ac6d39 /intern/audaspace/FX/AUD_ButterworthFactory.cpp | |
parent | 14d3870c2fc51670856a56d2ed3aa759d91d4d78 (diff) | |
parent | 245345fba5f44b5c8562c826ca0053fa5bbd2ff3 (diff) |
Merging r52851 through r52858 from trunk into blender-2.65-release
Diffstat (limited to 'intern/audaspace/FX/AUD_ButterworthFactory.cpp')
-rw-r--r-- | intern/audaspace/FX/AUD_ButterworthFactory.cpp | 37 |
1 files changed, 2 insertions, 35 deletions
diff --git a/intern/audaspace/FX/AUD_ButterworthFactory.cpp b/intern/audaspace/FX/AUD_ButterworthFactory.cpp index 97d85c8122f..12c4306c2f7 100644 --- a/intern/audaspace/FX/AUD_ButterworthFactory.cpp +++ b/intern/audaspace/FX/AUD_ButterworthFactory.cpp @@ -29,44 +29,11 @@ #include "AUD_ButterworthFactory.h" #include "AUD_IIRFilterReader.h" - -#include <cmath> - -#ifndef M_PI -#define M_PI 3.14159265358979323846 -#endif - -#define BWPB41 0.76536686473 -#define BWPB42 1.84775906502 +#include "AUD_ButterworthCalculator.h" AUD_ButterworthFactory::AUD_ButterworthFactory(boost::shared_ptr<AUD_IFactory> factory, float frequency) : - AUD_DynamicIIRFilterFactory(factory), - m_frequency(frequency) + AUD_DynamicIIRFilterFactory(factory, boost::shared_ptr<AUD_IDynamicIIRFilterCalculator>(new AUD_ButterworthCalculator(frequency))) { } -void AUD_ButterworthFactory::recalculateCoefficients(AUD_SampleRate rate, - std::vector<float> &b, - std::vector<float> &a) -{ - float omega = 2 * tan(m_frequency * M_PI / rate); - float o2 = omega * omega; - float o4 = o2 * o2; - float x1 = o2 + 2.0f * (float)BWPB41 * omega + 4.0f; - float x2 = o2 + 2.0f * (float)BWPB42 * omega + 4.0f; - float y1 = o2 - 2.0f * (float)BWPB41 * omega + 4.0f; - float y2 = o2 - 2.0f * (float)BWPB42 * omega + 4.0f; - float o228 = 2.0f * o2 - 8.0f; - float norm = x1 * x2; - a.push_back(1); - a.push_back((x1 + x2) * o228 / norm); - a.push_back((x1 * y2 + x2 * y1 + o228 * o228) / norm); - a.push_back((y1 + y2) * o228 / norm); - a.push_back(y1 * y2 / norm); - b.push_back(o4 / norm); - b.push_back(4 * o4 / norm); - b.push_back(6 * o4 / norm); - b.push_back(b[1]); - b.push_back(b[0]); -} |