开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > DVD > 查看源码
downmix.c
资源名称：dvd-munitions.tar.gz　[点击查看]
上传用户：aoeyumen
上传日期：2007-01-06
资源大小：3329k
文件大小：8k
源码类别：

DVD

开发平台：

Unix_Linux

downmix.c：源码内容
  1. /*
  2.  *
  3.  *  downmix.c
  4.  *    
  5.  * Copyright (C) Aaron Holtzman - Sept 1999
  6.  *
  7.  * Originally based on code by Yuqing Deng.
  8.  *
  9.  *  This file is part of ac3dec, a free Dolby AC-3 stream decoder.
  10.  *
  11.  *  ac3dec is free software; you can redistribute it and/or modify
  12.  *  it under the terms of the GNU General Public License as published by
  13.  *  the Free Software Foundation; either version 2, or (at your option)
  14.  *  any later version.
  15.  *   
  16.  *  ac3dec is distributed in the hope that it will be useful,
  17.  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
  18.  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  19.  *  GNU General Public License for more details.
  20.  *   
  21.  *  You should have received a copy of the GNU General Public License
  22.  *  along with GNU Make; see the file COPYING.  If not, write to
  23.  *  the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. 
  24.  *
  25.  *
  26.  */
  28. #include <stdlib.h>
  29. #include <stdio.h>
  30. #include <math.h>
  31. #include "ac3.h"
  32. #include "decode.h"
  33. #include "downmix.h"
  34. #include "debug.h"
  37. //Pre-scaled downmix coefficients
  38. static float cmixlev_lut[4] = { 0.2928, 0.2468, 0.2071, 0.2468 };
  39. static float smixlev_lut[4] = { 0.2928, 0.2071, 0.0   , 0.2071 };
  42. /* Downmix into _two_ channels...other downmix modes aren't implemented
  43.  * to reduce complexity. Realistically, there aren't many machines around
  44.  * with > 2 channel output anyways */
  46. void downmix(bsi_t* bsi, stream_samples_t* stream_sample)
  47. {
  48. int j;
  49. float right_tmp;
  50. float left_tmp;
  51. float clev,slev;
  52. float *centre = 0, *left = 0, *right = 0, *left_sur = 0, *right_sur = 0;
  54. if(bsi->acmod > 7)
  55. dprintf("(downmix) invalid acmod numbern"); 
  57. //There are two main cases, with or without Dolby Surround
  58. if(global_prefs.use_dolby_surround)
  59. {
  60. switch(bsi->acmod)
  61. {
  62. // 3/2
  63. case 7:
  64. left      = stream_sample->channel[0];
  65. centre    = stream_sample->channel[1];
  66. right     = stream_sample->channel[2];
  67. left_sur  = stream_sample->channel[3];
  68. right_sur = stream_sample->channel[4];
  70. for (j = 0; j < 256; j++) 
  71. {
  72. right_tmp = 0.2265 * *left_sur++ + 0.2265 * *right_sur++;
  73. left_tmp  = -1 * right_tmp;
  74. right_tmp += 0.3204 * *right++ + 0.2265 * *centre;
  75. left_tmp  += 0.3204 * *left++  + 0.2265 * *centre++;
  77. stream_sample->channel[1][j] = right_tmp;
  78. stream_sample->channel[0][j] = left_tmp;
  79. }
  81. break;
  83. // 2/2
  84. case 6:
  85. left      = stream_sample->channel[0];
  86. right     = stream_sample->channel[1];
  87. left_sur  = stream_sample->channel[2];
  88. right_sur = stream_sample->channel[3];
  90. for (j = 0; j < 256; j++) 
  91. {
  92. right_tmp = 0.2265 * *left_sur++ + 0.2265 * *right_sur++;
  93. left_tmp  = -1 * right_tmp;
  94. right_tmp += 0.3204 * *right++;
  95. left_tmp  += 0.3204 * *left++ ;
  97. stream_sample->channel[1][j] = right_tmp;
  98. stream_sample->channel[0][j] = left_tmp;
  99. }
  100. break;
  102. // 3/1
  103. case 5:
  104. left      = stream_sample->channel[0];
  105. centre    = stream_sample->channel[1];
  106. right     = stream_sample->channel[2];
  107. //Mono surround
  108. right_sur = stream_sample->channel[3];
  110. for (j = 0; j < 256; j++) 
  111. {
  112. right_tmp =  0.2265 * *right_sur++;
  113. left_tmp  = -1 * right_tmp;
  114. right_tmp += 0.3204 * *right++ + 0.2265 * *centre;
  115. left_tmp  += 0.3204 * *left++  + 0.2265 * *centre++;
  117. stream_sample->channel[1][j] = right_tmp;
  118. stream_sample->channel[0][j] = left_tmp;
  119. }
  120. break;
  122. // 2/1
  123. case 4:
  124. left      = stream_sample->channel[0];
  125. right     = stream_sample->channel[1];
  126. //Mono surround
  127. right_sur = stream_sample->channel[2];
  129. for (j = 0; j < 256; j++) 
  130. {
  131. right_tmp =  0.2265 * *right_sur++;
  132. left_tmp  = -1 * right_tmp;
  133. right_tmp += 0.3204 * *right++; 
  134. left_tmp  += 0.3204 * *left++;
  136. stream_sample->channel[1][j] = right_tmp;
  137. stream_sample->channel[0][j] = left_tmp;
  138. }
  139. break;
  141. // 3/0
  142. case 3:
  143. left      = stream_sample->channel[0];
  144. centre    = stream_sample->channel[1];
  145. right     = stream_sample->channel[2];
  147. for (j = 0; j < 256; j++) 
  148. {
  149. right_tmp = 0.3204 * *right++ + 0.2265 * *centre;
  150. left_tmp  = 0.3204 * *left++  + 0.2265 * *centre++;
  152. stream_sample->channel[1][j] = right_tmp;
  153. stream_sample->channel[0][j] = left_tmp;
  154. }
  155. break;
  157. // 2/0
  158. case 2:
  159. //Do nothing!
  160. break;
  162. // 1/0
  163. case 1:
  164. //Mono program!
  165. right = stream_sample->channel[0];
  167. for (j = 0; j < 256; j++) 
  168. {
  169. right_tmp = 0.7071 * *right++;
  171. stream_sample->channel[1][j] = right_tmp;
  172. stream_sample->channel[0][j] = right_tmp;
  173. }
  175. break;
  177. // 1+1
  178. case 0:
  179. //Dual mono, output selected by user
  180. right = stream_sample->channel[global_prefs.dual_mono_channel_select];
  182. for (j = 0; j < 256; j++) 
  183. {
  184. right_tmp = 0.7071 * *right++;
  186. stream_sample->channel[1][j] = right_tmp;
  187. stream_sample->channel[0][j] = right_tmp;
  188. }
  189. break;
  190. }
  191. }
  192. else
  193. {
  194. //Non-Dolby surround downmixes
  195. switch(bsi->acmod)
  196. {
  197. // 3/2
  198. case 7:
  199. left      = stream_sample->channel[0];
  200. centre    = stream_sample->channel[1];
  201. right     = stream_sample->channel[2];
  202. left_sur  = stream_sample->channel[3];
  203. right_sur = stream_sample->channel[4];
  205. clev = cmixlev_lut[bsi->cmixlev];
  206. slev = smixlev_lut[bsi->surmixlev];
  208. for (j = 0; j < 256; j++) 
  209. {
  210. right_tmp= 0.4142 * *right++ + clev * *centre   + slev * *right_sur++;
  211. left_tmp = 0.4142 * *left++  + clev * *centre++ + slev * *left_sur++;
  213. stream_sample->channel[1][j] = right_tmp;
  214. stream_sample->channel[0][j] = left_tmp;
  215. }
  217. break;
  219. // 2/2
  220. case 6:
  221. left      = stream_sample->channel[0];
  222. right     = stream_sample->channel[1];
  223. left_sur  = stream_sample->channel[2];
  224. right_sur = stream_sample->channel[3];
  226. slev = smixlev_lut[bsi->surmixlev];
  228. for (j = 0; j < 256; j++) 
  229. {
  230. right_tmp= 0.4142 * *right++ + slev * *right_sur++;
  231. left_tmp = 0.4142 * *left++  + slev * *left_sur++;
  233. stream_sample->channel[1][j] = right_tmp;
  234. stream_sample->channel[0][j] = left_tmp;
  235. }
  237. // 3/1
  238. case 5:
  239. left      = stream_sample->channel[0];
  240. centre    = stream_sample->channel[1];
  241. right     = stream_sample->channel[2];
  242. //Mono surround
  243. right_sur = stream_sample->channel[3];
  245. clev = cmixlev_lut[bsi->cmixlev];
  246. slev = smixlev_lut[bsi->surmixlev];
  248. for (j = 0; j < 256; j++) 
  249. {
  250. right_tmp= 0.4142 * *right++ + clev * *centre   + slev * *right_sur;
  251. left_tmp = 0.4142 * *left++  + clev * *centre++ + slev * *right_sur++;
  253. stream_sample->channel[1][j] = right_tmp;
  254. stream_sample->channel[0][j] = left_tmp;
  255. }
  257. break;
  259. // 2/1
  260. case 4:
  261. left      = stream_sample->channel[0];
  262. right     = stream_sample->channel[1];
  263. //Mono surround
  264. right_sur = stream_sample->channel[2];
  266. slev = smixlev_lut[bsi->surmixlev];
  268. for (j = 0; j < 256; j++) 
  269. {
  270. right_tmp= 0.4142 * *right++ + slev * *right_sur;
  271. left_tmp = 0.4142 * *left++  + slev * *right_sur++;
  273. stream_sample->channel[1][j] = right_tmp;
  274. stream_sample->channel[0][j] = left_tmp;
  275. }
  277. // 3/0
  278. case 3:
  279. left      = stream_sample->channel[0];
  280. centre    = stream_sample->channel[1];
  281. right     = stream_sample->channel[2];
  283. clev = cmixlev_lut[bsi->cmixlev];
  285. for (j = 0; j < 256; j++) 
  286. {
  287. right_tmp= 0.4142 * *right++ + clev * *centre;   
  288. left_tmp = 0.4142 * *left++  + clev * *centre++; 
  290. stream_sample->channel[1][j] = right_tmp;
  291. stream_sample->channel[0][j] = left_tmp;
  292. }
  294. break;
  296. case 2:
  297. //Do nothing!
  298. break;
  300. // 1/0
  301. case 1:
  302. //Mono program!
  303. right = stream_sample->channel[0];
  305. for (j = 0; j < 256; j++) 
  306. {
  307. right_tmp = 0.7071 * *right++;
  309. stream_sample->channel[1][j] = right_tmp;
  310. stream_sample->channel[0][j] = right_tmp;
  311. }
  313. break;
  315. // 1+1
  316. case 0:
  317. //Dual mono, output selected by user
  318. right = stream_sample->channel[global_prefs.dual_mono_channel_select];
  320. for (j = 0; j < 256; j++) 
  321. {
  322. right_tmp = 0.7071 * *right++;
  324. stream_sample->channel[1][j] = right_tmp;
  325. stream_sample->channel[0][j] = right_tmp;
  326. }
  327. break;
  329. }
  330. }
  331. }