From 1f9bdb6f5d9699493f6bce572cbbb58ef3b4ad27 Mon Sep 17 00:00:00 2001 From: cinap_lenrek Date: Sun, 28 Aug 2016 16:46:32 +0200 Subject: libmp: mpdiv(): fix divisor==quotient case (again) --- sys/src/libmp/port/mpdiv.c | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) (limited to 'sys/src') diff --git a/sys/src/libmp/port/mpdiv.c b/sys/src/libmp/port/mpdiv.c index 2c2ce73eb..ea68acb43 100644 --- a/sys/src/libmp/port/mpdiv.c +++ b/sys/src/libmp/port/mpdiv.c @@ -24,10 +24,11 @@ mpdiv(mpint *dividend, mpint *divisor, mpint *quotient, mpint *remainder) if(divisor->top == 1 && (divisor->p[0] & divisor->p[0]-1) == 0){ vlong r = (vlong)dividend->sign * (dividend->p[0] & divisor->p[0]-1); if(quotient != nil){ + sign = divisor->sign; for(s = 0; ((divisor->p[0] >> s) & 1) == 0; s++) ; mpright(dividend, s, quotient); - if(divisor->sign < 0) + if(sign < 0) quotient->sign ^= (-mpmagcmp(quotient, mpzero) >> 31) << 1; } if(remainder != nil){ -- cgit v1.2.3