mpdiv: negative divisor has to flip sign of quotient

This commit is contained in:
aiju 2016-08-28 12:00:25 +02:00
parent 344ff48512
commit 43bb71c8cc

View file

@ -24,9 +24,12 @@ 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(sign < 0 && quotient->top != 0)
quotient->sign ^= -2;
}
if(remainder != nil){
remainder->flags |= dividend->flags & MPtimesafe;