/*============================================================================= This file is part of ARB. ARB is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. ARB is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with ARB; if not, write to the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA =============================================================================*/ /****************************************************************************** Copyright (C) 2012 Fredrik Johansson ******************************************************************************/ #include "fmpr.h" long fmpr_exp(fmpr_t y, const fmpr_t x, slong prec, fmpr_rnd_t rnd) { if (fmpr_is_special(x)) { if (fmpr_is_zero(x)) fmpr_one(y); else if (fmpr_is_pos_inf(x)) fmpr_pos_inf(y); else if (fmpr_is_neg_inf(x)) fmpr_zero(y); else fmpr_nan(y); return FMPR_RESULT_EXACT; } else { slong r; CALL_MPFR_FUNC(r, mpfr_exp, y, x, prec, rnd); return r; } } long fmpr_expm1(fmpr_t y, const fmpr_t x, slong prec, fmpr_rnd_t rnd) { if (fmpr_is_special(x)) { if (fmpr_is_zero(x)) fmpr_zero(y); else if (fmpr_is_pos_inf(x)) fmpr_pos_inf(y); else if (fmpr_is_neg_inf(x)) fmpr_set_si(y, -1L); else fmpr_nan(y); return FMPR_RESULT_EXACT; } else { slong r; CALL_MPFR_FUNC(r, mpfr_expm1, y, x, prec, rnd); return r; } }