mymath.hxx

Go to the documentation of this file.

#ifndef CH_MATRIX_CLASSES_MYMATH_HXX
#define CH_MATRIX_CLASSES_MYMATH_HXX

#include <cmath>
#include <cstdlib>

namespace CH_Matrix_Classes {

  
inline double abs(double d)
{
  return std::fabs(d);
}
  
inline int abs(int d)
{
  return std::abs(d);
}

inline long abs(long d)
{
  return std::labs(d);
}
  
inline double max(double a,double b)
{
 return (a>=b)?a:b;
}

inline double min(double a,double b)
{
 return (a<=b)?a:b;
}

inline int max(int a,int b)
{
 return (a>=b)?a:b;
}

inline int min(int a,int b)
{
 return (a<=b)?a:b;
}

inline long max(long a,long b)
{
 return (a>=b)?a:b;
}

inline long min(long a,long b)
{
 return (a<=b)?a:b;
}

inline void swap(double& a,double& b)
{
 double h=a; a=b; b=h;
}

inline void swap(long& a,long& b)
{
 long h=a; a=b; b=h;
}

inline void swap(int& a,int& b)
{
 int h=a; a=b; b=h;
}

inline void swap(bool& a,bool& b)
{
 bool h=a; a=b; b=h;
}

inline int sqr(int a)
{
 return a*a;
}

inline long sqr(long a)
{
 return a*a;
}

inline double sqr(double a)
{
 return a*a;
}

inline double sqrt(int a)
{
  return std::sqrt(double(a));
}

inline double sqrt(long a)
{
  return std::sqrt(double(a));
}

inline int sign(int a)
{
 if (a>0) return 1;
 if (a<0) return -1;
 return 0;
}

inline long sign(long a)
{
 if (a>0) return 1;
 if (a<0) return -1;
 return 0;
}

inline double sign(double a)
{
 return (a>=0)?1.:-1.;
}

inline double sign(double a,double tol)
{
 if (a>tol) return 1.;
 if (a<-tol) return -1.;
 return 0.;
}

//---- functions for f2c translations

inline double d_sign(double a,double b)
{
 return (((b>=0)&&(a>=0))||((b<0)&&(a<0)))?a:-a;
}


}

#endif