tacy.c File Reference

Working with y values (concentrations) in TAC struct. 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	tacYRange (TAC *d, int i, double *ymin, double *ymax, int *smin, int *smax, int *imin, int *imax)
	Get the range of y values (concentrations) in TAC struct.
int	tacYRangeInXRange (TAC *d, int i, const double xmin, const double xmax, double *ymin, double *ymax, int *smin, int *smax, int *imin, int *imax)
	Get the range of y values (concentrations) in TAC struct, including only samples with x (times) inside given range.

Detailed Description

Working with y values (concentrations) in TAC struct.

Definition in file tacy.c.

Function Documentation

tacYRange()

int tacYRange	(	TAC *	d,
		int	i,
		double *	ymin,
		double *	ymax,
		int *	smin,
		int *	smax,
		int *	imin,
		int *	imax )

Get the range of y values (concentrations) in TAC struct.

Data is not modified. Data does not need to be sorted. Data can contain NaNs, but y values omitted if corresponding x is NaN.

See also

tacYRangeInXRange, tacXRange, tacYUnitConvert, tacCopyTacc

Author

Vesa Oikonen

Returns

Returns <>0 in case of failure.

Parameters

d	Pointer to TAC structure.
i	Index of TAC, [0, ..., tacNr-1] ; enter <0 to search for max and/or min from all TACs.
ymin	Pointer to variable for min y value (NULL if not needed).
ymax	Pointer to variable for max y value (NULL if not needed).
smin	Pointer to variable where sample index of min y value will be stored (NULL if not needed).
smax	Pointer to variable where sample index of max y value will be stored (NULL if not needed).
imin	Pointer to variable where tac index of min y value will be stored (NULL if not needed).
imax	Pointer to variable where tac index of max y value will be stored (NULL if not needed).

Definition at line 26 of file tacy.c.

  {
  if(ymin!=NULL) *ymin=nan("");
  if(ymax!=NULL) *ymax=nan("");
  if(smin!=NULL) *smin=-1;
  if(smax!=NULL) *smax=-1;
  if(imin!=NULL) *imin=-1;
  if(imax!=NULL) *imax=-1;
  /* Check the data */
  if(d==NULL || d->sampleNr<1) return(1);
  /* Check the index */
  if(i>=d->tacNr) return(2);
 
  /* Find the min and max time */
  double *xi, *xa;
  if(d->isframe) {xi=d->x1; xa=d->x2;} else xi=xa=d->x;
  double mi, ma;
  mi=ma=nan("");
  for(int j=0; j<d->sampleNr; j++) {
    if(isnan(xi[j]) || isnan(xa[j])) continue;
    for(int k=0; k<d->tacNr; k++) {
      if(i>=0 && i!=k) continue;
      if(isnan(mi) || isnan(ma)) {
        mi=d->c[k].y[j]; ma=d->c[k].y[j]; 
        if(smin!=NULL) *smin=j;
        if(smax!=NULL) *smax=j;
        if(imin!=NULL) *imin=k;
        if(imax!=NULL) *imax=k;
        continue;
      }
      if(d->c[k].y[j]<mi) {
        mi=d->c[k].y[j]; 
        if(smin!=NULL) *smin=j;
        if(imin!=NULL) *imin=k;
      } else if(d->c[k].y[j]>ma) {
        ma=d->c[k].y[j];
        if(smax!=NULL) *smax=j;
        if(imax!=NULL) *imax=k;
      }
    }
  }
  if(ymin!=NULL) *ymin=mi;
  if(ymax!=NULL) *ymax=ma;
  if(isnan(mi) || isnan(ma)) return(2);
  return(0);
}

Referenced by mtgaPlotSVG().

tacYRangeInXRange()

int tacYRangeInXRange	(	TAC *	d,
		int	i,
		const double	xmin,
		const double	xmax,
		double *	ymin,
		double *	ymax,
		int *	smin,
		int *	smax,
		int *	imin,
		int *	imax )

Get the range of y values (concentrations) in TAC struct, including only samples with x (times) inside given range.

Data is not modified. Data does not need to be sorted. Data can contain NaNs, but y values omitted if corresponding x is NaN.

See also

tacYRange, tacXRange, tacYUnitConvert

Author

Vesa Oikonen

Returns

Returns <>0 in case of failure.

Parameters

d	Pointer to TAC structure.
i	Index of TAC, [0, ..., tacNr-1] ; enter <0 to search for max and/or min from all TACs.
xmin	Min x value.
xmax	Max x value.
ymin	Pointer to variable for min y value (NULL if not needed).
ymax	Pointer to variable for max y value (NULL if not needed).
smin	Pointer to variable where sample index of min y value will be stored (NULL if not needed).
smax	Pointer to variable where sample index of max y value will be stored (NULL if not needed).
imin	Pointer to variable where tac index of min y value will be stored (NULL if not needed).
imax	Pointer to variable where tac index of max y value will be stored (NULL if not needed).

Definition at line 99 of file tacy.c.

  {
  if(ymin!=NULL) *ymin=nan("");
  if(ymax!=NULL) *ymax=nan("");
  if(smin!=NULL) *smin=-1;
  if(smax!=NULL) *smax=-1;
  if(imin!=NULL) *imin=-1;
  if(imax!=NULL) *imax=-1;
  /* Check the data */
  if(d==NULL || d->sampleNr<1) return(1);
  /* Check the index */
  if(i>=d->tacNr) return(2);
  /* Check x range */
  if(xmin>xmax) return(3);
 
  /* Find the min and max time */
  double *xi, *xa, x;
  if(d->isframe) {xi=d->x1; xa=d->x2;} else xi=xa=d->x;
  double mi, ma;
  mi=ma=nan("");
  for(int j=0; j<d->sampleNr; j++) {
    if(isnan(xi[j]) || isnan(xa[j])) continue;
    x=0.5*(xi[j]+xa[j]); if(x<xmin || x>xmax) continue;
    for(int k=0; k<d->tacNr; k++) {
      if(i>=0 && i!=k) continue;
      if(isnan(mi) || isnan(ma)) {
        mi=d->c[k].y[j]; ma=d->c[k].y[j]; 
        if(smin!=NULL) *smin=j;
        if(smax!=NULL) *smax=j;
        if(imin!=NULL) *imin=k;
        if(imax!=NULL) *imax=k;
        continue;
      }
      if(d->c[k].y[j]<mi) {
        mi=d->c[k].y[j]; 
        if(smin!=NULL) *smin=j;
        if(imin!=NULL) *imin=k;
      } else if(d->c[k].y[j]>ma) {
        ma=d->c[k].y[j];
        if(smax!=NULL) *smax=j;
        if(imax!=NULL) *imax=k;
      }
    }
  }
  if(ymin!=NULL) *ymin=mi;
  if(ymax!=NULL) *ymax=ma;
  if(isnan(mi) || isnan(ma)) return(2);
  return(0);
}

Referenced by tacPlotFitSVG(), tacPlotHistogramSVG(), and tacPlotLineSVG().