---

CTBmath.cxx

Go to the documentation of this file.

/*----------------------------------------------------------------------------*/
/* C Tool Box: Designed and implemented by:                                   */
/*    Walter F.J. Mueller   Gesellschaft fuer Schwerionenforschung (GSI)      */
/*                          Planckstrasse 1, D-64291 Darmstadt, Germany       */
/*                  Email:  W.F.J.Mueller@gsi.de                              */
/*                  WWW:    http://www-kp3.gsi.de/www/kp3/people/mueller.html */
/*------------------------------------------+---------------------------------*/
 
#include <errno.h>
#include <math.h>

#include "CTB.hxx"
#include "CTBnum.hxx"
#include "CTBexceptionMath.hxx"

//------------------------------------------+-----------------------------------
//------------------------------------------+-----------------------------------

double CTBsqrt(double x)
{
  double f;
  
  errno = 0;
  f = sqrt(x);
  if (errno == EDOM) throw CTBexceptionMathDomain(x,"CTBsqrt()");

  return f;
}
 
//------------------------------------------+-----------------------------------

float CTBsqrt(float x)
{
  float f;
  
  errno = 0;
  f = sqrtf(x);
  if (errno == EDOM) throw CTBexceptionMathDomain(x,"CTBsqrt()");

  return f;
}
 
//------------------------------------------+-----------------------------------

double CTBexp(double x)
{
  double f;
  
  errno = 0;
  f = exp(x);
  if (errno == ERANGE) throw CTBexceptionMathRange(x,"CTBexp()");

  return f;
}
 
//------------------------------------------+-----------------------------------

float CTBexp(float x)
{
  float f;
  
  errno = 0;
  f = expf(x);
  if (errno == EDOM) throw CTBexceptionMathRange(x,"CTBexp()");

  return f;
}
 
//------------------------------------------+-----------------------------------

double CTBlog(double x)
{
  double f;
  
  errno = 0;
  f = log(x);
  if (errno == EDOM)   throw CTBexceptionMathDomain(x,"CTBlog()");
  if (errno == ERANGE) throw CTBexceptionMathRange(x,"CTBlog()");

  return f;
}
 
//------------------------------------------+-----------------------------------

float CTBlog(float x)
{
  float f;
  
  errno = 0;
  f = logf(x);
  if (errno == EDOM)   throw CTBexceptionMathDomain(x,"CTBlog()");
  if (errno == ERANGE) throw CTBexceptionMathRange(x,"CTBlog()");

  return f;
}
 
//------------------------------------------+-----------------------------------

double CTBlog10(double x)
{
  double f;
  
  errno = 0;
  f = log10(x);
  if (errno == EDOM)   throw CTBexceptionMathDomain(x,"CTBlog10()");
  if (errno == ERANGE) throw CTBexceptionMathRange(x,"CTBlog10()");

  return f;
}
 
//------------------------------------------+-----------------------------------

float CTBlog10(float x)
{
  float f;
  
  errno = 0;
  f = log10f(x);
  if (errno == EDOM)   throw CTBexceptionMathDomain(x,"CTBlog10()");
  if (errno == ERANGE) throw CTBexceptionMathRange(x,"CTBlog10()");

  return f;
}
 
//------------------------------------------+-----------------------------------

double CTBasin(double x)
{
  double f;
  
  errno = 0;
  f = asin(x);
  if (errno == EDOM)   throw CTBexceptionMathDomain(x,"CTBasin()");

  return f;
}
 
//------------------------------------------+-----------------------------------

float CTBasin(float x)
{
  float f;
  
  errno = 0;
  f = asinf(x);
  if (errno == EDOM)   throw CTBexceptionMathDomain(x,"CTBasin()");

  return f;
}
 
//------------------------------------------+-----------------------------------

double CTBacos(double x)
{
  double f;
  
  errno = 0;
  f = acos(x);
  if (errno == EDOM)   throw CTBexceptionMathDomain(x,"CTBacos()");

  return f;
}
 
//------------------------------------------+-----------------------------------

float CTBacos(float x)
{
  float f;
  
  errno = 0;
  f = acosf(x);
  if (errno == EDOM)   throw CTBexceptionMathDomain(x,"CTBacos()");

  return f;
}
 
//------------------------------------------+-----------------------------------

double CTBpow(double x, int i_y)
{
  double f      = 1.;
  int    i_nmul = CTBabs(i_y);

  for (int i = 0; i < i_nmul; i++) f *= x;
  if (i_y < 0) f = 1./f;

  return f;
}
 
//------------------------------------------+-----------------------------------

float CTBpow(float x, int i_y)
{
  float f      = 1.;
  int   i_nmul = CTBabs(i_y);

  for (int i = 0; i < i_nmul; i++) f *= x;
  if (i_y < 0) f = 1./f;

  return f;
}
 
//------------------------------------------+-----------------------------------

CTBint CTBilog2(CTBint ix)
{
  CTBint ir = 0;
  if (ix <= 0) return -1;
  while (ix > 1) {
    if (ix & 1) return -1;
    ix >>= 1;
    ir  += 1;
  }
  return ir;
}

//------------------------------------------+-----------------------------------

double CTBpythag(double x, double y)
{
  double hypo;

  if (x == 0.)      hypo = fabs(y);
  else if (y == 0.) hypo = fabs(x);

  else {
    double absx = fabs(x);
    double absy = fabs(y);

    if (absx > absy) {
      double fac = absy / absx;
      hypo = absx * sqrt(1. + fac*fac);
    } else {
      double fac = absx / absy;
      hypo = absy * sqrt(1. + fac*fac);
    }
  }

  return hypo;
}
 
//------------------------------------------+-----------------------------------

float CTBpythag(float x, float y)
{
  float hypo;

  if (x == 0.)      hypo = fabs(y);
  else if (y == 0.) hypo = fabs(x);

  else {
    float absx = fabs(x);
    float absy = fabs(y);

    if (absx > absy) {
      float fac = absy / absx;
      hypo = absx * sqrt(1. + fac*fac);
    } else {
      float fac = absx / absy;
      hypo = absy * sqrt(1. + fac*fac);
    }
  }

  return hypo;
}
 
//------------------------------------------+-----------------------------------

---