MUSICOUT.C
上传用户:dshsh2009
上传日期:2007-01-07
资源大小:155k
文件大小:24k
源码类别:

mpeg/mp3

开发平台:

Unix_Linux

  1. /**********************************************************************
  2. Copyright (c) 1991 MPEG/audio software simulation group, All Rights Reserved
  3. musicout.c
  4. **********************************************************************/
  5. /**********************************************************************
  6.  * MPEG/audio coding/decoding software, work in progress              *
  7.  *   NOT for public distribution until verified and approved by the   *
  8.  *   MPEG/audio committee.  For further information, please contact   *
  9.  *   Davis Pan, 708-538-5671, e-mail: pan@ukraine.corp.mot.com        *
  10.  *                                                                    *
  11.  * VERSION 4.3                                                        *
  12.  *   changes made since last update:                                  *
  13.  *   date   programmers                comment                        *
  14.  * 2/25/91  Douglas Wong        start of version 1.0 records          *
  15.  * 3/06/91  Douglas Wong        rename setup.h to dedef.h             *
  16.  *                              removed extraneous variables          *
  17.  *                              removed window_samples (now part of   *
  18.  *                              filter_samples)                       *
  19.  * 3/07/91  Davis Pan           changed output file to "codmusic"     *
  20.  * 5/10/91  Vish (PRISM)        Ported to Macintosh and Unix.         *
  21.  *                              Incorporated new "out_fifo()" which   *
  22.  *                              writes out last incomplete buffer.    *
  23.  *                              Incorporated all AIFF routines which  *
  24.  *                              are also compatible with SUN.         *
  25.  *                              Incorporated user interface for       *
  26.  *                              specifying sound file names.          *
  27.  *                              Also incorporated user interface for  *
  28.  *                              writing AIFF compatible sound files.  *
  29.  * 27jun91  dpwe (Aware)        Added musicout and &sample_frames as  *
  30.  *                              args to out_fifo (were glob refs).    *
  31.  *                              Used new 'frame_params' struct.       *
  32.  *                              Clean,simplify, track clipped output  *
  33.  *                              and total bits/frame received.        *
  34.  * 7/10/91  Earle Jennings      changed to floats to FLOAT            *
  35.  *10/ 1/91  S.I. Sudharsanan,   Ported to IBM AIX platform.           *
  36.  *          Don H. Lee,                                               *
  37.  *          Peter W. Farrett                                          *
  38.  *10/ 3/91  Don H. Lee          implemented CRC-16 error protection   *
  39.  *                              newly introduced functions are        *
  40.  *                              buffer_CRC and recover_CRC_error      *
  41.  *                              Additions and revisions are marked    *
  42.  *                              with "dhl" for clarity                *
  43.  * 2/11/92  W. Joseph Carter    Ported new code to Macintosh.  Most   *
  44.  *                              important fixes involved changing     *
  45.  *                              16-bit ints to long or unsigned in    *
  46.  *                              bit alloc routines for quant of 65535 *
  47.  *                              and passing proper function args.     *
  48.  *                              Removed "Other Joint Stereo" option   *
  49.  *                              and made bitrate be total channel     *
  50.  *                              bitrate, irrespective of the mode.    *
  51.  *                              Fixed many small bugs & reorganized.  *
  52.  *19 aug 92 Soren H. Nielsen    Changed MS-DOS file name extensions.  *
  53.  * 8/27/93 Seymour Shlien,      Fixes in Unix and MSDOS ports,        *
  54.  *         Daniel Lauzon, and                                         *
  55.  *         Bill Truerniet                                             *
  56.  *--------------------------------------------------------------------*
  57.  * 4/23/92  J. Pineda           Added code for layer III.  LayerIII   *
  58.  *          Amit Gulati         decoding is currently performed in    *
  59.  *                              two-passes for ease of sideinfo and   *
  60.  *                              maindata buffering and decoding.      *
  61.  *                              The second (computation) pass is      *
  62.  *                              activated with "decode -3 <outfile>"  *
  63.  * 10/25/92 Amit Gulati         Modified usage() for layerIII         *
  64.  * 12/10/92 Amit Gulati         Changed processing order of re-order- *
  65.  *                              -ing step.  Fixed adjustment of       *
  66.  *                              main_data_end pointer to exclude      *
  67.  *                              side information.                     *
  68.  *  9/07/93 Toshiyuki Ishino    Integrated Layer III with Ver 3.9.    *
  69.  *--------------------------------------------------------------------*
  70.  * 11/20/93 Masahiro Iwadare    Integrated Layer III with Ver 4.0.    *
  71.  *--------------------------------------------------------------------*
  72.  *  7/14/94 Juergen Koller      Bug fixes in Layer III code           *
  73.  *--------------------------------------------------------------------*
  74.  * 08/11/94 IIS                 Bug fixes in Layer III code           *
  75.  *--------------------------------------------------------------------*
  76.  * 11/04/94 Jon Rowlands        Prototype fixes                       *
  77.  **********************************************************************/
  78. #include        "common.h"
  79. #include        "decoder.h"
  80. /********************************************************************
  81. /*
  82. /*        This part contains the MPEG I decoder for Layers I & II.
  83. /*
  84. /*********************************************************************/
  85. /****************************************************************
  86. /*
  87. /*        For MS-DOS user (Turbo c) change all instance of malloc
  88. /*        to _farmalloc and free to _farfree. Compiler model hugh
  89. /*        Also make sure all the pointer specified are changed to far.
  90. /*
  91. /*****************************************************************/
  92. /*********************************************************************
  93. /*
  94. /* Core of the Layer II decoder.  Default layer is Layer II.
  95. /*
  96. /*********************************************************************/
  97. /* Global variable definitions for "musicout.c" */
  98. char *programName;
  99. int main_data_slots();
  100. int side_info_slots();
  101. /* Implementations */
  102. main(argc, argv)
  103. int argc;
  104. char **argv;
  105. {
  106. /*typedef short PCM[2][3][SBLIMIT];*/
  107. typedef short PCM[2][SSLIMIT][SBLIMIT];
  108.     PCM FAR *pcm_sample;
  109. typedef unsigned int SAM[2][3][SBLIMIT];
  110.     SAM FAR *sample;
  111. typedef double FRA[2][3][SBLIMIT];
  112.     FRA FAR *fraction;
  113. typedef double VE[2][HAN_SIZE];
  114.     VE FAR *w;
  115.     Bit_stream_struc  bs;
  116.     frame_params      fr_ps;
  117.     layer             info;
  118.     FILE              *musicout;
  119.     unsigned long     sample_frames;
  120.     int               i, j, k, x, stereo, done=FALSE, clip, sync;
  121.     int               error_protection, crc_error_count, total_error_count;
  122.     unsigned int      old_crc, new_crc;
  123.     unsigned int      bit_alloc[2][SBLIMIT], scfsi[2][SBLIMIT],
  124.                       scale_index[2][3][SBLIMIT];
  125.     unsigned long     bitsPerSlot, samplesPerFrame, frameNum = 0;
  126.     unsigned long     frameBits, gotBits = 0;
  127.     IFF_AIFF          pcm_aiff_data;
  128.     char              encoded_file_name[MAX_NAME_SIZE];
  129.     char              decoded_file_name[MAX_NAME_SIZE];
  130.     char       default_file_name[MAX_NAME_SIZE];
  131.     char              t[50];
  132.     int               need_aiff;
  133.     int               need_esps;        /* MI */
  134.     int topSb = 0;
  135. III_scalefac_t III_scalefac;
  136. III_side_info_t III_side_info;
  137. #ifdef  MACINTOSH
  138.     console_options.nrows = MAC_WINDOW_SIZE;
  139.     argc = ccommand(&argv);
  140. #endif
  141.     /* Most large variables are declared dynamically to ensure
  142.        compatibility with smaller machines */
  143.     pcm_sample = (PCM FAR *) mem_alloc((long) sizeof(PCM), "PCM Samp");
  144.     sample = (SAM FAR *) mem_alloc((long) sizeof(SAM), "Sample");
  145.     fraction = (FRA FAR *) mem_alloc((long) sizeof(FRA), "fraction");
  146.     w = (VE FAR *) mem_alloc((long) sizeof(VE), "w");
  147.     fr_ps.header = &info;
  148.     fr_ps.tab_num = -1;                /* no table loaded */
  149.     fr_ps.alloc = NULL;
  150.     for (i=0;i<HAN_SIZE;i++) for (j=0;j<2;j++) (*w)[j][i] = 0.0;
  151.     programName = argv[0];
  152.     if(argc==1) {        /* no command line args -> interact */
  153.        do {
  154.           printf ("Enter encoded file name <required>: ");
  155.           gets (encoded_file_name);
  156.           if (encoded_file_name[0] == NULL_CHAR)
  157.              printf ("Encoded file name is required. n");
  158.        } while (encoded_file_name[0] == NULL_CHAR);
  159.        printf (">>> Encoded file name is: %s n", encoded_file_name);
  160.     x=0;
  161.     while (x <= MAX_NAME_SIZE)
  162.     {
  163. default_file_name[x] = NULL_CHAR;
  164. ++x;
  165. }
  166.     x=0;
  167.     while (x <= 8)
  168.     {
  169.    default_file_name[x] = encoded_file_name[x];
  170.     if (encoded_file_name[++x] == '.') 
  171. x = 9;
  172.   
  173.     }
  174.    
  175.     strcat(default_file_name,DFLT_OPEXT);
  176.     printf("Enter MPEG decoded output file name <%s>: ",
  177.    default_file_name); /* 92-08-19 shn */
  178.     gets(decoded_file_name);
  179.     if (decoded_file_name[0] == NULL_CHAR) {
  180.  
  181. /* replace old extension with new one, 92-08-19 shn */
  182.        strcpy(decoded_file_name,default_file_name);
  183.  
  184.     }
  185.        printf (">>> MPEG decoded file name is: %s n", decoded_file_name);
  186.        printf(
  187.           "Do you wish to write an AIFF compatible sound file ? (y/<n>) : ");
  188.        gets(t);
  189.        if (*t == 'y' || *t == 'Y') need_aiff = TRUE;
  190.        else                        need_aiff = FALSE;
  191.        if (need_aiff)
  192.             printf(">>> An AIFF compatible sound file will be writtenn");
  193.        else printf(">>> A non-headered PCM sound file will be writtenn");
  194.        printf(
  195.           "Do you wish to exit (last chance before decoding) ? (y/<n>) : ");
  196.        gets(t);
  197.        if (*t == 'y' || *t == 'Y') exit(0);
  198.     }
  199.     else {        /* interpret CL Args */
  200.        int i=0, err=0;
  201.        need_aiff = FALSE;
  202.        need_esps = FALSE; /* MI */
  203.        encoded_file_name[0] = '';
  204.        decoded_file_name[0] = '';
  205.        while(++i<argc && err == 0) {
  206.           char c, *token, *arg, *nextArg;
  207.           int  argUsed;
  208.           token = argv[i];
  209.           if(*token++ == '-') {
  210.              if(i+1 < argc) nextArg = argv[i+1];
  211.              else           nextArg = "";
  212.              argUsed = 0;
  213.              while(c = *token++) {
  214.                 if(*token /* NumericQ(token) */) arg = token;
  215.                 else                             arg = nextArg;
  216.                 switch(c) {
  217.                    case 's':  topSb = atoi(arg); argUsed = 1;
  218.                       if(topSb<1 || topSb>SBLIMIT) {
  219.                          fprintf(stderr, "%s: -s band %s not %d..%dn",
  220.                                  programName, arg, 1, SBLIMIT);
  221.                          err = 1;
  222.                       }
  223.                       break;
  224.                    case 'A':  need_aiff = TRUE; break;
  225.                    case 'E':  need_esps = TRUE; break; /* MI */
  226.                    default:   fprintf(stderr,"%s: unrecognized option %cn",
  227.                                       programName, c);
  228.                       err = 1; break;
  229.                 }
  230.                 if(argUsed) {
  231.                    if(arg == token) token = ""; /* no more from token */
  232.                    else             ++i; /* skip arg we used */
  233.                    arg = ""; argUsed = 0;
  234.                 }
  235.              }
  236.           }
  237.           else {
  238.              if(encoded_file_name[0] == '')
  239.                 strcpy(encoded_file_name, argv[i]);
  240.              else
  241.                 if(decoded_file_name[0] == '')
  242.                    strcpy(decoded_file_name, argv[i]);
  243.                 else {
  244.                    fprintf(stderr,
  245.                            "%s: excess arg %sn", programName, argv[i]);
  246.                    err = 1;
  247.                 }
  248.           }
  249.        }
  250.        if(err || encoded_file_name[0] == '') usage();  /* never returns */
  251.        if(decoded_file_name[0] == '') {
  252.           strcpy(decoded_file_name, encoded_file_name);
  253.           strcat(decoded_file_name, DFLT_OPEXT);
  254.        }
  255.     }
  256.     /* report results of dialog / command line */
  257.     printf("Input file = '%s'  output file = '%s'n",
  258.            encoded_file_name, decoded_file_name);
  259.     if(need_aiff) printf("Output file written in AIFF formatn");
  260.     if(need_esps) printf("Output file written in ESPS formatn"); /* MI */
  261.     if ((musicout = fopen(decoded_file_name, "w+b")) == NULL) {
  262.        printf ("Could not create "%s".n", decoded_file_name);
  263.        exit(1);
  264.     }
  265.     open_bit_stream_r(&bs, encoded_file_name, BUFFER_SIZE);
  266.     if (need_aiff)
  267.        if (aiff_seek_to_sound_data(musicout) == -1) {
  268.           printf("Could not seek to PCM sound data in "%s".n",
  269.                  decoded_file_name);
  270.           exit(1);
  271.        }
  272.     sample_frames = 0;
  273.     while (!end_bs(&bs)) {
  274.        sync = seek_sync(&bs, SYNC_WORD, SYNC_WORD_LNGTH);
  275.        frameBits = sstell(&bs) - gotBits;
  276.        if(frameNum > 0)        /* don't want to print on 1st loop; no lay */
  277.           if(frameBits%bitsPerSlot)
  278.              fprintf(stderr,"Got %ld bits = %ld slots plus %ldn",
  279.                      frameBits, frameBits/bitsPerSlot, frameBits%bitsPerSlot);
  280.        gotBits += frameBits;
  281.        if (!sync) {
  282.           printf("Frame cannot be locatedn");
  283.           printf("Input stream may be emptyn");
  284.           done = TRUE;
  285.           /* finally write out the buffer */
  286.           if (info.lay != 1) out_fifo(*pcm_sample, 3, &fr_ps, done,
  287.                                       musicout, &sample_frames);
  288.           else               out_fifo(*pcm_sample, 1, &fr_ps, done,
  289.                                       musicout, &sample_frames);
  290.           break;
  291.        }
  292.        decode_info(&bs, &fr_ps);
  293.        hdr_to_frps(&fr_ps);
  294.        stereo = fr_ps.stereo;
  295.        error_protection = info.error_protection;
  296.        crc_error_count = 0;
  297.        total_error_count = 0;
  298.        if(frameNum == 0) WriteHdr(&fr_ps, stdout);  /* printout layer/mode */
  299. #ifdef ESPS
  300. if (frameNum == 0 && need_esps) {
  301. esps_write_header(musicout,(long) sample_frames, (double)
  302. s_freq[info.sampling_frequency] * 1000,
  303. (int) stereo, decoded_file_name );
  304. } /* MI */
  305. #endif
  306.        fprintf(stderr, "{%4lu}", frameNum++); fflush(stderr);
  307.        if (error_protection) buffer_CRC(&bs, &old_crc);
  308.        switch (info.lay) {
  309.           case 1: {
  310.              bitsPerSlot = 32;        samplesPerFrame = 384;
  311.              I_decode_bitalloc(&bs,bit_alloc,&fr_ps);
  312.              I_decode_scale(&bs, bit_alloc, scale_index, &fr_ps);
  313.              if (error_protection) {
  314.                 I_CRC_calc(&fr_ps, bit_alloc, &new_crc);
  315.                 if (new_crc != old_crc) {
  316.                    crc_error_count++;
  317.                    total_error_count++;
  318.                    recover_CRC_error(*pcm_sample, crc_error_count,
  319.                                      &fr_ps, musicout, &sample_frames);
  320.                    break;
  321.                 }
  322.                 else crc_error_count = 0;
  323.              }
  324.              clip = 0;
  325.              for (i=0;i<SCALE_BLOCK;i++) {
  326.                 I_buffer_sample(&bs,(*sample),bit_alloc,&fr_ps);
  327.                 I_dequantize_sample(*sample,*fraction,bit_alloc,&fr_ps);
  328.                 I_denormalize_sample((*fraction),scale_index,&fr_ps);
  329.                 if(topSb>0)        /* clear channels to 0 */
  330.                    for(j=topSb; j<fr_ps.sblimit; ++j)
  331.                       for(k=0; k<stereo; ++k)
  332.                          (*fraction)[k][0][j] = 0;
  333.                 for (j=0;j<stereo;j++) {
  334.                    clip += SubBandSynthesis (&((*fraction)[j][0][0]), j,
  335.                                              &((*pcm_sample)[j][0][0]));
  336.                 }
  337.                 out_fifo(*pcm_sample, 1, &fr_ps, done,
  338.                          musicout, &sample_frames);
  339.              }
  340.              if(clip > 0) printf("%d output samples clippedn", clip);
  341.              break;
  342.           }
  343.           case 2: {
  344.              bitsPerSlot = 8;        samplesPerFrame = 1152;
  345.              II_decode_bitalloc(&bs, bit_alloc, &fr_ps);
  346.              II_decode_scale(&bs, scfsi, bit_alloc, scale_index, &fr_ps);
  347.              if (error_protection) { 
  348.                 II_CRC_calc(&fr_ps, bit_alloc, scfsi, &new_crc);
  349.                 if (new_crc != old_crc) {
  350.                    crc_error_count++;
  351.                    total_error_count++;
  352.                    recover_CRC_error(*pcm_sample, crc_error_count,
  353.                                      &fr_ps, musicout, &sample_frames);
  354.                    break;
  355.                 }
  356.                 else crc_error_count = 0;
  357.              }
  358.              clip = 0;
  359.              for (i=0;i<SCALE_BLOCK;i++) {
  360.                 II_buffer_sample(&bs,(*sample),bit_alloc,&fr_ps);
  361.                 II_dequantize_sample((*sample),bit_alloc,(*fraction),&fr_ps);
  362.                 II_denormalize_sample((*fraction),scale_index,&fr_ps,i>>2);
  363.                 if(topSb>0)        /* debug : clear channels to 0 */
  364.                    for(j=topSb; j<fr_ps.sblimit; ++j)
  365.                       for(k=0; k<stereo; ++k)
  366.                          (*fraction)[k][0][j] =
  367.                          (*fraction)[k][1][j] =
  368.                          (*fraction)[k][2][j] = 0;
  369.                 for (j=0;j<3;j++) for (k=0;k<stereo;k++) {
  370.                    clip += SubBandSynthesis (&((*fraction)[k][j][0]), k,
  371.                                              &((*pcm_sample)[k][j][0]));
  372.                 }
  373.                 out_fifo(*pcm_sample, 3, &fr_ps, done, musicout,
  374.                          &sample_frames);
  375.              }
  376.              if(clip > 0) printf("%d samples clippedn", clip);
  377.              break;
  378.           }
  379.           case 3: {
  380.              int nSlots;
  381.              int gr, ch, ss, sb, main_data_end, flush_main ;
  382.      int  bytes_to_discard ;
  383.      static int frame_start = 0;
  384.              bitsPerSlot = 8;        samplesPerFrame = 1152;
  385.              III_get_side_info(&bs, &III_side_info, &fr_ps);
  386.              nSlots = main_data_slots(fr_ps);
  387.              for (; nSlots > 0; nSlots--)  /* read main data. */
  388.                 hputbuf((unsigned int) getbits(&bs,8), 8);
  389.      main_data_end = hsstell() / 8; /*of privious frame*/
  390.              if ( flush_main=(hsstell() % bitsPerSlot) ) { 
  391.                 hgetbits((int)(bitsPerSlot - flush_main));
  392. main_data_end ++;
  393.      }
  394.              bytes_to_discard = frame_start - main_data_end
  395.               - III_side_info.main_data_begin ;
  396.              if( main_data_end > 4096 )
  397.              {   frame_start -= 4096;
  398.                  rewindNbytes( 4096 );
  399.              }
  400.              frame_start += main_data_slots(fr_ps);
  401.              if (bytes_to_discard < 0) {
  402.          printf("Not enough main data to decode frame %d.  Frame discarded.n", 
  403.                         frameNum - 1); break;
  404.              }
  405.              for (; bytes_to_discard > 0; bytes_to_discard--) hgetbits(8);
  406.              clip = 0;
  407.              for (gr=0;gr<2;gr++) {
  408.                double lr[2][SBLIMIT][SSLIMIT],ro[2][SBLIMIT][SSLIMIT];
  409.                for (ch=0; ch<stereo; ch++) {
  410.                  long int is[SBLIMIT][SSLIMIT];   /* Quantized samples. */
  411.                  int part2_start;
  412.                  part2_start = hsstell();
  413.                  III_get_scale_factors(III_scalefac,&III_side_info,gr,ch,
  414. &fr_ps);
  415.                  III_hufman_decode(is, &III_side_info, ch, gr, part2_start,
  416.                                    &fr_ps);
  417.                  III_dequantize_sample(is, ro[ch], &III_scalefac,
  418.                                    &(III_side_info.ch[ch].gr[gr]), ch, &fr_ps);
  419.                }
  420.                III_stereo(ro,lr,III_scalefac,
  421.                             &(III_side_info.ch[0].gr[gr]), &fr_ps);
  422.                for (ch=0; ch<stereo; ch++) {
  423.                     double re[SBLIMIT][SSLIMIT];
  424.                     double hybridIn[SBLIMIT][SSLIMIT];/* Hybrid filter input */
  425.                     double hybridOut[SBLIMIT][SSLIMIT];/* Hybrid filter out */
  426.                     double polyPhaseIn[SBLIMIT];     /* PolyPhase Input. */
  427.                     III_reorder (lr[ch],re,&(III_side_info.ch[ch].gr[gr]),
  428.                                   &fr_ps);
  429.                     III_antialias(re, hybridIn, /* Antialias butterflies. */
  430.                                   &(III_side_info.ch[ch].gr[gr]), &fr_ps);
  431.                     for (sb=0; sb<SBLIMIT; sb++) { /* Hybrid synthesis. */
  432.                         III_hybrid(hybridIn[sb], hybridOut[sb], sb, ch,
  433.                                    &(III_side_info.ch[ch].gr[gr]), &fr_ps);
  434.                     }
  435.                     for (ss=0;ss<18;ss++) /*Frequency inversion for polyphase.*/
  436.                        for (sb=0; sb<SBLIMIT; sb++)
  437.                           if ((ss%2) && (sb%2))
  438.                              hybridOut[sb][ss] = -hybridOut[sb][ss];
  439.                     for (ss=0;ss<18;ss++) { /* Polyphase synthesis */
  440.                         for (sb=0; sb<SBLIMIT; sb++)
  441.                             polyPhaseIn[sb] = hybridOut[sb][ss];
  442.                         clip += SubBandSynthesis (polyPhaseIn, ch,
  443.                                                   &((*pcm_sample)[ch][ss][0]));
  444.                         }
  445.                     }
  446.                 /* Output PCM sample points for one granule. */
  447.                 out_fifo(*pcm_sample, 18, &fr_ps, done, musicout,
  448.                          &sample_frames);
  449.              }
  450.              if(clip > 0) printf("%d samples clipped.n", clip);
  451.              break;
  452.           }
  453.        }
  454.     }
  455.     if (need_aiff) {
  456.        pcm_aiff_data.numChannels       = stereo;
  457.        pcm_aiff_data.numSampleFrames   = sample_frames;
  458.        pcm_aiff_data.sampleSize        = 16;
  459.        pcm_aiff_data.sampleRate        = s_freq[info.sampling_frequency]*1000;
  460. #ifdef IFF_LONG
  461.        pcm_aiff_data.sampleType        = IFF_ID_SSND;
  462. #else
  463.        strncpy(&pcm_aiff_data.sampleType,IFF_ID_SSND,4);
  464. #endif
  465.        pcm_aiff_data.blkAlgn.offset    = 0;
  466.        pcm_aiff_data.blkAlgn.blockSize = 0;
  467.        if (aiff_write_headers(musicout, &pcm_aiff_data) == -1) {
  468.           printf("Could not write AIFF headers to "%s"n",
  469.                  decoded_file_name);
  470.           exit(2);
  471.        }
  472.     }
  473.     printf("Avg slots/frame = %.3f; b/smp = %.2f; br = %.3f kbpsn",
  474.            (FLOAT) gotBits / (frameNum * bitsPerSlot),
  475.            (FLOAT) gotBits / (frameNum * samplesPerFrame),
  476.            (FLOAT) gotBits / (frameNum * samplesPerFrame) *
  477.            s_freq[info.sampling_frequency]);
  478.     close_bit_stream_r(&bs);
  479.     fclose(musicout);
  480.     /* for the correct AIFF header information */
  481.     /*             on the Macintosh            */
  482.     /* the file type and the file creator for  */
  483.     /* Macintosh compatible Digidesign is set  */
  484. #ifdef  MACINTOSH
  485.     if (need_aiff) set_mac_file_attr(decoded_file_name, VOL_REF_NUM,
  486.                                      CREATR_DEC_AIFF, FILTYP_DEC_AIFF);
  487.     else           set_mac_file_attr(decoded_file_name, VOL_REF_NUM,
  488.                                      CREATR_DEC_BNRY, FILTYP_DEC_BNRY);
  489. #endif
  490.     printf("Decoding of "%s" is finishedn", encoded_file_name);
  491.     printf("The decoded PCM output file name is "%s"n", decoded_file_name);
  492.     if (need_aiff)
  493.        printf(""%s" has been written with AIFF header informationn",
  494.               decoded_file_name);
  495.     exit( 0 );
  496. }
  497. static void usage()  /* print syntax & exit */
  498. {
  499.    fprintf(stderr,
  500.       "usage: %s                         queries for all arguments, orn",
  501.        programName);
  502.    fprintf(stderr,
  503.       "       %s [-A][-s sb] inputBS [outPCM]n", programName);
  504.    fprintf(stderr,"wheren");
  505.    fprintf(stderr," -A       write an AIFF output PCM sound filen");
  506.    fprintf(stderr," -s sb    resynth only up to this sb (debugging only)n");
  507.    fprintf(stderr," inputBS  input bit stream of encoded audion");
  508.    fprintf(stderr," outPCM   output PCM sound file (dflt inName+%s)n",
  509.            DFLT_OPEXT);
  510.    exit(1);
  511. }