PageRenderTime 76ms CodeModel.GetById 15ms RepoModel.GetById 1ms app.codeStats 0ms

/AVS_Transcoder_SDK/kernel/audio/AudioTransCoder_Lame/av3enc/src/sam_encode_cbc.cpp

https://gitlab.com/jchwei/avs-transcoder
C++ | 1518 lines | 1162 code | 176 blank | 180 comment | 241 complexity | 1fb56c1f348652ddfd0166b9045a02a3 MD5 | raw file
    

  1. /*

    2. 

  2. ***********************************************************************

    3. 

  3. * COPYRIGHT AND WARRANTY INFORMATION

    4. 

  4. *

    5. 

  5. * Copyright 2004,  Audio Video Coding Standard, Part III

    6. 

  6. *

    7. 

  7. * This software module was originally developed by

    8. 

  8. *

    9. 

  9. * JungHoe Kim (kjh94@samsung.com), Samsung AIT

    10. 

  10. *

    11. 

  11. * edited by

    12. 

  12. *

    13. 

  13. * Lei Miao (win.miaolei@samsung.com), Samsung AIT

    14. 

  14. *

    15. 

  15. * Lei Miao, fill data embedded, 2005-03-16

    16. 

  16. *

    17. 

  17. * Lei Miao, coding band bug fixed, 2005-06-15

    18. 

  18. * Lei Miao, flexible frame length/MC_present, 2005-06-15

    19. 

  19. * Lei Miao, CBC Multi-channel extension, 2005-09-19

    20. 

  20. * Lei Miao, CBC codebook reduction, 2005-12-22

    21. 

  21. * Xun Li  , fill bits for buffer control, 2006-4-16

    22. 

  22. *

    23. 

  23. * DISCLAIMER OF WARRANTY

    24. 

  24. *

    25. 

  25. * These software programs are available to the users without any

    26. 

  26. * license fee or royalty on an "as is" basis. The AVS disclaims

    27. 

  27. * any and all warranties, whether express, implied, or statutory,

    28. 

  28. * including any implied warranties of merchantability or of fitness

    29. 

  29. * for a particular purpose. In no event shall the contributors or 

    30. 

  30. * the AVS be liable for any incidental, punitive, or consequential

    31. 

  31. * damages of any kind whatsoever arising from the use of this program.

    32. 

  32. *

    33. 

  33. * This disclaimer of warranty extends to the user of this program

    34. 

  34. * and user's customers, employees, agents, transferees, successors,

    35. 

  35. * and assigns.

    36. 

  36. *

    37. 

  37. * The AVS does not represent or warrant that the program furnished

    38. 

  38. * hereunder are free of infringement of any third-party patents.

    39. 

  39. * Commercial implementations of AVS, including shareware, may be

    40. 

  40. * subject to royalty fees to patent holders. Information regarding

    41. 

  41. * the AVS patent policy is available from the AVS Web site at

    42. 

  42. * http://www.avs.org.cn

    43. 

  43. *

    44. 

  44. * THIS IS NOT A GRANT OF PATENT RIGHTS - SEE THE AVS PATENT POLICY.

    45. 

  45. ************************************************************************

    46. 

  46. */

    47. 

  47. 

    48. 

  48. #include <sys/types.h>

    49. 

  49. #include <sys/stat.h>

    50. 

  50. #include <stdio.h>

    51. 

  51. #include <stdlib.h>

    52. 

  52. #include <math.h>

    53. 

  53. 

    54. 

  54. #include "av3enc.h"

    55. 

  55. 

    56. 

  56. extern void BSHCInit(FILE *name);

    57. 

  57. extern void BSHCClose(void);

    58. 

  58. extern int EncodeBSHCQuant(int model,int upper,int cur);

    59. 

  59. extern int EncodeBSHCSi(int siVal);

    60. 

  60. extern int EncodeBSHCSf(int sfVal);

    61. 

  61. extern int EncodeBSHCBin(int ms);

    62. 

  62. extern void EncodeBSHCStart(void);

    63. 

  63. extern void EncodeBSHCEnd(void);

    64. 

  64. extern void EncodeBSHCPutBits(int val,int len);

    65. 

  65. extern int EncodeBSHCGetSize();

    66. 

  66. extern void EncodeBSHCHeader(int ch,int freq);

    67. 

  67. extern int EncodeBSHCByteAlign(void);

    68. 

  68. extern int BSHCModelSelect(int bshcModel,int bpl);

    69. 

  69. extern void EncodeBSHCFlush(void);

    70. 

  70. extern int WriteAASFHeader(int ch,int freq,int bitrate,int aasf_flag);		//for AASF

    71. 

  71. 

    72. 

  72. 

    73. 

  73. extern void Register_buffer(void);

    74. 

  74. 

    75. 

  75. extern void Restore_buffer(void);

    76. 

  76. 

    77. 

  77. extern int EncodeBSHCGetAATFSize(void);

    78. 

  78. 

    79. 

  79. int encode_cband_si(

    80. 

  80. 					int	*model_index[][8],

    81. 

  81. 					int	start_cband[8],

    82. 

  82. 					int end_cband[8],

    83. 

  83. 					int g,

    84. 

  84. 					int nch);

    85. 

  85. 

    86. 

  86. int encode_scfband_si(

    87. 

  87. 					  int	*scf[][8],

    88. 

  88. 					  int	start_sfb[8],

    89. 

  89. 					  int end_sfb[8],

    90. 

  90. 					  int stereo_mode,

    91. 

  91. 					  int *stereo_info[8],

    92. 

  92. 					  int stereo_si_coded[8][MAX_SCFAC_BANDS],

    93. 

  93. 					  int g,

    94. 

  94. 					  int nch);

    95. 

  95. 

    96. 

  96. int encode_spectra(

    97. 

  97. 				   int	*sample[][8],

    98. 

  98. 				   int s_reg,

    99. 

  99. 				   int e_reg,

    100. 

  100. 				   int s_freq[8],

    101. 

  101. 				   int e_freq[8],

    102. 

  102. 				   int	min_bpl,

    103. 

  103. 				   int available_len,

    104. 

  104. 				   int *model_index[][8],

    105. 

  105. 				   int *cur_bpl[][8],

    106. 

  106. 				   int *coded_samp_bit[][8],

    107. 

  107. 				   int *sign_coded[][8],

    108. 

  108. 				   int *enc_vec[][8],

    109. 

  109. 				   int	nch);

    110. 

  110. 

    111. 

  111. #define min(x,y) ((x) > (y) ? (y) : (x))

    112. 

  112. #define max(x,y) ((x) > (y) ? (x) : (y))

    113. 

  113. 

    114. 

  114. double  fs_tbl[12]=

    115. 

  115. {96000, 88200, 64000, /* zhanjie 0708 */ 48000., 44100., 32000., 24000., 22050, 16000, 12000, 11025, 8000};

    116. 

  116. 

    117. 

  117. int channel_mapping_51[3][2]=

    118. 

  118. {{2, 0}, {4, 2}, {3, 3}};

    119. 

  119. 

    120. 

  120. static int sampling_rate;

    121. 

  121. 

    122. 

  122. void init_layer_variables (

    123. 

  123.     int signal_type,

    124. 

  124.     int num_window_groups,

    125. 

  125.     int window_group_length[],

    126. 

  126.     int maxSfb,

    127. 

  127.     int base_band,

    128. 

  128.     short	*swb_offset[],

    129. 

  129.     int top_layer,

    130. 

  130.     int *slayer_size,

    131. 

  131.     int layer_reg[],

    132. 

  132.     int layer_max_freq[],

    133. 

  133.     int layer_max_cband[],

    134. 

  134.     int layer_max_qband[]    

    135. 

  135. );

    136. 

  136. 

    137. 

  137. /* function prototypes */

    138. 

  138. void sam_scale_bits_init(int fsidx)

    139. 

  139. {

    140. 

  140.   sampling_rate = (int)fs_tbl[fsidx];

    141. 

  141. }

    142. 

  142. 

    143. 

  143. int get_base_band() 

    144. 

  144. {

    145. 

  145. 	int base_band;

    146. 

  146. 	switch(sampling_rate) {

    147. 

  147. 		case 96000:

    148. 

  148. 		case 88200:

    149. 

  149. 			base_band = 256;			//12000Hz zhanjie 0708

    150. 

  150. 			break;

    151. 

  151. 		case 64000:

    152. 

  152. 			base_band = 320;			//10000Hz zhanjie 0708

    153. 

  153. 			break;

    154. 

  154. 		case 48000:

    155. 

  155. 		case 44100:

    156. 

  156. 			base_band = 288 + 32+64;	//9000Hz

    157. 

  157. 			break;

    158. 

  158. 		case 32000:

    159. 

  159. 			base_band = 416;		// 6500 Hz

    160. 

  160. 			break;

    161. 

  161. 		case 24000:

    162. 

  162. 		case 22050:

    163. 

  163. 			base_band = 512;	//6000Hz

    164. 

  164. 			break;

    165. 

  165. 		case 16000:

    166. 

  166. 			base_band = 640;	//5000Hz

    167. 

  167. 			break;

    168. 

  168. 		case 12000:

    169. 

  169. 		case 11025:

    170. 

  170. 			base_band = 768;	//4500Hz

    171. 

  171. 			break;

    172. 

  172. 		case  8000:

    173. 

  173. 			base_band = 832;	//3250Hz

    174. 

  174. 			break;

    175. 

  175. 	}

    176. 

  176. 	return base_band;

    177. 

  177. }

    178. 

  178. 

    179. 

  179. int enc_wflag = 0;

    180. 

  180. 

    181. 

  181. int sam_encode_cbc(

    182. 

  182.     ChannelInfo channelInfo[],				   

    183. 

  183. 	int outputformat,

    184. 

  184. 	int target_bitrate,

    185. 

  185. 	int	signal_type,

    186. 

  186. 	int	windowShape,

    187. 

  187. 	int	*sample[][8],

    188. 

  188. 	int *scalefac[][8],

    189. 

  189. 	int	maxSfb,

    190. 

  190. 	int minSfb,

    191. 

  191. 	int	num_window_groups,

    192. 

  192. 	int window_group_length[],

    193. 

  193. 	short	*swb_offset[],

    194. 

  194. 	int	*model_index[][8],

    195. 

  195. 	int	stereo_mode,

    196. 

  196. 	int	*stereo_info[],

    197. 

  197. 	int abits,

    198. 

  198. 	int nch,

    199. 

  199. 	int wflag,

    200. 

  200. 	int isExtended,

    201. 

  201. 	int channelIdx,

    202. 

  202. 	int mc_present,

    203. 

  203. 	int extended_bytes,

    204. 

  204. 	int *fill_bits)

    205. 

  205. {

    206. 

  206. 	int		i, ch, g, m, k;

    207. 

  207. 	int     top_layer;

    208. 

  208. 	int		sfb, qband, cband;

    209. 

  209. 	int		base_snf;

    210. 

  210. 	int		base_band;

    211. 

  211. 	int		layer_reg[100];       /* Array Size : 64 -> 100  shpark 2000.04.20 */  

    212. 

  212. 	int		layer_max_freq[100];  /* Array Size : 64 -> 100  shpark 2000.04.20 */  

    213. 

  213. 	int		layer_max_cband[100]; /* Array Size : 64 -> 100  shpark 2000.04.20 */  

    214. 

  214. 	int		layer_max_qband[100]; /* Array Size : 64 -> 100  shpark 2000.04.20 */  

    215. 

  215. 	//int		layer_bit_flush[100];    /* Array Size : 64 -> 100  shpark 2000.04.20 */  

    216. 

  216. 	int		layer_buf_offset[100];   /* Array Size : 64 -> 100  shpark 2000.04.20 */  

    217. 

  217. 	int		layer_si_maxlen[100];    /* Array Size : 64 -> 100  shpark 2000.04.20 */  

    218. 

  218. 	int		layer_extra_len[100];    /* Array Size : 64 -> 100  shpark 2000.04.20 */  

    219. 

  219. 	int		layer_cw_offset[100];    /* Array Size : 64 -> 100  shpark 2000.04.20 */  

    220. 

  220. 

    221. 

  221. 	int		start_freq[8];    /* Array Size : 64 -> 8  shpark 2000.04.20 */

    222. 

  222. 	int		start_qband[8];   /* Array Size : 64 -> 8  shpark 2000.04.20 */

    223. 

  223. 	int		start_cband[8];   /* Array Size : 64 -> 8  shpark 2000.04.20 */

    224. 

  224. 	int		end_freq[8];      /* Array Size : 64 -> 8  shpark 2000.04.20 */

    225. 

  225. 	int		end_cband[8];     /* Array Size : 64 -> 8  shpark 2000.04.20 */

    226. 

  226. 	int		end_qband[8];     /* Array Size : 64 -> 8  shpark 2000.04.20 */

    227. 

  227. 	int		s_freq[8];        /* Array Size : 64 -> 8  shpark 2000.04.20 */

    228. 

  228. 	int		e_freq[8];        /* Array Size : 64 -> 8  shpark 2000.04.20 */

    229. 

  229. 

    230. 

  230. 	int		cw_len;

    231. 

  231. 	int		layer;

    232. 

  232. 	int		slayer_size;

    233. 

  233. 	int		si_offset;

    234. 

  234. 	//int		est_total_len;

    235. 

  235. 	int		available_len;

    236. 

  236. 	//int		used_bits;

    237. 

  237. 	int		cur_snf_buf[2][1024];

    238. 

  238. 	int		*cur_snf[2][8];

    239. 

  239. 	int		band_snf[2][8][32];

    240. 

  240. 	int		sign_is_coded_buf[2][1024];

    241. 

  241. 	int		*sign_is_coded[2][8];

    242. 

  242. 	int		*coded_samp_bit[2][8];

    243. 

  243. 	int		coded_samp_bit_buf[2][1024];

    244. 

  244. 	int		*enc_vec[2][8];

    245. 

  245. 	int		enc_vec_buf[2][256];

    246. 

  246.     int     stereo_si_coded[8][MAX_SCFAC_BANDS];

    247. 

  247. 	int		code_len;

    248. 

  248. 	

    249. 

  249. 	int layerBit;

    250. 

  250. 	int headerLen;

    251. 

  251. 	int specBits=0;

    252. 

  252. 

    253. 

  253.     //for AATF 

    254. 

  254. 	int bit_count0;

    255. 

  255. 	int	bit_count1;

    256. 

  256. 

    257. 

  257. 	/* fill element variables */

    258. 

  258. 	int fill_enable=1;

    259. 

  259. 	int fill_length=0;

    260. 

  260. 	int fill_ele_number=0;	

    261. 

  261. 	int fill_cnt[128];

    262. 

  262. 	/* default fill */

    263. 

  263. 	int fill_type=FILL_DFT;

    264. 

  264. 	

    265. 

  265. 	int enhance_FLCnt;

    266. 

  266. 	int FL_fieldBits;

    267. 

  267. 	int enhance_enable;

    268. 

  268. 	int FL_Byte;

    269. 

  269. 	int FL_Bit;

    270. 

  270. 	int headerAlignBits, headerAlignBytes;

    271. 

  271. 

    272. 

  272. 	int channel_configuration_index;

    273. 

  273. 	int test=0;

    274. 

  274. 

    275. 

  275. 	// 20060116 Miao

    276. 

  276. 	int scalable_header_reuse=0;

    277. 

  277. 

    278. 

  278. 	enc_wflag = wflag; 

    279. 

  279. 	

    280. 

  280. 	// 2006-4-16 xun

    281. 

  281. 	*fill_bits = ((*fill_bits + 7) >> 3) << 3;

    282. 

  282. 

    283. 

  283. 	/* ***************************************************** */

    284. 

  284. 	/*                  Initialize variables                 */

    285. 

  285. 	/* ***************************************************** */

    286. 

  286. 	for(ch = 0; ch < nch; ch++) {

    287. 

  287. 		int s;

    288. 

  288. 

    289. 

  289. 		s = 0;

    290. 

  290. 		for(g = 0; g < num_window_groups; g++) { 

    291. 

  291. 			sign_is_coded[ch][g] = &(sign_is_coded_buf[ch][1024*s/8]);

    292. 

  292. 			coded_samp_bit[ch][g] = &(coded_samp_bit_buf[ch][1024*s/8]);

    293. 

  293. 			enc_vec[ch][g] = &(enc_vec_buf[ch][(1024*s/8)/4]);

    294. 

  294. 			cur_snf[ch][g] = &(cur_snf_buf	[ch][1024*s/8]);

    295. 

  295. 			s += window_group_length[g];

    296. 

  296. 

    297. 

  297.            for (i=0; i<MAX_SCFAC_BANDS; i++)

    298. 

  298.               stereo_si_coded[g][i] = 0;

    299. 

  299. 		}

    300. 

  300. 

    301. 

  301. 		for (i=0; i<1024; i++) {

    302. 

  302. 			sign_is_coded_buf[ch][i] = 0;

    303. 

  303. 			coded_samp_bit_buf[ch][i] = 0;

    304. 

  304. 			cur_snf_buf[ch][i] = 0;

    305. 

  305. 		}

    306. 

  306. 		for (i=0; i<256; i++) {

    307. 

  307. 			enc_vec_buf[ch][i] = 0;

    308. 

  308. 		}

    309. 

  309. 	}

    310. 

  310. 

    311. 

  311. 	top_layer = ((target_bitrate/nch) / 1000) - 16;

    312. 

  312. 	if (top_layer > 63) top_layer = 63; /* top_layer limit inserted by shpark 2000. 09. 25 */

    313. 

  313. 

    314. 

  314. 

    315. 

  315. 	/******************************/

    316. 

  316. 	/* Adjust necessary variables */

    317. 

  317. 	/******************************/

    318. 

  318. 	base_snf = 1;

    319. 

  319.     base_band = get_base_band();

    320. 

  320. 

    321. 

  321. 	/* ************************************************

    322. 

  322. 		Initialize

    323. 

  323. 			layer_max_freq

    324. 

  324. 			layer_max_cband

    325. 

  325. 			layer_max_qband

    326. 

  326. 			layer_bit_flush

    327. 

  327. 	************************************************ */

    328. 

  328. 

    329. 

  329.    init_layer_variables (signal_type, num_window_groups, 

    330. 

  330. 	   window_group_length, maxSfb,

    331. 

  331.        base_band, swb_offset, top_layer, &slayer_size, layer_reg,

    332. 

  332.        layer_max_freq,  layer_max_cband, layer_max_qband);

    333. 

  333. 

    334. 

  334.    /* maxSfb limitation */

    335. 

  335.    if (layer_max_qband[top_layer+slayer_size] < maxSfb)

    336. 

  336. 		maxSfb = layer_max_qband[top_layer+slayer_size-1];

    337. 

  337. 

    338. 

  338.    for(g = 0; g < num_window_groups; g++) 

    339. 

  339. 	  end_cband[g] = (swb_offset[g][maxSfb]+31)  / 32;

    340. 

  340. 

    341. 

  341.    if(signal_type == TRANSIENT_TYPE) {

    342. 

  342. 	  for (layer = 0; layer < slayer_size; layer++) {

    343. 

  343. 	 	 g = layer_reg[layer];

    344. 

  344. 		 end_qband[g] = layer_max_qband[layer];

    345. 

  345. 	  }

    346. 

  346.    }

    347. 

  347.    else

    348. 

  348. 	  end_qband[0] = layer_max_qband[slayer_size-1];

    349. 

  349. 

    350. 

  350. 	/* ***************************************************** */

    351. 

  351. 	/*                  syntax of cbc_header()              */

    352. 

  352. 	/* ***************************************************** */

    353. 

  353. 	

    354. 

  354. 	//HUFF start

    355. 

  355.    if(outputformat == 2)		

    356. 

  356. 		Register_buffer();

    357. 

  357.    else

    358. 

  358. 		EncodeBSHCStart();

    359. 

  359. 

    360. 

  360.    bit_count0 = EncodeBSHCGetSize();

    361. 

  361. 

    362. 

  362. 	/* frame length */

    363. 

  363. 	enhance_enable = 0;

    364. 

  364. 	/* base frame length */

    365. 

  365. 	if(abits == 20000)  // lixun 2008-7-28, max bits 6144 * 2=12288, so use 20000 as esmitated bits number

    366. 

  366. 	{

    367. 

  367. 		EncodeBSHCPutBits(0,8);

    368. 

  368. 	}

    369. 

  369. 	else

    370. 

  370. 	{	

    371. 

  371. 		enhance_enable = 1;

    372. 

  372. 		enhance_FLCnt = 0;

    373. 

  373. 		FL_Byte = abits >> 3;

    374. 

  374. 		FL_Byte += extended_bytes;

    375. 

  375. 		//bug fix: 20006-4-16 Xun

    376. 

  376. 		if (!isExtended)

    377. 

  377. 			FL_Byte += (*fill_bits>>3);

    378. 

  378. 

    379. 

  379. 		/* 4 bytes: minimum header length */

    380. 

  380. 		while( (FL_Byte-4) > ((1<<(7+enhance_FLCnt*3))-1) )

    381. 

  381. 			enhance_FLCnt ++;

    382. 

  382. 		FL_fieldBits = 7 + enhance_FLCnt*3;

    383. 

  383. 		/* base frame length */

    384. 

  384. 		EncodeBSHCPutBits((FL_Byte-4)>>(enhance_FLCnt*3),7);

    385. 

  385. 

    386. 

  386. 		while(enhance_FLCnt) {

    387. 

  387. 			EncodeBSHCPutBits(1, 1);

    388. 

  388. 			EncodeBSHCPutBits(((FL_Byte-4)>>((enhance_FLCnt-1)*3))&0x7,3);

    389. 

  389. 			enhance_FLCnt --;

    390. 

  390. 		}

    391. 

  391. 		/* frame length end flag */

    392. 

  392. 		EncodeBSHCPutBits(0, 1);

    393. 

  393. 	}

    394. 

  394. 

    395. 

  395. 	/* MC_present */

    396. 

  396. 	if(!isExtended) {

    397. 

  397. 		EncodeBSHCPutBits(mc_present,1);

    398. 

  398. 		if(mc_present) {

    399. 

  399. 			/* FL/FR length */

    400. 

  400. 			/* base frame length */

    401. 

  401. 			if(abits == 20000)

    402. 

  402. 			{

    403. 

  403. 				EncodeBSHCPutBits(0,8);

    404. 

  404. 			}

    405. 

  405. 			else {

    406. 

  406. 				enhance_enable = 1;

    407. 

  407. 				enhance_FLCnt = 0;

    408. 

  408. 				FL_Byte = abits >> 3;

    409. 

  409. 				/* 4 bytes: minimum header length */

    410. 

  410. 				while( (FL_Byte-4) > ((1<<(7+enhance_FLCnt*3))-1) )

    411. 

  411. 					enhance_FLCnt ++;

    412. 

  412. 				FL_fieldBits = 7 + enhance_FLCnt*3;

    413. 

  413. 				/* base frame length */

    414. 

  414. 				EncodeBSHCPutBits((FL_Byte-4)>>(enhance_FLCnt*3),7);

    415. 

  415. 				while(enhance_FLCnt) {

    416. 

  416. 					EncodeBSHCPutBits(1, 1);

    417. 

  417. 					EncodeBSHCPutBits(((FL_Byte-4)>>((enhance_FLCnt-1)*3))&0x7,3);

    418. 

  418. 					enhance_FLCnt --;

    419. 

  419. 				}

    420. 

  420. 				/* frame length end flag */

    421. 

  421. 				EncodeBSHCPutBits(0, 1);

    422. 

  422. 			}

    423. 

  423. 		}

    424. 

  424. 	}

    425. 

  425. 	else {

    426. 

  426. 		ch = 0;

    427. 

  427. 		while(channel_mapping_51[ch][0] != channelIdx)

    428. 

  428. 			ch ++;

    429. 

  429. 

    430. 

  430. 		channel_configuration_index = channel_mapping_51[ch][1];

    431. 

  431. 

    432. 

  432. 		/* channel configuration index: 4 */

    433. 

  433. 		EncodeBSHCPutBits(channel_configuration_index,4);

    434. 

  434. 

    435. 

  435. 		// 20060116 Miao

    436. 

  436. 		/* scalable_header_reuse: 1 bit */

    437. 

  437. 		EncodeBSHCPutBits(scalable_header_reuse, 1);

    438. 

  438. 

    439. 

  439. 		/* reserved_bit : 1 */

    440. 

  440. 		EncodeBSHCPutBits(0, 1);

    441. 

  441. 	}

    442. 

  442. 	

    443. 

  443.    	/* top layer index */

    444. 

  444. 	EncodeBSHCPutBits(top_layer,6);

    445. 

  445. 

    446. 

  446.    	/* base snf */

    447. 

  447. 	EncodeBSHCPutBits(base_snf-1,2);

    448. 

  448. 

    449. 

  449.    	/* base band */

    450. 

  450. 	EncodeBSHCPutBits(base_band/32,5);

    451. 

  451. 

    452. 

  452. 	/* ***************************************************** */

    453. 

  453. 	/*               syntax of general_header()              */

    454. 

  454. 	/* ***************************************************** */

    455. 

  455. 

    456. 

  456. 	/* reserved_bit : 1 */

    457. 

  457. 	EncodeBSHCPutBits(0,1);

    458. 

  458. 

    459. 

  459. 	for(ch = 0; ch < nch; ch++) {

    460. 

  460. 		EncodeBSHCPutBits(scalefac[ch][0][0],8);

    461. 

  461. 	}

    462. 

  462. 

    463. 

  463. 	EncodeBSHCPutBits((int)signal_type,1);

    464. 

  464. 

    465. 

  465. 	EncodeBSHCPutBits((int)windowShape,1);

    466. 

  466. 

    467. 

  467. 	if(signal_type == TRANSIENT_TYPE) {

    468. 

  468. 		extern int siCodeMode;

    469. 

  469. 		EncodeBSHCPutBits(maxSfb,4);

    470. 

  470. 

    471. 

  471. 		for (i=1; i<window_group_length[0]; i++)

    472. 

  472. 		{

    473. 

  473. 			EncodeBSHCPutBits(1,1);

    474. 

  474. 		}

    475. 

  475. 		for(g = 1; g < num_window_groups; g++) {

    476. 

  476. 			EncodeBSHCPutBits(0,1);

    477. 

  477. 			for (i=1; i<window_group_length[g]; i++)

    478. 

  478. 			{

    479. 

  479. 				EncodeBSHCPutBits(1,1);

    480. 

  480. 			}			

    481. 

  481. 		}

    482. 

  482. 		EncodeBSHCPutBits(siCodeMode,1);			

    483. 

  483. 		

    484. 

  484. 	} else {

    485. 

  485. 		EncodeBSHCPutBits(maxSfb,6);

    486. 

  486. 	}

    487. 

  487. 	

    488. 

  488. 	if(nch == 2) {

    489. 

  489. 		EncodeBSHCPutBits(stereo_mode,2);

    490. 

  490. 	}

    491. 

  491. 	

    492. 

  492. 	/* fill enable, 1 bit */

    493. 

  493. 	// xun 2006-4-16 fill bits

    494. 

  494. 	if (!isExtended) {

    495. 

  495. 		if (*fill_bits == 0)

    496. 

  496. 			fill_enable = 0;

    497. 

  497. 		fill_length = *fill_bits >> 3;

    498. 

  498. 		

    499. 

  499. 		EncodeBSHCPutBits(fill_enable,1);

    500. 

  500. 

    501. 

  501. 		if (*fill_bits > 0)

    502. 

  502. 			fill_ele_number = (*fill_bits + MAX_FILL_BITS - 1)/MAX_FILL_BITS - 1;

    503. 

  503. 		else

    504. 

  504. 			fill_ele_number = 0;

    505. 

  505. 		

    506. 

  506. 		if(fill_enable)

    507. 

  507. 		{			

    508. 

  508. 			int tmp = *fill_bits;

    509. 

  509. 

    510. 

  510. 			EncodeBSHCPutBits(fill_ele_number, 7);

    511. 

  511. 			

    512. 

  512. 			for (i = 0; i < fill_ele_number+1; i++)

    513. 

  513. 			{

    514. 

  514. 				fill_cnt[i] = min(tmp, MAX_FILL_BITS);

    515. 

  515. 				tmp -= fill_cnt[i];

    516. 

  516. 				fill_cnt[i] >>= 3;

    517. 

  517. 								

    518. 

  518. 				/* fill length field */

    519. 

  519. 				if(fill_cnt[i] < 15)

    520. 

  520. 				{

    521. 

  521. 					EncodeBSHCPutBits(fill_cnt[i],4);

    522. 

  522. 				}

    523. 

  523. 				else

    524. 

  524. 				{

    525. 

  525. 					EncodeBSHCPutBits(15,4);

    526. 

  526. 					EncodeBSHCPutBits(fill_cnt[i]-15,8);

    527. 

  527. 				}

    528. 

  528. 			}

    529. 

  529. 		}

    530. 

  530. 	}

    531. 

  531. 

    532. 

  532. 	/* ***************************************************** */

    533. 

  533. 	/*                  extension_coding_element()           */

    534. 

  534. 	/* ***************************************************** */

    535. 

  535. 	{		

    536. 

  536. 		int chIndex = 0;

    537. 

  537. 		FLPVQInfo *flpVQInfo;

    538. 

  538. 		chIndex  = channelIdx;			

    539. 

  539. 		for (ch = 0; ch < nch; ch++) {	

    540. 

  540. 			flpVQInfo = &channelInfo[chIndex + ch].flpInfo; // wlei modified [2005-12-22]

    541. 

  541. 			EncodeBSHCPutBits(flpVQInfo->nFlpPresent, 1);			

    542. 

  542. 			if (flpVQInfo->nFlpPresent) {

    543. 

  543. 				EncodeBSHCPutBits(flpVQInfo->nPredDirection, 1);

    544. 

  544. 				EncodeBSHCPutBits(flpVQInfo->nVQCode1Index, 10);

    545. 

  545. 				EncodeBSHCPutBits(flpVQInfo->nVQCode2Index, 10);

    546. 

  546. 				EncodeBSHCPutBits(flpVQInfo->nVQCode3Index, 9);

    547. 

  547. 				EncodeBSHCPutBits(flpVQInfo->nVQCode4Index, 8);

    548. 

  548. 				chIndex ++;

    549. 

  549. 			}			

    550. 

  550. 		}	

    551. 

  551. 	}

    552. 

  552. 

    553. 

  553. 	

    554. 

  554. 

    555. 

  555. 	if(!enhance_enable) {

    556. 

  556. 		headerAlignBits = EncodeBSHCGetSize() & 0x7;

    557. 

  557. 		if(headerAlignBits)

    558. 

  558. 			headerAlignBits = 8 - headerAlignBits;

    559. 

  559. 	}

    560. 

  560. 

    561. 

  561. 	/*********** byte align *********/

    562. 

  562.  	EncodeBSHCByteAlign();

    563. 

  563. 

    564. 

  564. 	if (outputformat == 2) 

    565. 

  565. 		headerLen = EncodeBSHCGetAATFSize();

    566. 

  566. 	else

    567. 

  567. 		headerLen = EncodeBSHCGetSize();

    568. 

  568. 

    569. 

  569. 

    570. 

  570. 	if(maxSfb == 0) {

    571. 

  571. 		if(nch == 1) return headerLen;

    572. 

  572. 		if(nch == 2 && maxSfb == 0) return headerLen;

    573. 

  573. 	}

    574. 

  574. 	

    575. 

  575. 	for(ch = 0; ch < nch; ch++) {

    576. 

  576. 		for(g = 0; g < num_window_groups; g++) {

    577. 

  577. 			for(cband = 0; cband < end_cband[g]; cband++) {

    578. 

  578. 				if(model_index[ch][g][cband]==0)

    579. 

  579. 				{

    580. 

  580. 					band_snf[ch][g][cband]=0;

    581. 

  581. 					continue;

    582. 

  582. 				}

    583. 

  583. 				if(model_index[ch][g][cband] >= 9)

    584. 

  584. 					band_snf[ch][g][cband] = model_index[ch][g][cband] -4;

    585. 

  585. 				else

    586. 

  586. 				{

    587. 

  587. 					band_snf[ch][g][cband] = (model_index[ch][g][cband]+1)/2;

    588. 

  588. 				}

    589. 

  589. 			}

    590. 

  590. 		}

    591. 

  591. 	}

    592. 

  592. 

    593. 

  593. 	/* ##################################################### */

    594. 

  594. 	/* 	                CBC MAIN ROUTINE                     */

    595. 

  595. 	/* ##################################################### */

    596. 

  596. 

    597. 

  597. 	/* ***************************************************** */

    598. 

  598. 	/*                  CBC_layer_stream()                   */

    599. 

  599. 	/* ***************************************************** */

    600. 

  600. 	for(g = 0; g < num_window_groups; g++)

    601. 

  601. 		end_cband[g] = end_qband[g] = 0;

    602. 

  602. 

    603. 

  603. 	for (layer=0; layer<(top_layer+slayer_size); layer++) {

    604. 

  604. 		g = layer_reg[layer];

    605. 

  605.  	 	start_qband[g] = end_qband[g];

    606. 

  606. 		start_cband[g] = end_cband[g];

    607. 

  607. 		end_qband[g] = layer_max_qband[layer];

    608. 

  608. 		end_cband[g] = layer_max_cband[layer];

    609. 

  609. 	}

    610. 

  610. 

    611. 

  611. 	/* to calculate core CBC data size */

    612. 

  612. 	layer_buf_offset[0] = (double)16.*(abits-headerLen)

    613. 

  613. 		/(double)(top_layer+16)/(double)slayer_size;

    614. 

  614. 	for (layer = 1; layer < slayer_size-1; layer++)

    615. 

  615. 	{

    616. 

  616. 		layer_buf_offset[layer] = layer_buf_offset[layer-1] + 

    617. 

  617. 			(double)16.*(abits-headerLen)/(top_layer+16)/slayer_size;

    618. 

  618. 	}

    619. 

  619. 	layer_buf_offset[slayer_size-1] = 

    620. 

  620. 		(double)16.*(abits-headerLen)/(top_layer+16);

    621. 

  621. 	layerBit = (abits-headerLen) - layer_buf_offset[slayer_size-1];

    622. 

  622. 	for (layer = slayer_size; layer < top_layer+slayer_size-1; layer++)

    623. 

  623. 	{

    624. 

  624. 		layer_buf_offset[layer] = layer_buf_offset[layer-1] +

    625. 

  625. 			layerBit/(top_layer);

    626. 

  626. 	}

    627. 

  627. 	layer_buf_offset[top_layer+slayer_size-1] = (abits-headerLen);

    628. 

  628. 

    629. 

  629. 	for (layer = top_layer+slayer_size; layer < 100; layer++)

    630. 

  630. 		layer_buf_offset[layer] = (abits-headerLen);

    631. 

  631. 

    632. 

  632. 	for(i = 0; i < 100; i++) {

    633. 

  633. 		layer_extra_len[i] = 0;

    634. 

  634. 		layer_cw_offset[i] = 0;

    635. 

  635. 	}

    636. 

  636. 

    637. 

  637. 	//est_total_len = used_bits;

    638. 

  638. 

    639. 

  639. 	for(g = 0; g < num_window_groups; g++) {

    640. 

  640. 	   	start_qband[g] = 0;

    641. 

  641. 		start_cband[g] = 0;

    642. 

  642. 		start_freq[g] = 0;

    643. 

  643. 	   	end_qband[g] = 0;

    644. 

  644. 		end_cband[g] = 0;

    645. 

  645. 		end_freq[g] = 0;

    646. 

  646. 	}

    647. 

  647. 

    648. 

  648. 	for(g = 0; g < num_window_groups; g++)

    649. 

  649. 		s_freq[g] = 0;

    650. 

  650. 

    651. 

  651. 	if (signal_type == TRANSIENT_TYPE) {

    652. 

  652. 		for (layer = 0; layer < slayer_size; layer++) {

    653. 

  653. 			g = layer_reg[layer];

    654. 

  654. 			e_freq[g]  = layer_max_freq[layer];

    655. 

  655. 		}

    656. 

  656. 	}

    657. 

  657. 

    658. 

  658. 	code_len = 0;

    659. 

  659. 	

    660. 

  660. 	available_len = 0;

    661. 

  661. 	/* ##################################################### */

    662. 

  662. 	/* 	                CBC MAIN LOOP                        */

    663. 

  663. 	/* ##################################################### */

    664. 

  664. 	for (layer = 0; layer < top_layer+slayer_size; layer++) {

    665. 

  665. 		int min_snf;

    666. 

  666. 		

    667. 

  667. 		g = layer_reg[layer];

    668. 

  668. 

    669. 

  669. 		start_freq[g]  = end_freq[g];

    670. 

  670. 	   	start_qband[g] = end_qband[g];

    671. 

  671. 		start_cband[g] = end_cband[g];

    672. 

  672. 

    673. 

  673. 		end_qband[g] = layer_max_qband[layer];

    674. 

  674. 		end_cband[g] = layer_max_cband[layer];

    675. 

  675. 		end_freq[g]  = layer_max_freq[layer];

    676. 

  676. 

    677. 

  677. 		for(ch = 0; ch < nch; ch++) {

    678. 

  678. 			for(i = start_freq[g]; i < end_freq[g]; i++)  

    679. 

  679. 				cur_snf[ch][g][i] = band_snf[ch][g][i/32];

    680. 

  680. 		}

    681. 

  681. 

    682. 

  682. 		if (signal_type == TRANSIENT_TYPE) {

    683. 

  683. 			if(layer >= slayer_size) {

    684. 

  684. 				g = layer_reg[layer];

    685. 

  685. 				e_freq[g]  = layer_max_freq[layer];

    686. 

  686. 			}

    687. 

  687. 		} else {

    688. 

  688. 			if (layer >= slayer_size) 

    689. 

  689. 				e_freq[0] = layer_max_freq[layer];

    690. 

  690. 			else

    691. 

  691.    				e_freq[0] = layer_max_freq[slayer_size-1];

    692. 

  692. 		}

    693. 

  693. 		if(layer==0)

    694. 

  694. 			available_len = layer_buf_offset[0];

    695. 

  695. 		else

    696. 

  696. 			available_len += layer_buf_offset[layer]-

    697. 

  697. 			layer_buf_offset[layer-1];

    698. 

  698. 		

    699. 

  699. 		//est_total_len += available_len;

    700. 

  700. 

    701. 

  701. 		cw_len =  encode_cband_si(model_index, start_cband, end_cband, g, nch);

    702. 

  702. 		code_len += cw_len;

    703. 

  703. 		

    704. 

  704. 		available_len -= cw_len;

    705. 

  705. 		g = layer_reg[layer];

    706. 

  706. 		

    707. 

  707. 		/* side infomation : scalefactor */

    708. 

  708. 		cw_len =  encode_scfband_si(scalefac, start_qband, 

    709. 

  709. 			end_qband, stereo_mode, stereo_info,

    710. 

  710. 						stereo_si_coded, g, nch);

    711. 

  711. 

    712. 

  712. 		code_len += cw_len;

    713. 

  713. 		available_len -= cw_len;

    714. 

  714. 		

    715. 

  715. 		min_snf = layer < slayer_size ? base_snf : 1;

    716. 

  716. 

    717. 

  717. 		/* Bit-Sliced Spectral Data Coding */

    718. 

  718. 		cw_len = encode_spectra(sample, g, g+1, start_freq, end_freq, 

    719. 

  719. 				min_snf, available_len, model_index, cur_snf, coded_samp_bit, 

    720. 

  720. 				sign_is_coded, enc_vec, nch);

    721. 

  721. 

    722. 

  722. 		test += cw_len;

    723. 

  723. 		

    724. 

  724. 		code_len += cw_len;

    725. 

  725. 		available_len -= cw_len;

    726. 

  726. 		specBits+= cw_len;

    727. 

  727. 

    728. 

  728. 		if(available_len > 0) {

    729. 

  729. 			cw_len = encode_spectra(sample, 0, num_window_groups, s_freq, e_freq,

    730. 

  730. 				1, available_len, model_index, cur_snf, coded_samp_bit, 

    731. 

  731. 				sign_is_coded, enc_vec, nch);

    732. 

  732. 			code_len += cw_len;

    733. 

  733. 			available_len -= cw_len;

    734. 

  734. 			specBits+= cw_len;

    735. 

  735. 

    736. 

  736. 			test += cw_len;

    737. 

  737. 		}

    738. 

  738. 		//est_total_len -= available_len;

    739. 

  739. 

    740. 

  740. 	}

    741. 

  741. 

    742. 

  742. 	if(!enhance_enable) {

    743. 

  743. 		int byte_incr = 0;

    744. 

  744. 		int FL_Bit_cur = 0;

    745. 

  745. 		int enhance_FLCnt_last = 0;

    746. 

  746. 		// 20060116 Miao

    747. 

  747. 		if(outputformat != 2)

    748. 

  748. 			FL_Bit = EncodeBSHCGetSize();

    749. 

  749. 		else

    750. 

  750. 			FL_Bit = EncodeBSHCGetSize() - bit_count0;

    751. 

  751. 

    752. 

  752. 		FL_Byte = (FL_Bit+7) >> 3;

    753. 

  753. 

    754. 

  754. 		FL_Bit_cur = FL_Bit;

    755. 

  755. 		while (1)

    756. 

  756. 		{

    757. 

  757. 			enhance_FLCnt = 0;

    758. 

  758. 			FL_Byte = (FL_Bit_cur+7) >> 3;

    759. 

  759. 			// 2006-4-16 Xun bug fixed

    760. 

  760. 			if (!isExtended)

    761. 

  761. 			{

    762. 

  762. 				FL_Byte += extended_bytes;

    763. 

  763. 				FL_Byte += fill_enable * fill_length;

    764. 

  764. 			}

    765. 

  765. 			/* 4 bytes: minimum header length */

    766. 

  766. 			while((FL_Byte-4) > ((1<<(7+enhance_FLCnt*3))-1))

    767. 

  767. 				enhance_FLCnt ++;

    768. 

  768. 			

    769. 

  769. 			// 2006-4-16 Xun

    770. 

  770. 			if(!isExtended && mc_present) {

    771. 

  771. 				int enhance_FLCnt2;

    772. 

  772. 				FL_Byte -= extended_bytes;

    773. 

  773. 				FL_Byte -= fill_enable * fill_length;

    774. 

  774. 

    775. 

  775. 				enhance_FLCnt2 = 0;

    776. 

  776. 				while((FL_Byte-4) > ((1<<(7+enhance_FLCnt2*3))-1))

    777. 

  777. 					enhance_FLCnt2 ++;

    778. 

  778. 				enhance_FLCnt += enhance_FLCnt2;

    779. 

  779. 			}	

    780. 

  780. 			if(enhance_FLCnt > enhance_FLCnt_last) {

    781. 

  781. 				if(headerAlignBits < (enhance_FLCnt<<2)) {

    782. 

  782. 					byte_incr = (enhance_FLCnt<<2) - headerAlignBits;					

    783. 

  783. 					byte_incr = (byte_incr+7)>>3;

    784. 

  784. 					// 2006-4-16 xun

    785. 

  785. 					FL_Bit_cur = FL_Bit + (byte_incr<<3);

    786. 

  786. 				}

    787. 

  787. 				else

    788. 

  788. 					break;

    789. 

  789. 			}

    790. 

  790. 			else

    791. 

  791. 				break;

    792. 

  792. 			enhance_FLCnt_last = enhance_FLCnt;

    793. 

  793. 		}

    794. 

  794. 		FL_Bit = FL_Bit_cur;

    795. 

  795. 		FL_Byte = (FL_Bit+7) >> 3;	

    796. 

  796. 		// 20060116 li xun

    797. 

  797. 		for(i = 0; i < byte_incr; i ++)

    798. 

  798. 			EncodeBSHCPutBits(0, 8);

    799. 

  799. 	}

    800. 

  800. 

    801. 

  801. 	EncodeBSHCByteAlign();

    802. 

  802. 

    803. 

  803. 	/* fill data */

    804. 

  804. 	// xun 2006-4-16 fill buffer

    805. 

  805. 	if(!isExtended && fill_enable)

    806. 

  806. 	{

    807. 

  807. 		for (k = 0; k < fill_ele_number + 1; k++)

    808. 

  808. 		{		

    809. 

  809. 			EncodeBSHCPutBits(fill_type,4);

    810. 

  810. 			switch(fill_type)

    811. 

  811. 			{			

    812. 

  812. 			case FILL_DFT:

    813. 

  813. 				for(i = 0; i < 8*(fill_cnt[k]-1)+4; i ++)

    814. 

  814. 					EncodeBSHCPutBits(1,1);

    815. 

  815. 				break;

    816. 

  816. 			default: 

    817. 

  817. 				;

    818. 

  818. 			}

    819. 

  819. 		}

    820. 

  820. 	}

    821. 

  821. 

    822. 

  822. 	//EncodeBSHCEnd();

    823. 

  823. 	// 20060116 Miao

    824. 

  824. 	bit_count1 = EncodeBSHCGetSize()-*fill_bits;

    825. 

  825. 

    826. 

  826. 	if(outputformat == 2 && !wflag)

    827. 

  827. 		Restore_buffer();		

    828. 

  828.  	return (bit_count1 - bit_count0);

    829. 

  829. 

    830. 

  830. }

    831. 

  831. 

    832. 

  832. void init_layer_variables (

    833. 

  833. int signal_type,

    834. 

  834. int num_window_groups,

    835. 

  835. int window_group_length[],

    836. 

  836. int maxSfb,

    837. 

  837. int base_band,

    838. 

  838. short	*swb_offset[],

    839. 

  839. int top_layer,

    840. 

  840. int *slayer_size,

    841. 

  841. int layer_reg[],

    842. 

  842. int layer_max_freq[],

    843. 

  843. int layer_max_cband[],

    844. 

  844. int layer_max_qband[]

    845. 

  845. )

    846. 

  846. {

    847. 

  847. 	int i, reg;

    848. 

  848.     int layer, g;

    849. 

  849. 	int sfb, cband;

    850. 

  850. 	int init_max_freq[8];

    851. 

  851. 	int end_cband[8];

    852. 

  852. 	int end_freq[8];

    853. 

  853. 	int end_layer;

    854. 

  854. 

    855. 

  855. 	/************************************************/

    856. 

  856. 	/* calculate slayer_size                        */

    857. 

  857. 	/************************************************/

    858. 

  858. 	if(signal_type == TRANSIENT_TYPE) {

    859. 

  859. 		*slayer_size = 0;

    860. 

  860. 		for(g = 0; g < num_window_groups; g++) {

    861. 

  861. 			init_max_freq[g] = (base_band * window_group_length[g]) / 32 * 4;

    862. 

  862. 			switch(sampling_rate) {

    863. 

  863. 			case 96000:

    864. 

  864. 			case 88200:

    865. 

  865. 			case 64000: /* zhanjie 0708 */

    866. 

  866. 			case 48000:	case 44100:

    867. 

  867. 				if((init_max_freq[g] % 32) >= 16)

    868. 

  868. 					init_max_freq[g] = (init_max_freq[g]/32) * 32 + 20;

    869. 

  869. 				else if((init_max_freq[g] % 32) >= 4)

    870. 

  870. 					init_max_freq[g] = (init_max_freq[g]/32) * 32 + 8;

    871. 

  871. 				break;

    872. 

  872. 			case 32000:	case 24000:	case 22050:

    873. 

  873. 				init_max_freq[g] = (init_max_freq[g]/16) * 16;

    874. 

  874. 				break;

    875. 

  875. 			case 16000:	case 12000:	case 11025:

    876. 

  876. 				init_max_freq[g] = (init_max_freq[g]/32) * 32;

    877. 

  877. 				break;

    878. 

  878. 			case 8000:

    879. 

  879. 				init_max_freq[g] = (init_max_freq[g]/64) * 64;

    880. 

  880. 				break;

    881. 

  881. 			}

    882. 

  882. 

    883. 

  883. 			init_max_freq[g] = min(init_max_freq[g], swb_offset[g][maxSfb]);

    884. 

  884. 

    885. 

  885. 			end_freq[g] = init_max_freq[g];

    886. 

  886. 			end_cband[g] = (end_freq[g] + 31) / 32;

    887. 

  887. 			*slayer_size += (init_max_freq[g]+31)/32;

    888. 

  888. 		}

    889. 

  889. 	} else { /* if(signal_type == TRANSIENT_TYPE) */

    890. 

  890. 		g = 0;

    891. 

  891. 		init_max_freq[g] = base_band;

    892. 

  892. 

    893. 

  893. 		if (init_max_freq[g] > swb_offset[g][maxSfb])

    894. 

  894. 			init_max_freq[g] = swb_offset[g][maxSfb];

    895. 

  895. 		*slayer_size = (init_max_freq[g]+31)/32;

    896. 

  896. 		end_freq[g] = init_max_freq[g];

    897. 

  897. 		end_cband[g] = (end_freq[g] + 31) / 32;

    898. 

  898. 	} /* if(signal_type == TRANSIENT_TYPE) ... else  */

    899. 

  899. 

    900. 

  900. 	/************************************************/

    901. 

  901. 	/* calculate layer_max_cband and layer_max_freq */

    902. 

  902. 	/************************************************/

    903. 

  903. 	for(g = 0, layer=0; g < num_window_groups; g++)

    904. 

  904. 	for (cband = 1; cband <= end_cband[g]; layer++, cband++) {

    905. 

  905. 		layer_reg[layer] = g;

    906. 

  906. 		layer_max_freq[layer] = cband * 32;

    907. 

  907. 		if (layer_max_freq[layer] > init_max_freq[g])

    908. 

  908. 			layer_max_freq[layer] = init_max_freq[g];

    909. 

  909. 		layer_max_cband[layer] = cband;

    910. 

  910. 	}

    911. 

  911. 	if(signal_type == TRANSIENT_TYPE) {	

    912. 

  912. 	   layer = *slayer_size;

    913. 

  913. 	   for(g = 0; g < num_window_groups; g++) {

    914. 

  914. 	  	   for (i=0; i<window_group_length[g]; i++) { 

    915. 

  915. 		  	   layer_reg[layer] = g;

    916. 

  916. 			   layer++;

    917. 

  917. 		   }

    918. 

  918. 	   }

    919. 

  919. 	   for (layer=*slayer_size+8; layer<100; layer++) 

    920. 

  920. 	      layer_reg[layer] = layer_reg[layer-8];

    921. 

  921.     }

    922. 

  922. 	else {

    923. 

  923. 	   for (layer=*slayer_size; layer<100; layer++) 

    924. 

  924. 	      layer_reg[layer] = 0;

    925. 

  925. 	}

    926. 

  926.    

    927. 

  927. 	for (layer = *slayer_size; layer <100; layer++) {

    928. 

  928. 		reg = layer_reg[layer];

    929. 

  929. 		switch(sampling_rate) {

    930. 

  930. 			/* zhanjie 0708 begin */

    931. 

  931. 		case 96000:

    932. 

  932. 		case 88200:

    933. 

  933. 			end_freq[reg] += 16;			

    934. 

  934. 			break;

    935. 

  935. 		case 64000:

    936. 

  936. 			end_freq[reg] += 12;

    937. 

  937. 			break;

    938. 

  938. 			/* zhanjie 0708 ends */

    939. 

  939. 

    940. 

  940. 	    case 48000: case 44100:

    941. 

  941. 		    end_freq[reg] += 8;

    942. 

  942. 		    if ((end_freq[reg]%32) != 0)

    943. 

  943. 				end_freq[reg] += 4;

    944. 

  944. 		    break;

    945. 

  945. 		case 32000:	case 24000:	case 22050:

    946. 

  946. 			end_freq[reg] += 16;

    947. 

  947. 			break;

    948. 

  948. 		case 16000:	case 12000:	case 11025:

    949. 

  949. 			end_freq[reg] += 32;

    950. 

  950. 		    break;

    951. 

  951. 		default:

    952. 

  952. 			end_freq[reg] += 64;

    953. 

  953. 			break;

    954. 

  954. 		}

    955. 

  955. 		if (end_freq[reg] > swb_offset[reg][maxSfb] )

    956. 

  956. 			end_freq[reg] = swb_offset[reg][maxSfb];

    957. 

  957. 		//if( (top_layer+*slayer_size) <= layer)//

    958. 

  958. 		//	end_freq[reg] = swb_offset[reg][maxSfb];

    959. 

  959. 

    960. 

  960. 

    961. 

  961. 		layer_max_freq[layer] = end_freq[reg];

    962. 

  962. 		layer_max_cband[layer] = (layer_max_freq[layer]+31)/32;

    963. 

  963. 	}

    964. 

  964. 

    965. 

  965. 	/*****************************/

    966. 

  966. 	/* calculate layer_bit_flush */

    967. 

  967. 	/*****************************/

    968. 

  968. #if 0

    969. 

  969. 	for (layer = 0; layer < (top_layer+*slayer_size-1); layer++) {

    970. 

  970. 		if (layer_max_cband[layer] != layer_max_cband[layer+1] ||

    971. 

  971. 			layer_reg[layer] != layer_reg[layer+1] )

    972. 

  972. 			layer_bit_flush[layer] = 1;

    973. 

  973. 		else

    974. 

  974. 			layer_bit_flush[layer] = 0;

    975. 

  975. 	}

    976. 

  976. 	for ( ; layer < 100; layer++)

    977. 

  977. 		layer_bit_flush[layer] = 1;

    978. 

  978. #endif	// this part is conflit with standard spec. -- lixun 08/07/21

    979. 

  979. 

    980. 

  980. 	/*****************************/

    981. 

  981. 	/* calculate layer_max_qband */

    982. 

  982. 	/*****************************/

    983. 

  983. 	for (layer = 0; layer < 100; layer++) {

    984. 

  984. 		end_layer = layer;

    985. 

  985. #if 0

    986. 

  986. 		while ( !layer_bit_flush[end_layer] ) 

    987. 

  987. 			end_layer++;

    988. 

  988. #endif		// this part is conflict with spec. -- lixun 08/07/21

    989. 

  989. 		reg = layer_reg[layer];

    990. 

  990. 		for (sfb=0; sfb<maxSfb; sfb++) {

    991. 

  991. 			if (layer_max_freq[end_layer] <= swb_offset[reg][sfb+1]) {

    992. 

  992. 				layer_max_qband[layer] = sfb+1;

    993. 

  993. 				break;

    994. 

  994. 			}

    995. 

  995. 		}

    996. 

  996. 		if (layer_max_qband[layer] > maxSfb)

    997. 

  997. 			layer_max_qband[layer] = maxSfb;

    998. 

  998. 	}

    999. 

  999. }

    1000. 

  1000. 

    1001. 

  1001. void get_end_band (

    1002. 

  1002. int signal_type,

    1003. 

  1003. int num_window_groups,

    1004. 

  1004. int window_group_length[],

    1005. 

  1005. int maxSfb,

    1006. 

  1006. short	*swb_offset[],

    1007. 

  1007. int top_layer,

    1008. 

  1008. int base_maxSfb[],

    1009. 

  1009. int end_sfb[],

    1010. 

  1010. int end_freq[],

    1011. 

  1011. int end_cband[]

    1012. 

  1012. )

    1013. 

  1013. {

    1014. 

  1014. 	int i, reg;

    1015. 

  1015.     int layer, g;

    1016. 

  1016. 	int sfb, cband;

    1017. 

  1017. 	int init_max_freq[8];

    1018. 

  1018.     int slayer_size;

    1019. 

  1019.     int layer_reg[100];

    1020. 

  1020.     int base_band = get_base_band();

    1021. 

  1021. 

    1022. 

  1022. 	/************************************************/

    1023. 

  1023. 	/* calculate slayer_size                        */

    1024. 

  1024. 	/************************************************/

    1025. 

  1025. 	if(signal_type == TRANSIENT_TYPE) {

    1026. 

  1026. 		slayer_size = 0;

    1027. 

  1027. 		for(g = 0; g < num_window_groups; g++) {

    1028. 

  1028. 			init_max_freq[g] = (base_band * window_group_length[g]) / 32 * 4;

    1029. 

  1029. 			switch(sampling_rate) {

    1030. 

  1030. 			case 96000: 

    1031. 

  1031. 			case 88200:

    1032. 

  1032. 			case 64000: // zhanjie 0708

    1033. 

  1033. 			case 48000:	case 44100:

    1034. 

  1034. 				if((init_max_freq[g] % 32) >= 16)

    1035. 

  1035. 					init_max_freq[g] = (init_max_freq[g]/32) * 32 + 20;

    1036. 

  1036. 				else if((init_max_freq[g] % 32) >= 4)

    1037. 

  1037. 					init_max_freq[g] = (init_max_freq[g]/32) * 32 + 8;

    1038. 

  1038. 				break;

    1039. 

  1039. 			case 32000:	case 24000:	case 22050:

    1040. 

  1040. 				init_max_freq[g] = (init_max_freq[g]/16) * 16;

    1041. 

  1041. 				break;

    1042. 

  1042. 			case 16000:	case 12000:	case 11025:

    1043. 

  1043. 				init_max_freq[g] = (init_max_freq[g]/32) * 32;

    1044. 

  1044. 				break;

    1045. 

  1045. 			case 8000:

    1046. 

  1046. 				init_max_freq[g] = (init_max_freq[g]/64) * 64;

    1047. 

  1047. 				break;

    1048. 

  1048. 			}

    1049. 

  1049. 

    1050. 

  1050. 			end_freq[g] = init_max_freq[g];

    1051. 

  1051. 			end_cband[g] = (end_freq[g] + 31) / 32;

    1052. 

  1052. 			slayer_size += (init_max_freq[g]+31)/32;

    1053. 

  1053. 		}

    1054. 

  1054. 	} else { /* if(signal_type == TRANSIENT_TYPE) */

    1055. 

  1055. 		g = 0;

    1056. 

  1056. 		init_max_freq[g] = base_band;

    1057. 

  1057. 

    1058. 

  1058. 		if (init_max_freq[g] > swb_offset[g][maxSfb])

    1059. 

  1059. 			init_max_freq[g] = swb_offset[g][maxSfb];

    1060. 

  1060. 		slayer_size = (init_max_freq[g]+31)/32;

    1061. 

  1061. 		end_freq[g] = init_max_freq[g];

    1062. 

  1062. 		end_cband[g] = (end_freq[g] + 31) / 32;

    1063. 

  1063. 	} /* if(signal_type == TRANSIENT_TYPE) ... else  */

    1064. 

  1064. 

    1065. 

  1065. 	/*****************************/

    1066. 

  1066. 	/* calculate base_maxSfb */

    1067. 

  1067. 	/*****************************/

    1068. 

  1068. 	for (g = 0; g < num_window_groups; g++) {

    1069. 

  1069. 		for (sfb=0; sfb<maxSfb; sfb++) {

    1070. 

  1070. 			if (init_max_freq[g] <= swb_offset[g][sfb+1]) {

    1071. 

  1071. 				base_maxSfb[g] = sfb+1;

    1072. 

  1072. 				break;

    1073. 

  1073. 			}

    1074. 

  1074. 		}

    1075. 

  1075. 		if (base_maxSfb[g] > maxSfb) base_maxSfb[g] = maxSfb;

    1076. 

  1076. 	}

    1077. 

  1077. 

    1078. 

  1078. 	/************************************************/

    1079. 

  1079. 	/* calculate layer_max_cband and layer_max_freq */

    1080. 

  1080. 	/************************************************/

    1081. 

  1081. 	for(g = 0, layer=0; g < num_window_groups; g++)

    1082. 

  1082. 	for (cband = 1; cband <= end_cband[g]; layer++, cband++) {

    1083. 

  1083. 		layer_reg[layer] = g;

    1084. 

  1084. 	}

    1085. 

  1085. 	if(signal_type == TRANSIENT_TYPE) {

    1086. 

  1086. 	   layer = slayer_size;

    1087. 

  1087. 	   for(g = 0; g < num_window_groups; g++) {

    1088. 

  1088. 	  	   for (i=0; i<window_group_length[g]; i++) { 

    1089. 

  1089. 		  	   layer_reg[layer] = g;

    1090. 

  1090. 			   layer++;

    1091. 

  1091. 		   }

    1092. 

  1092. 	   }

    1093. 

  1093. 	   for (layer=slayer_size+8; layer<(top_layer+slayer_size); layer++) 

    1094. 

  1094. 	      layer_reg[layer] = layer_reg[layer-8];

    1095. 

  1095.     }

    1096. 

  1096. 	else {

    1097. 

  1097. 	   for (layer=slayer_size; layer<(top_layer+slayer_size); layer++) 

    1098. 

  1098. 	      layer_reg[layer] = 0;

    1099. 

  1099. 	}

    1100. 

  1100.    

    1101. 

  1101. 	for (layer = slayer_size; layer <(top_layer+slayer_size); layer++) {

    1102. 

  1102. 		reg = layer_reg[layer];

    1103. 

  1103. 		switch(sampling_rate) {

    1104. 

  1104. 			/* zhanjie 0708 begin */

    1105. 

  1105. 		case 96000:

    1106. 

  1106. 		case 88200:

    1107. 

  1107. 			end_freq[reg] += 16;			

    1108. 

  1108. 			break;

    1109. 

  1109. 		case 64000:

    1110. 

  1110. 			end_freq[reg] += 12;

    1111. 

  1111. 			break;

    1112. 

  1112. 			/* zhanjie 0708 end */

    1113. 

  1113. 

    1114. 

  1114. 	    case 48000: case 44100:

    1115. 

  1115. 		    end_freq[reg] += 8;

    1116. 

  1116. 		    if ((end_freq[reg]%32) != 0) end_freq[reg] += 4;

    1117. 

  1117. 		    break;

    1118. 

  1118. 		case 32000:	case 24000:	case 22050:

    1119. 

  1119. 			end_freq[reg] += 16;

    1120. 

  1120. 			break;

    1121. 

  1121. 		case 16000:	case 12000:	case 11025:

    1122. 

  1122. 			end_freq[reg] += 32;

    1123. 

  1123. 		    break;

    1124. 

  1124. 		default:

    1125. 

  1125. 			end_freq[reg] += 64;

    1126. 

  1126. 			break;

    1127. 

  1127. 		}

    1128. 

  1128. 	}

    1129. 

  1129. 

    1130. 

  1130. 

    1131. 

  1131. 	/*****************************/

    1132. 

  1132. 	/* calculate layer_max_qband */

    1133. 

  1133. 	/*****************************/

    1134. 

  1134. 	for (reg=0; reg < num_window_groups; reg++) {

    1135. 

  1135. 	    if (end_freq[reg] > swb_offset[reg][maxSfb] )

    1136. 

  1136. 		    end_freq[reg] = swb_offset[reg][maxSfb];

    1137. 

  1137. 		end_cband[reg] = (end_freq[reg]+31)/32;

    1138. 

  1138. 		for (sfb=0; sfb<maxSfb; sfb++) {

    1139. 

  1139. 			if (end_freq[reg] <= swb_offset[reg][sfb+1]) {

    1140. 

  1140. 				end_sfb[reg] = sfb+1;

    1141. 

  1141. 				break;

    1142. 

  1142. 			}

    1143. 

  1143. 		}

    1144. 

  1144. 		if (end_sfb[reg] > maxSfb)

    1145. 

  1145. 			end_sfb[reg] = maxSfb;

    1146. 

  1146. 	}

    1147. 

  1147. }

    1148. 

  1148. 

    1149. 

  1149. void get_base_maxSfb (

    1150. 

  1150. int signal_type,

    1151. 

  1151. int num_window_groups,

    1152. 

  1152. int window_group_length[],

    1153. 

  1153. int maxSfb,

    1154. 

  1154. short	*swb_offset[],

    1155. 

  1155. int base_maxSfb[]

    1156. 

  1156. )

    1157. 

  1157. {

    1158. 

  1158.     int g;

    1159. 

  1159. 	int sfb;

    1160. 

  1160.     int base_band = get_base_band();

    1161. 

  1161. 	int init_max_freq[8];

    1162. 

  1162. 

    1163. 

  1163. 	/************************************************/

    1164. 

  1164. 	/* calculate slayer_size                        */

    1165. 

  1165. 	/************************************************/

    1166. 

  1166. 	if(signal_type == TRANSIENT_TYPE) {

    1167. 

  1167. 		for(g = 0; g < num_window_groups; g++) {

    1168. 

  1168. 			init_max_freq[g] = (base_band * window_group_length[g]) / 32 * 4;

    1169. 

  1169. 			switch(sampling_rate) {

    1170. 

  1170. 			case 96000:

    1171. 

  1171. 			case 88200:

    1172. 

  1172. 			case 64000: /* zhanjie 0708 */

    1173. 

  1173. 			case 48000:	case 44100:

    1174. 

  1174. 				if((init_max_freq[g] % 32) >= 16)

    1175. 

  1175. 					init_max_freq[g] = (init_max_freq[g]/32) * 32 + 20;

    1176. 

  1176. 				else if((init_max_freq[g] % 32) >= 4)

    1177. 

  1177. 					init_max_freq[g] = (init_max_freq[g]/32) * 32 + 8;

    1178. 

  1178. 				break;

    1179. 

  1179. 			case 32000:	case 24000:	case 22050:

    1180. 

  1180. 				init_max_freq[g] = (init_max_freq[g]/16) * 16;

    1181. 

  1181. 				break;

    1182. 

  1182. 			case 16000:	case 12000:	case 11025:

    1183. 

  1183. 				init_max_freq[g] = (init_max_freq[g]/32) * 32;

    1184. 

  1184. 				break;

    1185. 

  1185. 			case 8000:

    1186. 

  1186. 				init_max_freq[g] = (init_max_freq[g]/64) * 64;

    1187. 

  1187. 				break;

    1188. 

  1188. 			}

    1189. 

  1189. 		}

    1190. 

  1190. 	} else { /* if(signal_type == TRANSIENT_TYPE) */

    1191. 

  1191. 		init_max_freq[0] = base_band;

    1192. 

  1192. 	} /* if(signal_type == TRANSIENT_TYPE) ... else  */

    1193. 

  1193. 

    1194. 

  1194. 	/*****************************/

    1195. 

  1195. 	/* calculate layer_max_qband */

    1196. 

  1196. 	/*****************************/

    1197. 

  1197. 	for (g = 0; g < num_window_groups; g++) {

    1198. 

  1198. 		for (sfb=0; sfb<maxSfb; sfb++) {

    1199. 

  1199. 			if (init_max_freq[g] <= swb_offset[g][sfb+1]) {

    1200. 

  1200. 				base_maxSfb[g] = sfb+1;

    1201. 

  1201. 				break;

    1202. 

  1202. 			}

    1203. 

  1203. 		}

    1204. 

  1204. 		if (base_maxSfb[g] > maxSfb) base_maxSfb[g] = maxSfb;

    1205. 

  1205. 	}

    1206. 

  1206. }

    1207. 

  1207. 

    1208. 

  1208. int cband_si_esc[9] = 

    1209. 

  1209. {-4, -3, -2, -1, 30, 31, 32, 33, 34};

    1210. 

  1210. 

    1211. 

  1211. /*--------------------------------------------------------------*/

    1212. 

  1212. /***************   coding band side infomation   ****************/

    1213. 

  1213. /*--------------------------------------------------------------*/

    1214. 

  1214. int encode_cband_si(

    1215. 

  1215. 	int	*model_index[][8],

    1216. 

  1216. 	int	start_cband[8],

    1217. 

  1217. 	int end_cband[8],

    1218. 

  1218. 	int g,

    1219. 

  1219. 	int nch)

    1220. 

  1220. {

    1221. 

  1221. 	int ch;

    1222. 

  1222. 	int cband;

    1223. 

  1223. 	int	m;

    1224. 

  1224. 	int si_cw_len;

    1225. 

  1225. 	int cband_model;

    1226. 

  1226. 	

    1227. 

  1227. 	static int prevSi[2][8];

    1228. 

  1228. 	int t1,t2;

    1229. 

  1229. 	extern int siCodeMode;

    1230. 

  1230. 

    1231. 

  1231. 	si_cw_len = 0;

    1232. 

  1232. 	for(ch = 0; ch < nch; ch++) {

    1233. 

  1233. 		if(siCodeMode==0)

    1234. 

  1234. 		{

    1235. 

  1235. 			if(start_cband[g]==0 && g==0)

    1236. 

  1236. 				prevSi[ch][0] = 0;

    1237. 

  1237. 			for (cband = start_cband[g]; cband < end_cband[g]; cband++) {

    1238. 

  1238. 				t2 = model_index[ch][g][cband];

    1239. 

  1239. 				if(t2<-10000)

    1240. 

  1240. 					t2 = t2;

    1241. 

  1241. 		

    1242. 

  1242. 				if((prevSi[ch][0]-t2+15)<0 || (prevSi[ch][0]-t2+15)>29 )

    1243. 

  1243. 				{

    1244. 

  1244. 					/* coding escape code */

    1245. 

  1245. 					si_cw_len += EncodeBSHCSi(30);

    1246. 

  1246. 					/* encoding 9 escape value with fixed 4 bit */

    1247. 

  1247. 					t1 = prevSi[ch][0]-t2+15;

    1248. 

  1248. 					m = 0;

    1249. 

  1249. 					while(cband_si_esc[m] != t1)

    1250. 

  1250. 						m ++;

    1251. 

  1251. 					if(m & 0x8)

    1252. 

  1252. 						si_cw_len += EncodeBSHCBin(1);

    1253. 

  1253. 					else

    1254. 

  1254. 						si_cw_len += EncodeBSHCBin(0);

    1255. 

  1255. 					if(m & 0x4)

    1256. 

  1256. 						si_cw_len += EncodeBSHCBin(1);

    1257. 

  1257. 					else

    1258. 

  1258. 						si_cw_len += EncodeBSHCBin(0);

    1259. 

  1259. 					if(m & 0x2)

    1260. 

  1260. 						si_cw_len += EncodeBSHCBin(1);

    1261. 

  1261. 					else

    1262. 

  1262. 						si_cw_len += EncodeBSHCBin(0);

    1263. 

  1263. 					if(m & 0x1)

    1264. 

  1264. 						si_cw_len += EncodeBSHCBin(1);

    1265. 

  1265. 					else

    1266. 

  1266. 						si_cw_len += EncodeBSHCBin(0);

    1267. 

  1267. 				}

    1268. 

  1268. 				else

    1269. 

  1269. 					si_cw_len += EncodeBSHCSi(prevSi[ch][0]-t2+15);

    1270. 

  1270. 				

    1271. 

  1271. 				prevSi[ch][0] = t2;

    1272. 

  1272. 			}

    1273. 

  1273. 		}

    1274. 

  1274. 		else

    1275. 

  1275. 		{

    1276. 

  1276. 			if(start_cband[g]==0)

    1277. 

  1277. 				prevSi[ch][g] = 0;

    1278. 

  1278. 			for (cband = start_cband[g]; cband < end_cband[g]; cband++) {

    1279. 

  1279. 				t2 = model_index[ch][g][cband];

    1280. 

  1280. 				

    1281. 

  1281. 				if((prevSi[ch][g]-t2+15)<0 || (prevSi[ch][g]-t2+15)>29 )

    1282. 

  1282. 				{

    1283. 

  1283. 					/* coding escape code */

    1284. 

  1284. 					si_cw_len += EncodeBSHCSi(30);

    1285. 

  1285. 					/* encoding 9 escape value with fixed 4 bit */

    1286. 

  1286. 					t1 = prevSi[ch][g]-t2+15;

    1287. 

  1287. 					m = 0;

    1288. 

  1288. 					while(cband_si_esc[m] != t1)

    1289. 

  1289. 						m ++;

    1290. 

  1290. 					if(m & 0x8)

    1291. 

  1291. 						si_cw_len += EncodeBSHCBin(1);

    1292. 

  1292. 					else

    1293. 

  1293. 						si_cw_len += EncodeBSHCBin(0);

    1294. 

  1294. 					if(m & 0x4)

    1295. 

  1295. 						si_cw_len += EncodeBSHCBin(1);

    1296. 

  1296. 					else

    1297. 

  1297. 						si_cw_len += EncodeBSHCBin(0);

    1298. 

  1298. 					if(m & 0x2)

    1299. 

  1299. 						si_cw_len += EncodeBSHCBin(1);

    1300. 

  1300. 					else

    1301. 

  1301. 						si_cw_len += EncodeBSHCBin(0);

    1302. 

  1302. 					if(m & 0x1)

    1303. 

  1303. 						si_cw_len += EncodeBSHCBin(1);

    1304. 

  1304. 					else

    1305. 

  1305. 						si_cw_len += EncodeBSHCBin(0);

    1306. 

  1306. 				}

    1307. 

  1307. 				else

    1308. 

  1308. 					si_cw_len += EncodeBSHCSi(prevSi[ch][g]-t2+15);

    1309. 

  1309. 				

    1310. 

  1310. 				prevSi[ch][g] = t2;

    1311. 

  1311. 			}		

    1312. 

  1312. 		}

    1313. 

  1313. 	}

    1314. 

  1314. 

    1315. 

  1315. 	return si_cw_len;

    1316. 

  1316. }

    1317. 

  1317. 	

    1318. 

  1318. /*--------------------------------------------------------------*/

    1319. 

  1319. /***********  scalefactor band side infomation   ****************/

    1320. 

  1320. /*--------------------------------------------------------------*/

    1321. 

  1321. int encode_scfband_si(

    1322. 

  1322. 	int	*scf[][8],

    1323. 

  1323. 	int	start_sfb[8],

    1324. 

  1324. 	int end_sfb[8],

    1325. 

  1325. 	int stereo_mode,

    1326. 

  1326. 	int *stereo_info[8],

    1327. 

  1327. 	int stereo_si_coded[8][MAX_SCFAC_BANDS],

    1328. 

  1328. 	int g,

    1329. 

  1329. 	int nch)

    1330. 

  1330. {

    1331. 

  1331. 	int ch;

    1332. 

  1332. 	int	m;

    1333. 

  1333. 	int sfb;

    1334. 

  1334. 	int si_cw_len;

    1335. 

  1335. 	

    1336. 

  1336. 	static int prevSf[2][8];

    1337. 

  1337. 	if(start_sfb[g]==0)

    1338. 

  1338. 		for(ch = 0; ch < nch; ch++)

    1339. 

  1339. 			prevSf[ch][g] = scf[ch][0][0];

    1340. 

  1340. 	si_cw_len = 0;

    1341. 

  1341. 	for(ch = 0; ch < nch; ch++) {

    1342. 

  1342. 		

    1343. 

  1343. 		for(sfb = start_sfb[g]; sfb < end_sfb[g]; sfb++) {

    1344. 

  1344. 			if(stereo_si_coded[g][sfb] == 0) {

    1345. 

  1345. 				if(stereo_mode != 2) {

    1346. 

  1346. 					m = stereo_info[g][sfb];

    1347. 

  1347. 					if(stereo_mode == 1) {						

    1348. 

  1348. 						EncodeBSHCBin(m);

    1349. 

  1349. 						si_cw_len++;

    1350. 

  1350. 					}

    1351. 

  1351. 				}

    1352. 

  1352. 				stereo_si_coded[g][sfb] = 1;

    1353. 

  1353. 			}

    1354. 

  1354. 			

    1355. 

  1355. 			if(!(g==0 && sfb==0))

    1356. 

  1356. 			{

    1357. 

  1357.                                           

    1358. 

  1358. 				if((scf[ch][g][sfb]-prevSf[ch][g]+64)<0 || (scf[ch][g][sfb]-prevSf[ch][g]+64)>127)

    1359. 

  1359. 				{

    1360. 

  1360. 					si_cw_len += EncodeBSHCSf(64);

    1361. 

  1361. 				}

    1362. 

  1362. 				else

    1363. 

  1363.                                            

    1364. 

  1364. 				{

    1365. 

  1365. 					if ((stereo_mode == 2 || (stereo_mode == 1 && stereo_info[g][sfb] == 1)) && ch!=0)

    1366. 

  1366. 						si_cw_len += EncodeBSHCSf(scf[ch][g][sfb] - scf[0][g][sfb] + 64);

    1367. 

  1367. 					else

    1368. 

  1368. 						si_cw_len += EncodeBSHCSf(scf[ch][g][sfb]-prevSf[ch][g]+64);

    1369. 

  1369. 					prevSf[ch][g] = scf[ch][g][sfb];

    1370. 

  1370. 				}			

    1371. 

  1371. 			}			

    1372. 

  1372. 		}

    1373. 

  1373. 	}

    1374. 

  1374. 

    1375. 

  1375. 	return si_cw_len;

    1376. 

  1376. }

    1377. 

  1377. 

    1378. 

  1378. int cb_tab[16] = {0,1,2,3,4,3,3,5,6,3,3,5,3,5,5,5};

    1379. 

  1379. 

    1380. 

  1380. /*--------------------------------------------------------------*/

    1381. 

  1381. /***********     Bit-Sliced Spectral Data        ****************/

    1382. 

  1382. /*--------------------------------------------------------------*/

    1383. 

  1383. int encode_spectra(

    1384. 

  1384. 	int	*sample[][8],

    1385. 

  1385. 	int s_reg,

    1386. 

  1386. 	int e_reg,

    1387. 

  1387. 	int s_freq[8],

    1388. 

  1388. 	int e_freq[8],

    1389. 

  1389. 	int	min_bpl,

    1390. 

  1390. 	int available_len,

    1391. 

  1391. 	int *model_index[][8],

    1392. 

  1392. 	int *cur_bpl[][8],

    1393. 

  1393. 	int *coded_samp_bit[][8],

    1394. 

  1394. 	int *sign_coded[][8],

    1395. 

  1395. 	int *enc_vec[][8],

    1396. 

  1396. 	int	nch)

    1397. 

  1397. {

    1398. 

  1398. 	int ch, i, m, k, i4;

    1399. 

  1399. 	int maxbpl;

    1400. 

  1400. 	int	bpl;

    1401. 

  1401. 	int	cband;

    1402. 

  1402. 	int	maxhalf;

    1403. 

  1403. 	int	cw_len;

    1404. 

  1404. 	int	reg;

    1405. 

  1405. 

    1406. 

  1406. 	

    1407. 

  1407. 	int modelIdx;

    1408. 

  1408. 	int j;

    1409. 

  1409. 

    1410. 

  1410. 	int bit_mask[16] = { 0x0000, 0x0001, 0x0002, 0x0004, 0x0008, 0x0010, 0x0020, 0x0040, 

    1411. 

  1411. 						 0x0080, 0x0100, 0x0200, 0x0400, 0x0800, 0x1000, 0x2000, 0x4000 };  

    1412. 

  1412. 

    1413. 

  1413. 	maxbpl = 0;

    1414. 

  1414. 	for(ch = 0; ch < nch; ch++)

    1415. 

  1415. 	{

    1416. 

  1416. 		for(reg=s_reg; reg<e_reg; reg++)

    1417. 

  1417. 		for(i = s_freq[reg]; i < e_freq[reg]; i++)

    1418. 

  1418. 			if (maxbpl < cur_bpl[ch][reg][i]) 

    1419. 

  1419. 				maxbpl = cur_bpl[ch][reg][i];

    1420. 

  1420. 	}

    1421. 

  1421. 

    1422. 

  1422. 	cw_len = available_len;

    1423. 

  1423. 	for (bpl = maxbpl; bpl >= min_bpl; bpl--) {

    1424. 

  1424. 

    1425. 

  1425. 		if (cw_len <= 0 ) 

    1426. 

  1426. 		{

    1427. 

  1427. //			printf("\nbpl:%d, maxbpl:%d	", bpl, maxbpl);

    1428. 

  1428. 			return (available_len-cw_len);

    1429. 

  1429. 		}

    1430. 

  1430. 		maxhalf = bit_mask[bpl]; 

    1431. 

  1431. 		

    1432. 

  1432. 		for(reg = s_reg; reg < e_reg; reg++)

    1433. 

  1433. 	    for (i = s_freq[reg]; i < e_freq[reg]; i+=4) {

    1434. 

  1434. 		for(ch = 0; ch < nch; ch++) {

    1435. 

  1435. 

    1436. 

  1436. 			if(model_index[ch][reg][i/32]==0)

    1437. 

  1437. 				continue;

    1438. 

  1438. 	 		if ( cur_bpl[ch][reg][i] < bpl )

    1439. 

  1439. 				continue; 

    1440. 

  1440. 

    1441. 

  1441. 			if (coded_samp_bit[ch][reg][i]==0 || sign_coded[ch][reg][i]==1) {

    1442. 

  1442. 

    1443. 

  1443. 				cband = i/32;

    1444. 

  1444. 				i4 = i / 4;

    1445. 

  1445. 				if (i%4==0) {

    1446. 

  1446. 					enc_vec[ch][reg][i4] = 0;

    1447. 

  1447. 				}

    1448. 

  1448. 			

    1449. 

  1449. 				for(j=0;j<4;j++)

    1450. 

  1450. 				{

    1451. 

  1451. 					if( abs(sample[ch][reg][i+j]) & maxhalf) 	m = 1;

    1452. 

  1452. 					else 								    m = 0;

    1453. 

  1453. 

    1454. 

  1454. 					enc_vec[ch][reg][i4] = (enc_vec[ch][reg][i4]<<1) | m;

    1455. 

  1455. 					coded_samp_bit[ch][reg][i+j] = 

    1456. 

  1456. 						(coded_samp_bit[ch][reg][i+j]<<1) | m;

    1457. 

  1457. 					cur_bpl[ch][reg][i+j]-=1;

    1458. 

  1458. 				}

    1459. 

  1459. 				/* model selection */

    1460. 

  1460. 				{

    1461. 

  1461. 					int t[4];

    1462. 

  1462. 					t[0] = (coded_samp_bit[ch][reg][i]&0x2)==0x2 ? 1:0;

    1463. 

  1463. 					t[1] = (coded_samp_bit[ch][reg][i+1]&0x2)==0x2 ? 1:0;

    1464. 

  1464. 					t[2] = (coded_samp_bit[ch][reg][i+2]&0x2)==0x2 ? 1:0;

    1465. 

  1465. 					t[3] = (coded_samp_bit[ch][reg][i+3]&0x2)==0x2 ? 1:0;

    1466. 

  1466. 					t[0] = (t[0]<<3) + (t[1]<<2) + (t[2]<<1) + (t[3]);

    1467. 

  1467. 

    1468. 

  1468. 					modelIdx = BSHCModelSelect(model_index[ch][reg][cband],bpl);

    1469. 

  1469. 					/* upper vector mapping, codebook reduction */

    1470. 

  1470. 					t[0] = cb_tab[t[0]];

    1471. 

  1471. 					k = EncodeBSHCQuant(modelIdx,t[0],enc_vec[ch][reg][i4]);

    1472. 

  1472. 					

    1473. 

  1473. 					cw_len -= k;

    1474. 

  1474. 					if (coded_samp_bit[ch][reg][i] && sign_coded[ch][reg][i]==0)

    1475. 

  1475. 					{

    1476. 

  1476. 						if(sample[ch][reg][i] < 0)

    1477. 

  1477. 							EncodeBSHCBin(1);

    1478. 

  1478. 						else

    1479. 

  1479. 							EncodeBSHCBin(0);

    1480. 

  1480. 						cw_len--;

    1481. 

  1481. 						sign_coded[ch][reg][i] = 0x01;

    1482. 

  1482. 					}

    1483. 

  1483. 					if (coded_samp_bit[ch][reg][i+1] && sign_coded[ch][reg][i+1]==0)

    1484. 

  1484. 					{

    1485. 

  1485. 						if(sample[ch][reg][i+1] < 0)

    1486. 

  1486. 							EncodeBSHCBin(1);

    1487. 

  1487. 						else

    1488. 

  1488. 							EncodeBSHCBin(0);

    1489. 

  1489. 						cw_len--;

    1490. 

  1490. 						sign_coded[ch][reg][i+1] = 0x01;

    1491. 

  1491. 					}

    1492. 

  1492. 					if (coded_samp_bit[ch][reg][i+2] && sign_coded[ch][reg][i+2]==0)

    1493. 

  1493. 					{

    1494. 

  1494. 						if(sample[ch][reg][i+2] < 0)

    1495. 

  1495. 							EncodeBSHCBin(1);

    1496. 

  1496. 						else

    1497. 

  1497. 							EncodeBSHCBin(0);

    1498. 

  1498. 						cw_len--;

    1499. 

  1499. 						sign_coded[ch][reg][i+2] = 0x01;

    1500. 

  1500. 					}

    1501. 

  1501. 					if (coded_samp_bit[ch][reg][i+3] && sign_coded[ch][reg][i+3]==0)

    1502. 

  1502. 					{

    1503. 

  1503. 						if(sample[ch][reg][i+3] < 0)

    1504. 

  1504. 							EncodeBSHCBin(1);

    1505. 

  1505. 						else

    1506. 

  1506. 							EncodeBSHCBin(0);

    1507. 

  1507. 						cw_len--;

    1508. 

  1508. 						sign_coded[ch][reg][i+3] = 0x01;

    1509. 

  1509. 					}

    1510. 

  1510. 				}

    1511. 

  1511. 			}

    1512. 

  1512. 			if (cw_len <= 0 ) 

    1513. 

  1513. 				return (available_len-cw_len);

    1514. 

  1514. 		}

    1515. 

  1515. 	} 

    1516. 

  1516. }

    1517. 

  1517. 	return (available_len-cw_len);

    1518. 

  1518. }

    1519. 

  1519.