}
DEBUG_FUNCTION void
-debug (sreal &ref)
+debug (const sreal &ref)
{
ref.dump (stderr);
}
DEBUG_FUNCTION void
-debug (sreal *ptr)
+debug (const sreal *ptr)
{
if (ptr)
debug (*ptr);
m_exp += s;
- m_sig += (uint64_t) 1 << (s - 1);
+ m_sig += (int64_t) 1 << (s - 1);
m_sig >>= s;
}
void
sreal::normalize ()
{
- if (m_sig == 0)
+ int64_t s = m_sig < 0 ? -1 : 1;
+ unsigned HOST_WIDE_INT sig = absu_hwi (m_sig);
+
+ if (sig == 0)
{
- m_negative = 0;
m_exp = -SREAL_MAX_EXP;
}
- else if (m_sig < SREAL_MIN_SIG)
+ else if (sig < SREAL_MIN_SIG)
{
do
{
- m_sig <<= 1;
+ sig <<= 1;
m_exp--;
}
- while (m_sig < SREAL_MIN_SIG);
+ while (sig < SREAL_MIN_SIG);
/* Check underflow. */
if (m_exp < -SREAL_MAX_EXP)
{
m_exp = -SREAL_MAX_EXP;
- m_sig = 0;
+ sig = 0;
}
}
- else if (m_sig > SREAL_MAX_SIG)
+ else if (sig > SREAL_MAX_SIG)
{
int last_bit;
do
{
- last_bit = m_sig & 1;
- m_sig >>= 1;
+ last_bit = sig & 1;
+ sig >>= 1;
m_exp++;
}
- while (m_sig > SREAL_MAX_SIG);
+ while (sig > SREAL_MAX_SIG);
/* Round the number. */
- m_sig += last_bit;
- if (m_sig > SREAL_MAX_SIG)
+ sig += last_bit;
+ if (sig > SREAL_MAX_SIG)
{
- m_sig >>= 1;
+ sig >>= 1;
m_exp++;
}
if (m_exp > SREAL_MAX_EXP)
{
m_exp = SREAL_MAX_EXP;
- m_sig = SREAL_MAX_SIG;
+ sig = SREAL_MAX_SIG;
}
}
+
+ m_sig = s * sig;
}
/* Return integer value of *this. */
int64_t
sreal::to_int () const
{
- int64_t sign = m_negative ? -1 : 1;
+ int64_t sign = m_sig < 0 ? -1 : 1;
if (m_exp <= -SREAL_BITS)
return 0;
if (m_exp >= SREAL_PART_BITS)
return sign * INTTYPE_MAXIMUM (int64_t);
if (m_exp > 0)
- return sign * (m_sig << m_exp);
+ return m_sig << m_exp;
if (m_exp < 0)
- return sign * (m_sig >> -m_exp);
- return sign * m_sig;
+ return m_sig >> -m_exp;
+ return m_sig;
}
/* Return *this + other. */
sreal
sreal::operator+ (const sreal &other) const
{
- const sreal *a_p = this, *b_p = &other;
-
- if (a_p->m_negative && !b_p->m_negative)
- std::swap (a_p, b_p);
-
- /* a + -b => a - b. */
- if (!a_p->m_negative && b_p->m_negative)
- {
- sreal tmp = -(*b_p);
- if (*a_p < tmp)
- return signedless_minus (tmp, *a_p, true);
- else
- return signedless_minus (*a_p, tmp, false);
- }
-
- gcc_checking_assert (a_p->m_negative == b_p->m_negative);
-
- sreal r = signedless_plus (*a_p, *b_p, a_p->m_negative);
-
- return r;
-}
-
-sreal
-sreal::signedless_plus (const sreal &a, const sreal &b, bool negative)
-{
- const sreal *bb;
- sreal r, tmp;
int dexp;
- const sreal *a_p = &a;
- const sreal *b_p = &b;
+ sreal tmp, r;
+
+ const sreal *a_p = this, *b_p = &other, *bb;
if (a_p->m_exp < b_p->m_exp)
std::swap (a_p, b_p);
if (dexp > SREAL_BITS)
{
r.m_sig = a_p->m_sig;
- r.m_negative = negative;
return r;
}
r.m_sig = a_p->m_sig + bb->m_sig;
r.normalize ();
-
- r.m_negative = negative;
return r;
}
+
/* Return *this - other. */
sreal
sreal::operator- (const sreal &other) const
-{
- /* -a - b => -a + (-b). */
- if (m_negative && !other.m_negative)
- return signedless_plus (*this, -other, true);
-
- /* a - (-b) => a + b. */
- if (!m_negative && other.m_negative)
- return signedless_plus (*this, -other, false);
-
- gcc_checking_assert (m_negative == other.m_negative);
-
- /* We want to substract a smaller number from bigger
- for nonegative numbers. */
- if (!m_negative && *this < other)
- return signedless_minus (other, *this, true);
-
- /* Example: -2 - (-3) => 3 - 2 */
- if (m_negative && *this > other)
- return signedless_minus (-other, -(*this), false);
-
- sreal r = signedless_minus (*this, other, m_negative);
-
- return r;
-}
-
-sreal
-sreal::signedless_minus (const sreal &a, const sreal &b, bool negative)
{
int dexp;
sreal tmp, r;
const sreal *bb;
- const sreal *a_p = &a;
- const sreal *b_p = &b;
+ const sreal *a_p = this, *b_p = &other;
- dexp = a_p->m_exp - b_p->m_exp;
+ int64_t sign = 1;
+ if (a_p->m_exp < b_p->m_exp)
+ {
+ sign = -1;
+ std::swap (a_p, b_p);
+ }
+ dexp = a_p->m_exp - b_p->m_exp;
r.m_exp = a_p->m_exp;
if (dexp > SREAL_BITS)
{
- r.m_sig = a_p->m_sig;
- r.m_negative = negative;
+ r.m_sig = sign * a_p->m_sig;
return r;
}
if (dexp == 0)
bb = &tmp;
}
- r.m_sig = a_p->m_sig - bb->m_sig;
+ r.m_sig = sign * (a_p->m_sig - bb->m_sig);
r.normalize ();
-
- r.m_negative = negative;
return r;
}
sreal::operator* (const sreal &other) const
{
sreal r;
- if (m_sig < SREAL_MIN_SIG || other.m_sig < SREAL_MIN_SIG)
+ if (std::abs (m_sig) < SREAL_MIN_SIG || std::abs (other.m_sig) < SREAL_MIN_SIG)
{
r.m_sig = 0;
r.m_exp = -SREAL_MAX_EXP;
r.normalize ();
}
- r.m_negative = m_negative ^ other.m_negative;
return r;
}
sreal r;
r.m_sig = (m_sig << SREAL_PART_BITS) / other.m_sig;
r.m_exp = m_exp - other.m_exp - SREAL_PART_BITS;
- r.m_negative = m_negative ^ other.m_negative;
r.normalize ();
return r;
}
#define UINT64_BITS 64
-#define SREAL_MIN_SIG ((uint64_t) 1 << (SREAL_PART_BITS - 1))
-#define SREAL_MAX_SIG (((uint64_t) 1 << SREAL_PART_BITS) - 1)
+#define SREAL_MIN_SIG ((int64_t) 1 << (SREAL_PART_BITS - 2))
+#define SREAL_MAX_SIG (((int64_t) 1 << (SREAL_PART_BITS - 1)) - 1)
#define SREAL_MAX_EXP (INT_MAX / 4)
#define SREAL_BITS SREAL_PART_BITS
{
public:
/* Construct an uninitialized sreal. */
- sreal () : m_sig (-1), m_exp (-1), m_negative (0) {}
+ sreal () : m_sig (-1), m_exp (-1) {}
/* Construct a sreal. */
- sreal (int64_t sig, int exp = 0) : m_exp (exp)
+ sreal (int64_t sig, int exp = 0) : m_sig (sig), m_exp (exp)
{
- m_negative = sig < 0;
-
- if (sig < 0)
- sig = -sig;
-
- m_sig = (uint64_t) sig;
-
normalize ();
}
bool operator< (const sreal &other) const
{
- /* We negate result in case of negative numbers and
- it would return true for equal negative numbers. */
- if (*this == other)
- return false;
-
- if (m_negative != other.m_negative)
- return m_negative > other.m_negative;
-
- bool r = m_exp < other.m_exp
- || (m_exp == other.m_exp && m_sig < other.m_sig);
-
- return m_negative ? !r : r;
+ if (m_exp == other.m_exp)
+ return m_sig < other.m_sig;
+ else
+ {
+ bool negative = m_sig < 0;
+ bool other_negative = other.m_sig < 0;
+
+ if (negative != other_negative)
+ return negative > other_negative;
+
+ bool r = m_exp < other.m_exp;
+ return negative ? !r : r;
+ }
}
bool operator== (const sreal &other) const
{
- return m_exp == other.m_exp && m_sig == other.m_sig
- && m_negative == other.m_negative;
+ return m_exp == other.m_exp && m_sig == other.m_sig;
}
sreal operator- () const
{
- if (m_sig == 0)
- return *this;
-
sreal tmp = *this;
- tmp.m_negative = !tmp.m_negative;
+ tmp.m_sig *= -1;
return tmp;
}
private:
void normalize ();
void shift_right (int amount);
-
static sreal signedless_plus (const sreal &a, const sreal &b, bool negative);
static sreal signedless_minus (const sreal &a, const sreal &b, bool negative);
- uint64_t m_sig; /* Significant. */
+ int64_t m_sig; /* Significant. */
signed int m_exp; /* Exponent. */
- bool m_negative; /* Negative sign. */
};
-extern void debug (sreal &ref);
-extern void debug (sreal *ptr);
+extern void debug (const sreal &ref);
+extern void debug (const sreal *ptr);
inline sreal &operator+= (sreal &a, const sreal &b)
{