hiro/arb
1
0
Fork 0
mirror of https://github.com/vale981/arb synced 2025-03-05 17:31:38 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions

add utility function

Browse source
This commit is contained in:
Fredrik Johansson 2012-04-11 10:36:09 +02:00
parent e025325b3d
commit ccc0b034ed
1 changed files with 33 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

33
ufloat.h
View file

@ -261,6 +261,39 @@ ufloat_get_fmpz(fmpz_t z, const ufloat_t u)
} }
} }
static __inline__ void
ufloat_set_mpfr(ufloat_t u, const mpfr_t x)
{
if (mpfr_zero_p(x))
{
u->man = 0;
u->exp = 0;
}
else
{
mpz_t t;
long exp, bits;
mpz_init(t);
exp = mpfr_get_z_2exp(t, x);
mpz_abs(t, t);
bits = mpz_sizeinbase(t, 2);
if (bits > UFLOAT_PREC)
{
mpz_cdiv_q_2exp(t, t, bits - UFLOAT_PREC);
exp += bits - UFLOAT_PREC;
}
u->exp = exp;
u->man = mpz_get_ui(t);
ufloat_normalise(u);
mpz_clear(t);
}
}
void ufloat_log(ufloat_t z, const ufloat_t x); void ufloat_log(ufloat_t z, const ufloat_t x);
void ufloat_log1p(ufloat_t z, const ufloat_t x); void ufloat_log1p(ufloat_t z, const ufloat_t x);