aboutsummaryrefslogtreecommitdiffstats
path: root/alc/filters/splitter.cpp
diff options
context:
space:
mode:
authorChris Robinson <[email protected]>2020-12-04 09:42:13 -0800
committerChris Robinson <[email protected]>2020-12-04 11:15:50 -0800
commit69d55d7e03996484cc899de1e21172a7a4532d6b (patch)
treedf23284a2f4d6d01cc2c9cf8c4fb26f62652ad24 /alc/filters/splitter.cpp
parent84d47f7d4c2d1355a6eb914dd091b39683f83c15 (diff)
Move the filters to core
Diffstat (limited to 'alc/filters/splitter.cpp')
-rw-r--r--alc/filters/splitter.cpp113
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>;