/*  $Id: testCedarMemIO.c 3304 2011-01-17 15:25:59Z brideout $ */

/*
modification history
--------------------
00a,21Feb00         original
*/

/* 
* USAGE: testCedarMemIO
* This program test Cedar Memory IO.
* 
*/

#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <sys/types.h>
#include <cedarIO.h>
#include <cedar.h>

int
main (argc, argv)
     int     argc;
     char    *argv[];
{
    int i=0, ind1=0, ind2=0, nrec=0, debug=0;
    FILE *finp, *foutp;
    char madfile[128] = "";
    Int16 *cedarRecordp=(Int16 *)NULL;
    int fileType=0;
    void timestart(), timestop();
    char command[128];
    int recordpInMem = 0;

    /* madrigal state variables */
    int madrigalBlockSize=6720;
    int sigWords=0, blockIndex=0, prevRec=0, thisRec=0;
    Int16 *blockp=(Int16 *)NULL;

    /* Cbf state parameters */
    int cbfBlockSize=4096;

    /* Memory file state parameters */
    Int16 *cedarFile1p=(Int16 *)NULL, *cedarFile2p=(Int16 *)NULL;
    int mempos1, mempos2;
    int fileSize1, fileSize2;

    /* Get the name of the input madrigal file */
    cedarGetMadroot(madfile);
    if (argc == 1) {
        strcat(madfile, "/experiments/1998/mlh/20jan98/mil980120g.002");
    } else if (argc == 2) {
        strcpy(madfile, argv[1]);
    } else {
        return(-1);
    }

    fileType = cedarFileType(madfile, madrigalBlockSize, cbfBlockSize);
    if (debug) 
        (void) printf("%s - fileType = %d\n", madfile,fileType);

    /*** Input file to Madrigal ***/
    if ((finp = fopen(madfile, "r")) == (FILE *)NULL) {
        printf("Input file %s does not exist\n", madfile);
        return(1);
    }

    foutp = fopen("file1", "w");
    sigWords = 0;
    blockIndex = 0;
    blockp=(Int16 *)NULL;
    prevRec = 0;
    thisRec = 0;
    nrec = 0;
    while((ind1 = getNextMadrigalRecord(finp, &cedarRecordp,
                                        madrigalBlockSize,
                                        &sigWords)) == 0) {
        ind2 = putNextMadrigalRecord(foutp, &cedarRecordp,
                                     madrigalBlockSize,
                                     &blockIndex,
                                     &blockp,
                                     &prevRec,
                                     &thisRec);
        nrec++;
    }
    if (debug) 
        (void) printf("Madrigal to Madrigal - %d records copied\n", nrec);
    if (cedarRecordp != (Int16 *)NULL) {
        if (debug) 
            (void) printf("  Freeing cedarRecordp\n");
        free(cedarRecordp);
        cedarRecordp = (Int16 *)NULL;
    } 
    if (blockp != (Int16 *)NULL) {
        if (debug) 
            (void) printf("  Freeing blockp\n");
        free(blockp);
        blockp = (Int16 *)NULL;
    } 
    (void) fclose(finp);
    (void) fclose(foutp);
    fileType = cedarFileType("file1", madrigalBlockSize, cbfBlockSize);
    if (debug) 
        (void) printf("file1 - fileType = %d\n", fileType);


    /*** Madrigal to Memory ***/
    finp = fopen("file1", "r");
    sigWords = 0;
    mempos1 = 0;
    /*
    (void) fseek(finp, 0L, SEEK_END);
    fileSize = ftell(finp);
    if (debug) 
        printf("fileSize = %d\n", fileSize);
    (void) fseek(finp, 0L, SEEK_SET);
    cedarFile1p = malloc(fileSize+2);
    */
    fileSize1 = 0;
    nrec = 0;
    while((ind1 = getNextMadrigalRecord(finp, &cedarRecordp,
                                        madrigalBlockSize,
                                        &sigWords)) == 0) {
        ind2 = putMemNextCedarUnblockedRecord(&cedarFile1p, &cedarRecordp,
                                              &fileSize1, &mempos1);
        nrec++;
    }
    if (debug) 
        (void) printf("Madrigal to Memory - %d records copied\n", nrec);
    if (cedarRecordp != (Int16 *)NULL) {
        if (debug) 
            (void) printf("  Freeing cedarRecordp\n");
        free(cedarRecordp);
        cedarRecordp = (Int16 *)NULL;
    } 
    (void) fclose(finp);


    /*** Memory to Memory ***/
    if (debug)
        timestart();
    for (i=0; i<100; i++) {
        if (cedarFile2p != (Int16 *)NULL) free(cedarFile2p);
        mempos1 = 0;
        mempos2 = 0;
        /*
        cedarFile2p = malloc(fileSize+2);
        */
        fileSize2 = 0;
        nrec = 0;
        while((ind1 = getMemNextCedarUnblockedRecord(&cedarFile1p, 
                                                     &cedarRecordp,
                                                     &mempos1,
                                                     &recordpInMem)) == 0) {
            ind2 = putMemNextCedarUnblockedRecord(&cedarFile2p, &cedarRecordp,
                                                  &fileSize2,
                                                  &mempos2);
            nrec++;
        }
        /* (void) printf("Memory to Memory - %d records copied\n", nrec); */
        if (cedarRecordp != (Int16 *)NULL) {
            if (debug) 
                (void) printf("  Freeing cedarRecordp\n");
            free(cedarRecordp);
            cedarRecordp = (Int16 *)NULL;
        }
    }
    if (debug) 
        timestop("100 memory to memory copies");


    /*** Memory to Madrigal ***/
    foutp = fopen("file2", "w");
    blockIndex = 0;
    blockp = (Int16 *)NULL;
    prevRec = 0;
    thisRec = 0;
    mempos2 = 0;
    nrec = 0;
    while((ind1 = getMemNextCedarUnblockedRecord(&cedarFile2p, &cedarRecordp,
                                                 &mempos2, &recordpInMem)) == 0) {
        ind2 = putNextMadrigalRecord(foutp, &cedarRecordp,
                                     madrigalBlockSize,
                                     &blockIndex,
                                     &blockp,
                                     &prevRec,
                                     &thisRec);
        nrec++;
    }
    if (cedarFile2p != (Int16 *)NULL) free(cedarFile2p);
    if (debug) 
        (void) printf("Memory to Madrigal- %d records copied\n", nrec);
    if (cedarRecordp != (Int16 *)NULL) {
        if (debug) 
            (void) printf("  Freeing cedarRecordp\n");
        free(cedarRecordp);
        cedarRecordp = (Int16 *)NULL;
    } 
    if (blockp != (Int16 *)NULL) {
        if (debug) 
          (void) printf("  Freeing blockp\n");
        free(blockp);
        blockp = (Int16 *)NULL;
    } 
    (void) fclose(foutp);
    fileType = cedarFileType("file2", madrigalBlockSize, cbfBlockSize);
    if (debug) {
        (void) printf("file2 - fileType = %d\n", fileType);
        (void) printf("Diff %s and file2 - They could be different\n", madfile);
        sprintf(command, "diff %s file2", madfile);
        printf("Executing %s\n", command);
        system(command);
        (void) printf("Diff file1 and file2 - They should be the same\n");
        sprintf(command, "diff file1 file2");
        printf("Executing %s\n", command);
        system(command);
        printf("If there was no error message after the second diff command, the test ran successfully\n");
    }

    free(cedarRecordp);
    if (cedarFile1p != NULL)
        free(cedarFile1p);

    return(0);
}

struct tms {
    time_t tms_utime;
    time_t tms_stime;
    time_t tms_cutime;
    time_t tms_cstime;
};
long times();
#define TICKS 60.
static struct tms tbuf1;
static long real1;

void timestart()    /* start timer */
{
    real1 = times(&tbuf1);
}

void timestop(char *msg)    /* Stop timer */
{
    struct tms tbuf2;
    long real2;

    real2 = times(&tbuf2);

    fprintf(stderr, "%s: real %.2f; user %.2f; sys %.2f\n",msg,
      (real2-real1)/TICKS,
      (tbuf2.tms_utime-tbuf1.tms_utime)/TICKS,
      (tbuf2.tms_stime-tbuf1.tms_stime)/TICKS);
}