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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > mpeg/mp3 > 查看源码
MUSICOUT.C
资源名称：mpeg.zip　[点击查看]
上传用户：njqiyou
上传日期：2007-01-08
资源大小：574k
文件大小：50k
源码类别：

mpeg/mp3

开发平台：

C/C++

MUSICOUT.C：源码内容
  1. /**********************************************************************
  2.  * ISO MPEG Audio Subgroup Software Simulation Group (1996)
  3.  * ISO 13818-3 MPEG-2 Audio Multichannel Decoder
  4.  *
  5.  * $Id: musicout.c 1.11 1996/04/18 05:37:23 rowlands Exp $
  6.  *
  7.  * $Log: musicout.c $
  8.  * Revision 1.11  1996/04/18 05:37:23  rowlands
  9.  * Release following Florence meeting
  10.  *
  11.  * Revision 1.10  1996/02/12 07:13:06  rowlands
  12.  * Release following Munich meeting
  13.  *
  14.  * Revision 1.9.1.1  1996/01/20  17:28:45  rowlands
  15.  * Received from Ralf Schwalbe (Telekom FTZ) - includes prediction
  16.  *
  17.  * Revision 1.7.1.3  1995/08/14  08:13:37  tenkate
  18.  * ML-LSF added Warner ten Kate 7/8/95 (Philips)
  19.  * debugging
  20.  * change "alloc" and "sblimit" into "alloc_ml" and
  21.  * "sblimit_ml" where appropriate.
  22.  * adapt sample-loops to 6 or 12 depending on half or full ML rate.
  23.  *
  24.  * Revision 1.7.1.1  1995/07/14  06:12:46  rowlands
  25.  * Updated dynamic crosstalk from FTZ: revision FTZ_03
  26.  *
  27.  * Revision 1.3.3.1  1995/06/16  08:00:46  rowlands
  28.  * Input from Sang Wook Kim (Samsung AIT)
  29.  *
  30.  * Revision 1.3.1.1  1995/06/14  04:36:30  rowlands
  31.  * Added support for dematrixing procedure 2.
  32.  *
  33.  **********************************************************************/
  35. /**********************************************************************
  36.  *   date   programmers                comment                        *
  37.  * 2/25/91  Douglas Wong        start of version 1.0 records          *
  38.  * 3/06/91  Douglas Wong        rename setup.h to dedef.h             *
  39.  *                              removed extraneous variables          *
  40.  *                              removed window_samples (now part of   *
  41.  *                              filter_samples)                       *
  42.  * 3/07/91  Davis Pan           changed output file to "codmusic"     *
  43.  * 5/10/91  Vish (PRISM)        Ported to Macintosh and Unix.         *
  44.  *                              Incorporated new "out_fifo()" which   *
  45.  *                              writes out last incomplete buffer.    *
  46.  *                              Incorporated all AIFF routines which  *
  47.  *                              are also compatible with SUN.         *
  48.  *                              Incorporated user interface for       *
  49.  *                              specifying sound file names.          *
  50.  *                              Also incorporated user interface for  *
  51.  *                              writing AIFF compatible sound files.  *
  52.  * 27jun91  dpwe (Aware)        Added musicout and &sample_frames as  *
  53.  *                              args to out_fifo (were glob refs).    *
  54.  *                              Used new 'frame_params' struct.       *
  55.  *                              Clean,simplify, track clipped output  *
  56.  *                              and total bits/frame received.        *
  57.  * 7/10/91  Earle Jennings      changed to floats to FLOAT            *
  58.  *10/ 1/91  S.I. Sudharsanan,   Ported to IBM AIX platform.           *
  59.  *          Don H. Lee,                                               *
  60.  *          Peter W. Farrett                                          *
  61.  *10/ 3/91  Don H. Lee          implemented CRC-16 error protection   *
  62.  *                              newly introduced functions are        *
  63.  *                              buffer_CRC and recover_CRC_error      *
  64.  *                              Additions and revisions are marked    *
  65.  *                              with "dhl" for clarity                *
  66.  * 2/11/92  W. Joseph Carter    Ported new code to Macintosh.  Most   *
  67.  *                              important fixes involved changing     *
  68.  *                              16-bit ints to long or unsigned in    *
  69.  *                              bit alloc routines for quant of 65535 *
  70.  *                              and passing proper function args.     *
  71.  *                              Removed "Other Joint Stereo" option   *
  72.  *                              and made bitrate be total channel     *
  73.  *                              bitrate, irrespective of the mode.    *
  74.  *                              Fixed many small bugs & reorganized.  *
  75.  **********************************************************************
  76.  *                                                                    *
  77.  *                                                                    *
  78.  *  MPEG/audio Phase 2 coding/decoding multichannel                   *
  79.  *                                                                    *
  80.  *  Version 1.0                                                       *
  81.  *                                                                    *
  82.  *  7/27/93        Susanne Ritscher,  IRT Munich                      *
  83.  *                                                                    *
  84.  *                  thanks to                                         *
  85.  *                  Ralf Schwalbe,    Telekom FTZ Berlin              *
  86.  *                  Heiko Purnhagen,  Uni Hannover                    *
  87.  *                                                                    *
  88.  *  Version 2.0                                                       *
  89.  *                                                                    *
  90.  *  8/27/93        Susanne Ritscher, IRT Munich                       *
  91.  *                 Channel-Switching is working                       *
  92.  *                                                                    *
  93.  *  Version 2.1                                                       *
  94.  *                                                                    *
  95.  *  9/1/93         Susanne Ritscher,  IRT Munich                      *
  96.  *                 all channels normalized                            *
  97.  *                                                                    *
  98.  *  Version 3.0                                                       *
  99.  *                                                                    *
  100.  *  06/16/94       Ralf Schwalbe, Telekom FTZ Berlin                  *
  101.  *                 all sources and variables adapted due to MPEG-2 -  *
  102.  *                 DIS from March 1994                                *
  103.  *                  - dematrix and denormalize procedure              *
  104.  *                  - new tc - allocation (0-7)                       *
  105.  *                  - some new structures and variables as a basis    *
  106.  *                    for further decoding modes                      *
  107.  **********************************************************************
  108.  *       *
  109.  *  Version 1.0                                                       *
  110.  *                                                                    *
  111.  *  11/04/94       Ralf Schwalbe,  Telekom FTZ Berlin                 *
  112.  *                  - decoding tc-allocation                          *
  113.  *                  - some new subroutines, globale variables and     *
  114.  *                    structures (important to handle the ext. bitst.)*
  115.  *                  - changed all functions to ANSI-C funktion header *
  116.  *     - corrected some bugs to decode bitstreams > 512kB*
  117.  *       *
  118.  *  Version 1.1                                                       *
  119.  *                                                                    *
  120.  *  12/07/94       Ralf Schwalbe,  Telekom FTZ Berlin                 *
  121.  *                  - decoding extension bitstream                    *
  122.  *                                                                    *
  123.  *  Version 1.1.1                                                     *
  124.  *                 Ralf Schwalbe, Telekom FTZ Berlin                  *
  125.  *                  - fixed some bugs                                 *
  126.  *                                                                    *
  127.  *  Version 1.2                                                       *
  128.  *                                                                    *
  129.  *  6/21/95       Ralf Schwalbe, Deutsche Telekom FTZ Berlin          *
  130.  *                 - decoding dynamic crosstalk                      *
  131.  *                  - decoding phantom center                         * 
  132.  *                  - decoding MPEG1 compatible part only (stereo)    *
  133.  *                  - corrected some settings and bugs                *
  134.  *        *
  135.  *  7/12/95   Ralf Schwalbe, Deutsche Telekom FTZ Berlin          *
  136.  *     - corrected dynamic crosstalk        *
  137.  *     - 3/2,3/1,3/0,2/1, channel configurations         *
  138.  *                    are working                                     *
  139.  *          *
  140.  * 10/31/95   Ralf Schwalbe, Deutsche Telekom FTZ Berlin       *
  141.  *     - decoding of LFE-channel is working         *
  142.  *     - corrected any settings and bugs       *
  143.  *     - corrected table-switch for channel mode < 3/2   *
  144.  *       and tc-allocation / dyn-crosstalk       *
  145.  *     - still a problem with compl. bitstream 18/19     *
  146.  * 01/12/96   Ralf Schwalbe, Deutsche Telekom TZ Darmstadt       *
  147.  *     - decoder prediction installed       *
  148.  *     - problem with compl. bitstream 19 solved,        *
  149.  * 03/27/96   Ralf Schwalbe, Deutsche Telekom TZ Darmstadt        *
  150.  *     - Multilingual mode adapted for max. 7 channels   *
  151.  *                                                                    *       *     
  152.  *  Version 2.0                                                       *
  153.  *                                                                    *
  154.  * 01/28/97       Frans de Bont, Philips Sound & Vision, Eindhoven    *
  155.  *     - simultaneous use of tc allocation and dynamic   *
  156.  *       crosstalk working for all configurations       *
  157.  *     - prediction working for all configurations       *
  158.  *     - layer 1 MC working       *
  159.  *     - variable bit rate and extension bitstreams      *
  160.  *     - fully compliant to 13818-3.2                    *
  161.  *                                                                    *       *     
  162.  **********************************************************************/
  163. /**********************************************************************
  164.  *                                                                    *
  165.  *  06/06/95       Yeon Bae Thomas Kim,  Samsung AIT                  *
  166.  *                 ancillary data is working                          *
  167.  *       *
  168.  *  06/06/95       Sang Wook Kim,  Samsung AIT                        *
  169.  *                 corrected some bugs                                *
  170.  *       *
  171.  **********************************************************************/
  173. #include        "common.h"
  174. #include        "decoder.h"
  176. /********************************************************************
  177. *
  178. *        This part contains the MPEG I / II decoder for Layers II.
  179. *
  180. *        Core of the Layer II decoder.  Default layer is Layer II.
  181. *
  182. *********************************************************************/
  184. /* Global variable definitions for "musicout.c" */
  186. char     *programName;
  187. int     Frame_Bits,Bitrate;
  188. double    S_freq;
  190. /* Implementations */
  192. main (int argc, char **argv)  /* R.S. 7 channels for ML */
  193. {
  194. typedef long PCM[7][3][SBLIMIT];
  195. PCM  *pcm_sample;
  196. PCM  *pcm_sample_ml; /* 10/03/1995 JMZ Multilingual*/
  197. typedef unsigned int SAM[7][3][SBLIMIT];
  198. SAM  *sample;
  199. SAM  *sample_ml; /* 10/03/1995 JMZ Multilingual*/
  200. typedef double FRA[7][SBLIMIT][3][12]; /* 7.10.93 R.S. mem - alloc for DOS */ 
  201. FRA  *fraction;
  202. FRA  *fraction_ml; /*JMZ 09/03/1995 Multilingual */
  203. typedef double FRA_HELP[12][7][3][SBLIMIT]; /* 10/03/1995 JMZ Multilingual*/
  204. FRA_HELP *fraction_help;
  205. FRA_HELP *fraction_help_ml;
  206. typedef double FRA_BUF[2][8][36+PREDDEL];   
  207. FRA_BUF  *pred_buf; 
  208. typedef double VE[7][HAN_SIZE];
  209. VE  *w;
  211. frame_params      fr_ps;
  212. layer    info;
  213. Bit_stream        bs_mpg, bs_ext, bs_mc;
  214. FILE              *musicout, *lfe;
  215. FILE              *musicout_ml; /* 10/03/1995 JMZ Multilingual*/
  216. unsigned long     sample_frames;
  217. unsigned long     sample_frames_ml;
  218. int   i,jj, j, k, ii, stereo, done=FALSE, clip, sync,f;
  219. int               error_protection, crc_check = 1;
  220. int   crc_error_count, total_error_count;
  221. int               crc_error_count_mc, total_error_count_mc;
  222. #ifdef Augmentation_7ch
  223. int               crc_error_count_aug, total_error_count_aug;
  224. #endif
  225. unsigned int      old_crc, new_crc;
  226. unsigned int      I_bit_alloc[3][7][SBLIMIT];
  227. unsigned int      bit_alloc[7][SBLIMIT], scfsi[7][SBLIMIT],
  228.   scale_index[7][3][SBLIMIT]; 
  229. unsigned int      bit_alloc_ml[7][SBLIMIT], scfsi_ml[7][SBLIMIT],
  230.   scale_index_ml[7][3][SBLIMIT]; /* 09/03/1995 JMZ Multilingual */
  231. unsigned long     bitsPerSlot, samplesPerFrame;
  232. IFF_AIFF          pcm_aiff_data;
  233. char           encoded_file_name[MAX_NAME_SIZE];
  234. char           encoded_file_name1[MAX_NAME_SIZE]; /* 8/11/92.sr*/
  235. char              decoded_file_name[MAX_NAME_SIZE];
  236. char              decoded_file_name_ml[MAX_NAME_SIZE];/* 10/03/1995 JMZ Multilingual*/
  237. char   ext_bitstream_name[MAX_NAME_SIZE];
  238. char   frame_name[MAX_NAME_SIZE];
  239. char   lfe_file_name[MAX_NAME_SIZE];
  240. char              t[50];
  241. int               need_aiff;
  242. int   l, m, print_out = 0, ml =0;
  243. int   ch_start,lfe_init = 1;
  244. int   tca_log = 0, dynx_log = 0, bits_log = 0, rate_log = 0, scfsi_log = 0, any_log = 0;
  245. short int   lfe_tmp[12];
  246. unsigned int   sym_sample_ml[12][7][3][SBLIMIT]; /* 18/03/1996 FDB Multilingual*/
  247. int   part, layer_I_frames;
  248. int   mpeg = 2; /* R.S. mpeg default for mc */
  249. int   channels=2, mc_channel=0;
  250. unsigned long   frameNum = 0L;
  251. unsigned long   frameMod = 1L;
  252. unsigned long   frameBits;
  253. /*****************************************************************************/
  254. int hi, hu, ho;
  255. FILE *fp1;
  256. FILE *fp2;
  257. FILE *fp3;
  258. FILE *fp4;
  259. FILE *fp5;
  260. FILE *fp6;
  261. /******************************************************************************/
  263. #ifdef  MACINTOSH
  264. console_options.nrows = MAC_WINDOW_SIZE;
  265. argc = ccommand(&argv);
  266. #endif
  267. /* Most large variables are declared dynamically to ensure
  268.    compatibility with smaller machines */
  270. pcm_sample = (PCM *) mem_alloc((long) sizeof(PCM), "PCM Samp");
  271. pcm_sample_ml = (PCM *) mem_alloc((long) sizeof(PCM), "PCM Samp");
  272. sample     = (SAM *) mem_alloc((long) sizeof(SAM), "Sample");
  273. sample_ml   = (SAM *) mem_alloc((long) sizeof(SAM), "Sample");
  274. fraction   = (FRA *) mem_alloc((long) sizeof(FRA), "fraction");  /* R.S. */
  275. fraction_ml   = (FRA *) mem_alloc((long) sizeof(FRA), "fraction");  /* R.S. */
  276. fraction_help = (FRA_HELP *) mem_alloc((long) sizeof(FRA_HELP), "fraction");
  277. fraction_help_ml = (FRA_HELP *) mem_alloc((long) sizeof(FRA_HELP), "fraction");
  278. pred_buf  = (FRA_BUF *) mem_alloc((long) sizeof(FRA_BUF), "pred_buf");
  280. w = (VE *) mem_alloc((long) sizeof(VE), "w");
  282. bs_mpg.header_size = 0;
  283. bs_mpg.bits = (unsigned char *) mem_alloc ((long) 13824, "MPG Bits");
  284. bs_ext.header_size = 0;
  285. bs_ext.bits = (unsigned char *) mem_alloc ((long) 16376, "EXT Bits");
  286. bs_mc.bits  = (unsigned char *) mem_alloc ((long) 30208, "MC Bits");
  288. fr_ps.header = &info;
  289. fr_ps.bs_mpg = &bs_mpg;
  290. fr_ps.bs_ext = &bs_ext;
  291. fr_ps.bs_mc  = &bs_mc;
  292. fr_ps.tab_num = -1;                /* no table loaded */
  293. fr_ps.tab_num_mc = -1;                /* no table loaded */
  294. fr_ps.tab_num_ml = -1;
  295. fr_ps.alloc = NULL;
  296. fr_ps.alloc_mc = NULL;
  297. fr_ps.alloc_ml = NULL;
  299. info.mode_ext = 0;
  300. info.version = MPEG_AUDIO_ID;
  301. info.bitrate_index = 0;
  302. info.lfe = 0;         /* no low frequency effect channel present! */
  304. info.no_of_multi_lingual_ch = 0;
  305. info.multi_lingual_fs = 0;
  306. info.multi_lingual_layer = 0;
  308. info.ext_bit_stream_present = 0;
  309. info.ext_length = 0;
  310. info.n_ad_bytes = 0;    /* 7.12.94 R.S. */
  312. info.lfe = 0;
  314. for (i=0;i<HAN_SIZE;i++) for (j=0;j<5;j++) (*w)[j][i] = 0.0;
  316. program_information();
  317. programName = argv[0];
  318. if(argc==1)
  319. {        /* no command line args -> interact */
  320.    do
  321.    {
  322.   printf ("Enter encoded file name <required>: ");
  323.   gets (encoded_file_name);
  324.   f = strlen(encoded_file_name)-4;   /*cut off extension.8/11/92.sr*/
  325.   if (encoded_file_name[0] == NULL_CHAR)
  326.  printf ("Encoded file name is required. n");
  327.    } while (encoded_file_name[0] == NULL_CHAR);
  329.    printf (">>> Encoded file name is: %s n", encoded_file_name);
  330.    strcpy(encoded_file_name1, encoded_file_name);  /*8/11/92.sr*/
  331.    strcpy(&encoded_file_name1[f], DFLT_OPEXT_DEC);   /*.dec-extension.8/11/92.sr*/
  332.    printf ("Enter MPEG decoded file name <%s>: ", encoded_file_name1);
  333.    gets (decoded_file_name);
  334.    if (decoded_file_name[0] == NULL_CHAR)
  335. strcpy(decoded_file_name, encoded_file_name1);
  337. /* JMZ 10/03/1995 Multilingual */
  339.    printf (">>> MPEG decoded file name is: %s n", decoded_file_name);
  341.            /* encoded_file_name1 only for help */
  342.    strcpy(encoded_file_name1, encoded_file_name);
  343.    strcpy(&encoded_file_name1[f], DFLT_IPEXT_EXT);   /* .ext */
  344.    printf (">>> Enter MPEG 2 decoded extension filename <%s>: ", encoded_file_name1);
  345.    gets(ext_bitstream_name);
  346.    if( ext_bitstream_name[0] == NULL_CHAR )
  347. strcpy(ext_bitstream_name, encoded_file_name1);
  348.    printf("Extension bitstream <%s> will be decoded n",ext_bitstream_name);
  349.    /* 6/21/95 Ralf Schwalbe: stand a chance to decode MPEG1 compatible part */
  350.    if ((open_bit_stream_r (&bs_ext, ext_bitstream_name, BUFFER_SIZE)) == 0)
  351.    { 
  352. printf("No extension bitstream <%s> present n",ext_bitstream_name);
  353. printf("Do you want to decode an MPEG 2 bitstream ? (<y>/n) : ");
  354. gets(t);
  355.     if (*t == 'N' || *t == 'n') 
  356. {
  357. mpeg = 1;
  358. printf("An MPEG-1 Layer II bitstream will be decoded n");
  359. }
  360.     else
  361. {
  362. mpeg = 2;
  363. printf("An MPEG-2 Layer II bitstream will be decoded n");
  364. }
  365.    } 
  366.    printf("Do you wish to write an AIFF compatible sound file ? (<y>/n) : ");
  367.    gets(t);
  368.    if (*t == 'N' || *t == 'n') need_aiff = FALSE;
  369.    else                        need_aiff = TRUE;
  370.    if (need_aiff)
  371. printf(">>> An AIFF compatible sound file will be writtenn");
  372.    else printf(">>> A non-headered PCM sound file will be writtenn");
  373.    
  374.    printf("Do you want to print out all decoding information ? (y/<n>) :");
  375.    gets(t);
  376.    if( *t == 'y' || *t == 'Y' ) print_out = 1;
  377.    else        print_out = 0;
  379.    printf(
  380.   "Do you wish to exit (last chance before decoding) ? (y/<n>) : ");
  381.    gets(t);
  382.    if (*t == 'y' || *t == 'Y') exit(0);
  383. }
  385. else /************** interpret CL Args *****************/
  386. /* default: -write an AIFF file
  387.       -print encoded bitstream information  */
  388. {        
  389.    int i=0, err=0;
  391. #ifdef UNIX /* default values */
  392.    need_aiff = TRUE;
  393.    print_out = 0;   /* 10/31/95 Ralf Schwalbe: set print_out off in CL Args */
  394.    mpeg = 2;
  395. #else
  396.    need_aiff = FALSE;
  397. #endif
  398.    encoded_file_name[0] = NULL_CHAR;
  399.    decoded_file_name[0] = NULL_CHAR;
  400.    ext_bitstream_name[0] = NULL_CHAR;
  402.    while(++i<argc && err == 0)
  403.    {
  404.   char c, *token, *arg, *nextArg;
  405.   int  argUsed;
  407.   token = argv[i];
  408.   if(*token++ == '-')
  409.   {
  410.  if(i+1 < argc) nextArg = argv[i+1];
  411.  else           nextArg = "";
  412.  argUsed = 0;
  413.  while((c = *token++) != 0)
  414.  {
  415. if(*token /* NumericQ(token) */) arg = token;
  416. else                             arg = nextArg;
  417. switch(c)
  418. {
  419.    case 'a':  need_aiff = FALSE;
  420.       break;
  422.    case 'b':  bits_log = 1;
  423.       break;
  424.     
  425.    case 'd':  dynx_log = 1;
  426.       break;
  427.     
  428.    case 'h':  bs_mpg.header_size = 2304;
  429.       bs_ext.header_size = 2304;
  430.       break;
  432.    case 'i':  print_out = 1;
  433.       break;
  434.     
  435.    case 'm':  mpeg = 1;
  436.       break;
  437. #ifdef Augmentation_7ch
  438.    case 'n':  mpeg = 3;
  439.       break;
  440. #endif
  441.    case 'r':  rate_log = 1;
  442.       break;
  443.     
  444.    case 's':  scfsi_log = 1;
  445.       break;
  446.     
  447.    case 't':  tca_log = 1;
  448.       break;
  449.     
  450.    default:   fprintf(stderr,"%s: unrecognized option %cn",
  451.   programName, c);
  452.       err = 1;
  453.       break;
  454. }
  455. if(argUsed)
  456. {
  457. if(arg == token) token = ""; /* no more from token */
  458. else  ++i; /* skip arg we used */
  459. arg = ""; argUsed = 0;
  460. }
  461.  }
  462.   }
  463.   else
  464.   {
  465.      if (encoded_file_name[0] == NULL_CHAR)
  466.      {
  467. strcpy (encoded_file_name, argv[i]);
  468. f = strlen (encoded_file_name);
  469. strcpy (&encoded_file_name[f], DFLT_IPEXT_MPG);
  470. strcpy (ext_bitstream_name, encoded_file_name);
  471. strcpy (&ext_bitstream_name[f], DFLT_IPEXT_EXT);
  472.      }
  473.      else if (decoded_file_name[0] == NULL_CHAR)
  474.      {
  475. strcpy (decoded_file_name, argv[i]);
  476. f = strlen (decoded_file_name);
  477. strcpy (&decoded_file_name[f], DFLT_OPEXT_DEC);
  478.      }
  479.      else
  480.      {
  481. fprintf (stderr, "%s: excess arg %sn", programName, argv[i]);
  482. err = 1;
  483.      }
  484.   }
  485.        } /* end while argc */
  487.        if (bits_log || dynx_log || rate_log || scfsi_log || tca_log)
  488.   any_log = 1;
  491.        if (err || encoded_file_name[0] == NULL_CHAR)
  492.           usage ();  /* never returns */
  494.        if ((open_bit_stream_r (&bs_ext, ext_bitstream_name, BUFFER_SIZE)) == 0)
  495.        {
  496. printf("No extension bitstream <%s> present n", ext_bitstream_name);
  497. ext_bitstream_name[0] = NULL_CHAR;
  498.        }
  499.        
  500.        if (decoded_file_name[0] == NULL_CHAR)
  501.        {
  502.   f = strlen (encoded_file_name) - 4;
  503.   strcpy (decoded_file_name, encoded_file_name);
  504.   strcpy (&decoded_file_name[f], DFLT_OPEXT_DEC);
  505.        }
  507.     }  /* else arg-line end */
  509.     /* JMZ 10/03/1995 Multilingual */
  510.     strcpy(decoded_file_name_ml, decoded_file_name);
  511.     strcat(decoded_file_name_ml, DFLT_OPEXT_ML);
  513. /* report results of dialog / command line */
  514.     printf("Input file = '%s' nOutput file = '%s' ",
  515.    encoded_file_name, decoded_file_name);
  516.     if (mpeg >= 2 && ext_bitstream_name[0] != NULL_CHAR)
  517.      printf("Ext_bitstream = '%s'n", ext_bitstream_name);
  518.     if(mpeg == 1)
  519.     printf(">>>nA MPEG 1 Layer-1/2 bitstream will be decodedn");
  520.    else
  521. printf(">>>nA MPEG 2 Layer-1/2 bitstream will be decodedn");  
  522.   
  523.     if(need_aiff) printf("Output file written in AIFF formatn");
  525.     if ((musicout = fopen(decoded_file_name, "w+b")) == NULL)
  526.     {
  527.    fprintf (stderr, "Could not create "%s".n", decoded_file_name);
  528.    exit (1);
  529.     }
  530.     if (open_bit_stream_r (&bs_mpg, encoded_file_name, BUFFER_SIZE) == 0)
  531.     {
  532. fprintf (stderr, "Could not open "%s".n", encoded_file_name);
  533. exit (1);
  534.     }
  536.     if (need_aiff)
  537.   if (aiff_seek_to_sound_data(musicout) == -1)
  538.   {
  539.   fprintf (stderr, "Could not seek to PCM sound data in "%s".n",
  540.  decoded_file_name);
  541.   exit (1);
  542.    }
  544.     sample_frames = 0;
  545.     sample_frames_ml = 0;
  546.     total_error_count = total_error_count_mc = 0;
  547. #ifdef Augmentation_7ch
  548.     total_error_count_aug = 0;
  549. #endif
  551.     while (!end_bs (&bs_mpg))
  552.     {
  553.     /****************************************************************/
  554.     /*                                                              */
  555.     /*                  START OF FRAME LOOP                         */
  556.     /*                                                              */
  557.     /****************************************************************/
  559.    sync = seek_sync_mpg (&bs_mpg);
  560.    frameBits = bs_mpg.totbits;
  561.    if (!sync)
  562.    {
  563.   printf("nFrame cannot be located, end of input stream n");
  564.   done = TRUE;
  565.   /* finally write out the buffer */
  566.   out_fifo (*pcm_sample, 3, &fr_ps, done, musicout, &sample_frames, channels);
  567.   out_fifo_ml (*pcm_sample_ml, 3, &fr_ps, done, musicout_ml, &sample_frames_ml);
  568.   break;
  569.    }
  571.    info.lfe = 0;   /* R.S. ! */
  572.    decode_info (&bs_mpg, &fr_ps);
  573.    hdr_to_frps (&fr_ps);
  575.    stereo = fr_ps.stereo;
  576.    error_protection = info.error_protection;
  577.    crc_error_count = 0;
  578.    if (frameNum == 0) WriteHdr (&fr_ps, stdout);  /* printout layer/mode */
  580.    if (frameNum % frameMod == 0)
  581.    {
  582. fprintf(stderr, "{%4lu}r", frameNum++); fflush(stderr);
  583. if (frameNum / frameMod >= 100) frameMod *= 10;
  584.    }
  585.    else
  586. frameNum++;
  588.    Frame_Bits = 144*Bitrate/S_freq; 
  589.    if(info.padding)
  590. Frame_Bits++;    /* R.Schwalbe add a padding byte */
  591.    Frame_Bits*=8;  /* R.Schwalbe to get bits per frame *8 */
  592.    if (error_protection) buffer_CRC(&bs_mpg, &old_crc);
  593.           
  595.            /* 7/8/95 added message WtK 7/8/95 */
  596.    if (info.lay == 3)
  597.    {
  598.      fprintf (stderr, "nLayer = %d, not supported.nexit.nn",info.lay); break;
  599.      exit (1);
  600.    }
  601.    
  602.    bs_mc.totbits = bs_mc.curpos = 0;
  604.    if (info.lay == 1)
  605.    {
  606.       if (mpeg > 1)
  607.  layer_I_frames = 3;
  608.       else
  609.         layer_I_frames = 1;
  611.       for (part = 0; part < layer_I_frames; part++)
  612.       {
  613.  bitsPerSlot = 32;
  614.  samplesPerFrame = 384;
  615.  ch_start = 0;
  616.  channels = stereo;
  618.  if (part > 0)
  619.  {
  620.     sync = seek_sync_mpg (&bs_mpg);
  621.     frameBits = bs_mpg.totbits;
  622.     if (!sync)
  623.     {
  624.    fprintf (stderr, "2nd or 3th Frame cannot be located, end of input streamn");
  625.    exit (1);
  626.     }
  628.     decode_info (&bs_mpg, &fr_ps);
  629.     hdr_to_frps (&fr_ps);
  630.  
  631.     stereo = fr_ps.stereo;
  632.     error_protection = info.error_protection;
  633.     crc_error_count = 0;
  635.     Frame_Bits = 12*Bitrate/S_freq; 
  636.     if (info.padding)
  637.  Frame_Bits++;    /* R.Schwalbe add a padding byte */
  638.     Frame_Bits *=32;  /* R.Schwalbe to get bits per frame *32 */
  639.     if (error_protection)
  640.        buffer_CRC (&bs_mpg, &old_crc);
  641.           
  642.     /* 7/8/95 added message WtK 7/8/95 */
  643.     if (info.lay != 1)
  644.     {
  645.       fprintf (stderr, "Change of Layer = %d, not supported.nexit.nn", info.lay);
  646.       exit (1);
  647.     }
  648.  }
  650.  I_decode_bitalloc (&bs_mpg, &fr_ps, I_bit_alloc[part], bits_log);
  651.  I_decode_scale (&bs_mpg, &fr_ps, I_bit_alloc[part], scale_index, part, scfsi_log);
  653.  if (error_protection)
  654.  {
  655. I_CRC_calc (&fr_ps, I_bit_alloc[part], &new_crc);
  656. if (new_crc != old_crc)
  657. {
  658.    printf(" n ERROR in LAYER 1 - CRC! n");
  659.    crc_error_count++;
  660.    total_error_count++;
  661.    recover_CRC_error (*pcm_sample, crc_error_count,
  662.  &fr_ps, musicout, &sample_frames, channels);
  663. }
  664. else
  665.    crc_error_count = 0;
  666.  }
  668.  clip = 0;
  670.  for (i=0;i<12;i++)
  671.  {
  672. I_buffer_sample (&bs_mpg, &fr_ps, (*sample), I_bit_alloc[part]);
  673. if (!any_log)
  674. {
  675. I_dequantize_sample ((*sample), part, I_bit_alloc[part], *fraction, &fr_ps, &i);
  676. I_denormalize_sample (*fraction, scale_index, part, &fr_ps, &i);
  677. }
  678.  }   /* end of for loop */
  680.  if (mpeg > 1)
  681.  {
  682.     while (bs_mpg.curpos < bs_mpg.totbits-1)
  683.        bs_mc.bits[bs_mc.totbits++] = bs_mpg.bits[bs_mpg.curpos++];
  685.     if (part == 0)
  686.     {
  687.        mc_header (&bs_mc, &fr_ps);
  688.        mc_hdr_to_frps (&fr_ps);
  689.     }
  691.     bs_mc.totbits -= fr_ps.header->n_ad_bytes * 8;
  692.  }
  693.       }
  694.       if (mpeg > 1)
  695.       {
  696.  if (fr_ps.header->ext_bit_stream_present)
  697.  {
  698.     sync = seek_sync_ext (&bs_ext, &fr_ps);
  699.     frameBits += bs_ext.totbits;
  700.     if (!sync)
  701.     {
  702.        fprintf(stderr, "Extension frame cannot be located, end of input streamn");
  703.        exit (1);
  704.     }
  706.     while (bs_ext.curpos < bs_ext.totbits)
  707.        bs_mc.bits[bs_mc.totbits++] = bs_ext.bits[bs_ext.curpos++];
  708.  }
  710.  for (i=0; i < channels; i++)
  711.     for (j=0; j < SBLIMIT; j++)
  712.        if (I_bit_alloc[0][i][j] > 0)
  713.           bit_alloc[i][j] = I_bit_alloc[0][i][j];
  714.        else if (I_bit_alloc[1][i][j] > 0)
  715.           bit_alloc[i][j] = I_bit_alloc[1][i][j];
  716.        else if (I_bit_alloc[2][i][j] > 0)
  717.           bit_alloc[i][j] = I_bit_alloc[2][i][j];
  718.        else
  719.           bit_alloc[i][j] = 0;
  720.       }
  721.    }
  722.    else if (info.lay == 2 )
  723.    {
  724.       layer_I_frames = 3;
  725.       bitsPerSlot = 8;
  726.       samplesPerFrame = 1152;
  727.       ch_start = 0;
  728.       channels = stereo;
  729.       II_decode_bitalloc (&bs_mpg, &fr_ps, bit_alloc, bits_log);
  730.       II_decode_scale (&bs_mpg, &fr_ps, scfsi, bit_alloc, scale_index, &ch_start, &channels, scfsi_log);
  732.       if (error_protection)
  733.       {
  734.  II_CRC_calc(&fr_ps, bit_alloc, scfsi, &new_crc);
  735.  if (new_crc != old_crc)
  736.  {
  737.     printf(" n ERROR in LAYER 2 - CRC! n");
  738.     crc_error_count++;
  739.     total_error_count++;
  740.     recover_CRC_error (*pcm_sample, crc_error_count,
  741.        &fr_ps, musicout, &sample_frames, channels);
  742.  }
  743.  else crc_error_count = 0;
  744.       }
  746.       clip = 0;
  748.       for (i=0;i<12;i++)
  749.       {
  750.  II_buffer_sample (&bs_mpg, &fr_ps, (*sample), bit_alloc);
  751. if (!any_log)
  752. {
  753.  II_dequantize_sample ((*sample), bit_alloc, *fraction, &fr_ps, &i);
  754.  II_denormalize_sample (*fraction, scale_index, &fr_ps, i>>2, &i);
  755. }
  756.       }   /* end of for loop */
  758.       if (mpeg > 1)
  759.       {
  760.  while (bs_mpg.curpos < bs_mpg.totbits)
  761.     bs_mc.bits[bs_mc.totbits++] = bs_mpg.bits[bs_mpg.curpos++];
  763.  mc_header (&bs_mc, &fr_ps);
  764.  mc_hdr_to_frps (&fr_ps);
  766.  bs_mc.totbits -= fr_ps.header->n_ad_bytes * 8;
  768.  if (fr_ps.header->ext_bit_stream_present)
  769.  {
  770.     sync = seek_sync_ext (&bs_ext, &fr_ps);
  771.     frameBits += bs_ext.totbits;
  772.     if (!sync)
  773.     {
  774.        fprintf(stderr, "Extension frame cannot be located, end of input streamn");
  775.        exit (1);
  776.     }
  778.     while (bs_ext.curpos < bs_ext.totbits)
  779.        bs_mc.bits[bs_mc.totbits++] = bs_ext.bits[bs_ext.curpos++];
  780.  }
  781.       }
  782.    }
  784.    /**********************************************************/
  785.    /*      multichannel - decoding                */
  786.    /*             7.07.93 Susanne Ritscher                   */
  787.    /*            13.10.93 Ralf Schwalbe                      */
  788.    /*        30.05.94 Ralf Schwalbe    */
  789.    /**********************************************************/
  791.    if (rate_log)
  792.      printf ("frame: %5d bitrate: %8.3f kbpsn", frameNum, frameBits * S_freq / samplesPerFrame);
  793.    else if (mpeg >= 2)
  794.    {
  795.      if ( (frameNum==1) && (fr_ps.header->no_of_multi_lingual_ch>0) )
  796.      {
  797.        /*frameNum==1 NOT frameNum==0*/
  798.        if (fr_ps.header->no_of_multi_lingual_ch>7)
  799.        {
  800.  printf("nn***Warning. There are %d multilingual channels.n", fr_ps.header->no_of_multi_lingual_ch);
  801.  printf    ("            This decoder supports only upto 7 channels.nn");
  802.        }
  803.        printf("There are %d multilingual channels, coded at ",fr_ps.header->no_of_multi_lingual_ch);
  804.        if (fr_ps.header->multi_lingual_fs==0) printf("full"); else printf("half");
  805.        printf(" sampling rate.n They are outputted in "%s".n",decoded_file_name_ml);
  806.        if ((musicout_ml = fopen(decoded_file_name_ml, "w+b")) == NULL)
  807.        {
  808.  fprintf (stderr, "***Error. Could not create "%s".nExitn", decoded_file_name_ml);
  809.  exit (1);
  810.        }
  811.        if(need_aiff)
  812.        {
  813.  printf("Multilingual file written in AIFF formatn");
  814.  if (aiff_seek_to_sound_data (musicout_ml) == -1)
  815.  {
  816.     fprintf (stderr, "Could not seek to ML PCM sound data in "%s".n",
  817.     decoded_file_name_ml);
  818.     exit (1);
  819.  }  
  820.        }
  821.      }
  823.      mc_channel = fr_ps.mc_channel;
  824.      crc_error_count_mc = 0;
  825.      buffer_CRC (&bs_mc, &old_crc);   /* read CRC - check from header */
  826.      mc_composite_status_info (&bs_mc, &fr_ps, tca_log, dynx_log);
  827.      
  828.      if (frameNum == 1) 
  829.     if (info.lfe)
  830.     printf ("mc-channel=%d, MPEG2-ext-length=%d kbit/s lfe-channel presentn",
  831.      mc_channel, (int) (info.ext_length * s_freq[info.sampling_frequency] / 144));
  832.   else
  833.     printf ("mc-channel=%d, MPEG2-ext-length=%d kbit/sn",
  834.      mc_channel, (int) (info.ext_length * s_freq[info.sampling_frequency] / 144));
  836.      /* Achtung: fr_ps.jsbound = 27; wird in Abh鋘gikeit von fs gesetzt*/
  837.      ch_start = stereo;
  838.      channels = stereo + mc_channel;
  840.      /* decoding bitallocation's */
  841.    
  842.      II_decode_bitalloc_mc (&bs_mc, &fr_ps, bit_alloc, &ch_start, &channels, bits_log);
  843.      II_decode_scale (&bs_mc, &fr_ps, scfsi, bit_alloc, scale_index, &ch_start, &channels, scfsi_log);
  845.      /* mandatory CRC check is set */
  846.      if (crc_check)
  847.      {
  848. mc_error_check (&fr_ps, bit_alloc, scfsi, &new_crc, ch_start, channels);
  849. if (new_crc != old_crc)
  850. {
  851.    crc_error_count_mc++;
  852.    total_error_count_mc++;
  853.    for (i = 0; i < SBLIMIT; i ++)
  854.       for (ii = ch_start; ii < channels; ++ii)
  855.  bit_alloc[ii][i] = 0;
  856.    printf ("nERROR in MC-CRC -> mc frame can't be decoded !n");
  857. }
  858. else
  859.    crc_error_count_mc = 0;
  860.      }  /* endif crc_check */
  862.      clip = 0;
  863.      for (i = 0; i < 12; i++)
  864.      {
  865. II_buffer_sample_mc (&bs_mc, &fr_ps, (*sample), bit_alloc, ch_start, channels, i);
  866. if (!any_log)
  867. {
  868. II_dequantize_sample_mc ((*sample), bit_alloc, *fraction, &fr_ps, ch_start, channels, &i);
  869. II_denormalize_sample_mc (*fraction, scale_index, &fr_ps, i>>2, ch_start, channels, &i);
  870. }
  871.      }
  873.      /* 10/31/95 Ralf Schwalbe LFE */
  874.      if( info.lfe ) 
  875.      { 
  876. if(lfe_init)
  877. {
  878.    f = strlen (encoded_file_name) - 4;
  879.    strcpy (lfe_file_name, encoded_file_name);
  880.    strcpy (&lfe_file_name[f], DFLT_OPEXT_LFE);
  881.    if ((lfe = fopen (lfe_file_name, "w+b")) == NULL)
  882.    {
  883.       fprintf (stderr, "Could not create "%s".n",lfe_file_name );
  884.       exit (1);
  885.    }
  886.    if (need_aiff)
  887.    {
  888.       printf ("LFE file written in AIFF formatn");
  889.       if (aiff_seek_to_sound_data (lfe) == -1)
  890.       {
  891.  fprintf (stderr, "Could not seek to LFE PCM sound data in "%s".n",
  892.  lfe_file_name);
  893.  exit (1);
  894.       }  
  895.    }
  896. }
  897. lfe_init = 0; /* R.S. only one time */
  898. II_lfe_calc(&fr_ps);  
  899.      }
  900. /*********************************************************************/
  901. /* JMZ 09/03/1995 Multilingual , WtK 7/8/95 */
  903.      if (info.no_of_multi_lingual_ch>0) 
  904.      {
  905. /* SWKim060695 add the condion with 32, 44.1kHz case */
  906. if (info.sampling_frequency == 1)
  907.    fr_ps.jsbound = 27;
  908. else
  909.    fr_ps.jsbound = 30;
  911. m = info.no_of_multi_lingual_ch;
  913. II_decode_bitalloc_ml (&bs_mc, &fr_ps, bit_alloc_ml, &m);
  914. II_decode_scale_ml (&bs_mc, &fr_ps, scfsi_ml, bit_alloc_ml, scale_index_ml, &m);
  917.    clip = 0;
  918.    for (i=0; i< ( (info.multi_lingual_fs==0) ? 12 : 6 ); i++) /* ngr : n_of_granules loop */
  919.    {
  920.       II_buffer_sample_ml (&bs_mc, &fr_ps, (*sample_ml), bit_alloc_ml, &m);
  921. if (!any_log)
  922. {
  923.       II_dequantize_sample_ml ((*sample_ml), bit_alloc_ml, *fraction_ml, &fr_ps, &m, &i);
  924.       II_denormalize_sample_ml (*fraction_ml, scale_index_ml, &fr_ps,
  925.       ( (info.multi_lingual_fs==0) ? (i>>2) : (i>>1) ), &m, &i);
  926. }
  927.       for (k = 0; k < info.no_of_multi_lingual_ch; k++)
  928.  for (ii = 0; ii < 3; ii++)
  929.     for (j = 0; j < 32; j++)
  930.        sym_sample_ml[i][k][ii][j] = (*sample_ml)[k][ii][j];
  931.    }
  932.      }
  934. /* JMZ 09/03/1995 Multilingual */
  935. /*********************************************************************/
  937. #ifdef Augmentation_7ch
  938.      if (mpeg == 3)
  939.        if (stereo != 2 || fr_ps.header->surround != 2 || fr_ps.header->center == 0)
  940.        {
  941.  /* no 3/2 mode */
  942.  printf ("WARNING: Stream does not contain 5/2 infon");
  943.  printf ("         Continue decoding of MPEG2 MC partn");
  944.  mpeg = 2;
  945.        }
  947.      if (mpeg == 3)
  948.      {
  949.        crc_error_count_aug = 0;
  950.        mc_aug_composite_status_info (&bs_mc, &fr_ps, tca_log, dynx_log);
  952.        buffer_CRC (&bs_mc, &old_crc);   /* read CRC - check from header */
  953.        
  954.        /* Achtung: fr_ps.jsbound = 27; wird in Abh鋘gikeit von fs gesetzt*/
  955.        ch_start = 5;
  956.        channels = 7;
  958.        /* decoding bitallocation's */
  959.      
  960.        II_decode_bitalloc_aug (&bs_mc, &fr_ps, bit_alloc, &ch_start, &channels, bits_log);
  961.        II_decode_scale (&bs_mc, &fr_ps, scfsi, bit_alloc, scale_index, &ch_start, &channels, scfsi_log);
  963.        /* mandatory CRC check is set */
  964.        if (crc_check)
  965.        {
  966.      mc_aug_error_check(&fr_ps, bit_alloc, scfsi, &new_crc);
  967.      if (new_crc != old_crc)
  968.      {
  969.      crc_error_count_aug++;
  970.      total_error_count_aug++;
  971.      for( i = 0; i < SBLIMIT; i ++)
  972.      for(ii = ch_start; ii < channels; ++ii)
  973.      bit_alloc[ii][i] = 0;
  975.      printf("ERROR in MC-AUG-CRC -> mc aug frame can't be decoded !n");
  976.      }
  977.      else
  978.      crc_error_count_aug = 0;
  979.        }  /* endif crc_check */
  981.        clip = 0;
  982.        for (i=0;i<12;i++)
  983.        {
  984.      II_buffer_sample_aug (&bs_mc, &fr_ps, (*sample), bit_alloc, i);
  985. if (!any_log)
  986. {
  987.      II_dequantize_sample_aug ((*sample), bit_alloc, *fraction, &fr_ps, &i);
  988.      II_denormalize_sample_aug (*fraction, scale_index, &fr_ps, i>>2, &i);
  989. }
  990. }
  991.      }   /* endif mpeg == 3 */
  992. #endif
  994.      if (info.mc_prediction_on)
  995. for (k=0; k<2; k++)
  996.    for (ii=0; ii<8; ii++) 
  997.    {
  998.       for (j=0; j<PREDDEL; j++)
  999.  (*pred_buf)[k][ii][j] = (*pred_buf)[k][ii][j+36];
  1000.       for (j =0; j <3; j++)
  1001.  for (jj=0;jj<12;jj++)
  1002.     (*pred_buf)[k][ii][PREDDEL+j+3*jj] = (*fraction)[k][ii][j][jj];
  1003.    }
  1005.      /* prediction & dematricing MPEG2 part */
  1006. if (!any_log)
  1007.      if (crc_error_count_mc == 0)
  1008. dematricing_mc (*fraction, &fr_ps, *pred_buf);
  1010. #ifdef Augmentation_7ch
  1011.      if (mpeg == 3)
  1012.      {
  1013. if (!any_log)
  1014. if (crc_error_count_aug == 0)
  1015. {
  1016.    dematricing_aug (*fraction, &fr_ps);
  1017.    denormalizing_aug (*fraction, &fr_ps);
  1018. }
  1019.       }
  1020.       else
  1021. #endif
  1022.       if (crc_error_count_mc == 0)
  1023. denormalizing_mc (*fraction, &fr_ps, channels);
  1024.    }   /* endif mpeg >= 2 */
  1026. if (!any_log)
  1027.    for (jj=0;jj<4*layer_I_frames;jj++)
  1028.       for (k=0; k < channels; k++)
  1029.   for (j =0; j <3; j++)
  1030.      for (ii=0; ii< SBLIMIT; ii++)
  1031. (*fraction_help)[jj][k][j][ii] = (*fraction)[k][ii][j][jj];
  1032. if (!any_log)
  1033.    for (i=0;i<4*layer_I_frames;i++)
  1034.    {   
  1035.       for (j=0;j<3;j++) 
  1036.          for (k=0; k < channels; k++)
  1037.     clip += SubBandSynthesis (&((*fraction_help)[i][k][j][0]), k,
  1038.     &((*pcm_sample)[k][j][0]));
  1039.    
  1040.       out_fifo (*pcm_sample, 3, &fr_ps, done, musicout,
  1041.     &sample_frames, channels);
  1042.    }
  1044. /* 31/10/95 Ralf Schwalbe LFE output */
  1045. if (!any_log)
  1046.    if (info.lfe) 
  1047.    {
  1048.       for (jj = 0; jj < 12; jj++)
  1049.       {
  1050. /* lfe_tmp[jj] = (short int) (info.lfe_spl_fraction[jj] * SCALE); */
  1051. double foo = floor (info.lfe_spl_fraction[jj] * SCALE + 0.5);
  1053. if (foo >= (long) SCALE)      {lfe_tmp[jj] = SCALE-1; clip++;}
  1054. else if (foo < (long) -SCALE) {lfe_tmp[jj] =-SCALE;   clip++;}
  1055. else                           lfe_tmp[jj] = foo;
  1056.       }
  1057.       fwrite (&lfe_tmp[0], 2, 12, lfe); 
  1058.    }
  1059.    if (clip > 0)
  1060.       printf("n%d samples clippedn", clip);
  1062. /***************************************************************/
  1063. /* 10/03/1995 JMZ Multilingual , WtK 7/8/95 */
  1064.    if (info.no_of_multi_lingual_ch>0) 
  1065.    {
  1066.       for(jj=0;jj<( (info.multi_lingual_fs==0) ? 12 : 6 );jj++)
  1067.  for(k=0; k<info.no_of_multi_lingual_ch; k++)
  1068.      for(j =0; j <3; j++)
  1069. for(ii=0; ii< SBLIMIT; ii++)
  1070.       (*fraction_help_ml)[jj][k][j][ii] = (*fraction_ml)[k][ii][j][jj];
  1072.       clip = 0;
  1073.       for (i=0;i<( (info.multi_lingual_fs==0) ? 12 : 6 );i++)
  1074.       {   
  1075. for (j=0;j<3;j++) 
  1076.   for (k=0; k<info.no_of_multi_lingual_ch; k++)
  1077.       clip += SubBandSynthesis_ml(&((*fraction_help_ml)[i][k][j][0]), k,
  1078.        &((*pcm_sample_ml)[k][j][0]));
  1080. out_fifo_ml(*pcm_sample_ml, 3, &fr_ps, done, musicout_ml, &sample_frames_ml);
  1081.       }
  1082.       if(clip > 0) printf("n%d MultiLingual samples clippedn", clip);
  1083.    }
  1084. /* 10/03/1995 JMZ Multilingual */
  1085. /**************************************************************/
  1087. if( print_out == 1)
  1088. {
  1089. if(frameNum == 1) 
  1090. {
  1091. fp6 = fopen("SUB_ML","w+");
  1092. fp5 = fopen("DYN_CROSS","w+");
  1093. fp4 = fopen("INFODEC", "w+");
  1094. fp2 = fopen("BALDEC", "w+");
  1095. fp3 = fopen("SCFSIDEC", "w+");
  1096. fp1 = fopen("SCFDEC", "w+");
  1097. }
  1098. /*********************Ausgabe***************************************************/
  1099. fprintf(fp1, "******************** FRAME %d *****************************n", frameNum-1);
  1100. fprintf(fp2, "n******************** FRAME %d *****************************n", frameNum-1);
  1101. fprintf(fp3, "n******************** FRAME %d *****************************n", frameNum-1);
  1102. fprintf(fp4, "******************** FRAME %d *****************************n", frameNum-1);
  1103. fprintf(fp6, "******************** FRAME %d *****************************n", frameNum-1);
  1105. fprintf(fp4, "version = %d, ", info.version);
  1106. fprintf(fp4, "lay = %d, ", info.lay);
  1107. fprintf(fp4, "error_protection = %dn", info.error_protection);
  1109. fprintf(fp4, "bitrate_index = %d, ", info.bitrate_index);
  1111. fprintf(fp4, "sampling_frequency = %d,", info.sampling_frequency);
  1112. fprintf(fp4, "padding = %d, ", info.padding);
  1113. fprintf(fp4, "extension = %dn", info.extension);
  1115. fprintf(fp4, "mode = %d, ", info.mode);
  1116. fprintf(fp4, "mode_ext = %d,", info.mode_ext);
  1117. fprintf(fp4, "copyright = %d, ", info.copyright);
  1118. fprintf(fp4, "original = %d, ", info.original);
  1119. fprintf(fp4, "emphasis = %dn", info.emphasis);
  1121. fprintf(fp4, "center = %d, ", info.center);
  1122. fprintf(fp4, "surround = %d, ", info.surround);
  1123. fprintf(fp4, "dematrix_procedure = %d, ", info.dematrix_procedure);
  1124. fprintf(fp4, "lfe = %dn", info.lfe);
  1126. fprintf(fp4, "no_of_multi_lingual_ch = %d, ", info.no_of_multi_lingual_ch);
  1127. fprintf(fp4, "multi_lingual_fs = %d, ", info.multi_lingual_fs);
  1128. fprintf(fp4, "multi_lingual_layer = %dn", info.multi_lingual_layer);
  1130. fprintf(fp4, "dyn_cross_on = %dn", info.dyn_cross_on);
  1131. fprintf(fp4, "mc_prediction_on = %dn", info.mc_prediction_on);
  1133. #ifdef DEBUG_PREDICTION
  1134. if(info.mc_prediction_on == 1)
  1135.    for(hu = 0; hu < 8; ++ hu)
  1136.    {
  1137. fprintf(fp4, "mc_prediction[%d] = %d ", hu, info.mc_prediction[hu]);
  1138. for(ho = 0; ho < 4; ho++)
  1139. {
  1140. fprintf(fp4, "mc_predsi[%d][%d] = %d  ", hu,ho, info.mc_predsi[hu][ho]);
  1141. fprintf(fp4, "mc_delay[%d][%d] = %d  ", hu,ho, info.mc_delay_comp[hu][ho]);
  1142. fprintf(fp4, " n");
  1143. for(hi = 0;hi <3; hi++)
  1144. fprintf(fp4, "pred_coeff[%d][%d][%d] = %d ",hu,ho,hi,info.mc_pred_coeff[hu][ho][hi]);
  1145. fprintf(fp4, " n");
  1146. }
  1147.    }
  1148. #endif
  1150. if(info.dyn_cross_on == 1)
  1151. {
  1152. fprintf(fp5, "n******************** FRAME %d *****************************n", frameNum-1);
  1153. fprintf(fp5,"dyn_cross_LR = %dn",info.dyn_cross_LR);
  1154. for(ho = 0; ho < 12; ho ++)
  1155. {
  1156. fprintf(fp5,"dyn_cross[%d] = %d   ",ho,info.dyn_cross_mode[ho]);
  1157. if((ho+1)%4 == 0) fprintf(fp5,"n");
  1158. }
  1159. if (info.surround == 3)
  1160. for(ho = 0; ho < 12; ho ++)
  1161. {
  1162. fprintf(fp5,"dyn_second_stereo[%d] = %d   ",ho,info.dyn_second_stereo[ho]);
  1163. if((ho+1)%4 == 0) fprintf(fp5,"n");
  1164. }
  1165. }
  1167. fprintf(fp4, "tc_sbgr_select = %dn", info.tc_sbgr_select);
  1168. if(info.tc_sbgr_select == 1)
  1169.   fprintf(fp4, "tc_allocation = %dn", info.tc_allocation);
  1170. else
  1171. {
  1172. for(ho = 0; ho < 12; ho ++)
  1173. {
  1174. if(info.mc_prediction[ho] == 0)
  1175. fprintf(fp4, "tc_alloc[%d] = %dn", ho, info.tc_alloc[ho]);
  1176. else{
  1177. for(hu = 0; hu < 4; hu++)
  1178. fprintf(fp4, "predsi[%d][%d] = %dn", ho, hu, info.mc_predsi[ho][hu]);
  1179. }
  1180. }
  1181. fprintf(fp4, "n");
  1182. }
  1184. for(hu = 0; hu < SBLIMIT; ++hu)
  1185. {
  1186. fprintf(fp2, "n");
  1187. fprintf(fp3, "n");
  1189. for(hi = 0; hi < channels; hi ++)
  1190. {
  1191. fprintf(fp1, "n");
  1192. if (layer_I_frames == 1)
  1193.    fprintf(fp2, "BAL[%d][%02d] = %d   ", hi, hu, I_bit_alloc[0][hi][hu]);
  1194. else
  1195.    fprintf(fp2, "BAL[%d][%02d] = %d   ", hi, hu, bit_alloc[hi][hu]);
  1196. fprintf(fp3, "SCFSI[%d][%02d] = %d   ", hi, hu, scfsi[hi][hu]);
  1198. for( ho = 0; ho < layer_I_frames; ho ++)
  1199. {
  1200.    fprintf(fp1, "SCF[%d][%d][%02d] = %d   ", hi, ho, hu, scale_index[hi][ho][hu]);
  1201. }
  1202. }
  1203. for(hi = 0; hi < info.no_of_multi_lingual_ch; hi ++)
  1204. {
  1205.    fprintf(fp1, "n");
  1206.    /* 960315 added FDB Multi-Lingual bitalloc */
  1207.    fprintf(fp2, "BAL_ML[%d][%02d] = %d   ", hi, hu, bit_alloc_ml[hi][hu]);
  1208.    /* 960315 added FDB Multi-Lingual scf select info */
  1209.    fprintf(fp3, "SCFSI_ML[%d][%02d] = %d   ", hi, hu, scfsi_ml[hi][hu]);
  1210.    /* 960315 added FDB Multi-Lingual scfs */
  1211.    for( ho = 0; ho < 3; ho ++)
  1212.       fprintf(fp1, "SCF_ML[%d][%d][%02d] = %d   ", hi, ho, hu, scale_index_ml[hi][ho][hu]);
  1213.         }
  1214. fprintf(fp1, "n");
  1215. }
  1217. /* 960318 added FDB Multi-Lingual subband samples symbols*/
  1218. for(k=0; k<info.no_of_multi_lingual_ch; k++)
  1219.    for(ii=0; ii< SBLIMIT; ii++)
  1220.    {
  1221.       fprintf(fp6, "SB_SYM_ML[%d][%02d]: ", k, ii);
  1222.       for(jj=0;jj<( (info.multi_lingual_fs==0) ? 12 : 6 );jj++)
  1223.          for(j =0; j <3; j++)
  1224.             fprintf(fp6, " %6d", sym_sample_ml[jj][k][j][ii]);
  1225.       fprintf(fp6, "n");
  1226.    }
  1229. /* 960318 added FDB Multi-Lingual subband samples */
  1230. for(k=0; k<info.no_of_multi_lingual_ch; k++)
  1231.    for(ii=0; ii< SBLIMIT; ii++)
  1232.    {
  1233.       fprintf(fp6, "SB_ML[%d][%02d]: ", k, ii);
  1234.       for(jj=0;jj<( (info.multi_lingual_fs==0) ? 12 : 6 );jj++)
  1235.          for(j =0; j <3; j++)
  1236.             fprintf(fp6, " %12.5e", (*fraction_help_ml)[jj][k][j][ii]);
  1237.       fprintf(fp6, "n");
  1238.    }
  1239. fflush (fp1);
  1240. fflush (fp2);
  1241. fflush (fp3);
  1242. fflush (fp4);
  1243. fflush (fp5);
  1244. fflush (fp6);
  1245. /*******************************************************************************/
  1247. } /* endif of print_out == 1 */
  1249.     /****************************************************************/
  1250.     /*                                                              */
  1251.     /*                   END OF FRAME LOOP                          */
  1252.     /*                                                              */
  1253.     /****************************************************************/
  1255. }/*end of while(!endof(bs)) - loop */
  1257. if (need_aiff)
  1258. {
  1259. #ifdef Augmentation_7ch
  1260.    if (mpeg == 3)
  1261.       pcm_aiff_data.numChannels       = 7;
  1262.    else
  1263. #endif
  1264.       pcm_aiff_data.numChannels       = stereo + mc_channel;
  1265.    pcm_aiff_data.numSampleFrames   = sample_frames;
  1266.    pcm_aiff_data.sampleSize        = 16;
  1267.    pcm_aiff_data.sampleRate        = s_freq[info.sampling_frequency]*1000;
  1268.    strcpy(pcm_aiff_data.sampleType, IFF_ID_SSND);
  1269.    pcm_aiff_data.blkAlgn.offset    = 0;
  1270.    pcm_aiff_data.blkAlgn.blockSize = 0;
  1272.    if (aiff_write_headers(musicout, &pcm_aiff_data) == -1)
  1273.    {
  1274.   fprintf (stderr, "Could not write AIFF headers to "%s"n",
  1275.  decoded_file_name);
  1276.   exit (1);
  1277.            }
  1279.            if (fr_ps.header->no_of_multi_lingual_ch>0)
  1280.    {
  1281.      pcm_aiff_data.numChannels       = fr_ps.header->no_of_multi_lingual_ch;
  1282.      pcm_aiff_data.numSampleFrames   = sample_frames;
  1283.      pcm_aiff_data.sampleSize        = 16;
  1284.      pcm_aiff_data.sampleRate        = s_freq[info.sampling_frequency]*1000;
  1285.      if (fr_ps.header->multi_lingual_fs==1) pcm_aiff_data.sampleRate *= 0.5;
  1286.      strcpy(pcm_aiff_data.sampleType, IFF_ID_SSND);
  1287.      pcm_aiff_data.blkAlgn.offset    = 0;
  1288.      pcm_aiff_data.blkAlgn.blockSize = 0;
  1289.   
  1290.      if (aiff_write_headers(musicout_ml, &pcm_aiff_data) == -1)
  1291.      {
  1292.     fprintf (stderr, "Could not write AIFF headers to "%s"n",
  1293.    decoded_file_name_ml);
  1294.     exit (1);
  1295.      }
  1296.    }
  1297.    if (lfe_init == 0)
  1298.    {
  1299.      pcm_aiff_data.numChannels       = 1;
  1300.      pcm_aiff_data.numSampleFrames   = sample_frames / 96;
  1301.      pcm_aiff_data.sampleSize        = 16;
  1302.      pcm_aiff_data.sampleRate        = (s_freq[info.sampling_frequency]*1000) / 96;
  1303.      strcpy(pcm_aiff_data.sampleType, IFF_ID_SSND);
  1304.      pcm_aiff_data.blkAlgn.offset    = 0;
  1305.      pcm_aiff_data.blkAlgn.blockSize = 0;
  1307.      if (aiff_write_headers(lfe, &pcm_aiff_data) == -1)
  1308.      {
  1309.     fprintf (stderr, "Could not write AIFF headers to "%s"n",
  1310.    lfe_file_name);
  1311.     exit (1);
  1312.      }
  1313.    }
  1314.  }
  1316. close_bit_stream_r(&bs_mpg);
  1318. if (info.ext_bit_stream_present)
  1319. {
  1320.      close_bit_stream_r(&bs_ext);
  1322. }
  1323. fclose(musicout);
  1324. if(info.lfe) fclose(lfe);
  1327. #ifdef  MACINTOSH
  1328. if (need_aiff) set_mac_file_attr(decoded_file_name, VOL_REF_NUM,
  1329.  CREATR_DEC_AIFF, FILTYP_DEC_AIFF);
  1330. else           set_mac_file_attr(decoded_file_name, VOL_REF_NUM,
  1331.  CREATR_DEC_BNRY, FILTYP_DEC_BNRY);
  1332. #endif
  1334. printf("Decoding of "%s" with %lu frames is finishedn", encoded_file_name, frameNum);
  1335. printf("The decoded PCM output file name is "%s"n", decoded_file_name);
  1336. if (need_aiff)
  1337.    printf(""%s" has been written with AIFF header informationn",
  1338.   decoded_file_name);
  1340. if (fr_ps.header->no_of_multi_lingual_ch>0) {
  1341.           printf("There are %d multilingual channels, coded at ",fr_ps.header->no_of_multi_lingual_ch);
  1342.           if (fr_ps.header->multi_lingual_fs==0) printf("full"); else printf("half");
  1343.           printf(" sampling rate.n They are outputted in "%s".n",decoded_file_name_ml);
  1344.           if (need_aiff)
  1345.    printf("Multilingual file has been written with AIFF header informationn");
  1346. }
  1348. if (total_error_count_mc != 0)
  1349.     printf ("There were %d frames,  which were not in multichannel!!n", total_error_count_mc);
  1350. #ifdef Augmentation_7ch
  1351. if (total_error_count_aug != 0)
  1352.     printf ("There were %d frames,  which were not in augmented multichannel!!n", total_error_count_aug);
  1353. #endif
  1354.     if (print_out == 1)
  1355.     {
  1356. fclose (fp1);
  1357. fclose (fp2);
  1358. fclose (fp3);
  1359. fclose (fp4);
  1360. fclose (fp5);
  1361. fclose (fp6);
  1362.     }
  1363.     return (1);
  1364. }
  1367. void usage (void)  /* print syntax & exit */
  1368. {
  1369.    fprintf (stderr,
  1370.       "usage: %s                         queries for all arguments, orn",
  1371.    programName);
  1372.    fprintf (stderr,
  1373.       "       %s [-a] [-h] [-i] [-m] inputBS [outPCM]n", programName);
  1374.    fprintf (stderr,"wheren");
  1375.    fprintf (stderr," -a       write a RAW data sound file  (default: write an AIFF sound file)n");
  1376.    fprintf (stderr," -b       write alloc bits per subband to stdout (default: off)n");
  1377.    fprintf (stderr," -d       write dynx modes to stdout (default: off)n");
  1378.    fprintf (stderr," -h       bitstreams contain a 2304 byte header  (default: no header)n");
  1379.    fprintf (stderr," -i       write decoded bitstream-information in a filen"); 
  1380.       /* default: write encoded information */
  1381.    fprintf (stderr," -m       encoded bitstream is a MPEG 1 Layer-2 filen"); 
  1382.       /* default: decode a MPEG 2 Layer-2 bitstream */
  1383. #ifdef Augmentation_7ch
  1384.    fprintf (stderr," -n       encoded bitstream is a MPEG 2 Layer-2 file with 7.1 augmentationn"); 
  1385.       /* default: decode a MPEG 2 Layer-2 bitstream */
  1386. #endif
  1387.    fprintf (stderr," -s       write scf bits per subband to stdout (default: off)n");
  1388.    fprintf (stderr," -t       write tca modes to stdout (default: off)n");
  1389.    fprintf (stderr," inputBS  input bit stream of encoded audion");
  1390.    fprintf (stderr," outPCM   output PCM sound file (dflt inName)n");
  1391.    exit (1);
  1392. }