ScalarTraits.hpp

#pragma once

#include <cmath>
#include <limits>


namespace LSLOpt {

template<typename Scalar>
struct scalar_traits {

};

template<>
struct scalar_traits<double> {
    static constexpr double epsilon() {return std::numeric_limits<double>::epsilon();}
    static constexpr double infinity() {return std::numeric_limits<double>::infinity();}
    static constexpr double nan() {return std::numeric_limits<double>::quiet_NaN();}

    static bool is_nan(const double& d) {return std::isnan(d);}
};

}