norms.hpp

Go to the documentation of this file.

 
/*****************************************************************************
** Ifdefs
*****************************************************************************/
 
#ifndef ECL_MATH_NORMS_HPP_
#define ECL_MATH_NORMS_HPP_
 
/*****************************************************************************
** Includes
*****************************************************************************/
 
#include <cmath>
 
/*****************************************************************************
** Namespaces
*****************************************************************************/
 
namespace ecl {
 
/*****************************************************************************
** Functors
*****************************************************************************/
class EuclideanNorm {
public:
        template <typename T>
        T operator()(const T& x1, const T& x2) {
                return std::sqrt(x1*x1 + x2*x2);
        }
 
        template <typename T>
        T operator()(const T& x1, const T& x2, const T& x3) {
                return std::sqrt(x1*x1 + x2*x2 + x3*x3);
        }
 
        // n dimensional cases using vectors, valarrays and eigen vectors here.
 
        // maybe also a norm(vectorx, vectory) which is sqrt((x0-y0)^2+(x1-y1)^2.....)
};
 
/*****************************************************************************
** Functions
*****************************************************************************/
template <typename T>
T euclidean_norm(const T& x1, const T& x2) {
        return std::sqrt(x1*x1 + x2*x2);
}
 
template <typename T>
T euclidean_norm(const T& x1, const T& x2, const T& x3) {
        return std::sqrt(x1*x1 + x2*x2 + x3*x3);
}
 
} // namespace ecl
 
#endif /* ECL_MATH_NORMS_HPP_ */