home *** CD-ROM | disk | FTP | other *** search
/ Program Metropolis - Software Boutique 95 / SOFTWARECD.iso / stat45 / disk1 / develop.exe / DEMO.C next >
Encoding:
C/C++ Source or Header  |  1992-12-02  |  11.6 KB  |  352 lines
  1. /*===========================================================================*/
  2. /*                                                                           */ 
  3. /* This program demonstrates basic operations on a CSS data file. It will    */
  4. /* read a CSS raw data file, print its basic characteristics (including      */
  5. /* alphanumeric values for the second variable in the data file), and then   */
  6. /* create a new CSS data file with one additional variable.  That variable   */
  7. /* will contain the sum of all variables in the input file.  Arguments (input*/
  8. /* file name, output file name) are passed by command line (from DOS) in the */
  9. /* following manner:                                                         */
  10. /*                                                                           */
  11. /*                       demo input output                                   */
  12. /*                                                                           */
  13. /* This program was developed using a standard subset of the C language, and */
  14. /* it is ANSI compatible.                                                    */
  15. /* Enclosed demo.exe was created with Microsoft C 7.00 using command         */
  16. /*   cl demo.c                                                               */
  17. /*===========================================================================*/
  18.  
  19.  
  20. #include   <stdio.h>
  21. #include   <stdlib.h>
  22.  
  23. #include "demo.h"
  24.  
  25. FILE   *in,   /* input stream */
  26.       *out;   /* output stream */
  27.  
  28. struct   HeadTag   InHeader,   /* header of input file */
  29.                   OutHeader;   /* header of output file */
  30. struct   InfoTag   InArray,    /* information array about input file */
  31.                   OutArray;    /* information array about output file */
  32.  
  33. main(argc,argv)
  34. int argc;
  35. char * argv[];
  36. { 
  37.    /*-- get arguments from command line --*/
  38.    if ( argc <= 2 ) {
  39.       printf("\n To use program enter \ndemo inp_file out_file\n");
  40.       exit(1);
  41.    }
  42.    in = fopen( argv[1], "rb");
  43.    if ( in == NULL ){
  44.       printf("\n Cannot open input file.\n");
  45.       exit(1);
  46.    }
  47.  
  48.    /*-- initializing data structures for the input file --*/
  49.    InitInputFile();
  50.  
  51.    /*-- printing some information about input file --*/
  52.    PrintInputInfo();
  53.  
  54.    /*-- creating second (output) file --*/
  55.    out = fopen( argv[2], "wb");
  56.    if ( in == NULL ){
  57.       printf("\n Cannot create output file.\n");
  58.       exit(1);
  59.    }
  60.  
  61.    /*-- creating information header for output file; writing it to disk --*/
  62.    CreateOutputHeader();
  63.  
  64.    /*-- processing data for output file --*/
  65.    ProcessData();
  66.  
  67.    /*-- closing files and exiting --*/
  68.    fclose(in);
  69.    fclose(out);
  70.    printf("\nNormal program termination.\n");
  71.    exit(0);
  72. }
  73.  
  74. /*------------------------------- end of main ----------------------------------*/
  75.  
  76. void   InitInputFile(void)
  77. /*-- read from disk basic information about input file and store it in
  78.    InHeader and InArray structures --*/
  79. {
  80.    fread( &InHeader, SizeFH, 1, in);
  81.    if ( memcmp( InHeader.Title, "CSS ", 4 ) ||
  82.         (InHeader.Code != -9876L && InHeader.Code != -9877L) ) {
  83.       printf(" This is not a CSS raw data file.");
  84.       exit(1);
  85.    }
  86.    if ( InHeader.NV <= 1 ){
  87.       printf(" The input file for this program must contain at least two variables.");
  88.       exit(1);
  89.    }
  90.    if ( InHeader.NV >= 300 ){
  91.       printf(" Number of variables must not exceed 300. ");
  92.       exit(1);
  93.    }
  94.    switch ((int) InHeader.NArr ){
  95.       case 5 :
  96.          InHeader.SL1Len = 0;
  97.          InHeader.SL2Len = 0;
  98.       case 7 :
  99.          InHeader.SL3Len = 0;
  100.    }
  101.    fseek( in, 24 + InHeader.NArr * 2, SEEK_SET);
  102.    fread( &InArray, SizeFI, 1, in);
  103. }
  104.  
  105. /*------------------------------- end of InitInputFile -------------------------*/
  106.  
  107. void   PrintInputInfo(void)
  108. {
  109.    char Name[8];
  110.    union ValueTag MD;
  111.    long Offset;
  112.    struct FmTag Format;
  113.    int NbAlpha,NbRead,i;
  114.    long OffsAlpha;
  115.    struct LabTag Alpha[10];
  116.  
  117.    printf("\nInput file characteristics:\n Variables: %ld\n Cases: %ld\n Precision: %ldB\n",
  118.           InHeader.NV,InHeader.NC,InHeader.Prec);
  119.    if ( InHeader.Code == -9877L )
  120.       printf("\n This is a STATISTICA/W data file\n");
  121.    /*-- now print specifications for the second variable --*/
  122.    Offset = 24L + 2 * InHeader.NArr + InHeader.NCLen + InHeader.HLen;
  123.    /*-- find name --*/
  124.    fseek( in, Offset + 8, SEEK_SET);
  125.    fread( Name, 8, 1, in);
  126.    Offset += InHeader.VLen;
  127.    fseek( in, Offset + 2, SEEK_SET);
  128.    fread( &Format, 2, 1,in);
  129.    Offset += InHeader.FLen;
  130.    fseek( in, Offset + InHeader.Prec, SEEK_SET);
  131.    fread( &MD, (int)InHeader.Prec, 1, in);
  132.    printf("\nSpecifications for second variable:\n");
  133.    printf(" Name: %.8s\n",Name);
  134.    if ( InHeader.Prec == 4 )
  135.       printf(" MD Code: %f\n",MD.f);
  136.    else
  137.       printf(" MD Code: %lf\n",MD.d);
  138.    printf(" Format: width = %d, decimal = %d\n",(int)Format.Nz,(int)Format.Nd);
  139.    if ( InHeader.NArr >= 7 ){
  140.       /*-- if alphanumeric value labels exist, find out how many--*/
  141.       Offset += InHeader.MDLen;
  142.       fseek( in, Offset + 2, SEEK_SET);
  143.       fread( &NbAlpha, 2, 1, in);
  144.       if ( NbAlpha ){
  145.          /*-- find where they begin --*/
  146.          Offset += InHeader.SL1Len;
  147.          fseek( in, Offset + 4, SEEK_SET);
  148.          fread( &OffsAlpha, 4, 1, in);
  149.          Offset += InHeader.SL2Len;
  150.          /*-- here the offsets to long value labels begin --*/
  151.          Offset += InHeader.SL3Len;
  152.          /*-- here the long variable names begin --*/
  153.          Offset += InArray.MemVar;
  154.          /*-- here alphanumerics starts --*/
  155.          fseek( in, Offset + OffsAlpha, SEEK_SET);
  156.          /*-- read first 10 alphanumeric specification --*/
  157.          fread( Alpha, SizeLT, NbRead = min( NbAlpha, 10), in);
  158.          printf(" Number of alphanumerics: %d\n",NbAlpha);
  159.          for ( i = 0 ; i < NbRead; i++ )
  160.             printf(" Alpha: %.8s Value: %f\n",Alpha[i].T,Alpha[i].V);
  161.       }
  162.       else
  163.          printf(" Second variable contains now alphanumeric value labels\n");
  164.    }
  165.    else
  166.       printf(" There are no alphanumeric value labels in this file.\n");
  167. }
  168.  
  169. /*------------------------------- end of PrintInputInfo ------------------------*/
  170.  
  171. void   CreateOutputHeader(void)
  172. /*-create output file header, make changes to accomodate additional variable-*/
  173. {
  174.    union ValueTag MD;
  175.    int NbAlpha;
  176.    long OffsAlpha,OffsLongAlpha;
  177.  
  178.    /*-- copy basic information --*/
  179.    memcpy( &OutHeader, &InHeader, SizeFH );
  180.    /* we are going to create CSS (not STAT/W) data file;
  181.       STAT/W will read CSS data file */
  182.    OutHeader.Code = -9876L;
  183.    OutHeader.NV++;
  184.    OutHeader.VLen += 8;
  185.    OutHeader.FLen += 2;
  186.    OutHeader.MDLen += OutHeader.Prec;
  187.    if ( OutHeader.SL1Len )
  188.       OutHeader.SL1Len += 2;
  189.    if ( OutHeader.SL2Len )
  190.       OutHeader.SL2Len += 4;
  191.    if ( OutHeader.SL3Len )
  192.       OutHeader.SL3Len += 4;
  193.    fwrite( &OutHeader, (unsigned int)(24 + OutHeader.NArr * 2), 1, out);
  194.  
  195.    /*-- copy information array --*/
  196.    memcpy( &OutArray, &InArray, SizeFI );
  197.    OutArray.ExtraInfo = 0L;
  198.    fwrite( &OutArray, SizeFI, 1, out);
  199.  
  200.    /*-- copy text header --*/
  201.    fseek( in, 24 + OutHeader.NArr * 2 + SizeFI, SEEK_SET );
  202.    fcopy( out, in, 80);
  203.  
  204.    /*-- variable names --*/
  205.    fwrite( "     SUM", 8, 1, out);
  206.    fcopy( out, in, InHeader.VLen);
  207.  
  208.    /*-- variable formats --*/
  209.    fwrite( "", 2, 1, out); /*-- create default 8.3 --*/
  210.    fcopy( out, in, InHeader.FLen);
  211.  
  212.    /*-- MD codes --*/
  213.    if ( InHeader.Prec == 4 )
  214.       MD.f = (float)-9999.0;
  215.    else
  216.       MD.d = -9999.0;
  217.    fwrite( &MD, (int)InHeader.Prec, 1, out);
  218.    fcopy( out, in, InHeader.MDLen );
  219.  
  220.    /*-- number of alphanumeric value labels --*/
  221.    if ( InHeader.SL1Len ){
  222.       NbAlpha = 0;
  223.       fwrite( &NbAlpha, 2, 1, out);
  224.       fcopy( out, in, InHeader.SL1Len );
  225.    }
  226.  
  227.    /*-- offset to alphanumeric value labels --*/
  228.    if ( InHeader.SL2Len ){
  229.       OffsAlpha = 0L;
  230.       fwrite( &OffsAlpha, 4, 1, out);
  231.       fcopy( out, in, InHeader.SL2Len );
  232.    }
  233.  
  234.    /*-- offset to long alphanumeric value labels --*/
  235.    if ( InHeader.SL3Len ){
  236.       OffsLongAlpha = 0L;
  237.       fwrite( &OffsLongAlpha, 4, 1, out);
  238.       fcopy( out, in, InHeader.SL3Len );
  239.    }
  240.  
  241.    /*-- long variable names (labels/formulas) --*/
  242.    fcopy( out, in, (unsigned int)OutArray.MemVar);
  243.  
  244.    /*-- short alphanumeric values for all variables;
  245.       NOTE: make sure that third argument in fcopy is not larger than
  246.       the max unsigned integer value; otherwise write your own fcopy routine --*/
  247.    fcopy( out, in, (unsigned int)OutArray.MemSvl);
  248.  
  249.    /*-- long alphanumerics - same as above --*/
  250.    fcopy( out, in, (unsigned int)OutArray.MemLvl);
  251. }
  252.  
  253. /*----------------------------- end of CreateOutputHeader ----------------------*/
  254.  
  255. void   ProcessData(void)
  256. /*-- copy data, process each case and compute the sum of variables;
  257.    we assume that the file pointers are at the end of header specifications 
  258.    and at the bginning of the data for both files --*/
  259. {
  260.    double * pD,sumD;
  261.    float * pF,sumF;
  262.    long Offset,i;
  263.    int RecordLen,j;
  264.    BYTE * pData;
  265.    char CaseName[20];
  266.  
  267.    /* get offset on input file */
  268.    Offset = ftell( in );
  269.    /* in case of STATISTICA/W we must skip some part of input data file */
  270.    if ( InHeader.Code == -9877L )
  271.       Offset += InArray.ExtraInfo;
  272.    /*-- first move pointers to nearest multiply of 512B --*/
  273.    Ceil512( &Offset );
  274.    fseek( in, Offset, SEEK_SET );
  275.    Offset = ftell( out );
  276.    Ceil512( &Offset );
  277.    fseek( out, Offset, SEEK_SET );
  278.  
  279.    /*-- next allocate memory for one record of data --*/
  280.    if ( (pData = malloc( (unsigned int)(InHeader.Prec*InHeader.NV) )) == NULL ){
  281.       printf("\n Memory allocation error.");
  282.       exit(1);
  283.    }
  284.  
  285.    /*-- process data --*/
  286.    for ( i = 0, RecordLen = (int)(InHeader.Prec * InHeader.NV);
  287.          i < InHeader.NC;
  288.          i++ ){
  289.       if ( InArray.LCName > 0 ){
  290.          fread( CaseName, (int)InArray.LCName, 1, in);
  291.          fwrite( CaseName, (int)InArray.LCName, 1, out);
  292.       }
  293.       fread( pData, RecordLen, 1, in);
  294.       if (InHeader.Prec == 4) {
  295.          /*-- compute sum of variables --*/
  296.          for ( j = 0, sumF = (float)0, pF = (float *)pData;
  297.                j < InHeader.NV;
  298.                j++, pF++ )
  299.             sumF += *pF;
  300.          /*-- write results to output file --*/
  301.          fwrite( &sumF, 4, 1, out);
  302.       }
  303.       else {
  304.          /*-- compute sum of variables --*/
  305.          for ( j = 0, sumD = (double)0, pD = (double *)pData;
  306.                j < InHeader.NV;
  307.                j++, pD++ )
  308.             sumD += *pD;
  309.          /*-- write results to output file --*/
  310.          fwrite( &sumD, 8, 1, out);
  311.       }
  312.       /*-- write remaining part of record to the output --*/
  313.       fwrite( pData, RecordLen, 1, out);
  314.    }
  315.    /*-- free allocated memory buffer --*/
  316.    free(pData);
  317. }
  318.  
  319. /*-------------------------------- end of ProcessData --------------------------*/
  320.  
  321. void fcopy( FILE * out, FILE * in, unsigned int size )
  322. /*-- copy SIZE bytes from input stream to output stream using dynamically
  323.    allocated buffer --*/
  324. {
  325.    BYTE * p;
  326.  
  327.    if ( size == 0 )
  328.       return;
  329.    if ( (p = malloc(size)) == NULL ){
  330.       printf("\n Memory allocation error.");
  331.       exit(1);
  332.    }
  333.    fread( p, size, 1, in);
  334.    fwrite( p, size, 1, out);
  335.    free(p);
  336. }
  337.  
  338. /*--------------------------------- end of fcopy -------------------------------*/
  339.  
  340. void Ceil512( long * pL )
  341. /*-- adjust number to nearest multiple of 512 --*/
  342. {
  343.   ldiv_t S;
  344.  
  345.   S = ldiv( *pL, 512L );
  346.   if ( S.rem != 0L )
  347.      *pL = (S.quot + 1L) * 512L;
  348. }
  349.  
  350. /*--------------------------------- end of Ceil512 ------------------------------*/
  351.  
  352.