core: Avoid Floating point exception dividing negative numbers.
This fixes kaem --verbose --strict bin/mes-m2 -c '(display (modulo -1 -1))' * src/math.c (modulo): Use unsigned division.
This commit is contained in:
parent
7c3d55323d
commit
87e73e5b01
28
src/math.c
28
src/math.c
|
@ -188,16 +188,26 @@ modulo (struct scm *a, struct scm *b)
|
|||
{
|
||||
assert_number ("modulo", a);
|
||||
assert_number ("modulo", b);
|
||||
long x = a->value;
|
||||
long y = b->value;
|
||||
if (y == 0)
|
||||
long n = a->value;
|
||||
long v = b->value;
|
||||
if (v == 0)
|
||||
error (cstring_to_symbol ("divide-by-zero"), a);
|
||||
while (x < 0)
|
||||
x = x + y;
|
||||
if (x != 0)
|
||||
x = x % y;
|
||||
|
||||
return make_number (x);
|
||||
int sign_p = 0;
|
||||
size_t w = v;
|
||||
if (v < 0)
|
||||
{
|
||||
sign_p = 1;
|
||||
w = -v;
|
||||
}
|
||||
while (n < 0)
|
||||
n = n + w;
|
||||
size_t u = n;
|
||||
if (u != 0)
|
||||
u = u % w;
|
||||
n = u;
|
||||
if (sign_p)
|
||||
n = -n;
|
||||
return make_number (n);
|
||||
}
|
||||
|
||||
struct scm *
|
||||
|
|
Loading…
Reference in New Issue