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 : * phase_dispersion()
49 : *
50 : * Post-processing to enhance noise at low bitrate
51 : *-----------------------------------------------------------------------*/
52 :
53 1630211 : void phase_dispersion(
54 : const float gain_code, /* i : gain of code */
55 : const float gain_pit, /* i : gain of pitch */
56 : float code[], /* i/o: code vector */
57 : const int16_t mode, /* i : level, 0=hi, 1=lo, 2=off */
58 : float disp_mem[] /* i/o: static memory (size = 8) */
59 : )
60 : {
61 : int16_t i, j, state;
62 : float *prev_gain_pit, *prev_gain_code, *prev_state;
63 : float code2[2 * L_SUBFR];
64 : float h_disp[L_SUBFR], *code2_real, *code2_imag, *code_real, *code_imag, *h_real, *h_imag;
65 :
66 1630211 : prev_state = disp_mem;
67 1630211 : prev_gain_code = disp_mem + 1;
68 1630211 : prev_gain_pit = disp_mem + 2;
69 :
70 1630211 : state = 2;
71 1630211 : if ( gain_pit < 0.6f )
72 : {
73 650074 : state = 0;
74 : }
75 980137 : else if ( gain_pit < 0.9f )
76 : {
77 529609 : state = 1;
78 : }
79 :
80 9781266 : for ( i = 5; i > 0; i-- )
81 : {
82 8151055 : prev_gain_pit[i] = prev_gain_pit[i - 1];
83 : }
84 1630211 : prev_gain_pit[0] = gain_pit;
85 :
86 1630211 : if ( gain_code - 3.0f * *prev_gain_code > 0.0f )
87 : {
88 52008 : if ( state < 2 )
89 : {
90 39326 : state++;
91 : }
92 : }
93 : else
94 : {
95 1578203 : j = 0;
96 11047421 : for ( i = 0; i < 6; i++ )
97 : {
98 9469218 : if ( prev_gain_pit[i] < 0.6f )
99 : {
100 3844962 : j++;
101 : }
102 : }
103 :
104 1578203 : if ( j > 2 )
105 : {
106 760119 : state = 0;
107 : }
108 :
109 1578203 : if ( ( state - (int16_t) *prev_state ) > 1 )
110 : {
111 58965 : state--;
112 : }
113 : }
114 :
115 1630211 : *prev_gain_code = gain_code;
116 1630211 : *prev_state = (float) state;
117 :
118 : /*-----------------------------------------------------------------*
119 : * Circular convolution
120 : *-----------------------------------------------------------------*/
121 :
122 1630211 : state += mode; /* level of dispersion */
123 1630211 : if ( state < 2 )
124 : {
125 108312 : fft_rel( code, L_SUBFR, 6 );
126 :
127 108312 : if ( state == 0 )
128 : {
129 29619 : mvr2r( low_H, h_disp, L_SUBFR );
130 : }
131 :
132 108312 : if ( state == 1 )
133 : {
134 78693 : mvr2r( mid_H, h_disp, L_SUBFR );
135 : }
136 :
137 108312 : code2_real = code2;
138 108312 : code2_imag = code2 + L_SUBFR - 1;
139 108312 : code_real = code;
140 108312 : code_imag = code + L_SUBFR - 1;
141 108312 : h_real = h_disp;
142 108312 : h_imag = h_disp + L_SUBFR - 1;
143 108312 : *code2_real++ = *code_real++ * *h_real++;
144 :
145 3465984 : for ( i = 1; i < L_SUBFR / 2; i++ )
146 : {
147 3357672 : *code2_real++ = *code_real * *h_real - *code_imag * *h_imag;
148 3357672 : *code2_imag-- = *code_real++ * *h_imag-- + *code_imag-- * *h_real++;
149 : }
150 :
151 108312 : *code2_real++ = *code_real++ * *h_real++;
152 108312 : ifft_rel( code2, L_SUBFR, 6 );
153 :
154 108312 : mvr2r( code2, code, L_SUBFR );
155 : }
156 :
157 1630211 : return;
158 : }
|
