ratio.h

Go to the documentation of this file.

 
 
#ifndef THEORETICA_RATIO_H
#define THEORETICA_RATIO_H
 
#ifndef THEORETICA_NO_PRINT
#include <sstream>
#include <ostream>
#endif
 
#include "real_analysis.h"
 
 
namespace theoretica {
 
    template <typename T1, typename T2 = T1>
    class ratio {
        public:
 
            T1 num;
 
            T2 den;
 
            ratio() : num(T1()), den(T2()) {}
 
            ratio(const T1& n, const T2& d) : num(n), den(d) {}
 
            ~ratio() = default;
 
            inline ratio operator*(const ratio& r) const {
                return ratio(num * r.num, den * r.den);
            }
 
            inline ratio operator/(const ratio& r) const {
                return ratio(num * r.den, den * r.num);
            }
 
            inline ratio operator+(const ratio& r) const {
                return ratio(num * r.den + r.num * den, den * r.den);
            }
 
            inline ratio operator-(const ratio& r) const {
                return ratio(num * r.den - r.num * den, den * r.den);
            }
 
            inline ratio operator*(const T1& a) const {
                return ratio(num * a, den);
            }
 
            inline ratio operator/(const T2& b) const {
                return ratio(num, den * b);
            }
 
            template<typename T>
            inline T eval_as() {
                return static_cast<T>(num) / static_cast<T>(den);
            }
 
 
            inline T2 eval() {
                return static_cast<T2>(num) / static_cast<T2>(den);
            }
 
 
            inline operator T2() {
                return eval();
            }
 
 
            template <
                typename = std::enable_if_t<std::is_integral<T1>::value>,
                typename = std::enable_if_t<std::is_integral<T2>::value>
            >
            inline ratio<T1, T2>& reduce() {
 
                if (num == 0) {
                    den = 1;
                    return *this;
                }
 
                const int g = gcd(num, den);
                num /= g;
                den /= g;
 
                return *this;
            }
 
 
#ifndef THEORETICA_NO_PRINT
 
            inline std::string to_string() const {
 
                std::stringstream res;
                res << num << "/" << den;
                return res.str();
            }
 
 
            inline operator std::string() {
                return to_string();
            }
 
 
            inline friend std::ostream& operator<<(std::ostream& out, const ratio<T1, T2>& obj) {
                return out << obj.to_string();
            }
 
#endif
 
    };
}
 
 
#endif