Line data    Source code

       1             : /******************************************************************************************************
       2             : 
       3             :    (C) 2022-2025 IVAS codec Public Collaboration with portions copyright Dolby International AB, Ericsson AB,
       4             :    Fraunhofer-Gesellschaft zur Foerderung der angewandten Forschung e.V., Huawei Technologies Co. LTD.,
       5             :    Koninklijke Philips N.V., Nippon Telegraph and Telephone Corporation, Nokia Technologies Oy, Orange,
       6             :    Panasonic Holdings Corporation, Qualcomm Technologies, Inc., VoiceAge Corporation, and other
       7             :    contributors to this repository. All Rights Reserved.
       8             : 
       9             :    This software is protected by copyright law and by international treaties.
      10             :    The IVAS codec Public Collaboration consisting of Dolby International AB, Ericsson AB,
      11             :    Fraunhofer-Gesellschaft zur Foerderung der angewandten Forschung e.V., Huawei Technologies Co. LTD.,
      12             :    Koninklijke Philips N.V., Nippon Telegraph and Telephone Corporation, Nokia Technologies Oy, Orange,
      13             :    Panasonic Holdings Corporation, Qualcomm Technologies, Inc., VoiceAge Corporation, and other
      14             :    contributors to this repository retain full ownership rights in their respective contributions in
      15             :    the software. This notice grants no license of any kind, including but not limited to patent
      16             :    license, nor is any license granted by implication, estoppel or otherwise.
      17             : 
      18             :    Contributors are required to enter into the IVAS codec Public Collaboration agreement before making
      19             :    contributions.
      20             : 
      21             :    This software is provided "AS IS", without any express or implied warranties. The software is in the
      22             :    development stage. It is intended exclusively for experts who have experience with such software and
      23             :    solely for the purpose of inspection. All implied warranties of non-infringement, merchantability
      24             :    and fitness for a particular purpose are hereby disclaimed and excluded.
      25             : 
      26             :    Any dispute, controversy or claim arising under or in relation to providing this software shall be
      27             :    submitted to and settled by the final, binding jurisdiction of the courts of Munich, Germany in
      28             :    accordance with the laws of the Federal Republic of Germany excluding its conflict of law rules and
      29             :    the United Nations Convention on Contracts on the International Sales of Goods.
      30             : 
      31             : *******************************************************************************************************/
      32             : 
      33             : /*====================================================================================
      34             :     EVS Codec 3GPP TS26.443 Nov 04, 2021. Version 12.14.0 / 13.10.0 / 14.6.0 / 15.4.0 / 16.3.0
      35             :   ====================================================================================*/
      36             : 
      37             : #include <stdint.h>
      38             : #include "options.h"
      39             : #ifdef DEBUGGING
      40             : #include "debug.h"
      41             : #endif
      42             : #include "cnst.h"
      43             : #include "prot.h"
      44             : #include "rom_com.h"
      45             : #include "wmc_auto.h"
      46             : 
      47             : /*-------------------------------------------------------------------*
      48             :  * Interpol_delay()
      49             :  *
      50             :  * Interpolation of pitch lag
      51             :  *--------------------------------------------------------------------*/
      52             : 
      53           0 : void Interpol_delay(
      54             :     float *out,       /* o  : pitch interpolation output   */
      55             :     float *last,      /* i  : last frame pitch lag         */
      56             :     float *current,   /* i  : current frame pitch lag      */
      57             :     int16_t SubNum,   /* i  : subframe number              */
      58             :     const float *frac /* i  : interpolation constant       */
      59             : )
      60             : {
      61           0 :     out[0] = ( 1.0f - frac[SubNum] ) * *last + frac[SubNum] * *current;
      62           0 :     out[1] = ( 1.0f - frac[SubNum + 1] ) * *last + frac[SubNum + 1] * *current;
      63           0 :     out[2] = ( 1.0f - frac[SubNum + 2] ) * *last + frac[SubNum + 2] * *current;
      64             : 
      65           0 :     return;
      66             : }
      67             : 
      68             : 
      69             : /*-------------------------------------------------------------------*
      70             :  * deemph_lpc()
      71             :  *
      72             :  * De-emphasis of LP coefficients
      73             :  * convolve LPC with [1 -PREEMPH_FAC] to de-emphasise LPC
      74             :  *--------------------------------------------------------------------*/
      75             : 
      76           0 : void deemph_lpc(
      77             :     float *p_Aq_curr,   /* i  : LP coefficients current frame                 */
      78             :     float *p_Aq_old,    /* i  : LP coefficients previous frame                */
      79             :     float *LPC_de_curr, /* o  : De-emphasized LP coefficients current frame   */
      80             :     float *LPC_de_old,  /* o  : De-emphasized LP coefficients previous frame  */
      81             :     const int16_t deemph_old )
      82             : {
      83             :     int16_t k;
      84             :     float b[M + 2];
      85             :     const float a[2] = { -PREEMPH_FAC, 1.0 };
      86             : 
      87           0 :     b[0] = 1.0;
      88           0 :     for ( k = 0; k < M; k++ )
      89             :     {
      90           0 :         b[k + 1] = p_Aq_curr[k];
      91             :     }
      92           0 :     b[M + 1] = 0.0;
      93             : 
      94           0 :     for ( k = 0; k <= M; k++ )
      95             :     {
      96           0 :         LPC_de_curr[k] = a[0] * b[k] + a[1] * b[k + 1];
      97             :     }
      98             : 
      99           0 :     if ( deemph_old == 1 )
     100             :     {
     101             :         /* ignoring the 1st value which is 1.0 in this case */
     102           0 :         b[0] = 1.0;
     103           0 :         for ( k = 0; k < M; k++ )
     104             :         {
     105           0 :             b[k + 1] = p_Aq_old[k + 1];
     106             :         }
     107           0 :         b[M + 1] = 0.0;
     108             : 
     109           0 :         for ( k = 0; k <= M; k++ )
     110             :         {
     111           0 :             LPC_de_old[k] = a[0] * b[k] + a[1] * b[k + 1];
     112             :         }
     113             :     }
     114             : 
     115           0 :     return;
     116             : }