##// END OF EJS Templates
DATABASES » madrigal
RSS
Review updata app
Review updata app
jespinoza - - Load All Authors

File last commit:

r0:b84e1135c2c4
r8:9c1b6919a8a9
Show More
testCedarIOBoundaries.c
323 lines | 9.0 KiB | text/x-c | CLexer
/ source / madc / madrec / testCedarIOBoundaries.c
History | Source | Raw |Copy content |Copy permalink
jespinoza
Initial
 r0 /* $Id: testCedarIOBoundaries.c 3304 2011-01-17 15:25:59Z brideout $ */
/*
modification history
--------------------
00a,21Feb00 original
*/
/*
USAGE: testCedarIOBoundries madrigalFileName numberOfMadrigalCopies numberOfCbfCopies
e.g.: testCedarIOBoundries /opt/madrigal/experiments/1998/mlh/23mar98/mil980323g.003 3 3
The Madrigal and CBF formats are blocked and each block contains
control words in addition to CEDAR records. In practice, the block
lengths are fixed, but the read and write routines in cedarIO.c allow
an arbitrary block size (beyond the minimum possible size). This
program test the blocking by repeatedly doing two transformations:
Madrigal (standard block size) -> Madrigal(random block size) -> Madrigal(standard block size)
Madrigal (standard block size) -> CBF (random block size) -> Madrigal (standard block size)
The first and third files, file1 and file2, are then compared.
Relatively small random block sizes are used to enhance the liklihood
of encountering edge effects.
*/
#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#include <cedarIO.h>
int
main (argc, argv)
int argc;
char *argv[];
{
int i=0, nwords=0, ind1=0, ind2=0, nrec=0, nmad=0, ncbf=0;
char infile[100];
Int16 w1=0, w2=0;
FILE *finp, *foutp;
Int16 *cedarRecordp=(Int16 *)NULL;
/* madrigal state variables */
int madrigalBlockSize=6720, madrigalBlockSizex=0;
int sigWords=0, blockIndex=0, prevRec=0, thisRec=0;
Int16 *blockp=(Int16 *)NULL;
/* Cbf state parameters */
int cbfBlockSize=4096, cbfBlockSizex=0;
int forceCosRead=0, initPos8=0, initFwi=0, initPos=0, initLCosBlock=0,
lCosBlock=0, gpos=0, fwi=0;
Int16 *cosRecordp=(Int16 *) NULL;
int lbuf=2; /* Are two extra words in Cos block */
int ppos=1; /* First word of Cos block is it's length */
int blockNumber=0;
long lastControlWord=-1;
int previousFileIndex=-1, previousRecordIndex=-1;
if (argc == 4) {
strcpy(infile, argv[1]);
nmad = atoi(argv[2]);
ncbf = atoi(argv[3]);
} else {
(void) printf("Usage: testCedarIOBoundries madrigalFileName numberOfMadrigalCopies numberOfCbfCopies\n");
return(1);
}
/* Madrigal to Madrigal */
finp = fopen(infile, "r");
foutp = fopen("file1", "w");
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);
if (ind1 != 0 || ind2 != 0)
(void) printf("Error: ind1,ind2 = %d %d\n", ind1,ind2);
nrec++;
}
/* (void) printf("Initialization - %d records copied\n", nrec); */
if (cedarRecordp != (Int16 *)NULL) {
free(cedarRecordp);
cedarRecordp = (Int16 *)NULL;
}
if (blockp != (Int16 *)NULL) {
free(blockp);
blockp = (Int16 *)NULL;
}
(void) fclose(finp);
(void) fclose(foutp);
for (i=0; i<nmad; i++) {
madrigalBlockSizex = 20 + 1000*((float)rand())/RAND_MAX;
/* printf("madrigalBlockSizex = %d\n", madrigalBlockSizex); */
/* Madrigal to Madrigal */
finp = fopen("file1", "r");
foutp = fopen("file2", "w");
blockIndex = 0;
blockp=(Int16 *)NULL;
prevRec = 0;
thisRec = 0;
nrec = 0;
while((ind1 = getNextMadrigalRecord(finp, &cedarRecordp,
madrigalBlockSize,
&sigWords)) == 0) {
ind2 = putNextMadrigalRecord(foutp, &cedarRecordp,
madrigalBlockSizex,
&blockIndex,
&blockp,
&prevRec,
&thisRec);
nrec++;
}
/* (void) printf(" %d records copied\n", nrec); */
if (cedarRecordp != (Int16 *)NULL) {
free(cedarRecordp);
cedarRecordp = (Int16 *)NULL;
}
if (blockp != (Int16 *)NULL) {
free(blockp);
blockp = (Int16 *)NULL;
}
(void) fclose(finp);
(void) fclose(foutp);
/* Madrigal to Madrigal */
finp = fopen("file2", "r");
foutp = fopen("file3", "w");
blockIndex = 0;
blockp=(Int16 *)NULL;
prevRec = 0;
thisRec = 0;
nrec = 0;
while((ind1 = getNextMadrigalRecord(finp, &cedarRecordp,
madrigalBlockSizex,
&sigWords)) == 0) {
ind2 = putNextMadrigalRecord(foutp, &cedarRecordp,
madrigalBlockSize,
&blockIndex,
&blockp,
&prevRec,
&thisRec);
nrec++;
}
/* (void) printf(" %d records copied\n", nrec); */
if (cedarRecordp != (Int16 *)NULL) {
free(cedarRecordp);
cedarRecordp = (Int16 *)NULL;
}
if (blockp != (Int16 *)NULL) {
free(blockp);
blockp = (Int16 *)NULL;
}
(void) fclose(finp);
(void) fclose(foutp);
finp = fopen("file1", "r");
foutp = fopen("file3", "r");
nwords = 0;
for (;;) {
if (fread(&w1, (size_t) 2, (size_t) 1, finp) != 1) break;
if (fread(&w2, (size_t) 2, (size_t) 1, foutp) != 1) break;
nwords++;
if (w1 != w2) {
printf("***> The files are different - %6d %6d %6d\n <***",
i,w1,w2);
return(0);
}
}
printf(" OK - madrigal - blockSize = %4d, %d words compared\n",
madrigalBlockSizex, nwords);
(void) fclose(finp);
(void) fclose(foutp);
}
for (i=0; i<ncbf; i++) {
cbfBlockSizex = 20 + 492*((float)rand())/RAND_MAX;
cbfBlockSizex = 8*cbfBlockSizex ; /* Must be multiple of 8 bytes */
/* printf("cbfBlockSizex = %d\n", cbfBlockSizex); */
/* Madrigal to Cbf */
finp = fopen("file1", "r");
foutp = fopen("file2", "w+");
lbuf = 2;
ppos = 1;
cosRecordp = (Int16 *) NULL;
blockNumber = 0;
previousFileIndex = -1;
previousRecordIndex = -1;
lastControlWord = -1;
nrec = 0;
while((ind1 = getNextMadrigalRecord(finp, &cedarRecordp,
madrigalBlockSize,
&sigWords)) == 0) {
ind2 = putNextCedarCbfRecord(foutp, &cedarRecordp,
cbfBlockSizex,
&lbuf,
&ppos,
&cosRecordp,
&blockNumber,
&previousFileIndex,
&previousRecordIndex,
&lastControlWord);
nrec++;
}
/* (void) printf("Madrigal to Cbf - %d records copied\n", nrec); */
flushCedarCbfRecord(foutp, &cedarRecordp,
cbfBlockSize,
&lbuf,
&ppos,
&cosRecordp,
&blockNumber,
&previousFileIndex,
&previousRecordIndex,
&lastControlWord);
endFileCedarCbfRecord(foutp, &cedarRecordp,
cbfBlockSize,
&previousFileIndex,
&lastControlWord);
endDataCedarCbfRecord(foutp, &cedarRecordp,
cbfBlockSize,
&lastControlWord);
if (cosRecordp != (Int16 *)NULL) {
free(cosRecordp);
cosRecordp = (Int16 *)NULL;
}
(void) fclose(finp);
(void) fclose(foutp);
/* Cbf to Madrigal */
/* N.B.: Only copies first file */
finp = fopen("file2", "r");
foutp = fopen("file3", "w");
forceCosRead = 0;
lCosBlock = 0;
gpos = 0;
fwi = 0;
cosRecordp = (Int16 *) NULL;
blockIndex = 0;
blockp=(Int16 *)NULL;
prevRec = 0;
thisRec = 0;
nrec = 0;
while((ind1 = getNextCedarCbfRecord(finp, &cedarRecordp,
forceCosRead,
cbfBlockSizex,
&initPos8,
&initFwi,
&initPos,
&initLCosBlock,
&lCosBlock,
&gpos,
&fwi,
&cosRecordp)) == 0) {
ind2 = putNextMadrigalRecord(foutp, &cedarRecordp,
madrigalBlockSize,
&blockIndex,
&blockp,
&prevRec,
&thisRec);
nrec++;
}
/* (void) printf("Cbf to Madrigal - %d records copied\n", nrec); */
if (cedarRecordp != (Int16 *)NULL) {
free(cedarRecordp);
cedarRecordp = (Int16 *)NULL;
}
if (blockp != (Int16 *)NULL) {
free(blockp);
blockp = (Int16 *)NULL;
}
if (cosRecordp != (Int16 *)NULL) {
free(cosRecordp);
cosRecordp = (Int16 *)NULL;
}
(void) fclose(finp);
(void) fclose(foutp);
finp = fopen("file1", "r");
foutp = fopen("file3", "r");
nwords = 0;
for (;;) {
if (fread(&w1, (size_t) 2, (size_t) 1, finp) != 1) break;
if (fread(&w2, (size_t) 2, (size_t) 1, foutp) != 1) break;
nwords++;
if (w1 != w2) {
printf("***> The files are different - %6d %6d %6d <***\n",
i,w1,w2);
return(0);
}
}
if (nwords > 0) {
printf(" OK - cbf - blockSize = %4d, %d words compared\n",
cbfBlockSizex, nwords);
} else {
printf(" ***> ERROR - %d words compared <***\n", nwords);
return(0);
}
(void) fclose(finp);
(void) fclose(foutp);
}
return(0);
}