EliasFleckenstein03
/
plan9front
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

libmp: mpdiv(): fix divisor==quotient case (again)

front
cinap_lenrek 2016-08-28 16:46:32 +02:00
parent bdc87e17bd
commit 1f9bdb6f5d
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
sys/src/libmp/port/mpdiv.c
View File

@ -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){