sif.c

Go to the documentation of this file.

/*****************************************************************************/
int SIF_TEST;
char siferrmsg[128];
/*****************************************************************************/
#include "libtpcimgio.h"
/*****************************************************************************/
 
/*****************************************************************************/
void sifInit(
  SIF *data
) {
  if(SIF_TEST) printf("sifInit()\n");
  if(data==NULL) return;
  memset(data, 0, sizeof(SIF));
  data->frameNr=data->colNr=0;
}
/*****************************************************************************/
 
/*****************************************************************************/
void sifEmpty(
  SIF *data
) {
  if(SIF_TEST) printf("sifEmpty()\n");
  if(data==NULL) return;
  if(data->frameNr>0) {
    free((char*)(data->x1)); free((char*)(data->x2));
    free((char*)(data->prompts)); free((char*)(data->randoms));
    free((char*)(data->trues)); free((char*)(data->weights));
    data->frameNr=data->colNr=0;
  }
  data->scantime=(time_t)0; data->version=0;
  strcpy(data->studynr, ""); strcpy(data->isotope_name, "");
}
/*****************************************************************************/
 
/*****************************************************************************/
int sifSetmem(
  SIF *data,
  int frameNr
) {
  if(SIF_TEST) printf("sifSetmem()\n");
  if(data==NULL) return(1);
  /* Clear previous data, if necessary */
  if(data->frameNr>0) sifEmpty(data);
  if(frameNr<1) return(0);
  
  /* Allocate memory */
  data->x1=(double*)calloc(frameNr, sizeof(double));
  data->x2=(double*)calloc(frameNr, sizeof(double));
  data->prompts=(double*)calloc(frameNr, sizeof(double));
  data->randoms=(double*)calloc(frameNr, sizeof(double));
  data->trues=(double*)calloc(frameNr, sizeof(double));
  data->weights=(double*)calloc(frameNr, sizeof(double));
  if(data->x1==NULL || data->x2==NULL || data->prompts==NULL ||
     data->randoms==NULL || data->trues==NULL || data->weights==NULL) {
    strcpy(siferrmsg, "out of memory"); return(1);}
  data->frameNr=frameNr;
 
  return(0);
}
/*****************************************************************************/
 
/*****************************************************************************/