arb/arf/set_mpfr.c
2016-04-26 17:20:05 +02:00

35 lines
932 B
C

/*
Copyright (C) 2014 Fredrik Johansson
This file is part of Arb.
Arb is free software: you can redistribute it and/or modify it under
the terms of the GNU Lesser General Public License (LGPL) as published
by the Free Software Foundation; either version 2.1 of the License, or
(at your option) any later version. See <http://www.gnu.org/licenses/>.
*/
#include "arf.h"
void
arf_set_mpfr(arf_t x, const mpfr_t y)
{
if (!mpfr_regular_p(y))
{
if (mpfr_zero_p(y))
arf_zero(x);
else if (mpfr_inf_p(y) && mpfr_sgn(y) > 0)
arf_pos_inf(x);
else if (mpfr_inf_p(y) && mpfr_sgn(y) < 0)
arf_neg_inf(x);
else
arf_nan(x);
}
else
{
mp_size_t n = (y->_mpfr_prec + FLINT_BITS - 1) / FLINT_BITS;
arf_set_mpn(x, y->_mpfr_d, n, y->_mpfr_sign < 0);
fmpz_set_si(ARF_EXPREF(x), y->_mpfr_exp);
}
}