##// END OF EJS Templates
DATABASES » madrigal
RSS
Update html images routes
Update html images routes
jespinoza - - Load All Authors

File last commit:

r0:b84e1135c2c4
r14:8f437cc4b32c
Show More
testCedarMemIO.c
257 lines | 7.4 KiB | text/x-c | CLexer
/ source / madc / madrec / testCedarMemIO.c
History | Source | Raw |Copy content |Copy permalink
jespinoza
Initial
 r0 /* $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);
}