Suga koubou / Mbed 2 deprecated i2s_audio_madplayer

streo mp3 player see: http://mbed.org/users/okini3939/notebook/I2S_AUDIO

Dependencies:   FatFileSystemCpp I2SSlave TLV320 mbed

Fork of madplayer by Andreas Grün

layer3.cpp@0:7627c79db971, 2010-11-26 (annotated)

Committer:
Gruenfrosch
Date:
Fri Nov 26 12:18:30 2010 +0000
Revision:
0:7627c79db971
Child:
2:f28cf0afd021
First Version

Who changed what in which revision?

UserRevisionLine numberNew contents of line
Gruenfrosch 0:7627c79db971 1 /*
Gruenfrosch 0:7627c79db971 2 * libmad - MPEG audio decoder library
Gruenfrosch 0:7627c79db971 3 * Copyright (C) 2000-2004 Underbit Technologies, Inc.
Gruenfrosch 0:7627c79db971 4 *
Gruenfrosch 0:7627c79db971 5 * This program is free software; you can redistribute it and/or modify
Gruenfrosch 0:7627c79db971 6 * it under the terms of the GNU General Public License as published by
Gruenfrosch 0:7627c79db971 7 * the Free Software Foundation; either version 2 of the License, or
Gruenfrosch 0:7627c79db971 8 * (at your option) any later version.
Gruenfrosch 0:7627c79db971 9 *
Gruenfrosch 0:7627c79db971 10 * This program is distributed in the hope that it will be useful,
Gruenfrosch 0:7627c79db971 11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
Gruenfrosch 0:7627c79db971 12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
Gruenfrosch 0:7627c79db971 13 * GNU General Public License for more details.
Gruenfrosch 0:7627c79db971 14 *
Gruenfrosch 0:7627c79db971 15 * You should have received a copy of the GNU General Public License
Gruenfrosch 0:7627c79db971 16 * along with this program; if not, write to the Free Software
Gruenfrosch 0:7627c79db971 17 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
Gruenfrosch 0:7627c79db971 18 *
Gruenfrosch 0:7627c79db971 19 * $Id: layer3.c,v 1.1 2010/11/23 20:12:57 andy Exp $
Gruenfrosch 0:7627c79db971 20 */
Gruenfrosch 0:7627c79db971 21
Gruenfrosch 0:7627c79db971 22 # include "config.h"
Gruenfrosch 0:7627c79db971 23
Gruenfrosch 0:7627c79db971 24 # include "global.h"
Gruenfrosch 0:7627c79db971 25
Gruenfrosch 0:7627c79db971 26 # include <stdlib.h>
Gruenfrosch 0:7627c79db971 27 # include <string.h>
Gruenfrosch 0:7627c79db971 28
Gruenfrosch 0:7627c79db971 29 # ifdef HAVE_ASSERT_H
Gruenfrosch 0:7627c79db971 30 # include <assert.h>
Gruenfrosch 0:7627c79db971 31 # endif
Gruenfrosch 0:7627c79db971 32
Gruenfrosch 0:7627c79db971 33 # ifdef HAVE_LIMITS_H
Gruenfrosch 0:7627c79db971 34 # include <limits.h>
Gruenfrosch 0:7627c79db971 35 # else
Gruenfrosch 0:7627c79db971 36 # define CHAR_BIT 8
Gruenfrosch 0:7627c79db971 37 # endif
Gruenfrosch 0:7627c79db971 38
Gruenfrosch 0:7627c79db971 39 # include "fixed.h"
Gruenfrosch 0:7627c79db971 40 # include "bit.h"
Gruenfrosch 0:7627c79db971 41 # include "stream.h"
Gruenfrosch 0:7627c79db971 42 # include "frame.h"
Gruenfrosch 0:7627c79db971 43 # include "huffman.h"
Gruenfrosch 0:7627c79db971 44 # include "layer3.h"
Gruenfrosch 0:7627c79db971 45
Gruenfrosch 0:7627c79db971 46 /* --- Layer III ----------------------------------------------------------- */
Gruenfrosch 0:7627c79db971 47
Gruenfrosch 0:7627c79db971 48 enum {
Gruenfrosch 0:7627c79db971 49 count1table_select = 0x01,
Gruenfrosch 0:7627c79db971 50 scalefac_scale = 0x02,
Gruenfrosch 0:7627c79db971 51 preflag = 0x04,
Gruenfrosch 0:7627c79db971 52 mixed_block_flag = 0x08
Gruenfrosch 0:7627c79db971 53 };
Gruenfrosch 0:7627c79db971 54
Gruenfrosch 0:7627c79db971 55 enum {
Gruenfrosch 0:7627c79db971 56 I_STEREO = 0x1,
Gruenfrosch 0:7627c79db971 57 MS_STEREO = 0x2
Gruenfrosch 0:7627c79db971 58 };
Gruenfrosch 0:7627c79db971 59
Gruenfrosch 0:7627c79db971 60
Gruenfrosch 0:7627c79db971 61 struct channel {
Gruenfrosch 0:7627c79db971 62 /* from side info */
Gruenfrosch 0:7627c79db971 63 unsigned short part2_3_length;
Gruenfrosch 0:7627c79db971 64 unsigned short big_values;
Gruenfrosch 0:7627c79db971 65 unsigned short global_gain;
Gruenfrosch 0:7627c79db971 66 unsigned short scalefac_compress;
Gruenfrosch 0:7627c79db971 67
Gruenfrosch 0:7627c79db971 68 unsigned char flags;
Gruenfrosch 0:7627c79db971 69 unsigned char block_type;
Gruenfrosch 0:7627c79db971 70 unsigned char table_select[3];
Gruenfrosch 0:7627c79db971 71 unsigned char subblock_gain[3];
Gruenfrosch 0:7627c79db971 72 unsigned char region0_count;
Gruenfrosch 0:7627c79db971 73 unsigned char region1_count;
Gruenfrosch 0:7627c79db971 74
Gruenfrosch 0:7627c79db971 75 /* from main_data */
Gruenfrosch 0:7627c79db971 76 unsigned char scalefac[39]; /* scalefac_l and/or scalefac_s */
Gruenfrosch 0:7627c79db971 77 };
Gruenfrosch 0:7627c79db971 78 struct granule {
Gruenfrosch 0:7627c79db971 79 struct channel ch[2];
Gruenfrosch 0:7627c79db971 80 };
Gruenfrosch 0:7627c79db971 81
Gruenfrosch 0:7627c79db971 82
Gruenfrosch 0:7627c79db971 83 struct sideinfo {
Gruenfrosch 0:7627c79db971 84 unsigned int main_data_begin;
Gruenfrosch 0:7627c79db971 85 unsigned int private_bits;
Gruenfrosch 0:7627c79db971 86
Gruenfrosch 0:7627c79db971 87 unsigned char scfsi[2];
Gruenfrosch 0:7627c79db971 88
Gruenfrosch 0:7627c79db971 89 struct granule gr[2];
Gruenfrosch 0:7627c79db971 90 };
Gruenfrosch 0:7627c79db971 91
Gruenfrosch 0:7627c79db971 92 /*
Gruenfrosch 0:7627c79db971 93 * scalefactor bit lengths
Gruenfrosch 0:7627c79db971 94 * derived from section 2.4.2.7 of ISO/IEC 11172-3
Gruenfrosch 0:7627c79db971 95 */
Gruenfrosch 0:7627c79db971 96 static
Gruenfrosch 0:7627c79db971 97 struct {
Gruenfrosch 0:7627c79db971 98 unsigned char slen1;
Gruenfrosch 0:7627c79db971 99 unsigned char slen2;
Gruenfrosch 0:7627c79db971 100 } const sflen_table[16] = {
Gruenfrosch 0:7627c79db971 101 { 0, 0 }, { 0, 1 }, { 0, 2 }, { 0, 3 },
Gruenfrosch 0:7627c79db971 102 { 3, 0 }, { 1, 1 }, { 1, 2 }, { 1, 3 },
Gruenfrosch 0:7627c79db971 103 { 2, 1 }, { 2, 2 }, { 2, 3 }, { 3, 1 },
Gruenfrosch 0:7627c79db971 104 { 3, 2 }, { 3, 3 }, { 4, 2 }, { 4, 3 }
Gruenfrosch 0:7627c79db971 105 };
Gruenfrosch 0:7627c79db971 106
Gruenfrosch 0:7627c79db971 107 /*
Gruenfrosch 0:7627c79db971 108 * number of LSF scalefactor band values
Gruenfrosch 0:7627c79db971 109 * derived from section 2.4.3.2 of ISO/IEC 13818-3
Gruenfrosch 0:7627c79db971 110 */
Gruenfrosch 0:7627c79db971 111 static
Gruenfrosch 0:7627c79db971 112 unsigned char const nsfb_table[6][3][4] = {
Gruenfrosch 0:7627c79db971 113 { { 6, 5, 5, 5 },
Gruenfrosch 0:7627c79db971 114 { 9, 9, 9, 9 },
Gruenfrosch 0:7627c79db971 115 { 6, 9, 9, 9 } },
Gruenfrosch 0:7627c79db971 116
Gruenfrosch 0:7627c79db971 117 { { 6, 5, 7, 3 },
Gruenfrosch 0:7627c79db971 118 { 9, 9, 12, 6 },
Gruenfrosch 0:7627c79db971 119 { 6, 9, 12, 6 } },
Gruenfrosch 0:7627c79db971 120
Gruenfrosch 0:7627c79db971 121 { { 11, 10, 0, 0 },
Gruenfrosch 0:7627c79db971 122 { 18, 18, 0, 0 },
Gruenfrosch 0:7627c79db971 123 { 15, 18, 0, 0 } },
Gruenfrosch 0:7627c79db971 124
Gruenfrosch 0:7627c79db971 125 { { 7, 7, 7, 0 },
Gruenfrosch 0:7627c79db971 126 { 12, 12, 12, 0 },
Gruenfrosch 0:7627c79db971 127 { 6, 15, 12, 0 } },
Gruenfrosch 0:7627c79db971 128
Gruenfrosch 0:7627c79db971 129 { { 6, 6, 6, 3 },
Gruenfrosch 0:7627c79db971 130 { 12, 9, 9, 6 },
Gruenfrosch 0:7627c79db971 131 { 6, 12, 9, 6 } },
Gruenfrosch 0:7627c79db971 132
Gruenfrosch 0:7627c79db971 133 { { 8, 8, 5, 0 },
Gruenfrosch 0:7627c79db971 134 { 15, 12, 9, 0 },
Gruenfrosch 0:7627c79db971 135 { 6, 18, 9, 0 } }
Gruenfrosch 0:7627c79db971 136 };
Gruenfrosch 0:7627c79db971 137
Gruenfrosch 0:7627c79db971 138 /*
Gruenfrosch 0:7627c79db971 139 * MPEG-1 scalefactor band widths
Gruenfrosch 0:7627c79db971 140 * derived from Table B.8 of ISO/IEC 11172-3
Gruenfrosch 0:7627c79db971 141 */
Gruenfrosch 0:7627c79db971 142 static
Gruenfrosch 0:7627c79db971 143 unsigned char const sfb_48000_long[] = {
Gruenfrosch 0:7627c79db971 144 4, 4, 4, 4, 4, 4, 6, 6, 6, 8, 10,
Gruenfrosch 0:7627c79db971 145 12, 16, 18, 22, 28, 34, 40, 46, 54, 54, 192
Gruenfrosch 0:7627c79db971 146 };
Gruenfrosch 0:7627c79db971 147
Gruenfrosch 0:7627c79db971 148 static
Gruenfrosch 0:7627c79db971 149 unsigned char const sfb_44100_long[] = {
Gruenfrosch 0:7627c79db971 150 4, 4, 4, 4, 4, 4, 6, 6, 8, 8, 10,
Gruenfrosch 0:7627c79db971 151 12, 16, 20, 24, 28, 34, 42, 50, 54, 76, 158
Gruenfrosch 0:7627c79db971 152 };
Gruenfrosch 0:7627c79db971 153
Gruenfrosch 0:7627c79db971 154 static
Gruenfrosch 0:7627c79db971 155 unsigned char const sfb_32000_long[] = {
Gruenfrosch 0:7627c79db971 156 4, 4, 4, 4, 4, 4, 6, 6, 8, 10, 12,
Gruenfrosch 0:7627c79db971 157 16, 20, 24, 30, 38, 46, 56, 68, 84, 102, 26
Gruenfrosch 0:7627c79db971 158 };
Gruenfrosch 0:7627c79db971 159
Gruenfrosch 0:7627c79db971 160 static
Gruenfrosch 0:7627c79db971 161 unsigned char const sfb_48000_short[] = {
Gruenfrosch 0:7627c79db971 162 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 6,
Gruenfrosch 0:7627c79db971 163 6, 6, 6, 6, 6, 10, 10, 10, 12, 12, 12, 14, 14,
Gruenfrosch 0:7627c79db971 164 14, 16, 16, 16, 20, 20, 20, 26, 26, 26, 66, 66, 66
Gruenfrosch 0:7627c79db971 165 };
Gruenfrosch 0:7627c79db971 166
Gruenfrosch 0:7627c79db971 167 static
Gruenfrosch 0:7627c79db971 168 unsigned char const sfb_44100_short[] = {
Gruenfrosch 0:7627c79db971 169 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 6,
Gruenfrosch 0:7627c79db971 170 6, 6, 8, 8, 8, 10, 10, 10, 12, 12, 12, 14, 14,
Gruenfrosch 0:7627c79db971 171 14, 18, 18, 18, 22, 22, 22, 30, 30, 30, 56, 56, 56
Gruenfrosch 0:7627c79db971 172 };
Gruenfrosch 0:7627c79db971 173
Gruenfrosch 0:7627c79db971 174 static
Gruenfrosch 0:7627c79db971 175 unsigned char const sfb_32000_short[] = {
Gruenfrosch 0:7627c79db971 176 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 6,
Gruenfrosch 0:7627c79db971 177 6, 6, 8, 8, 8, 12, 12, 12, 16, 16, 16, 20, 20,
Gruenfrosch 0:7627c79db971 178 20, 26, 26, 26, 34, 34, 34, 42, 42, 42, 12, 12, 12
Gruenfrosch 0:7627c79db971 179 };
Gruenfrosch 0:7627c79db971 180
Gruenfrosch 0:7627c79db971 181 static
Gruenfrosch 0:7627c79db971 182 unsigned char const sfb_48000_mixed[] = {
Gruenfrosch 0:7627c79db971 183 /* long */ 4, 4, 4, 4, 4, 4, 6, 6,
Gruenfrosch 0:7627c79db971 184 /* short */ 4, 4, 4, 6, 6, 6, 6, 6, 6, 10,
Gruenfrosch 0:7627c79db971 185 10, 10, 12, 12, 12, 14, 14, 14, 16, 16,
Gruenfrosch 0:7627c79db971 186 16, 20, 20, 20, 26, 26, 26, 66, 66, 66
Gruenfrosch 0:7627c79db971 187 };
Gruenfrosch 0:7627c79db971 188
Gruenfrosch 0:7627c79db971 189 static
Gruenfrosch 0:7627c79db971 190 unsigned char const sfb_44100_mixed[] = {
Gruenfrosch 0:7627c79db971 191 /* long */ 4, 4, 4, 4, 4, 4, 6, 6,
Gruenfrosch 0:7627c79db971 192 /* short */ 4, 4, 4, 6, 6, 6, 8, 8, 8, 10,
Gruenfrosch 0:7627c79db971 193 10, 10, 12, 12, 12, 14, 14, 14, 18, 18,
Gruenfrosch 0:7627c79db971 194 18, 22, 22, 22, 30, 30, 30, 56, 56, 56
Gruenfrosch 0:7627c79db971 195 };
Gruenfrosch 0:7627c79db971 196
Gruenfrosch 0:7627c79db971 197 static
Gruenfrosch 0:7627c79db971 198 unsigned char const sfb_32000_mixed[] = {
Gruenfrosch 0:7627c79db971 199 /* long */ 4, 4, 4, 4, 4, 4, 6, 6,
Gruenfrosch 0:7627c79db971 200 /* short */ 4, 4, 4, 6, 6, 6, 8, 8, 8, 12,
Gruenfrosch 0:7627c79db971 201 12, 12, 16, 16, 16, 20, 20, 20, 26, 26,
Gruenfrosch 0:7627c79db971 202 26, 34, 34, 34, 42, 42, 42, 12, 12, 12
Gruenfrosch 0:7627c79db971 203 };
Gruenfrosch 0:7627c79db971 204
Gruenfrosch 0:7627c79db971 205 /*
Gruenfrosch 0:7627c79db971 206 * MPEG-2 scalefactor band widths
Gruenfrosch 0:7627c79db971 207 * derived from Table B.2 of ISO/IEC 13818-3
Gruenfrosch 0:7627c79db971 208 */
Gruenfrosch 0:7627c79db971 209 static
Gruenfrosch 0:7627c79db971 210 unsigned char const sfb_24000_long[] = {
Gruenfrosch 0:7627c79db971 211 6, 6, 6, 6, 6, 6, 8, 10, 12, 14, 16,
Gruenfrosch 0:7627c79db971 212 18, 22, 26, 32, 38, 46, 54, 62, 70, 76, 36
Gruenfrosch 0:7627c79db971 213 };
Gruenfrosch 0:7627c79db971 214
Gruenfrosch 0:7627c79db971 215 static
Gruenfrosch 0:7627c79db971 216 unsigned char const sfb_22050_long[] = {
Gruenfrosch 0:7627c79db971 217 6, 6, 6, 6, 6, 6, 8, 10, 12, 14, 16,
Gruenfrosch 0:7627c79db971 218 20, 24, 28, 32, 38, 46, 52, 60, 68, 58, 54
Gruenfrosch 0:7627c79db971 219 };
Gruenfrosch 0:7627c79db971 220
Gruenfrosch 0:7627c79db971 221 # define sfb_16000_long sfb_22050_long
Gruenfrosch 0:7627c79db971 222
Gruenfrosch 0:7627c79db971 223 static
Gruenfrosch 0:7627c79db971 224 unsigned char const sfb_24000_short[] = {
Gruenfrosch 0:7627c79db971 225 4, 4, 4, 4, 4, 4, 4, 4, 4, 6, 6, 6, 8,
Gruenfrosch 0:7627c79db971 226 8, 8, 10, 10, 10, 12, 12, 12, 14, 14, 14, 18, 18,
Gruenfrosch 0:7627c79db971 227 18, 24, 24, 24, 32, 32, 32, 44, 44, 44, 12, 12, 12
Gruenfrosch 0:7627c79db971 228 };
Gruenfrosch 0:7627c79db971 229
Gruenfrosch 0:7627c79db971 230 static
Gruenfrosch 0:7627c79db971 231 unsigned char const sfb_22050_short[] = {
Gruenfrosch 0:7627c79db971 232 4, 4, 4, 4, 4, 4, 4, 4, 4, 6, 6, 6, 6,
Gruenfrosch 0:7627c79db971 233 6, 6, 8, 8, 8, 10, 10, 10, 14, 14, 14, 18, 18,
Gruenfrosch 0:7627c79db971 234 18, 26, 26, 26, 32, 32, 32, 42, 42, 42, 18, 18, 18
Gruenfrosch 0:7627c79db971 235 };
Gruenfrosch 0:7627c79db971 236
Gruenfrosch 0:7627c79db971 237 static
Gruenfrosch 0:7627c79db971 238 unsigned char const sfb_16000_short[] = {
Gruenfrosch 0:7627c79db971 239 4, 4, 4, 4, 4, 4, 4, 4, 4, 6, 6, 6, 8,
Gruenfrosch 0:7627c79db971 240 8, 8, 10, 10, 10, 12, 12, 12, 14, 14, 14, 18, 18,
Gruenfrosch 0:7627c79db971 241 18, 24, 24, 24, 30, 30, 30, 40, 40, 40, 18, 18, 18
Gruenfrosch 0:7627c79db971 242 };
Gruenfrosch 0:7627c79db971 243
Gruenfrosch 0:7627c79db971 244 static
Gruenfrosch 0:7627c79db971 245 unsigned char const sfb_24000_mixed[] = {
Gruenfrosch 0:7627c79db971 246 /* long */ 6, 6, 6, 6, 6, 6,
Gruenfrosch 0:7627c79db971 247 /* short */ 6, 6, 6, 8, 8, 8, 10, 10, 10, 12,
Gruenfrosch 0:7627c79db971 248 12, 12, 14, 14, 14, 18, 18, 18, 24, 24,
Gruenfrosch 0:7627c79db971 249 24, 32, 32, 32, 44, 44, 44, 12, 12, 12
Gruenfrosch 0:7627c79db971 250 };
Gruenfrosch 0:7627c79db971 251
Gruenfrosch 0:7627c79db971 252 static
Gruenfrosch 0:7627c79db971 253 unsigned char const sfb_22050_mixed[] = {
Gruenfrosch 0:7627c79db971 254 /* long */ 6, 6, 6, 6, 6, 6,
Gruenfrosch 0:7627c79db971 255 /* short */ 6, 6, 6, 6, 6, 6, 8, 8, 8, 10,
Gruenfrosch 0:7627c79db971 256 10, 10, 14, 14, 14, 18, 18, 18, 26, 26,
Gruenfrosch 0:7627c79db971 257 26, 32, 32, 32, 42, 42, 42, 18, 18, 18
Gruenfrosch 0:7627c79db971 258 };
Gruenfrosch 0:7627c79db971 259
Gruenfrosch 0:7627c79db971 260 static
Gruenfrosch 0:7627c79db971 261 unsigned char const sfb_16000_mixed[] = {
Gruenfrosch 0:7627c79db971 262 /* long */ 6, 6, 6, 6, 6, 6,
Gruenfrosch 0:7627c79db971 263 /* short */ 6, 6, 6, 8, 8, 8, 10, 10, 10, 12,
Gruenfrosch 0:7627c79db971 264 12, 12, 14, 14, 14, 18, 18, 18, 24, 24,
Gruenfrosch 0:7627c79db971 265 24, 30, 30, 30, 40, 40, 40, 18, 18, 18
Gruenfrosch 0:7627c79db971 266 };
Gruenfrosch 0:7627c79db971 267
Gruenfrosch 0:7627c79db971 268 /*
Gruenfrosch 0:7627c79db971 269 * MPEG 2.5 scalefactor band widths
Gruenfrosch 0:7627c79db971 270 * derived from public sources
Gruenfrosch 0:7627c79db971 271 */
Gruenfrosch 0:7627c79db971 272 # define sfb_12000_long sfb_16000_long
Gruenfrosch 0:7627c79db971 273 # define sfb_11025_long sfb_12000_long
Gruenfrosch 0:7627c79db971 274
Gruenfrosch 0:7627c79db971 275 static
Gruenfrosch 0:7627c79db971 276 unsigned char const sfb_8000_long[] = {
Gruenfrosch 0:7627c79db971 277 12, 12, 12, 12, 12, 12, 16, 20, 24, 28, 32,
Gruenfrosch 0:7627c79db971 278 40, 48, 56, 64, 76, 90, 2, 2, 2, 2, 2
Gruenfrosch 0:7627c79db971 279 };
Gruenfrosch 0:7627c79db971 280
Gruenfrosch 0:7627c79db971 281 # define sfb_12000_short sfb_16000_short
Gruenfrosch 0:7627c79db971 282 # define sfb_11025_short sfb_12000_short
Gruenfrosch 0:7627c79db971 283
Gruenfrosch 0:7627c79db971 284 static
Gruenfrosch 0:7627c79db971 285 unsigned char const sfb_8000_short[] = {
Gruenfrosch 0:7627c79db971 286 8, 8, 8, 8, 8, 8, 8, 8, 8, 12, 12, 12, 16,
Gruenfrosch 0:7627c79db971 287 16, 16, 20, 20, 20, 24, 24, 24, 28, 28, 28, 36, 36,
Gruenfrosch 0:7627c79db971 288 36, 2, 2, 2, 2, 2, 2, 2, 2, 2, 26, 26, 26
Gruenfrosch 0:7627c79db971 289 };
Gruenfrosch 0:7627c79db971 290
Gruenfrosch 0:7627c79db971 291 # define sfb_12000_mixed sfb_16000_mixed
Gruenfrosch 0:7627c79db971 292 # define sfb_11025_mixed sfb_12000_mixed
Gruenfrosch 0:7627c79db971 293
Gruenfrosch 0:7627c79db971 294 /* the 8000 Hz short block scalefactor bands do not break after
Gruenfrosch 0:7627c79db971 295 the first 36 frequency lines, so this is probably wrong */
Gruenfrosch 0:7627c79db971 296 static
Gruenfrosch 0:7627c79db971 297 unsigned char const sfb_8000_mixed[] = {
Gruenfrosch 0:7627c79db971 298 /* long */ 12, 12, 12,
Gruenfrosch 0:7627c79db971 299 /* short */ 4, 4, 4, 8, 8, 8, 12, 12, 12, 16, 16, 16,
Gruenfrosch 0:7627c79db971 300 20, 20, 20, 24, 24, 24, 28, 28, 28, 36, 36, 36,
Gruenfrosch 0:7627c79db971 301 2, 2, 2, 2, 2, 2, 2, 2, 2, 26, 26, 26
Gruenfrosch 0:7627c79db971 302 };
Gruenfrosch 0:7627c79db971 303
Gruenfrosch 0:7627c79db971 304 static
Gruenfrosch 0:7627c79db971 305 struct {
Gruenfrosch 0:7627c79db971 306 unsigned char const *l;
Gruenfrosch 0:7627c79db971 307 unsigned char const *s;
Gruenfrosch 0:7627c79db971 308 unsigned char const *m;
Gruenfrosch 0:7627c79db971 309 } const sfbwidth_table[9] = {
Gruenfrosch 0:7627c79db971 310 { sfb_48000_long, sfb_48000_short, sfb_48000_mixed },
Gruenfrosch 0:7627c79db971 311 { sfb_44100_long, sfb_44100_short, sfb_44100_mixed },
Gruenfrosch 0:7627c79db971 312 { sfb_32000_long, sfb_32000_short, sfb_32000_mixed },
Gruenfrosch 0:7627c79db971 313 { sfb_24000_long, sfb_24000_short, sfb_24000_mixed },
Gruenfrosch 0:7627c79db971 314 { sfb_22050_long, sfb_22050_short, sfb_22050_mixed },
Gruenfrosch 0:7627c79db971 315 { sfb_16000_long, sfb_16000_short, sfb_16000_mixed },
Gruenfrosch 0:7627c79db971 316 { sfb_12000_long, sfb_12000_short, sfb_12000_mixed },
Gruenfrosch 0:7627c79db971 317 { sfb_11025_long, sfb_11025_short, sfb_11025_mixed },
Gruenfrosch 0:7627c79db971 318 { sfb_8000_long, sfb_8000_short, sfb_8000_mixed }
Gruenfrosch 0:7627c79db971 319 };
Gruenfrosch 0:7627c79db971 320
Gruenfrosch 0:7627c79db971 321 /*
Gruenfrosch 0:7627c79db971 322 * scalefactor band preemphasis (used only when preflag is set)
Gruenfrosch 0:7627c79db971 323 * derived from Table B.6 of ISO/IEC 11172-3
Gruenfrosch 0:7627c79db971 324 */
Gruenfrosch 0:7627c79db971 325 static
Gruenfrosch 0:7627c79db971 326 unsigned char const pretab[22] = {
Gruenfrosch 0:7627c79db971 327 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 3, 3, 3, 2, 0
Gruenfrosch 0:7627c79db971 328 };
Gruenfrosch 0:7627c79db971 329
Gruenfrosch 0:7627c79db971 330 /*
Gruenfrosch 0:7627c79db971 331 * table for requantization
Gruenfrosch 0:7627c79db971 332 *
Gruenfrosch 0:7627c79db971 333 * rq_table[x].mantissa * 2^(rq_table[x].exponent) = x^(4/3)
Gruenfrosch 0:7627c79db971 334 */
Gruenfrosch 0:7627c79db971 335 static
Gruenfrosch 0:7627c79db971 336 struct fixedfloat {
Gruenfrosch 0:7627c79db971 337 unsigned long mantissa : 27;
Gruenfrosch 0:7627c79db971 338 unsigned short exponent : 5;
Gruenfrosch 0:7627c79db971 339 } const rq_table[8207] = {
Gruenfrosch 0:7627c79db971 340 # include "rq_table.h"
Gruenfrosch 0:7627c79db971 341 };
Gruenfrosch 0:7627c79db971 342
Gruenfrosch 0:7627c79db971 343 /*
Gruenfrosch 0:7627c79db971 344 * fractional powers of two
Gruenfrosch 0:7627c79db971 345 * used for requantization and joint stereo decoding
Gruenfrosch 0:7627c79db971 346 *
Gruenfrosch 0:7627c79db971 347 * root_table[3 + x] = 2^(x/4)
Gruenfrosch 0:7627c79db971 348 */
Gruenfrosch 0:7627c79db971 349 static
Gruenfrosch 0:7627c79db971 350 mad_fixed_t const root_table[7] = {
Gruenfrosch 0:7627c79db971 351 MAD_F(0x09837f05) /* 2^(-3/4) == 0.59460355750136 */,
Gruenfrosch 0:7627c79db971 352 MAD_F(0x0b504f33) /* 2^(-2/4) == 0.70710678118655 */,
Gruenfrosch 0:7627c79db971 353 MAD_F(0x0d744fcd) /* 2^(-1/4) == 0.84089641525371 */,
Gruenfrosch 0:7627c79db971 354 MAD_F(0x10000000) /* 2^( 0/4) == 1.00000000000000 */,
Gruenfrosch 0:7627c79db971 355 MAD_F(0x1306fe0a) /* 2^(+1/4) == 1.18920711500272 */,
Gruenfrosch 0:7627c79db971 356 MAD_F(0x16a09e66) /* 2^(+2/4) == 1.41421356237310 */,
Gruenfrosch 0:7627c79db971 357 MAD_F(0x1ae89f99) /* 2^(+3/4) == 1.68179283050743 */
Gruenfrosch 0:7627c79db971 358 };
Gruenfrosch 0:7627c79db971 359
Gruenfrosch 0:7627c79db971 360 /*
Gruenfrosch 0:7627c79db971 361 * coefficients for aliasing reduction
Gruenfrosch 0:7627c79db971 362 * derived from Table B.9 of ISO/IEC 11172-3
Gruenfrosch 0:7627c79db971 363 *
Gruenfrosch 0:7627c79db971 364 * c[] = { -0.6, -0.535, -0.33, -0.185, -0.095, -0.041, -0.0142, -0.0037 }
Gruenfrosch 0:7627c79db971 365 * cs[i] = 1 / sqrt(1 + c[i]^2)
Gruenfrosch 0:7627c79db971 366 * ca[i] = c[i] / sqrt(1 + c[i]^2)
Gruenfrosch 0:7627c79db971 367 */
Gruenfrosch 0:7627c79db971 368 static
Gruenfrosch 0:7627c79db971 369 mad_fixed_t const cs[8] = {
Gruenfrosch 0:7627c79db971 370 +MAD_F(0x0db84a81) /* +0.857492926 */, +MAD_F(0x0e1b9d7f) /* +0.881741997 */,
Gruenfrosch 0:7627c79db971 371 +MAD_F(0x0f31adcf) /* +0.949628649 */, +MAD_F(0x0fbba815) /* +0.983314592 */,
Gruenfrosch 0:7627c79db971 372 +MAD_F(0x0feda417) /* +0.995517816 */, +MAD_F(0x0ffc8fc8) /* +0.999160558 */,
Gruenfrosch 0:7627c79db971 373 +MAD_F(0x0fff964c) /* +0.999899195 */, +MAD_F(0x0ffff8d3) /* +0.999993155 */
Gruenfrosch 0:7627c79db971 374 };
Gruenfrosch 0:7627c79db971 375
Gruenfrosch 0:7627c79db971 376 static
Gruenfrosch 0:7627c79db971 377 mad_fixed_t const ca[8] = {
Gruenfrosch 0:7627c79db971 378 -MAD_F(0x083b5fe7) /* -0.514495755 */, -MAD_F(0x078c36d2) /* -0.471731969 */,
Gruenfrosch 0:7627c79db971 379 -MAD_F(0x05039814) /* -0.313377454 */, -MAD_F(0x02e91dd1) /* -0.181913200 */,
Gruenfrosch 0:7627c79db971 380 -MAD_F(0x0183603a) /* -0.094574193 */, -MAD_F(0x00a7cb87) /* -0.040965583 */,
Gruenfrosch 0:7627c79db971 381 -MAD_F(0x003a2847) /* -0.014198569 */, -MAD_F(0x000f27b4) /* -0.003699975 */
Gruenfrosch 0:7627c79db971 382 };
Gruenfrosch 0:7627c79db971 383
Gruenfrosch 0:7627c79db971 384 /*
Gruenfrosch 0:7627c79db971 385 * IMDCT coefficients for short blocks
Gruenfrosch 0:7627c79db971 386 * derived from section 2.4.3.4.10.2 of ISO/IEC 11172-3
Gruenfrosch 0:7627c79db971 387 *
Gruenfrosch 0:7627c79db971 388 * imdct_s[i/even][k] = cos((PI / 24) * (2 * (i / 2) + 7) * (2 * k + 1))
Gruenfrosch 0:7627c79db971 389 * imdct_s[i /odd][k] = cos((PI / 24) * (2 * (6 + (i-1)/2) + 7) * (2 * k + 1))
Gruenfrosch 0:7627c79db971 390 */
Gruenfrosch 0:7627c79db971 391 static
Gruenfrosch 0:7627c79db971 392 mad_fixed_t const imdct_s[6][6] = {
Gruenfrosch 0:7627c79db971 393 # include "imdct_s.h"
Gruenfrosch 0:7627c79db971 394 };
Gruenfrosch 0:7627c79db971 395
Gruenfrosch 0:7627c79db971 396 # if !defined(ASO_IMDCT)
Gruenfrosch 0:7627c79db971 397 /*
Gruenfrosch 0:7627c79db971 398 * windowing coefficients for long blocks
Gruenfrosch 0:7627c79db971 399 * derived from section 2.4.3.4.10.3 of ISO/IEC 11172-3
Gruenfrosch 0:7627c79db971 400 *
Gruenfrosch 0:7627c79db971 401 * window_l[i] = sin((PI / 36) * (i + 1/2))
Gruenfrosch 0:7627c79db971 402 */
Gruenfrosch 0:7627c79db971 403 static
Gruenfrosch 0:7627c79db971 404 mad_fixed_t const window_l[36] = {
Gruenfrosch 0:7627c79db971 405 MAD_F(0x00b2aa3e) /* 0.043619387 */, MAD_F(0x0216a2a2) /* 0.130526192 */,
Gruenfrosch 0:7627c79db971 406 MAD_F(0x03768962) /* 0.216439614 */, MAD_F(0x04cfb0e2) /* 0.300705800 */,
Gruenfrosch 0:7627c79db971 407 MAD_F(0x061f78aa) /* 0.382683432 */, MAD_F(0x07635284) /* 0.461748613 */,
Gruenfrosch 0:7627c79db971 408 MAD_F(0x0898c779) /* 0.537299608 */, MAD_F(0x09bd7ca0) /* 0.608761429 */,
Gruenfrosch 0:7627c79db971 409 MAD_F(0x0acf37ad) /* 0.675590208 */, MAD_F(0x0bcbe352) /* 0.737277337 */,
Gruenfrosch 0:7627c79db971 410 MAD_F(0x0cb19346) /* 0.793353340 */, MAD_F(0x0d7e8807) /* 0.843391446 */,
Gruenfrosch 0:7627c79db971 411
Gruenfrosch 0:7627c79db971 412 MAD_F(0x0e313245) /* 0.887010833 */, MAD_F(0x0ec835e8) /* 0.923879533 */,
Gruenfrosch 0:7627c79db971 413 MAD_F(0x0f426cb5) /* 0.953716951 */, MAD_F(0x0f9ee890) /* 0.976296007 */,
Gruenfrosch 0:7627c79db971 414 MAD_F(0x0fdcf549) /* 0.991444861 */, MAD_F(0x0ffc19fd) /* 0.999048222 */,
Gruenfrosch 0:7627c79db971 415 MAD_F(0x0ffc19fd) /* 0.999048222 */, MAD_F(0x0fdcf549) /* 0.991444861 */,
Gruenfrosch 0:7627c79db971 416 MAD_F(0x0f9ee890) /* 0.976296007 */, MAD_F(0x0f426cb5) /* 0.953716951 */,
Gruenfrosch 0:7627c79db971 417 MAD_F(0x0ec835e8) /* 0.923879533 */, MAD_F(0x0e313245) /* 0.887010833 */,
Gruenfrosch 0:7627c79db971 418
Gruenfrosch 0:7627c79db971 419 MAD_F(0x0d7e8807) /* 0.843391446 */, MAD_F(0x0cb19346) /* 0.793353340 */,
Gruenfrosch 0:7627c79db971 420 MAD_F(0x0bcbe352) /* 0.737277337 */, MAD_F(0x0acf37ad) /* 0.675590208 */,
Gruenfrosch 0:7627c79db971 421 MAD_F(0x09bd7ca0) /* 0.608761429 */, MAD_F(0x0898c779) /* 0.537299608 */,
Gruenfrosch 0:7627c79db971 422 MAD_F(0x07635284) /* 0.461748613 */, MAD_F(0x061f78aa) /* 0.382683432 */,
Gruenfrosch 0:7627c79db971 423 MAD_F(0x04cfb0e2) /* 0.300705800 */, MAD_F(0x03768962) /* 0.216439614 */,
Gruenfrosch 0:7627c79db971 424 MAD_F(0x0216a2a2) /* 0.130526192 */, MAD_F(0x00b2aa3e) /* 0.043619387 */,
Gruenfrosch 0:7627c79db971 425 };
Gruenfrosch 0:7627c79db971 426 # endif /* ASO_IMDCT */
Gruenfrosch 0:7627c79db971 427
Gruenfrosch 0:7627c79db971 428 /*
Gruenfrosch 0:7627c79db971 429 * windowing coefficients for short blocks
Gruenfrosch 0:7627c79db971 430 * derived from section 2.4.3.4.10.3 of ISO/IEC 11172-3
Gruenfrosch 0:7627c79db971 431 *
Gruenfrosch 0:7627c79db971 432 * window_s[i] = sin((PI / 12) * (i + 1/2))
Gruenfrosch 0:7627c79db971 433 */
Gruenfrosch 0:7627c79db971 434 static
Gruenfrosch 0:7627c79db971 435 mad_fixed_t const window_s[12] = {
Gruenfrosch 0:7627c79db971 436 MAD_F(0x0216a2a2) /* 0.130526192 */, MAD_F(0x061f78aa) /* 0.382683432 */,
Gruenfrosch 0:7627c79db971 437 MAD_F(0x09bd7ca0) /* 0.608761429 */, MAD_F(0x0cb19346) /* 0.793353340 */,
Gruenfrosch 0:7627c79db971 438 MAD_F(0x0ec835e8) /* 0.923879533 */, MAD_F(0x0fdcf549) /* 0.991444861 */,
Gruenfrosch 0:7627c79db971 439 MAD_F(0x0fdcf549) /* 0.991444861 */, MAD_F(0x0ec835e8) /* 0.923879533 */,
Gruenfrosch 0:7627c79db971 440 MAD_F(0x0cb19346) /* 0.793353340 */, MAD_F(0x09bd7ca0) /* 0.608761429 */,
Gruenfrosch 0:7627c79db971 441 MAD_F(0x061f78aa) /* 0.382683432 */, MAD_F(0x0216a2a2) /* 0.130526192 */,
Gruenfrosch 0:7627c79db971 442 };
Gruenfrosch 0:7627c79db971 443
Gruenfrosch 0:7627c79db971 444 /*
Gruenfrosch 0:7627c79db971 445 * coefficients for intensity stereo processing
Gruenfrosch 0:7627c79db971 446 * derived from section 2.4.3.4.9.3 of ISO/IEC 11172-3
Gruenfrosch 0:7627c79db971 447 *
Gruenfrosch 0:7627c79db971 448 * is_ratio[i] = tan(i * (PI / 12))
Gruenfrosch 0:7627c79db971 449 * is_table[i] = is_ratio[i] / (1 + is_ratio[i])
Gruenfrosch 0:7627c79db971 450 */
Gruenfrosch 0:7627c79db971 451 static
Gruenfrosch 0:7627c79db971 452 mad_fixed_t const is_table[7] = {
Gruenfrosch 0:7627c79db971 453 MAD_F(0x00000000) /* 0.000000000 */,
Gruenfrosch 0:7627c79db971 454 MAD_F(0x0361962f) /* 0.211324865 */,
Gruenfrosch 0:7627c79db971 455 MAD_F(0x05db3d74) /* 0.366025404 */,
Gruenfrosch 0:7627c79db971 456 MAD_F(0x08000000) /* 0.500000000 */,
Gruenfrosch 0:7627c79db971 457 MAD_F(0x0a24c28c) /* 0.633974596 */,
Gruenfrosch 0:7627c79db971 458 MAD_F(0x0c9e69d1) /* 0.788675135 */,
Gruenfrosch 0:7627c79db971 459 MAD_F(0x10000000) /* 1.000000000 */
Gruenfrosch 0:7627c79db971 460 };
Gruenfrosch 0:7627c79db971 461
Gruenfrosch 0:7627c79db971 462 /*
Gruenfrosch 0:7627c79db971 463 * coefficients for LSF intensity stereo processing
Gruenfrosch 0:7627c79db971 464 * derived from section 2.4.3.2 of ISO/IEC 13818-3
Gruenfrosch 0:7627c79db971 465 *
Gruenfrosch 0:7627c79db971 466 * is_lsf_table[0][i] = (1 / sqrt(sqrt(2)))^(i + 1)
Gruenfrosch 0:7627c79db971 467 * is_lsf_table[1][i] = (1 / sqrt(2)) ^(i + 1)
Gruenfrosch 0:7627c79db971 468 */
Gruenfrosch 0:7627c79db971 469 static
Gruenfrosch 0:7627c79db971 470 mad_fixed_t const is_lsf_table[2][15] = {
Gruenfrosch 0:7627c79db971 471 {
Gruenfrosch 0:7627c79db971 472 MAD_F(0x0d744fcd) /* 0.840896415 */,
Gruenfrosch 0:7627c79db971 473 MAD_F(0x0b504f33) /* 0.707106781 */,
Gruenfrosch 0:7627c79db971 474 MAD_F(0x09837f05) /* 0.594603558 */,
Gruenfrosch 0:7627c79db971 475 MAD_F(0x08000000) /* 0.500000000 */,
Gruenfrosch 0:7627c79db971 476 MAD_F(0x06ba27e6) /* 0.420448208 */,
Gruenfrosch 0:7627c79db971 477 MAD_F(0x05a8279a) /* 0.353553391 */,
Gruenfrosch 0:7627c79db971 478 MAD_F(0x04c1bf83) /* 0.297301779 */,
Gruenfrosch 0:7627c79db971 479 MAD_F(0x04000000) /* 0.250000000 */,
Gruenfrosch 0:7627c79db971 480 MAD_F(0x035d13f3) /* 0.210224104 */,
Gruenfrosch 0:7627c79db971 481 MAD_F(0x02d413cd) /* 0.176776695 */,
Gruenfrosch 0:7627c79db971 482 MAD_F(0x0260dfc1) /* 0.148650889 */,
Gruenfrosch 0:7627c79db971 483 MAD_F(0x02000000) /* 0.125000000 */,
Gruenfrosch 0:7627c79db971 484 MAD_F(0x01ae89fa) /* 0.105112052 */,
Gruenfrosch 0:7627c79db971 485 MAD_F(0x016a09e6) /* 0.088388348 */,
Gruenfrosch 0:7627c79db971 486 MAD_F(0x01306fe1) /* 0.074325445 */
Gruenfrosch 0:7627c79db971 487 }, {
Gruenfrosch 0:7627c79db971 488 MAD_F(0x0b504f33) /* 0.707106781 */,
Gruenfrosch 0:7627c79db971 489 MAD_F(0x08000000) /* 0.500000000 */,
Gruenfrosch 0:7627c79db971 490 MAD_F(0x05a8279a) /* 0.353553391 */,
Gruenfrosch 0:7627c79db971 491 MAD_F(0x04000000) /* 0.250000000 */,
Gruenfrosch 0:7627c79db971 492 MAD_F(0x02d413cd) /* 0.176776695 */,
Gruenfrosch 0:7627c79db971 493 MAD_F(0x02000000) /* 0.125000000 */,
Gruenfrosch 0:7627c79db971 494 MAD_F(0x016a09e6) /* 0.088388348 */,
Gruenfrosch 0:7627c79db971 495 MAD_F(0x01000000) /* 0.062500000 */,
Gruenfrosch 0:7627c79db971 496 MAD_F(0x00b504f3) /* 0.044194174 */,
Gruenfrosch 0:7627c79db971 497 MAD_F(0x00800000) /* 0.031250000 */,
Gruenfrosch 0:7627c79db971 498 MAD_F(0x005a827a) /* 0.022097087 */,
Gruenfrosch 0:7627c79db971 499 MAD_F(0x00400000) /* 0.015625000 */,
Gruenfrosch 0:7627c79db971 500 MAD_F(0x002d413d) /* 0.011048543 */,
Gruenfrosch 0:7627c79db971 501 MAD_F(0x00200000) /* 0.007812500 */,
Gruenfrosch 0:7627c79db971 502 MAD_F(0x0016a09e) /* 0.005524272 */
Gruenfrosch 0:7627c79db971 503 }
Gruenfrosch 0:7627c79db971 504 };
Gruenfrosch 0:7627c79db971 505
Gruenfrosch 0:7627c79db971 506 /*
Gruenfrosch 0:7627c79db971 507 * NAME: III_sideinfo()
Gruenfrosch 0:7627c79db971 508 * DESCRIPTION: decode frame side information from a bitstream
Gruenfrosch 0:7627c79db971 509 */
Gruenfrosch 0:7627c79db971 510 static
Gruenfrosch 0:7627c79db971 511 enum mad_error III_sideinfo(struct mad_bitptr *ptr, unsigned int nch,
Gruenfrosch 0:7627c79db971 512 int lsf, struct sideinfo *si,
Gruenfrosch 0:7627c79db971 513 unsigned int *data_bitlen,
Gruenfrosch 0:7627c79db971 514 unsigned int *priv_bitlen)
Gruenfrosch 0:7627c79db971 515 {
Gruenfrosch 0:7627c79db971 516 unsigned int ngr, gr, ch, i;
Gruenfrosch 0:7627c79db971 517 enum mad_error result = MAD_ERROR_NONE;
Gruenfrosch 0:7627c79db971 518
Gruenfrosch 0:7627c79db971 519 *data_bitlen = 0;
Gruenfrosch 0:7627c79db971 520 *priv_bitlen = lsf ? ((nch == 1) ? 1 : 2) : ((nch == 1) ? 5 : 3);
Gruenfrosch 0:7627c79db971 521
Gruenfrosch 0:7627c79db971 522 si->main_data_begin = mad_bit_read(ptr, lsf ? 8 : 9);
Gruenfrosch 0:7627c79db971 523 si->private_bits = mad_bit_read(ptr, *priv_bitlen);
Gruenfrosch 0:7627c79db971 524
Gruenfrosch 0:7627c79db971 525 ngr = 1;
Gruenfrosch 0:7627c79db971 526 if (!lsf) {
Gruenfrosch 0:7627c79db971 527 ngr = 2;
Gruenfrosch 0:7627c79db971 528
Gruenfrosch 0:7627c79db971 529 for (ch = 0; ch < nch; ++ch)
Gruenfrosch 0:7627c79db971 530 si->scfsi[ch] = mad_bit_read(ptr, 4);
Gruenfrosch 0:7627c79db971 531 }
Gruenfrosch 0:7627c79db971 532
Gruenfrosch 0:7627c79db971 533 for (gr = 0; gr < ngr; ++gr) {
Gruenfrosch 0:7627c79db971 534 struct granule *granule = &si->gr[gr];
Gruenfrosch 0:7627c79db971 535
Gruenfrosch 0:7627c79db971 536 for (ch = 0; ch < nch; ++ch) {
Gruenfrosch 0:7627c79db971 537 struct channel *channel = &granule->ch[ch];
Gruenfrosch 0:7627c79db971 538
Gruenfrosch 0:7627c79db971 539 channel->part2_3_length = mad_bit_read(ptr, 12);
Gruenfrosch 0:7627c79db971 540 channel->big_values = mad_bit_read(ptr, 9);
Gruenfrosch 0:7627c79db971 541 channel->global_gain = mad_bit_read(ptr, 8);
Gruenfrosch 0:7627c79db971 542 channel->scalefac_compress = mad_bit_read(ptr, lsf ? 9 : 4);
Gruenfrosch 0:7627c79db971 543
Gruenfrosch 0:7627c79db971 544 *data_bitlen += channel->part2_3_length;
Gruenfrosch 0:7627c79db971 545
Gruenfrosch 0:7627c79db971 546 if (channel->big_values > 288 && result == 0)
Gruenfrosch 0:7627c79db971 547 result = MAD_ERROR_BADBIGVALUES;
Gruenfrosch 0:7627c79db971 548
Gruenfrosch 0:7627c79db971 549 channel->flags = 0;
Gruenfrosch 0:7627c79db971 550
Gruenfrosch 0:7627c79db971 551 /* window_switching_flag */
Gruenfrosch 0:7627c79db971 552 if (mad_bit_read(ptr, 1)) {
Gruenfrosch 0:7627c79db971 553 channel->block_type = mad_bit_read(ptr, 2);
Gruenfrosch 0:7627c79db971 554
Gruenfrosch 0:7627c79db971 555 if (channel->block_type == 0 && result == 0)
Gruenfrosch 0:7627c79db971 556 result = MAD_ERROR_BADBLOCKTYPE;
Gruenfrosch 0:7627c79db971 557
Gruenfrosch 0:7627c79db971 558 if (!lsf && channel->block_type == 2 && si->scfsi[ch] && result == 0)
Gruenfrosch 0:7627c79db971 559 result = MAD_ERROR_BADSCFSI;
Gruenfrosch 0:7627c79db971 560
Gruenfrosch 0:7627c79db971 561 channel->region0_count = 7;
Gruenfrosch 0:7627c79db971 562 channel->region1_count = 36;
Gruenfrosch 0:7627c79db971 563
Gruenfrosch 0:7627c79db971 564 if (mad_bit_read(ptr, 1))
Gruenfrosch 0:7627c79db971 565 channel->flags |= mixed_block_flag;
Gruenfrosch 0:7627c79db971 566 else if (channel->block_type == 2)
Gruenfrosch 0:7627c79db971 567 channel->region0_count = 8;
Gruenfrosch 0:7627c79db971 568
Gruenfrosch 0:7627c79db971 569 for (i = 0; i < 2; ++i)
Gruenfrosch 0:7627c79db971 570 channel->table_select[i] = mad_bit_read(ptr, 5);
Gruenfrosch 0:7627c79db971 571
Gruenfrosch 0:7627c79db971 572 # if defined(DEBUG)
Gruenfrosch 0:7627c79db971 573 channel->table_select[2] = 4; /* not used */
Gruenfrosch 0:7627c79db971 574 # endif
Gruenfrosch 0:7627c79db971 575
Gruenfrosch 0:7627c79db971 576 for (i = 0; i < 3; ++i)
Gruenfrosch 0:7627c79db971 577 channel->subblock_gain[i] = mad_bit_read(ptr, 3);
Gruenfrosch 0:7627c79db971 578 }
Gruenfrosch 0:7627c79db971 579 else {
Gruenfrosch 0:7627c79db971 580 channel->block_type = 0;
Gruenfrosch 0:7627c79db971 581
Gruenfrosch 0:7627c79db971 582 for (i = 0; i < 3; ++i)
Gruenfrosch 0:7627c79db971 583 channel->table_select[i] = mad_bit_read(ptr, 5);
Gruenfrosch 0:7627c79db971 584
Gruenfrosch 0:7627c79db971 585 channel->region0_count = mad_bit_read(ptr, 4);
Gruenfrosch 0:7627c79db971 586 channel->region1_count = mad_bit_read(ptr, 3);
Gruenfrosch 0:7627c79db971 587 }
Gruenfrosch 0:7627c79db971 588
Gruenfrosch 0:7627c79db971 589 /* [preflag,] scalefac_scale, count1table_select */
Gruenfrosch 0:7627c79db971 590 channel->flags |= mad_bit_read(ptr, lsf ? 2 : 3);
Gruenfrosch 0:7627c79db971 591 }
Gruenfrosch 0:7627c79db971 592 }
Gruenfrosch 0:7627c79db971 593
Gruenfrosch 0:7627c79db971 594 return result;
Gruenfrosch 0:7627c79db971 595 }
Gruenfrosch 0:7627c79db971 596
Gruenfrosch 0:7627c79db971 597 /*
Gruenfrosch 0:7627c79db971 598 * NAME: III_scalefactors_lsf()
Gruenfrosch 0:7627c79db971 599 * DESCRIPTION: decode channel scalefactors for LSF from a bitstream
Gruenfrosch 0:7627c79db971 600 */
Gruenfrosch 0:7627c79db971 601 static
Gruenfrosch 0:7627c79db971 602 unsigned int III_scalefactors_lsf(struct mad_bitptr *ptr,
Gruenfrosch 0:7627c79db971 603 struct channel *channel,
Gruenfrosch 0:7627c79db971 604 struct channel *gr1ch, int mode_extension)
Gruenfrosch 0:7627c79db971 605 {
Gruenfrosch 0:7627c79db971 606 struct mad_bitptr start;
Gruenfrosch 0:7627c79db971 607 unsigned int scalefac_compress, index, slen[4], part, n, i;
Gruenfrosch 0:7627c79db971 608 unsigned char const *nsfb;
Gruenfrosch 0:7627c79db971 609
Gruenfrosch 0:7627c79db971 610 start = *ptr;
Gruenfrosch 0:7627c79db971 611
Gruenfrosch 0:7627c79db971 612 scalefac_compress = channel->scalefac_compress;
Gruenfrosch 0:7627c79db971 613 index = (channel->block_type == 2) ?
Gruenfrosch 0:7627c79db971 614 ((channel->flags & mixed_block_flag) ? 2 : 1) : 0;
Gruenfrosch 0:7627c79db971 615
Gruenfrosch 0:7627c79db971 616 if (!((mode_extension & I_STEREO) && gr1ch)) {
Gruenfrosch 0:7627c79db971 617 if (scalefac_compress < 400) {
Gruenfrosch 0:7627c79db971 618 slen[0] = (scalefac_compress >> 4) / 5;
Gruenfrosch 0:7627c79db971 619 slen[1] = (scalefac_compress >> 4) % 5;
Gruenfrosch 0:7627c79db971 620 slen[2] = (scalefac_compress % 16) >> 2;
Gruenfrosch 0:7627c79db971 621 slen[3] = scalefac_compress % 4;
Gruenfrosch 0:7627c79db971 622
Gruenfrosch 0:7627c79db971 623 nsfb = nsfb_table[0][index];
Gruenfrosch 0:7627c79db971 624 }
Gruenfrosch 0:7627c79db971 625 else if (scalefac_compress < 500) {
Gruenfrosch 0:7627c79db971 626 scalefac_compress -= 400;
Gruenfrosch 0:7627c79db971 627
Gruenfrosch 0:7627c79db971 628 slen[0] = (scalefac_compress >> 2) / 5;
Gruenfrosch 0:7627c79db971 629 slen[1] = (scalefac_compress >> 2) % 5;
Gruenfrosch 0:7627c79db971 630 slen[2] = scalefac_compress % 4;
Gruenfrosch 0:7627c79db971 631 slen[3] = 0;
Gruenfrosch 0:7627c79db971 632
Gruenfrosch 0:7627c79db971 633 nsfb = nsfb_table[1][index];
Gruenfrosch 0:7627c79db971 634 }
Gruenfrosch 0:7627c79db971 635 else {
Gruenfrosch 0:7627c79db971 636 scalefac_compress -= 500;
Gruenfrosch 0:7627c79db971 637
Gruenfrosch 0:7627c79db971 638 slen[0] = scalefac_compress / 3;
Gruenfrosch 0:7627c79db971 639 slen[1] = scalefac_compress % 3;
Gruenfrosch 0:7627c79db971 640 slen[2] = 0;
Gruenfrosch 0:7627c79db971 641 slen[3] = 0;
Gruenfrosch 0:7627c79db971 642
Gruenfrosch 0:7627c79db971 643 channel->flags |= preflag;
Gruenfrosch 0:7627c79db971 644
Gruenfrosch 0:7627c79db971 645 nsfb = nsfb_table[2][index];
Gruenfrosch 0:7627c79db971 646 }
Gruenfrosch 0:7627c79db971 647
Gruenfrosch 0:7627c79db971 648 n = 0;
Gruenfrosch 0:7627c79db971 649 for (part = 0; part < 4; ++part) {
Gruenfrosch 0:7627c79db971 650 for (i = 0; i < nsfb[part]; ++i)
Gruenfrosch 0:7627c79db971 651 channel->scalefac[n++] = mad_bit_read(ptr, slen[part]);
Gruenfrosch 0:7627c79db971 652 }
Gruenfrosch 0:7627c79db971 653
Gruenfrosch 0:7627c79db971 654 while (n < 39)
Gruenfrosch 0:7627c79db971 655 channel->scalefac[n++] = 0;
Gruenfrosch 0:7627c79db971 656 }
Gruenfrosch 0:7627c79db971 657 else { /* (mode_extension & I_STEREO) && gr1ch (i.e. ch == 1) */
Gruenfrosch 0:7627c79db971 658 scalefac_compress >>= 1;
Gruenfrosch 0:7627c79db971 659
Gruenfrosch 0:7627c79db971 660 if (scalefac_compress < 180) {
Gruenfrosch 0:7627c79db971 661 slen[0] = scalefac_compress / 36;
Gruenfrosch 0:7627c79db971 662 slen[1] = (scalefac_compress % 36) / 6;
Gruenfrosch 0:7627c79db971 663 slen[2] = (scalefac_compress % 36) % 6;
Gruenfrosch 0:7627c79db971 664 slen[3] = 0;
Gruenfrosch 0:7627c79db971 665
Gruenfrosch 0:7627c79db971 666 nsfb = nsfb_table[3][index];
Gruenfrosch 0:7627c79db971 667 }
Gruenfrosch 0:7627c79db971 668 else if (scalefac_compress < 244) {
Gruenfrosch 0:7627c79db971 669 scalefac_compress -= 180;
Gruenfrosch 0:7627c79db971 670
Gruenfrosch 0:7627c79db971 671 slen[0] = (scalefac_compress % 64) >> 4;
Gruenfrosch 0:7627c79db971 672 slen[1] = (scalefac_compress % 16) >> 2;
Gruenfrosch 0:7627c79db971 673 slen[2] = scalefac_compress % 4;
Gruenfrosch 0:7627c79db971 674 slen[3] = 0;
Gruenfrosch 0:7627c79db971 675
Gruenfrosch 0:7627c79db971 676 nsfb = nsfb_table[4][index];
Gruenfrosch 0:7627c79db971 677 }
Gruenfrosch 0:7627c79db971 678 else {
Gruenfrosch 0:7627c79db971 679 scalefac_compress -= 244;
Gruenfrosch 0:7627c79db971 680
Gruenfrosch 0:7627c79db971 681 slen[0] = scalefac_compress / 3;
Gruenfrosch 0:7627c79db971 682 slen[1] = scalefac_compress % 3;
Gruenfrosch 0:7627c79db971 683 slen[2] = 0;
Gruenfrosch 0:7627c79db971 684 slen[3] = 0;
Gruenfrosch 0:7627c79db971 685
Gruenfrosch 0:7627c79db971 686 nsfb = nsfb_table[5][index];
Gruenfrosch 0:7627c79db971 687 }
Gruenfrosch 0:7627c79db971 688
Gruenfrosch 0:7627c79db971 689 n = 0;
Gruenfrosch 0:7627c79db971 690 for (part = 0; part < 4; ++part) {
Gruenfrosch 0:7627c79db971 691 unsigned int max, is_pos;
Gruenfrosch 0:7627c79db971 692
Gruenfrosch 0:7627c79db971 693 max = (1 << slen[part]) - 1;
Gruenfrosch 0:7627c79db971 694
Gruenfrosch 0:7627c79db971 695 for (i = 0; i < nsfb[part]; ++i) {
Gruenfrosch 0:7627c79db971 696 is_pos = mad_bit_read(ptr, slen[part]);
Gruenfrosch 0:7627c79db971 697
Gruenfrosch 0:7627c79db971 698 channel->scalefac[n] = is_pos;
Gruenfrosch 0:7627c79db971 699 gr1ch->scalefac[n++] = (is_pos == max);
Gruenfrosch 0:7627c79db971 700 }
Gruenfrosch 0:7627c79db971 701 }
Gruenfrosch 0:7627c79db971 702
Gruenfrosch 0:7627c79db971 703 while (n < 39) {
Gruenfrosch 0:7627c79db971 704 channel->scalefac[n] = 0;
Gruenfrosch 0:7627c79db971 705 gr1ch->scalefac[n++] = 0; /* apparently not illegal */
Gruenfrosch 0:7627c79db971 706 }
Gruenfrosch 0:7627c79db971 707 }
Gruenfrosch 0:7627c79db971 708
Gruenfrosch 0:7627c79db971 709 return mad_bit_length(&start, ptr);
Gruenfrosch 0:7627c79db971 710 }
Gruenfrosch 0:7627c79db971 711
Gruenfrosch 0:7627c79db971 712 /*
Gruenfrosch 0:7627c79db971 713 * NAME: III_scalefactors()
Gruenfrosch 0:7627c79db971 714 * DESCRIPTION: decode channel scalefactors of one granule from a bitstream
Gruenfrosch 0:7627c79db971 715 */
Gruenfrosch 0:7627c79db971 716 static
Gruenfrosch 0:7627c79db971 717 unsigned int III_scalefactors(struct mad_bitptr *ptr, struct channel *channel,
Gruenfrosch 0:7627c79db971 718 struct channel const *gr0ch, unsigned int scfsi)
Gruenfrosch 0:7627c79db971 719 {
Gruenfrosch 0:7627c79db971 720 struct mad_bitptr start;
Gruenfrosch 0:7627c79db971 721 unsigned int slen1, slen2, sfbi;
Gruenfrosch 0:7627c79db971 722
Gruenfrosch 0:7627c79db971 723 start = *ptr;
Gruenfrosch 0:7627c79db971 724
Gruenfrosch 0:7627c79db971 725 slen1 = sflen_table[channel->scalefac_compress].slen1;
Gruenfrosch 0:7627c79db971 726 slen2 = sflen_table[channel->scalefac_compress].slen2;
Gruenfrosch 0:7627c79db971 727
Gruenfrosch 0:7627c79db971 728 if (channel->block_type == 2) {
Gruenfrosch 0:7627c79db971 729 unsigned int nsfb;
Gruenfrosch 0:7627c79db971 730
Gruenfrosch 0:7627c79db971 731 sfbi = 0;
Gruenfrosch 0:7627c79db971 732
Gruenfrosch 0:7627c79db971 733 nsfb = (channel->flags & mixed_block_flag) ? 8 + 3 * 3 : 6 * 3;
Gruenfrosch 0:7627c79db971 734 while (nsfb--)
Gruenfrosch 0:7627c79db971 735 channel->scalefac[sfbi++] = mad_bit_read(ptr, slen1);
Gruenfrosch 0:7627c79db971 736
Gruenfrosch 0:7627c79db971 737 nsfb = 6 * 3;
Gruenfrosch 0:7627c79db971 738 while (nsfb--)
Gruenfrosch 0:7627c79db971 739 channel->scalefac[sfbi++] = mad_bit_read(ptr, slen2);
Gruenfrosch 0:7627c79db971 740
Gruenfrosch 0:7627c79db971 741 nsfb = 1 * 3;
Gruenfrosch 0:7627c79db971 742 while (nsfb--)
Gruenfrosch 0:7627c79db971 743 channel->scalefac[sfbi++] = 0;
Gruenfrosch 0:7627c79db971 744 }
Gruenfrosch 0:7627c79db971 745 else { /* channel->block_type != 2 */
Gruenfrosch 0:7627c79db971 746 if (scfsi & 0x8) {
Gruenfrosch 0:7627c79db971 747 for (sfbi = 0; sfbi < 6; ++sfbi)
Gruenfrosch 0:7627c79db971 748 channel->scalefac[sfbi] = gr0ch->scalefac[sfbi];
Gruenfrosch 0:7627c79db971 749 }
Gruenfrosch 0:7627c79db971 750 else {
Gruenfrosch 0:7627c79db971 751 for (sfbi = 0; sfbi < 6; ++sfbi)
Gruenfrosch 0:7627c79db971 752 channel->scalefac[sfbi] = mad_bit_read(ptr, slen1);
Gruenfrosch 0:7627c79db971 753 }
Gruenfrosch 0:7627c79db971 754
Gruenfrosch 0:7627c79db971 755 if (scfsi & 0x4) {
Gruenfrosch 0:7627c79db971 756 for (sfbi = 6; sfbi < 11; ++sfbi)
Gruenfrosch 0:7627c79db971 757 channel->scalefac[sfbi] = gr0ch->scalefac[sfbi];
Gruenfrosch 0:7627c79db971 758 }
Gruenfrosch 0:7627c79db971 759 else {
Gruenfrosch 0:7627c79db971 760 for (sfbi = 6; sfbi < 11; ++sfbi)
Gruenfrosch 0:7627c79db971 761 channel->scalefac[sfbi] = mad_bit_read(ptr, slen1);
Gruenfrosch 0:7627c79db971 762 }
Gruenfrosch 0:7627c79db971 763
Gruenfrosch 0:7627c79db971 764 if (scfsi & 0x2) {
Gruenfrosch 0:7627c79db971 765 for (sfbi = 11; sfbi < 16; ++sfbi)
Gruenfrosch 0:7627c79db971 766 channel->scalefac[sfbi] = gr0ch->scalefac[sfbi];
Gruenfrosch 0:7627c79db971 767 }
Gruenfrosch 0:7627c79db971 768 else {
Gruenfrosch 0:7627c79db971 769 for (sfbi = 11; sfbi < 16; ++sfbi)
Gruenfrosch 0:7627c79db971 770 channel->scalefac[sfbi] = mad_bit_read(ptr, slen2);
Gruenfrosch 0:7627c79db971 771 }
Gruenfrosch 0:7627c79db971 772
Gruenfrosch 0:7627c79db971 773 if (scfsi & 0x1) {
Gruenfrosch 0:7627c79db971 774 for (sfbi = 16; sfbi < 21; ++sfbi)
Gruenfrosch 0:7627c79db971 775 channel->scalefac[sfbi] = gr0ch->scalefac[sfbi];
Gruenfrosch 0:7627c79db971 776 }
Gruenfrosch 0:7627c79db971 777 else {
Gruenfrosch 0:7627c79db971 778 for (sfbi = 16; sfbi < 21; ++sfbi)
Gruenfrosch 0:7627c79db971 779 channel->scalefac[sfbi] = mad_bit_read(ptr, slen2);
Gruenfrosch 0:7627c79db971 780 }
Gruenfrosch 0:7627c79db971 781
Gruenfrosch 0:7627c79db971 782 channel->scalefac[21] = 0;
Gruenfrosch 0:7627c79db971 783 }
Gruenfrosch 0:7627c79db971 784
Gruenfrosch 0:7627c79db971 785 return mad_bit_length(&start, ptr);
Gruenfrosch 0:7627c79db971 786 }
Gruenfrosch 0:7627c79db971 787
Gruenfrosch 0:7627c79db971 788 /*
Gruenfrosch 0:7627c79db971 789 * The Layer III formula for requantization and scaling is defined by
Gruenfrosch 0:7627c79db971 790 * section 2.4.3.4.7.1 of ISO/IEC 11172-3, as follows:
Gruenfrosch 0:7627c79db971 791 *
Gruenfrosch 0:7627c79db971 792 * long blocks:
Gruenfrosch 0:7627c79db971 793 * xr[i] = sign(is[i]) * abs(is[i])^(4/3) *
Gruenfrosch 0:7627c79db971 794 * 2^((1/4) * (global_gain - 210)) *
Gruenfrosch 0:7627c79db971 795 * 2^-(scalefac_multiplier *
Gruenfrosch 0:7627c79db971 796 * (scalefac_l[sfb] + preflag * pretab[sfb]))
Gruenfrosch 0:7627c79db971 797 *
Gruenfrosch 0:7627c79db971 798 * short blocks:
Gruenfrosch 0:7627c79db971 799 * xr[i] = sign(is[i]) * abs(is[i])^(4/3) *
Gruenfrosch 0:7627c79db971 800 * 2^((1/4) * (global_gain - 210 - 8 * subblock_gain[w])) *
Gruenfrosch 0:7627c79db971 801 * 2^-(scalefac_multiplier * scalefac_s[sfb][w])
Gruenfrosch 0:7627c79db971 802 *
Gruenfrosch 0:7627c79db971 803 * where:
Gruenfrosch 0:7627c79db971 804 * scalefac_multiplier = (scalefac_scale + 1) / 2
Gruenfrosch 0:7627c79db971 805 *
Gruenfrosch 0:7627c79db971 806 * The routines III_exponents() and III_requantize() facilitate this
Gruenfrosch 0:7627c79db971 807 * calculation.
Gruenfrosch 0:7627c79db971 808 */
Gruenfrosch 0:7627c79db971 809
Gruenfrosch 0:7627c79db971 810 /*
Gruenfrosch 0:7627c79db971 811 * NAME: III_exponents()
Gruenfrosch 0:7627c79db971 812 * DESCRIPTION: calculate scalefactor exponents
Gruenfrosch 0:7627c79db971 813 */
Gruenfrosch 0:7627c79db971 814 static
Gruenfrosch 0:7627c79db971 815 void III_exponents(struct channel const *channel,
Gruenfrosch 0:7627c79db971 816 unsigned char const *sfbwidth, signed int exponents[39])
Gruenfrosch 0:7627c79db971 817 {
Gruenfrosch 0:7627c79db971 818 signed int gain;
Gruenfrosch 0:7627c79db971 819 unsigned int scalefac_multiplier, sfbi;
Gruenfrosch 0:7627c79db971 820
Gruenfrosch 0:7627c79db971 821 gain = (signed int) channel->global_gain - 210;
Gruenfrosch 0:7627c79db971 822 scalefac_multiplier = (channel->flags & scalefac_scale) ? 2 : 1;
Gruenfrosch 0:7627c79db971 823
Gruenfrosch 0:7627c79db971 824 if (channel->block_type == 2) {
Gruenfrosch 0:7627c79db971 825 unsigned int l;
Gruenfrosch 0:7627c79db971 826 signed int gain0, gain1, gain2;
Gruenfrosch 0:7627c79db971 827
Gruenfrosch 0:7627c79db971 828 sfbi = l = 0;
Gruenfrosch 0:7627c79db971 829
Gruenfrosch 0:7627c79db971 830 if (channel->flags & mixed_block_flag) {
Gruenfrosch 0:7627c79db971 831 unsigned int premask;
Gruenfrosch 0:7627c79db971 832
Gruenfrosch 0:7627c79db971 833 premask = (channel->flags & preflag) ? ~0 : 0;
Gruenfrosch 0:7627c79db971 834
Gruenfrosch 0:7627c79db971 835 /* long block subbands 0-1 */
Gruenfrosch 0:7627c79db971 836
Gruenfrosch 0:7627c79db971 837 while (l < 36) {
Gruenfrosch 0:7627c79db971 838 exponents[sfbi] = gain -
Gruenfrosch 0:7627c79db971 839 (signed int) ((channel->scalefac[sfbi] + (pretab[sfbi] & premask)) <<
Gruenfrosch 0:7627c79db971 840 scalefac_multiplier);
Gruenfrosch 0:7627c79db971 841
Gruenfrosch 0:7627c79db971 842 l += sfbwidth[sfbi++];
Gruenfrosch 0:7627c79db971 843 }
Gruenfrosch 0:7627c79db971 844 }
Gruenfrosch 0:7627c79db971 845
Gruenfrosch 0:7627c79db971 846 /* this is probably wrong for 8000 Hz short/mixed blocks */
Gruenfrosch 0:7627c79db971 847
Gruenfrosch 0:7627c79db971 848 gain0 = gain - 8 * (signed int) channel->subblock_gain[0];
Gruenfrosch 0:7627c79db971 849 gain1 = gain - 8 * (signed int) channel->subblock_gain[1];
Gruenfrosch 0:7627c79db971 850 gain2 = gain - 8 * (signed int) channel->subblock_gain[2];
Gruenfrosch 0:7627c79db971 851
Gruenfrosch 0:7627c79db971 852 while (l < 576) {
Gruenfrosch 0:7627c79db971 853 exponents[sfbi + 0] = gain0 -
Gruenfrosch 0:7627c79db971 854 (signed int) (channel->scalefac[sfbi + 0] << scalefac_multiplier);
Gruenfrosch 0:7627c79db971 855 exponents[sfbi + 1] = gain1 -
Gruenfrosch 0:7627c79db971 856 (signed int) (channel->scalefac[sfbi + 1] << scalefac_multiplier);
Gruenfrosch 0:7627c79db971 857 exponents[sfbi + 2] = gain2 -
Gruenfrosch 0:7627c79db971 858 (signed int) (channel->scalefac[sfbi + 2] << scalefac_multiplier);
Gruenfrosch 0:7627c79db971 859
Gruenfrosch 0:7627c79db971 860 l += 3 * sfbwidth[sfbi];
Gruenfrosch 0:7627c79db971 861 sfbi += 3;
Gruenfrosch 0:7627c79db971 862 }
Gruenfrosch 0:7627c79db971 863 }
Gruenfrosch 0:7627c79db971 864 else { /* channel->block_type != 2 */
Gruenfrosch 0:7627c79db971 865 if (channel->flags & preflag) {
Gruenfrosch 0:7627c79db971 866 for (sfbi = 0; sfbi < 22; ++sfbi) {
Gruenfrosch 0:7627c79db971 867 exponents[sfbi] = gain -
Gruenfrosch 0:7627c79db971 868 (signed int) ((channel->scalefac[sfbi] + pretab[sfbi]) <<
Gruenfrosch 0:7627c79db971 869 scalefac_multiplier);
Gruenfrosch 0:7627c79db971 870 }
Gruenfrosch 0:7627c79db971 871 }
Gruenfrosch 0:7627c79db971 872 else {
Gruenfrosch 0:7627c79db971 873 for (sfbi = 0; sfbi < 22; ++sfbi) {
Gruenfrosch 0:7627c79db971 874 exponents[sfbi] = gain -
Gruenfrosch 0:7627c79db971 875 (signed int) (channel->scalefac[sfbi] << scalefac_multiplier);
Gruenfrosch 0:7627c79db971 876 }
Gruenfrosch 0:7627c79db971 877 }
Gruenfrosch 0:7627c79db971 878 }
Gruenfrosch 0:7627c79db971 879 }
Gruenfrosch 0:7627c79db971 880
Gruenfrosch 0:7627c79db971 881 /*
Gruenfrosch 0:7627c79db971 882 * NAME: III_requantize()
Gruenfrosch 0:7627c79db971 883 * DESCRIPTION: requantize one (positive) value
Gruenfrosch 0:7627c79db971 884 */
Gruenfrosch 0:7627c79db971 885 static
Gruenfrosch 0:7627c79db971 886 mad_fixed_t III_requantize(unsigned int value, signed int exp)
Gruenfrosch 0:7627c79db971 887 {
Gruenfrosch 0:7627c79db971 888 mad_fixed_t requantized;
Gruenfrosch 0:7627c79db971 889 signed int frac;
Gruenfrosch 0:7627c79db971 890 struct fixedfloat const *power;
Gruenfrosch 0:7627c79db971 891
Gruenfrosch 0:7627c79db971 892 frac = exp % 4; /* assumes sign(frac) == sign(exp) */
Gruenfrosch 0:7627c79db971 893 exp /= 4;
Gruenfrosch 0:7627c79db971 894
Gruenfrosch 0:7627c79db971 895 power = &rq_table[value];
Gruenfrosch 0:7627c79db971 896 requantized = power->mantissa;
Gruenfrosch 0:7627c79db971 897 exp += power->exponent;
Gruenfrosch 0:7627c79db971 898
Gruenfrosch 0:7627c79db971 899 if (exp < 0) {
Gruenfrosch 0:7627c79db971 900 if (-exp >= sizeof(mad_fixed_t) * CHAR_BIT) {
Gruenfrosch 0:7627c79db971 901 /* underflow */
Gruenfrosch 0:7627c79db971 902 requantized = 0;
Gruenfrosch 0:7627c79db971 903 }
Gruenfrosch 0:7627c79db971 904 else {
Gruenfrosch 0:7627c79db971 905 requantized += 1L << (-exp - 1);
Gruenfrosch 0:7627c79db971 906 requantized >>= -exp;
Gruenfrosch 0:7627c79db971 907 }
Gruenfrosch 0:7627c79db971 908 }
Gruenfrosch 0:7627c79db971 909 else {
Gruenfrosch 0:7627c79db971 910 if (exp >= 5) {
Gruenfrosch 0:7627c79db971 911 /* overflow */
Gruenfrosch 0:7627c79db971 912 # if defined(DEBUG)
Gruenfrosch 0:7627c79db971 913 fprintf(stderr, "requantize overflow (%f * 2^%d)\n",
Gruenfrosch 0:7627c79db971 914 mad_f_todouble(requantized), exp);
Gruenfrosch 0:7627c79db971 915 # endif
Gruenfrosch 0:7627c79db971 916 requantized = MAD_F_MAX;
Gruenfrosch 0:7627c79db971 917 }
Gruenfrosch 0:7627c79db971 918 else
Gruenfrosch 0:7627c79db971 919 requantized <<= exp;
Gruenfrosch 0:7627c79db971 920 }
Gruenfrosch 0:7627c79db971 921
Gruenfrosch 0:7627c79db971 922 return frac ? mad_f_mul(requantized, root_table[3 + frac]) : requantized;
Gruenfrosch 0:7627c79db971 923 }
Gruenfrosch 0:7627c79db971 924
Gruenfrosch 0:7627c79db971 925 /* we must take care that sz >= bits and sz < sizeof(long) lest bits == 0 */
Gruenfrosch 0:7627c79db971 926 # define MASK(cache, sz, bits) \
Gruenfrosch 0:7627c79db971 927 (((cache) >> ((sz) - (bits))) & ((1 << (bits)) - 1))
Gruenfrosch 0:7627c79db971 928 # define MASK1BIT(cache, sz) \
Gruenfrosch 0:7627c79db971 929 ((cache) & (1 << ((sz) - 1)))
Gruenfrosch 0:7627c79db971 930
Gruenfrosch 0:7627c79db971 931 /*
Gruenfrosch 0:7627c79db971 932 * NAME: III_huffdecode()
Gruenfrosch 0:7627c79db971 933 * DESCRIPTION: decode Huffman code words of one channel of one granule
Gruenfrosch 0:7627c79db971 934 */
Gruenfrosch 0:7627c79db971 935 static
Gruenfrosch 0:7627c79db971 936 enum mad_error III_huffdecode(struct mad_bitptr *ptr, mad_fixed_t xr[576],
Gruenfrosch 0:7627c79db971 937 struct channel *channel,
Gruenfrosch 0:7627c79db971 938 unsigned char const *sfbwidth,
Gruenfrosch 0:7627c79db971 939 unsigned int part2_length)
Gruenfrosch 0:7627c79db971 940 {
Gruenfrosch 0:7627c79db971 941 signed int exponents[39], exp;
Gruenfrosch 0:7627c79db971 942 signed int const *expptr;
Gruenfrosch 0:7627c79db971 943 struct mad_bitptr peek;
Gruenfrosch 0:7627c79db971 944 signed int bits_left, cachesz;
Gruenfrosch 0:7627c79db971 945 register mad_fixed_t *xrptr;
Gruenfrosch 0:7627c79db971 946 mad_fixed_t const *sfbound;
Gruenfrosch 0:7627c79db971 947 register unsigned long bitcache;
Gruenfrosch 0:7627c79db971 948
Gruenfrosch 0:7627c79db971 949 bits_left = (signed) channel->part2_3_length - (signed) part2_length;
Gruenfrosch 0:7627c79db971 950 if (bits_left < 0)
Gruenfrosch 0:7627c79db971 951 return MAD_ERROR_BADPART3LEN;
Gruenfrosch 0:7627c79db971 952
Gruenfrosch 0:7627c79db971 953 III_exponents(channel, sfbwidth, exponents);
Gruenfrosch 0:7627c79db971 954
Gruenfrosch 0:7627c79db971 955 peek = *ptr;
Gruenfrosch 0:7627c79db971 956 mad_bit_skip(ptr, bits_left);
Gruenfrosch 0:7627c79db971 957
Gruenfrosch 0:7627c79db971 958 /* align bit reads to byte boundaries */
Gruenfrosch 0:7627c79db971 959 cachesz = mad_bit_bitsleft(&peek);
Gruenfrosch 0:7627c79db971 960 cachesz += ((32 - 1 - 24) + (24 - cachesz)) & ~7;
Gruenfrosch 0:7627c79db971 961
Gruenfrosch 0:7627c79db971 962 bitcache = mad_bit_read(&peek, cachesz);
Gruenfrosch 0:7627c79db971 963 bits_left -= cachesz;
Gruenfrosch 0:7627c79db971 964
Gruenfrosch 0:7627c79db971 965 xrptr = &xr[0];
Gruenfrosch 0:7627c79db971 966
Gruenfrosch 0:7627c79db971 967 /* big_values */
Gruenfrosch 0:7627c79db971 968 {
Gruenfrosch 0:7627c79db971 969 unsigned int region, rcount;
Gruenfrosch 0:7627c79db971 970 struct hufftable const *entry;
Gruenfrosch 0:7627c79db971 971 union huffpair const *table;
Gruenfrosch 0:7627c79db971 972 unsigned int linbits, startbits, big_values, reqhits;
Gruenfrosch 0:7627c79db971 973 mad_fixed_t reqcache[16];
Gruenfrosch 0:7627c79db971 974
Gruenfrosch 0:7627c79db971 975 sfbound = xrptr + *sfbwidth++;
Gruenfrosch 0:7627c79db971 976 rcount = channel->region0_count + 1;
Gruenfrosch 0:7627c79db971 977
Gruenfrosch 0:7627c79db971 978 entry = &mad_huff_pair_table[channel->table_select[region = 0]];
Gruenfrosch 0:7627c79db971 979 table = entry->table;
Gruenfrosch 0:7627c79db971 980 linbits = entry->linbits;
Gruenfrosch 0:7627c79db971 981 startbits = entry->startbits;
Gruenfrosch 0:7627c79db971 982
Gruenfrosch 0:7627c79db971 983 if (table == 0)
Gruenfrosch 0:7627c79db971 984 return MAD_ERROR_BADHUFFTABLE;
Gruenfrosch 0:7627c79db971 985
Gruenfrosch 0:7627c79db971 986 expptr = &exponents[0];
Gruenfrosch 0:7627c79db971 987 exp = *expptr++;
Gruenfrosch 0:7627c79db971 988 reqhits = 0;
Gruenfrosch 0:7627c79db971 989
Gruenfrosch 0:7627c79db971 990 big_values = channel->big_values;
Gruenfrosch 0:7627c79db971 991
Gruenfrosch 0:7627c79db971 992 while (big_values-- && cachesz + bits_left > 0) {
Gruenfrosch 0:7627c79db971 993 union huffpair const *pair;
Gruenfrosch 0:7627c79db971 994 unsigned int clumpsz, value;
Gruenfrosch 0:7627c79db971 995 register mad_fixed_t requantized;
Gruenfrosch 0:7627c79db971 996
Gruenfrosch 0:7627c79db971 997 if (xrptr == sfbound) {
Gruenfrosch 0:7627c79db971 998 sfbound += *sfbwidth++;
Gruenfrosch 0:7627c79db971 999
Gruenfrosch 0:7627c79db971 1000 /* change table if region boundary */
Gruenfrosch 0:7627c79db971 1001
Gruenfrosch 0:7627c79db971 1002 if (--rcount == 0) {
Gruenfrosch 0:7627c79db971 1003 if (region == 0)
Gruenfrosch 0:7627c79db971 1004 rcount = channel->region1_count + 1;
Gruenfrosch 0:7627c79db971 1005 else
Gruenfrosch 0:7627c79db971 1006 rcount = 0; /* all remaining */
Gruenfrosch 0:7627c79db971 1007
Gruenfrosch 0:7627c79db971 1008 entry = &mad_huff_pair_table[channel->table_select[++region]];
Gruenfrosch 0:7627c79db971 1009 table = entry->table;
Gruenfrosch 0:7627c79db971 1010 linbits = entry->linbits;
Gruenfrosch 0:7627c79db971 1011 startbits = entry->startbits;
Gruenfrosch 0:7627c79db971 1012
Gruenfrosch 0:7627c79db971 1013 if (table == 0)
Gruenfrosch 0:7627c79db971 1014 return MAD_ERROR_BADHUFFTABLE;
Gruenfrosch 0:7627c79db971 1015 }
Gruenfrosch 0:7627c79db971 1016
Gruenfrosch 0:7627c79db971 1017 if (exp != *expptr) {
Gruenfrosch 0:7627c79db971 1018 exp = *expptr;
Gruenfrosch 0:7627c79db971 1019 reqhits = 0;
Gruenfrosch 0:7627c79db971 1020 }
Gruenfrosch 0:7627c79db971 1021
Gruenfrosch 0:7627c79db971 1022 ++expptr;
Gruenfrosch 0:7627c79db971 1023 }
Gruenfrosch 0:7627c79db971 1024
Gruenfrosch 0:7627c79db971 1025 if (cachesz < 21) {
Gruenfrosch 0:7627c79db971 1026 unsigned int bits;
Gruenfrosch 0:7627c79db971 1027
Gruenfrosch 0:7627c79db971 1028 bits = ((32 - 1 - 21) + (21 - cachesz)) & ~7;
Gruenfrosch 0:7627c79db971 1029 bitcache = (bitcache << bits) | mad_bit_read(&peek, bits);
Gruenfrosch 0:7627c79db971 1030 cachesz += bits;
Gruenfrosch 0:7627c79db971 1031 bits_left -= bits;
Gruenfrosch 0:7627c79db971 1032 }
Gruenfrosch 0:7627c79db971 1033
Gruenfrosch 0:7627c79db971 1034 /* hcod (0..19) */
Gruenfrosch 0:7627c79db971 1035
Gruenfrosch 0:7627c79db971 1036 clumpsz = startbits;
Gruenfrosch 0:7627c79db971 1037 pair = &table[MASK(bitcache, cachesz, clumpsz)];
Gruenfrosch 0:7627c79db971 1038
Gruenfrosch 0:7627c79db971 1039 while (!pair->final) {
Gruenfrosch 0:7627c79db971 1040 cachesz -= clumpsz;
Gruenfrosch 0:7627c79db971 1041
Gruenfrosch 0:7627c79db971 1042 clumpsz = pair->ptr.bits;
Gruenfrosch 0:7627c79db971 1043 pair = &table[pair->ptr.offset + MASK(bitcache, cachesz, clumpsz)];
Gruenfrosch 0:7627c79db971 1044 }
Gruenfrosch 0:7627c79db971 1045
Gruenfrosch 0:7627c79db971 1046 cachesz -= pair->value.hlen;
Gruenfrosch 0:7627c79db971 1047
Gruenfrosch 0:7627c79db971 1048 if (linbits) {
Gruenfrosch 0:7627c79db971 1049 /* x (0..14) */
Gruenfrosch 0:7627c79db971 1050
Gruenfrosch 0:7627c79db971 1051 value = pair->value.x;
Gruenfrosch 0:7627c79db971 1052
Gruenfrosch 0:7627c79db971 1053 switch (value) {
Gruenfrosch 0:7627c79db971 1054 case 0:
Gruenfrosch 0:7627c79db971 1055 xrptr[0] = 0;
Gruenfrosch 0:7627c79db971 1056 break;
Gruenfrosch 0:7627c79db971 1057
Gruenfrosch 0:7627c79db971 1058 case 15:
Gruenfrosch 0:7627c79db971 1059 if (cachesz < linbits + 2) {
Gruenfrosch 0:7627c79db971 1060 bitcache = (bitcache << 16) | mad_bit_read(&peek, 16);
Gruenfrosch 0:7627c79db971 1061 cachesz += 16;
Gruenfrosch 0:7627c79db971 1062 bits_left -= 16;
Gruenfrosch 0:7627c79db971 1063 }
Gruenfrosch 0:7627c79db971 1064
Gruenfrosch 0:7627c79db971 1065 value += MASK(bitcache, cachesz, linbits);
Gruenfrosch 0:7627c79db971 1066 cachesz -= linbits;
Gruenfrosch 0:7627c79db971 1067
Gruenfrosch 0:7627c79db971 1068 requantized = III_requantize(value, exp);
Gruenfrosch 0:7627c79db971 1069 goto x_final;
Gruenfrosch 0:7627c79db971 1070
Gruenfrosch 0:7627c79db971 1071 default:
Gruenfrosch 0:7627c79db971 1072 if (reqhits & (1 << value))
Gruenfrosch 0:7627c79db971 1073 requantized = reqcache[value];
Gruenfrosch 0:7627c79db971 1074 else {
Gruenfrosch 0:7627c79db971 1075 reqhits |= (1 << value);
Gruenfrosch 0:7627c79db971 1076 requantized = reqcache[value] = III_requantize(value, exp);
Gruenfrosch 0:7627c79db971 1077 }
Gruenfrosch 0:7627c79db971 1078
Gruenfrosch 0:7627c79db971 1079 x_final:
Gruenfrosch 0:7627c79db971 1080 xrptr[0] = MASK1BIT(bitcache, cachesz--) ?
Gruenfrosch 0:7627c79db971 1081 -requantized : requantized;
Gruenfrosch 0:7627c79db971 1082 }
Gruenfrosch 0:7627c79db971 1083
Gruenfrosch 0:7627c79db971 1084 /* y (0..14) */
Gruenfrosch 0:7627c79db971 1085
Gruenfrosch 0:7627c79db971 1086 value = pair->value.y;
Gruenfrosch 0:7627c79db971 1087
Gruenfrosch 0:7627c79db971 1088 switch (value) {
Gruenfrosch 0:7627c79db971 1089 case 0:
Gruenfrosch 0:7627c79db971 1090 xrptr[1] = 0;
Gruenfrosch 0:7627c79db971 1091 break;
Gruenfrosch 0:7627c79db971 1092
Gruenfrosch 0:7627c79db971 1093 case 15:
Gruenfrosch 0:7627c79db971 1094 if (cachesz < linbits + 1) {
Gruenfrosch 0:7627c79db971 1095 bitcache = (bitcache << 16) | mad_bit_read(&peek, 16);
Gruenfrosch 0:7627c79db971 1096 cachesz += 16;
Gruenfrosch 0:7627c79db971 1097 bits_left -= 16;
Gruenfrosch 0:7627c79db971 1098 }
Gruenfrosch 0:7627c79db971 1099
Gruenfrosch 0:7627c79db971 1100 value += MASK(bitcache, cachesz, linbits);
Gruenfrosch 0:7627c79db971 1101 cachesz -= linbits;
Gruenfrosch 0:7627c79db971 1102
Gruenfrosch 0:7627c79db971 1103 requantized = III_requantize(value, exp);
Gruenfrosch 0:7627c79db971 1104 goto y_final;
Gruenfrosch 0:7627c79db971 1105
Gruenfrosch 0:7627c79db971 1106 default:
Gruenfrosch 0:7627c79db971 1107 if (reqhits & (1 << value))
Gruenfrosch 0:7627c79db971 1108 requantized = reqcache[value];
Gruenfrosch 0:7627c79db971 1109 else {
Gruenfrosch 0:7627c79db971 1110 reqhits |= (1 << value);
Gruenfrosch 0:7627c79db971 1111 requantized = reqcache[value] = III_requantize(value, exp);
Gruenfrosch 0:7627c79db971 1112 }
Gruenfrosch 0:7627c79db971 1113
Gruenfrosch 0:7627c79db971 1114 y_final:
Gruenfrosch 0:7627c79db971 1115 xrptr[1] = MASK1BIT(bitcache, cachesz--) ?
Gruenfrosch 0:7627c79db971 1116 -requantized : requantized;
Gruenfrosch 0:7627c79db971 1117 }
Gruenfrosch 0:7627c79db971 1118 }
Gruenfrosch 0:7627c79db971 1119 else {
Gruenfrosch 0:7627c79db971 1120 /* x (0..1) */
Gruenfrosch 0:7627c79db971 1121
Gruenfrosch 0:7627c79db971 1122 value = pair->value.x;
Gruenfrosch 0:7627c79db971 1123
Gruenfrosch 0:7627c79db971 1124 if (value == 0)
Gruenfrosch 0:7627c79db971 1125 xrptr[0] = 0;
Gruenfrosch 0:7627c79db971 1126 else {
Gruenfrosch 0:7627c79db971 1127 if (reqhits & (1 << value))
Gruenfrosch 0:7627c79db971 1128 requantized = reqcache[value];
Gruenfrosch 0:7627c79db971 1129 else {
Gruenfrosch 0:7627c79db971 1130 reqhits |= (1 << value);
Gruenfrosch 0:7627c79db971 1131 requantized = reqcache[value] = III_requantize(value, exp);
Gruenfrosch 0:7627c79db971 1132 }
Gruenfrosch 0:7627c79db971 1133
Gruenfrosch 0:7627c79db971 1134 xrptr[0] = MASK1BIT(bitcache, cachesz--) ?
Gruenfrosch 0:7627c79db971 1135 -requantized : requantized;
Gruenfrosch 0:7627c79db971 1136 }
Gruenfrosch 0:7627c79db971 1137
Gruenfrosch 0:7627c79db971 1138 /* y (0..1) */
Gruenfrosch 0:7627c79db971 1139
Gruenfrosch 0:7627c79db971 1140 value = pair->value.y;
Gruenfrosch 0:7627c79db971 1141
Gruenfrosch 0:7627c79db971 1142 if (value == 0)
Gruenfrosch 0:7627c79db971 1143 xrptr[1] = 0;
Gruenfrosch 0:7627c79db971 1144 else {
Gruenfrosch 0:7627c79db971 1145 if (reqhits & (1 << value))
Gruenfrosch 0:7627c79db971 1146 requantized = reqcache[value];
Gruenfrosch 0:7627c79db971 1147 else {
Gruenfrosch 0:7627c79db971 1148 reqhits |= (1 << value);
Gruenfrosch 0:7627c79db971 1149 requantized = reqcache[value] = III_requantize(value, exp);
Gruenfrosch 0:7627c79db971 1150 }
Gruenfrosch 0:7627c79db971 1151
Gruenfrosch 0:7627c79db971 1152 xrptr[1] = MASK1BIT(bitcache, cachesz--) ?
Gruenfrosch 0:7627c79db971 1153 -requantized : requantized;
Gruenfrosch 0:7627c79db971 1154 }
Gruenfrosch 0:7627c79db971 1155 }
Gruenfrosch 0:7627c79db971 1156
Gruenfrosch 0:7627c79db971 1157 xrptr += 2;
Gruenfrosch 0:7627c79db971 1158 }
Gruenfrosch 0:7627c79db971 1159 }
Gruenfrosch 0:7627c79db971 1160
Gruenfrosch 0:7627c79db971 1161 if (cachesz + bits_left < 0)
Gruenfrosch 0:7627c79db971 1162 return MAD_ERROR_BADHUFFDATA; /* big_values overrun */
Gruenfrosch 0:7627c79db971 1163
Gruenfrosch 0:7627c79db971 1164 /* count1 */
Gruenfrosch 0:7627c79db971 1165 {
Gruenfrosch 0:7627c79db971 1166 union huffquad const *table;
Gruenfrosch 0:7627c79db971 1167 register mad_fixed_t requantized;
Gruenfrosch 0:7627c79db971 1168
Gruenfrosch 0:7627c79db971 1169 table = mad_huff_quad_table[channel->flags & count1table_select];
Gruenfrosch 0:7627c79db971 1170
Gruenfrosch 0:7627c79db971 1171 requantized = III_requantize(1, exp);
Gruenfrosch 0:7627c79db971 1172
Gruenfrosch 0:7627c79db971 1173 while (cachesz + bits_left > 0 && xrptr <= &xr[572]) {
Gruenfrosch 0:7627c79db971 1174 union huffquad const *quad;
Gruenfrosch 0:7627c79db971 1175
Gruenfrosch 0:7627c79db971 1176 /* hcod (1..6) */
Gruenfrosch 0:7627c79db971 1177
Gruenfrosch 0:7627c79db971 1178 if (cachesz < 10) {
Gruenfrosch 0:7627c79db971 1179 bitcache = (bitcache << 16) | mad_bit_read(&peek, 16);
Gruenfrosch 0:7627c79db971 1180 cachesz += 16;
Gruenfrosch 0:7627c79db971 1181 bits_left -= 16;
Gruenfrosch 0:7627c79db971 1182 }
Gruenfrosch 0:7627c79db971 1183
Gruenfrosch 0:7627c79db971 1184 quad = &table[MASK(bitcache, cachesz, 4)];
Gruenfrosch 0:7627c79db971 1185
Gruenfrosch 0:7627c79db971 1186 /* quad tables guaranteed to have at most one extra lookup */
Gruenfrosch 0:7627c79db971 1187 if (!quad->final) {
Gruenfrosch 0:7627c79db971 1188 cachesz -= 4;
Gruenfrosch 0:7627c79db971 1189
Gruenfrosch 0:7627c79db971 1190 quad = &table[quad->ptr.offset +
Gruenfrosch 0:7627c79db971 1191 MASK(bitcache, cachesz, quad->ptr.bits)];
Gruenfrosch 0:7627c79db971 1192 }
Gruenfrosch 0:7627c79db971 1193
Gruenfrosch 0:7627c79db971 1194 cachesz -= quad->value.hlen;
Gruenfrosch 0:7627c79db971 1195
Gruenfrosch 0:7627c79db971 1196 if (xrptr == sfbound) {
Gruenfrosch 0:7627c79db971 1197 sfbound += *sfbwidth++;
Gruenfrosch 0:7627c79db971 1198
Gruenfrosch 0:7627c79db971 1199 if (exp != *expptr) {
Gruenfrosch 0:7627c79db971 1200 exp = *expptr;
Gruenfrosch 0:7627c79db971 1201 requantized = III_requantize(1, exp);
Gruenfrosch 0:7627c79db971 1202 }
Gruenfrosch 0:7627c79db971 1203
Gruenfrosch 0:7627c79db971 1204 ++expptr;
Gruenfrosch 0:7627c79db971 1205 }
Gruenfrosch 0:7627c79db971 1206
Gruenfrosch 0:7627c79db971 1207 /* v (0..1) */
Gruenfrosch 0:7627c79db971 1208
Gruenfrosch 0:7627c79db971 1209 xrptr[0] = quad->value.v ?
Gruenfrosch 0:7627c79db971 1210 (MASK1BIT(bitcache, cachesz--) ? -requantized : requantized) : 0;
Gruenfrosch 0:7627c79db971 1211
Gruenfrosch 0:7627c79db971 1212 /* w (0..1) */
Gruenfrosch 0:7627c79db971 1213
Gruenfrosch 0:7627c79db971 1214 xrptr[1] = quad->value.w ?
Gruenfrosch 0:7627c79db971 1215 (MASK1BIT(bitcache, cachesz--) ? -requantized : requantized) : 0;
Gruenfrosch 0:7627c79db971 1216
Gruenfrosch 0:7627c79db971 1217 xrptr += 2;
Gruenfrosch 0:7627c79db971 1218
Gruenfrosch 0:7627c79db971 1219 if (xrptr == sfbound) {
Gruenfrosch 0:7627c79db971 1220 sfbound += *sfbwidth++;
Gruenfrosch 0:7627c79db971 1221
Gruenfrosch 0:7627c79db971 1222 if (exp != *expptr) {
Gruenfrosch 0:7627c79db971 1223 exp = *expptr;
Gruenfrosch 0:7627c79db971 1224 requantized = III_requantize(1, exp);
Gruenfrosch 0:7627c79db971 1225 }
Gruenfrosch 0:7627c79db971 1226
Gruenfrosch 0:7627c79db971 1227 ++expptr;
Gruenfrosch 0:7627c79db971 1228 }
Gruenfrosch 0:7627c79db971 1229
Gruenfrosch 0:7627c79db971 1230 /* x (0..1) */
Gruenfrosch 0:7627c79db971 1231
Gruenfrosch 0:7627c79db971 1232 xrptr[0] = quad->value.x ?
Gruenfrosch 0:7627c79db971 1233 (MASK1BIT(bitcache, cachesz--) ? -requantized : requantized) : 0;
Gruenfrosch 0:7627c79db971 1234
Gruenfrosch 0:7627c79db971 1235 /* y (0..1) */
Gruenfrosch 0:7627c79db971 1236
Gruenfrosch 0:7627c79db971 1237 xrptr[1] = quad->value.y ?
Gruenfrosch 0:7627c79db971 1238 (MASK1BIT(bitcache, cachesz--) ? -requantized : requantized) : 0;
Gruenfrosch 0:7627c79db971 1239
Gruenfrosch 0:7627c79db971 1240 xrptr += 2;
Gruenfrosch 0:7627c79db971 1241 }
Gruenfrosch 0:7627c79db971 1242
Gruenfrosch 0:7627c79db971 1243 if (cachesz + bits_left < 0) {
Gruenfrosch 0:7627c79db971 1244 # if 0 && defined(DEBUG)
Gruenfrosch 0:7627c79db971 1245 fprintf(stderr, "huffman count1 overrun (%d bits)\n",
Gruenfrosch 0:7627c79db971 1246 -(cachesz + bits_left));
Gruenfrosch 0:7627c79db971 1247 # endif
Gruenfrosch 0:7627c79db971 1248
Gruenfrosch 0:7627c79db971 1249 /* technically the bitstream is misformatted, but apparently
Gruenfrosch 0:7627c79db971 1250 some encoders are just a bit sloppy with stuffing bits */
Gruenfrosch 0:7627c79db971 1251
Gruenfrosch 0:7627c79db971 1252 xrptr -= 4;
Gruenfrosch 0:7627c79db971 1253 }
Gruenfrosch 0:7627c79db971 1254 }
Gruenfrosch 0:7627c79db971 1255
Gruenfrosch 0:7627c79db971 1256 assert(-bits_left <= MAD_BUFFER_GUARD * CHAR_BIT);
Gruenfrosch 0:7627c79db971 1257
Gruenfrosch 0:7627c79db971 1258 # if 0 && defined(DEBUG)
Gruenfrosch 0:7627c79db971 1259 if (bits_left < 0)
Gruenfrosch 0:7627c79db971 1260 fprintf(stderr, "read %d bits too many\n", -bits_left);
Gruenfrosch 0:7627c79db971 1261 else if (cachesz + bits_left > 0)
Gruenfrosch 0:7627c79db971 1262 fprintf(stderr, "%d stuffing bits\n", cachesz + bits_left);
Gruenfrosch 0:7627c79db971 1263 # endif
Gruenfrosch 0:7627c79db971 1264
Gruenfrosch 0:7627c79db971 1265 /* rzero */
Gruenfrosch 0:7627c79db971 1266 while (xrptr < &xr[576]) {
Gruenfrosch 0:7627c79db971 1267 xrptr[0] = 0;
Gruenfrosch 0:7627c79db971 1268 xrptr[1] = 0;
Gruenfrosch 0:7627c79db971 1269
Gruenfrosch 0:7627c79db971 1270 xrptr += 2;
Gruenfrosch 0:7627c79db971 1271 }
Gruenfrosch 0:7627c79db971 1272
Gruenfrosch 0:7627c79db971 1273 return MAD_ERROR_NONE;
Gruenfrosch 0:7627c79db971 1274 }
Gruenfrosch 0:7627c79db971 1275
Gruenfrosch 0:7627c79db971 1276 # undef MASK
Gruenfrosch 0:7627c79db971 1277 # undef MASK1BIT
Gruenfrosch 0:7627c79db971 1278
Gruenfrosch 0:7627c79db971 1279 /*
Gruenfrosch 0:7627c79db971 1280 * NAME: III_reorder()
Gruenfrosch 0:7627c79db971 1281 * DESCRIPTION: reorder frequency lines of a short block into subband order
Gruenfrosch 0:7627c79db971 1282 */
Gruenfrosch 0:7627c79db971 1283 static
Gruenfrosch 0:7627c79db971 1284 void III_reorder(mad_fixed_t xr[576], struct channel const *channel,
Gruenfrosch 0:7627c79db971 1285 unsigned char const sfbwidth[39])
Gruenfrosch 0:7627c79db971 1286 {
Gruenfrosch 0:7627c79db971 1287 mad_fixed_t tmp[32][3][6];
Gruenfrosch 0:7627c79db971 1288 unsigned int sb, l, f, w, sbw[3], sw[3];
Gruenfrosch 0:7627c79db971 1289
Gruenfrosch 0:7627c79db971 1290 /* this is probably wrong for 8000 Hz mixed blocks */
Gruenfrosch 0:7627c79db971 1291
Gruenfrosch 0:7627c79db971 1292 sb = 0;
Gruenfrosch 0:7627c79db971 1293 if (channel->flags & mixed_block_flag) {
Gruenfrosch 0:7627c79db971 1294 sb = 2;
Gruenfrosch 0:7627c79db971 1295
Gruenfrosch 0:7627c79db971 1296 l = 0;
Gruenfrosch 0:7627c79db971 1297 while (l < 36)
Gruenfrosch 0:7627c79db971 1298 l += *sfbwidth++;
Gruenfrosch 0:7627c79db971 1299 }
Gruenfrosch 0:7627c79db971 1300
Gruenfrosch 0:7627c79db971 1301 for (w = 0; w < 3; ++w) {
Gruenfrosch 0:7627c79db971 1302 sbw[w] = sb;
Gruenfrosch 0:7627c79db971 1303 sw[w] = 0;
Gruenfrosch 0:7627c79db971 1304 }
Gruenfrosch 0:7627c79db971 1305
Gruenfrosch 0:7627c79db971 1306 f = *sfbwidth++;
Gruenfrosch 0:7627c79db971 1307 w = 0;
Gruenfrosch 0:7627c79db971 1308
Gruenfrosch 0:7627c79db971 1309 for (l = 18 * sb; l < 576; ++l) {
Gruenfrosch 0:7627c79db971 1310 if (f-- == 0) {
Gruenfrosch 0:7627c79db971 1311 f = *sfbwidth++ - 1;
Gruenfrosch 0:7627c79db971 1312 w = (w + 1) % 3;
Gruenfrosch 0:7627c79db971 1313 }
Gruenfrosch 0:7627c79db971 1314
Gruenfrosch 0:7627c79db971 1315 tmp[sbw[w]][w][sw[w]++] = xr[l];
Gruenfrosch 0:7627c79db971 1316
Gruenfrosch 0:7627c79db971 1317 if (sw[w] == 6) {
Gruenfrosch 0:7627c79db971 1318 sw[w] = 0;
Gruenfrosch 0:7627c79db971 1319 ++sbw[w];
Gruenfrosch 0:7627c79db971 1320 }
Gruenfrosch 0:7627c79db971 1321 }
Gruenfrosch 0:7627c79db971 1322
Gruenfrosch 0:7627c79db971 1323 memcpy(&xr[18 * sb], &tmp[sb], (576 - 18 * sb) * sizeof(mad_fixed_t));
Gruenfrosch 0:7627c79db971 1324 }
Gruenfrosch 0:7627c79db971 1325
Gruenfrosch 0:7627c79db971 1326 /*
Gruenfrosch 0:7627c79db971 1327 * NAME: III_stereo()
Gruenfrosch 0:7627c79db971 1328 * DESCRIPTION: perform joint stereo processing on a granule
Gruenfrosch 0:7627c79db971 1329 */
Gruenfrosch 0:7627c79db971 1330 static
Gruenfrosch 0:7627c79db971 1331 enum mad_error III_stereo(mad_fixed_t xr[2][576],
Gruenfrosch 0:7627c79db971 1332 struct granule const *granule,
Gruenfrosch 0:7627c79db971 1333 struct mad_header *header,
Gruenfrosch 0:7627c79db971 1334 unsigned char const *sfbwidth)
Gruenfrosch 0:7627c79db971 1335 {
Gruenfrosch 0:7627c79db971 1336 short modes[39];
Gruenfrosch 0:7627c79db971 1337 unsigned int sfbi, l, n, i;
Gruenfrosch 0:7627c79db971 1338
Gruenfrosch 0:7627c79db971 1339 if (granule->ch[0].block_type !=
Gruenfrosch 0:7627c79db971 1340 granule->ch[1].block_type ||
Gruenfrosch 0:7627c79db971 1341 (granule->ch[0].flags & mixed_block_flag) !=
Gruenfrosch 0:7627c79db971 1342 (granule->ch[1].flags & mixed_block_flag))
Gruenfrosch 0:7627c79db971 1343 return MAD_ERROR_BADSTEREO;
Gruenfrosch 0:7627c79db971 1344
Gruenfrosch 0:7627c79db971 1345 for (i = 0; i < 39; ++i)
Gruenfrosch 0:7627c79db971 1346 modes[i] = header->mode_extension;
Gruenfrosch 0:7627c79db971 1347
Gruenfrosch 0:7627c79db971 1348 /* intensity stereo */
Gruenfrosch 0:7627c79db971 1349
Gruenfrosch 0:7627c79db971 1350 if (header->mode_extension & I_STEREO) {
Gruenfrosch 0:7627c79db971 1351 struct channel const *right_ch = &granule->ch[1];
Gruenfrosch 0:7627c79db971 1352 mad_fixed_t const *right_xr = xr[1];
Gruenfrosch 0:7627c79db971 1353 unsigned int is_pos;
Gruenfrosch 0:7627c79db971 1354
Gruenfrosch 0:7627c79db971 1355 header->flags |= MAD_FLAG_I_STEREO;
Gruenfrosch 0:7627c79db971 1356
Gruenfrosch 0:7627c79db971 1357 /* first determine which scalefactor bands are to be processed */
Gruenfrosch 0:7627c79db971 1358
Gruenfrosch 0:7627c79db971 1359 if (right_ch->block_type == 2) {
Gruenfrosch 0:7627c79db971 1360 unsigned int lower, start, max, bound[3], w;
Gruenfrosch 0:7627c79db971 1361
Gruenfrosch 0:7627c79db971 1362 lower = start = max = bound[0] = bound[1] = bound[2] = 0;
Gruenfrosch 0:7627c79db971 1363
Gruenfrosch 0:7627c79db971 1364 sfbi = l = 0;
Gruenfrosch 0:7627c79db971 1365
Gruenfrosch 0:7627c79db971 1366 if (right_ch->flags & mixed_block_flag) {
Gruenfrosch 0:7627c79db971 1367 while (l < 36) {
Gruenfrosch 0:7627c79db971 1368 n = sfbwidth[sfbi++];
Gruenfrosch 0:7627c79db971 1369
Gruenfrosch 0:7627c79db971 1370 for (i = 0; i < n; ++i) {
Gruenfrosch 0:7627c79db971 1371 if (right_xr[i]) {
Gruenfrosch 0:7627c79db971 1372 lower = sfbi;
Gruenfrosch 0:7627c79db971 1373 break;
Gruenfrosch 0:7627c79db971 1374 }
Gruenfrosch 0:7627c79db971 1375 }
Gruenfrosch 0:7627c79db971 1376
Gruenfrosch 0:7627c79db971 1377 right_xr += n;
Gruenfrosch 0:7627c79db971 1378 l += n;
Gruenfrosch 0:7627c79db971 1379 }
Gruenfrosch 0:7627c79db971 1380
Gruenfrosch 0:7627c79db971 1381 start = sfbi;
Gruenfrosch 0:7627c79db971 1382 }
Gruenfrosch 0:7627c79db971 1383
Gruenfrosch 0:7627c79db971 1384 w = 0;
Gruenfrosch 0:7627c79db971 1385 while (l < 576) {
Gruenfrosch 0:7627c79db971 1386 n = sfbwidth[sfbi++];
Gruenfrosch 0:7627c79db971 1387
Gruenfrosch 0:7627c79db971 1388 for (i = 0; i < n; ++i) {
Gruenfrosch 0:7627c79db971 1389 if (right_xr[i]) {
Gruenfrosch 0:7627c79db971 1390 max = bound[w] = sfbi;
Gruenfrosch 0:7627c79db971 1391 break;
Gruenfrosch 0:7627c79db971 1392 }
Gruenfrosch 0:7627c79db971 1393 }
Gruenfrosch 0:7627c79db971 1394
Gruenfrosch 0:7627c79db971 1395 right_xr += n;
Gruenfrosch 0:7627c79db971 1396 l += n;
Gruenfrosch 0:7627c79db971 1397 w = (w + 1) % 3;
Gruenfrosch 0:7627c79db971 1398 }
Gruenfrosch 0:7627c79db971 1399
Gruenfrosch 0:7627c79db971 1400 if (max)
Gruenfrosch 0:7627c79db971 1401 lower = start;
Gruenfrosch 0:7627c79db971 1402
Gruenfrosch 0:7627c79db971 1403 /* long blocks */
Gruenfrosch 0:7627c79db971 1404
Gruenfrosch 0:7627c79db971 1405 for (i = 0; i < lower; ++i)
Gruenfrosch 0:7627c79db971 1406 modes[i] = header->mode_extension & ~I_STEREO;
Gruenfrosch 0:7627c79db971 1407
Gruenfrosch 0:7627c79db971 1408 /* short blocks */
Gruenfrosch 0:7627c79db971 1409
Gruenfrosch 0:7627c79db971 1410 w = 0;
Gruenfrosch 0:7627c79db971 1411 for (i = start; i < max; ++i) {
Gruenfrosch 0:7627c79db971 1412 if (i < bound[w])
Gruenfrosch 0:7627c79db971 1413 modes[i] = header->mode_extension & ~I_STEREO;
Gruenfrosch 0:7627c79db971 1414
Gruenfrosch 0:7627c79db971 1415 w = (w + 1) % 3;
Gruenfrosch 0:7627c79db971 1416 }
Gruenfrosch 0:7627c79db971 1417 }
Gruenfrosch 0:7627c79db971 1418 else { /* right_ch->block_type != 2 */
Gruenfrosch 0:7627c79db971 1419 unsigned int bound;
Gruenfrosch 0:7627c79db971 1420
Gruenfrosch 0:7627c79db971 1421 bound = 0;
Gruenfrosch 0:7627c79db971 1422 for (sfbi = l = 0; l < 576; l += n) {
Gruenfrosch 0:7627c79db971 1423 n = sfbwidth[sfbi++];
Gruenfrosch 0:7627c79db971 1424
Gruenfrosch 0:7627c79db971 1425 for (i = 0; i < n; ++i) {
Gruenfrosch 0:7627c79db971 1426 if (right_xr[i]) {
Gruenfrosch 0:7627c79db971 1427 bound = sfbi;
Gruenfrosch 0:7627c79db971 1428 break;
Gruenfrosch 0:7627c79db971 1429 }
Gruenfrosch 0:7627c79db971 1430 }
Gruenfrosch 0:7627c79db971 1431
Gruenfrosch 0:7627c79db971 1432 right_xr += n;
Gruenfrosch 0:7627c79db971 1433 }
Gruenfrosch 0:7627c79db971 1434
Gruenfrosch 0:7627c79db971 1435 for (i = 0; i < bound; ++i)
Gruenfrosch 0:7627c79db971 1436 modes[i] = header->mode_extension & ~I_STEREO;
Gruenfrosch 0:7627c79db971 1437 }
Gruenfrosch 0:7627c79db971 1438
Gruenfrosch 0:7627c79db971 1439 /* now do the actual processing */
Gruenfrosch 0:7627c79db971 1440
Gruenfrosch 0:7627c79db971 1441 if (header->flags & MAD_FLAG_LSF_EXT) {
Gruenfrosch 0:7627c79db971 1442 unsigned char const *illegal_pos = granule[1].ch[1].scalefac;
Gruenfrosch 0:7627c79db971 1443 mad_fixed_t const *lsf_scale;
Gruenfrosch 0:7627c79db971 1444
Gruenfrosch 0:7627c79db971 1445 /* intensity_scale */
Gruenfrosch 0:7627c79db971 1446 lsf_scale = is_lsf_table[right_ch->scalefac_compress & 0x1];
Gruenfrosch 0:7627c79db971 1447
Gruenfrosch 0:7627c79db971 1448 for (sfbi = l = 0; l < 576; ++sfbi, l += n) {
Gruenfrosch 0:7627c79db971 1449 n = sfbwidth[sfbi];
Gruenfrosch 0:7627c79db971 1450
Gruenfrosch 0:7627c79db971 1451 if (!(modes[sfbi] & I_STEREO))
Gruenfrosch 0:7627c79db971 1452 continue;
Gruenfrosch 0:7627c79db971 1453
Gruenfrosch 0:7627c79db971 1454 if (illegal_pos[sfbi]) {
Gruenfrosch 0:7627c79db971 1455 modes[sfbi] &= ~I_STEREO;
Gruenfrosch 0:7627c79db971 1456 continue;
Gruenfrosch 0:7627c79db971 1457 }
Gruenfrosch 0:7627c79db971 1458
Gruenfrosch 0:7627c79db971 1459 is_pos = right_ch->scalefac[sfbi];
Gruenfrosch 0:7627c79db971 1460
Gruenfrosch 0:7627c79db971 1461 for (i = 0; i < n; ++i) {
Gruenfrosch 0:7627c79db971 1462 register mad_fixed_t left;
Gruenfrosch 0:7627c79db971 1463
Gruenfrosch 0:7627c79db971 1464 left = xr[0][l + i];
Gruenfrosch 0:7627c79db971 1465
Gruenfrosch 0:7627c79db971 1466 if (is_pos == 0)
Gruenfrosch 0:7627c79db971 1467 xr[1][l + i] = left;
Gruenfrosch 0:7627c79db971 1468 else {
Gruenfrosch 0:7627c79db971 1469 register mad_fixed_t opposite;
Gruenfrosch 0:7627c79db971 1470
Gruenfrosch 0:7627c79db971 1471 opposite = mad_f_mul(left, lsf_scale[(is_pos - 1) / 2]);
Gruenfrosch 0:7627c79db971 1472
Gruenfrosch 0:7627c79db971 1473 if (is_pos & 1) {
Gruenfrosch 0:7627c79db971 1474 xr[0][l + i] = opposite;
Gruenfrosch 0:7627c79db971 1475 xr[1][l + i] = left;
Gruenfrosch 0:7627c79db971 1476 }
Gruenfrosch 0:7627c79db971 1477 else
Gruenfrosch 0:7627c79db971 1478 xr[1][l + i] = opposite;
Gruenfrosch 0:7627c79db971 1479 }
Gruenfrosch 0:7627c79db971 1480 }
Gruenfrosch 0:7627c79db971 1481 }
Gruenfrosch 0:7627c79db971 1482 }
Gruenfrosch 0:7627c79db971 1483 else { /* !(header->flags & MAD_FLAG_LSF_EXT) */
Gruenfrosch 0:7627c79db971 1484 for (sfbi = l = 0; l < 576; ++sfbi, l += n) {
Gruenfrosch 0:7627c79db971 1485 n = sfbwidth[sfbi];
Gruenfrosch 0:7627c79db971 1486
Gruenfrosch 0:7627c79db971 1487 if (!(modes[sfbi] & I_STEREO))
Gruenfrosch 0:7627c79db971 1488 continue;
Gruenfrosch 0:7627c79db971 1489
Gruenfrosch 0:7627c79db971 1490 is_pos = right_ch->scalefac[sfbi];
Gruenfrosch 0:7627c79db971 1491
Gruenfrosch 0:7627c79db971 1492 if (is_pos >= 7) { /* illegal intensity position */
Gruenfrosch 0:7627c79db971 1493 modes[sfbi] &= ~I_STEREO;
Gruenfrosch 0:7627c79db971 1494 continue;
Gruenfrosch 0:7627c79db971 1495 }
Gruenfrosch 0:7627c79db971 1496
Gruenfrosch 0:7627c79db971 1497 for (i = 0; i < n; ++i) {
Gruenfrosch 0:7627c79db971 1498 register mad_fixed_t left;
Gruenfrosch 0:7627c79db971 1499
Gruenfrosch 0:7627c79db971 1500 left = xr[0][l + i];
Gruenfrosch 0:7627c79db971 1501
Gruenfrosch 0:7627c79db971 1502 xr[0][l + i] = mad_f_mul(left, is_table[ is_pos]);
Gruenfrosch 0:7627c79db971 1503 xr[1][l + i] = mad_f_mul(left, is_table[6 - is_pos]);
Gruenfrosch 0:7627c79db971 1504 }
Gruenfrosch 0:7627c79db971 1505 }
Gruenfrosch 0:7627c79db971 1506 }
Gruenfrosch 0:7627c79db971 1507 }
Gruenfrosch 0:7627c79db971 1508
Gruenfrosch 0:7627c79db971 1509 /* middle/side stereo */
Gruenfrosch 0:7627c79db971 1510
Gruenfrosch 0:7627c79db971 1511 if (header->mode_extension & MS_STEREO) {
Gruenfrosch 0:7627c79db971 1512 register mad_fixed_t invsqrt2;
Gruenfrosch 0:7627c79db971 1513
Gruenfrosch 0:7627c79db971 1514 header->flags |= MAD_FLAG_MS_STEREO;
Gruenfrosch 0:7627c79db971 1515
Gruenfrosch 0:7627c79db971 1516 invsqrt2 = root_table[3 + -2];
Gruenfrosch 0:7627c79db971 1517
Gruenfrosch 0:7627c79db971 1518 for (sfbi = l = 0; l < 576; ++sfbi, l += n) {
Gruenfrosch 0:7627c79db971 1519 n = sfbwidth[sfbi];
Gruenfrosch 0:7627c79db971 1520
Gruenfrosch 0:7627c79db971 1521 if (modes[sfbi] != MS_STEREO)
Gruenfrosch 0:7627c79db971 1522 continue;
Gruenfrosch 0:7627c79db971 1523
Gruenfrosch 0:7627c79db971 1524 for (i = 0; i < n; ++i) {
Gruenfrosch 0:7627c79db971 1525 register mad_fixed_t m, s;
Gruenfrosch 0:7627c79db971 1526
Gruenfrosch 0:7627c79db971 1527 m = xr[0][l + i];
Gruenfrosch 0:7627c79db971 1528 s = xr[1][l + i];
Gruenfrosch 0:7627c79db971 1529
Gruenfrosch 0:7627c79db971 1530 xr[0][l + i] = mad_f_mul(m + s, invsqrt2); /* l = (m + s) / sqrt(2) */
Gruenfrosch 0:7627c79db971 1531 xr[1][l + i] = mad_f_mul(m - s, invsqrt2); /* r = (m - s) / sqrt(2) */
Gruenfrosch 0:7627c79db971 1532 }
Gruenfrosch 0:7627c79db971 1533 }
Gruenfrosch 0:7627c79db971 1534 }
Gruenfrosch 0:7627c79db971 1535
Gruenfrosch 0:7627c79db971 1536 return MAD_ERROR_NONE;
Gruenfrosch 0:7627c79db971 1537 }
Gruenfrosch 0:7627c79db971 1538
Gruenfrosch 0:7627c79db971 1539 /*
Gruenfrosch 0:7627c79db971 1540 * NAME: III_aliasreduce()
Gruenfrosch 0:7627c79db971 1541 * DESCRIPTION: perform frequency line alias reduction
Gruenfrosch 0:7627c79db971 1542 */
Gruenfrosch 0:7627c79db971 1543 static
Gruenfrosch 0:7627c79db971 1544 void III_aliasreduce(mad_fixed_t xr[576], int lines)
Gruenfrosch 0:7627c79db971 1545 {
Gruenfrosch 0:7627c79db971 1546 mad_fixed_t const *bound;
Gruenfrosch 0:7627c79db971 1547 int i;
Gruenfrosch 0:7627c79db971 1548
Gruenfrosch 0:7627c79db971 1549 bound = &xr[lines];
Gruenfrosch 0:7627c79db971 1550 for (xr += 18; xr < bound; xr += 18) {
Gruenfrosch 0:7627c79db971 1551 for (i = 0; i < 8; ++i) {
Gruenfrosch 0:7627c79db971 1552 register mad_fixed_t a, b;
Gruenfrosch 0:7627c79db971 1553 register mad_fixed64hi_t hi;
Gruenfrosch 0:7627c79db971 1554 register mad_fixed64lo_t lo;
Gruenfrosch 0:7627c79db971 1555
Gruenfrosch 0:7627c79db971 1556 a = xr[-1 - i];
Gruenfrosch 0:7627c79db971 1557 b = xr[ i];
Gruenfrosch 0:7627c79db971 1558
Gruenfrosch 0:7627c79db971 1559 # if defined(ASO_ZEROCHECK)
Gruenfrosch 0:7627c79db971 1560 if (a | b) {
Gruenfrosch 0:7627c79db971 1561 # endif
Gruenfrosch 0:7627c79db971 1562 MAD_F_ML0(hi, lo, a, cs[i]);
Gruenfrosch 0:7627c79db971 1563 MAD_F_MLA(hi, lo, -b, ca[i]);
Gruenfrosch 0:7627c79db971 1564
Gruenfrosch 0:7627c79db971 1565 xr[-1 - i] = MAD_F_MLZ(hi, lo);
Gruenfrosch 0:7627c79db971 1566
Gruenfrosch 0:7627c79db971 1567 MAD_F_ML0(hi, lo, b, cs[i]);
Gruenfrosch 0:7627c79db971 1568 MAD_F_MLA(hi, lo, a, ca[i]);
Gruenfrosch 0:7627c79db971 1569
Gruenfrosch 0:7627c79db971 1570 xr[ i] = MAD_F_MLZ(hi, lo);
Gruenfrosch 0:7627c79db971 1571 # if defined(ASO_ZEROCHECK)
Gruenfrosch 0:7627c79db971 1572 }
Gruenfrosch 0:7627c79db971 1573 # endif
Gruenfrosch 0:7627c79db971 1574 }
Gruenfrosch 0:7627c79db971 1575 }
Gruenfrosch 0:7627c79db971 1576 }
Gruenfrosch 0:7627c79db971 1577
Gruenfrosch 0:7627c79db971 1578 # if defined(ASO_IMDCT)
Gruenfrosch 0:7627c79db971 1579 void III_imdct_l(mad_fixed_t const [18], mad_fixed_t [36], unsigned int);
Gruenfrosch 0:7627c79db971 1580 # else
Gruenfrosch 0:7627c79db971 1581 # if 1
Gruenfrosch 0:7627c79db971 1582 static
Gruenfrosch 0:7627c79db971 1583 void fastsdct(mad_fixed_t const x[9], mad_fixed_t y[18])
Gruenfrosch 0:7627c79db971 1584 {
Gruenfrosch 0:7627c79db971 1585 mad_fixed_t a0, a1, a2, a3, a4, a5, a6, a7, a8, a9, a10, a11, a12;
Gruenfrosch 0:7627c79db971 1586 mad_fixed_t a13, a14, a15, a16, a17, a18, a19, a20, a21, a22, a23, a24, a25;
Gruenfrosch 0:7627c79db971 1587 mad_fixed_t m0, m1, m2, m3, m4, m5, m6, m7;
Gruenfrosch 0:7627c79db971 1588
Gruenfrosch 0:7627c79db971 1589 enum {
Gruenfrosch 0:7627c79db971 1590 c0 = MAD_F(0x1f838b8d), /* 2 * cos( 1 * PI / 18) */
Gruenfrosch 0:7627c79db971 1591 c1 = MAD_F(0x1bb67ae8), /* 2 * cos( 3 * PI / 18) */
Gruenfrosch 0:7627c79db971 1592 c2 = MAD_F(0x18836fa3), /* 2 * cos( 4 * PI / 18) */
Gruenfrosch 0:7627c79db971 1593 c3 = MAD_F(0x1491b752), /* 2 * cos( 5 * PI / 18) */
Gruenfrosch 0:7627c79db971 1594 c4 = MAD_F(0x0af1d43a), /* 2 * cos( 7 * PI / 18) */
Gruenfrosch 0:7627c79db971 1595 c5 = MAD_F(0x058e86a0), /* 2 * cos( 8 * PI / 18) */
Gruenfrosch 0:7627c79db971 1596 c6 = -MAD_F(0x1e11f642) /* 2 * cos(16 * PI / 18) */
Gruenfrosch 0:7627c79db971 1597 };
Gruenfrosch 0:7627c79db971 1598
Gruenfrosch 0:7627c79db971 1599 a0 = x[3] + x[5];
Gruenfrosch 0:7627c79db971 1600 a1 = x[3] - x[5];
Gruenfrosch 0:7627c79db971 1601 a2 = x[6] + x[2];
Gruenfrosch 0:7627c79db971 1602 a3 = x[6] - x[2];
Gruenfrosch 0:7627c79db971 1603 a4 = x[1] + x[7];
Gruenfrosch 0:7627c79db971 1604 a5 = x[1] - x[7];
Gruenfrosch 0:7627c79db971 1605 a6 = x[8] + x[0];
Gruenfrosch 0:7627c79db971 1606 a7 = x[8] - x[0];
Gruenfrosch 0:7627c79db971 1607
Gruenfrosch 0:7627c79db971 1608 a8 = a0 + a2;
Gruenfrosch 0:7627c79db971 1609 a9 = a0 - a2;
Gruenfrosch 0:7627c79db971 1610 a10 = a0 - a6;
Gruenfrosch 0:7627c79db971 1611 a11 = a2 - a6;
Gruenfrosch 0:7627c79db971 1612 a12 = a8 + a6;
Gruenfrosch 0:7627c79db971 1613 a13 = a1 - a3;
Gruenfrosch 0:7627c79db971 1614 a14 = a13 + a7;
Gruenfrosch 0:7627c79db971 1615 a15 = a3 + a7;
Gruenfrosch 0:7627c79db971 1616 a16 = a1 - a7;
Gruenfrosch 0:7627c79db971 1617 a17 = a1 + a3;
Gruenfrosch 0:7627c79db971 1618
Gruenfrosch 0:7627c79db971 1619 m0 = mad_f_mul(a17, -c3);
Gruenfrosch 0:7627c79db971 1620 m1 = mad_f_mul(a16, -c0);
Gruenfrosch 0:7627c79db971 1621 m2 = mad_f_mul(a15, -c4);
Gruenfrosch 0:7627c79db971 1622 m3 = mad_f_mul(a14, -c1);
Gruenfrosch 0:7627c79db971 1623 m4 = mad_f_mul(a5, -c1);
Gruenfrosch 0:7627c79db971 1624 m5 = mad_f_mul(a11, -c6);
Gruenfrosch 0:7627c79db971 1625 m6 = mad_f_mul(a10, -c5);
Gruenfrosch 0:7627c79db971 1626 m7 = mad_f_mul(a9, -c2);
Gruenfrosch 0:7627c79db971 1627
Gruenfrosch 0:7627c79db971 1628 a18 = x[4] + a4;
Gruenfrosch 0:7627c79db971 1629 a19 = 2 * x[4] - a4;
Gruenfrosch 0:7627c79db971 1630 a20 = a19 + m5;
Gruenfrosch 0:7627c79db971 1631 a21 = a19 - m5;
Gruenfrosch 0:7627c79db971 1632 a22 = a19 + m6;
Gruenfrosch 0:7627c79db971 1633 a23 = m4 + m2;
Gruenfrosch 0:7627c79db971 1634 a24 = m4 - m2;
Gruenfrosch 0:7627c79db971 1635 a25 = m4 + m1;
Gruenfrosch 0:7627c79db971 1636
Gruenfrosch 0:7627c79db971 1637 /* output to every other slot for convenience */
Gruenfrosch 0:7627c79db971 1638
Gruenfrosch 0:7627c79db971 1639 y[ 0] = a18 + a12;
Gruenfrosch 0:7627c79db971 1640 y[ 2] = m0 - a25;
Gruenfrosch 0:7627c79db971 1641 y[ 4] = m7 - a20;
Gruenfrosch 0:7627c79db971 1642 y[ 6] = m3;
Gruenfrosch 0:7627c79db971 1643 y[ 8] = a21 - m6;
Gruenfrosch 0:7627c79db971 1644 y[10] = a24 - m1;
Gruenfrosch 0:7627c79db971 1645 y[12] = a12 - 2 * a18;
Gruenfrosch 0:7627c79db971 1646 y[14] = a23 + m0;
Gruenfrosch 0:7627c79db971 1647 y[16] = a22 + m7;
Gruenfrosch 0:7627c79db971 1648 }
Gruenfrosch 0:7627c79db971 1649
Gruenfrosch 0:7627c79db971 1650 static inline
Gruenfrosch 0:7627c79db971 1651 void sdctII(mad_fixed_t const x[18], mad_fixed_t X[18])
Gruenfrosch 0:7627c79db971 1652 {
Gruenfrosch 0:7627c79db971 1653 mad_fixed_t tmp[9];
Gruenfrosch 0:7627c79db971 1654 int i;
Gruenfrosch 0:7627c79db971 1655
Gruenfrosch 0:7627c79db971 1656 /* scale[i] = 2 * cos(PI * (2 * i + 1) / (2 * 18)) */
Gruenfrosch 0:7627c79db971 1657 static mad_fixed_t const scale[9] = {
Gruenfrosch 0:7627c79db971 1658 MAD_F(0x1fe0d3b4), MAD_F(0x1ee8dd47), MAD_F(0x1d007930),
Gruenfrosch 0:7627c79db971 1659 MAD_F(0x1a367e59), MAD_F(0x16a09e66), MAD_F(0x125abcf8),
Gruenfrosch 0:7627c79db971 1660 MAD_F(0x0d8616bc), MAD_F(0x08483ee1), MAD_F(0x02c9fad7)
Gruenfrosch 0:7627c79db971 1661 };
Gruenfrosch 0:7627c79db971 1662
Gruenfrosch 0:7627c79db971 1663 /* divide the 18-point SDCT-II into two 9-point SDCT-IIs */
Gruenfrosch 0:7627c79db971 1664
Gruenfrosch 0:7627c79db971 1665 /* even input butterfly */
Gruenfrosch 0:7627c79db971 1666
Gruenfrosch 0:7627c79db971 1667 for (i = 0; i < 9; i += 3) {
Gruenfrosch 0:7627c79db971 1668 tmp[i + 0] = x[i + 0] + x[18 - (i + 0) - 1];
Gruenfrosch 0:7627c79db971 1669 tmp[i + 1] = x[i + 1] + x[18 - (i + 1) - 1];
Gruenfrosch 0:7627c79db971 1670 tmp[i + 2] = x[i + 2] + x[18 - (i + 2) - 1];
Gruenfrosch 0:7627c79db971 1671 }
Gruenfrosch 0:7627c79db971 1672
Gruenfrosch 0:7627c79db971 1673 fastsdct(tmp, &X[0]);
Gruenfrosch 0:7627c79db971 1674
Gruenfrosch 0:7627c79db971 1675 /* odd input butterfly and scaling */
Gruenfrosch 0:7627c79db971 1676
Gruenfrosch 0:7627c79db971 1677 for (i = 0; i < 9; i += 3) {
Gruenfrosch 0:7627c79db971 1678 tmp[i + 0] = mad_f_mul(x[i + 0] - x[18 - (i + 0) - 1], scale[i + 0]);
Gruenfrosch 0:7627c79db971 1679 tmp[i + 1] = mad_f_mul(x[i + 1] - x[18 - (i + 1) - 1], scale[i + 1]);
Gruenfrosch 0:7627c79db971 1680 tmp[i + 2] = mad_f_mul(x[i + 2] - x[18 - (i + 2) - 1], scale[i + 2]);
Gruenfrosch 0:7627c79db971 1681 }
Gruenfrosch 0:7627c79db971 1682
Gruenfrosch 0:7627c79db971 1683 fastsdct(tmp, &X[1]);
Gruenfrosch 0:7627c79db971 1684
Gruenfrosch 0:7627c79db971 1685 /* output accumulation */
Gruenfrosch 0:7627c79db971 1686
Gruenfrosch 0:7627c79db971 1687 for (i = 3; i < 18; i += 8) {
Gruenfrosch 0:7627c79db971 1688 X[i + 0] -= X[(i + 0) - 2];
Gruenfrosch 0:7627c79db971 1689 X[i + 2] -= X[(i + 2) - 2];
Gruenfrosch 0:7627c79db971 1690 X[i + 4] -= X[(i + 4) - 2];
Gruenfrosch 0:7627c79db971 1691 X[i + 6] -= X[(i + 6) - 2];
Gruenfrosch 0:7627c79db971 1692 }
Gruenfrosch 0:7627c79db971 1693 }
Gruenfrosch 0:7627c79db971 1694
Gruenfrosch 0:7627c79db971 1695 static inline
Gruenfrosch 0:7627c79db971 1696 void dctIV(mad_fixed_t const y[18], mad_fixed_t X[18])
Gruenfrosch 0:7627c79db971 1697 {
Gruenfrosch 0:7627c79db971 1698 mad_fixed_t tmp[18];
Gruenfrosch 0:7627c79db971 1699 int i;
Gruenfrosch 0:7627c79db971 1700
Gruenfrosch 0:7627c79db971 1701 /* scale[i] = 2 * cos(PI * (2 * i + 1) / (4 * 18)) */
Gruenfrosch 0:7627c79db971 1702 static mad_fixed_t const scale[18] = {
Gruenfrosch 0:7627c79db971 1703 MAD_F(0x1ff833fa), MAD_F(0x1fb9ea93), MAD_F(0x1f3dd120),
Gruenfrosch 0:7627c79db971 1704 MAD_F(0x1e84d969), MAD_F(0x1d906bcf), MAD_F(0x1c62648b),
Gruenfrosch 0:7627c79db971 1705 MAD_F(0x1afd100f), MAD_F(0x1963268b), MAD_F(0x1797c6a4),
Gruenfrosch 0:7627c79db971 1706 MAD_F(0x159e6f5b), MAD_F(0x137af940), MAD_F(0x11318ef3),
Gruenfrosch 0:7627c79db971 1707 MAD_F(0x0ec6a507), MAD_F(0x0c3ef153), MAD_F(0x099f61c5),
Gruenfrosch 0:7627c79db971 1708 MAD_F(0x06ed12c5), MAD_F(0x042d4544), MAD_F(0x0165547c)
Gruenfrosch 0:7627c79db971 1709 };
Gruenfrosch 0:7627c79db971 1710
Gruenfrosch 0:7627c79db971 1711 /* scaling */
Gruenfrosch 0:7627c79db971 1712
Gruenfrosch 0:7627c79db971 1713 for (i = 0; i < 18; i += 3) {
Gruenfrosch 0:7627c79db971 1714 tmp[i + 0] = mad_f_mul(y[i + 0], scale[i + 0]);
Gruenfrosch 0:7627c79db971 1715 tmp[i + 1] = mad_f_mul(y[i + 1], scale[i + 1]);
Gruenfrosch 0:7627c79db971 1716 tmp[i + 2] = mad_f_mul(y[i + 2], scale[i + 2]);
Gruenfrosch 0:7627c79db971 1717 }
Gruenfrosch 0:7627c79db971 1718
Gruenfrosch 0:7627c79db971 1719 /* SDCT-II */
Gruenfrosch 0:7627c79db971 1720
Gruenfrosch 0:7627c79db971 1721 sdctII(tmp, X);
Gruenfrosch 0:7627c79db971 1722
Gruenfrosch 0:7627c79db971 1723 /* scale reduction and output accumulation */
Gruenfrosch 0:7627c79db971 1724
Gruenfrosch 0:7627c79db971 1725 X[0] /= 2;
Gruenfrosch 0:7627c79db971 1726 for (i = 1; i < 17; i += 4) {
Gruenfrosch 0:7627c79db971 1727 X[i + 0] = X[i + 0] / 2 - X[(i + 0) - 1];
Gruenfrosch 0:7627c79db971 1728 X[i + 1] = X[i + 1] / 2 - X[(i + 1) - 1];
Gruenfrosch 0:7627c79db971 1729 X[i + 2] = X[i + 2] / 2 - X[(i + 2) - 1];
Gruenfrosch 0:7627c79db971 1730 X[i + 3] = X[i + 3] / 2 - X[(i + 3) - 1];
Gruenfrosch 0:7627c79db971 1731 }
Gruenfrosch 0:7627c79db971 1732 X[17] = X[17] / 2 - X[16];
Gruenfrosch 0:7627c79db971 1733 }
Gruenfrosch 0:7627c79db971 1734
Gruenfrosch 0:7627c79db971 1735 /*
Gruenfrosch 0:7627c79db971 1736 * NAME: imdct36
Gruenfrosch 0:7627c79db971 1737 * DESCRIPTION: perform X[18]->x[36] IMDCT using Szu-Wei Lee's fast algorithm
Gruenfrosch 0:7627c79db971 1738 */
Gruenfrosch 0:7627c79db971 1739 static inline
Gruenfrosch 0:7627c79db971 1740 void imdct36(mad_fixed_t const x[18], mad_fixed_t y[36])
Gruenfrosch 0:7627c79db971 1741 {
Gruenfrosch 0:7627c79db971 1742 mad_fixed_t tmp[18];
Gruenfrosch 0:7627c79db971 1743 int i;
Gruenfrosch 0:7627c79db971 1744
Gruenfrosch 0:7627c79db971 1745 /* DCT-IV */
Gruenfrosch 0:7627c79db971 1746
Gruenfrosch 0:7627c79db971 1747 dctIV(x, tmp);
Gruenfrosch 0:7627c79db971 1748
Gruenfrosch 0:7627c79db971 1749 /* convert 18-point DCT-IV to 36-point IMDCT */
Gruenfrosch 0:7627c79db971 1750
Gruenfrosch 0:7627c79db971 1751 for (i = 0; i < 9; i += 3) {
Gruenfrosch 0:7627c79db971 1752 y[i + 0] = tmp[9 + (i + 0)];
Gruenfrosch 0:7627c79db971 1753 y[i + 1] = tmp[9 + (i + 1)];
Gruenfrosch 0:7627c79db971 1754 y[i + 2] = tmp[9 + (i + 2)];
Gruenfrosch 0:7627c79db971 1755 }
Gruenfrosch 0:7627c79db971 1756 for (i = 9; i < 27; i += 3) {
Gruenfrosch 0:7627c79db971 1757 y[i + 0] = -tmp[36 - (9 + (i + 0)) - 1];
Gruenfrosch 0:7627c79db971 1758 y[i + 1] = -tmp[36 - (9 + (i + 1)) - 1];
Gruenfrosch 0:7627c79db971 1759 y[i + 2] = -tmp[36 - (9 + (i + 2)) - 1];
Gruenfrosch 0:7627c79db971 1760 }
Gruenfrosch 0:7627c79db971 1761 for (i = 27; i < 36; i += 3) {
Gruenfrosch 0:7627c79db971 1762 y[i + 0] = -tmp[(i + 0) - 27];
Gruenfrosch 0:7627c79db971 1763 y[i + 1] = -tmp[(i + 1) - 27];
Gruenfrosch 0:7627c79db971 1764 y[i + 2] = -tmp[(i + 2) - 27];
Gruenfrosch 0:7627c79db971 1765 }
Gruenfrosch 0:7627c79db971 1766 }
Gruenfrosch 0:7627c79db971 1767 # else
Gruenfrosch 0:7627c79db971 1768 /*
Gruenfrosch 0:7627c79db971 1769 * NAME: imdct36
Gruenfrosch 0:7627c79db971 1770 * DESCRIPTION: perform X[18]->x[36] IMDCT
Gruenfrosch 0:7627c79db971 1771 */
Gruenfrosch 0:7627c79db971 1772 static inline
Gruenfrosch 0:7627c79db971 1773 void imdct36(mad_fixed_t const X[18], mad_fixed_t x[36])
Gruenfrosch 0:7627c79db971 1774 {
Gruenfrosch 0:7627c79db971 1775 mad_fixed_t t0, t1, t2, t3, t4, t5, t6, t7;
Gruenfrosch 0:7627c79db971 1776 mad_fixed_t t8, t9, t10, t11, t12, t13, t14, t15;
Gruenfrosch 0:7627c79db971 1777 register mad_fixed64hi_t hi;
Gruenfrosch 0:7627c79db971 1778 register mad_fixed64lo_t lo;
Gruenfrosch 0:7627c79db971 1779
Gruenfrosch 0:7627c79db971 1780 MAD_F_ML0(hi, lo, X[4], MAD_F(0x0ec835e8));
Gruenfrosch 0:7627c79db971 1781 MAD_F_MLA(hi, lo, X[13], MAD_F(0x061f78aa));
Gruenfrosch 0:7627c79db971 1782
Gruenfrosch 0:7627c79db971 1783 t6 = MAD_F_MLZ(hi, lo);
Gruenfrosch 0:7627c79db971 1784
Gruenfrosch 0:7627c79db971 1785 MAD_F_MLA(hi, lo, (t14 = X[1] - X[10]), -MAD_F(0x061f78aa));
Gruenfrosch 0:7627c79db971 1786 MAD_F_MLA(hi, lo, (t15 = X[7] + X[16]), -MAD_F(0x0ec835e8));
Gruenfrosch 0:7627c79db971 1787
Gruenfrosch 0:7627c79db971 1788 t0 = MAD_F_MLZ(hi, lo);
Gruenfrosch 0:7627c79db971 1789
Gruenfrosch 0:7627c79db971 1790 MAD_F_MLA(hi, lo, (t8 = X[0] - X[11] - X[12]), MAD_F(0x0216a2a2));
Gruenfrosch 0:7627c79db971 1791 MAD_F_MLA(hi, lo, (t9 = X[2] - X[9] - X[14]), MAD_F(0x09bd7ca0));
Gruenfrosch 0:7627c79db971 1792 MAD_F_MLA(hi, lo, (t10 = X[3] - X[8] - X[15]), -MAD_F(0x0cb19346));
Gruenfrosch 0:7627c79db971 1793 MAD_F_MLA(hi, lo, (t11 = X[5] - X[6] - X[17]), -MAD_F(0x0fdcf549));
Gruenfrosch 0:7627c79db971 1794
Gruenfrosch 0:7627c79db971 1795 x[7] = MAD_F_MLZ(hi, lo);
Gruenfrosch 0:7627c79db971 1796 x[10] = -x[7];
Gruenfrosch 0:7627c79db971 1797
Gruenfrosch 0:7627c79db971 1798 MAD_F_ML0(hi, lo, t8, -MAD_F(0x0cb19346));
Gruenfrosch 0:7627c79db971 1799 MAD_F_MLA(hi, lo, t9, MAD_F(0x0fdcf549));
Gruenfrosch 0:7627c79db971 1800 MAD_F_MLA(hi, lo, t10, MAD_F(0x0216a2a2));
Gruenfrosch 0:7627c79db971 1801 MAD_F_MLA(hi, lo, t11, -MAD_F(0x09bd7ca0));
Gruenfrosch 0:7627c79db971 1802
Gruenfrosch 0:7627c79db971 1803 x[19] = x[34] = MAD_F_MLZ(hi, lo) - t0;
Gruenfrosch 0:7627c79db971 1804
Gruenfrosch 0:7627c79db971 1805 t12 = X[0] - X[3] + X[8] - X[11] - X[12] + X[15];
Gruenfrosch 0:7627c79db971 1806 t13 = X[2] + X[5] - X[6] - X[9] - X[14] - X[17];
Gruenfrosch 0:7627c79db971 1807
Gruenfrosch 0:7627c79db971 1808 MAD_F_ML0(hi, lo, t12, -MAD_F(0x0ec835e8));
Gruenfrosch 0:7627c79db971 1809 MAD_F_MLA(hi, lo, t13, MAD_F(0x061f78aa));
Gruenfrosch 0:7627c79db971 1810
Gruenfrosch 0:7627c79db971 1811 x[22] = x[31] = MAD_F_MLZ(hi, lo) + t0;
Gruenfrosch 0:7627c79db971 1812
Gruenfrosch 0:7627c79db971 1813 MAD_F_ML0(hi, lo, X[1], -MAD_F(0x09bd7ca0));
Gruenfrosch 0:7627c79db971 1814 MAD_F_MLA(hi, lo, X[7], MAD_F(0x0216a2a2));
Gruenfrosch 0:7627c79db971 1815 MAD_F_MLA(hi, lo, X[10], -MAD_F(0x0fdcf549));
Gruenfrosch 0:7627c79db971 1816 MAD_F_MLA(hi, lo, X[16], MAD_F(0x0cb19346));
Gruenfrosch 0:7627c79db971 1817
Gruenfrosch 0:7627c79db971 1818 t1 = MAD_F_MLZ(hi, lo) + t6;
Gruenfrosch 0:7627c79db971 1819
Gruenfrosch 0:7627c79db971 1820 MAD_F_ML0(hi, lo, X[0], MAD_F(0x03768962));
Gruenfrosch 0:7627c79db971 1821 MAD_F_MLA(hi, lo, X[2], MAD_F(0x0e313245));
Gruenfrosch 0:7627c79db971 1822 MAD_F_MLA(hi, lo, X[3], -MAD_F(0x0ffc19fd));
Gruenfrosch 0:7627c79db971 1823 MAD_F_MLA(hi, lo, X[5], -MAD_F(0x0acf37ad));
Gruenfrosch 0:7627c79db971 1824 MAD_F_MLA(hi, lo, X[6], MAD_F(0x04cfb0e2));
Gruenfrosch 0:7627c79db971 1825 MAD_F_MLA(hi, lo, X[8], -MAD_F(0x0898c779));
Gruenfrosch 0:7627c79db971 1826 MAD_F_MLA(hi, lo, X[9], MAD_F(0x0d7e8807));
Gruenfrosch 0:7627c79db971 1827 MAD_F_MLA(hi, lo, X[11], MAD_F(0x0f426cb5));
Gruenfrosch 0:7627c79db971 1828 MAD_F_MLA(hi, lo, X[12], -MAD_F(0x0bcbe352));
Gruenfrosch 0:7627c79db971 1829 MAD_F_MLA(hi, lo, X[14], MAD_F(0x00b2aa3e));
Gruenfrosch 0:7627c79db971 1830 MAD_F_MLA(hi, lo, X[15], -MAD_F(0x07635284));
Gruenfrosch 0:7627c79db971 1831 MAD_F_MLA(hi, lo, X[17], -MAD_F(0x0f9ee890));
Gruenfrosch 0:7627c79db971 1832
Gruenfrosch 0:7627c79db971 1833 x[6] = MAD_F_MLZ(hi, lo) + t1;
Gruenfrosch 0:7627c79db971 1834 x[11] = -x[6];
Gruenfrosch 0:7627c79db971 1835
Gruenfrosch 0:7627c79db971 1836 MAD_F_ML0(hi, lo, X[0], -MAD_F(0x0f426cb5));
Gruenfrosch 0:7627c79db971 1837 MAD_F_MLA(hi, lo, X[2], -MAD_F(0x00b2aa3e));
Gruenfrosch 0:7627c79db971 1838 MAD_F_MLA(hi, lo, X[3], MAD_F(0x0898c779));
Gruenfrosch 0:7627c79db971 1839 MAD_F_MLA(hi, lo, X[5], MAD_F(0x0f9ee890));
Gruenfrosch 0:7627c79db971 1840 MAD_F_MLA(hi, lo, X[6], MAD_F(0x0acf37ad));
Gruenfrosch 0:7627c79db971 1841 MAD_F_MLA(hi, lo, X[8], -MAD_F(0x07635284));
Gruenfrosch 0:7627c79db971 1842 MAD_F_MLA(hi, lo, X[9], -MAD_F(0x0e313245));
Gruenfrosch 0:7627c79db971 1843 MAD_F_MLA(hi, lo, X[11], -MAD_F(0x0bcbe352));
Gruenfrosch 0:7627c79db971 1844 MAD_F_MLA(hi, lo, X[12], -MAD_F(0x03768962));
Gruenfrosch 0:7627c79db971 1845 MAD_F_MLA(hi, lo, X[14], MAD_F(0x0d7e8807));
Gruenfrosch 0:7627c79db971 1846 MAD_F_MLA(hi, lo, X[15], MAD_F(0x0ffc19fd));
Gruenfrosch 0:7627c79db971 1847 MAD_F_MLA(hi, lo, X[17], MAD_F(0x04cfb0e2));
Gruenfrosch 0:7627c79db971 1848
Gruenfrosch 0:7627c79db971 1849 x[23] = x[30] = MAD_F_MLZ(hi, lo) + t1;
Gruenfrosch 0:7627c79db971 1850
Gruenfrosch 0:7627c79db971 1851 MAD_F_ML0(hi, lo, X[0], -MAD_F(0x0bcbe352));
Gruenfrosch 0:7627c79db971 1852 MAD_F_MLA(hi, lo, X[2], MAD_F(0x0d7e8807));
Gruenfrosch 0:7627c79db971 1853 MAD_F_MLA(hi, lo, X[3], -MAD_F(0x07635284));
Gruenfrosch 0:7627c79db971 1854 MAD_F_MLA(hi, lo, X[5], MAD_F(0x04cfb0e2));
Gruenfrosch 0:7627c79db971 1855 MAD_F_MLA(hi, lo, X[6], MAD_F(0x0f9ee890));
Gruenfrosch 0:7627c79db971 1856 MAD_F_MLA(hi, lo, X[8], -MAD_F(0x0ffc19fd));
Gruenfrosch 0:7627c79db971 1857 MAD_F_MLA(hi, lo, X[9], -MAD_F(0x00b2aa3e));
Gruenfrosch 0:7627c79db971 1858 MAD_F_MLA(hi, lo, X[11], MAD_F(0x03768962));
Gruenfrosch 0:7627c79db971 1859 MAD_F_MLA(hi, lo, X[12], -MAD_F(0x0f426cb5));
Gruenfrosch 0:7627c79db971 1860 MAD_F_MLA(hi, lo, X[14], MAD_F(0x0e313245));
Gruenfrosch 0:7627c79db971 1861 MAD_F_MLA(hi, lo, X[15], MAD_F(0x0898c779));
Gruenfrosch 0:7627c79db971 1862 MAD_F_MLA(hi, lo, X[17], -MAD_F(0x0acf37ad));
Gruenfrosch 0:7627c79db971 1863
Gruenfrosch 0:7627c79db971 1864 x[18] = x[35] = MAD_F_MLZ(hi, lo) - t1;
Gruenfrosch 0:7627c79db971 1865
Gruenfrosch 0:7627c79db971 1866 MAD_F_ML0(hi, lo, X[4], MAD_F(0x061f78aa));
Gruenfrosch 0:7627c79db971 1867 MAD_F_MLA(hi, lo, X[13], -MAD_F(0x0ec835e8));
Gruenfrosch 0:7627c79db971 1868
Gruenfrosch 0:7627c79db971 1869 t7 = MAD_F_MLZ(hi, lo);
Gruenfrosch 0:7627c79db971 1870
Gruenfrosch 0:7627c79db971 1871 MAD_F_MLA(hi, lo, X[1], -MAD_F(0x0cb19346));
Gruenfrosch 0:7627c79db971 1872 MAD_F_MLA(hi, lo, X[7], MAD_F(0x0fdcf549));
Gruenfrosch 0:7627c79db971 1873 MAD_F_MLA(hi, lo, X[10], MAD_F(0x0216a2a2));
Gruenfrosch 0:7627c79db971 1874 MAD_F_MLA(hi, lo, X[16], -MAD_F(0x09bd7ca0));
Gruenfrosch 0:7627c79db971 1875
Gruenfrosch 0:7627c79db971 1876 t2 = MAD_F_MLZ(hi, lo);
Gruenfrosch 0:7627c79db971 1877
Gruenfrosch 0:7627c79db971 1878 MAD_F_MLA(hi, lo, X[0], MAD_F(0x04cfb0e2));
Gruenfrosch 0:7627c79db971 1879 MAD_F_MLA(hi, lo, X[2], MAD_F(0x0ffc19fd));
Gruenfrosch 0:7627c79db971 1880 MAD_F_MLA(hi, lo, X[3], -MAD_F(0x0d7e8807));
Gruenfrosch 0:7627c79db971 1881 MAD_F_MLA(hi, lo, X[5], MAD_F(0x03768962));
Gruenfrosch 0:7627c79db971 1882 MAD_F_MLA(hi, lo, X[6], -MAD_F(0x0bcbe352));
Gruenfrosch 0:7627c79db971 1883 MAD_F_MLA(hi, lo, X[8], -MAD_F(0x0e313245));
Gruenfrosch 0:7627c79db971 1884 MAD_F_MLA(hi, lo, X[9], MAD_F(0x07635284));
Gruenfrosch 0:7627c79db971 1885 MAD_F_MLA(hi, lo, X[11], -MAD_F(0x0acf37ad));
Gruenfrosch 0:7627c79db971 1886 MAD_F_MLA(hi, lo, X[12], MAD_F(0x0f9ee890));
Gruenfrosch 0:7627c79db971 1887 MAD_F_MLA(hi, lo, X[14], MAD_F(0x0898c779));
Gruenfrosch 0:7627c79db971 1888 MAD_F_MLA(hi, lo, X[15], MAD_F(0x00b2aa3e));
Gruenfrosch 0:7627c79db971 1889 MAD_F_MLA(hi, lo, X[17], MAD_F(0x0f426cb5));
Gruenfrosch 0:7627c79db971 1890
Gruenfrosch 0:7627c79db971 1891 x[5] = MAD_F_MLZ(hi, lo);
Gruenfrosch 0:7627c79db971 1892 x[12] = -x[5];
Gruenfrosch 0:7627c79db971 1893
Gruenfrosch 0:7627c79db971 1894 MAD_F_ML0(hi, lo, X[0], MAD_F(0x0acf37ad));
Gruenfrosch 0:7627c79db971 1895 MAD_F_MLA(hi, lo, X[2], -MAD_F(0x0898c779));
Gruenfrosch 0:7627c79db971 1896 MAD_F_MLA(hi, lo, X[3], MAD_F(0x0e313245));
Gruenfrosch 0:7627c79db971 1897 MAD_F_MLA(hi, lo, X[5], -MAD_F(0x0f426cb5));
Gruenfrosch 0:7627c79db971 1898 MAD_F_MLA(hi, lo, X[6], -MAD_F(0x03768962));
Gruenfrosch 0:7627c79db971 1899 MAD_F_MLA(hi, lo, X[8], MAD_F(0x00b2aa3e));
Gruenfrosch 0:7627c79db971 1900 MAD_F_MLA(hi, lo, X[9], -MAD_F(0x0ffc19fd));
Gruenfrosch 0:7627c79db971 1901 MAD_F_MLA(hi, lo, X[11], MAD_F(0x0f9ee890));
Gruenfrosch 0:7627c79db971 1902 MAD_F_MLA(hi, lo, X[12], -MAD_F(0x04cfb0e2));
Gruenfrosch 0:7627c79db971 1903 MAD_F_MLA(hi, lo, X[14], MAD_F(0x07635284));
Gruenfrosch 0:7627c79db971 1904 MAD_F_MLA(hi, lo, X[15], MAD_F(0x0d7e8807));
Gruenfrosch 0:7627c79db971 1905 MAD_F_MLA(hi, lo, X[17], -MAD_F(0x0bcbe352));
Gruenfrosch 0:7627c79db971 1906
Gruenfrosch 0:7627c79db971 1907 x[0] = MAD_F_MLZ(hi, lo) + t2;
Gruenfrosch 0:7627c79db971 1908 x[17] = -x[0];
Gruenfrosch 0:7627c79db971 1909
Gruenfrosch 0:7627c79db971 1910 MAD_F_ML0(hi, lo, X[0], -MAD_F(0x0f9ee890));
Gruenfrosch 0:7627c79db971 1911 MAD_F_MLA(hi, lo, X[2], -MAD_F(0x07635284));
Gruenfrosch 0:7627c79db971 1912 MAD_F_MLA(hi, lo, X[3], -MAD_F(0x00b2aa3e));
Gruenfrosch 0:7627c79db971 1913 MAD_F_MLA(hi, lo, X[5], MAD_F(0x0bcbe352));
Gruenfrosch 0:7627c79db971 1914 MAD_F_MLA(hi, lo, X[6], MAD_F(0x0f426cb5));
Gruenfrosch 0:7627c79db971 1915 MAD_F_MLA(hi, lo, X[8], MAD_F(0x0d7e8807));
Gruenfrosch 0:7627c79db971 1916 MAD_F_MLA(hi, lo, X[9], MAD_F(0x0898c779));
Gruenfrosch 0:7627c79db971 1917 MAD_F_MLA(hi, lo, X[11], -MAD_F(0x04cfb0e2));
Gruenfrosch 0:7627c79db971 1918 MAD_F_MLA(hi, lo, X[12], -MAD_F(0x0acf37ad));
Gruenfrosch 0:7627c79db971 1919 MAD_F_MLA(hi, lo, X[14], -MAD_F(0x0ffc19fd));
Gruenfrosch 0:7627c79db971 1920 MAD_F_MLA(hi, lo, X[15], -MAD_F(0x0e313245));
Gruenfrosch 0:7627c79db971 1921 MAD_F_MLA(hi, lo, X[17], -MAD_F(0x03768962));
Gruenfrosch 0:7627c79db971 1922
Gruenfrosch 0:7627c79db971 1923 x[24] = x[29] = MAD_F_MLZ(hi, lo) + t2;
Gruenfrosch 0:7627c79db971 1924
Gruenfrosch 0:7627c79db971 1925 MAD_F_ML0(hi, lo, X[1], -MAD_F(0x0216a2a2));
Gruenfrosch 0:7627c79db971 1926 MAD_F_MLA(hi, lo, X[7], -MAD_F(0x09bd7ca0));
Gruenfrosch 0:7627c79db971 1927 MAD_F_MLA(hi, lo, X[10], MAD_F(0x0cb19346));
Gruenfrosch 0:7627c79db971 1928 MAD_F_MLA(hi, lo, X[16], MAD_F(0x0fdcf549));
Gruenfrosch 0:7627c79db971 1929
Gruenfrosch 0:7627c79db971 1930 t3 = MAD_F_MLZ(hi, lo) + t7;
Gruenfrosch 0:7627c79db971 1931
Gruenfrosch 0:7627c79db971 1932 MAD_F_ML0(hi, lo, X[0], MAD_F(0x00b2aa3e));
Gruenfrosch 0:7627c79db971 1933 MAD_F_MLA(hi, lo, X[2], MAD_F(0x03768962));
Gruenfrosch 0:7627c79db971 1934 MAD_F_MLA(hi, lo, X[3], -MAD_F(0x04cfb0e2));
Gruenfrosch 0:7627c79db971 1935 MAD_F_MLA(hi, lo, X[5], -MAD_F(0x07635284));
Gruenfrosch 0:7627c79db971 1936 MAD_F_MLA(hi, lo, X[6], MAD_F(0x0898c779));
Gruenfrosch 0:7627c79db971 1937 MAD_F_MLA(hi, lo, X[8], MAD_F(0x0acf37ad));
Gruenfrosch 0:7627c79db971 1938 MAD_F_MLA(hi, lo, X[9], -MAD_F(0x0bcbe352));
Gruenfrosch 0:7627c79db971 1939 MAD_F_MLA(hi, lo, X[11], -MAD_F(0x0d7e8807));
Gruenfrosch 0:7627c79db971 1940 MAD_F_MLA(hi, lo, X[12], MAD_F(0x0e313245));
Gruenfrosch 0:7627c79db971 1941 MAD_F_MLA(hi, lo, X[14], MAD_F(0x0f426cb5));
Gruenfrosch 0:7627c79db971 1942 MAD_F_MLA(hi, lo, X[15], -MAD_F(0x0f9ee890));
Gruenfrosch 0:7627c79db971 1943 MAD_F_MLA(hi, lo, X[17], -MAD_F(0x0ffc19fd));
Gruenfrosch 0:7627c79db971 1944
Gruenfrosch 0:7627c79db971 1945 x[8] = MAD_F_MLZ(hi, lo) + t3;
Gruenfrosch 0:7627c79db971 1946 x[9] = -x[8];
Gruenfrosch 0:7627c79db971 1947
Gruenfrosch 0:7627c79db971 1948 MAD_F_ML0(hi, lo, X[0], -MAD_F(0x0e313245));
Gruenfrosch 0:7627c79db971 1949 MAD_F_MLA(hi, lo, X[2], MAD_F(0x0bcbe352));
Gruenfrosch 0:7627c79db971 1950 MAD_F_MLA(hi, lo, X[3], MAD_F(0x0f9ee890));
Gruenfrosch 0:7627c79db971 1951 MAD_F_MLA(hi, lo, X[5], -MAD_F(0x0898c779));
Gruenfrosch 0:7627c79db971 1952 MAD_F_MLA(hi, lo, X[6], -MAD_F(0x0ffc19fd));
Gruenfrosch 0:7627c79db971 1953 MAD_F_MLA(hi, lo, X[8], MAD_F(0x04cfb0e2));
Gruenfrosch 0:7627c79db971 1954 MAD_F_MLA(hi, lo, X[9], MAD_F(0x0f426cb5));
Gruenfrosch 0:7627c79db971 1955 MAD_F_MLA(hi, lo, X[11], -MAD_F(0x00b2aa3e));
Gruenfrosch 0:7627c79db971 1956 MAD_F_MLA(hi, lo, X[12], -MAD_F(0x0d7e8807));
Gruenfrosch 0:7627c79db971 1957 MAD_F_MLA(hi, lo, X[14], -MAD_F(0x03768962));
Gruenfrosch 0:7627c79db971 1958 MAD_F_MLA(hi, lo, X[15], MAD_F(0x0acf37ad));
Gruenfrosch 0:7627c79db971 1959 MAD_F_MLA(hi, lo, X[17], MAD_F(0x07635284));
Gruenfrosch 0:7627c79db971 1960
Gruenfrosch 0:7627c79db971 1961 x[21] = x[32] = MAD_F_MLZ(hi, lo) + t3;
Gruenfrosch 0:7627c79db971 1962
Gruenfrosch 0:7627c79db971 1963 MAD_F_ML0(hi, lo, X[0], -MAD_F(0x0d7e8807));
Gruenfrosch 0:7627c79db971 1964 MAD_F_MLA(hi, lo, X[2], MAD_F(0x0f426cb5));
Gruenfrosch 0:7627c79db971 1965 MAD_F_MLA(hi, lo, X[3], MAD_F(0x0acf37ad));
Gruenfrosch 0:7627c79db971 1966 MAD_F_MLA(hi, lo, X[5], -MAD_F(0x0ffc19fd));
Gruenfrosch 0:7627c79db971 1967 MAD_F_MLA(hi, lo, X[6], -MAD_F(0x07635284));
Gruenfrosch 0:7627c79db971 1968 MAD_F_MLA(hi, lo, X[8], MAD_F(0x0f9ee890));
Gruenfrosch 0:7627c79db971 1969 MAD_F_MLA(hi, lo, X[9], MAD_F(0x03768962));
Gruenfrosch 0:7627c79db971 1970 MAD_F_MLA(hi, lo, X[11], -MAD_F(0x0e313245));
Gruenfrosch 0:7627c79db971 1971 MAD_F_MLA(hi, lo, X[12], MAD_F(0x00b2aa3e));
Gruenfrosch 0:7627c79db971 1972 MAD_F_MLA(hi, lo, X[14], MAD_F(0x0bcbe352));
Gruenfrosch 0:7627c79db971 1973 MAD_F_MLA(hi, lo, X[15], -MAD_F(0x04cfb0e2));
Gruenfrosch 0:7627c79db971 1974 MAD_F_MLA(hi, lo, X[17], -MAD_F(0x0898c779));
Gruenfrosch 0:7627c79db971 1975
Gruenfrosch 0:7627c79db971 1976 x[20] = x[33] = MAD_F_MLZ(hi, lo) - t3;
Gruenfrosch 0:7627c79db971 1977
Gruenfrosch 0:7627c79db971 1978 MAD_F_ML0(hi, lo, t14, -MAD_F(0x0ec835e8));
Gruenfrosch 0:7627c79db971 1979 MAD_F_MLA(hi, lo, t15, MAD_F(0x061f78aa));
Gruenfrosch 0:7627c79db971 1980
Gruenfrosch 0:7627c79db971 1981 t4 = MAD_F_MLZ(hi, lo) - t7;
Gruenfrosch 0:7627c79db971 1982
Gruenfrosch 0:7627c79db971 1983 MAD_F_ML0(hi, lo, t12, MAD_F(0x061f78aa));
Gruenfrosch 0:7627c79db971 1984 MAD_F_MLA(hi, lo, t13, MAD_F(0x0ec835e8));
Gruenfrosch 0:7627c79db971 1985
Gruenfrosch 0:7627c79db971 1986 x[4] = MAD_F_MLZ(hi, lo) + t4;
Gruenfrosch 0:7627c79db971 1987 x[13] = -x[4];
Gruenfrosch 0:7627c79db971 1988
Gruenfrosch 0:7627c79db971 1989 MAD_F_ML0(hi, lo, t8, MAD_F(0x09bd7ca0));
Gruenfrosch 0:7627c79db971 1990 MAD_F_MLA(hi, lo, t9, -MAD_F(0x0216a2a2));
Gruenfrosch 0:7627c79db971 1991 MAD_F_MLA(hi, lo, t10, MAD_F(0x0fdcf549));
Gruenfrosch 0:7627c79db971 1992 MAD_F_MLA(hi, lo, t11, -MAD_F(0x0cb19346));
Gruenfrosch 0:7627c79db971 1993
Gruenfrosch 0:7627c79db971 1994 x[1] = MAD_F_MLZ(hi, lo) + t4;
Gruenfrosch 0:7627c79db971 1995 x[16] = -x[1];
Gruenfrosch 0:7627c79db971 1996
Gruenfrosch 0:7627c79db971 1997 MAD_F_ML0(hi, lo, t8, -MAD_F(0x0fdcf549));
Gruenfrosch 0:7627c79db971 1998 MAD_F_MLA(hi, lo, t9, -MAD_F(0x0cb19346));
Gruenfrosch 0:7627c79db971 1999 MAD_F_MLA(hi, lo, t10, -MAD_F(0x09bd7ca0));
Gruenfrosch 0:7627c79db971 2000 MAD_F_MLA(hi, lo, t11, -MAD_F(0x0216a2a2));
Gruenfrosch 0:7627c79db971 2001
Gruenfrosch 0:7627c79db971 2002 x[25] = x[28] = MAD_F_MLZ(hi, lo) + t4;
Gruenfrosch 0:7627c79db971 2003
Gruenfrosch 0:7627c79db971 2004 MAD_F_ML0(hi, lo, X[1], -MAD_F(0x0fdcf549));
Gruenfrosch 0:7627c79db971 2005 MAD_F_MLA(hi, lo, X[7], -MAD_F(0x0cb19346));
Gruenfrosch 0:7627c79db971 2006 MAD_F_MLA(hi, lo, X[10], -MAD_F(0x09bd7ca0));
Gruenfrosch 0:7627c79db971 2007 MAD_F_MLA(hi, lo, X[16], -MAD_F(0x0216a2a2));
Gruenfrosch 0:7627c79db971 2008
Gruenfrosch 0:7627c79db971 2009 t5 = MAD_F_MLZ(hi, lo) - t6;
Gruenfrosch 0:7627c79db971 2010
Gruenfrosch 0:7627c79db971 2011 MAD_F_ML0(hi, lo, X[0], MAD_F(0x0898c779));
Gruenfrosch 0:7627c79db971 2012 MAD_F_MLA(hi, lo, X[2], MAD_F(0x04cfb0e2));
Gruenfrosch 0:7627c79db971 2013 MAD_F_MLA(hi, lo, X[3], MAD_F(0x0bcbe352));
Gruenfrosch 0:7627c79db971 2014 MAD_F_MLA(hi, lo, X[5], MAD_F(0x00b2aa3e));
Gruenfrosch 0:7627c79db971 2015 MAD_F_MLA(hi, lo, X[6], MAD_F(0x0e313245));
Gruenfrosch 0:7627c79db971 2016 MAD_F_MLA(hi, lo, X[8], -MAD_F(0x03768962));
Gruenfrosch 0:7627c79db971 2017 MAD_F_MLA(hi, lo, X[9], MAD_F(0x0f9ee890));
Gruenfrosch 0:7627c79db971 2018 MAD_F_MLA(hi, lo, X[11], -MAD_F(0x07635284));
Gruenfrosch 0:7627c79db971 2019 MAD_F_MLA(hi, lo, X[12], MAD_F(0x0ffc19fd));
Gruenfrosch 0:7627c79db971 2020 MAD_F_MLA(hi, lo, X[14], -MAD_F(0x0acf37ad));
Gruenfrosch 0:7627c79db971 2021 MAD_F_MLA(hi, lo, X[15], MAD_F(0x0f426cb5));
Gruenfrosch 0:7627c79db971 2022 MAD_F_MLA(hi, lo, X[17], -MAD_F(0x0d7e8807));
Gruenfrosch 0:7627c79db971 2023
Gruenfrosch 0:7627c79db971 2024 x[2] = MAD_F_MLZ(hi, lo) + t5;
Gruenfrosch 0:7627c79db971 2025 x[15] = -x[2];
Gruenfrosch 0:7627c79db971 2026
Gruenfrosch 0:7627c79db971 2027 MAD_F_ML0(hi, lo, X[0], MAD_F(0x07635284));
Gruenfrosch 0:7627c79db971 2028 MAD_F_MLA(hi, lo, X[2], MAD_F(0x0acf37ad));
Gruenfrosch 0:7627c79db971 2029 MAD_F_MLA(hi, lo, X[3], MAD_F(0x03768962));
Gruenfrosch 0:7627c79db971 2030 MAD_F_MLA(hi, lo, X[5], MAD_F(0x0d7e8807));
Gruenfrosch 0:7627c79db971 2031 MAD_F_MLA(hi, lo, X[6], -MAD_F(0x00b2aa3e));
Gruenfrosch 0:7627c79db971 2032 MAD_F_MLA(hi, lo, X[8], MAD_F(0x0f426cb5));
Gruenfrosch 0:7627c79db971 2033 MAD_F_MLA(hi, lo, X[9], -MAD_F(0x04cfb0e2));
Gruenfrosch 0:7627c79db971 2034 MAD_F_MLA(hi, lo, X[11], MAD_F(0x0ffc19fd));
Gruenfrosch 0:7627c79db971 2035 MAD_F_MLA(hi, lo, X[12], -MAD_F(0x0898c779));
Gruenfrosch 0:7627c79db971 2036 MAD_F_MLA(hi, lo, X[14], MAD_F(0x0f9ee890));
Gruenfrosch 0:7627c79db971 2037 MAD_F_MLA(hi, lo, X[15], -MAD_F(0x0bcbe352));
Gruenfrosch 0:7627c79db971 2038 MAD_F_MLA(hi, lo, X[17], MAD_F(0x0e313245));
Gruenfrosch 0:7627c79db971 2039
Gruenfrosch 0:7627c79db971 2040 x[3] = MAD_F_MLZ(hi, lo) + t5;
Gruenfrosch 0:7627c79db971 2041 x[14] = -x[3];
Gruenfrosch 0:7627c79db971 2042
Gruenfrosch 0:7627c79db971 2043 MAD_F_ML0(hi, lo, X[0], -MAD_F(0x0ffc19fd));
Gruenfrosch 0:7627c79db971 2044 MAD_F_MLA(hi, lo, X[2], -MAD_F(0x0f9ee890));
Gruenfrosch 0:7627c79db971 2045 MAD_F_MLA(hi, lo, X[3], -MAD_F(0x0f426cb5));
Gruenfrosch 0:7627c79db971 2046 MAD_F_MLA(hi, lo, X[5], -MAD_F(0x0e313245));
Gruenfrosch 0:7627c79db971 2047 MAD_F_MLA(hi, lo, X[6], -MAD_F(0x0d7e8807));
Gruenfrosch 0:7627c79db971 2048 MAD_F_MLA(hi, lo, X[8], -MAD_F(0x0bcbe352));
Gruenfrosch 0:7627c79db971 2049 MAD_F_MLA(hi, lo, X[9], -MAD_F(0x0acf37ad));
Gruenfrosch 0:7627c79db971 2050 MAD_F_MLA(hi, lo, X[11], -MAD_F(0x0898c779));
Gruenfrosch 0:7627c79db971 2051 MAD_F_MLA(hi, lo, X[12], -MAD_F(0x07635284));
Gruenfrosch 0:7627c79db971 2052 MAD_F_MLA(hi, lo, X[14], -MAD_F(0x04cfb0e2));
Gruenfrosch 0:7627c79db971 2053 MAD_F_MLA(hi, lo, X[15], -MAD_F(0x03768962));
Gruenfrosch 0:7627c79db971 2054 MAD_F_MLA(hi, lo, X[17], -MAD_F(0x00b2aa3e));
Gruenfrosch 0:7627c79db971 2055
Gruenfrosch 0:7627c79db971 2056 x[26] = x[27] = MAD_F_MLZ(hi, lo) + t5;
Gruenfrosch 0:7627c79db971 2057 }
Gruenfrosch 0:7627c79db971 2058 # endif
Gruenfrosch 0:7627c79db971 2059
Gruenfrosch 0:7627c79db971 2060 /*
Gruenfrosch 0:7627c79db971 2061 * NAME: III_imdct_l()
Gruenfrosch 0:7627c79db971 2062 * DESCRIPTION: perform IMDCT and windowing for long blocks
Gruenfrosch 0:7627c79db971 2063 */
Gruenfrosch 0:7627c79db971 2064 static
Gruenfrosch 0:7627c79db971 2065 void III_imdct_l(mad_fixed_t const X[18], mad_fixed_t z[36],
Gruenfrosch 0:7627c79db971 2066 unsigned int block_type)
Gruenfrosch 0:7627c79db971 2067 {
Gruenfrosch 0:7627c79db971 2068 unsigned int i;
Gruenfrosch 0:7627c79db971 2069
Gruenfrosch 0:7627c79db971 2070 /* IMDCT */
Gruenfrosch 0:7627c79db971 2071
Gruenfrosch 0:7627c79db971 2072 imdct36(X, z);
Gruenfrosch 0:7627c79db971 2073
Gruenfrosch 0:7627c79db971 2074 /* windowing */
Gruenfrosch 0:7627c79db971 2075
Gruenfrosch 0:7627c79db971 2076 switch (block_type) {
Gruenfrosch 0:7627c79db971 2077 case 0: /* normal window */
Gruenfrosch 0:7627c79db971 2078 # if defined(ASO_INTERLEAVE1)
Gruenfrosch 0:7627c79db971 2079 {
Gruenfrosch 0:7627c79db971 2080 register mad_fixed_t tmp1, tmp2;
Gruenfrosch 0:7627c79db971 2081
Gruenfrosch 0:7627c79db971 2082 tmp1 = window_l[0];
Gruenfrosch 0:7627c79db971 2083 tmp2 = window_l[1];
Gruenfrosch 0:7627c79db971 2084
Gruenfrosch 0:7627c79db971 2085 for (i = 0; i < 34; i += 2) {
Gruenfrosch 0:7627c79db971 2086 z[i + 0] = mad_f_mul(z[i + 0], tmp1);
Gruenfrosch 0:7627c79db971 2087 tmp1 = window_l[i + 2];
Gruenfrosch 0:7627c79db971 2088 z[i + 1] = mad_f_mul(z[i + 1], tmp2);
Gruenfrosch 0:7627c79db971 2089 tmp2 = window_l[i + 3];
Gruenfrosch 0:7627c79db971 2090 }
Gruenfrosch 0:7627c79db971 2091
Gruenfrosch 0:7627c79db971 2092 z[34] = mad_f_mul(z[34], tmp1);
Gruenfrosch 0:7627c79db971 2093 z[35] = mad_f_mul(z[35], tmp2);
Gruenfrosch 0:7627c79db971 2094 }
Gruenfrosch 0:7627c79db971 2095 # elif defined(ASO_INTERLEAVE2)
Gruenfrosch 0:7627c79db971 2096 {
Gruenfrosch 0:7627c79db971 2097 register mad_fixed_t tmp1, tmp2;
Gruenfrosch 0:7627c79db971 2098
Gruenfrosch 0:7627c79db971 2099 tmp1 = z[0];
Gruenfrosch 0:7627c79db971 2100 tmp2 = window_l[0];
Gruenfrosch 0:7627c79db971 2101
Gruenfrosch 0:7627c79db971 2102 for (i = 0; i < 35; ++i) {
Gruenfrosch 0:7627c79db971 2103 z[i] = mad_f_mul(tmp1, tmp2);
Gruenfrosch 0:7627c79db971 2104 tmp1 = z[i + 1];
Gruenfrosch 0:7627c79db971 2105 tmp2 = window_l[i + 1];
Gruenfrosch 0:7627c79db971 2106 }
Gruenfrosch 0:7627c79db971 2107
Gruenfrosch 0:7627c79db971 2108 z[35] = mad_f_mul(tmp1, tmp2);
Gruenfrosch 0:7627c79db971 2109 }
Gruenfrosch 0:7627c79db971 2110 # elif 1
Gruenfrosch 0:7627c79db971 2111 for (i = 0; i < 36; i += 4) {
Gruenfrosch 0:7627c79db971 2112 z[i + 0] = mad_f_mul(z[i + 0], window_l[i + 0]);
Gruenfrosch 0:7627c79db971 2113 z[i + 1] = mad_f_mul(z[i + 1], window_l[i + 1]);
Gruenfrosch 0:7627c79db971 2114 z[i + 2] = mad_f_mul(z[i + 2], window_l[i + 2]);
Gruenfrosch 0:7627c79db971 2115 z[i + 3] = mad_f_mul(z[i + 3], window_l[i + 3]);
Gruenfrosch 0:7627c79db971 2116 }
Gruenfrosch 0:7627c79db971 2117 # else
Gruenfrosch 0:7627c79db971 2118 for (i = 0; i < 36; ++i) z[i] = mad_f_mul(z[i], window_l[i]);
Gruenfrosch 0:7627c79db971 2119 # endif
Gruenfrosch 0:7627c79db971 2120 break;
Gruenfrosch 0:7627c79db971 2121
Gruenfrosch 0:7627c79db971 2122 case 1: /* start block */
Gruenfrosch 0:7627c79db971 2123 for (i = 0; i < 18; i += 3) {
Gruenfrosch 0:7627c79db971 2124 z[i + 0] = mad_f_mul(z[i + 0], window_l[i + 0]);
Gruenfrosch 0:7627c79db971 2125 z[i + 1] = mad_f_mul(z[i + 1], window_l[i + 1]);
Gruenfrosch 0:7627c79db971 2126 z[i + 2] = mad_f_mul(z[i + 2], window_l[i + 2]);
Gruenfrosch 0:7627c79db971 2127 }
Gruenfrosch 0:7627c79db971 2128 /* (i = 18; i < 24; ++i) z[i] unchanged */
Gruenfrosch 0:7627c79db971 2129 for (i = 24; i < 30; ++i) z[i] = mad_f_mul(z[i], window_s[i - 18]);
Gruenfrosch 0:7627c79db971 2130 for (i = 30; i < 36; ++i) z[i] = 0;
Gruenfrosch 0:7627c79db971 2131 break;
Gruenfrosch 0:7627c79db971 2132
Gruenfrosch 0:7627c79db971 2133 case 3: /* stop block */
Gruenfrosch 0:7627c79db971 2134 for (i = 0; i < 6; ++i) z[i] = 0;
Gruenfrosch 0:7627c79db971 2135 for (i = 6; i < 12; ++i) z[i] = mad_f_mul(z[i], window_s[i - 6]);
Gruenfrosch 0:7627c79db971 2136 /* (i = 12; i < 18; ++i) z[i] unchanged */
Gruenfrosch 0:7627c79db971 2137 for (i = 18; i < 36; i += 3) {
Gruenfrosch 0:7627c79db971 2138 z[i + 0] = mad_f_mul(z[i + 0], window_l[i + 0]);
Gruenfrosch 0:7627c79db971 2139 z[i + 1] = mad_f_mul(z[i + 1], window_l[i + 1]);
Gruenfrosch 0:7627c79db971 2140 z[i + 2] = mad_f_mul(z[i + 2], window_l[i + 2]);
Gruenfrosch 0:7627c79db971 2141 }
Gruenfrosch 0:7627c79db971 2142 break;
Gruenfrosch 0:7627c79db971 2143 }
Gruenfrosch 0:7627c79db971 2144 }
Gruenfrosch 0:7627c79db971 2145 # endif /* ASO_IMDCT */
Gruenfrosch 0:7627c79db971 2146
Gruenfrosch 0:7627c79db971 2147 /*
Gruenfrosch 0:7627c79db971 2148 * NAME: III_imdct_s()
Gruenfrosch 0:7627c79db971 2149 * DESCRIPTION: perform IMDCT and windowing for short blocks
Gruenfrosch 0:7627c79db971 2150 */
Gruenfrosch 0:7627c79db971 2151 static
Gruenfrosch 0:7627c79db971 2152 void III_imdct_s(mad_fixed_t const X[18], mad_fixed_t z[36])
Gruenfrosch 0:7627c79db971 2153 {
Gruenfrosch 0:7627c79db971 2154 mad_fixed_t y[36], *yptr;
Gruenfrosch 0:7627c79db971 2155 mad_fixed_t const *wptr;
Gruenfrosch 0:7627c79db971 2156 int w, i;
Gruenfrosch 0:7627c79db971 2157 register mad_fixed64hi_t hi;
Gruenfrosch 0:7627c79db971 2158 register mad_fixed64lo_t lo;
Gruenfrosch 0:7627c79db971 2159
Gruenfrosch 0:7627c79db971 2160 /* IMDCT */
Gruenfrosch 0:7627c79db971 2161
Gruenfrosch 0:7627c79db971 2162 yptr = &y[0];
Gruenfrosch 0:7627c79db971 2163
Gruenfrosch 0:7627c79db971 2164 for (w = 0; w < 3; ++w) {
Gruenfrosch 0:7627c79db971 2165 register mad_fixed_t const (*s)[6];
Gruenfrosch 0:7627c79db971 2166
Gruenfrosch 0:7627c79db971 2167 s = imdct_s;
Gruenfrosch 0:7627c79db971 2168
Gruenfrosch 0:7627c79db971 2169 for (i = 0; i < 3; ++i) {
Gruenfrosch 0:7627c79db971 2170 MAD_F_ML0(hi, lo, X[0], (*s)[0]);
Gruenfrosch 0:7627c79db971 2171 MAD_F_MLA(hi, lo, X[1], (*s)[1]);
Gruenfrosch 0:7627c79db971 2172 MAD_F_MLA(hi, lo, X[2], (*s)[2]);
Gruenfrosch 0:7627c79db971 2173 MAD_F_MLA(hi, lo, X[3], (*s)[3]);
Gruenfrosch 0:7627c79db971 2174 MAD_F_MLA(hi, lo, X[4], (*s)[4]);
Gruenfrosch 0:7627c79db971 2175 MAD_F_MLA(hi, lo, X[5], (*s)[5]);
Gruenfrosch 0:7627c79db971 2176
Gruenfrosch 0:7627c79db971 2177 yptr[i + 0] = MAD_F_MLZ(hi, lo);
Gruenfrosch 0:7627c79db971 2178 yptr[5 - i] = -yptr[i + 0];
Gruenfrosch 0:7627c79db971 2179
Gruenfrosch 0:7627c79db971 2180 ++s;
Gruenfrosch 0:7627c79db971 2181
Gruenfrosch 0:7627c79db971 2182 MAD_F_ML0(hi, lo, X[0], (*s)[0]);
Gruenfrosch 0:7627c79db971 2183 MAD_F_MLA(hi, lo, X[1], (*s)[1]);
Gruenfrosch 0:7627c79db971 2184 MAD_F_MLA(hi, lo, X[2], (*s)[2]);
Gruenfrosch 0:7627c79db971 2185 MAD_F_MLA(hi, lo, X[3], (*s)[3]);
Gruenfrosch 0:7627c79db971 2186 MAD_F_MLA(hi, lo, X[4], (*s)[4]);
Gruenfrosch 0:7627c79db971 2187 MAD_F_MLA(hi, lo, X[5], (*s)[5]);
Gruenfrosch 0:7627c79db971 2188
Gruenfrosch 0:7627c79db971 2189 yptr[ i + 6] = MAD_F_MLZ(hi, lo);
Gruenfrosch 0:7627c79db971 2190 yptr[11 - i] = yptr[i + 6];
Gruenfrosch 0:7627c79db971 2191
Gruenfrosch 0:7627c79db971 2192 ++s;
Gruenfrosch 0:7627c79db971 2193 }
Gruenfrosch 0:7627c79db971 2194
Gruenfrosch 0:7627c79db971 2195 yptr += 12;
Gruenfrosch 0:7627c79db971 2196 X += 6;
Gruenfrosch 0:7627c79db971 2197 }
Gruenfrosch 0:7627c79db971 2198
Gruenfrosch 0:7627c79db971 2199 /* windowing, overlapping and concatenation */
Gruenfrosch 0:7627c79db971 2200
Gruenfrosch 0:7627c79db971 2201 yptr = &y[0];
Gruenfrosch 0:7627c79db971 2202 wptr = &window_s[0];
Gruenfrosch 0:7627c79db971 2203
Gruenfrosch 0:7627c79db971 2204 for (i = 0; i < 6; ++i) {
Gruenfrosch 0:7627c79db971 2205 z[i + 0] = 0;
Gruenfrosch 0:7627c79db971 2206 z[i + 6] = mad_f_mul(yptr[ 0 + 0], wptr[0]);
Gruenfrosch 0:7627c79db971 2207
Gruenfrosch 0:7627c79db971 2208 MAD_F_ML0(hi, lo, yptr[ 0 + 6], wptr[6]);
Gruenfrosch 0:7627c79db971 2209 MAD_F_MLA(hi, lo, yptr[12 + 0], wptr[0]);
Gruenfrosch 0:7627c79db971 2210
Gruenfrosch 0:7627c79db971 2211 z[i + 12] = MAD_F_MLZ(hi, lo);
Gruenfrosch 0:7627c79db971 2212
Gruenfrosch 0:7627c79db971 2213 MAD_F_ML0(hi, lo, yptr[12 + 6], wptr[6]);
Gruenfrosch 0:7627c79db971 2214 MAD_F_MLA(hi, lo, yptr[24 + 0], wptr[0]);
Gruenfrosch 0:7627c79db971 2215
Gruenfrosch 0:7627c79db971 2216 z[i + 18] = MAD_F_MLZ(hi, lo);
Gruenfrosch 0:7627c79db971 2217
Gruenfrosch 0:7627c79db971 2218 z[i + 24] = mad_f_mul(yptr[24 + 6], wptr[6]);
Gruenfrosch 0:7627c79db971 2219 z[i + 30] = 0;
Gruenfrosch 0:7627c79db971 2220
Gruenfrosch 0:7627c79db971 2221 ++yptr;
Gruenfrosch 0:7627c79db971 2222 ++wptr;
Gruenfrosch 0:7627c79db971 2223 }
Gruenfrosch 0:7627c79db971 2224 }
Gruenfrosch 0:7627c79db971 2225
Gruenfrosch 0:7627c79db971 2226 /*
Gruenfrosch 0:7627c79db971 2227 * NAME: III_overlap()
Gruenfrosch 0:7627c79db971 2228 * DESCRIPTION: perform overlap-add of windowed IMDCT outputs
Gruenfrosch 0:7627c79db971 2229 */
Gruenfrosch 0:7627c79db971 2230 static
Gruenfrosch 0:7627c79db971 2231 void III_overlap(mad_fixed_t const output[36], mad_fixed_t overlap[18],
Gruenfrosch 0:7627c79db971 2232 mad_fixed_t sample[18][32], unsigned int sb)
Gruenfrosch 0:7627c79db971 2233 {
Gruenfrosch 0:7627c79db971 2234 unsigned int i;
Gruenfrosch 0:7627c79db971 2235
Gruenfrosch 0:7627c79db971 2236 # if defined(ASO_INTERLEAVE2)
Gruenfrosch 0:7627c79db971 2237 {
Gruenfrosch 0:7627c79db971 2238 register mad_fixed_t tmp1, tmp2;
Gruenfrosch 0:7627c79db971 2239
Gruenfrosch 0:7627c79db971 2240 tmp1 = overlap[0];
Gruenfrosch 0:7627c79db971 2241 tmp2 = overlap[1];
Gruenfrosch 0:7627c79db971 2242
Gruenfrosch 0:7627c79db971 2243 for (i = 0; i < 16; i += 2) {
Gruenfrosch 0:7627c79db971 2244 sample[i + 0][sb] = output[i + 0 + 0] + tmp1;
Gruenfrosch 0:7627c79db971 2245 overlap[i + 0] = output[i + 0 + 18];
Gruenfrosch 0:7627c79db971 2246 tmp1 = overlap[i + 2];
Gruenfrosch 0:7627c79db971 2247
Gruenfrosch 0:7627c79db971 2248 sample[i + 1][sb] = output[i + 1 + 0] + tmp2;
Gruenfrosch 0:7627c79db971 2249 overlap[i + 1] = output[i + 1 + 18];
Gruenfrosch 0:7627c79db971 2250 tmp2 = overlap[i + 3];
Gruenfrosch 0:7627c79db971 2251 }
Gruenfrosch 0:7627c79db971 2252
Gruenfrosch 0:7627c79db971 2253 sample[16][sb] = output[16 + 0] + tmp1;
Gruenfrosch 0:7627c79db971 2254 overlap[16] = output[16 + 18];
Gruenfrosch 0:7627c79db971 2255 sample[17][sb] = output[17 + 0] + tmp2;
Gruenfrosch 0:7627c79db971 2256 overlap[17] = output[17 + 18];
Gruenfrosch 0:7627c79db971 2257 }
Gruenfrosch 0:7627c79db971 2258 # elif 0
Gruenfrosch 0:7627c79db971 2259 for (i = 0; i < 18; i += 2) {
Gruenfrosch 0:7627c79db971 2260 sample[i + 0][sb] = output[i + 0 + 0] + overlap[i + 0];
Gruenfrosch 0:7627c79db971 2261 overlap[i + 0] = output[i + 0 + 18];
Gruenfrosch 0:7627c79db971 2262
Gruenfrosch 0:7627c79db971 2263 sample[i + 1][sb] = output[i + 1 + 0] + overlap[i + 1];
Gruenfrosch 0:7627c79db971 2264 overlap[i + 1] = output[i + 1 + 18];
Gruenfrosch 0:7627c79db971 2265 }
Gruenfrosch 0:7627c79db971 2266 # else
Gruenfrosch 0:7627c79db971 2267 for (i = 0; i < 18; ++i) {
Gruenfrosch 0:7627c79db971 2268 sample[i][sb] = output[i + 0] + overlap[i];
Gruenfrosch 0:7627c79db971 2269 overlap[i] = output[i + 18];
Gruenfrosch 0:7627c79db971 2270 }
Gruenfrosch 0:7627c79db971 2271 # endif
Gruenfrosch 0:7627c79db971 2272 }
Gruenfrosch 0:7627c79db971 2273
Gruenfrosch 0:7627c79db971 2274 /*
Gruenfrosch 0:7627c79db971 2275 * NAME: III_overlap_z()
Gruenfrosch 0:7627c79db971 2276 * DESCRIPTION: perform "overlap-add" of zero IMDCT outputs
Gruenfrosch 0:7627c79db971 2277 */
Gruenfrosch 0:7627c79db971 2278 static inline
Gruenfrosch 0:7627c79db971 2279 void III_overlap_z(mad_fixed_t overlap[18],
Gruenfrosch 0:7627c79db971 2280 mad_fixed_t sample[18][32], unsigned int sb)
Gruenfrosch 0:7627c79db971 2281 {
Gruenfrosch 0:7627c79db971 2282 unsigned int i;
Gruenfrosch 0:7627c79db971 2283
Gruenfrosch 0:7627c79db971 2284 # if defined(ASO_INTERLEAVE2)
Gruenfrosch 0:7627c79db971 2285 {
Gruenfrosch 0:7627c79db971 2286 register mad_fixed_t tmp1, tmp2;
Gruenfrosch 0:7627c79db971 2287
Gruenfrosch 0:7627c79db971 2288 tmp1 = overlap[0];
Gruenfrosch 0:7627c79db971 2289 tmp2 = overlap[1];
Gruenfrosch 0:7627c79db971 2290
Gruenfrosch 0:7627c79db971 2291 for (i = 0; i < 16; i += 2) {
Gruenfrosch 0:7627c79db971 2292 sample[i + 0][sb] = tmp1;
Gruenfrosch 0:7627c79db971 2293 overlap[i + 0] = 0;
Gruenfrosch 0:7627c79db971 2294 tmp1 = overlap[i + 2];
Gruenfrosch 0:7627c79db971 2295
Gruenfrosch 0:7627c79db971 2296 sample[i + 1][sb] = tmp2;
Gruenfrosch 0:7627c79db971 2297 overlap[i + 1] = 0;
Gruenfrosch 0:7627c79db971 2298 tmp2 = overlap[i + 3];
Gruenfrosch 0:7627c79db971 2299 }
Gruenfrosch 0:7627c79db971 2300
Gruenfrosch 0:7627c79db971 2301 sample[16][sb] = tmp1;
Gruenfrosch 0:7627c79db971 2302 overlap[16] = 0;
Gruenfrosch 0:7627c79db971 2303 sample[17][sb] = tmp2;
Gruenfrosch 0:7627c79db971 2304 overlap[17] = 0;
Gruenfrosch 0:7627c79db971 2305 }
Gruenfrosch 0:7627c79db971 2306 # else
Gruenfrosch 0:7627c79db971 2307 for (i = 0; i < 18; ++i) {
Gruenfrosch 0:7627c79db971 2308 sample[i][sb] = overlap[i];
Gruenfrosch 0:7627c79db971 2309 overlap[i] = 0;
Gruenfrosch 0:7627c79db971 2310 }
Gruenfrosch 0:7627c79db971 2311 # endif
Gruenfrosch 0:7627c79db971 2312 }
Gruenfrosch 0:7627c79db971 2313
Gruenfrosch 0:7627c79db971 2314 /*
Gruenfrosch 0:7627c79db971 2315 * NAME: III_freqinver()
Gruenfrosch 0:7627c79db971 2316 * DESCRIPTION: perform subband frequency inversion for odd sample lines
Gruenfrosch 0:7627c79db971 2317 */
Gruenfrosch 0:7627c79db971 2318 static
Gruenfrosch 0:7627c79db971 2319 void III_freqinver(mad_fixed_t sample[18][32], unsigned int sb)
Gruenfrosch 0:7627c79db971 2320 {
Gruenfrosch 0:7627c79db971 2321 unsigned int i;
Gruenfrosch 0:7627c79db971 2322
Gruenfrosch 0:7627c79db971 2323 # if 1 || defined(ASO_INTERLEAVE1) || defined(ASO_INTERLEAVE2)
Gruenfrosch 0:7627c79db971 2324 {
Gruenfrosch 0:7627c79db971 2325 register mad_fixed_t tmp1, tmp2;
Gruenfrosch 0:7627c79db971 2326
Gruenfrosch 0:7627c79db971 2327 tmp1 = sample[1][sb];
Gruenfrosch 0:7627c79db971 2328 tmp2 = sample[3][sb];
Gruenfrosch 0:7627c79db971 2329
Gruenfrosch 0:7627c79db971 2330 for (i = 1; i < 13; i += 4) {
Gruenfrosch 0:7627c79db971 2331 sample[i + 0][sb] = -tmp1;
Gruenfrosch 0:7627c79db971 2332 tmp1 = sample[i + 4][sb];
Gruenfrosch 0:7627c79db971 2333 sample[i + 2][sb] = -tmp2;
Gruenfrosch 0:7627c79db971 2334 tmp2 = sample[i + 6][sb];
Gruenfrosch 0:7627c79db971 2335 }
Gruenfrosch 0:7627c79db971 2336
Gruenfrosch 0:7627c79db971 2337 sample[13][sb] = -tmp1;
Gruenfrosch 0:7627c79db971 2338 tmp1 = sample[17][sb];
Gruenfrosch 0:7627c79db971 2339 sample[15][sb] = -tmp2;
Gruenfrosch 0:7627c79db971 2340 sample[17][sb] = -tmp1;
Gruenfrosch 0:7627c79db971 2341 }
Gruenfrosch 0:7627c79db971 2342 # else
Gruenfrosch 0:7627c79db971 2343 for (i = 1; i < 18; i += 2)
Gruenfrosch 0:7627c79db971 2344 sample[i][sb] = -sample[i][sb];
Gruenfrosch 0:7627c79db971 2345 # endif
Gruenfrosch 0:7627c79db971 2346 }
Gruenfrosch 0:7627c79db971 2347
Gruenfrosch 0:7627c79db971 2348 /*
Gruenfrosch 0:7627c79db971 2349 * NAME: III_decode()
Gruenfrosch 0:7627c79db971 2350 * DESCRIPTION: decode frame main_data
Gruenfrosch 0:7627c79db971 2351 */
Gruenfrosch 0:7627c79db971 2352 static
Gruenfrosch 0:7627c79db971 2353 enum mad_error III_decode(struct mad_bitptr *ptr, struct mad_frame *frame,
Gruenfrosch 0:7627c79db971 2354 struct sideinfo *si, unsigned int nch)
Gruenfrosch 0:7627c79db971 2355 {
Gruenfrosch 0:7627c79db971 2356 struct mad_header *header = &frame->header;
Gruenfrosch 0:7627c79db971 2357 unsigned int sfreqi, ngr, gr;
Gruenfrosch 0:7627c79db971 2358
Gruenfrosch 0:7627c79db971 2359 {
Gruenfrosch 0:7627c79db971 2360 unsigned int sfreq;
Gruenfrosch 0:7627c79db971 2361
Gruenfrosch 0:7627c79db971 2362 sfreq = header->samplerate;
Gruenfrosch 0:7627c79db971 2363 if (header->flags & MAD_FLAG_MPEG_2_5_EXT)
Gruenfrosch 0:7627c79db971 2364 sfreq *= 2;
Gruenfrosch 0:7627c79db971 2365
Gruenfrosch 0:7627c79db971 2366 /* 48000 => 0, 44100 => 1, 32000 => 2,
Gruenfrosch 0:7627c79db971 2367 24000 => 3, 22050 => 4, 16000 => 5 */
Gruenfrosch 0:7627c79db971 2368 sfreqi = ((sfreq >> 7) & 0x000f) +
Gruenfrosch 0:7627c79db971 2369 ((sfreq >> 15) & 0x0001) - 8;
Gruenfrosch 0:7627c79db971 2370
Gruenfrosch 0:7627c79db971 2371 if (header->flags & MAD_FLAG_MPEG_2_5_EXT)
Gruenfrosch 0:7627c79db971 2372 sfreqi += 3;
Gruenfrosch 0:7627c79db971 2373 }
Gruenfrosch 0:7627c79db971 2374
Gruenfrosch 0:7627c79db971 2375 /* scalefactors, Huffman decoding, requantization */
Gruenfrosch 0:7627c79db971 2376
Gruenfrosch 0:7627c79db971 2377 ngr = (header->flags & MAD_FLAG_LSF_EXT) ? 1 : 2;
Gruenfrosch 0:7627c79db971 2378
Gruenfrosch 0:7627c79db971 2379 for (gr = 0; gr < ngr; ++gr) {
Gruenfrosch 0:7627c79db971 2380 struct granule *granule = &si->gr[gr];
Gruenfrosch 0:7627c79db971 2381 unsigned char const *sfbwidth[2];
Gruenfrosch 0:7627c79db971 2382 mad_fixed_t xr[2][576];
Gruenfrosch 0:7627c79db971 2383 unsigned int ch;
Gruenfrosch 0:7627c79db971 2384 enum mad_error error;
Gruenfrosch 0:7627c79db971 2385
Gruenfrosch 0:7627c79db971 2386 for (ch = 0; ch < nch; ++ch) {
Gruenfrosch 0:7627c79db971 2387 struct channel *channel = &granule->ch[ch];
Gruenfrosch 0:7627c79db971 2388 unsigned int part2_length;
Gruenfrosch 0:7627c79db971 2389
Gruenfrosch 0:7627c79db971 2390 sfbwidth[ch] = sfbwidth_table[sfreqi].l;
Gruenfrosch 0:7627c79db971 2391 if (channel->block_type == 2) {
Gruenfrosch 0:7627c79db971 2392 sfbwidth[ch] = (channel->flags & mixed_block_flag) ?
Gruenfrosch 0:7627c79db971 2393 sfbwidth_table[sfreqi].m : sfbwidth_table[sfreqi].s;
Gruenfrosch 0:7627c79db971 2394 }
Gruenfrosch 0:7627c79db971 2395
Gruenfrosch 0:7627c79db971 2396 if (header->flags & MAD_FLAG_LSF_EXT) {
Gruenfrosch 0:7627c79db971 2397 part2_length = III_scalefactors_lsf(ptr, channel,
Gruenfrosch 0:7627c79db971 2398 ch == 0 ? 0 : &si->gr[1].ch[1],
Gruenfrosch 0:7627c79db971 2399 header->mode_extension);
Gruenfrosch 0:7627c79db971 2400 }
Gruenfrosch 0:7627c79db971 2401 else {
Gruenfrosch 0:7627c79db971 2402 part2_length = III_scalefactors(ptr, channel, &si->gr[0].ch[ch],
Gruenfrosch 0:7627c79db971 2403 gr == 0 ? 0 : si->scfsi[ch]);
Gruenfrosch 0:7627c79db971 2404 }
Gruenfrosch 0:7627c79db971 2405
Gruenfrosch 0:7627c79db971 2406 error = III_huffdecode(ptr, xr[ch], channel, sfbwidth[ch], part2_length);
Gruenfrosch 0:7627c79db971 2407 if (error)
Gruenfrosch 0:7627c79db971 2408 return error;
Gruenfrosch 0:7627c79db971 2409 }
Gruenfrosch 0:7627c79db971 2410
Gruenfrosch 0:7627c79db971 2411 /* joint stereo processing */
Gruenfrosch 0:7627c79db971 2412
Gruenfrosch 0:7627c79db971 2413 if (header->mode == MAD_MODE_JOINT_STEREO && header->mode_extension) {
Gruenfrosch 0:7627c79db971 2414 error = III_stereo(xr, granule, header, sfbwidth[0]);
Gruenfrosch 0:7627c79db971 2415 if (error)
Gruenfrosch 0:7627c79db971 2416 return error;
Gruenfrosch 0:7627c79db971 2417 }
Gruenfrosch 0:7627c79db971 2418
Gruenfrosch 0:7627c79db971 2419 /* reordering, alias reduction, IMDCT, overlap-add, frequency inversion */
Gruenfrosch 0:7627c79db971 2420
Gruenfrosch 0:7627c79db971 2421 for (ch = 0; ch < nch; ++ch) {
Gruenfrosch 0:7627c79db971 2422 struct channel const *channel = &granule->ch[ch];
Gruenfrosch 0:7627c79db971 2423 mad_fixed_t (*sample)[32] = &frame->sbsample[ch][18 * gr];
Gruenfrosch 0:7627c79db971 2424 unsigned int sb, l, i, sblimit;
Gruenfrosch 0:7627c79db971 2425 mad_fixed_t output[36];
Gruenfrosch 0:7627c79db971 2426
Gruenfrosch 0:7627c79db971 2427 if (channel->block_type == 2) {
Gruenfrosch 0:7627c79db971 2428 III_reorder(xr[ch], channel, sfbwidth[ch]);
Gruenfrosch 0:7627c79db971 2429
Gruenfrosch 0:7627c79db971 2430 # if !defined(OPT_STRICT)
Gruenfrosch 0:7627c79db971 2431 /*
Gruenfrosch 0:7627c79db971 2432 * According to ISO/IEC 11172-3, "Alias reduction is not applied for
Gruenfrosch 0:7627c79db971 2433 * granules with block_type == 2 (short block)." However, other
Gruenfrosch 0:7627c79db971 2434 * sources suggest alias reduction should indeed be performed on the
Gruenfrosch 0:7627c79db971 2435 * lower two subbands of mixed blocks. Most other implementations do
Gruenfrosch 0:7627c79db971 2436 * this, so by default we will too.
Gruenfrosch 0:7627c79db971 2437 */
Gruenfrosch 0:7627c79db971 2438 if (channel->flags & mixed_block_flag)
Gruenfrosch 0:7627c79db971 2439 III_aliasreduce(xr[ch], 36);
Gruenfrosch 0:7627c79db971 2440 # endif
Gruenfrosch 0:7627c79db971 2441 }
Gruenfrosch 0:7627c79db971 2442 else
Gruenfrosch 0:7627c79db971 2443 III_aliasreduce(xr[ch], 576);
Gruenfrosch 0:7627c79db971 2444
Gruenfrosch 0:7627c79db971 2445 l = 0;
Gruenfrosch 0:7627c79db971 2446
Gruenfrosch 0:7627c79db971 2447 /* subbands 0-1 */
Gruenfrosch 0:7627c79db971 2448
Gruenfrosch 0:7627c79db971 2449 if (channel->block_type != 2 || (channel->flags & mixed_block_flag)) {
Gruenfrosch 0:7627c79db971 2450 unsigned int block_type;
Gruenfrosch 0:7627c79db971 2451
Gruenfrosch 0:7627c79db971 2452 block_type = channel->block_type;
Gruenfrosch 0:7627c79db971 2453 if (channel->flags & mixed_block_flag)
Gruenfrosch 0:7627c79db971 2454 block_type = 0;
Gruenfrosch 0:7627c79db971 2455
Gruenfrosch 0:7627c79db971 2456 /* long blocks */
Gruenfrosch 0:7627c79db971 2457 for (sb = 0; sb < 2; ++sb, l += 18) {
Gruenfrosch 0:7627c79db971 2458 III_imdct_l(&xr[ch][l], output, block_type);
Gruenfrosch 0:7627c79db971 2459 III_overlap(output, (*frame->overlap)[ch][sb], sample, sb);
Gruenfrosch 0:7627c79db971 2460 }
Gruenfrosch 0:7627c79db971 2461 }
Gruenfrosch 0:7627c79db971 2462 else {
Gruenfrosch 0:7627c79db971 2463 /* short blocks */
Gruenfrosch 0:7627c79db971 2464 for (sb = 0; sb < 2; ++sb, l += 18) {
Gruenfrosch 0:7627c79db971 2465 III_imdct_s(&xr[ch][l], output);
Gruenfrosch 0:7627c79db971 2466 III_overlap(output, (*frame->overlap)[ch][sb], sample, sb);
Gruenfrosch 0:7627c79db971 2467 }
Gruenfrosch 0:7627c79db971 2468 }
Gruenfrosch 0:7627c79db971 2469
Gruenfrosch 0:7627c79db971 2470 III_freqinver(sample, 1);
Gruenfrosch 0:7627c79db971 2471
Gruenfrosch 0:7627c79db971 2472 /* (nonzero) subbands 2-31 */
Gruenfrosch 0:7627c79db971 2473
Gruenfrosch 0:7627c79db971 2474 i = 576;
Gruenfrosch 0:7627c79db971 2475 while (i > 36 && xr[ch][i - 1] == 0)
Gruenfrosch 0:7627c79db971 2476 --i;
Gruenfrosch 0:7627c79db971 2477
Gruenfrosch 0:7627c79db971 2478 sblimit = 32 - (576 - i) / 18;
Gruenfrosch 0:7627c79db971 2479
Gruenfrosch 0:7627c79db971 2480 if (channel->block_type != 2) {
Gruenfrosch 0:7627c79db971 2481 /* long blocks */
Gruenfrosch 0:7627c79db971 2482 for (sb = 2; sb < sblimit; ++sb, l += 18) {
Gruenfrosch 0:7627c79db971 2483 III_imdct_l(&xr[ch][l], output, channel->block_type);
Gruenfrosch 0:7627c79db971 2484 III_overlap(output, (*frame->overlap)[ch][sb], sample, sb);
Gruenfrosch 0:7627c79db971 2485
Gruenfrosch 0:7627c79db971 2486 if (sb & 1)
Gruenfrosch 0:7627c79db971 2487 III_freqinver(sample, sb);
Gruenfrosch 0:7627c79db971 2488 }
Gruenfrosch 0:7627c79db971 2489 }
Gruenfrosch 0:7627c79db971 2490 else {
Gruenfrosch 0:7627c79db971 2491 /* short blocks */
Gruenfrosch 0:7627c79db971 2492 for (sb = 2; sb < sblimit; ++sb, l += 18) {
Gruenfrosch 0:7627c79db971 2493 III_imdct_s(&xr[ch][l], output);
Gruenfrosch 0:7627c79db971 2494 III_overlap(output, (*frame->overlap)[ch][sb], sample, sb);
Gruenfrosch 0:7627c79db971 2495
Gruenfrosch 0:7627c79db971 2496 if (sb & 1)
Gruenfrosch 0:7627c79db971 2497 III_freqinver(sample, sb);
Gruenfrosch 0:7627c79db971 2498 }
Gruenfrosch 0:7627c79db971 2499 }
Gruenfrosch 0:7627c79db971 2500
Gruenfrosch 0:7627c79db971 2501 /* remaining (zero) subbands */
Gruenfrosch 0:7627c79db971 2502
Gruenfrosch 0:7627c79db971 2503 for (sb = sblimit; sb < 32; ++sb) {
Gruenfrosch 0:7627c79db971 2504 III_overlap_z((*frame->overlap)[ch][sb], sample, sb);
Gruenfrosch 0:7627c79db971 2505
Gruenfrosch 0:7627c79db971 2506 if (sb & 1)
Gruenfrosch 0:7627c79db971 2507 III_freqinver(sample, sb);
Gruenfrosch 0:7627c79db971 2508 }
Gruenfrosch 0:7627c79db971 2509 }
Gruenfrosch 0:7627c79db971 2510 }
Gruenfrosch 0:7627c79db971 2511
Gruenfrosch 0:7627c79db971 2512 return MAD_ERROR_NONE;
Gruenfrosch 0:7627c79db971 2513 }
Gruenfrosch 0:7627c79db971 2514
Gruenfrosch 0:7627c79db971 2515 /*
Gruenfrosch 0:7627c79db971 2516 * NAME: layer->III()
Gruenfrosch 0:7627c79db971 2517 * DESCRIPTION: decode a single Layer III frame
Gruenfrosch 0:7627c79db971 2518 */
Gruenfrosch 0:7627c79db971 2519 int mad_layer_III(struct mad_stream *stream, struct mad_frame *frame)
Gruenfrosch 0:7627c79db971 2520 {
Gruenfrosch 0:7627c79db971 2521 struct mad_header *header = &frame->header;
Gruenfrosch 0:7627c79db971 2522 unsigned int nch, priv_bitlen, next_md_begin = 0;
Gruenfrosch 0:7627c79db971 2523 unsigned int si_len, data_bitlen, md_len;
Gruenfrosch 0:7627c79db971 2524 unsigned int frame_space, frame_used, frame_free;
Gruenfrosch 0:7627c79db971 2525 struct mad_bitptr ptr;
Gruenfrosch 0:7627c79db971 2526 struct sideinfo si;
Gruenfrosch 0:7627c79db971 2527 enum mad_error error;
Gruenfrosch 0:7627c79db971 2528 int result = 0;
Gruenfrosch 0:7627c79db971 2529
Gruenfrosch 0:7627c79db971 2530 /* allocate Layer III dynamic structures */
Gruenfrosch 0:7627c79db971 2531
Gruenfrosch 0:7627c79db971 2532 if (stream->main_data == 0) {
Gruenfrosch 0:7627c79db971 2533 stream->main_data = (unsigned char (*)[MAD_BUFFER_MDLEN])malloc(MAD_BUFFER_MDLEN);
Gruenfrosch 0:7627c79db971 2534 if (stream->main_data == 0) {
Gruenfrosch 0:7627c79db971 2535 stream->error = MAD_ERROR_NOMEM;
Gruenfrosch 0:7627c79db971 2536 return -1;
Gruenfrosch 0:7627c79db971 2537 }
Gruenfrosch 0:7627c79db971 2538 }
Gruenfrosch 0:7627c79db971 2539
Gruenfrosch 0:7627c79db971 2540 if (frame->overlap == 0) {
Gruenfrosch 0:7627c79db971 2541 frame->overlap = (mad_fixed_t (*)[2][32][18])calloc(2 * 32 * 18, sizeof(mad_fixed_t));
Gruenfrosch 0:7627c79db971 2542 if (frame->overlap == 0) {
Gruenfrosch 0:7627c79db971 2543 stream->error = MAD_ERROR_NOMEM;
Gruenfrosch 0:7627c79db971 2544 return -1;
Gruenfrosch 0:7627c79db971 2545 }
Gruenfrosch 0:7627c79db971 2546 }
Gruenfrosch 0:7627c79db971 2547
Gruenfrosch 0:7627c79db971 2548 nch = MAD_NCHANNELS(header);
Gruenfrosch 0:7627c79db971 2549 si_len = (header->flags & MAD_FLAG_LSF_EXT) ?
Gruenfrosch 0:7627c79db971 2550 (nch == 1 ? 9 : 17) : (nch == 1 ? 17 : 32);
Gruenfrosch 0:7627c79db971 2551
Gruenfrosch 0:7627c79db971 2552 /* check frame sanity */
Gruenfrosch 0:7627c79db971 2553
Gruenfrosch 0:7627c79db971 2554 if (stream->next_frame - mad_bit_nextbyte(&stream->ptr) <
Gruenfrosch 0:7627c79db971 2555 (signed int) si_len) {
Gruenfrosch 0:7627c79db971 2556 stream->error = MAD_ERROR_BADFRAMELEN;
Gruenfrosch 0:7627c79db971 2557 stream->md_len = 0;
Gruenfrosch 0:7627c79db971 2558 return -1;
Gruenfrosch 0:7627c79db971 2559 }
Gruenfrosch 0:7627c79db971 2560
Gruenfrosch 0:7627c79db971 2561 /* check CRC word */
Gruenfrosch 0:7627c79db971 2562
Gruenfrosch 0:7627c79db971 2563 if (header->flags & MAD_FLAG_PROTECTION) {
Gruenfrosch 0:7627c79db971 2564 header->crc_check =
Gruenfrosch 0:7627c79db971 2565 mad_bit_crc(stream->ptr, si_len * CHAR_BIT, header->crc_check);
Gruenfrosch 0:7627c79db971 2566
Gruenfrosch 0:7627c79db971 2567 if (header->crc_check != header->crc_target &&
Gruenfrosch 0:7627c79db971 2568 !(frame->options & MAD_OPTION_IGNORECRC)) {
Gruenfrosch 0:7627c79db971 2569 stream->error = MAD_ERROR_BADCRC;
Gruenfrosch 0:7627c79db971 2570 result = -1;
Gruenfrosch 0:7627c79db971 2571 }
Gruenfrosch 0:7627c79db971 2572 }
Gruenfrosch 0:7627c79db971 2573
Gruenfrosch 0:7627c79db971 2574 /* decode frame side information */
Gruenfrosch 0:7627c79db971 2575
Gruenfrosch 0:7627c79db971 2576 error = III_sideinfo(&stream->ptr, nch, header->flags & MAD_FLAG_LSF_EXT,
Gruenfrosch 0:7627c79db971 2577 &si, &data_bitlen, &priv_bitlen);
Gruenfrosch 0:7627c79db971 2578 if (error && result == 0) {
Gruenfrosch 0:7627c79db971 2579 stream->error = error;
Gruenfrosch 0:7627c79db971 2580 result = -1;
Gruenfrosch 0:7627c79db971 2581 }
Gruenfrosch 0:7627c79db971 2582
Gruenfrosch 0:7627c79db971 2583 header->flags |= priv_bitlen;
Gruenfrosch 0:7627c79db971 2584 header->private_bits |= si.private_bits;
Gruenfrosch 0:7627c79db971 2585
Gruenfrosch 0:7627c79db971 2586 /* find main_data of next frame */
Gruenfrosch 0:7627c79db971 2587
Gruenfrosch 0:7627c79db971 2588 {
Gruenfrosch 0:7627c79db971 2589 struct mad_bitptr peek;
Gruenfrosch 0:7627c79db971 2590 unsigned long header;
Gruenfrosch 0:7627c79db971 2591
Gruenfrosch 0:7627c79db971 2592 mad_bit_init(&peek, stream->next_frame);
Gruenfrosch 0:7627c79db971 2593
Gruenfrosch 0:7627c79db971 2594 header = mad_bit_read(&peek, 32);
Gruenfrosch 0:7627c79db971 2595 if ((header & 0xffe60000L) /* syncword | layer */ == 0xffe20000L) {
Gruenfrosch 0:7627c79db971 2596 if (!(header & 0x00010000L)) /* protection_bit */
Gruenfrosch 0:7627c79db971 2597 mad_bit_skip(&peek, 16); /* crc_check */
Gruenfrosch 0:7627c79db971 2598
Gruenfrosch 0:7627c79db971 2599 next_md_begin =
Gruenfrosch 0:7627c79db971 2600 mad_bit_read(&peek, (header & 0x00080000L) /* ID */ ? 9 : 8);
Gruenfrosch 0:7627c79db971 2601 }
Gruenfrosch 0:7627c79db971 2602
Gruenfrosch 0:7627c79db971 2603 mad_bit_finish(&peek);
Gruenfrosch 0:7627c79db971 2604 }
Gruenfrosch 0:7627c79db971 2605
Gruenfrosch 0:7627c79db971 2606 /* find main_data of this frame */
Gruenfrosch 0:7627c79db971 2607
Gruenfrosch 0:7627c79db971 2608 frame_space = stream->next_frame - mad_bit_nextbyte(&stream->ptr);
Gruenfrosch 0:7627c79db971 2609
Gruenfrosch 0:7627c79db971 2610 if (next_md_begin > si.main_data_begin + frame_space)
Gruenfrosch 0:7627c79db971 2611 next_md_begin = 0;
Gruenfrosch 0:7627c79db971 2612
Gruenfrosch 0:7627c79db971 2613 md_len = si.main_data_begin + frame_space - next_md_begin;
Gruenfrosch 0:7627c79db971 2614
Gruenfrosch 0:7627c79db971 2615 frame_used = 0;
Gruenfrosch 0:7627c79db971 2616
Gruenfrosch 0:7627c79db971 2617 if (si.main_data_begin == 0) {
Gruenfrosch 0:7627c79db971 2618 ptr = stream->ptr;
Gruenfrosch 0:7627c79db971 2619 stream->md_len = 0;
Gruenfrosch 0:7627c79db971 2620
Gruenfrosch 0:7627c79db971 2621 frame_used = md_len;
Gruenfrosch 0:7627c79db971 2622 }
Gruenfrosch 0:7627c79db971 2623 else {
Gruenfrosch 0:7627c79db971 2624 if (si.main_data_begin > stream->md_len) {
Gruenfrosch 0:7627c79db971 2625 if (result == 0) {
Gruenfrosch 0:7627c79db971 2626 stream->error = MAD_ERROR_BADDATAPTR;
Gruenfrosch 0:7627c79db971 2627 result = -1;
Gruenfrosch 0:7627c79db971 2628 }
Gruenfrosch 0:7627c79db971 2629 }
Gruenfrosch 0:7627c79db971 2630 else {
Gruenfrosch 0:7627c79db971 2631 mad_bit_init(&ptr,
Gruenfrosch 0:7627c79db971 2632 *stream->main_data + stream->md_len - si.main_data_begin);
Gruenfrosch 0:7627c79db971 2633
Gruenfrosch 0:7627c79db971 2634 if (md_len > si.main_data_begin) {
Gruenfrosch 0:7627c79db971 2635 assert(stream->md_len + md_len -
Gruenfrosch 0:7627c79db971 2636 si.main_data_begin <= MAD_BUFFER_MDLEN);
Gruenfrosch 0:7627c79db971 2637
Gruenfrosch 0:7627c79db971 2638 memcpy(*stream->main_data + stream->md_len,
Gruenfrosch 0:7627c79db971 2639 mad_bit_nextbyte(&stream->ptr),
Gruenfrosch 0:7627c79db971 2640 frame_used = md_len - si.main_data_begin);
Gruenfrosch 0:7627c79db971 2641 stream->md_len += frame_used;
Gruenfrosch 0:7627c79db971 2642 }
Gruenfrosch 0:7627c79db971 2643 }
Gruenfrosch 0:7627c79db971 2644 }
Gruenfrosch 0:7627c79db971 2645
Gruenfrosch 0:7627c79db971 2646 frame_free = frame_space - frame_used;
Gruenfrosch 0:7627c79db971 2647
Gruenfrosch 0:7627c79db971 2648 /* decode main_data */
Gruenfrosch 0:7627c79db971 2649
Gruenfrosch 0:7627c79db971 2650 if (result == 0) {
Gruenfrosch 0:7627c79db971 2651 error = III_decode(&ptr, frame, &si, nch);
Gruenfrosch 0:7627c79db971 2652 if (error) {
Gruenfrosch 0:7627c79db971 2653 stream->error = error;
Gruenfrosch 0:7627c79db971 2654 result = -1;
Gruenfrosch 0:7627c79db971 2655 }
Gruenfrosch 0:7627c79db971 2656
Gruenfrosch 0:7627c79db971 2657 /* designate ancillary bits */
Gruenfrosch 0:7627c79db971 2658
Gruenfrosch 0:7627c79db971 2659 stream->anc_ptr = ptr;
Gruenfrosch 0:7627c79db971 2660 stream->anc_bitlen = md_len * CHAR_BIT - data_bitlen;
Gruenfrosch 0:7627c79db971 2661 }
Gruenfrosch 0:7627c79db971 2662
Gruenfrosch 0:7627c79db971 2663 # if 0 && defined(DEBUG)
Gruenfrosch 0:7627c79db971 2664 fprintf(stderr,
Gruenfrosch 0:7627c79db971 2665 "main_data_begin:%u, md_len:%u, frame_free:%u, "
Gruenfrosch 0:7627c79db971 2666 "data_bitlen:%u, anc_bitlen: %u\n",
Gruenfrosch 0:7627c79db971 2667 si.main_data_begin, md_len, frame_free,
Gruenfrosch 0:7627c79db971 2668 data_bitlen, stream->anc_bitlen);
Gruenfrosch 0:7627c79db971 2669 # endif
Gruenfrosch 0:7627c79db971 2670
Gruenfrosch 0:7627c79db971 2671 /* preload main_data buffer with up to 511 bytes for next frame(s) */
Gruenfrosch 0:7627c79db971 2672
Gruenfrosch 0:7627c79db971 2673 if (frame_free >= next_md_begin) {
Gruenfrosch 0:7627c79db971 2674 memcpy(*stream->main_data,
Gruenfrosch 0:7627c79db971 2675 stream->next_frame - next_md_begin, next_md_begin);
Gruenfrosch 0:7627c79db971 2676 stream->md_len = next_md_begin;
Gruenfrosch 0:7627c79db971 2677 }
Gruenfrosch 0:7627c79db971 2678 else {
Gruenfrosch 0:7627c79db971 2679 if (md_len < si.main_data_begin) {
Gruenfrosch 0:7627c79db971 2680 unsigned int extra;
Gruenfrosch 0:7627c79db971 2681
Gruenfrosch 0:7627c79db971 2682 extra = si.main_data_begin - md_len;
Gruenfrosch 0:7627c79db971 2683 if (extra + frame_free > next_md_begin)
Gruenfrosch 0:7627c79db971 2684 extra = next_md_begin - frame_free;
Gruenfrosch 0:7627c79db971 2685
Gruenfrosch 0:7627c79db971 2686 if (extra < stream->md_len) {
Gruenfrosch 0:7627c79db971 2687 memmove(*stream->main_data,
Gruenfrosch 0:7627c79db971 2688 *stream->main_data + stream->md_len - extra, extra);
Gruenfrosch 0:7627c79db971 2689 stream->md_len = extra;
Gruenfrosch 0:7627c79db971 2690 }
Gruenfrosch 0:7627c79db971 2691 }
Gruenfrosch 0:7627c79db971 2692 else
Gruenfrosch 0:7627c79db971 2693 stream->md_len = 0;
Gruenfrosch 0:7627c79db971 2694
Gruenfrosch 0:7627c79db971 2695 memcpy(*stream->main_data + stream->md_len,
Gruenfrosch 0:7627c79db971 2696 stream->next_frame - frame_free, frame_free);
Gruenfrosch 0:7627c79db971 2697 stream->md_len += frame_free;
Gruenfrosch 0:7627c79db971 2698 }
Gruenfrosch 0:7627c79db971 2699
Gruenfrosch 0:7627c79db971 2700 return result;
Gruenfrosch 0:7627c79db971 2701 }