FFmpeg coverage


Directory: ../../../ffmpeg/
File: src/libavcodec/aacps_tablegen.h
Date: 2024-03-28 14:59:00
Exec Total Coverage
Lines: 75 75 100.0%
Functions: 2 2 100.0%
Branches: 29 30 96.7%

Line Branch Exec Source
1 /*
2 * Header file for hardcoded Parametric Stereo tables
3 *
4 * Copyright (c) 2010 Alex Converse <alex.converse@gmail.com>
5 *
6 * This file is part of FFmpeg.
7 *
8 * FFmpeg is free software; you can redistribute it and/or
9 * modify it under the terms of the GNU Lesser General Public
10 * License as published by the Free Software Foundation; either
11 * version 2.1 of the License, or (at your option) any later version.
12 *
13 * FFmpeg is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 * Lesser General Public License for more details.
17 *
18 * You should have received a copy of the GNU Lesser General Public
19 * License along with FFmpeg; if not, write to the Free Software
20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
21 */
22
23 #ifndef AVCODEC_AACPS_TABLEGEN_H
24 #define AVCODEC_AACPS_TABLEGEN_H
25
26 #include <math.h>
27 #include <stdint.h>
28
29 #if CONFIG_HARDCODED_TABLES
30 #define ps_tableinit()
31 #define TABLE_CONST const
32 #include "libavcodec/aacps_tables.h"
33 #else
34 #include "libavutil/common.h"
35 #include "libavutil/libm.h"
36 #include "libavutil/mathematics.h"
37 #include "libavutil/mem_internal.h"
38 #define NR_ALLPASS_BANDS20 30
39 #define NR_ALLPASS_BANDS34 50
40 #define PS_AP_LINKS 3
41 #define TABLE_CONST
42 static float pd_re_smooth[8*8*8];
43 static float pd_im_smooth[8*8*8];
44 static float HA[46][8][4];
45 static float HB[46][8][4];
46 static DECLARE_ALIGNED(16, float, f20_0_8) [ 8][8][2];
47 static DECLARE_ALIGNED(16, float, f34_0_12)[12][8][2];
48 static DECLARE_ALIGNED(16, float, f34_1_8) [ 8][8][2];
49 static DECLARE_ALIGNED(16, float, f34_2_4) [ 4][8][2];
50 static TABLE_CONST DECLARE_ALIGNED(16, float, Q_fract_allpass)[2][50][3][2];
51 static DECLARE_ALIGNED(16, float, phi_fract)[2][50][2];
52
53 static const float g0_Q8[] = {
54 0.00746082949812f, 0.02270420949825f, 0.04546865930473f, 0.07266113929591f,
55 0.09885108575264f, 0.11793710567217f, 0.125f
56 };
57
58 static const float g0_Q12[] = {
59 0.04081179924692f, 0.03812810994926f, 0.05144908135699f, 0.06399831151592f,
60 0.07428313801106f, 0.08100347892914f, 0.08333333333333f
61 };
62
63 static const float g1_Q8[] = {
64 0.01565675600122f, 0.03752716391991f, 0.05417891378782f, 0.08417044116767f,
65 0.10307344158036f, 0.12222452249753f, 0.125f
66 };
67
68 static const float g2_Q4[] = {
69 -0.05908211155639f, -0.04871498374946f, 0.0f, 0.07778723915851f,
70 0.16486303567403f, 0.23279856662996f, 0.25f
71 };
72
73 644 static av_cold void make_filters_from_proto(float (*filter)[8][2], const float *proto, int bands)
74 {
75 int q, n;
76
2/2
✓ Branch 0 taken 5152 times.
✓ Branch 1 taken 644 times.
5796 for (q = 0; q < bands; q++) {
77
2/2
✓ Branch 0 taken 36064 times.
✓ Branch 1 taken 5152 times.
41216 for (n = 0; n < 7; n++) {
78 36064 double theta = 2 * M_PI * (q + 0.5) * (n - 6) / bands;
79 36064 filter[q][n][0] = proto[n] * cos(theta);
80 36064 filter[q][n][1] = proto[n] * -sin(theta);
81 }
82 }
83 644 }
84
85 161 static av_cold void ps_tableinit(void)
86 {
87 static const float ipdopd_sin[] = { 0, M_SQRT1_2, 1, M_SQRT1_2, 0, -M_SQRT1_2, -1, -M_SQRT1_2 };
88 static const float ipdopd_cos[] = { 1, M_SQRT1_2, 0, -M_SQRT1_2, -1, -M_SQRT1_2, 0, M_SQRT1_2 };
89 int pd0, pd1, pd2;
90
91 static const float iid_par_dequant[] = {
92 //iid_par_dequant_default
93 0.05623413251903, 0.12589254117942, 0.19952623149689, 0.31622776601684,
94 0.44668359215096, 0.63095734448019, 0.79432823472428, 1,
95 1.25892541179417, 1.58489319246111, 2.23872113856834, 3.16227766016838,
96 5.01187233627272, 7.94328234724282, 17.7827941003892,
97 //iid_par_dequant_fine
98 0.00316227766017, 0.00562341325190, 0.01, 0.01778279410039,
99 0.03162277660168, 0.05623413251903, 0.07943282347243, 0.11220184543020,
100 0.15848931924611, 0.22387211385683, 0.31622776601684, 0.39810717055350,
101 0.50118723362727, 0.63095734448019, 0.79432823472428, 1,
102 1.25892541179417, 1.58489319246111, 1.99526231496888, 2.51188643150958,
103 3.16227766016838, 4.46683592150963, 6.30957344480193, 8.91250938133745,
104 12.5892541179417, 17.7827941003892, 31.6227766016838, 56.2341325190349,
105 100, 177.827941003892, 316.227766016837,
106 };
107 static const float icc_invq[] = {
108 1, 0.937, 0.84118, 0.60092, 0.36764, 0, -0.589, -1
109 };
110 static const float acos_icc_invq[] = {
111 0, 0.35685527, 0.57133466, 0.92614472, 1.1943263, M_PI/2, 2.2006171, M_PI
112 };
113 int iid, icc;
114
115 int k, m;
116 static const int8_t f_center_20[] = {
117 -3, -1, 1, 3, 5, 7, 10, 14, 18, 22,
118 };
119 static const int8_t f_center_34[] = {
120 2, 6, 10, 14, 18, 22, 26, 30,
121 34,-10, -6, -2, 51, 57, 15, 21,
122 27, 33, 39, 45, 54, 66, 78, 42,
123 102, 66, 78, 90,102,114,126, 90,
124 };
125 static const float fractional_delay_links[] = { 0.43f, 0.75f, 0.347f };
126 161 const float fractional_delay_gain = 0.39f;
127
128
2/2
✓ Branch 0 taken 1288 times.
✓ Branch 1 taken 161 times.
1449 for (pd0 = 0; pd0 < 8; pd0++) {
129 1288 float pd0_re = ipdopd_cos[pd0];
130 1288 float pd0_im = ipdopd_sin[pd0];
131
2/2
✓ Branch 0 taken 10304 times.
✓ Branch 1 taken 1288 times.
11592 for (pd1 = 0; pd1 < 8; pd1++) {
132 10304 float pd1_re = ipdopd_cos[pd1];
133 10304 float pd1_im = ipdopd_sin[pd1];
134
2/2
✓ Branch 0 taken 82432 times.
✓ Branch 1 taken 10304 times.
92736 for (pd2 = 0; pd2 < 8; pd2++) {
135 82432 float pd2_re = ipdopd_cos[pd2];
136 82432 float pd2_im = ipdopd_sin[pd2];
137 82432 float re_smooth = 0.25f * pd0_re + 0.5f * pd1_re + pd2_re;
138 82432 float im_smooth = 0.25f * pd0_im + 0.5f * pd1_im + pd2_im;
139 82432 float pd_mag = 1 / hypot(im_smooth, re_smooth);
140 82432 pd_re_smooth[pd0*64+pd1*8+pd2] = re_smooth * pd_mag;
141 82432 pd_im_smooth[pd0*64+pd1*8+pd2] = im_smooth * pd_mag;
142 }
143 }
144 }
145
146
2/2
✓ Branch 0 taken 7406 times.
✓ Branch 1 taken 161 times.
7567 for (iid = 0; iid < 46; iid++) {
147 7406 float c = iid_par_dequant[iid]; ///< Linear Inter-channel Intensity Difference
148 7406 float c1 = (float)M_SQRT2 / sqrtf(1.0f + c*c);
149 7406 float c2 = c * c1;
150
2/2
✓ Branch 0 taken 59248 times.
✓ Branch 1 taken 7406 times.
66654 for (icc = 0; icc < 8; icc++) {
151 /*if (PS_BASELINE || ps->icc_mode < 3)*/ {
152 59248 float alpha = 0.5f * acos_icc_invq[icc];
153 59248 float beta = alpha * (c1 - c2) * (float)M_SQRT1_2;
154 59248 HA[iid][icc][0] = c2 * cosf(beta + alpha);
155 59248 HA[iid][icc][1] = c1 * cosf(beta - alpha);
156 59248 HA[iid][icc][2] = c2 * sinf(beta + alpha);
157 59248 HA[iid][icc][3] = c1 * sinf(beta - alpha);
158 } /* else */ {
159 float alpha, gamma, mu, rho;
160 float alpha_c, alpha_s, gamma_c, gamma_s;
161
2/2
✓ Branch 0 taken 37030 times.
✓ Branch 1 taken 22218 times.
59248 rho = FFMAX(icc_invq[icc], 0.05f);
162 59248 alpha = 0.5f * atan2f(2.0f * c * rho, c*c - 1.0f);
163 59248 mu = c + 1.0f / c;
164 59248 mu = sqrtf(1 + (4 * rho * rho - 4)/(mu * mu));
165 59248 gamma = atanf(sqrtf((1.0f - mu)/(1.0f + mu)));
166
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 59248 times.
59248 if (alpha < 0) alpha += M_PI/2;
167 59248 alpha_c = cosf(alpha);
168 59248 alpha_s = sinf(alpha);
169 59248 gamma_c = cosf(gamma);
170 59248 gamma_s = sinf(gamma);
171 59248 HB[iid][icc][0] = M_SQRT2 * alpha_c * gamma_c;
172 59248 HB[iid][icc][1] = M_SQRT2 * alpha_s * gamma_c;
173 59248 HB[iid][icc][2] = -M_SQRT2 * alpha_s * gamma_s;
174 59248 HB[iid][icc][3] = M_SQRT2 * alpha_c * gamma_s;
175 }
176 }
177 }
178
179
2/2
✓ Branch 0 taken 4830 times.
✓ Branch 1 taken 161 times.
4991 for (k = 0; k < NR_ALLPASS_BANDS20; k++) {
180 double f_center, theta;
181
2/2
✓ Branch 0 taken 1610 times.
✓ Branch 1 taken 3220 times.
4830 if (k < FF_ARRAY_ELEMS(f_center_20))
182 1610 f_center = f_center_20[k] * 0.125;
183 else
184 3220 f_center = k - 6.5f;
185
2/2
✓ Branch 0 taken 14490 times.
✓ Branch 1 taken 4830 times.
19320 for (m = 0; m < PS_AP_LINKS; m++) {
186 14490 theta = -M_PI * fractional_delay_links[m] * f_center;
187 14490 Q_fract_allpass[0][k][m][0] = cos(theta);
188 14490 Q_fract_allpass[0][k][m][1] = sin(theta);
189 }
190 4830 theta = -M_PI*fractional_delay_gain*f_center;
191 4830 phi_fract[0][k][0] = cos(theta);
192 4830 phi_fract[0][k][1] = sin(theta);
193 }
194
2/2
✓ Branch 0 taken 8050 times.
✓ Branch 1 taken 161 times.
8211 for (k = 0; k < NR_ALLPASS_BANDS34; k++) {
195 double f_center, theta;
196
2/2
✓ Branch 0 taken 5152 times.
✓ Branch 1 taken 2898 times.
8050 if (k < FF_ARRAY_ELEMS(f_center_34))
197 5152 f_center = f_center_34[k] / 24.0;
198 else
199 2898 f_center = k - 26.5f;
200
2/2
✓ Branch 0 taken 24150 times.
✓ Branch 1 taken 8050 times.
32200 for (m = 0; m < PS_AP_LINKS; m++) {
201 24150 theta = -M_PI * fractional_delay_links[m] * f_center;
202 24150 Q_fract_allpass[1][k][m][0] = cos(theta);
203 24150 Q_fract_allpass[1][k][m][1] = sin(theta);
204 }
205 8050 theta = -M_PI*fractional_delay_gain*f_center;
206 8050 phi_fract[1][k][0] = cos(theta);
207 8050 phi_fract[1][k][1] = sin(theta);
208 }
209
210 161 make_filters_from_proto(f20_0_8, g0_Q8, 8);
211 161 make_filters_from_proto(f34_0_12, g0_Q12, 12);
212 161 make_filters_from_proto(f34_1_8, g1_Q8, 8);
213 161 make_filters_from_proto(f34_2_4, g2_Q4, 4);
214 161 }
215 #endif /* CONFIG_HARDCODED_TABLES */
216
217 #endif /* AVCODEC_AACPS_TABLEGEN_H */
218