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 "prot.h"
      40             : #include "wmc_auto.h"
      41             : 
      42             : /*-------------------------------------------------------------------*
      43             :  * gauss_L2:
      44             :  *
      45             :  * encode an additional Gaussian excitation for unvoiced subframes and compute
      46             :  * associated xcorrelations for gains computation
      47             :  *
      48             :  * Gaussian excitation is generated by a white noise and shapes it with LPC-derived filter
      49             :  *-------------------------------------------------------------------*/
      50             : 
      51        8148 : void gauss_L2(
      52             :     const float h[],       /* i  : weighted LP filter impulse response    */
      53             :     float code[],          /* o  : gaussian excitation                    */
      54             :     float y2[],            /* i  : zero-memory filtered code. excitation  */
      55             :     float y11[],           /* o  : zero-memory filtered gauss. excitation */
      56             :     float *gain,           /* o  : excitation gain                        */
      57             :     float g_corr[],        /* i/o: correlation structure for gain coding  */
      58             :     float tilt_code,       /* i  : tilt of code                           */
      59             :     const float *Aq,       /* i  : quantized LPCs                         */
      60             :     float formant_enh_num, /* i  : formant enhancement factor             */
      61             :     int16_t *seed_acelp    /* i/o: random seed                            */
      62             : )
      63             : {
      64             :     int16_t i;
      65             : 
      66             :     /*-----------------------------------------------------------------*
      67             :      * Find new target for the Gaussian codebook
      68             :      *-----------------------------------------------------------------*/
      69             : 
      70             :     /*Generate white gaussian noise using central limit theorem method (N only 4 as E_util_random is not purely uniform)*/
      71      529620 :     for ( i = 0; i < L_SUBFR; i++ )
      72             :     {
      73      521472 :         code[i] = (float) ( own_random( seed_acelp ) ) / ( 1 << 15 );
      74      521472 :         code[i] += (float) ( own_random( seed_acelp ) ) / ( 1 << 15 );
      75      521472 :         code[i] += (float) ( own_random( seed_acelp ) ) / ( 1 << 15 );
      76             :     }
      77             : 
      78             :     /*Shape the gaussian excitation*/
      79        8148 :     cb_shape( 1, 0, 0, 1, 0, formant_enh_num, FORMANT_SHARPENING_G2, Aq, code, tilt_code, 0, L_SUBFR );
      80             : 
      81             :     /*compute 0s memory weighted synthesis contribution and find gain*/
      82        8148 :     conv( code, h, y11, L_SUBFR );
      83        8148 :     *gain = 0.f;
      84             : 
      85             :     /*Update correlations for gains coding */
      86        8148 :     g_corr[0] = 0.01F + dotp( y11, y11, L_SUBFR );
      87        8148 :     g_corr[4] = 0.01F + dotp( y11, y2, L_SUBFR );
      88             : 
      89        8148 :     return;
      90             : }