mirror/guile
mirror/guile
1
Fork 0
mirror of https://git.savannah.gnu.org/git/guile.git synced 2025-05-20 11:40:18 +02:00
Code Activity

Expose frexp from integers lib

Browse source 
* libguile/integers.h:
* libguile/integers.c (scm_integer_frexp_z): Expose internally.
(scm_integer_to_double_z):
(scm_integer_inexact_sqrt_z): Adapt uses.
* libguile/numbers.c (log_of_exact_integer): Use exposed
scm_integer_frexp_z.
This commit is contained in:
Andy Wingo 2022-01-07 10:11:45 +01:00
parent 124d889227
commit f167627bb0
3 changed files with 7 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

8
libguile/integers.c
View file

@ -2696,8 +2696,8 @@ mpz_get_d_2exp (long *exp, mpz_srcptr z)
}
#endif
static double
bignum_frexp (struct scm_bignum *x, long *exp)
double
scm_integer_frexp_z (struct scm_bignum *x, long *exp)
{
mpz_t zx;
alias_bignum_to_mpz (x, zx);
@ -2730,7 +2730,7 @@ double
scm_integer_to_double_z (struct scm_bignum *x)
{
long exponent;
double significand = bignum_frexp (x, &exponent);
double significand = scm_integer_frexp_z (x, &exponent);
return ldexp (significand, exponent);
}
@ -3134,7 +3134,7 @@ scm_integer_inexact_sqrt_z (struct scm_bignum *k)
mpz_init (zs);
long expon;
double signif = bignum_frexp (k, &expon);
double signif = scm_integer_frexp_z (k, &expon);
int negative = signif < 0;
if (negative)
signif = -signif;

1
libguile/integers.h
View file

@ -169,6 +169,7 @@ SCM_INTERNAL int scm_is_integer_less_than_rz (double y, struct scm_bignum *x);
SCM_INTERNAL int scm_is_integer_positive_z (struct scm_bignum *x);
SCM_INTERNAL int scm_is_integer_negative_z (struct scm_bignum *x);
SCM_INTERNAL double scm_integer_frexp_z (struct scm_bignum *x, long *exp);
SCM_INTERNAL double scm_integer_to_double_z (struct scm_bignum *x);
SCM_INTERNAL SCM scm_integer_from_double (double val);

4
libguile/numbers.c
View file

@ -7209,11 +7209,11 @@ log_of_exact_integer (SCM n)
else if (SCM_BIGP (n))
{
long expon;
double signif = scm_i_big2dbl_2exp (n, &expon);
double signif = scm_integer_frexp_z (scm_bignum (n), &expon);
return log_of_shifted_double (signif, expon);
}
else
scm_wrong_type_arg ("log_of_exact_integer", SCM_ARG1, n);
abort ();
}
/* Returns log(n/d), for exact non-zero integers n and d */