blob: d63365c6090472cad0bca411a5227b7c0c497931 [file] [log] [blame]
/*
* Copyright (C) 2004-2010 NXP Software
* Copyright (C) 2010 The Android Open Source Project
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
#include "BIQUAD.h"
#include "BQ_2I_D32F32Cll_TRC_WRA_01_Private.h"
#include "LVM_Macros.h"
/**************************************************************************
ASSUMPTIONS:
COEFS-
pBiquadState->coefs[0] is A2, pBiquadState->coefs[1] is A1
pBiquadState->coefs[2] is A0, pBiquadState->coefs[3] is -B2
pBiquadState->coefs[4] is -B1, these are in Q30 format
DELAYS-
pBiquadState->pDelays[0] is x(n-1)L in Q0 format
pBiquadState->pDelays[1] is x(n-1)R in Q0 format
pBiquadState->pDelays[2] is x(n-2)L in Q0 format
pBiquadState->pDelays[3] is x(n-2)R in Q0 format
pBiquadState->pDelays[4] is y(n-1)L in Q0 format
pBiquadState->pDelays[5] is y(n-1)R in Q0 format
pBiquadState->pDelays[6] is y(n-2)L in Q0 format
pBiquadState->pDelays[7] is y(n-2)R in Q0 format
***************************************************************************/
#ifdef BUILD_FLOAT
void BQ_2I_D32F32C30_TRC_WRA_01 ( Biquad_FLOAT_Instance_t *pInstance,
LVM_FLOAT *pDataIn,
LVM_FLOAT *pDataOut,
LVM_INT16 NrSamples)
{
LVM_FLOAT ynL,ynR,templ,tempd;
LVM_INT16 ii;
PFilter_State_FLOAT pBiquadState = (PFilter_State_FLOAT) pInstance;
for (ii = NrSamples; ii != 0; ii--)
{
/**************************************************************************
PROCESSING OF THE LEFT CHANNEL
***************************************************************************/
/* ynL= ( A2 * x(n-2)L ) */
ynL = pBiquadState->coefs[0] * pBiquadState->pDelays[2];
/* ynL+= ( A1 * x(n-1)L )*/
templ = pBiquadState->coefs[1] * pBiquadState->pDelays[0];
ynL += templ;
/* ynL+= ( A0 * x(n)L ) */
templ = pBiquadState->coefs[2] * (*pDataIn);
ynL += templ;
/* ynL+= (-B2 * y(n-2)L ) */
templ = pBiquadState->coefs[3] * pBiquadState->pDelays[6];
ynL += templ;
/* ynL+= (-B1 * y(n-1)L )*/
templ = pBiquadState->coefs[4] * pBiquadState->pDelays[4];
ynL += templ;
/**************************************************************************
PROCESSING OF THE RIGHT CHANNEL
***************************************************************************/
/* ynR= ( A2 * x(n-2)R ) */
ynR = pBiquadState->coefs[0] * pBiquadState->pDelays[3];
/* ynR+= ( A1 * x(n-1)R ) */
templ = pBiquadState->coefs[1] * pBiquadState->pDelays[1];
ynR += templ;
/* ynR+= ( A0 * x(n)R ) */
tempd =* (pDataIn+1);
templ = pBiquadState->coefs[2] * tempd;
ynR += templ;
/* ynR+= (-B2 * y(n-2)R ) */
templ = pBiquadState->coefs[3] * pBiquadState->pDelays[7];
ynR += templ;
/* ynR+= (-B1 * y(n-1)R ) */
templ = pBiquadState->coefs[4] * pBiquadState->pDelays[5];
ynR += templ;
/**************************************************************************
UPDATING THE DELAYS
***************************************************************************/
pBiquadState->pDelays[7] = pBiquadState->pDelays[5]; /* y(n-2)R=y(n-1)R*/
pBiquadState->pDelays[6] = pBiquadState->pDelays[4]; /* y(n-2)L=y(n-1)L*/
pBiquadState->pDelays[3] = pBiquadState->pDelays[1]; /* x(n-2)R=x(n-1)R*/
pBiquadState->pDelays[2] = pBiquadState->pDelays[0]; /* x(n-2)L=x(n-1)L*/
pBiquadState->pDelays[5] = (LVM_FLOAT)ynR; /* Update y(n-1)R */
pBiquadState->pDelays[4] = (LVM_FLOAT)ynL; /* Update y(n-1)L */
pBiquadState->pDelays[0] = (*pDataIn); /* Update x(n-1)L */
pDataIn++;
pBiquadState->pDelays[1] = (*pDataIn); /* Update x(n-1)R */
pDataIn++;
/**************************************************************************
WRITING THE OUTPUT
***************************************************************************/
*pDataOut = (LVM_FLOAT)ynL; /* Write Left output */
pDataOut++;
*pDataOut = (LVM_FLOAT)ynR; /* Write Right ouput */
pDataOut++;
}
}
#ifdef SUPPORT_MC
/**************************************************************************
ASSUMPTIONS:
COEFS-
pBiquadState->coefs[0] is A2, pBiquadState->coefs[1] is A1
pBiquadState->coefs[2] is A0, pBiquadState->coefs[3] is -B2
pBiquadState->coefs[4] is -B1
DELAYS-
pBiquadState->pDelays[0] to
pBiquadState->pDelays[NrChannels - 1] is x(n-1) for all NrChannels
pBiquadState->pDelays[NrChannels] to
pBiquadState->pDelays[2*NrChannels - 1] is x(n-2) for all NrChannels
pBiquadState->pDelays[2*NrChannels] to
pBiquadState->pDelays[3*NrChannels - 1] is y(n-1) for all NrChannels
pBiquadState->pDelays[3*NrChannels] to
pBiquadState->pDelays[4*NrChannels - 1] is y(n-2) for all NrChannels
***************************************************************************/
void BQ_MC_D32F32C30_TRC_WRA_01 ( Biquad_FLOAT_Instance_t *pInstance,
LVM_FLOAT *pDataIn,
LVM_FLOAT *pDataOut,
LVM_INT16 NrFrames,
LVM_INT16 NrChannels)
{
LVM_FLOAT yn, temp;
LVM_INT16 ii, jj;
PFilter_State_FLOAT pBiquadState = (PFilter_State_FLOAT) pInstance;
for (ii = NrFrames; ii != 0; ii--)
{
/**************************************************************************
PROCESSING CHANNEL-WISE
***************************************************************************/
for (jj = 0; jj < NrChannels; jj++)
{
/* yn= (A2 * x(n-2)) */
yn = pBiquadState->coefs[0] * pBiquadState->pDelays[NrChannels + jj];
/* yn+= (A1 * x(n-1)) */
temp = pBiquadState->coefs[1] * pBiquadState->pDelays[jj];
yn += temp;
/* yn+= (A0 * x(n)) */
temp = pBiquadState->coefs[2] * (*pDataIn);
yn += temp;
/* yn+= (-B2 * y(n-2)) */
temp = pBiquadState->coefs[3] * pBiquadState->pDelays[NrChannels*3 + jj];
yn += temp;
/* yn+= (-B1 * y(n-1)) */
temp = pBiquadState->coefs[4] * pBiquadState->pDelays[NrChannels*2 + jj];
yn += temp;
/**************************************************************************
UPDATING THE DELAYS
***************************************************************************/
pBiquadState->pDelays[NrChannels * 3 + jj] =
pBiquadState->pDelays[NrChannels * 2 + jj]; /* y(n-2)=y(n-1)*/
pBiquadState->pDelays[NrChannels * 1 + jj] =
pBiquadState->pDelays[jj]; /* x(n-2)=x(n-1)*/
pBiquadState->pDelays[NrChannels * 2 + jj] = (LVM_FLOAT)yn; /* Update y(n-1)*/
pBiquadState->pDelays[jj] = (*pDataIn); /* Update x(n-1)*/
pDataIn++;
/**************************************************************************
WRITING THE OUTPUT
***************************************************************************/
*pDataOut = (LVM_FLOAT)yn; /* Write jj Channel output */
pDataOut++;
}
}
}
#endif /*SUPPORT_MC*/
#else
void BQ_2I_D32F32C30_TRC_WRA_01 ( Biquad_Instance_t *pInstance,
LVM_INT32 *pDataIn,
LVM_INT32 *pDataOut,
LVM_INT16 NrSamples)
{
LVM_INT32 ynL,ynR,templ,tempd;
LVM_INT16 ii;
PFilter_State pBiquadState = (PFilter_State) pInstance;
for (ii = NrSamples; ii != 0; ii--)
{
/**************************************************************************
PROCESSING OF THE LEFT CHANNEL
***************************************************************************/
/* ynL= ( A2 (Q30) * x(n-2)L (Q0) ) >>30 in Q0*/
MUL32x32INTO32(pBiquadState->coefs[0],pBiquadState->pDelays[2],ynL,30)
/* ynL+= ( A1 (Q30) * x(n-1)L (Q0) ) >> 30 in Q0*/
MUL32x32INTO32(pBiquadState->coefs[1],pBiquadState->pDelays[0],templ,30)
ynL+=templ;
/* ynL+= ( A0 (Q30) * x(n)L (Q0) ) >> 30 in Q0*/
MUL32x32INTO32(pBiquadState->coefs[2],*pDataIn,templ,30)
ynL+=templ;
/* ynL+= (-B2 (Q30) * y(n-2)L (Q0) ) >> 30 in Q0*/
MUL32x32INTO32(pBiquadState->coefs[3],pBiquadState->pDelays[6],templ,30)
ynL+=templ;
/* ynL+= (-B1 (Q30) * y(n-1)L (Q0) ) >> 30 in Q0 */
MUL32x32INTO32(pBiquadState->coefs[4],pBiquadState->pDelays[4],templ,30)
ynL+=templ;
/**************************************************************************
PROCESSING OF THE RIGHT CHANNEL
***************************************************************************/
/* ynR= ( A2 (Q30) * x(n-2)R (Q0) ) >> 30 in Q0*/
MUL32x32INTO32(pBiquadState->coefs[0],pBiquadState->pDelays[3],ynR,30)
/* ynR+= ( A1 (Q30) * x(n-1)R (Q0) ) >> 30 in Q0*/
MUL32x32INTO32(pBiquadState->coefs[1],pBiquadState->pDelays[1],templ,30)
ynR+=templ;
/* ynR+= ( A0 (Q30) * x(n)R (Q0) ) >> 30 in Q0*/
tempd=*(pDataIn+1);
MUL32x32INTO32(pBiquadState->coefs[2],tempd,templ,30)
ynR+=templ;
/* ynR+= (-B2 (Q30) * y(n-2)R (Q0) ) >> 30 in Q0*/
MUL32x32INTO32(pBiquadState->coefs[3],pBiquadState->pDelays[7],templ,30)
ynR+=templ;
/* ynR+= (-B1 (Q30) * y(n-1)R (Q0) ) >> 30 in Q0 */
MUL32x32INTO32(pBiquadState->coefs[4],pBiquadState->pDelays[5],templ,30)
ynR+=templ;
/**************************************************************************
UPDATING THE DELAYS
***************************************************************************/
pBiquadState->pDelays[7]=pBiquadState->pDelays[5]; /* y(n-2)R=y(n-1)R*/
pBiquadState->pDelays[6]=pBiquadState->pDelays[4]; /* y(n-2)L=y(n-1)L*/
pBiquadState->pDelays[3]=pBiquadState->pDelays[1]; /* x(n-2)R=x(n-1)R*/
pBiquadState->pDelays[2]=pBiquadState->pDelays[0]; /* x(n-2)L=x(n-1)L*/
pBiquadState->pDelays[5]=(LVM_INT32)ynR; /* Update y(n-1)R in Q0*/
pBiquadState->pDelays[4]=(LVM_INT32)ynL; /* Update y(n-1)L in Q0*/
pBiquadState->pDelays[0]=(*pDataIn); /* Update x(n-1)L in Q0*/
pDataIn++;
pBiquadState->pDelays[1]=(*pDataIn); /* Update x(n-1)R in Q0*/
pDataIn++;
/**************************************************************************
WRITING THE OUTPUT
***************************************************************************/
*pDataOut=(LVM_INT32)ynL; /* Write Left output in Q0*/
pDataOut++;
*pDataOut=(LVM_INT32)ynR; /* Write Right ouput in Q0*/
pDataOut++;
}
}
#endif /*BUILD_FLOAT*/