distance.h

Go to the documentation of this file.

 
 
#ifndef THEORETICA_DISTANCE_H
#define THEORETICA_DISTANCE_H
 
#include "./vec.h"
#include "./algebra.h"
#include "../core/constants.h"
#include "../core/core_traits.h"
#include "../core/real_analysis.h"
#include "../complex/complex_analysis.h"
 
 
namespace theoretica {
 
    namespace algebra {
 
 
 
 
        template<typename Vector>
        inline real lp_norm(const Vector& v, unsigned int p) {
 
            real sum = 0;
 
            for (unsigned int i = 0; i < v.size(); ++i)
                sum += pow(abs(v[i]), p);
 
            return root(sum, p);
        }
 
 
        template<typename Vector>
        inline real l1_norm(const Vector& v) {
 
            real sum = 0;
 
            for (unsigned int i = 0; i < v.size(); ++i)
                sum += abs(v[i]);
 
            return sum;
        }
 
        template<typename Vector>
        inline real l2_norm(const Vector& v) {
 
            auto sum = pair_inner_product(v[0], v[0]);
 
            for (unsigned int i = 1; i < v.size(); ++i)
                sum += pair_inner_product(v[i], v[i]);
 
            return real(sqrt(sum));
        }
 
 
        template<typename Vector>
        inline real linf_norm(const Vector& v) {
 
            real res = abs(v[0]);
 
            for (unsigned int i = 1; i < v.size(); ++i)
                res = max(res, abs(v[i]));
 
            return res;
        }
 
 
 
 
        template<typename Vector>
        inline real euclidean_distance(const Vector& v1, const Vector& v2) {
 
            if(v1.size() != v2.size()) {
                TH_MATH_ERROR("euclidean_distance", v1.size(), INVALID_ARGUMENT);
                return nan();
            }
 
            return l2_norm(v1 - v2);
        }
 
 
        template<typename Vector>
        inline real distance(const Vector& v1, const Vector& v2) {
            return euclidean_distance(v1, v2);
        }
 
 
        inline real euclidean_distance(real a, real b) {
            return abs(a - b);
        }
        
        
        template<typename T>
        inline real euclidean_distance(complex<T> z1, complex<T> z2) {
            return abs(z1 - z2);
        }
 
 
        inline real distance(real a, real b) {
            return euclidean_distance(a, b);
        }
 
 
        template<typename T>
        inline real distance(complex<T> z1, complex<T> z2) {
            return euclidean_distance(z1, z2);
        }
 
 
        template<typename Vector>
        inline real minkowski_distance(const Vector& v1, const Vector& v2, unsigned int p) {
 
            if(v1.size() != v2.size()) {
                TH_MATH_ERROR("minkowski_distance", v1.size(), INVALID_ARGUMENT);
                return nan();
            }
 
            return lp_norm(v1 - v2, p);
        }
 
 
        inline real minkowski_distance(real a, real b, unsigned int p) {
            return root(pow(abs(b - a), p), p);
        }
        
        
        template<typename T>
        inline real minkowski_distance(complex<T> z1, complex<T> z2, unsigned int p) {
            return root(pow(abs(z2 - z1), p), p);
        }
 
 
        template<typename Vector>
        inline real manhattan_distance(const Vector& v1, const Vector& v2) {
 
            if(v1.size() != v2.size()) {
                TH_MATH_ERROR("manhattan_distance", v1.size(), INVALID_ARGUMENT);
                return nan();
            }
 
            return l1_norm(v1 - v2);
        }
 
 
        template<typename Vector>
        inline real chebyshev_distance(const Vector& v1, const Vector& v2) {
 
            if(v1.size() != v2.size()) {
                TH_MATH_ERROR("chebyshev_distance", v1.size(), INVALID_ARGUMENT);
                return nan();
            }
 
            return linf_norm(v1 - v2);
        }
 
 
        template<typename Vector>
        inline real discrete_distance(
            const Vector& v1, const Vector& v2, real tolerance = MACH_EPSILON) {
 
            if(v1.size() != v2.size()) {
                TH_MATH_ERROR("discrete_distance", v1.size(), INVALID_ARGUMENT);
                return nan();
            }
 
            bool diff = false;
 
            for (size_t i = 0; i < v1.size(); ++i) {
 
                // The vectors differ if a pair of elements differs
                // more than the given tolerance
                if(abs(v1[i] - v2[i]) > tolerance) {
                    diff = true;
                    break;
                }
            }
 
            return diff ? 1 : 0;
        }
 
 
        template<typename Vector>
        inline real canberra_distance(const Vector& v1, const Vector& v2) {
 
            if(v1.size() != v2.size()) {
                TH_MATH_ERROR("canberra_distance", v1.size(), INVALID_ARGUMENT);
                return nan();
            }
 
            real sum = 0;
 
            for (size_t i = 0; i < v1.size(); ++i)
                sum += abs(v1[i] - v2[i]) / (abs(v1[i]) + abs(v2[i]));
 
            return sum;
        }
 
 
        template<typename Vector>
        inline real cosine_distance(const Vector& v1, const Vector& v2) {
 
            if(v1.size() != v2.size()) {
                TH_MATH_ERROR("cosine_distance", v1.size(), INVALID_ARGUMENT);
                return nan();
            }
 
            real product = 0;
            real sum_sqr_x = 0;
            real sum_sqr_y = 0;
 
            for (size_t i = 0; i < v1.size(); ++i) {
                product += v1[i] * v2[i];
                sum_sqr_x += square(v1[i]);
                sum_sqr_y += square(v2[i]);
            }
 
            return product / sqrt(sum_sqr_x * sum_sqr_y);
        }
 
 
        template<typename Vector>
        inline real hamming_distance(
            const Vector& v1, const Vector& v2, real tolerance = MACH_EPSILON) {
 
            if(v1.size() != v2.size()) {
                TH_MATH_ERROR("hamming_distance", v1.size(), INVALID_ARGUMENT);
                return nan();
            }
 
            unsigned int count = 0;
 
            // Count how many elements differ to some tolerance
            for (unsigned int i = 0; i < v1.size(); ++i)
                if(abs(v1[i] - v2[i]) > tolerance)
                    count++;
 
            return count;
        }
    }
}
 
#endif