dt/fx/scfx_rep.cc

817     double result = id;
820         result += m_sign * scfx_pow2(exp - SCFX_IEEE_DOUBLE_M_SIZE);
822     return result;
840     uint64 result = 0;
847         result += static_cast<uint64>(m_mant[idx]) << shift;
852     return m_sign > 0 ? result : -result;
1170 //  result has the same size
1175 add_mants(int size, scfx_mant &result, const scfx_mant &a, const scfx_mant &b)
1189         result[index] = y;
1196 sub_mants(int size, scfx_mant &result, const scfx_mant &a, const scfx_mant &b)
1210         result[index] = y;
1219     scfx_rep &result = *new scfx_rep;
1226         result.set_nan();
1227         return &result;
1231         result.set_inf(lhs.m_sign);
1232         return &result;
1236         result.set_inf(rhs.m_sign);
1237         return &result;
1252     // size the result mantissa
1254     result.resize_to(len_mant);
1255     result.m_wp = new_wp;
1261         add_mants(len_mant, result.m_mant, lhs_mant, rhs_mant);
1262         result.m_sign = lhs.m_sign;
1267             sub_mants(len_mant, result.m_mant, lhs_mant, rhs_mant);
1268             result.m_sign = lhs.m_sign;
1270             sub_mants(len_mant, result.m_mant, rhs_mant, lhs_mant);
1271             result.m_sign = rhs.m_sign;
1273             result.m_mant.clear();
1274             result.m_sign = 1;
1278     result.find_sw();
1279     result.round(max_wl);
1281     return &result;
1289 //  result has the same size
1328     scfx_rep &result = *new scfx_rep;
1335         result.set_nan();
1336         return &result;
1340         result.set_inf(lhs.m_sign);
1341         return &result;
1345         result.set_inf(-1 * rhs.m_sign);
1346         return &result;
1361     // size the result mantissa
1363     result.resize_to(len_mant);
1364     result.m_wp = new_wp;
1370         add_mants(len_mant, result.m_mant, lhs_mant, rhs_mant);
1371         result.m_sign = lhs.m_sign;
1376             sub_mants(len_mant, result.m_mant, lhs_mant, rhs_mant);
1377             result.m_sign = lhs.m_sign;
1379             sub_mants(len_mant, result.m_mant, rhs_mant, lhs_mant);
1380             result.m_sign = -rhs.m_sign;
1382             result.m_mant.clear();
1383             result.m_sign = 1;
1387     result.find_sw();
1388     result.round(max_wl);
1390     return &result;
1420 multiply(scfx_rep &result, const scfx_rep &lhs, const scfx_rep &rhs,
1429         result.set_nan();
1434         result.set_inf(lhs.m_sign * rhs.m_sign);
1439         result.set_zero(lhs.m_sign * rhs.m_sign);
1453     result.resize_to(new_size);
1454     result.m_mant.clear();
1455     result.m_wp = new_wp;
1456     result.m_sign = new_sign;
1457     result.m_state = scfx_rep::normal;
1462     half_word *t = result.m_mant.half_addr();
1485     result.find_sw();
1486     result.round(max_wl);
1497     scfx_rep &result = *new scfx_rep;
1504         result.set_nan();
1505         return &result;
1509         result.set_inf(lhs.m_sign * rhs.m_sign);
1510         return &result;
1514         result.set_zero(lhs.m_sign * rhs.m_sign);
1515         return &result;
1525     result.resize_to(sc_max(n_word(div_wl) + 1, min_mant));
1526     result.m_mant.clear();
1527     result.m_sign = lhs.m_sign * rhs.m_sign;
1540     result_index = (result.size() - c) * bits_in_word + msb_res % bits_in_word;
1541     result.m_wp = (result.size() - c) - msb_res / bits_in_word;
1561             result.set_bin(result_index);
1572         int index = result_index + 1 - result.m_wp * bits_in_word;
1574         scfx_index x = result.calc_indices(index);
1575         scfx_index x1 = result.calc_indices(index + 1);
1577         if (result.o_bit_at(x) && result.o_bit_at(x1))
1578             result.q_incr(x);
1580         result.m_r_flag = true;
1583     result.find_sw();
1585     return &result;
2233     scfx_mant *result = new scfx_mant(new_size);
2235     result->clear();
2240         (*result)[j + shift] = m_mant[j];
2243     return result;
2375     bool result = (m_mant[x.wi()] & (1 << x.bi())) != 0;
2379     return result;