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

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

mpeg/mp3

开发平台：

C/C++

COMMON.C：源码内容
  1. /**********************************************************************
  2.  * ISO MPEG Audio Subgroup Software Simulation Group (1996)
  3.  * ISO 13818-3 MPEG-2 Audio Multichannel Decoder
  4.  *
  5.  * $Id: common.c 1.12 1997/01/17 04:25:12 rowlands Exp $
  6.  *
  7.  * $Log: common.c $
  8.  * Revision 1.12  1997/01/17 04:25:12  rowlands
  9.  * Fixed base bitstream joint stereo bugs.
  10.  *
  11.  * Removed incorrect dependence on jsbound in the calculation of the
  12.  * mc crc.
  13.  *
  14.  * Revision 1.11  1996/06/18 22:47:29  rowlands
  15.  * Added extra debug output
  16.  *
  17.  * Revision 1.10  1996/04/18 05:37:23  rowlands
  18.  * Release following Florence meeting
  19.  *
  20.  * Revision 1.9  1996/02/12 07:13:06  rowlands
  21.  * Release following Munich meeting
  22.  *
  23.  * Revision 1.8.1.1  1996/01/20  17:28:45  rowlands
  24.  * Received from Ralf Schwalbe (Telekom FTZ) - includes prediction
  25.  *
  26.  * Revision 1.6.1.3  1995/08/14  08:10:19  tenkate
  27.  * ML-LSF added Warner ten Kate 7/8/95 (Philips)
  28.  * "alloc_4" table copied from LSF-directory.
  29.  * SmpFrqIndex() extended.
  30.  *
  31.  * Revision 1.6.1.1  1995/07/14  06:12:46  rowlands
  32.  * Updated dynamic crosstalk from FTZ: revision FTZ_03
  33.  *
  34.  * Revision 1.6  1995/06/22  01:23:00  rowlands
  35.  * Merged FTZ dynamic crosstalk
  36.  *
  37.  * Revision 1.3.3.1  1995/06/16  08:40:12  rowlands
  38.  * Input from Sang Wook Kim (Samsung AIT)
  39.  *
  40.  **********************************************************************/
  42. /**********************************************************************
  43.  *   date   programmers                comment                        *
  44.  * 2/25/91  Douglas Wong        start of version 1.0 records          *
  45.  * 3/06/91  Douglas Wong        rename setup.h to dedef.h             *
  46.  *                              removed extraneous variables          *
  47.  *                              removed window_samples (now part of   *
  48.  *                              filter_samples)                       *
  49.  * 3/07/91  Davis Pan           changed output file to "codmusic"     *
  50.  * 5/10/91  Vish (PRISM)        Ported to Macintosh and Unix.         *
  51.  *                              Incorporated new "out_fifo()" which   *
  52.  *                              writes out last incomplete buffer.    *
  53.  *                              Incorporated all AIFF routines which  *
  54.  *                              are also compatible with SUN.         *
  55.  *                              Incorporated user interface for       *
  56.  *                              specifying sound file names.          *
  57.  *                              Also incorporated user interface for  *
  58.  *                              writing AIFF compatible sound files.  *
  59.  * 27jun91  dpwe (Aware)        Added musicout and &sample_frames as  *
  60.  *                              args to out_fifo (were glob refs).    *
  61.  *                              Used new 'frame_params' struct.       *
  62.  *                              Clean,simplify, track clipped output  *
  63.  *                              and total bits/frame received.        *
  64.  * 7/10/91  Earle Jennings      changed to floats to FLOAT            *
  65.  *10/ 1/91  S.I. Sudharsanan,   Ported to IBM AIX platform.           *
  66.  *          Don H. Lee,                                               *
  67.  *          Peter W. Farrett                                          *
  68.  *10/ 3/91  Don H. Lee          implemented CRC-16 error protection   *
  69.  *                              newly introduced functions are        *
  70.  *                              buffer_CRC and recover_CRC_error      *
  71.  *                              Additions and revisions are marked    *
  72.  *                              with "dhl" for clarity                *
  73.  * 2/11/92  W. Joseph Carter    Ported new code to Macintosh.  Most   *
  74.  *                              important fixes involved changing     *
  75.  *                              16-bit ints to long or unsigned in    *
  76.  *                              bit alloc routines for quant of 65535 *
  77.  *                              and passing proper function args.     *
  78.  *                              Removed "Other Joint Stereo" option   *
  79.  *                              and made bitrate be total channel     *
  80.  *                              bitrate, irrespective of the mode.    *
  81.  *                              Fixed many small bugs & reorganized.  *
  82.  **********************************************************************
  83.  *                                                                    *
  84.  *                                                                    *
  85.  *  MPEG/audio Phase 2 coding/decoding multichannel                   *
  86.  *                                                                    *
  87.  *  Version 1.0                                                       *
  88.  *                                                                    *
  89.  *  7/27/93        Susanne Ritscher,  IRT Munich                      *
  90.  *                                                                    *
  91.  *                  thanks to                                         *
  92.  *                  Ralf Schwalbe,    Telekom FTZ Berlin              *
  93.  *                  Heiko Purnhagen,  Uni Hannover                    *
  94.  *                                                                    *
  95.  *  Version 2.0                                                       *
  96.  *                                                                    *
  97.  *  8/27/93        Susanne Ritscher, IRT Munich                       *
  98.  *                 Channel-Switching is working                       *
  99.  *                                                                    *
  100.  *  Version 2.1                                                       *
  101.  *                                                                    *
  102.  *  9/1/93         Susanne Ritscher,  IRT Munich                      *
  103.  *                 all channels normalized                            *
  104.  *                                                                    *
  105.  *  Version 3.0                                                       *
  106.  *                                                                    *
  107.  *  06/16/94       Ralf Schwalbe, Telekom FTZ Berlin                  *
  108.  *                 all sources and variables adapted due to MPEG-2 -  *
  109.  *                 DIS from March 1994                                *
  110.  *                  - dematrix and denormalize procedure              *
  111.  *                  - new tc - allocation (0-7)                       *
  112.  *                  - some new structures and variables as a basis    *
  113.  *                    for further decoding modes                      *
  114.  **********************************************************************
  115.  *       *
  116.  *  Version 1.0                                                       *
  117.  *                                                                    *
  118.  *  11/04/94       Ralf Schwalbe,  Telekom FTZ Berlin                 *
  119.  *                  - decoding tc-allocation                          *
  120.  *                  - some new subroutines, globale variables and     *
  121.  *                    structures (important to handle the ext. bitst.)*
  122.  *                  - changed all functions to ANSI-C funktion header *
  123.  *     - corrected some bugs to decode bitstreams > 512kB*
  124.  *       *
  125.  *  Version 1.1                                                       *
  126.  *                                                                    *
  127.  *  12/07/94       Ralf Schwalbe,  Telekom FTZ Berlin                 *
  128.  *                  - decoding extension bitstream                    *
  129.  *                                                                    *
  130.  *  Version 1.1.1                                                     *
  131.  *                 Ralf Schwalbe, Telekom FTZ Berlin                  *
  132.  *                  - fixed some bugs                                 *
  133.  *                                                                    *
  134.  *  Version 1.2                                                       *
  135.  *                                                                    *
  136.  *  6/21/95       Ralf Schwalbe, Deutsche Telekom FTZ Berlin          *
  137.  *                 - decoding dynamic crosstalk                      *
  138.  *                  - decoding phantom center                         * 
  139.  *                  - decoding MPEG1 compatible part only (stereo)    *
  140.  *                  - corrected some settings and bugs                *
  141.  *                                                                    *  
  142.  *  7/12/95   Ralf Schwalbe, Deutsche Telekom FTZ Berlin          *
  143.  *     - corrected dynamic crosstalk        *
  144.  *     - 3/2,3/1,3/0,2/1, channel configurations         *
  145.  *                    are working                                     *
  146.  *          *
  147.  * 10/31/95   Ralf Schwalbe, Deutsche Telekom FTZ Berlin       *
  148.  *     - decoding of LFE-channel is working         *
  149.  *     - corrected any settings and bugs       *
  150.  *     - corrected table-switch for channel mode < 3/2   *
  151.  *       and tc-allocation / dyn-crosstalk       *
  152.  *     - still a problem with compl. bitstream 18/19     *
  153.  *       *
  154.  * 01/12/96   Ralf Schwalbe, Deutsche Telekom TZ Darmstadt       *
  155.  *     - decoder prediction installed       *
  156.  *     - problem with compl. bitstream 19 solved,        *
  157.  *       *  
  158.  *  Version 2.0                                                       *
  159.  *                                                                    *
  160.  * 01/28/97       Frans de Bont, Philips Sound & Vision, Eindhoven    *
  161.  *     - simultaneous use of tc allocation and dynamic   *
  162.  *       crosstalk working for all configurations       *
  163.  *     - prediction working for all configurations       *
  164.  *     - layer 1 MC working       *
  165.  *     - variable bit rate and extension bitstreams      *
  166.  *     - fully compliant to 13818-3.2                    *
  167.  *                                                                    *       *     
  168.  **********************************************************************/
  170. /**********************************************************************
  171.  *                                                                    *
  172.  *  06/06/95       Sang Wook Kim,  Samsung AIT                        *
  173.  *                 Bug corrected in pick_table                        *
  174.  *       *
  175.  **********************************************************************/
  177. /***********************************************************************
  178. *
  179. *  Global Include Files
  180. *
  181. ***********************************************************************/
  183. #include        "common.h"
  184. #include        "decoder.h" /* 23/03/1995 JMZ Test_bit_stream */
  186. #ifdef  MACINTOSH
  188. #include        <SANE.h>
  189. #include        <pascal.h>
  191. #endif
  193. #include <string.h>
  194. #include <ctype.h>
  196. /***********************************************************************
  197. *
  198. *  Global Variable Definitions
  199. *
  200. ***********************************************************************/
  201. char *mode_names[4] = { "stereo", "j-stereo", "dual-ch", "single-ch" };
  202. char *layer_names[3] = { "I", "II", "III" };
  204. double  s_freq[4] = {44.1, 48, 32, 0};
  206. int   bitrate[3][15] = {
  207.  {1000,32,64,96,128,160,192,224,256,288,320,352,384,416,448},
  208.  {1000,32,48,56,64,80,96,112,128,160,192,224,256,320,384},
  209.  {1000,32,40,48,56,64,80,96,112,128,160,192,224,256,320}
  210.         };
  212. double multiple[64] = {
  213. 2.00000000000000, 1.58740105196820, 1.25992104989487,
  214. 1.00000000000000, 0.79370052598410, 0.62996052494744, 0.50000000000000,
  215. 0.39685026299205, 0.31498026247372, 0.25000000000000, 0.19842513149602,
  216. 0.15749013123686, 0.12500000000000, 0.09921256574801, 0.07874506561843,
  217. 0.06250000000000, 0.04960628287401, 0.03937253280921, 0.03125000000000,
  218. 0.02480314143700, 0.01968626640461, 0.01562500000000, 0.01240157071850,
  219. 0.00984313320230, 0.00781250000000, 0.00620078535925, 0.00492156660115,
  220. 0.00390625000000, 0.00310039267963, 0.00246078330058, 0.00195312500000,
  221. 0.00155019633981, 0.00123039165029, 0.00097656250000, 0.00077509816991,
  222. 0.00061519582514, 0.00048828125000, 0.00038754908495, 0.00030759791257,
  223. 0.00024414062500, 0.00019377454248, 0.00015379895629, 0.00012207031250,
  224. 0.00009688727124, 0.00007689947814, 0.00006103515625, 0.00004844363562,
  225. 0.00003844973907, 0.00003051757813, 0.00002422181781, 0.00001922486954,
  226. 0.00001525878906, 0.00001211090890, 0.00000961243477, 0.00000762939453,
  227. 0.00000605545445, 0.00000480621738, 0.00000381469727, 0.00000302772723,
  228. 0.00000240310869, 0.00000190734863, 0.00000151386361, 0.00000120155435,
  229. 1E-20
  230. };
  232. int sb_groups[12] = { 0, 1, 2, 3, 4, 5, 6, 7, 9, 11, 15, 26};
  234. /*int transmission_channel[8][5] = {{0, 1, 2, 3, 4},
  235.     {0, 1, 5, 3, 4},
  236.     {0, 1, 6, 3, 4},
  237.     {0, 1, 2, 5, 4},
  238.     {0, 1, 2, 3, 6},
  239.     {0, 1, 2, 5, 6},
  240.     {0, 1, 6, 5, 4},
  241.     {0, 1, 5, 3, 6}};
  242. */
  245. /***********************************************************************
  246. *
  247. *  Global Function Definitions
  248. *
  249. ***********************************************************************/
  251. /* The system uses a variety of data files.  By opening them via this
  252.    function, we can accommodate various locations. */
  254. FILE *OpenTableFile(char *name)
  255. {
  256. char fulname[80];
  257. char *envdir;
  258. FILE *f;
  260.  fulname[0] = '';
  262. #ifdef TABLES_PATH
  263.    strcpy(fulname, TABLES_PATH);   /* default relative path for tables */
  264. #endif /* TABLES_PATH */          /* (includes terminal path seperator */
  266. #ifdef UNIX                       /* envir. variables for UNIX only */
  267.    {
  268. char *getenv(const char *);
  270.         envdir = getenv(MPEGTABENV);   /* check for environment */
  271.         if(envdir != NULL)
  272.             strcpy(fulname, envdir);
  273.         strcat(fulname, PATH_SEPARATOR);  /* add a "/" on the end */
  274.       }
  275. #endif /* UNIX */
  277. strcat(fulname, name);
  278.     if( (f=fopen(fulname,"r"))==NULL ) {
  279.         fprintf(stderr,"OpenTable: could not find %sn", fulname);
  281. #ifdef UNIX
  282.           if(envdir != NULL)
  283. fprintf(stderr,"Check %s directory '%s'n",MPEGTABENV, envdir);
  284.           else
  285.             fprintf(stderr,"Check local directory './%s' or setenv %sn",
  286.                     TABLES_PATH, MPEGTABENV);
  287. #else /* not unix : no environment variables */
  289. #ifdef TABLES_PATH
  290.             fprintf(stderr,"Check local directory './%s'n",TABLES_PATH);
  291. #endif /* TABLES_PATH */
  293. #endif /* UNIX */
  295.     }
  296.     return f;
  297. }
  299. /***********************************************************************
  300. /*
  301. /* Read one of the data files ("alloc_*") specifying the bit allocation/
  302. /* quatization parameters for each subband in layer II encoding
  303. /*
  304. /**********************************************************************/
  305. /* read in table, return # subbands */
  306. int read_bit_alloc(int table, al_table *alloc)  
  308. {
  309.         long a, b, c, d, i, j;
  310. FILE *fp;
  311.         char name[16], t[80];
  312.         int sblim;
  314.         strcpy(name, "alloc_0");
  316. switch (table) {
  317.                 case 0 : name[6] = '0';         break;
  318.                 case 1 : name[6] = '1';         break;
  319.                 case 2 : name[6] = '2';         break;
  320.                 case 3 : name[6] = '3';         break;
  321.                 case 4 : name[6] = '4';         break; /* for MultiLingual-LSF , WtK 7/8/95 */
  322.                 default : name[6] = '0';
  323. }
  325.         if (!(fp = OpenTableFile(name))) {
  326.                 fprintf (stderr, "Please check bit allocation table %sn", name);
  327.                 exit (1);
  328.         }
  330. #ifdef PrintBitDebug
  331. printf("using bit allocation table %sn", name);
  332. #endif
  334.         fgets(t, 80, fp);
  335.         sscanf(t, "%dn", &sblim);
  336.         while (!feof(fp)) {
  337. fgets(t, 80, fp);
  338.                 sscanf(t, "%ld %ld %ld %ld %ld %ldn", &i, &j, &a, &b, &c, &d);
  339. (*alloc)[i][j].steps = a;
  340.                         (*alloc)[i][j].bits  = b;
  341.                         (*alloc)[i][j].group = c;
  342.                         (*alloc)[i][j].quant = d;
  343.         }
  344.         fclose(fp);
  345. return sblim;
  346. }
  348. /***********************************************************************
  349. /*
  350. /* Using the decoded info the appropriate possible quantization per
  351. /* subband table is loaded
  352. /*
  353. /**********************************************************************/
  354. /* choose table, load if necess, return # sb's */
  355. int pick_table (frame_params *fr_ps)   
  356. {
  357.    int table, tableml, lay, ws, bsp, br_per_ch, sfrq;
  358.    int sblim = fr_ps->sblimit;     /* return current value if no load */
  360.    lay = fr_ps->header->lay - 1;
  361.    bsp = fr_ps->header->bitrate_index;
  363.    /* BUG corrected SWKim060695 */
  364.    br_per_ch = bitrate[lay][bsp] / (fr_ps->stereo /* + fr_ps->mc_channel*/);
  365.    ws = fr_ps->header->sampling_frequency;
  366.    sfrq = s_freq[ws];
  368.    /* decision rules refer to per-channel bitrates (kbits/sec/chan) */
  369.    if ((sfrq == 48 && br_per_ch >= 56) ||
  370.        (br_per_ch >= 56 && br_per_ch <= 80)) table = 0;
  371.    else if (sfrq != 48 && br_per_ch >= 96) table = 1;
  372.    else if (sfrq != 32 && br_per_ch <= 48) table = 2;
  373.    else table = 3;
  374.    if (fr_ps->tab_num != table)
  375.    {
  376.       if (fr_ps->tab_num >= 0)
  377.         mem_free ((void **) &(fr_ps->alloc));
  378.       fr_ps->alloc = (al_table *) mem_alloc (sizeof(al_table), "alloc");
  379.       sblim = read_bit_alloc (fr_ps->tab_num = table, fr_ps->alloc);
  380.    }
  381.    return sblim;
  382. }
  385. /* choose MC and ML tables, load if necess, return # sb's */
  386. void mc_pick_table (frame_params *fr_ps)   
  387. {
  388.    int tablemc, tableml, ws, sfrq;
  390.    /* BUG corrected SWKim060695 */
  391.    ws = fr_ps->header->sampling_frequency;
  392.    sfrq = s_freq[ws];
  394.    /* decision rules refer to sampling frequency */
  395.    if (sfrq == 48)
  396.       tablemc = 0;
  397.    else
  398.       tablemc = 1;
  400.    if (fr_ps->tab_num_mc != tablemc)
  401.    {
  402.       if (fr_ps->tab_num >= 0)
  403.  mem_free ((void **)& (fr_ps->alloc_mc));
  404.       fr_ps->alloc_mc = (al_table *) mem_alloc (sizeof(al_table), "alloc_mc");
  405.       fr_ps->sblimit_mc = read_bit_alloc (fr_ps->tab_num_mc = tablemc, fr_ps->alloc_mc);
  406.    }
  408.    if (fr_ps->header->no_of_multi_lingual_ch > 0)
  409.    {
  410.       /* LSF MultiLingual 7/8/95 WtK */
  411.       if (fr_ps->header->multi_lingual_fs == 0)
  412.       {
  413.  fr_ps->tab_num_ml = fr_ps->tab_num_mc;
  414.  fr_ps->alloc_ml   = fr_ps->alloc_mc;
  415.  fr_ps->sblimit_ml = fr_ps->sblimit_mc;
  416.       }
  417.       else
  418.       {
  419.  tableml = 4;
  420.  if (fr_ps->tab_num_ml != tableml)
  421.  {
  422.     if (fr_ps->tab_num_ml >= 0)
  423.        mem_free ((void **)&(fr_ps->alloc_ml));
  424.     fr_ps->alloc_ml = (al_table *) mem_alloc (sizeof(al_table), "alloc_ml");
  425.     fr_ps->sblimit_ml = read_bit_alloc (fr_ps->tab_num_ml = tableml, fr_ps->alloc_ml);
  426.  }
  427.       } /* ml_fs==0 */
  428.    } /* ml_ch>0 */
  429. }
  432. int js_bound(int lay, int m_ext)
  433. {
  434. static int jsb_table[3][4] =  { { 4, 8, 12, 16 }, { 4, 8, 12, 16},
  435. { 0, 4, 8, 16} };  /* lay+m_e -> jsbound */
  437.     if(lay<1 || lay >3 || m_ext<0 || m_ext>3) {
  438.         fprintf (stderr, "js_bound bad layer/modext (%d/%d)n", lay, m_ext);
  439.         exit (1);
  440.     }
  441.     return(jsb_table[lay-1][m_ext]);
  442. }
  444. /* interpret data in hdr str to fields in fr_ps */
  445. void hdr_to_frps (frame_params *fr_ps)
  446. {
  447. layer *hdr = fr_ps->header;     /* (or pass in as arg?) */
  449. fr_ps->actual_mode = hdr->mode;
  450. if (hdr->mode != MPG_MD_NONE)
  451. fr_ps->stereo = (hdr->mode == MPG_MD_MONO) ? 1 : 2;
  452. else
  453. fr_ps->stereo = 0;
  455. Bitrate = bitrate[hdr->lay-1][hdr->bitrate_index];
  456. S_freq = s_freq[hdr->sampling_frequency];
  458. if (hdr->lay == 2)
  459. fr_ps->sblimit = pick_table (fr_ps);
  460. else
  461. fr_ps->sblimit = SBLIMIT;
  463. if (hdr->mode == MPG_MD_JOINT_STEREO)
  464. fr_ps->jsbound = js_bound (hdr->lay, hdr->mode_ext);
  465. else
  466. fr_ps->jsbound = fr_ps->sblimit;
  467. }
  469. /* interpret data in hdr str to fields in fr_ps */
  470. void mc_hdr_to_frps (frame_params *fr_ps)
  471. {
  472. layer *hdr = fr_ps->header;     /* (or pass in as arg?) */
  474. if (fr_ps->header->surround == 3) /* second stereo program */
  475. if (fr_ps->header->center == 1 || fr_ps->header->center == 3)
  476. {
  477. fr_ps->mc_channel = 3;     /* 3 + 2 */
  478. fr_ps->alloc_bits = 2;     /* 2 bits for tc_alloc */
  479. fr_ps->dyn_cross_bits = 1;
  480. fr_ps->pred_mode = 2; /* FdB for prediction  3/0 */
  481. }
  482. else
  483. {
  484. fr_ps->mc_channel = 2;     /* 2 + 2 or 1 + 2 */
  485. fr_ps->alloc_bits = 0;     /* 0 bits for tc_alloc */
  486. fr_ps->dyn_cross_bits = 0;
  487. fr_ps->pred_mode = 5; /* FdB for prediction  2/0 */
  488. }
  489. else if (fr_ps->header->surround == 2) /* stereo surround */
  490. if (fr_ps->header->center == 1 || fr_ps->header->center == 3)
  491. {
  492. fr_ps->mc_channel = 3;  /* 3 mc-channels 3/2 */
  493. fr_ps->alloc_bits = 3;  /* 3 bits for tc_alloc */
  494. fr_ps->dyn_cross_bits = 4; 
  495. fr_ps->pred_mode = 0; /* 7/4/95 R.S. for prediction  3/2 */
  496. }
  497. else
  498. {
  499. fr_ps->mc_channel = 2;  /* 2 mc-channels 2/2 */
  500. fr_ps->alloc_bits = 2;  /* 2 bits for tc_alloc */
  501. fr_ps->dyn_cross_bits = 3;
  502. fr_ps->pred_mode = 3; /* 7/4/95 R.S. for prediction  2/2 */
  503. }
  504. else if (fr_ps->header->surround == 1) /* mono surround */
  505. if (fr_ps->header->center == 1 || fr_ps->header->center == 3)
  506. {
  507. fr_ps->mc_channel = 2;  /* 3/1 */
  508. fr_ps->alloc_bits = 3;  /* 3 bits for tc_alloc */
  509. fr_ps->dyn_cross_bits = 3;
  510. fr_ps->pred_mode = 1; /* 7/4/95 R.S. for prediction 3/1 */
  511. }
  512. else
  513. {
  514. fr_ps->mc_channel = 1; /* 2/1 */
  515. fr_ps->alloc_bits = 2;  /* 2 bits for tc_alloc */
  516. fr_ps->dyn_cross_bits = 1;
  517. fr_ps->pred_mode = 4; /* 7/4/95 R.S. for prediction 2/1 */
  518. }
  519. else /* no surround */
  520. if (fr_ps->header->center == 1 || fr_ps->header->center == 3)
  521. {
  522. fr_ps->mc_channel = 1; /* 3/0 */
  523. fr_ps->alloc_bits = 2; /* 2 bits for tc_alloc */ 
  524. fr_ps->dyn_cross_bits = 1;
  525. fr_ps->pred_mode = 2;         /* 3/0 */
  526. }
  527. else
  528. {
  529. fr_ps->mc_channel = 0; /* 2/0 or 1/0 */
  530. fr_ps->alloc_bits = 0;  /* 2 bits for tc_alloc */
  531. fr_ps->dyn_cross_bits = 0;
  532. fr_ps->pred_mode = 5; /* FdB for prediction  2/0 */
  533. }
  535. if ( (hdr->no_of_multi_lingual_ch>0) && (hdr->multi_lingual_layer>0) )
  536. {
  537. fprintf (stderr, "MultiLingual not in Layer 2!n exit.n");
  538. exit (1);
  539. }
  541. if (fr_ps->mc_channel > 0 || hdr->lfe || hdr->no_of_multi_lingual_ch > 0)
  542. mc_pick_table (fr_ps);
  543. }
  545. void WriteHdr(frame_params *fr_ps, FILE *s)
  546. {
  547. layer *info = fr_ps->header;
  549.    fprintf(s, "HDR:  s=FFF, id=%X, l=%X, ep=%X, br=%X, sf=%X, pd=%X, ",
  550.    info->version, info->lay, !info->error_protection,
  551.    info->bitrate_index, info->sampling_frequency, info->padding);
  552.    fprintf(s, "pr=%X, m=%X, js=%X, c=%X, o=%X, e=%Xn",
  553.    info->extension, info->mode, info->mode_ext,
  554.    info->copyright, info->original, info->emphasis);
  555.    fprintf(s, "layer=%s, MPEG1-bitrate=%d, sfrq=%.1f, mode=%s, ",
  556.    layer_names[info->lay-1],Bitrate=bitrate[info->lay-1][info->bitrate_index],
  557.    S_freq=s_freq[info->sampling_frequency], mode_names[info->mode]);
  558.    fprintf(s, "sblim=%d, jsbd=%d, MPEG1-ch=%dn",
  559.    fr_ps->sblimit, fr_ps->jsbound, fr_ps->stereo);
  560.    fflush(s);
  561. }
  563. /*******************************************************************************
  564. *
  565. *  Allocate number of bytes of memory equal to "block".
  566. *
  567. *******************************************************************************/
  569. void    *mem_alloc(unsigned long block, char *item)
  570. {
  571. void   *ptr;
  573. #ifdef  MACINTOSH
  574.     ptr = NewPtr(block);
  575. #endif
  577. #ifdef  MSDOS
  578.  ptr = (void *)_fmalloc((unsigned int)block); /* far memory, 92-07-08 sr */
  579. #endif
  581. #if ! defined (MACINTOSH) && ! defined (MSDOS)
  582.     ptr = (void *) malloc(block);
  583. #endif
  584.  if (ptr != NULL){
  585. #ifdef  MSDOS
  586. _fmemset(ptr, 0, (unsigned int)block); /* far memory, 92-07-08 sr */
  587. #else
  588. memset(ptr, 0, block);
  589. #endif
  590.     }
  591.     else{
  592.         fprintf (stderr, "Unable to allocate %sn", item);
  593.         exit (1);
  594.     }
  596.     return(ptr);
  598. }
  600. /****************************************************************************
  601. *
  602. *  Free memory pointed to by "*ptr_addr".
  603. *
  604. *****************************************************************************/
  606. void    mem_free(void **ptr_addr)
  607. {
  609.     if (*ptr_addr != NULL){
  610. #ifdef  MACINTOSH
  611.         DisposPtr(*ptr_addr);
  612. #else
  613.         free(*ptr_addr);
  614. #endif
  615.         *ptr_addr = NULL;
  616.     }
  618. }
  620. /****************************************************************************
  621. *
  622. *  Routines to convert between the Apple SANE extended floating point format
  623. *  and the IEEE double precision floating point format.  These routines are
  624. *  called from within the Audio Interchange File Format (AIFF) routines.
  625. *
  626. *****************************************************************************/
  628. /*
  629. *** Apple's 80-bit SANE extended has the following format:
  631.  1       15      1            63
  632. +-+-------------+-+-----------------------------+
  633. |s|       e     |i|            f                |
  634. +-+-------------+-+-----------------------------+
  635.   msb        lsb   msb                       lsb
  637. The value v of the number is determined by these fields as follows:
  638. If 0 <= e < 32767,              then v = (-1)^s * 2^(e-16383) * (i.f).
  639. If e == 32767 and f == 0,       then v = (-1)^s * (infinity), regardless of i.
  640. If e == 32767 and f != 0,       then v is a NaN, regardless of i.
  642. *** IEEE Draft Standard 754 Double Precision has the following format:
  644. MSB
  645. +-+---------+-----------------------------+
  646. |1| 11 Bits |           52 Bits           |
  647. +-+---------+-----------------------------+
  648.  ^     ^                ^
  649.  |     |                |
  650.  Sign  Exponent         Mantissa
  651. */
  653. /*****************************************************************************
  654. *
  655. *  double_to_extended()
  656. *
  657. *  Purpose:     Convert from IEEE double precision format to SANE extended
  658. *               format.
  659. *
  660. *  Passed:      Pointer to the double precision number and a pointer to what
  661. *               will hold the Apple SANE extended format value.
  662. *
  663. *  Outputs:     The SANE extended format pointer will be filled with the
  664. *               converted value.
  665. *
  666. *  Returned:    Nothing.
  667. *
  668. *****************************************************************************/
  671. void    double_to_extended(double *pd, char ps[10])
  672. {
  674. #ifdef  MACINTOSH
  676.         x96tox80(pd, (extended *) ps);
  678. #else
  680. /* fixed bus alignment error, HP 27-may-93 */
  682. register unsigned long  top2bits;
  684. register unsigned short *ps2;
  685. register IEEE_DBL       *p_dbl;
  686. register SANE_EXT       *p_ext;
  687. SANE_EXT ext_align;
  688. char *c_align;
  689. int i;
  691.    p_dbl = (IEEE_DBL *) pd;
  692.    p_ext = &ext_align;
  693.    top2bits = p_dbl->hi & 0xc0000000L;
  694.    p_ext->l1 = ((p_dbl->hi >> 4) & 0x3ff0000L) | top2bits;
  695.    p_ext->l1 |= ((p_dbl->hi >> 5) & 0x7fffL) | 0x8000L;
  696.    p_ext->l2 = (p_dbl->hi << 27) & 0xf8000000L;
  697.    p_ext->l2 |= ((p_dbl->lo >> 5) & 0x07ffffffL);
  698.    ps2 = (unsigned short *) & (p_dbl->lo);
  699.    ps2++;
  700.    p_ext->s1 = (*ps2 << 11) & 0xf800;
  702.    c_align = (char *) p_ext;
  703.    for (i=0;i<10;i++)
  704.  ps[i] = c_align[i];
  706. #endif
  708. }
  711. /*****************************************************************************
  712. *
  713. *  extended_to_double()
  714. *
  715. *  Purpose:     Convert from SANE extended format to IEEE double precision
  716. *               format.
  717. *
  718. *  Passed:      Pointer to the Apple SANE extended format value and a pointer
  719. *               to what will hold the the IEEE double precision number.
  720. *
  721. *  Outputs:     The IEEE double precision format pointer will be filled with
  722. *               the converted value.
  723. *
  724. *  Returned:    Nothing.
  725. *
  726. *****************************************************************************/
  728. void    extended_to_double(char ps[10], double *pd)
  729. {
  731. #ifdef  MACINTOSH
  733.    x80tox96((extended *) ps, pd);
  735. #else
  737. /* fixed bus alignment error, HP 27-may-93 */
  739. register unsigned long  top2bits;
  741. register IEEE_DBL       *p_dbl;
  742. register SANE_EXT       *p_ext;
  743. SANE_EXT ext_align;
  744. char *c_align;
  745. int i;
  747.    p_dbl = (IEEE_DBL *) pd;
  748.    p_ext = &ext_align;
  750.    c_align = (char *) p_ext;
  751.    for (i=0;i<10;i++)
  752.  c_align[i] = ps[i];
  754.    top2bits = p_ext->l1 & 0xc0000000L;
  755.    p_dbl->hi = ((p_ext->l1 << 4) & 0x3ff00000L) | top2bits;
  756.    p_dbl->hi |= (p_ext->l1 << 5) & 0xffff0L;
  757.    p_dbl->hi |= (p_ext->l2 >> 27) & 0x1f;
  758.    p_dbl->lo = (p_ext->l2 << 5) & 0xffffffe0L;
  759.    p_dbl->lo |= (unsigned long) ((p_ext->s1 >> 11) & 0x1f);
  761. #endif
  763. }
  765. /*****************************************************************************
  766. *
  767. *  Read Audio Interchange File Format (AIFF) headers.
  768. *
  769. *****************************************************************************/
  771. int aiff_read_headers(FILE *file_ptr,
  772.       IFF_AIFF *aiff_ptr, 
  773.       int *byte_per_sample)
  774. {
  776. register char   i;
  777. register long   seek_offset;
  779. char temp_sampleRate[10];
  780. char *dummy;
  781. char holder;
  782. Chunk           FormChunk;
  783. CommonChunk     CommChunk;
  784. SoundDataChunk  SndDChunk;
  785. identifier ident;
  787.    if (fseek(file_ptr, 0, SEEK_SET) != 0)
  788.   return(-1);
  790.    if (fread(&FormChunk, sizeof(Chunk), 1, file_ptr) != 1)
  791.   return(-1);
  793. #ifdef MSDOS
  795.   holder = FormChunk.ckID[0];
  796.   FormChunk.ckID[0] = FormChunk.ckID[3];
  797.   FormChunk.ckID[3] = holder;
  798.   holder = FormChunk.ckID[1];
  799.   FormChunk.ckID[1] = FormChunk.ckID[2];
  800.   FormChunk.ckID[2] = holder;
  802. /* fixed bug in next line, HP 27-may-93 */
  803.   holder = FormChunk.formType[0];
  804.   FormChunk.formType[0] = FormChunk.formType[3];
  805.   FormChunk.formType[3] = holder;
  806.   holder = FormChunk.formType[1];
  807.   FormChunk.formType[1] = FormChunk.formType[2];
  808.   FormChunk.formType[2] = holder;
  810.   FormChunk.ckSize = _lrotl(FormChunk.ckSize, 8);
  811. #endif
  814. /* fixed bug in next line, HP 27-may-93 */
  815. if (strncmp(FormChunk.ckID,IFF_ID_FORM,4) != 0 ||
  816.  strncmp(FormChunk.formType,IFF_ID_AIFF,4) != 0)
  817.       return(-1);  /* warning: different levels of indirection.7/8/92.sr*/
  819. if (fread(&ident,sizeof(identifier), 1, file_ptr) != 1)
  820. return(-1);
  822. #ifdef MSDOS
  824. holder = ident.name[0];
  825. ident.name[0] = ident.name[3];
  826. ident.name[3] = holder;
  827. holder = ident.name[1];
  828. ident.name[1] = ident.name[2];
  829. ident.name[2] = holder;
  831. ident.ck_length = _lrotl(ident.ck_length, 8);
  833. #endif
  835. /* fixed bug in next line, HP 27-may-93 */
  836. while(strncmp(ident.name,IFF_ID_COMM,4) != 0)
  837. {
  838. dummy = calloc( ident.ck_length, sizeof(char));
  839. /* changed "fread( &dummy," to "fread ( dummy,", HP 26-may-93 */
  840. if(fread( dummy, ident.ck_length, 1, file_ptr) != 1)
  841.   return(-1);
  842. free(dummy);
  843. /* fixed bug in next line, HP 27-may-93 */
  844. if(fread( &ident, sizeof(identifier),1, file_ptr) != 1)
  845.   return(-1);
  847. #ifdef MSDOS
  849. holder = ident.name[0];
  850. ident.name[0] = ident.name[3];
  851. ident.name[3] = holder;
  852. holder = ident.name[1];
  853. ident.name[1] = ident.name[2];
  854. ident.name[2] = holder;
  856. ident.ck_length = _lrotl(ident.ck_length, 8);
  858. #endif
  859. }
  861. for( i = 0; i < 4; ++i)
  862. CommChunk.ckID[i] = ident.name[i];
  864. CommChunk.ckSize = ident.ck_length;
  866. if (fread(&CommChunk.numChannels, sizeof(short), 1, file_ptr) != 1)
  867. return(-1);
  869. if (fread(&CommChunk.numSampleFrames, sizeof(unsigned long), 1,
  870.  file_ptr) != 1)
  871. return(-1);
  873. if (fread(&CommChunk.sampleSize, sizeof(short), 1, file_ptr) != 1)
  874. return(-1);
  876. if (fread(CommChunk.sampleRate, sizeof(char[10]), 1, file_ptr) != 1)
  877. return(-1);
  879. #ifdef MSDOS
  881.   CommChunk.sampleSize = _rotl(CommChunk.sampleSize, 8);
  882.   CommChunk.ckSize = _lrotl(CommChunk.ckSize, 8);
  883.   CommChunk.numChannels = _rotl(CommChunk.numChannels, 8);
  884.   CommChunk.numSampleFrames = _lrotl(CommChunk.numSampleFrames, 8);
  886. #endif
  889.   *byte_per_sample = ceil((double)CommChunk.sampleSize / 8);
  891.   for (i = 0; i < sizeof(char[10]); i++)
  892. temp_sampleRate[i] = CommChunk.sampleRate[i];
  894.   extended_to_double(temp_sampleRate, &aiff_ptr->sampleRate);
  897. /* to start the search again from the beginning, HP 27-may-93 */
  898.          fseek (file_ptr, sizeof(Chunk), SEEK_SET);
  900.  if (fread(&ident, sizeof(identifier), 1, file_ptr) != 1)
  901. return(-1);
  903. #ifdef MSDOS
  905. holder = ident.name[0];
  906. ident.name[0] = ident.name[3];
  907. ident.name[3] = holder;
  908. holder = ident.name[1];
  909. ident.name[1] = ident.name[2];
  910. ident.name[2] = holder;
  912. ident.ck_length = _lrotl(ident.ck_length, 8);
  914. #endif
  916. /* fixed bug in next line, HP 27-may-93 */
  917. while(strncmp(ident.name,IFF_ID_SSND,4) != 0)
  918. {
  919. dummy = calloc( ident.ck_length, sizeof(char));
  920. /* changed "fread( &dummy," to "fread ( dummy,", HP 26-may-93 */
  921. if(fread( dummy, ident.ck_length, 1, file_ptr) != 1)
  922.   return(-1);
  923. free(dummy);
  924. if(fread( &ident, sizeof(identifier),1, file_ptr) != 1)
  925.   return (-1);
  926. /* the following lines are not necessary, HP 27-may-93 */
  927. /*
  928.   {
  929.   fseek(file_ptr, 0, SEEK_SET);
  930.   if(fread( &ident, sizeof(identifier), 1, file_ptr) != 1)
  931.   return(-1);
  932.   }
  933. */
  935. #ifdef MSDOS
  938. holder = ident.name[0];
  939. ident.name[0] = ident.name[3];
  940. ident.name[3] = holder;
  941. holder = ident.name[1];
  942. ident.name[1] = ident.name[2];
  943. ident.name[2] = holder;
  945. ident.ck_length = _lrotl(ident.ck_length, 8);
  947. #endif
  948. }
  950. for(i = 0; i < 4; ++i)
  951. SndDChunk.ckID[i] = ident.name[i];
  953. SndDChunk.ckSize = ident.ck_length;
  955. if (fread(&SndDChunk.offset, sizeof(unsigned long), 1, file_ptr) != 1)
  956. return(-1);
  958. if (fread(&SndDChunk.blockSize, sizeof(unsigned long), 1,
  959.              file_ptr) != 1)
  960. return(-1);
  962. #ifdef MSDOS
  964.   SndDChunk.offset = _lrotl(SndDChunk.offset, 8);
  965.   SndDChunk.blockSize = _lrotl(SndDChunk.blockSize, 8);
  967. #endif
  969. /* why seek behinde the SSND Chunk ????, HP 27-may-93 */
  970. /*
  971.    seek_offset = SndDChunk.ckSize - sizeof(SoundDataChunk) +
  972.       sizeof(ChunkHeader);
  974.    if (fseek(file_ptr, seek_offset, SEEK_CUR) != 0)
  975.   return(-1);
  976. */
  978.    aiff_ptr->numChannels       = CommChunk.numChannels;
  979.    aiff_ptr->numSampleFrames   = CommChunk.numSampleFrames;
  980.    aiff_ptr->sampleSize        = CommChunk.sampleSize;
  981.    aiff_ptr->blkAlgn.offset    = SndDChunk.offset;
  982.    aiff_ptr->blkAlgn.blockSize = SndDChunk.blockSize;
  983.    strncpy(aiff_ptr->sampleType, SndDChunk.ckID, 4);
  985.    return(0);
  987. }
  990. /*****************************************************************************
  991. *
  992. *  Seek past some Audio Interchange File Format (AIFF) headers to sound data.
  993. *
  994. *****************************************************************************/
  996. int aiff_seek_to_sound_data(FILE *file_ptr)
  997. {
  999. /*   if (fseek(file_ptr, sizeof(Chunk) + sizeof(SoundDataChunk), SEEK_SET) != 0) */
  1000.    if (fseek(file_ptr, sizeof(Chunk) + sizeof(CommonChunk) + sizeof(SoundDataChunk), SEEK_SET) != 0)
  1001.       return(-1);
  1002.    else
  1003.       return(0);
  1005. }
  1007. /*******************************************************************************
  1008. *
  1009. *  Write Audio Interchange File Format (AIFF) headers.
  1010. *
  1011. *******************************************************************************/
  1013. int aiff_write_headers(FILE *file_ptr,
  1014.        IFF_AIFF *aiff_ptr)
  1015. {
  1017. register char   i;
  1018. register long   seek_offset;
  1020. char            temp_sampleRate[10];
  1022. Chunk           FormChunk;
  1023. CommonChunk     CommChunk;
  1024. SoundDataChunk  SndDChunk;
  1027. strcpy( FormChunk.ckID, IFF_ID_FORM);
  1028. strcpy( FormChunk.formType, IFF_ID_AIFF);
  1029. strcpy( CommChunk.ckID, IFF_ID_COMM);  /*7/7/93,SR,changed FormChunk to CommChunk*/
  1032.    double_to_extended(&aiff_ptr->sampleRate, temp_sampleRate);
  1034.    for (i = 0; i < sizeof(char[10]); i++)
  1035.       CommChunk.sampleRate[i] = temp_sampleRate[i];
  1037.    CommChunk.numChannels             = aiff_ptr->numChannels;
  1038.    CommChunk.numSampleFrames         = aiff_ptr->numSampleFrames;
  1039.    CommChunk.sampleSize              = aiff_ptr->sampleSize;
  1040.    SndDChunk.offset                  = aiff_ptr->blkAlgn.offset + 2;
  1041.    SndDChunk.blockSize               = aiff_ptr->blkAlgn.blockSize;
  1042.    strncpy(/*(unsigned long *)*/ SndDChunk.ckID, aiff_ptr->sampleType,  4);
  1044.    CommChunk.ckSize = sizeof(CommChunk.numChannels) +
  1045.       sizeof(CommChunk.numSampleFrames) + sizeof(CommChunk.sampleSize) +
  1046.       sizeof(CommChunk.sampleRate);
  1048.    SndDChunk.ckSize = sizeof(SoundDataChunk) - sizeof(ChunkHeader) +
  1049.       (CommChunk.sampleSize + BITS_IN_A_BYTE - 1) / BITS_IN_A_BYTE *
  1050.       CommChunk.numChannels * CommChunk.numSampleFrames;
  1052.    FormChunk.ckSize = sizeof(Chunk) + SndDChunk.ckSize + sizeof(ChunkHeader) +
  1053.   CommChunk.ckSize;
  1055.    if (fseek(file_ptr, 0, SEEK_SET) != 0)
  1056.       return(-1);
  1058.    if (fwrite(&FormChunk, sizeof(Chunk), 1, file_ptr) != 1)
  1059.       return(-1);
  1060. /*
  1061.    if (fwrite(&SndDChunk, sizeof(SoundDataChunk), 1, file_ptr) != 1)
  1062.       return(-1);
  1064.    seek_offset = SndDChunk.ckSize - sizeof(SoundDataChunk) +
  1065.       sizeof(ChunkHeader);
  1067.    if (fseek(file_ptr, seek_offset, SEEK_CUR) != 0)
  1068.       return(-1);
  1069. */
  1071.    if (fwrite(CommChunk.ckID, sizeof(ID), 1, file_ptr) != 1)
  1072.       return(-1);
  1074.    if (fwrite(&CommChunk.ckSize, sizeof(long), 1, file_ptr) != 1)
  1075.       return(-1);
  1077.    if (fwrite(&CommChunk.numChannels, sizeof(short), 1, file_ptr) != 1)
  1078.       return(-1);
  1080.    if (fwrite(&CommChunk.numSampleFrames, sizeof(unsigned long), 1,
  1081.       file_ptr) != 1)
  1082.   return(-1);
  1084.    if (fwrite(&CommChunk.sampleSize, sizeof(short), 1, file_ptr) != 1)
  1085.   return(-1);
  1087.    if (fwrite(CommChunk.sampleRate, sizeof(char[10]), 1, file_ptr) != 1)
  1088.       return(-1);
  1090.    /* 960815 FdB put the sound data chunk after the common chunk */
  1091.    if (fwrite(&SndDChunk, sizeof(SoundDataChunk), 1, file_ptr) != 1)
  1092.       return(-1);
  1094.    return(0);
  1096. }
  1098. /*****************************************************************************
  1099. *
  1100. *  bit_stream.c package
  1101. *  Author:  Jean-Georges Fritsch, C-Cube Microsystems
  1102. *
  1103. *****************************************************************************/
  1105. /********************************************************************
  1106.   This package provides functions to write (exclusive or read)
  1107.   information from (exclusive or to) the bit stream.
  1109.   If the bit stream is opened in read mode only the get functions are
  1110.   available. If the bit stream is opened in write mode only the put
  1111.   functions are available.
  1112. ********************************************************************/
  1114. static char *he = "0123456789ABCDEF";
  1116. /* open the device to read the bit stream from it */
  1117. int open_bit_stream_r (Bit_stream *bs,   /* bit stream structure */
  1118.        char *bs_filenam,       /* name of the bit stream file */
  1119.        int size                /* size of the buffer */)
  1120. {
  1121.    register unsigned long n;
  1122.    register int i=0,j=0;
  1123.    register unsigned char flag = 1;
  1124.    unsigned char val;
  1126.    if ((bs->pt = fopen (bs_filenam, "rb")) == NULL)
  1127.    {
  1128.   printf("Could not find "%s".n", bs_filenam);
  1129.   return (0);
  1130.    }
  1131.    fseek (bs->pt, bs->header_size, 0);
  1133.    do
  1134.    {
  1135.  n = fread (&val, sizeof (unsigned char), 1, bs->pt);
  1136.  switch (val)
  1137.  {
  1138.   case 0x30:
  1139.   case 0x31:
  1140.   case 0x32:
  1141.   case 0x33:
  1142.   case 0x34:
  1143.   case 0x35:
  1144.   case 0x36:
  1145.   case 0x37:
  1146.   case 0x38:
  1147.   case 0x39:
  1148.   case 0x41:
  1149.   case 0x42:
  1150.   case 0x43:
  1151.   case 0x44:
  1152.   case 0x45:
  1153.   case 0x46:  j++;
  1154.   case 0xa:  /* n */
  1155.   break;
  1157.   default: /* detection of an binary character */
  1158.   flag--;
  1159.   i = 300;
  1160.   break;
  1161.  }
  1163.    } while (flag & n);
  1165.    if (flag)
  1166.    {
  1167.   printf ("the bit stream file %s is an ASCII filen", bs_filenam);
  1168.   printf ("ASCII files are not supported in this versionn");
  1169.   exit (1);
  1170.   bs->format = ASCII;
  1171.    }
  1172.    else
  1173.    {
  1174.   bs->format = BINARY;
  1175.   printf ("the bit stream file %s is a BINARY filen", bs_filenam);
  1176.    }
  1178.    fclose (bs->pt);
  1180.    if ((bs->pt = fopen (bs_filenam, "rb")) == NULL)
  1181.    {
  1182.       printf("Could not find "%s".n", bs_filenam);
  1183.       return(0);
  1184.    }
  1185.    fseek (bs->pt, bs->header_size, 0);
  1187.    bs->curpos=0;
  1188.    bs->totbits=0;
  1189.    bs->mode = READ_MODE;
  1190.    bs->eobs = FALSE;
  1191.    return(1);
  1192. }
  1194. /*close the device containing the bit stream after a read process*/
  1195. void close_bit_stream_r(Bit_stream *bs)
  1196. {
  1197.    fclose(bs->pt);
  1198. }
  1200. unsigned int get1bit(Bit_stream *bs)
  1201. {
  1202.    if (bs->curpos < bs->totbits)
  1203.    {
  1204. #ifdef PrintBitDebug
  1205.       printf ("pos: %5d bits: %2d code: %4x val: %5dn",
  1206. bs->curpos,  1, bs->bits[bs->curpos], bs->bits[bs->curpos]);
  1207.       fflush (stdout);
  1208. #endif
  1209.       return (bs->bits[bs->curpos++]);
  1210.    }
  1211.    else
  1212.    {
  1213.       fprintf (stderr, "get1bits: bs->curpos = %4d bs->totbits = %4d bitsn",
  1214.        bs->curpos, bs->totbits);
  1215.       exit (1);
  1216.    }
  1217. }
  1219. unsigned long getbits(Bit_stream *bs, int n)
  1220. {
  1221.    unsigned long val;
  1222.    int i;
  1224.    if (bs->curpos + n <= bs->totbits)
  1225.    {
  1226.       for (i = val = 0; i < n; i++)
  1227.          val = 2*val + bs->bits[bs->curpos++];
  1229. #ifdef PrintBitDebug
  1230.       printf ("pos: %5d bits: %2d code: %4x val: %5dn", 
  1231.  bs->curpos - n,  n, val, val);
  1232.       fflush (stdout);
  1233. #endif
  1235.       return (val);
  1236.    }
  1237.    else
  1238.    {
  1239.       fprintf (stderr, "getbits: n: %4d bs->curpos: %4d bs->totbits = %4d bitsn",
  1240.        n, bs->curpos, bs->totbits);
  1241.       exit (1);
  1242.    }
  1243. }
  1245. void program_information(void)
  1246. {
  1247.  printf("ISO MPEG Audio Subgroup Software Simulation Group (1996)n");
  1248.  printf("ISO 13818-3 MPEG-2 Audio Multichannel Decodern");
  1249.  printf("%sn", VERSION);
  1250. }
  1252. /*******************************************************************/
  1254. /*return the status of the bit stream*/
  1255. /* returns 1 if end of bit stream was reached */
  1256. /* returns 0 if end of bit stream was not reached */
  1257. int end_bs(Bit_stream *bs)
  1258. {
  1259.   return(bs->eobs);
  1260. }
  1262. static void bytes_to_bits (char *w_code, unsigned char *in, int nbytes)
  1263. {
  1264.   int i, j, bpos, d;
  1266.   bpos = 0;
  1267.   for (i = 0; i < nbytes; i++) {
  1268.     d = in[i];
  1269.     for (j = 7; j >= 0; j--) {
  1270.       w_code[bpos+j] = d & 1;
  1271.       d = d / 2;
  1272.     }
  1273.     bpos += 8;
  1274.   }
  1275. }
  1277. /*this function seeks for a byte aligned sync word in the bit stream and
  1278.   places the bit stream pointer right after the sync.
  1279.   This function returns 1 if the sync was found otherwise it returns 0 */
  1281. int seek_sync_mpg (Bit_stream *bs)   /* bit stream structure */
  1282. {
  1283.    unsigned long val = 0;
  1284.    unsigned char byte;
  1285.    unsigned char bytes[2000]; /* bytes of an MPEG-1 frame */
  1286.    long sync = SYNC_WORD; /* sync word maximum 32 bits */
  1287.    int N = SYNC_WORD_LNGTH; /* sync word length */
  1288.    int i, sync_bytes = N / 8;
  1289.    int version, layer, br_index, bit_rate, c, f_sampl, padding, slots, nbytes;
  1291.    if (fread (bytes, 1, 4, bs->pt) != 4)
  1292.    {
  1293.       printf ("next mpg header not foundn");
  1294.       return (0);
  1295.    }
  1297.    for (i = 0; i < sync_bytes; i++)
  1298.    {
  1299.       val <<= 8;
  1300.       val |= bytes[i];
  1301.    }
  1302.    if ((N % 8) != 0)
  1303.    {
  1304.       val <<= N % 8;
  1305.       byte = bytes[sync_bytes] >> (8 - (N % 8));
  1306.       val |= byte;
  1307.    }
  1308.   
  1309.    if (val != sync)
  1310.    {
  1311.       printf ("no mpg sync found: %4xn", val);
  1312.       return (0);
  1313.    }
  1315.    bytes_to_bits (bs->bits, bytes, 4);
  1316.   
  1317.    version = bs->bits[12];
  1318.    if (version == 0)
  1319.    {
  1320.       fprintf (stderr, "Error: Illegal version bit.n");
  1321.       exit (1);
  1322.    }
  1323.    layer = 4 - 2*bs->bits[13] - bs->bits[14];
  1324.    for (i = 16, br_index = 0; i < 20; i++)
  1325.       br_index = 2*br_index + bs->bits[i];
  1326.    if (br_index == 0 || br_index == 15)
  1327.    {
  1328.       fprintf (stderr, "Error: Illegal bit_rate index (0 or 15).n");
  1329.       exit (1);
  1330.    }
  1331.    bit_rate = bitrate[layer-1][br_index] * 1000;
  1332.    c = 2*bs->bits[20] + bs->bits[21];
  1333.    switch (c) {
  1334.    case 0: f_sampl = 44100;
  1335.    break;
  1336.    case 1: f_sampl = 48000;
  1337.    break;
  1338.    case 2: f_sampl = 32000;
  1339.    break;
  1340.    case 3: fprintf (stderr, "Error: reserved sampling rate.n");
  1341.    exit (1);
  1342.    }
  1343.    padding = bs->bits[22];
  1345.    if (layer == 1)
  1346.       slots = bit_rate * 12 / f_sampl;
  1347.    else
  1348.       slots = bit_rate * 144 / f_sampl;
  1350.    slots += padding;
  1352.    if (layer == 1)
  1353.       nbytes = slots * 4;
  1354.    else
  1355.       nbytes = slots;
  1357.    if (fread (&bytes[4], 1, nbytes-4, bs->pt) != nbytes-4)
  1358.    {
  1359.       printf ("no n-4 bytesn");
  1360.       return (0);
  1361.    }
  1363.    bytes_to_bits (&bs->bits[32], &bytes[4], nbytes-4);
  1365.    bs->curpos = 12;
  1366.    bs->totbits = nbytes * 8;
  1367.    bs->eobs = 0;
  1369.    return (1);
  1370. }
  1372. int seek_sync_ext (Bit_stream *bs, frame_params *fr_ps)   /* bit stream structure */
  1373. {
  1374.    layer *info = fr_ps->header;
  1375.    unsigned long val = 0;
  1376.    unsigned char byte;
  1377.    unsigned char bytes[2048]; /* bytes of an MPEG-2 extention frame */
  1378.    long sync = EXT_SYNCWORD; /* sync word maximum 32 bits */
  1379.    int N = SYNC_WORD_LNGTH; /* sync word length */
  1380.    int i, sync_bytes = N / 8;
  1381.    int nbytes, lcrc, la_bytes, la_bits;
  1383.    if (fread (bytes, 1, 5, bs->pt) != 5)
  1384.    {
  1385.       printf ("next ext header not foundn");
  1386.       return (0);
  1387.    }
  1389.    for (i = 0; i < sync_bytes; i++)
  1390.    {
  1391.       val <<= 8;
  1392.       val |= bytes[i];
  1393.    }
  1394.    if ((N % 8) != 0)
  1395.    {
  1396.       val <<= N % 8;
  1397.       byte = bytes[sync_bytes] >> (8 - (N % 8));
  1398.       val |= byte;
  1399.    }
  1400.   
  1401.    if (val != sync)
  1402.    {
  1403.       printf ("no ext sync found: %4xn", val);
  1404.       return (0);
  1405.    }
  1407.    nbytes = (((int) (bytes[3] & 0x0f) * 256) + bytes[4]) / 2;
  1409.    if (fread (&bytes[5], 1, nbytes-5, bs->pt) != nbytes-5)
  1410.       return (0);
  1411.    bytes_to_bits (bs->bits, bytes, nbytes);
  1413.    bs->totbits = nbytes * 8;
  1414.    bs->curpos = 12;
  1415.    bs->eobs = 0;
  1417.    info->ext_crc_check = getbits (bs, 16);
  1418.    info->ext_length = getbits (bs, 11);
  1419.    info->reserved_bit = get1bit (bs);
  1421. #ifdef EXTENSION_CRC_CHECK
  1423.    lcrc = info->ext_length * 8 - 28;
  1424.    if (lcrc > 128) lcrc = 128;
  1425.    la_bytes = (lcrc - 12) / 8;
  1426.    la_bits = (lcrc - 12) % 8;
  1428.    for (i = 0; i < la_bytes; i++)
  1429.       info->ext_crc_bits[i] = getbits (bs, 8);
  1430.    info->ext_crc_bits[la_bytes] = getbits (bs, la_bits);
  1431.    bs->curpos = 40;
  1433.    if (mc_ext_error_check (fr_ps, la_bytes, la_bits) == 0)
  1434.       return (0);
  1435. #endif
  1437.    return (1);
  1438. }
  1440. /*****************************************************************************
  1441. *
  1442. *  End of bit_stream.c package
  1443. *
  1444. *****************************************************************************/
  1446. /*****************************************************************************
  1447. *
  1448. *  CRC error protection package
  1449. *
  1450. *****************************************************************************/
  1452. void update_CRC(unsigned int data,
  1453. unsigned int length,
  1454. unsigned int *crc)
  1455. {
  1456. unsigned int  masking, carry;
  1458. masking = 1 << length;
  1460. while((masking >>= 1)){
  1461. carry = *crc & 0x8000;
  1462. *crc <<= 1;
  1463. if (!carry ^ !(data & masking))
  1464. *crc ^= CRC16_POLYNOMIAL;
  1465. }
  1466. *crc &= 0xffff;
  1468. #ifdef PrintCRCDebug
  1469. printf ("crc_len: %2d code: %4x crc: %4xn", length, data, *crc);
  1470. fflush (stdout);
  1471. #endif
  1472. }
  1474. void I_CRC_calc (frame_params *fr_ps,
  1475.  unsigned int bit_alloc[7][SBLIMIT], 
  1476.  unsigned int *crc)
  1477. {
  1478.    int i, k;
  1479.    layer *info = fr_ps->header;
  1480.    int stereo  = fr_ps->stereo;
  1481.    int jsbound = fr_ps->jsbound;
  1483.    *crc = 0xffff; /* changed from '0' 92-08-11 shn */
  1484.    update_CRC (info->bitrate_index, 4, crc);
  1485.    update_CRC (info->sampling_frequency, 2, crc);
  1486.    update_CRC (info->padding, 1, crc);
  1487.    update_CRC (info->extension, 1, crc);
  1488.    update_CRC (info->mode, 2, crc);
  1489.    update_CRC (info->mode_ext, 2, crc);
  1490.    update_CRC (info->copyright, 1, crc);
  1491.    update_CRC (info->original, 1, crc);
  1492.    update_CRC (info->emphasis, 2, crc);
  1494.    for (i=0; i<SBLIMIT; i++)
  1495.       for (k=0; k<((i<jsbound)?stereo:1); k++)
  1496.  update_CRC (bit_alloc[k][i], 4, crc);
  1497. }
  1500. void II_CRC_calc(frame_params *fr_ps,
  1501.  unsigned int bit_alloc[7][SBLIMIT], 
  1502.  unsigned int scfsi[7][SBLIMIT],
  1503.  unsigned int *crc)
  1504. {
  1505.    int i, k;
  1506.    layer *info = fr_ps->header;
  1507.    int stereo  = fr_ps->stereo;
  1508.    int sblimit = fr_ps->sblimit;
  1509.    int jsbound = fr_ps->jsbound;
  1510.    al_table *alloc = fr_ps->alloc;
  1512.    *crc = 0xffff; /* changed from '0' 92-08-11 shn */
  1513.    update_CRC (info->bitrate_index, 4, crc);
  1514.    update_CRC (info->sampling_frequency, 2, crc);
  1515.    update_CRC (info->padding, 1, crc);
  1516.    update_CRC (info->extension, 1, crc);
  1517.    update_CRC (info->mode, 2, crc);
  1518.    update_CRC (info->mode_ext, 2, crc);
  1519.    update_CRC (info->copyright, 1, crc);
  1520.    update_CRC (info->original, 1, crc);
  1521.    update_CRC (info->emphasis, 2, crc);
  1523.    for (i=0; i<sblimit; i++)
  1524.       for (k=0; k<((i<jsbound)?stereo:1); k++)
  1525.  update_CRC (bit_alloc[k][i], (*alloc)[i][0].bits, crc);
  1527.    for (i=0; i<sblimit; i++)
  1528.       for (k=0; k<stereo; k++)
  1529.  if (bit_alloc[k][i])
  1530.     update_CRC (scfsi[k][i], 2, crc);
  1531. }
  1534. void mc_error_check(frame_params *fr_ps,
  1535.     unsigned int bit_alloc[7][SBLIMIT], 
  1536.     unsigned int scfsi[7][SBLIMIT],
  1537.     unsigned int *crc,
  1538.      int ch_start,
  1539.      int channels)
  1540. {
  1541.  int i, m, l, pci,c,sbgr;
  1542.  layer *info = fr_ps->header;
  1543.  int sblimit = fr_ps->sblimit_mc;
  1544.  al_table *alloc = fr_ps->alloc_mc;
  1546.  *crc = 0xffff; /* changed from '0' 92-08-11 shn */
  1547.  update_CRC(info->ext_bit_stream_present, 1, crc);
  1548.  if ( info->ext_bit_stream_present == 1)
  1549. update_CRC(info->n_ad_bytes, 8, crc);
  1551.  update_CRC(info->center, 2, crc);
  1552.  update_CRC(info->surround, 2, crc);
  1553.  update_CRC(info->lfe, 1, crc);
  1554.  update_CRC(info->audio_mix, 1,crc);
  1555.  update_CRC(info->dematrix_procedure, 2, crc);
  1556.  update_CRC(info->no_of_multi_lingual_ch, 3, crc);
  1557.  update_CRC(info->multi_lingual_fs, 1, crc);
  1558.  update_CRC(info->multi_lingual_layer, 1, crc);
  1559.  update_CRC(info->copyright_ident_bit, 1, crc);
  1560.  update_CRC(info->copyright_ident_start, 1, crc);
  1561.  
  1562.  update_CRC(info->tc_sbgr_select, 1, crc);
  1563.  update_CRC(info->dyn_cross_on, 1, crc);
  1564.  update_CRC(info->mc_prediction_on, 1, crc);
  1566.  if (info->tc_sbgr_select == 1)
  1567.     update_CRC (info->tc_allocation, fr_ps->alloc_bits, crc);
  1568.  else
  1569.     for (i = 0; i < 12; i++)
  1570.        update_CRC (info->tc_alloc[i], fr_ps->alloc_bits, crc);
  1571.    
  1572.  if (info->dyn_cross_on == 1)
  1573.  {
  1574. update_CRC(info->dyn_cross_LR, 1, crc);
  1575. for(i = 0; i < 12; i++)
  1576. {
  1577. update_CRC(info->dyn_cross_mode[i], fr_ps->dyn_cross_bits, crc);
  1578. if (info->surround == 3)
  1579. update_CRC(info->dyn_second_stereo[i], 1, crc);
  1580. }
  1581.  }
  1583.  if (info->mc_prediction_on == 1)
  1584.  {
  1585. for(i = 0; i < 8; i++)
  1586. {  /* new sbgr < 8 */
  1587. update_CRC(info->mc_prediction[i], 1, crc);
  1588. if( info->mc_prediction[i] == 1 )
  1589. {
  1590.    for(pci=0; pci< pred_coef_table[fr_ps->pred_mode][info->dyn_cross_mode[i]]; pci++)
  1591. update_CRC(info->mc_predsi[i][pci], 2, crc);
  1592. }
  1593. }
  1594.  }
  1596.  if( info->lfe )
  1597. update_CRC(info->lfe_alloc, (*alloc)[0][0].bits, crc);
  1599.  for (i=0;i<sblimit;i++)
  1600.    if( info->dyn_cross_on == 0)
  1601. for(m = ch_start; m < channels; ++m)
  1602. {
  1603. if((i < 12) || (m != 2) || (fr_ps->header->center != 3))
  1604. update_CRC(bit_alloc[m][i], (*alloc)[i][0].bits, crc);
  1605. /* to limit the CRC calc., if the phant. center valid */
  1606. }
  1607.    else
  1608.    {
  1609. if(i == 0) sbgr = 0;
  1610. else
  1611.        for(c = 1; c < 12; c++)
  1612.        {
  1613.        if((sb_groups[c-1] < i) && (i <= sb_groups[c]))
  1614.        {
  1615.       sbgr = c;  /* search the valid subband group */
  1616.       break;
  1617.           }
  1618.    }
  1620. /* 960816 FdB new setup for dyn. crosstalk modes */
  1621. if (info->dyn_cross_mode[sbgr] == 0)
  1622. {
  1623.    for (m = ch_start; m < channels; m++)
  1624.       if ((info->center != 3) || (i < 12) || (m !=2))
  1625.  if (info->surround == 3 && info->dyn_second_stereo[sbgr] == 1)
  1626.  {
  1627.     if ((info->center == 0 || m != 4) && (info->center != 0 || m != (ch_start+1)))
  1628.        update_CRC (bit_alloc[m][i], (*alloc)[i][0].bits, crc);
  1629.  }
  1630.  else
  1631.     update_CRC (bit_alloc[m][i], (*alloc)[i][0].bits, crc);
  1632. }
  1633. else if (fr_ps->dyn_cross_bits == 1)   /* for channel mode 3/0 and 2/1 */
  1634. {
  1635.    if (info->surround == 3) /* 3/0 + 2/0 */
  1636.    {
  1637.       update_CRC (bit_alloc[3][i], (*alloc)[i][0].bits, crc);
  1638.       if (info->dyn_second_stereo[sbgr] == 0)
  1639.          update_CRC (bit_alloc[4][i], (*alloc)[i][0].bits, crc);
  1640.    }
  1641. }
  1642. else if (fr_ps->dyn_cross_bits == 3)   /* for channel mode 3/1 and 2/2 */
  1643. {
  1644.    if ((info->center != 3) || (i < 12))    /* 3/2 + no phantom center */
  1645.       if ((info->dyn_cross_mode[sbgr] == 1) || (info->dyn_cross_mode[sbgr] == 4))
  1646.          update_CRC (bit_alloc[2][i], (*alloc)[i][0].bits, crc);
  1648.    if (info->dyn_cross_mode[sbgr] == 2)
  1649.       update_CRC (bit_alloc[3][i], (*alloc)[i][0].bits, crc);
  1650. } 
  1651. else if (fr_ps->dyn_cross_bits == 4)   /* for channel mode 3/2 */
  1652. {
  1653.    /* T2 */
  1654.    if ((info->center != 3) || (i < 12))    /* 3/2 + no phantom center */
  1655.       switch (info->dyn_cross_mode[sbgr])
  1656.    {
  1657.    case 1 :   /* T2 contains bitalloc info */
  1658.    case 2 :
  1659.    case 4 :
  1660.    case 8 :
  1661.    case 9 :
  1662.    case 10: 
  1663.    case 11: 
  1664.    case 12: 
  1665.    case 14: 
  1666.       update_CRC (bit_alloc[2][i], (*alloc)[i][0].bits, crc);
  1667.       break;
  1668.    }
  1670.    /* T3 */
  1671.    switch (info->dyn_cross_mode[sbgr])
  1672.    {
  1673.    case 1 :   /* T3 contains bitalloc info */
  1674.    case 3 :
  1675.    case 5 :
  1676.    case 8 :
  1677.    case 10: 
  1678.    case 13: 
  1679.       update_CRC (bit_alloc[3][i], (*alloc)[i][0].bits, crc);
  1680.       break;
  1681.    }
  1683.    /* T4 */
  1684.    switch (info->dyn_cross_mode[sbgr])
  1685.    {
  1686.    case 2 :   /* T4 contains bitalloc info */
  1687.    case 3 :
  1688.    case 6 :
  1689.    case 9 : 
  1690.       update_CRC (bit_alloc[4][i], (*alloc)[i][0].bits, crc);
  1691.       break;
  1692.    }
  1693. }
  1694.    }
  1696.      for (i = 0; i < sblimit; i++)
  1697. for (m = ch_start; m < channels; m++)
  1698. if (bit_alloc[m][i])
  1699. update_CRC (scfsi[m][i], 2, crc); 
  1700. }
  1703. #ifdef Augmentation_7ch
  1704. void mc_aug_error_check(frame_params *fr_ps,
  1705.     unsigned int bit_alloc[7][SBLIMIT], 
  1706.     unsigned int scfsi[7][SBLIMIT],
  1707.     unsigned int *crc)
  1708. {
  1709.    int i, m, l, pci,c,sbgr;
  1710.    layer *info = fr_ps->header;
  1711.    int sblimit = fr_ps->sblimit_mc;
  1712.    al_table *alloc = fr_ps->alloc_mc;
  1714.    *crc = 0xffff; /* changed from '0' 92-08-11 shn */
  1716.    update_CRC(info->aug_mtx_proc, 2, crc);
  1717.    update_CRC(info->aug_dyn_cross_on, 1, crc);
  1718.    update_CRC(info->aug_future_ext, 1, crc);
  1720.    if(info->aug_mtx_proc == 0)
  1721.       for(i = 0; i < 12; i++)
  1722.  update_CRC(info->tc_aug_alloc[i], 3, crc);
  1723.    else if(info->aug_mtx_proc == 1)
  1724.       for(i = 0; i < 12; i++)
  1725.  update_CRC(info->tc_aug_alloc[i], 2, crc);
  1726.      
  1727.    if(info->aug_dyn_cross_on == 1)
  1728.       for(i = 0; i < 12; i++)
  1729.  update_CRC(info->dyn_cross_aug_mode[i], 5, crc);
  1731.    for (i=0;i<sblimit;i++)
  1732.    {
  1733.       if(i == 0) sbgr = 0;
  1734.       else
  1735.  for(c = 1; c < 12; c++)
  1736.  {
  1737.     if((sb_groups[c-1] < i) && (i <= sb_groups[c]))
  1738.     {
  1739.       sbgr = c;  /* search the valid subband group */
  1740.       break;
  1741.     }
  1742.  }
  1743.  
  1744.       /* check bitalloc info */
  1745.       switch (info->dyn_cross_aug_mode[sbgr])
  1746.       {
  1747.       case  0: /* T5 and T6 contains bitalloc info */
  1748.        update_CRC(bit_alloc[5][i], (*alloc)[i][0].bits, crc);
  1749.        update_CRC(bit_alloc[6][i], (*alloc)[i][0].bits, crc);
  1750.        break;
  1752.       case  1: /* T5 contains bitalloc info */
  1753.       case  2:
  1754.       case  3:
  1755.       case  4:
  1756.        update_CRC(bit_alloc[5][i], (*alloc)[i][0].bits, crc);
  1757.        break;
  1759.       case  5:   /* T6 contains bitalloc info */
  1760.       case 10:
  1761.       case 14:
  1762.        update_CRC(bit_alloc[6][i], (*alloc)[i][0].bits, crc);
  1763.        break;
  1764.       }
  1765.    }
  1767.    for (i = 0; i < sblimit; i++)
  1768.       for (m = 5; m < 7; m++)
  1769.  if (bit_alloc[m][i])
  1770.  update_CRC (scfsi[m][i], 2, crc); 
  1771. }
  1772. #endif
  1775. int mc_ext_error_check (frame_params *fr_ps, int bytes, int bits)
  1776. {
  1777. layer *info = fr_ps->header;
  1778. int i;
  1779. unsigned int crc;
  1781. crc = 0xffff;
  1782. update_CRC(info->ext_length, 11, &crc);
  1783. update_CRC(info->reserved_bit, 1, &crc);
  1784. for( i = 0; i < bytes; i++)
  1785. update_CRC(info->ext_crc_bits[i], 8, &crc);
  1786. update_CRC(info->ext_crc_bits[bytes], bits, &crc);
  1788. if (crc != info->ext_crc_check)
  1789. {
  1790. printf ("nERROR in EXT.-CRC n");
  1791. return (0);
  1792. }
  1793. else
  1794. return(1);
  1795. }
  1797. /*****************************************************************************
  1798. *
  1799. *  End of CRC error protection package
  1800. *
  1801. *****************************************************************************/
  1803. #ifdef  MACINTOSH
  1804. /*****************************************************************************
  1805. *
  1806. *  Set Macintosh file attributes.
  1807. *
  1808. *****************************************************************************/
  1810. void    set_mac_file_attr(char fileName[MAX_NAME_SIZE],
  1811.   short vRefNum,
  1812.   OsType creator,
  1813.   OsType fileType)
  1814. {
  1816. short   theFile;
  1817. char    pascal_fileName[MAX_NAME_SIZE];
  1818. FInfo   fndrInfo;
  1820. CtoPstr(strcpy(pascal_fileName, fileName));
  1822.         FSOpen(pascal_fileName, vRefNum, &theFile);
  1823.         GetFInfo(pascal_fileName, vRefNum, &fndrInfo);
  1824. fndrInfo.fdCreator = creator;
  1825.         fndrInfo.fdType = fileType;
  1826.         SetFInfo(pascal_fileName, vRefNum, &fndrInfo);
  1827.         FSClose(theFile);
  1829. }
  1830. #endif