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             : #include <math.h>
      40             : #include "prot.h"
      41             : #include "wmc_auto.h"
      42             : 
      43             : /*-------------------------------------------------------------------*
      44             :  * tfaCalcEnv()
      45             :  *
      46             :  *
      47             :  *-------------------------------------------------------------------*/
      48             : 
      49         604 : void tfaCalcEnv(
      50             :     const float *shb_speech,
      51             :     float *enr )
      52             : {
      53             :     int16_t i, j, k;
      54       10268 :     for ( i = 0, k = 0; i < N_TEC_TFA_SUBFR; i++ )
      55             :     {
      56        9664 :         enr[i] = 1e-12f;
      57      202944 :         for ( j = 0; j < L_TEC_TFA_SUBFR16k; j++ ) /* XX/2 since Fs = 16kHz */
      58             :         {
      59      193280 :             enr[i] += shb_speech[k] * shb_speech[k];
      60      193280 :             k++;
      61             :         }
      62             :     }
      63             : 
      64         604 :     return;
      65             : }
      66             : 
      67             : /*-------------------------------------------------------------------*
      68             :  * tfaEnc_TBE()
      69             :  *
      70             :  *
      71             :  *-------------------------------------------------------------------*/
      72             : 
      73         604 : int16_t tfaEnc_TBE(
      74             :     const float *enr,
      75             :     const int16_t last_core,
      76             :     const float *voicing,
      77             :     const float *pitch_buf )
      78             : {
      79             :     int16_t i;
      80             :     float m_g, m_a; /* m_g: geometrical mean, m_a: arithmetical mean */
      81         604 :     float flatness = 0;
      82             :     float voicing_sum;
      83             :     int16_t tfa_flag;
      84             : 
      85             :     float pitch_buf_sum;
      86         604 :     const float m_a_bottom = 10000.0 / N_TEC_TFA_SUBFR;
      87             : 
      88         604 :     tfa_flag = 0;
      89             : 
      90         604 :     m_a = 1e-12f;
      91         604 :     m_g = 0.0;
      92       10268 :     for ( i = 0; i < N_TEC_TFA_SUBFR; i++ )
      93             :     {
      94        9664 :         m_a += enr[i];
      95        9664 :         m_g = (float) ( m_g + log10( enr[i] ) );
      96             :     }
      97         604 :     m_a /= N_TEC_TFA_SUBFR;
      98         604 :     m_g /= N_TEC_TFA_SUBFR;
      99         604 :     m_g = (float) pow( 10.0, m_g );
     100             : 
     101         604 :     flatness = m_g / m_a;
     102             : 
     103         604 :     voicing_sum = voicing[0] + voicing[1];
     104             : 
     105         604 :     pitch_buf_sum = pitch_buf[0] + pitch_buf[1] + pitch_buf[2] + pitch_buf[3];
     106             : 
     107         604 :     if ( ( flatness > 0.70 && pitch_buf_sum > 4.0 * 110 && voicing_sum > 2.0 * 0.70 ) ||
     108          25 :          ( last_core == TCX_20_CORE && flatness > 0.50 && voicing_sum < 2.0 * 0.70 ) )
     109             :     {
     110          73 :         tfa_flag = 1;
     111             :     }
     112             : 
     113             :     /* energy lower limit */
     114         604 :     if ( m_a < m_a_bottom )
     115             :     {
     116          74 :         tfa_flag = 0;
     117             :     }
     118             : 
     119         604 :     return tfa_flag;
     120             : }