1/*
2 * Header file for hardcoded QDM2 tables
3 *
4 * Copyright (c) 2010 Reimar Döffinger <Reimar.Doeffinger@gmx.de>
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_QDM2_TABLEGEN_H
24#define AVCODEC_QDM2_TABLEGEN_H
25
26#include <stdint.h>
27#include <math.h>
28#include "libavutil/attributes.h"
29#include "qdm2data.h"
30
31#define SOFTCLIP_THRESHOLD 27600
32#define HARDCLIP_THRESHOLD 35716
33
34#if CONFIG_HARDCODED_TABLES
35#define softclip_table_init()
36#define rnd_table_init()
37#define init_noise_samples()
38#define qdm2_init_vlc()
39#include "libavcodec/qdm2_tables.h"
40#else
41static uint16_t softclip_table[HARDCLIP_THRESHOLD - SOFTCLIP_THRESHOLD + 1];
42static float noise_table[4096 + 20];
43static uint8_t random_dequant_index[256][5];
44static uint8_t random_dequant_type24[128][3];
45static float noise_samples[128];
46
47static av_cold void softclip_table_init(void) {
48    int i;
49    double dfl = SOFTCLIP_THRESHOLD - 32767;
50    float delta = 1.0 / -dfl;
51    for (i = 0; i < HARDCLIP_THRESHOLD - SOFTCLIP_THRESHOLD + 1; i++)
52        softclip_table[i] = SOFTCLIP_THRESHOLD - ((int)(sin((float)i * delta) * dfl) & 0x0000FFFF);
53}
54
55
56// random generated table
57static av_cold void rnd_table_init(void) {
58    int i,j;
59    uint32_t ldw;
60    uint64_t random_seed = 0;
61    float delta = 1.0 / 16384.0;
62    for(i = 0; i < 4096 ;i++) {
63        random_seed = random_seed * 214013 + 2531011;
64        noise_table[i] = (delta * (float)(((int32_t)random_seed >> 16) & 0x00007FFF)- 1.0) * 1.3;
65    }
66
67    for (i = 0; i < 256 ;i++) {
68        random_seed = 81;
69        ldw = i;
70        for (j = 0; j < 5 ;j++) {
71            random_dequant_index[i][j] = ldw / random_seed;
72            ldw %= random_seed;
73            random_seed /= 3;
74        }
75    }
76    for (i = 0; i < 128 ;i++) {
77        random_seed = 25;
78        ldw = i;
79        for (j = 0; j < 3 ;j++) {
80            random_dequant_type24[i][j] = ldw / random_seed;
81            ldw %= random_seed;
82            random_seed /= 5;
83        }
84    }
85}
86
87
88static av_cold void init_noise_samples(void) {
89    int i;
90    unsigned random_seed = 0;
91    float delta = 1.0 / 16384.0;
92    for (i = 0; i < 128;i++) {
93        random_seed = random_seed * 214013 + 2531011;
94        noise_samples[i] = (delta * (float)((random_seed >> 16) & 0x00007fff) - 1.0);
95    }
96}
97
98static VLC vlc_tab_level;
99static VLC vlc_tab_diff;
100static VLC vlc_tab_run;
101static VLC fft_level_exp_alt_vlc;
102static VLC fft_level_exp_vlc;
103static VLC fft_stereo_exp_vlc;
104static VLC fft_stereo_phase_vlc;
105static VLC vlc_tab_tone_level_idx_hi1;
106static VLC vlc_tab_tone_level_idx_mid;
107static VLC vlc_tab_tone_level_idx_hi2;
108static VLC vlc_tab_type30;
109static VLC vlc_tab_type34;
110static VLC vlc_tab_fft_tone_offset[5];
111
112static VLCElem qdm2_table[3838];
113
114static av_cold void build_vlc(VLC *vlc, int nb_bits, int nb_codes,
115                              unsigned *offset, const uint8_t tab[][2])
116{
117    vlc->table           = &qdm2_table[*offset];
118    vlc->table_allocated = FF_ARRAY_ELEMS(qdm2_table) - *offset;
119    ff_init_vlc_from_lengths(vlc, nb_bits, nb_codes,
120                             &tab[0][1], 2, &tab[0][0], 2, 1,
121                             -1, INIT_VLC_STATIC_OVERLONG | INIT_VLC_LE, NULL);
122    *offset += vlc->table_size;
123}
124
125static av_cold void qdm2_init_vlc(void)
126{
127    const uint8_t (*tab)[2] = tab_fft_tone_offset;
128    unsigned offset = 0;
129
130    build_vlc(&vlc_tab_level, 8, 24, &offset, tab_level);
131    build_vlc(&vlc_tab_diff,  8, 33, &offset, tab_diff);
132    build_vlc(&vlc_tab_run,   5,  6, &offset, tab_run);
133
134    build_vlc(&fft_level_exp_alt_vlc, 8, 28, &offset, fft_level_exp_alt);
135    build_vlc(&fft_level_exp_vlc,     8, 20, &offset, fft_level_exp);
136
137    build_vlc(&fft_stereo_exp_vlc,   6, 7, &offset, fft_stereo_exp);
138    build_vlc(&fft_stereo_phase_vlc, 6, 9, &offset, fft_stereo_phase);
139
140    build_vlc(&vlc_tab_tone_level_idx_hi1, 8, 20, &offset, tab_tone_level_idx_hi1);
141    build_vlc(&vlc_tab_tone_level_idx_mid, 8, 13, &offset, tab_tone_level_idx_mid);
142    build_vlc(&vlc_tab_tone_level_idx_hi2, 8, 18, &offset, tab_tone_level_idx_hi2);
143
144    build_vlc(&vlc_tab_type30, 6,  9, &offset, tab_type30);
145    build_vlc(&vlc_tab_type34, 5, 10, &offset, tab_type34);
146
147    for (int i = 0; i < 5; i++) {
148        build_vlc(&vlc_tab_fft_tone_offset[i], 8, tab_fft_tone_offset_sizes[i],
149                  &offset, tab);
150        tab += tab_fft_tone_offset_sizes[i];
151    }
152}
153
154#endif /* CONFIG_HARDCODED_TABLES */
155
156#endif /* AVCODEC_QDM2_TABLEGEN_H */
157