cheng/mpir
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Small improvement to sb_div_q.

Browse Source
This commit is contained in:
William Hart 2014-02-25 15:58:30 +00:00
parent bcf545c219
commit c7cdb70393
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
mpn/generic/sb_div_q.c
View File

@ -31,7 +31,7 @@ along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */
#include "gmp-impl.h"
#include "longlong.h"
#define SB_DIVAPPR_Q_SMALL_THRESHOLD 30
#define SB_DIVAPPR_Q_SMALL_THRESHOLD 11
void __div_helper(mp_ptr qp, mp_ptr np, mp_srcptr dp, mp_size_t qn)
{
@ -121,12 +121,15 @@ mpn_sb_div_q (mp_ptr qp,
np--;
/* rare case where truncation ruins normalisation */
if (cy > dp[qn] || (cy == dp[qn] && mpn_cmp(np - qn + 1, dp, qn) >= 0))
if (cy >= dp[qn])
{
if (cy > dp[qn] || mpn_cmp(np - qn + 1, dp, qn) >= 0)
{
__div_helper(qp, np - qn, dp, qn);
flag = 0;
break;
}
}
mpir_divapprox32_preinv2(q, cy, np[0], d1inv);