summaryrefslogtreecommitdiffstats
path: root/media/libstagefright/codecs/aacdec/calc_sbr_anafilterbank.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'media/libstagefright/codecs/aacdec/calc_sbr_anafilterbank.cpp')
-rw-r--r--media/libstagefright/codecs/aacdec/calc_sbr_anafilterbank.cpp360
1 files changed, 360 insertions, 0 deletions
diff --git a/media/libstagefright/codecs/aacdec/calc_sbr_anafilterbank.cpp b/media/libstagefright/codecs/aacdec/calc_sbr_anafilterbank.cpp
new file mode 100644
index 0000000..5ec3f69
--- /dev/null
+++ b/media/libstagefright/codecs/aacdec/calc_sbr_anafilterbank.cpp
@@ -0,0 +1,360 @@
+/* ------------------------------------------------------------------
+ * Copyright (C) 1998-2009 PacketVideo
+ *
+ * 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.
+ * -------------------------------------------------------------------
+ */
+/*
+
+ Filename: calc_sbr_anafilterbank.c
+
+
+------------------------------------------------------------------------------
+ REVISION HISTORY
+
+
+ Who: Date: MM/DD/YYYY
+ Description:
+
+------------------------------------------------------------------------------
+ INPUT AND OUTPUT DEFINITIONS
+
+
+
+------------------------------------------------------------------------------
+ FUNCTION DESCRIPTION
+
+
+------------------------------------------------------------------------------
+ REQUIREMENTS
+
+
+------------------------------------------------------------------------------
+ REFERENCES
+
+SC 29 Software Copyright Licencing Disclaimer:
+
+This software module was originally developed by
+ Coding Technologies
+
+and edited by
+ -
+
+in the course of development of the ISO/IEC 13818-7 and ISO/IEC 14496-3
+standards for reference purposes and its performance may not have been
+optimized. This software module is an implementation of one or more tools as
+specified by the ISO/IEC 13818-7 and ISO/IEC 14496-3 standards.
+ISO/IEC gives users free license to this software module or modifications
+thereof for use in products claiming conformance to audiovisual and
+image-coding related ITU Recommendations and/or ISO/IEC International
+Standards. ISO/IEC gives users the same free license to this software module or
+modifications thereof for research purposes and further ISO/IEC standardisation.
+Those intending to use this software module in products are advised that its
+use may infringe existing patents. ISO/IEC have no liability for use of this
+software module or modifications thereof. Copyright is not released for
+products that do not conform to audiovisual and image-coding related ITU
+Recommendations and/or ISO/IEC International Standards.
+The original developer retains full right to modify and use the code for its
+own purpose, assign or donate the code to a third party and to inhibit third
+parties from using the code for products that do not conform to audiovisual and
+image-coding related ITU Recommendations and/or ISO/IEC International Standards.
+This copyright notice must be included in all copies or derivative works.
+Copyright (c) ISO/IEC 2002.
+
+------------------------------------------------------------------------------
+ PSEUDO-CODE
+
+------------------------------------------------------------------------------
+*/
+
+
+/*----------------------------------------------------------------------------
+; INCLUDES
+----------------------------------------------------------------------------*/
+
+#ifdef AAC_PLUS
+
+
+#include "calc_sbr_anafilterbank.h"
+#include "qmf_filterbank_coeff.h"
+#include "analysis_sub_band.h"
+
+#include "aac_mem_funcs.h"
+#include "fxp_mul32.h"
+
+
+
+/*----------------------------------------------------------------------------
+; MACROS
+; Define module specific macros here
+----------------------------------------------------------------------------*/
+
+
+/*----------------------------------------------------------------------------
+; DEFINES
+; Include all pre-processor statements here. Include conditional
+; compile variables also.
+----------------------------------------------------------------------------*/
+
+/*----------------------------------------------------------------------------
+; LOCAL FUNCTION DEFINITIONS
+; Function Prototype declaration
+----------------------------------------------------------------------------*/
+
+/*----------------------------------------------------------------------------
+; LOCAL STORE/BUFFER/POINTER DEFINITIONS
+; Variable declaration - defined here and used outside this module
+----------------------------------------------------------------------------*/
+
+/*----------------------------------------------------------------------------
+; EXTERNAL FUNCTION REFERENCES
+; Declare functions defined elsewhere and referenced in this module
+----------------------------------------------------------------------------*/
+
+/*----------------------------------------------------------------------------
+; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES
+; Declare variables used in this module but defined elsewhere
+----------------------------------------------------------------------------*/
+
+/*----------------------------------------------------------------------------
+; FUNCTION CODE
+----------------------------------------------------------------------------*/
+
+void calc_sbr_anafilterbank_LC(Int32 * Sr,
+ Int16 * X,
+ Int32 scratch_mem[][64],
+ Int32 maxBand)
+{
+
+ Int i;
+ Int32 *p_Y_1;
+ Int32 *p_Y_2;
+
+ Int16 * pt_X_1;
+ Int16 * pt_X_2;
+ Int32 realAccu1;
+ Int32 realAccu2;
+
+ Int32 tmp1;
+ Int32 tmp2;
+
+
+ const Int32 * pt_C;
+
+ p_Y_1 = scratch_mem[0];
+
+
+ p_Y_2 = p_Y_1 + 63;
+ pt_C = &sbrDecoderFilterbankCoefficients_an_filt_LC[0];
+
+ pt_X_1 = X;
+
+
+ realAccu1 = fxp_mul32_by_16(Qfmt27(-0.51075594183097F), pt_X_1[-192]);
+
+ realAccu1 = fxp_mac32_by_16(Qfmt27(-0.51075594183097F), -pt_X_1[-128], realAccu1);
+ realAccu1 = fxp_mac32_by_16(Qfmt27(-0.01876919066980F), pt_X_1[-256], realAccu1);
+ *(p_Y_1++) = fxp_mac32_by_16(Qfmt27(-0.01876919066980F), -pt_X_1[ -64], realAccu1);
+
+
+ /* create array Y */
+
+ pt_X_1 = &X[-1];
+ pt_X_2 = &X[-319];
+
+
+ for (i = 15; i != 0; i--)
+ {
+ tmp1 = *(pt_X_1--);
+ tmp2 = *(pt_X_2++);
+
+ realAccu1 = fxp_mul32_by_16(*(pt_C), tmp1);
+ realAccu2 = fxp_mul32_by_16(*(pt_C++), tmp2);
+ tmp1 = pt_X_1[ -63];
+ tmp2 = pt_X_2[ +63];
+ realAccu1 = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ realAccu2 = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ tmp1 = pt_X_1[ -127];
+ tmp2 = pt_X_2[ +127];
+ realAccu1 = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ realAccu2 = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ tmp1 = pt_X_1[ -191];
+ tmp2 = pt_X_2[ +191];
+ realAccu1 = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ realAccu2 = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ tmp1 = pt_X_1[ -255];
+ tmp2 = pt_X_2[ +255];
+ *(p_Y_1++) = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ *(p_Y_2--) = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+
+ tmp1 = *(pt_X_1--);
+ tmp2 = *(pt_X_2++);
+ realAccu1 = fxp_mul32_by_16(*(pt_C), tmp1);
+ realAccu2 = fxp_mul32_by_16(*(pt_C++), tmp2);
+
+ tmp1 = pt_X_1[ -63];
+ tmp2 = pt_X_2[ +63];
+ realAccu1 = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ realAccu2 = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ tmp1 = pt_X_1[ -127];
+ tmp2 = pt_X_2[ +127];
+ realAccu1 = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ realAccu2 = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ tmp1 = pt_X_1[ -191];
+ tmp2 = pt_X_2[ +191];
+ realAccu1 = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ realAccu2 = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ tmp1 = pt_X_1[ -255];
+ tmp2 = pt_X_2[ +255];
+ *(p_Y_1++) = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ *(p_Y_2--) = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+
+ }
+
+
+ tmp1 = *(pt_X_1--);
+ tmp2 = *(pt_X_2++);
+ realAccu1 = fxp_mul32_by_16(*(pt_C), tmp1);
+ realAccu2 = fxp_mul32_by_16(*(pt_C++), tmp2);
+
+ tmp1 = pt_X_1[ -63];
+ tmp2 = pt_X_2[ +63];
+ realAccu1 = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ realAccu2 = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ tmp1 = pt_X_1[ -127];
+ tmp2 = pt_X_2[ +127];
+ realAccu1 = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ realAccu2 = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ tmp1 = pt_X_1[ -191];
+ tmp2 = pt_X_2[ +191];
+ realAccu1 = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ realAccu2 = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ tmp1 = pt_X_1[ -255];
+ tmp2 = pt_X_2[ +255];
+ *(p_Y_1++) = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ *(p_Y_2--) = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+
+
+ pt_X_1 = X;
+
+ realAccu2 = fxp_mul32_by_16(Qfmt27(0.00370548843500F), X[ -32]);
+
+ realAccu2 = fxp_mac32_by_16(Qfmt27(0.00370548843500F), pt_X_1[-288], realAccu2);
+ realAccu2 = fxp_mac32_by_16(Qfmt27(0.09949460091720F), pt_X_1[ -96], realAccu2);
+ realAccu2 = fxp_mac32_by_16(Qfmt27(0.09949460091720F), pt_X_1[-224], realAccu2);
+ *(p_Y_1++) = fxp_mac32_by_16(Qfmt27(1.20736865027288F), pt_X_1[-160], realAccu2);
+
+
+ analysis_sub_band_LC(scratch_mem[0],
+ Sr,
+ maxBand,
+ (Int32(*)[64])scratch_mem[1]);
+
+}
+
+
+
+#ifdef HQ_SBR
+
+void calc_sbr_anafilterbank(Int32 * Sr,
+ Int32 * Si,
+ Int16 * X,
+ Int32 scratch_mem[][64],
+ Int32 maxBand)
+{
+ Int i;
+ Int32 *p_Y_1;
+ Int32 *p_Y_2;
+
+
+
+
+ const Int32 * pt_C;
+ Int16 * pt_X_1;
+ Int16 * pt_X_2;
+ Int32 realAccu1;
+ Int32 realAccu2;
+
+ Int32 tmp1;
+ Int32 tmp2;
+
+
+ p_Y_1 = scratch_mem[0];
+
+
+ p_Y_2 = p_Y_1 + 63;
+ pt_C = &sbrDecoderFilterbankCoefficients_an_filt[0];
+
+ realAccu1 = fxp_mul32_by_16(Qfmt27(-0.36115899F), X[-192]);
+
+
+ realAccu1 = fxp_mac32_by_16(Qfmt27(-0.36115899F), -X[-128], realAccu1);
+ realAccu1 = fxp_mac32_by_16(Qfmt27(-0.013271822F), X[-256], realAccu1);
+ *(p_Y_1++) = fxp_mac32_by_16(Qfmt27(-0.013271822F), -X[ -64], realAccu1);
+
+ /* create array Y */
+
+ pt_X_1 = &X[-1];
+ pt_X_2 = &X[-319];
+
+
+ for (i = 31; i != 0; i--)
+ {
+ tmp1 = *(pt_X_1--);
+ tmp2 = *(pt_X_2++);
+ realAccu1 = fxp_mul32_by_16(*(pt_C), tmp1);
+ realAccu2 = fxp_mul32_by_16(*(pt_C++), tmp2);
+ tmp1 = pt_X_1[ -63];
+ tmp2 = pt_X_2[ 63];
+ realAccu1 = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ realAccu2 = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ tmp1 = pt_X_1[ -127];
+ tmp2 = pt_X_2[ 127];
+ realAccu1 = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ realAccu2 = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ tmp1 = pt_X_1[ -191];
+ tmp2 = pt_X_2[ 191];
+ realAccu1 = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ realAccu2 = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ tmp1 = pt_X_1[ -255];
+ tmp2 = pt_X_2[ 255];
+ *(p_Y_1++) = fxp_mac32_by_16(*(pt_C), tmp1, realAccu1);
+ *(p_Y_2--) = fxp_mac32_by_16(*(pt_C++), tmp2, realAccu2);
+ }
+
+
+ realAccu2 = fxp_mul32_by_16(Qfmt27(0.002620176F), X[ -32]);
+ realAccu2 = fxp_mac32_by_16(Qfmt27(0.002620176F), X[-288], realAccu2);
+ realAccu2 = fxp_mac32_by_16(Qfmt27(0.070353307F), X[ -96], realAccu2);
+ realAccu2 = fxp_mac32_by_16(Qfmt27(0.070353307F), X[-224], realAccu2);
+
+
+ *(p_Y_1++) = fxp_mac32_by_16(Qfmt27(0.85373856F), (X[-160]), realAccu2);
+
+
+ analysis_sub_band(scratch_mem[0],
+ Sr,
+ Si,
+ maxBand,
+ (Int32(*)[64])scratch_mem[1]);
+
+}
+
+
+#endif
+
+
+
+#endif /* AAC_PLUS */
+
generated by cgit v1.1 at 2024-05-29 20:33:49 +0000