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 <assert.h> 38 : #include <stdint.h> 39 : #include "options.h" 40 : #include <math.h> 41 : #include "cnst.h" 42 : #include "prot.h" 43 : #include "wmc_auto.h" 44 : 45 : /*------------------------------------------------------------------* 46 : * lsf_weight_2st() 47 : * 48 : * 49 : *------------------------------------------------------------------*/ 50 : 51 540808 : void lsf_weight_2st( 52 : const float *lsfq, 53 : float *w, 54 : const int16_t mode, 55 : const int32_t sr_core ) 56 : { 57 : int16_t i; 58 : float d[M + 1]; 59 540808 : float freq_max = sr_core / 2.f; 60 540808 : float freq_div = freq_max / (float) M; 61 : 62 : /* compute lsf distance */ 63 540808 : d[0] = lsfq[0]; 64 8652928 : for ( i = 1; i < M; i++ ) 65 : { 66 8112120 : d[i] = lsfq[i] - lsfq[i - 1]; 67 : } 68 540808 : d[M] = freq_max - lsfq[M - 1]; 69 : 70 : /* weighting function */ 71 9193736 : for ( i = 0; i < M; i++ ) 72 : { 73 8652928 : assert( d[i] > 0 ); 74 : 75 8652928 : if ( mode == 0 ) 76 : { 77 6799184 : w[i] = (float) ( 60.0f / ( freq_div / sqrt( d[i] * d[i + 1] ) ) ); /* abs */ 78 : } 79 1853744 : else if ( mode == 1 ) 80 : { 81 1776832 : w[i] = (float) ( 65.0f / ( freq_div / sqrt( d[i] * d[i + 1] ) ) ); /* mid */ 82 : } 83 : else 84 : { 85 76912 : w[i] = (float) ( 63.0f / ( freq_div / sqrt( d[i] * d[i + 1] ) ) ); /* rel2 */ 86 : } 87 : } 88 : 89 540808 : return; 90 : }