diff options
author | Chris Robinson <[email protected]> | 2020-12-04 09:42:13 -0800 |
---|---|---|
committer | Chris Robinson <[email protected]> | 2020-12-04 11:15:50 -0800 |
commit | 69d55d7e03996484cc899de1e21172a7a4532d6b (patch) | |
tree | df23284a2f4d6d01cc2c9cf8c4fb26f62652ad24 /alc/filters/splitter.cpp | |
parent | 84d47f7d4c2d1355a6eb914dd091b39683f83c15 (diff) |
Move the filters to core
Diffstat (limited to 'alc/filters/splitter.cpp')
-rw-r--r-- | alc/filters/splitter.cpp | 113 |
1 files changed, 0 insertions, 113 deletions
diff --git a/alc/filters/splitter.cpp b/alc/filters/splitter.cpp deleted file mode 100644 index 5cc670b7..00000000 --- a/alc/filters/splitter.cpp +++ /dev/null @@ -1,113 +0,0 @@ - -#include "config.h" - -#include "splitter.h" - -#include <algorithm> -#include <cmath> -#include <limits> - -#include "math_defs.h" -#include "opthelpers.h" - - -template<typename Real> -void BandSplitterR<Real>::init(Real f0norm) -{ - const Real w{f0norm * al::MathDefs<Real>::Tau()}; - const Real cw{std::cos(w)}; - if(cw > std::numeric_limits<float>::epsilon()) - mCoeff = (std::sin(w) - 1.0f) / cw; - else - mCoeff = cw * -0.5f; - - mLpZ1 = 0.0f; - mLpZ2 = 0.0f; - mApZ1 = 0.0f; -} - -template<typename Real> -void BandSplitterR<Real>::process(const al::span<const Real> input, Real *hpout, Real *lpout) -{ - const Real ap_coeff{mCoeff}; - const Real lp_coeff{mCoeff*0.5f + 0.5f}; - Real lp_z1{mLpZ1}; - Real lp_z2{mLpZ2}; - Real ap_z1{mApZ1}; - auto proc_sample = [ap_coeff,lp_coeff,&lp_z1,&lp_z2,&ap_z1,&lpout](const Real in) noexcept -> Real - { - /* Low-pass sample processing. */ - Real d{(in - lp_z1) * lp_coeff}; - Real lp_y{lp_z1 + d}; - lp_z1 = lp_y + d; - - d = (lp_y - lp_z2) * lp_coeff; - lp_y = lp_z2 + d; - lp_z2 = lp_y + d; - - *(lpout++) = lp_y; - - /* All-pass sample processing. */ - Real ap_y{in*ap_coeff + ap_z1}; - ap_z1 = in - ap_y*ap_coeff; - - /* High-pass generated from removing low-passed output. */ - return ap_y - lp_y; - }; - std::transform(input.cbegin(), input.cend(), hpout, proc_sample); - mLpZ1 = lp_z1; - mLpZ2 = lp_z2; - mApZ1 = ap_z1; -} - -template<typename Real> -void BandSplitterR<Real>::processHfScale(const al::span<Real> samples, const Real hfscale) -{ - const Real ap_coeff{mCoeff}; - const Real lp_coeff{mCoeff*0.5f + 0.5f}; - Real lp_z1{mLpZ1}; - Real lp_z2{mLpZ2}; - Real ap_z1{mApZ1}; - auto proc_sample = [hfscale,ap_coeff,lp_coeff,&lp_z1,&lp_z2,&ap_z1](const Real in) noexcept -> Real - { - /* Low-pass sample processing. */ - Real d{(in - lp_z1) * lp_coeff}; - Real lp_y{lp_z1 + d}; - lp_z1 = lp_y + d; - - d = (lp_y - lp_z2) * lp_coeff; - lp_y = lp_z2 + d; - lp_z2 = lp_y + d; - - /* All-pass sample processing. */ - Real ap_y{in*ap_coeff + ap_z1}; - ap_z1 = in - ap_y*ap_coeff; - - /* High-pass generated by removing the low-passed signal, which is then - * scaled and added back to the low-passed signal. - */ - return (ap_y-lp_y)*hfscale + lp_y; - }; - std::transform(samples.begin(), samples.end(), samples.begin(), proc_sample); - mLpZ1 = lp_z1; - mLpZ2 = lp_z2; - mApZ1 = ap_z1; -} - -template<typename Real> -void BandSplitterR<Real>::applyAllpass(const al::span<Real> samples) const -{ - const Real coeff{mCoeff}; - Real z1{0.0f}; - auto proc_sample = [coeff,&z1](const Real in) noexcept -> Real - { - const Real out{in*coeff + z1}; - z1 = in - out*coeff; - return out; - }; - std::transform(samples.begin(), samples.end(), samples.begin(), proc_sample); -} - - -template class BandSplitterR<float>; -template class BandSplitterR<double>; |