X-Git-Url: https://git.ucc.asn.au/?p=ipdf%2Fcode.git;a=blobdiff_plain;f=src%2Frational.h;h=0fa5d9910f3b7fa8136717bc971c5a7b6f23487a;hp=f9157c766f6b41f0fbdf061f3369183c500f81cc;hb=888817a67a9d840be66b52811b01eb77f10ff3e6;hpb=8b3424a48d2d2e20de1a0e60ff6e1d84b9b5e226 diff --git a/src/rational.h b/src/rational.h index f9157c7..0fa5d99 100644 --- a/src/rational.h +++ b/src/rational.h @@ -23,6 +23,7 @@ template T Tabs(const T & a) template <> Arbint Tabs(const Arbint & a); template <> Gmpint Tabs(const Gmpint & a); + /* Recursive version of GCD template T gcd(const T & a, const T & b) @@ -172,6 +173,15 @@ struct Rational Rational & operator-=(const Rational & r) {this->operator=(*this-r); return *this;} Rational & operator*=(const Rational & r) {this->operator=(*this*r); return *this;} Rational & operator/=(const Rational & r) {this->operator=(*this/r); return *this;} + Rational Sqrt() const + { + return Rational(sqrt(ToDouble())); + } + + int64_t ToInt64() const + { + return (int64_t)ToDouble(); + } double ToDouble() const { @@ -207,6 +217,11 @@ struct Rational }; +template +Rational

Abs(const Rational

& a) +{ + return Rational

(Tabs(a.P), a.Q); +} }