Floating point number with eps
(Mylib/Geometry/Float/double_eps.cpp)
Operations
Requirements
Notes
Problems
References
Verified with
Code
#pragma once
#include <cmath>
#include <iostream>
#include <limits>
namespace haar_lib {
template <typename T, const T &eps>
struct double_eps {
using value_type = T;
private:
T value_;
public:
double_eps() : value_(0) {}
double_eps(T value_) : value_(value_) {}
auto &operator=(const double_eps &rhs) {
this->value_ = rhs.value_;
return *this;
}
auto &operator+=(const double_eps &rhs) {
this->value_ += rhs.value_;
return *this;
}
auto &operator-=(const double_eps &rhs) {
this->value_ -= rhs.value_;
return *this;
}
auto &operator*=(const double_eps &rhs) {
this->value_ *= rhs.value_;
return *this;
}
auto &operator/=(const double_eps &rhs) {
this->value_ /= rhs.value_;
return *this;
}
auto operator+(const double_eps &rhs) const { return double_eps(this->value_ + rhs.value_); }
auto operator-(const double_eps &rhs) const { return double_eps(this->value_ - rhs.value_); }
auto operator*(const double_eps &rhs) const { return double_eps(this->value_ * rhs.value_); }
auto operator/(const double_eps &rhs) const { return double_eps(this->value_ / rhs.value_); }
bool operator==(const double_eps &rhs) const { return std::abs(this->value_ - rhs.value_) < eps; }
bool operator!=(const double_eps &rhs) const { return !(*this == rhs); }
bool operator<(const double_eps &rhs) const { return this->value_ - rhs.value_ < -eps; }
bool operator<=(const double_eps &rhs) const { return this->value_ - rhs.value_ < eps; }
bool operator>(const double_eps &rhs) const { return !(*this <= rhs); }
bool operator>=(const double_eps &rhs) const { return !(*this < rhs); }
auto operator-() const { return double_eps(-(this->value_)); }
explicit operator double() const noexcept { return value_; }
explicit operator long double() const noexcept { return value_; }
friend std::ostream &operator<<(std::ostream &s, const double_eps &rhs) {
s << rhs.value_;
return s;
}
friend std::istream &operator>>(std::istream &s, double_eps &rhs) {
s >> rhs.value_;
return s;
}
friend double_eps sin(double_eps x) { return std::sin((T) x); }
friend double_eps cos(double_eps x) { return std::cos((T) x); }
friend double_eps tan(double_eps x) { return std::tan((T) x); }
friend double_eps acos(double_eps x) { return std::acos((T) x); }
friend double_eps atan2(double_eps y, double_eps x) { return std::atan2((T) y, (T) x); }
friend double_eps abs(double_eps x) { return std::abs((T) x); }
friend double_eps sqrt(double_eps x) { return std::sqrt(std::max<T>(0, (T) x)); }
};
} // namespace haar_lib
namespace std {
template <typename T, const T &eps>
class numeric_limits<haar_lib::double_eps<T, eps>> {
public:
static haar_lib::double_eps<T, eps> infinity() { return numeric_limits<T>::infinity(); }
static haar_lib::double_eps<T, eps> min() { return numeric_limits<T>::min(); }
static haar_lib::double_eps<T, eps> max() { return numeric_limits<T>::max(); }
static haar_lib::double_eps<T, eps> lowest() { return numeric_limits<T>::lowest(); }
};
} // namespace std
#line 2 "Mylib/Geometry/Float/double_eps.cpp"
#include <cmath>
#include <iostream>
#include <limits>
namespace haar_lib {
template <typename T, const T &eps>
struct double_eps {
using value_type = T;
private:
T value_;
public:
double_eps() : value_(0) {}
double_eps(T value_) : value_(value_) {}
auto &operator=(const double_eps &rhs) {
this->value_ = rhs.value_;
return *this;
}
auto &operator+=(const double_eps &rhs) {
this->value_ += rhs.value_;
return *this;
}
auto &operator-=(const double_eps &rhs) {
this->value_ -= rhs.value_;
return *this;
}
auto &operator*=(const double_eps &rhs) {
this->value_ *= rhs.value_;
return *this;
}
auto &operator/=(const double_eps &rhs) {
this->value_ /= rhs.value_;
return *this;
}
auto operator+(const double_eps &rhs) const { return double_eps(this->value_ + rhs.value_); }
auto operator-(const double_eps &rhs) const { return double_eps(this->value_ - rhs.value_); }
auto operator*(const double_eps &rhs) const { return double_eps(this->value_ * rhs.value_); }
auto operator/(const double_eps &rhs) const { return double_eps(this->value_ / rhs.value_); }
bool operator==(const double_eps &rhs) const { return std::abs(this->value_ - rhs.value_) < eps; }
bool operator!=(const double_eps &rhs) const { return !(*this == rhs); }
bool operator<(const double_eps &rhs) const { return this->value_ - rhs.value_ < -eps; }
bool operator<=(const double_eps &rhs) const { return this->value_ - rhs.value_ < eps; }
bool operator>(const double_eps &rhs) const { return !(*this <= rhs); }
bool operator>=(const double_eps &rhs) const { return !(*this < rhs); }
auto operator-() const { return double_eps(-(this->value_)); }
explicit operator double() const noexcept { return value_; }
explicit operator long double() const noexcept { return value_; }
friend std::ostream &operator<<(std::ostream &s, const double_eps &rhs) {
s << rhs.value_;
return s;
}
friend std::istream &operator>>(std::istream &s, double_eps &rhs) {
s >> rhs.value_;
return s;
}
friend double_eps sin(double_eps x) { return std::sin((T) x); }
friend double_eps cos(double_eps x) { return std::cos((T) x); }
friend double_eps tan(double_eps x) { return std::tan((T) x); }
friend double_eps acos(double_eps x) { return std::acos((T) x); }
friend double_eps atan2(double_eps y, double_eps x) { return std::atan2((T) y, (T) x); }
friend double_eps abs(double_eps x) { return std::abs((T) x); }
friend double_eps sqrt(double_eps x) { return std::sqrt(std::max<T>(0, (T) x)); }
};
} // namespace haar_lib
namespace std {
template <typename T, const T &eps>
class numeric_limits<haar_lib::double_eps<T, eps>> {
public:
static haar_lib::double_eps<T, eps> infinity() { return numeric_limits<T>::infinity(); }
static haar_lib::double_eps<T, eps> min() { return numeric_limits<T>::min(); }
static haar_lib::double_eps<T, eps> max() { return numeric_limits<T>::max(); }
static haar_lib::double_eps<T, eps> lowest() { return numeric_limits<T>::lowest(); }
};
} // namespace std
Back to top page
test/aoj/2136/main.test.cpp
test/aoj/CGL_5_A/main.test.cpp