/*! The file is opened.
\param filename (char *) the filename
\param mode (int) the i/o mode , possible values are O_RDONLY, O_WRONLY,
O_RDWR, O_APPEND, O_TMP, O_CREAT, O_TRUNC - see ccp4_sysdep.h
\return (void *) CMMFile structure */
void *ccp4_cmap_open(const char *filename, int mode)
{
CMMFile *mfile;
CCP4File *cfile;
const size_t stamp_offset = 212U;
if ((cfile = ccp4_file_open(filename, mode)) == NULL) {
ccp4_signal( CCP4_ERRLEVEL(3) | CMAP_ERRNO(CMERR_CantOpenFile),
"ccp4_cmap_open",NULL);
return (NULL); }
ccp4_file_raw_setstamp(cfile, stamp_offset);
/* read or write only */
if (cfile->read) {
if (!is_cmap(cfile) || cfile->length < 1025) {
ccp4_signal( CCP4_ERRLEVEL(3) | CMAP_ERRNO(CMERR_NoHeader),
"ccp4_cmap_open",NULL);
ccp4_file_close(cfile);
return NULL; }
ccp4_file_rarch(cfile);
mfile = init_cmap_read();
mfile->stream = cfile;
mfile->file_name = cfile->name;
parse_mapheader(mfile);
parse_maplabels(mfile);
} else if (cfile->write) {
mfile = init_cmap_write();
mfile->stream = cfile;
mfile->file_name = cfile->name;
write_mapheader(mfile);
write_maplabels(mfile);
} else {
ccp4_signal( CCP4_ERRLEVEL(3) | CMAP_ERRNO(CMERR_CantOpenFile),
"ccp4_cmap_open",NULL);
return (NULL); }
return (mfile);
}
/*! Close the file.
In write mode the header is output, along with the machine
stamp. In read mode the file is just closed.
Write mode supports ways of updating the map statistics (
only active for FLOAT32).
/param mfile (CMMFile *)
/return void */
void ccp4_cmap_close(CMMFile *mfile)
{
int i;
if ( mfile == NULL)
return;
if (ccp4_file_is_write(mfile->stream) ) {
if ( mfile->data_mode == FLOAT32) {
switch (mfile->close_mode) {
case 1:
break;
case 2:
mfile->stats.offset = 0.0f;
case 0:
default:
if (mfile->stats.total != 0) {
mfile->stats.mean /= mfile->stats.total;
mfile->stats.rms /= mfile->stats.total;
mfile->stats.rms -= mfile->stats.mean*mfile->stats.mean;
mfile->stats.rms = (mfile->stats.rms > 0) ? sqrt(mfile->stats.rms) : 0;
mfile->stats.mean += (double) mfile->stats.offset;
}
break;
}
}
write_mapheader(mfile);
write_maplabels(mfile);
ccp4_file_warch(mfile->stream);
}
ccp4_file_close(mfile->stream);
for (i=0 ; i != mfile->labels.number ; i++)
if (mfile->labels.labels[i] != NULL)
free(mfile->labels.labels[i]);
free(mfile);
}
