X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=flower%2Frational.cc;h=a0c0d22410e235726eeee037b760fee1f7dcaef9;hb=afa40986a9697cdf1cdb499f641fc80046ebe465;hp=2d3195514a555f2b0fc10d8da67cd49cb76eb4ff;hpb=c659cb200486c2f908703696a1b2873e78c8160a;p=lilypond.git diff --git a/flower/rational.cc b/flower/rational.cc index 2d3195514a..a0c0d22410 100644 --- a/flower/rational.cc +++ b/flower/rational.cc @@ -1,34 +1,39 @@ /* rational.cc -- implement Rational - + source file of the Flower Library - (c) 1997--2004 Han-Wen Nienhuys + (c) 1997--2006 Han-Wen Nienhuys */ #include "rational.hh" #include #include +using namespace std; -#include "string-convert.hh" +#include "string-convert.hh" #include "libc-extension.hh" Rational::operator double () const { - return (double)sign_ * num_ / den_; + return ((double)sign_) * num_ / den_; +} + +Rational::operator bool () const +{ + return sign_ && num_; } #ifdef STREAM_SUPPORT ostream & operator << (ostream &o, Rational r) { - o << r.string (); + o << r.string (); return o; } #endif - Rational Rational::trunc_rat () const { @@ -53,10 +58,10 @@ Rational::Rational (int n) { sign_ = ::sign (n); num_ = abs (n); - den_= 1; + den_ = 1; } -static +static inline int gcd (int a, int b) { int t; @@ -68,18 +73,10 @@ int gcd (int a, int b) return b; } -#if 0 -static -int lcm (int a, int b) -{ - return abs (a*b / gcd (a,b)); -} -#endif - void Rational::set_infinite (int s) { - sign_ = ::sign (s) * 2; + sign_ = ::sign (s) * 2; } Rational @@ -126,7 +123,7 @@ Rational::normalise () } else { - int g = gcd (num_ , den_); + int g = gcd (num_, den_); num_ /= g; den_ /= g; @@ -149,22 +146,19 @@ Rational::compare (Rational const &r, Rational const &s) return 0; else if (r.sign_ == 0) return 0; - else - { - return r.sign_ * ::sign (int (r.num_ * s.den_) - int (s.num_ * r.den_)); - } + return r.sign_ * ::sign (int (r.num_ * s.den_) - int (s.num_ * r.den_)); } int compare (Rational const &r, Rational const &s) { - return Rational::compare (r, s ); + return Rational::compare (r, s); } Rational & Rational::operator %= (Rational r) { - *this = r.mod_rat (r); + *this = mod_rat (r); return *this; } @@ -174,25 +168,23 @@ Rational::operator += (Rational r) if (is_infinity ()) ; else if (r.is_infinity ()) + *this = r; + else { - *this = r; - } - else - { - int n = sign_ * num_ *r.den_ + r.sign_ * den_ * r.num_; - int d = den_ * r.den_; - sign_ = ::sign (n) * ::sign (d); + int lcm = (den_ / gcd (r.den_, den_)) * r.den_; + int n = sign_ * num_ * (lcm / den_) + r.sign_ * r.num_ * (lcm / r.den_); + int d = lcm; + sign_ = ::sign (n) * ::sign (d); num_ = abs (n); den_ = abs (d); normalise (); } return *this; } - /* copied from libg++ 2.8.0 - */ +*/ Rational::Rational (double x) { if (x != 0.0) @@ -215,7 +207,7 @@ Rational::Rational (double x) num_ = (unsigned int) (mantissa * FACT); den_ = (unsigned int) FACT; - normalise (); + normalise (); if (expt < 0) den_ <<= -expt; else @@ -226,17 +218,16 @@ Rational::Rational (double x) { num_ = 0; den_ = 1; - sign_ =0; + sign_ = 0; normalise (); } } - void Rational::invert () { int r (num_); - num_ = den_; + num_ = den_; den_ = r; } @@ -245,7 +236,7 @@ Rational::operator *= (Rational r) { sign_ *= ::sign (r.sign_); if (r.is_infinity ()) - { + { sign_ = sign () * 2; goto exit_func; } @@ -257,7 +248,7 @@ Rational::operator *= (Rational r) exit_func: return *this; } - + Rational & Rational::operator /= (Rational r) { @@ -271,23 +262,23 @@ Rational::negate () sign_ *= -1; } -Rational& +Rational & Rational::operator -= (Rational r) { r.negate (); return (*this += r); } -String +string Rational::to_string () const { if (is_infinity ()) { - String s (sign_ > 0 ? "" : "-" ); - return String (s + "infinity"); + string s (sign_ > 0 ? "" : "-"); + return string (s + "infinity"); } - String s = ::to_string (num ()); + string s = ::to_string (num ()); if (den () != 1 && num ()) s += "/" + ::to_string (den ()); return s; @@ -296,7 +287,7 @@ Rational::to_string () const int Rational::to_int () const { - return num () / den (); + return (int) num () / den (); } int