tacunits.c File Reference

TAC data unit conversions. More...

#include "tpcclibConfig.h"
#include "tpcift.h"
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
#include "tpctac.h"

Go to the source code of this file.

Functions
int	tacXUnitConvert (TAC *tac, const int u, TPCSTATUS *status)
int	tacYUnitConvert (TAC *tac, const int u, TPCSTATUS *status)
int	tacYUnitVolume2Mass (TAC *tac, const double density, TPCSTATUS *status)
int	tacYUnitMass2Volume (TAC *tac, const double density, TPCSTATUS *status)

Detailed Description

TAC data unit conversions.

Definition in file tacunits.c.

Function Documentation

tacXUnitConvert()

int tacXUnitConvert	(	TAC *	tac,
		const int	u,
		TPCSTATUS *	status )

Convert X values (sample times) in TAC struct to the specified unit.

See also

tacYUnitConvert, unitConversionFactor

Returns

Returns TPCERROR_OK (0) when successful.

Parameters

tac	Pointer to TAC struct
u	Enum unit
status	Pointer to status data; enter NULL if not needed

Definition at line 23 of file tacunits.c.

  {
  int verbose=0; if(status!=NULL) verbose=status->verbose;
  if(verbose>0) {printf("%s(%s)\n", __func__, unitName(u)); fflush(stdout);}
  if(tac==NULL || tac->sampleNr<1) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_NO_DATA);
    return TPCERROR_NO_DATA;
  }
 
  /* If both units are unknown, then return suitable error code */
  if(u==UNIT_UNKNOWN && tac->tunit==UNIT_UNKNOWN) {
    if(verbose>1) {printf("  unknown x units\n"); fflush(stdout);}
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_UNKNOWN_UNIT);
    return TPCERROR_UNKNOWN_UNIT;
  }
 
  /* Determine the correction factor */
  double cf;
  cf=unitConversionFactor(tac->tunit, u);
  if(isnan(cf)) {
    if(verbose>1) {printf("  cannot make x conversion factor\n"); fflush(stdout);}
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_INCOMPATIBLE_UNIT);
    return TPCERROR_INCOMPATIBLE_UNIT;
  }
  
  /* Multiply x values */
  if(verbose>1) {printf("  converting x with factor %g\n", cf); fflush(stdout);}
  for(int i=0; i<tac->sampleNr; i++) {
    tac->x[i]*=cf; tac->x1[i]*=cf; tac->x2[i]*=cf;
  }
  /* Set new unit */
  tac->tunit=u;
  
  statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_OK);
  return(TPCERROR_OK);
}

Referenced by tacInput2sim(), tacReadModelingData(), tacReadModelingInput(), tacReadReference(), tacSetWeights(), tacSetX(), and tacWriteSIF().

tacYUnitConvert()

int tacYUnitConvert	(	TAC *	tac,
		const int	u,
		TPCSTATUS *	status )

Convert Y values (concentrations) in TAC struct to the specified unit.

See also

tacXUnitConvert, unitConversionFactor

Returns

Returns TPCERROR_OK (0) when successful.

Parameters

tac	Pointer to TAC struct
u	Enum unit
status	Pointer to status data; enter NULL if not needed

Definition at line 72 of file tacunits.c.

  {
  int verbose=0; if(status!=NULL) verbose=status->verbose;
  if(verbose>0) printf("%s(%s)\n", __func__, unitName(u));
  if(tac==NULL || tac->sampleNr<1 || tac->tacNr<1) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_NO_DATA);
    return TPCERROR_NO_DATA;
  }
  
  /* If both units are unknown, then return suitable error code */
  if(u==UNIT_UNKNOWN && tac->cunit==UNIT_UNKNOWN) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_UNKNOWN_UNIT);
    return TPCERROR_UNKNOWN_UNIT;
  }
 
  /* Determine the correction factor */
  double cf;
  cf=unitConversionFactor(tac->cunit, u);
  if(isnan(cf)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_INCOMPATIBLE_UNIT);
    return TPCERROR_INCOMPATIBLE_UNIT;
  }
  
  /* Multiply y values */
  for(int j=0; j<tac->tacNr; j++) {
    for(int i=0; i<tac->sampleNr; i++) {
      tac->c[j].y[i]*=cf;
    }
  }
  /* Set new unit */
  tac->cunit=u;
 
  statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_OK);
  return(TPCERROR_OK);
}

Referenced by tacReadModelingData(), tacReadModelingInput(), tacReadReference(), and tacSetWeights().

tacYUnitMass2Volume()

int tacYUnitMass2Volume	(	TAC *	tac,
		const double	density,
		TPCSTATUS *	status )

Convert Y values (per mass concentrations) in TAC struct to the corresponding per volume concentrations.

See also

tacYUnitVolume2Mass, tacXUnitConvert, tacYUnitConvert, unitConversionFactor

Returns

Returns TPCERROR_OK (0) when successful.

Parameters

tac	Pointer to TAC struct.
density	Density (g/mL).
status	Pointer to status data; enter NULL if not needed.

Definition at line 169 of file tacunits.c.

  {
  int verbose=0; if(status!=NULL) verbose=status->verbose;
  if(verbose>0) {printf("%s(%g)\n", __func__, density); fflush(stdout);}
  if(tac==NULL || tac->sampleNr<1 || tac->tacNr<1 || !(density>0.0)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_NO_DATA);
    return TPCERROR_NO_DATA;
  }
 
  /* Check that unit contains mass in divider */
  if(!unitDividerHasMass(tac->cunit)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_INCOMPATIBLE_UNIT);
    return TPCERROR_INCOMPATIBLE_UNIT;
  }
 
  /* Try to convert the unit */
  int uu=unitDividerMassVolumeConversion(tac->cunit);
  if(uu==UNIT_UNKNOWN) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_INCOMPATIBLE_UNIT);
    return TPCERROR_INCOMPATIBLE_UNIT;
  }
 
  /* Multiply y values with the density */
  for(int j=0; j<tac->tacNr; j++) {
    for(int i=0; i<tac->sampleNr; i++) {
      tac->c[j].y[i]*=density;
    }
  }
  /* Set new unit */
  tac->cunit=uu;
 
  statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_OK);
  return(TPCERROR_OK);
}

tacYUnitVolume2Mass()

int tacYUnitVolume2Mass	(	TAC *	tac,
		const double	density,
		TPCSTATUS *	status )

Convert Y values (per volume concentrations) in TAC struct to the corresponding per mass concentrations.

See also

tacYUnitMass2Volume, tacXUnitConvert, tacYUnitConvert, unitConversionFactor

Returns

Returns TPCERROR_OK (0) when successful.

Parameters

tac	Pointer to TAC struct.
density	Density (g/mL).
status	Pointer to status data; enter NULL if not needed.

Definition at line 121 of file tacunits.c.

  {
  int verbose=0; if(status!=NULL) verbose=status->verbose;
  if(verbose>0) {printf("%s(%g)\n", __func__, density); fflush(stdout);}
  if(tac==NULL || tac->sampleNr<1 || tac->tacNr<1 || !(density>0.0)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_NO_DATA);
    return TPCERROR_NO_DATA;
  }
 
  /* Check that unit contains volume in divider */
  if(!unitDividerHasVolume(tac->cunit)) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_INCOMPATIBLE_UNIT);
    return TPCERROR_INCOMPATIBLE_UNIT;
  }
 
  /* Try to convert the unit */
  int uu=unitDividerMassVolumeConversion(tac->cunit);
  if(uu==UNIT_UNKNOWN) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_INCOMPATIBLE_UNIT);
    return TPCERROR_INCOMPATIBLE_UNIT;
  }
 
  /* Divide y values with the density */
  for(int j=0; j<tac->tacNr; j++) {
    for(int i=0; i<tac->sampleNr; i++) {
      tac->c[j].y[i]/=density;
    }
  }
  /* Set new unit */
  tac->cunit=uu;
 
  statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_OK);
  return(TPCERROR_OK);
}