1/* mpfr_stack -- initialize a floating-point number with given allocation area 2 3Copyright 2005, 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2013 Free Software Foundation, Inc. 4Contributed by the AriC and Caramel projects, INRIA. 5 6This file is part of the GNU MPFR Library. 7 8The GNU MPFR Library is free software; you can redistribute it and/or modify 9it under the terms of the GNU Lesser General Public License as published by 10the Free Software Foundation; either version 3 of the License, or (at your 11option) any later version. 12 13The GNU MPFR Library is distributed in the hope that it will be useful, but 14WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 15or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public 16License for more details. 17 18You should have received a copy of the GNU Lesser General Public License 19along with the GNU MPFR Library; see the file COPYING.LESSER. If not, see 20http://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., 2151 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA. */ 22 23#include "mpfr-impl.h" 24 25#undef mpfr_custom_get_size 26size_t 27mpfr_custom_get_size (mpfr_prec_t prec) 28{ 29 return MPFR_PREC2LIMBS (prec) * BYTES_PER_MP_LIMB; 30} 31 32#undef mpfr_custom_init 33void 34mpfr_custom_init (void *mantissa, mpfr_prec_t prec) 35{ 36 return ; 37} 38 39#undef mpfr_custom_get_significand 40void * 41mpfr_custom_get_significand (mpfr_srcptr x) 42{ 43 return (void*) MPFR_MANT (x); 44} 45 46#undef mpfr_custom_get_exp 47mpfr_exp_t 48mpfr_custom_get_exp (mpfr_srcptr x) 49{ 50 return MPFR_EXP (x); 51} 52 53#undef mpfr_custom_move 54void 55mpfr_custom_move (mpfr_ptr x, void *new_position) 56{ 57 MPFR_MANT (x) = (mp_limb_t *) new_position; 58} 59 60#undef mpfr_custom_init_set 61void 62mpfr_custom_init_set (mpfr_ptr x, int kind, mpfr_exp_t exp, 63 mpfr_prec_t prec, void *mantissa) 64{ 65 mpfr_kind_t t; 66 int s; 67 mpfr_exp_t e; 68 69 if (kind >= 0) 70 { 71 t = (mpfr_kind_t) kind; 72 s = MPFR_SIGN_POS; 73 } 74 else 75 { 76 t = (mpfr_kind_t) -kind; 77 s = MPFR_SIGN_NEG; 78 } 79 MPFR_ASSERTD (t <= MPFR_REGULAR_KIND); 80 e = MPFR_LIKELY (t == MPFR_REGULAR_KIND) ? exp : 81 MPFR_UNLIKELY (t == MPFR_NAN_KIND) ? MPFR_EXP_NAN : 82 MPFR_UNLIKELY (t == MPFR_INF_KIND) ? MPFR_EXP_INF : MPFR_EXP_ZERO; 83 84 MPFR_PREC (x) = prec; 85 MPFR_SET_SIGN (x, s); 86 MPFR_EXP (x) = e; 87 MPFR_MANT (x) = (mp_limb_t*) mantissa; 88 return; 89} 90 91#undef mpfr_custom_get_kind 92int 93mpfr_custom_get_kind (mpfr_srcptr x) 94{ 95 if (MPFR_LIKELY (!MPFR_IS_SINGULAR (x))) 96 return (int) MPFR_REGULAR_KIND * MPFR_INT_SIGN (x); 97 if (MPFR_IS_INF (x)) 98 return (int) MPFR_INF_KIND * MPFR_INT_SIGN (x); 99 if (MPFR_IS_NAN (x)) 100 return (int) MPFR_NAN_KIND; 101 MPFR_ASSERTD (MPFR_IS_ZERO (x)); 102 return (int) MPFR_ZERO_KIND * MPFR_INT_SIGN (x); 103} 104 105