struct DataMap *DataMapReadBlock(int fid,int *s) {
char *buf=NULL;
struct DataMap *ptr;
int32 code,sze,*iptr;
int size=0,cnt=0,num=0,st=0;
st=ConvertReadInt(fid,&code);
if (st==-1) return NULL;
st=ConvertReadInt(fid,&sze);
if (st==-1) return NULL;
if (sze==0) return NULL;
size=sze;
buf=malloc(size);
if (buf==NULL) return NULL;
iptr=(int32 *) buf;
iptr[0]=code;
iptr[1]=sze;
cnt=0;
num=size-2*sizeof(int32);
while (cnt<num) {
st=read(fid,buf+2*sizeof(int32)+cnt,num-cnt);
if (st<=0) break;
cnt+=st;
}
if (cnt<num) {
free(buf);
return NULL;
}
ptr=DataMapDecodeBuffer(buf,size);
free(buf);
if (s !=NULL) *s=size+2*sizeof(int32);
return ptr;
}
int OldFitReadInx(struct OldFitFp *ptr,int32 *buffer,int recno) {
int i,status=0;
if (ptr->inxfp !=-1) {
if (recno !=-1) /* jump to this record */
lseek(ptr->inxfp,(recno+1)*ptr->inx_recl,SEEK_SET);
for (i=0;(i<4) && (status==0);i++)
status=ConvertReadInt(ptr->inxfp,&buffer[i]);
return status;
}
return -1;
}
int fit_read_inx(struct fitfp *ptr,int32 *buffer,int recno) {
int i,status=0;
if (ptr->inxfp !=-1) {
if (recno !=-1) /* jump to this record */
lseek(ptr->inxfp,(recno+1)*ptr->inx_recl,SEEK_SET);
for (i=0;(i<4) && (status==0);i++)
status=ConvertReadInt(ptr->inxfp,&buffer[i]);
return status;
} else if (ptr->ibuf !=NULL) {
if (recno !=-1) ptr->iptr=(recno+1)*ptr->inx_recl;
if ((ptr->istat.st_size-ptr->iptr)<(4*sizeof(int32))) return -1;
for (i=0;i<4;i++)
ConvertToInt(ptr->ibuf+ptr->iptr+sizeof(int32)*i,&buffer[i]);
ptr->iptr+=sizeof(int32)*4;
return 0;
}
return -1;
}
struct RawIndex *RawIndexLoad(int fid) {
void *tmp;
double tme;
int32 inx;
int sze=0;
int st;
struct RawIndex *ptr;
ptr=malloc(sizeof(struct RawIndex));
if (ptr==NULL) return NULL;
ptr->tme=malloc(sizeof(double)*FILE_STEP);
ptr->inx=malloc(sizeof(int)*FILE_STEP);
ptr->num=0;
sze=FILE_STEP;
if ((ptr->tme==NULL) || (ptr->inx==NULL)) {
if (ptr->tme !=NULL) free(ptr->tme);
if (ptr->inx !=NULL) free(ptr->inx);
free(ptr);
return NULL;
}
do {
st=ConvertReadDouble(fid,&tme);
if (st !=0) break;
st=ConvertReadInt(fid,&inx);
if (st !=0) break;
ptr->tme[ptr->num]=tme;
ptr->inx[ptr->num]=inx;
ptr->num++;
st=1;
if (ptr->num>=sze) {
sze+=FILE_STEP;
tmp=realloc(ptr->tme,sizeof(double)*sze);
if (tmp==NULL) break;
ptr->tme=(double *) tmp;
tmp=realloc(ptr->inx,sizeof(int)*sze);
if (tmp==NULL) break;
ptr->inx=(int *) tmp;
}
} while (1);
if (st==1) {
free(ptr->tme);
free(ptr->inx);
free(ptr);
return NULL;
}
tmp=realloc(ptr->tme,sizeof(double)*ptr->num);
if (tmp==NULL) {
free(ptr->tme);
free(ptr->inx);
free(ptr);
return NULL;
}
ptr->tme=tmp;
tmp=realloc(ptr->inx,sizeof(int)*ptr->num);
if (tmp==NULL) {
free(ptr->tme);
free(ptr->inx);
free(ptr);
return NULL;
}
ptr->inx=tmp;
return ptr;
}
int CFitRead(struct CFitfp *fptr,struct CFitdata *ptr) {
float val;
int i;
if (fptr->fp !=-1) {
if (ConvertReadDouble(fptr->fp,&ptr->time) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->stid) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->scan) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->cp) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->bmnum) !=0) return -1;
if (ConvertReadFloat(fptr->fp,&ptr->bmazm) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->channel) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->intt.sc) !=0) return -1;
if (ConvertReadInt(fptr->fp,&ptr->intt.us) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->frang) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->rsep) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->rxrise) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->tfreq) !=0) return -1;
if (ConvertReadInt(fptr->fp,&ptr->noise) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->atten) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->nave) !=0) return -1;
if (ConvertReadShort(fptr->fp,&ptr->nrang) !=0) return -1;
} else {
if ((fptr->fptr+sizeof(double)+7*sizeof(int16)+
sizeof(int32)+sizeof(float))>fptr->fsze) return -1;
ConvertToDouble(&fptr->fbuf[fptr->fptr],&ptr->time);
fptr->fptr+=sizeof(double);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->stid);
fptr->fptr+=sizeof(int16);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->scan);
fptr->fptr+=sizeof(int16);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->cp);
fptr->fptr+=sizeof(int16);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->bmnum);
fptr->fptr+=sizeof(int16);
ConvertToFloat(fptr->fbuf+fptr->fptr,&ptr->bmazm);
fptr->fptr+=sizeof(float);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->channel);
fptr->fptr+=sizeof(int16);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->intt.sc);
fptr->fptr+=sizeof(int16);
ConvertToInt(fptr->fbuf+fptr->fptr,&ptr->intt.us);
fptr->fptr+=sizeof(int32);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->frang);
fptr->fptr+=sizeof(int16);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->rsep);
fptr->fptr+=sizeof(int16);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->rxrise);
fptr->fptr+=sizeof(int16);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->tfreq);
fptr->fptr+=sizeof(int16);
ConvertToInt(fptr->fbuf+fptr->fptr,&ptr->noise);
fptr->fptr+=sizeof(int32);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->atten);
fptr->fptr+=sizeof(int16);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->nave);
fptr->fptr+=sizeof(int16);
ConvertToShort(fptr->fbuf+fptr->fptr,&ptr->nrang);
fptr->fptr+=sizeof(int16);
}
fptr->ctime=ptr->time;
if (fptr->fp !=-1) {
if (read(fptr->fp,&ptr->num,1) !=1) return -1;
if (read(fptr->fp,ptr->rng,ptr->num) !=ptr->num) return -1;
} else {
if (fptr->fptr>=fptr->fsze) return -1;
ptr->num=fptr->fbuf[fptr->fptr];
fptr->fptr+=1;
if ((fptr->fptr+ptr->num)>=fptr->fsze) return -1;
memcpy(ptr->rng,&fptr->fbuf[fptr->fptr],ptr->num);
fptr->fptr+=ptr->num;
}
for (i=0;i<ptr->num;i++) {
if (fptr->fp !=-1) {
if (read(fptr->fp,&ptr->data[i].gsct,1) !=1) return -1;
if (ptr->data[i].gsct==EOF) return -1;
if (ConvertReadFloat(fptr->fp,&val) !=0) return -1;
ptr->data[i].p_0=val;
if (ConvertReadFloat(fptr->fp,&val) !=0) return -1;
ptr->data[i].p_0_e=val;
if (ConvertReadFloat(fptr->fp,&val) !=0) return -1;
ptr->data[i].v=val;
if (ConvertReadFloat(fptr->fp,&val) !=0) return -1;
ptr->data[i].p_l=val;
if (ConvertReadFloat(fptr->fp,&val) !=0) return -1;
ptr->data[i].w_l=val;
if (ConvertReadFloat(fptr->fp,&val) !=0) return -1;
ptr->data[i].v_e=val;
if (ConvertReadFloat(fptr->fp,&val) !=0) return -1;
ptr->data[i].p_l_e=val;
if (ConvertReadFloat(fptr->fp,&val) !=0) return -1;
ptr->data[i].w_l_e=val;
} else {
if ((fptr->fptr+1+sizeof(float)*6)>=fptr->fsze) return -1;
ptr->data[i].gsct=fptr->fbuf[fptr->fptr];
fptr->fptr+=1;
ConvertToFloat(fptr->fbuf+fptr->fptr,&val);
fptr->fptr+=sizeof(float);
ptr->data[i].p_0=val;
ConvertToFloat(fptr->fbuf+fptr->fptr,&val);
fptr->fptr+=sizeof(float);
ptr->data[i].p_0_e=val;
ConvertToFloat(fptr->fbuf+fptr->fptr,&val);
fptr->fptr+=sizeof(float);
ptr->data[i].v=val;
ConvertToFloat(fptr->fbuf+fptr->fptr,&val);
fptr->fptr+=sizeof(float);
ptr->data[i].p_l=val;
ConvertToFloat(fptr->fbuf+fptr->fptr,&val);
fptr->fptr+=sizeof(float);
ptr->data[i].w_l=val;
ConvertToFloat(fptr->fbuf+fptr->fptr,&val);
fptr->fptr+=sizeof(float);
ptr->data[i].v_e=val;
ConvertToFloat(fptr->fbuf+fptr->fptr,&val);
fptr->fptr+=sizeof(float);
ptr->data[i].p_l_e=val;
ConvertToFloat(fptr->fbuf+fptr->fptr,&val);
fptr->fptr+=sizeof(float);
ptr->data[i].w_l_e=val;
}
}
return 0;
}
