AnimaL

Tutorial

Documentation

---

Documentation

• Overview
• Modules
• Namespace List
• Class Hierarchy
• Alphabetical List
• Compound List
• File List
• Namespace Members
• Compound Members
• File Members
• Examples

---

numerics.h

Go to the documentation of this file.

#ifndef ANIMAL_TRAITS_NUMERICS_H
#define ANIMAL_TRAITS_NUMERICS_H

#include <math.h>
#include <float.h>



namespace animal {

template <class T> inline T pow2(const T &f) { return f * f; }

template <class T> inline T pow3(const T &f) { return f * f * f; }

template <class T> inline T pow4(const T &f) { return f * f * f * f; }

typedef unsigned int Natural;

// -------------------------------------------------------

//template<class T>
//struct Value
//{
//  static inline T zero(){ return 0; }
//};

// -------------------------------------------------------
//  
//  Numerics class.
//  
// -------------------------------------------------------
template <class Number> class Numerics
{

public:
  
    typedef int Boolean;

    typedef double Real;

    
    static inline Real PI()
        { return 3.141592653589793238462643383;}
        
    static inline Real CONV_RAD_TO_DEG()
        { return 1.0/2.0/PI()*360.0; }
                
    static inline Real CONV_DEG_TO_RAD()
        { return 2.0*PI()/360.0; }
    
    static inline Real E()
        { return 2.7182818284590452354; }
    
    

    static Number fpabs(Number k) { return k<0 ? -k : k ; };

    static Real sqroot(Number k) { return sqrt(k); };

    static Real sine(Number k) { return sin(k); };

    static Real asine(Number k) { return ::asinf(k); };

    static Real cosine(Number k) { return cos(k); };

    static Real acosine(Number k) { return acos(k); };

    static Real tan(Number k) { return ::tanf(k); };

    static Real atan(Number k) { return ::atanf(k); };

    static Real e(Number k) { return ::expf(k); };

    static Real ln(Number k) { return ::logf(k); };

    static Real power(Number k1,Number k2) { return ::powf(k1,k2); };

    static Real prev_integer(Number a) { return ::floorf(a); };

    static Real next_integer(Number a) { return ::ceilf(a); };

    static Boolean isfinite(Number k) { return true; };

    static Real numthreshold() { return 0; };
    
  
}; // class Numerics<>


// -------------------------------------------------------
//  
//  Numerics<float> class.
//  
// -------------------------------------------------------

template <> 
class Numerics<float>
{
public:
  
    typedef bool Boolean;

    typedef float Real;

    
    static inline Real PI()
        { return 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998;}
        
    static inline Real CONV_RAD_TO_DEG()
        { return 1.0/2.0/PI()*360.0; }
                
    static inline Real CONV_DEG_TO_RAD()
        { return 2.0*PI()/360.0; }
    
    static inline Real E()
        { return 2.7182818284590452354; }
    
    

    static Real fpabs(Real k) { return ::fabsf(k); }

    static Real sqroot(Real k) { return ::sqrtf(k); }

    static Real sine(Real k) { return ::sinf(k); }

    static Real asine(Real k) { return ::asinf(k); }

    static Real cosine(Real k) { return ::cosf(k); }

    static Real acosine(Real k) { return ::acosf(k); }

    static Real tan(Real k) { return ::tanf(k); }

    static Real atan(Real k) { return ::atanf(k); }

    static Real e(Real k) { return ::expf(k); }

    static Real ln(Real k) { return ::logf(k); }
    
    static Real power(Real k1,Real k2) { return ::powf(k1,k2); }
    
    static Real prev_integer(Real a) { return ::floorf(a); };

    static Real next_integer(Real a) { return ::ceilf(a); };
    
    //static Boolean isfinite(Real k) { return ::finite( static_cast<double>(k) ); }

    //static Real numthreshold() { return FLT_EPSILON; }

    //static Real max() { return FLT_MAX; }

    //static Real min() { return FLT_MIN; }
    
  
}; // class Numerics<float>



// -------------------------------------------------------
//  
//  Numerics<double> class.
//  
// -------------------------------------------------------

template <> class Numerics<double>
{

public:

    typedef bool Boolean;

    typedef double Real;
    
    
    
    static inline Real PI()
        { return 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998;}
        
    static inline Real CONV_RAD_TO_DEG()
        { return 1.0/2.0/PI()*360.0; }
                
    static inline Real CONV_DEG_TO_RAD()
        { return 2.0*PI()/360.0; }
    
    static inline Real E()
        { return 2.7182818284590452354; }
    
    

    static Real fpabs(Real k) { return ::fabs(k); }

    static Real sqroot(Real k) { return ::sqrt(k); }

    static Real sine(Real k) { return ::sin(k); }

    static Real asine(Real k) { return ::asinf(k); }

    static Real cosine(Real k) { return ::cos(k); }

    static Real acosine(Real k) { return ::acos(k); }

    static Real tan(Real k) { return ::tan(k); }

    static Real atan(Real k) { return ::atan(k); }

    static Real e(Real k) { return ::exp(k); }

    static Real ln(Real k) { return ::log(k); }
    
    static Real power(Real k1,Real k2) { return ::pow(k1,k2); }
    
    static Real prev_integer(Real a) { return ::floor(a); };

    static Real next_integer(Real a) { return ::ceil(a); };
    
    //static Boolean isfinite(Real k) { return ::finite(k); }

    //static Real numthreshold() { return DBL_EPSILON; }

    //static Real max() { return DBL_MAX; }

    //static Real min() { return DBL_MIN; }
    //
  
}; // class Numerics<double>


} // namespace animal 



#endif // ANIMAL_TRAITS_NUMERICS_H

---