F5/unison
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fix 128 bit to/from double conversions

Browse Source
This commit is contained in:
Mathieu Lacage
2006-12-04 12:51:57 +01:00
parent a912bad39d
commit b6acf05f80
1 changed files with 18 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

24
src/simulator/high-precision-128.cc
View File

@@ -60,12 +60,15 @@ HighPrecision::GetInteger (void) const
double
HighPrecision::GetDouble (void) const
{
cairo_int128_t value = _cairo_int128_rsa (m_value, 64);
cairo_int128_t rem = _cairo_int128_sub (m_value, value);
double frem = _cairo_int128_to_int64 (rem);
frem /= MAX_64;
double retval = _cairo_int128_to_int64 (value);
retval += frem;
bool is_negative = _cairo_int128_negative (m_value);
cairo_int128_t value = is_negative?_cairo_int128_negate (m_value):m_value;
cairo_int128_t hi = _cairo_int128_rsa (value, 64);
cairo_uint128_t lo = _cairo_int128_sub (value, _cairo_uint128_lsl (hi, 64));
double flo = _cairo_uint128_to_uint64 (lo);
flo /= MAX_64;
double retval = _cairo_uint128_to_uint64 (hi);
retval += flo;
retval *= is_negative?-1.0:1.0;
return retval;
}
bool
@@ -90,10 +93,19 @@ HighPrecision::Mul (HighPrecision const &o)
bool
HighPrecision::Div (HighPrecision const &o)
{
#if 1
cairo_int128_t div = _cairo_int128_rsa (o.m_value, 64);
cairo_quorem128_t qr;
qr = _cairo_int128_divrem (m_value, div);
m_value = qr.quo;
#else
cairo_quorem128_t qr;
qr = _cairo_int128_divrem (m_value, o.m_value);
m_value = qr.quo;
m_value = _cairo_int128_lsl (m_value, 64);
cairo_int128_t rem = _cairo_int128_rsa (qr.rem, 64);
m_value = _cairo_int128_add (m_value, rem);
#endif
return false;
}
int