xref: /aosp_15_r20/external/oboe/src/flowgraph/Limiter.cpp (revision 05767d913155b055644481607e6fa1e35e2fe72c)
1*05767d91SRobert Wu /*
2*05767d91SRobert Wu  * Copyright 2022 The Android Open Source Project
3*05767d91SRobert Wu  *
4*05767d91SRobert Wu  * Licensed under the Apache License, Version 2.0 (the "License");
5*05767d91SRobert Wu  * you may not use this file except in compliance with the License.
6*05767d91SRobert Wu  * You may obtain a copy of the License at
7*05767d91SRobert Wu  *
8*05767d91SRobert Wu  *      http://www.apache.org/licenses/LICENSE-2.0
9*05767d91SRobert Wu  *
10*05767d91SRobert Wu  * Unless required by applicable law or agreed to in writing, software
11*05767d91SRobert Wu  * distributed under the License is distributed on an "AS IS" BASIS,
12*05767d91SRobert Wu  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13*05767d91SRobert Wu  * See the License for the specific language governing permissions and
14*05767d91SRobert Wu  * limitations under the License.
15*05767d91SRobert Wu  */
16*05767d91SRobert Wu 
17*05767d91SRobert Wu #include <algorithm>
18*05767d91SRobert Wu #include <math.h>
19*05767d91SRobert Wu #include <unistd.h>
20*05767d91SRobert Wu #include "FlowGraphNode.h"
21*05767d91SRobert Wu #include "Limiter.h"
22*05767d91SRobert Wu 
23*05767d91SRobert Wu using namespace FLOWGRAPH_OUTER_NAMESPACE::flowgraph;
24*05767d91SRobert Wu 
Limiter(int32_t channelCount)25*05767d91SRobert Wu Limiter::Limiter(int32_t channelCount)
26*05767d91SRobert Wu         : FlowGraphFilter(channelCount) {
27*05767d91SRobert Wu }
28*05767d91SRobert Wu 
onProcess(int32_t numFrames)29*05767d91SRobert Wu int32_t Limiter::onProcess(int32_t numFrames) {
30*05767d91SRobert Wu     const float *inputBuffer = input.getBuffer();
31*05767d91SRobert Wu     float *outputBuffer = output.getBuffer();
32*05767d91SRobert Wu 
33*05767d91SRobert Wu     int32_t numSamples = numFrames * output.getSamplesPerFrame();
34*05767d91SRobert Wu 
35*05767d91SRobert Wu     // Cache the last valid output to reduce memory read/write
36*05767d91SRobert Wu     float lastValidOutput = mLastValidOutput;
37*05767d91SRobert Wu 
38*05767d91SRobert Wu     for (int32_t i = 0; i < numSamples; i++) {
39*05767d91SRobert Wu         // Use the previous output if the input is NaN
40*05767d91SRobert Wu         if (!isnan(*inputBuffer)) {
41*05767d91SRobert Wu             lastValidOutput = processFloat(*inputBuffer);
42*05767d91SRobert Wu         }
43*05767d91SRobert Wu         inputBuffer++;
44*05767d91SRobert Wu         *outputBuffer++ = lastValidOutput;
45*05767d91SRobert Wu     }
46*05767d91SRobert Wu     mLastValidOutput = lastValidOutput;
47*05767d91SRobert Wu 
48*05767d91SRobert Wu     return numFrames;
49*05767d91SRobert Wu }
50*05767d91SRobert Wu 
processFloat(float in)51*05767d91SRobert Wu float Limiter::processFloat(float in)
52*05767d91SRobert Wu {
53*05767d91SRobert Wu     float in_abs = fabsf(in);
54*05767d91SRobert Wu     if (in_abs <= 1) {
55*05767d91SRobert Wu         return in;
56*05767d91SRobert Wu     }
57*05767d91SRobert Wu     float out;
58*05767d91SRobert Wu     if (in_abs < kXWhenYis3Decibels) {
59*05767d91SRobert Wu         out = (kPolynomialSplineA * in_abs + kPolynomialSplineB) * in_abs + kPolynomialSplineC;
60*05767d91SRobert Wu     } else {
61*05767d91SRobert Wu         out = M_SQRT2;
62*05767d91SRobert Wu     }
63*05767d91SRobert Wu     if (in < 0) {
64*05767d91SRobert Wu         out = -out;
65*05767d91SRobert Wu     }
66*05767d91SRobert Wu     return out;
67*05767d91SRobert Wu }
68