Loading...
Searching...
No Matches
juce_FastMathApproximations.h
/*
==============================================================================
This file is part of the JUCE library.
Copyright (c) 2022 - Raw Material Software Limited
JUCE is an open source library subject to commercial or open-source
licensing.
By using JUCE, you agree to the terms of both the JUCE 7 End-User License
Agreement and JUCE Privacy Policy.
End User License Agreement: www.juce.com/juce-7-licence
Privacy Policy: www.juce.com/juce-privacy-policy
Or: You may also use this code under the terms of the GPL v3 (see
www.gnu.org/licenses).
JUCE IS PROVIDED "AS IS" WITHOUT ANY WARRANTY, AND ALL WARRANTIES, WHETHER
EXPRESSED OR IMPLIED, INCLUDING MERCHANTABILITY AND FITNESS FOR PURPOSE, ARE
DISCLAIMED.
==============================================================================
*/
namespace juce::dsp
{
struct FastMathApproximations
{
template <typename FloatType>
static FloatType cosh (FloatType x) noexcept
{
auto x2 = x * x;
auto numerator = -(39251520 + x2 * (18471600 + x2 * (1075032 + 14615 * x2)));
auto denominator = -39251520 + x2 * (1154160 + x2 * (-16632 + 127 * x2));
return numerator / denominator;
}
template <typename FloatType>
static void cosh (FloatType* values, size_t numValues) noexcept
{
for (size_t i = 0; i < numValues; ++i)
values[i] = FastMathApproximations::cosh (values[i]);
}
template <typename FloatType>
static FloatType sinh (FloatType x) noexcept
{
auto x2 = x * x;
auto numerator = -x * (11511339840 + x2 * (1640635920 + x2 * (52785432 + x2 * 479249)));
auto denominator = -11511339840 + x2 * (277920720 + x2 * (-3177720 + x2 * 18361));
return numerator / denominator;
}
template <typename FloatType>
static void sinh (FloatType* values, size_t numValues) noexcept
{
for (size_t i = 0; i < numValues; ++i)
values[i] = FastMathApproximations::sinh (values[i]);
}
template <typename FloatType>
static FloatType tanh (FloatType x) noexcept
{
auto x2 = x * x;
auto numerator = x * (135135 + x2 * (17325 + x2 * (378 + x2)));
auto denominator = 135135 + x2 * (62370 + x2 * (3150 + 28 * x2));
return numerator / denominator;
}
template <typename FloatType>
static void tanh (FloatType* values, size_t numValues) noexcept
{
for (size_t i = 0; i < numValues; ++i)
values[i] = FastMathApproximations::tanh (values[i]);
}
//==============================================================================
template <typename FloatType>
static FloatType cos (FloatType x) noexcept
{
auto x2 = x * x;
auto numerator = -(-39251520 + x2 * (18471600 + x2 * (-1075032 + 14615 * x2)));
auto denominator = 39251520 + x2 * (1154160 + x2 * (16632 + x2 * 127));
return numerator / denominator;
}
template <typename FloatType>
static void cos (FloatType* values, size_t numValues) noexcept
{
for (size_t i = 0; i < numValues; ++i)
values[i] = FastMathApproximations::cos (values[i]);
}
template <typename FloatType>
static FloatType sin (FloatType x) noexcept
{
auto x2 = x * x;
auto numerator = -x * (-11511339840 + x2 * (1640635920 + x2 * (-52785432 + x2 * 479249)));
auto denominator = 11511339840 + x2 * (277920720 + x2 * (3177720 + x2 * 18361));
return numerator / denominator;
}
template <typename FloatType>
static void sin (FloatType* values, size_t numValues) noexcept
{
for (size_t i = 0; i < numValues; ++i)
values[i] = FastMathApproximations::sin (values[i]);
}
template <typename FloatType>
static FloatType tan (FloatType x) noexcept
{
auto x2 = x * x;
auto numerator = x * (-135135 + x2 * (17325 + x2 * (-378 + x2)));
auto denominator = -135135 + x2 * (62370 + x2 * (-3150 + 28 * x2));
return numerator / denominator;
}
template <typename FloatType>
static void tan (FloatType* values, size_t numValues) noexcept
{
for (size_t i = 0; i < numValues; ++i)
values[i] = FastMathApproximations::tan (values[i]);
}
//==============================================================================
template <typename FloatType>
static FloatType exp (FloatType x) noexcept
{
auto numerator = 1680 + x * (840 + x * (180 + x * (20 + x)));
auto denominator = 1680 + x *(-840 + x * (180 + x * (-20 + x)));
return numerator / denominator;
}
template <typename FloatType>
static void exp (FloatType* values, size_t numValues) noexcept
{
for (size_t i = 0; i < numValues; ++i)
values[i] = FastMathApproximations::exp (values[i]);
}
template <typename FloatType>
static FloatType logNPlusOne (FloatType x) noexcept
{
auto numerator = x * (7560 + x * (15120 + x * (9870 + x * (2310 + x * 137))));
auto denominator = 7560 + x * (18900 + x * (16800 + x * (6300 + x * (900 + 30 * x))));
return numerator / denominator;
}
template <typename FloatType>
static void logNPlusOne (FloatType* values, size_t numValues) noexcept
{
for (size_t i = 0; i < numValues; ++i)
values[i] = FastMathApproximations::logNPlusOne (values[i]);
}
};
} // namespace juce::dsp
juce::dsp::FastMathApproximations
Definition juce_FastMathApproximations.h:35
juce::dsp::FastMathApproximations::cos
static void cos(FloatType *values, size_t numValues) noexcept
Definition juce_FastMathApproximations.h:143
juce::dsp::FastMathApproximations::tanh
static FloatType tanh(FloatType x) noexcept
Definition juce_FastMathApproximations.h:99
juce::dsp::FastMathApproximations::tanh
static void tanh(FloatType *values, size_t numValues) noexcept
Definition juce_FastMathApproximations.h:114
juce::dsp::FastMathApproximations::sinh
static FloatType sinh(FloatType x) noexcept
Definition juce_FastMathApproximations.h:71
juce::dsp::FastMathApproximations::sin
static FloatType sin(FloatType x) noexcept
Definition juce_FastMathApproximations.h:156
juce::dsp::FastMathApproximations::logNPlusOne
static FloatType logNPlusOne(FloatType x) noexcept
Definition juce_FastMathApproximations.h:240
juce::dsp::FastMathApproximations::sin
static void sin(FloatType *values, size_t numValues) noexcept
Definition juce_FastMathApproximations.h:171
juce::dsp::FastMathApproximations::exp
static FloatType exp(FloatType x) noexcept
Definition juce_FastMathApproximations.h:213
juce::dsp::FastMathApproximations::logNPlusOne
static void logNPlusOne(FloatType *values, size_t numValues) noexcept
Definition juce_FastMathApproximations.h:254
juce::dsp::FastMathApproximations::cos
static FloatType cos(FloatType x) noexcept
Definition juce_FastMathApproximations.h:128
juce::dsp::FastMathApproximations::tan
static FloatType tan(FloatType x) noexcept
Definition juce_FastMathApproximations.h:184
juce::dsp::FastMathApproximations::cosh
static void cosh(FloatType *values, size_t numValues) noexcept
Definition juce_FastMathApproximations.h:58
juce::dsp::FastMathApproximations::sinh
static void sinh(FloatType *values, size_t numValues) noexcept
Definition juce_FastMathApproximations.h:86
juce::dsp::FastMathApproximations::exp
static void exp(FloatType *values, size_t numValues) noexcept
Definition juce_FastMathApproximations.h:227
juce::dsp::FastMathApproximations::cosh
static FloatType cosh(FloatType x) noexcept
Definition juce_FastMathApproximations.h:43
juce::dsp::FastMathApproximations::tan
static void tan(FloatType *values, size_t numValues) noexcept
Definition juce_FastMathApproximations.h:199