resift.c File Reference

Utility functions for working with RES and IFT struct. More...

#include "libtpccurveio.h"

Go to the source code of this file.

Functions
int	res2ift (RES *res, IFT *ift, int verbose)

Detailed Description

Utility functions for working with RES and IFT struct.

Author

Vesa Oikonen

Definition in file resift.c.

Function Documentation

res2ift()

int res2ift	(	RES *	res,
		IFT *	ift,
		int	verbose )

Copy results in RES structure into IFT structure.

Returns

Returns 0 when successful.

See also

fit2Result, dftToResult

Parameters

res	Pointer to RES structure
ift	Pointer to initiated IFT structure
verbose	Verbose level; if zero, then nothing is printed into stdout or stderr

Definition at line 14 of file resift.c.

  {
  int ri, pi, ret;
  char tmp[1024], tmp2[1024], tmp3[1024];
  struct tm st;
 
  if(verbose>0) printf("res2ift()\n");
 
  /* Check the input */
  if(res==NULL || ift==NULL) {
    if(verbose>0) fprintf(stderr, "Error: invalid input\n");
    return(1);
  }
 
  /* Delete any previous content */
  iftEmpty(ift);
  
  /* Set header info */
  if(strlen(res->program)>0) {
    ret=iftPut(ift, "program", res->program, NULL, 0);
    if(ret!=0) {
      if(verbose>0) fprintf(stderr, "Error: cannot set IFT content\n");
      return(3);
    }
  }
  if(gmtime_r(&res->time, &st)!=NULL) {
    strftime(tmp, 256, "%Y-%m-%d %H:%M:%S", &st);
    iftPut(ift, "date", tmp, NULL, 0);
  }
  if(res->studynr[0]) iftPut(ift, "studynr", res->studynr, NULL, 0);
  if(res->datafile[0]) iftPut(ift, "datafile", res->datafile, NULL, 0);
  if(res->plasmafile[0]) iftPut(ift, "plasmafile", res->plasmafile, NULL, 0);
  if(res->plasmafile2[0]) iftPut(ift, "plasmafile2", res->plasmafile2, NULL, 0);
  if(res->bloodfile[0]) iftPut(ift, "bloodfile", res->bloodfile, NULL, 0);
  if(res->reffile[0]) iftPut(ift, "reffile", res->reffile, NULL, 0);
  if(res->refroi[0]) iftPut(ift, "refroi", res->refroi, NULL, 0);
 
  if(res->datarange[0]) iftPut(ift, "datarange", res->datarange, NULL, 0);
  if(res->datanr>0) {
    snprintf(tmp, 1024, "%d", res->datanr); iftPut(ift, "datanr", tmp, NULL, 0);}
  if(res->fitmethod[0]) iftPut(ift, "fitmethod", res->fitmethod, NULL, 0);
 
  if(res->density>0) {
    snprintf(tmp, 1024, "%g", res->density); iftPut(ift, "density", tmp, NULL, 0);}
  if(res->lc>0) {
    snprintf(tmp, 1024, "%g", res->lc); iftPut(ift, "lc", tmp, NULL, 0);}
  if(res->concentration>0) {
    snprintf(tmp, 1024, "%g", res->concentration);
    iftPut(ift, "concentration", tmp, NULL, 0);
  }
  if(res->beta>0) {
    snprintf(tmp, 1024, "%g", res->beta); iftPut(ift, "beta", tmp, NULL, 0);}
  if(res->Vb>0) {
    snprintf(tmp, 1024, "%g", res->Vb); iftPut(ift, "Vb", tmp, NULL, 0);}
  if(res->fA>0) {
    snprintf(tmp, 1024, "%g", res->fA); iftPut(ift, "fA", tmp, NULL, 0);}
  if(res->E>0) {
    snprintf(tmp, 1024, "%g", res->E); iftPut(ift, "E", tmp, NULL, 0);}
  if(res->isweight>0) strcpy(tmp, "yes");
  else if(res->isweight==0) strcpy(tmp, "no");
  else strcpy(tmp, "unknown");
  iftPut(ift, "weighting", tmp, NULL, 0);
 
  /* Set each ROI and parameter value */
  for(ri=0; ri<res->voiNr; ri++) {
    if(res->voiNr>1) {
      strcpy(tmp, res->voi[ri].name); rnameRmDots(tmp, NULL);
      strcat(tmp, "_");
    } else {
      strcpy(tmp, "");
    }
    for(pi=0; pi<res->parNr; pi++) {
      snprintf(tmp2, 1024, "%s%s", tmp, res->parname[pi]);
      if(isnan(res->voi[ri].parameter[pi])) strcpy(tmp3, "");
      else sprintf(tmp3, "%g", res->voi[ri].parameter[pi]);
      if(strlen(res->parunit[pi])) {
        strcat(tmp3, " "); strcat(tmp3, res->parunit[pi]);}
      iftPut(ift, tmp2, tmp3, NULL, 0);
      if(!isnan(res->voi[ri].sd[pi])) {
        snprintf(tmp2, 1024, "%s%s_%s", tmp, res->parname[pi], "SD");
        snprintf(tmp3, 1024, "%g", res->voi[ri].sd[pi]);
        if(strlen(res->parunit[pi])) {
          strcat(tmp3, " "); strcat(tmp3, res->parunit[pi]);}
        iftPut(ift, tmp2, tmp3, NULL, 0);
      }
      if(!isnan(res->voi[ri].cl1[pi])) {
        snprintf(tmp2, 1024, "%s%s_%s", tmp, res->parname[pi], "CL1");
        snprintf(tmp3, 1024, "%g", res->voi[ri].cl1[pi]);
        if(strlen(res->parunit[pi])) {
          strcat(tmp3, " "); strcat(tmp3, res->parunit[pi]);}
        iftPut(ift, tmp2, tmp3, NULL, 0);
      }
      if(!isnan(res->voi[ri].cl2[pi])) {
        snprintf(tmp2, 1024, "%s%s_%s", tmp, res->parname[pi], "CL2");
        snprintf(tmp3, 1024, "%g", res->voi[ri].cl2[pi]);
        if(strlen(res->parunit[pi])) {
          strcat(tmp3, " "); strcat(tmp3, res->parunit[pi]);}
        iftPut(ift, tmp2, tmp3, NULL, 0);
      }
    }
  }
 
  return(0);
}