mirror of
https://github.com/aicodix/dsp.git
synced 2026-04-27 14:30:36 +00:00
50 lines
772 B
C++
50 lines
772 B
C++
/*
|
|
Some little helpers
|
|
|
|
Copyright 2018 Ahmet Inan <inan@aicodix.de>
|
|
*/
|
|
|
|
#pragma once
|
|
|
|
#include "const.hh"
|
|
|
|
namespace DSP {
|
|
|
|
template <typename TYPE>
|
|
int signum(TYPE v)
|
|
{
|
|
return (v > TYPE(0)) - (v < TYPE(0));
|
|
}
|
|
|
|
template <typename X, typename AB>
|
|
AB lerp(X x, AB a, AB b)
|
|
{
|
|
return (X(1) - x) * a + x * b;
|
|
}
|
|
|
|
template <typename TYPE>
|
|
TYPE normal_pdf(TYPE x, TYPE m, TYPE s)
|
|
{
|
|
return exp(-pow((x - m) / s, TYPE(2)) / TYPE(2)) / (Const<TYPE>::SqrtTwoPi() * s);
|
|
}
|
|
|
|
template <typename TYPE>
|
|
TYPE sinc(TYPE x)
|
|
{
|
|
return TYPE(0) == x ? TYPE(1) : sin(Const<TYPE>::Pi() * x) / (Const<TYPE>::Pi() * x);
|
|
}
|
|
|
|
template <typename TYPE>
|
|
TYPE delta(TYPE x)
|
|
{
|
|
return TYPE(0) == x ? TYPE(1) : TYPE(0);
|
|
}
|
|
|
|
template <typename TYPE>
|
|
TYPE decibel(TYPE v)
|
|
{
|
|
return TYPE(10) * log10(v);
|
|
}
|
|
|
|
}
|
|
|