mirrors/ecl
1
Fork 0
mirror of https://gitlab.com/embeddable-common-lisp/ecl.git synced 2026-01-16 22:32:05 -08:00
Code Issues Projects Releases Packages Wiki Activity

Fixed LOG of negative bignums and of complex numbers with equal positive and negative parts.

Browse source
This commit is contained in:
Juan Jose Garcia Ripoll 2009-12-08 23:26:16 +01:00
parent 19943c0a22
commit 7c050399fa
1 changed files with 9 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

16
src/c/num_sfun.d
View file

@ -235,7 +235,6 @@ ecl_log1_complex(cl_object r, cl_object i)
*/
a = ecl_times(a, a);
a = ecl_divide(ecl_log1(ecl_plus(a, a)), MAKE_FIXNUM(2));
p = MAKE_FIXNUM(0);
goto OUTPUT;
}
/* For the real part of the output we use the formula
@ -244,8 +243,8 @@ ecl_log1_complex(cl_object r, cl_object i)
a = ecl_divide(a, p);
a = ecl_plus(ecl_divide(ecl_log1p(ecl_times(a,a)), MAKE_FIXNUM(2)),
ecl_log1(p));
p = ecl_atan2(i, r);
OUTPUT:
p = ecl_atan2(i, r);
return ecl_make_complex(a, p);
}
@ -264,11 +263,14 @@ ecl_log1(cl_object x)
break;
}
case t_bignum: {
cl_fixnum l = ecl_integer_length(x) - 1;
cl_object r = ecl_make_ratio(x, ecl_ash(MAKE_FIXNUM(1), l));
float d = logf(number_to_float(r)) + l * logf(2.0);
if (d < 0) goto COMPLEX;
output = ecl_make_singlefloat(d);
if (ecl_minusp(x)) {
goto COMPLEX;
} else {
cl_fixnum l = ecl_integer_length(x) - 1;
cl_object r = ecl_make_ratio(x, ecl_ash(MAKE_FIXNUM(1), l));
float d = logf(number_to_float(r)) + l * logf(2.0);
output = ecl_make_singlefloat(d);
}
break;
}
#ifdef ECL_SHORT_FLOAT