floatutil.c

Go to the documentation of this file.

/*****************************************************************************/
 
/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <string.h>
#include <ctype.h>
#include <unistd.h>
#include <math.h>
#include "tpcextensions.h"
/*****************************************************************************/
 
/*****************************************************************************/
int floatMatch(
  const float v1,
  const float v2,
  const float lim
) {
  if(isnan(v1) && isnan(v2)) return 1;
  if(isnan(v1) || isnan(v2)) return 0;
  if(v1==v2) return 1;
  if(isnan(lim) || lim<0.0) return 0;
  if(fabsf(v1-v2)<=lim) return 1;
  return 0;
}
/*****************************************************************************/
 
/*****************************************************************************/
int floatMatchRel(
  const float v1,
  const float v2,
  const float lim
) {
  if(isnan(v1) && isnan(v2)) return 1;
  if(isnan(v1) || isnan(v2)) return 0;
  if(v1==v2) return 1;
  if(isnan(lim)) return 0;
  float mean;
  mean=0.5*(v1+v2); if(!isnormal(mean)) return 0;
  if(fabsf((v1-v2)/mean)<=lim) return 1;
  return 0;
}
/*****************************************************************************/
 
/*****************************************************************************/
float floatMachEps()
{
  float macheps=1.0;
  do {macheps/=2.0;} while((1.0+macheps/2.0)!=1.0);
  return(macheps);
}
/*****************************************************************************/
 
/*****************************************************************************/
void floatCopy(
  float *t,
  float *s,
  const unsigned int n
) {
  unsigned int i;
  if(t==NULL || s==NULL || n<1) return;
  for(i=0; i<n; i++) t[i]=s[i];
}
/*****************************************************************************/
 
/*****************************************************************************/
unsigned int floatMaxIndex(
  float *a,
  const unsigned int n
) {
  if(a==NULL || n<1) return(0);
  unsigned int i, mi=0;
  float mv=nanf("");
  for(i=0; i<n; i++) if(isnan(mv) || a[i]>mv) {mv=a[i]; mi=i;}
  return(mi);
}
/*****************************************************************************/
 
/*****************************************************************************/
float floatSum(
  float *a,
  const unsigned int n
) {
  float s=0.0;
  if(a==NULL || n<1) return(s);
  for(unsigned int i=0; i<n; i++) if(!isnan(a[i])) s+=a[i];
  return(s);
}
/*****************************************************************************/
 
/*****************************************************************************/
float floatMean(
  float *a,
  const unsigned int n
) {
  if(a==NULL || n<1) return(nanf(""));
  float s=0.0;
  unsigned int i, ci=0;
  for(i=0; i<n; i++) if(!isnan(a[i])) {ci++; s+=a[i];}
  if(ci<1) return(nanf(""));
  return(s/(float)ci);
}
/*****************************************************************************/
 
/*****************************************************************************/
int floatGetWithUnit(
  const char *s,
  float *v,
  int *u
) {
  if(v!=NULL) *v=nanf(""); 
  if(u!=NULL) *u=UNIT_UNKNOWN; 
  if(s==NULL) return 1;
  int n=strTokenNr(s, " \t"); if(n==0 || n>2) return 1;
  if(n==1) { // no unit
    double dv;
    if(atofCheck(s, &dv)) return 1;
    *v=(float)dv; return 0;
  }
  /* separate number and unit */
  char buf[64];
  n=strTokenNCpy(s, " \t", 1, buf, 64);
  double dv;
  if(atofCheck(buf, &dv)) return 1;
  *v=(float)dv;
  if(u==NULL) return 0;
  n=strTokenNCpy(s, " \t", 2, buf, 64);
  *u=unitIdentify(buf);
  return 0;
}
/*****************************************************************************/
 
/*****************************************************************************/
int floatSpanPositives(
  float *a,
  const int n
) {
  if(a==NULL || n<1) return(0);
  int i=0;
  for(i=0; i<n; i++) if(!(a[i]>0.0)) break;
  return(i);
}
/*****************************************************************************/
 
/*****************************************************************************/
int floatCSpanPositives(
  float *a,
  const int n
) {
  if(a==NULL || n<1) return(0);
  int i=0;
  for(i=0; i<n; i++) if(a[i]>0.0) break;
  return(i);
}
/*****************************************************************************/
 
/*****************************************************************************/
unsigned int floatNonzeroes(
  float *a,
  const unsigned int n
) {
  if(a==NULL || n<1) return(0);
  unsigned int i, m=0;
  for(i=0; i<n; i++) if(fabsf(a[i])>0.0) m++;
  return(m);
}
/*****************************************************************************/
 
/*****************************************************************************/