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 WuLimiter::Limiter(int32_t channelCount) 26*05767d91SRobert Wu : FlowGraphFilter(channelCount) { 27*05767d91SRobert Wu } 28*05767d91SRobert Wu onProcess(int32_t numFrames)29*05767d91SRobert Wuint32_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 Wufloat 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