AnimaL

Tutorial

Documentation

---

Documentation

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

---

utils.h

Go to the documentation of this file.

#ifndef ANIMAL_ALGO_UTILS
#define ANIMAL_ALGO_UTILS

#include <assert.h>
#include <iostream>
#include <vector>

namespace nr {

template <class T>
inline T 
MAX( const T& a, const T& b ) { return a > b ? a : b; }

template <class T1,class T2>
inline T2 
MAX( const T1& a, const T2& b ) { return a > b ? static_cast<T2>(a) : b; }


template <class T>
inline T 
SQR( const T& a ) { return a * a; }

template <class T>
inline T 
SIGN(const T& a, const T& b) { return ((b) >= 0.0 ? fabs(a) : -fabs(a)); }

inline void nrerror( const char* s ){
    std::cerr << s << std::endl;
    assert(false);
}



//=============================================================================
//=============================================================================
template 
<
    class RealT = double,
    class VectorT = RealT*,
    class MatrixT = RealT**
>
struct Old_Traits
{

    typedef RealT Real;
    typedef VectorT Vector;
    typedef MatrixT Matrix;
    

    static inline void alloc_vector( Vector& v, unsigned int size  )
    {
        v = new Real[ size ];
    }

    static inline void free_vector( Vector& vec )
    {
        delete vec;
    }

    static inline void alloc_matrix( Matrix& mat, unsigned int rows, unsigned int columns  )
    {
        mat = new Real*[rows];
        for( unsigned int i = 0; i<rows; ++i ){
            mat[i] = new Real[columns];
        }
    }

    static inline void free_matrix( Matrix& mat, unsigned int rows, unsigned int  )
    {
        for( unsigned int i = 0; i<rows; ++i ){
            delete mat[i];
        }
        delete mat;
    }

    
};




//=============================================================================
//=============================================================================
template 
<
    class RealT = double,
    class VectorT = std::vector< RealT >,
    class MatrixT = std::vector< VectorT >
>
struct STL_Traits
{

    typedef RealT Real;
    typedef VectorT Vector;
    typedef MatrixT Matrix;

    static void alloc_vector( Vector& v, unsigned int size )
    {
        v = Vector(size);
    }

    static inline void free_vector( Vector& )
    {
    }
    
    static inline void alloc_matrix( Matrix& mat, unsigned int rows, unsigned int columns  )
    {
        mat = Matrix(rows);
        for( unsigned int i = 0; i<rows; ++i ){
            mat[i] = Vector(columns);
        }
    }
    
    static inline void free_matrix( Matrix&, unsigned int, unsigned int  )
    {
    }

};

}
#endif

---