dftres.c

Go to the documentation of this file.

/*****************************************************************************/
#include "libtpccurveio.h"
/*****************************************************************************/
 
/*****************************************************************************/
int res_allocate_with_dft(
  RES *res,
  DFT *dft
) {
  int ri;
 
  //printf("res_allocate_with_dft()\n"); fflush(stdout);
  // Check the input data
  if(res==NULL || dft==NULL || dft->voiNr<1) return 1;
  // Allocate memory
  if(resSetmem(res, dft->voiNr)!=0) return 2;
  res->voiNr=dft->voiNr;
  // Copy header information
  strcpy(res->studynr, dft->studynr);
  res->Vb=-1.0;
  res->fA=-1.0;
  res->E=-1.0;
  res->time=time(NULL); // Set current time to results
  res->isweight=dft->isweight;
  /* Copy region names, etc */
  for(ri=0; ri<dft->voiNr; ri++) {
    strcpy(res->voi[ri].name, dft->voi[ri].name);
    strcpy(res->voi[ri].voiname, dft->voi[ri].voiname);
    strcpy(res->voi[ri].hemisphere, dft->voi[ri].hemisphere);
    strcpy(res->voi[ri].place, dft->voi[ri].place);
  }
  /* Set data range */
  if(dft->timetype==DFT_TIME_STARTEND)
    sprintf(res->datarange, "%g - %g %s",
      dft->x1[0], dft->x2[dft->frameNr-1], petTunit(dft->timeunit) );
  else
    sprintf(res->datarange, "%g - %g %s",
      dft->x[0], dft->x[dft->frameNr-1], petTunit(dft->timeunit) );
  res->datanr=dft->frameNr;
 
  return 0;
}
/*****************************************************************************/
 
/*****************************************************************************/
int dftToResult(
  DFT *dft,
  RES *res,
  char *status
) {
  int ri, fi;
 
  // Check the input data
  if(status!=NULL) sprintf(status, "program error");
  if(res==NULL || dft==NULL || dft->voiNr<1 || dft->frameNr<1) return 1;
  // Allocate memory and copy most of headers
  if(res_allocate_with_dft(res, dft) != 0) {
    if(status!=NULL) sprintf(status, "cannot setup results data");
    return 2;
  }
  res->parNr=dft->frameNr; if(res->parNr>MAX_RESPARAMS) {
    sprintf(status, "only %d frames can be copied to results", MAX_RESPARAMS);
    res->parNr=MAX_RESPARAMS;
  }
  // Set parameter titles and units
  for(fi=0; fi<res->parNr; fi++) {
    sprintf(res->parname[fi], "%d", fi+1);
    strcpy(res->parunit[fi], dft->unit);
  }
  // Copy regional values
  for(ri=0; ri<dft->voiNr; ri++) for(fi=0; fi<res->parNr; fi++)
    res->voi[ri].parameter[fi]=dft->voi[ri].y[fi];
 
  if(dft->frameNr>MAX_RESPARAMS) return -1;
  /* Set also deprecated parameter name and unit representations, for now */
  resFixParnames(res);
  return 0;
}
/*****************************************************************************/
 
/*****************************************************************************/