Home
last modified time | relevance | path

Searched refs:rnd_mode (Results 1 – 25 of 118) sorted by relevance

12345

/dragonfly/contrib/mpfr/src/
HDatan2.c27 pi_div_2ui (mpfr_ptr dest, int i, int neg, mpfr_rnd_t rnd_mode) in pi_div_2ui() argument
35 inexact = - mpfr_const_pi (dest, MPFR_INVERT_RND (rnd_mode)); in pi_div_2ui()
40 inexact = mpfr_const_pi (dest, rnd_mode); in pi_div_2ui()
42 mpfr_div_2ui (dest, dest, i, rnd_mode); /* exact */ in pi_div_2ui()
44 return mpfr_check_range (dest, inexact, rnd_mode); in pi_div_2ui()
48 mpfr_atan2 (mpfr_ptr dest, mpfr_srcptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode) in mpfr_atan2() argument
60 mpfr_get_prec (x), mpfr_log_prec, x, rnd_mode), in mpfr_atan2()
94 inexact = mpfr_const_pi (dest, MPFR_INVERT_RND (rnd_mode)); in mpfr_atan2()
99 return mpfr_const_pi (dest, rnd_mode); in mpfr_atan2()
111 return pi_div_2ui (dest, 1, MPFR_IS_NEG (y), rnd_mode); in mpfr_atan2()
[all …]
HDminmax.c33 mpfr_min (mpfr_ptr z, mpfr_srcptr x, mpfr_srcptr y, mpfr_rnd_t rnd_mode) in mpfr_min() argument
43 return mpfr_set(z, y, rnd_mode); in mpfr_min()
45 return mpfr_set(z, x, rnd_mode); in mpfr_min()
49 return mpfr_set(z, x, rnd_mode); in mpfr_min()
51 return mpfr_set(z, y, rnd_mode); in mpfr_min()
55 return mpfr_set(z, x, rnd_mode); in mpfr_min()
57 return mpfr_set(z, y, rnd_mode); in mpfr_min()
67 mpfr_max (mpfr_ptr z, mpfr_srcptr x, mpfr_srcptr y, mpfr_rnd_t rnd_mode) in mpfr_max() argument
77 return mpfr_set(z, y, rnd_mode); in mpfr_max()
79 return mpfr_set(z, x, rnd_mode); in mpfr_max()
[all …]
HDrint.c28 mpfr_rint (mpfr_ptr r, mpfr_srcptr u, mpfr_rnd_t rnd_mode) in mpfr_rint() argument
60 rnd_mode == MPFR_RNDD ? sign < 0 : in mpfr_rint()
61 rnd_mode == MPFR_RNDU ? sign > 0 : in mpfr_rint()
62 rnd_mode == MPFR_RNDZ ? 0 : in mpfr_rint()
63 rnd_mode == MPFR_RNDA ? 1 : in mpfr_rint()
77 (exp == 0 && (rnd_mode == MPFR_RNDNA || in mpfr_rint()
154 if (rnd_mode == MPFR_RNDN && in mpfr_rint()
239 else if (rnd_mode == MPFR_RNDN && in mpfr_rint()
289 return mpfr_overflow(r, rnd_mode, MPFR_SIGN(r)) >= 0 ? in mpfr_rint()
337 mpfr_rint_round (mpfr_ptr r, mpfr_srcptr u, mpfr_rnd_t rnd_mode) in mpfr_rint_round() argument
[all …]
HDurandom.c42 mpfr_urandom (mpfr_ptr rop, gmp_randstate_t rstate, mpfr_rnd_t rnd_mode) in mpfr_urandom() argument
61 if (rnd_mode == MPFR_RNDU || rnd_mode == MPFR_RNDA in mpfr_urandom()
62 || (emin == 1 && rnd_mode == MPFR_RNDN in mpfr_urandom()
65 mpfr_set_ui_2exp (rop, 1, emin - 1, rnd_mode); in mpfr_urandom()
101 if (rnd_mode == MPFR_RNDU || rnd_mode == MPFR_RNDA in mpfr_urandom()
102 || (rnd_mode == MPFR_RNDN && cnt == exp - emin - 1 in mpfr_urandom()
105 mpfr_set_ui_2exp (rop, 1, emin - 1, rnd_mode); in mpfr_urandom()
131 if (rnd_mode == MPFR_RNDU || rnd_mode == MPFR_RNDA in mpfr_urandom()
132 || (rnd_mode == MPFR_RNDN && random_rounding_bit (rstate))) in mpfr_urandom()
136 mpfr_set_ui_2exp (rop, 1, emin - 1, rnd_mode); in mpfr_urandom()
[all …]
HDget_flt.c35 mpfr_get_flt (mpfr_srcptr src, mpfr_rnd_t rnd_mode) in mpfr_get_flt() argument
44 return (float) mpfr_get_d (src, rnd_mode); in mpfr_get_flt()
49 if (MPFR_UNLIKELY(rnd_mode == MPFR_RNDA)) in mpfr_get_flt()
50 rnd_mode = negative ? MPFR_RNDD : MPFR_RNDU; in mpfr_get_flt()
61 (rnd_mode == MPFR_RNDD || in mpfr_get_flt()
62 (rnd_mode == MPFR_RNDN && mpfr_cmp_si_2exp (src, -1, -150) < 0) in mpfr_get_flt()
64 (rnd_mode == MPFR_RNDU || in mpfr_get_flt()
65 (rnd_mode == MPFR_RNDN && mpfr_cmp_si_2exp (src, 1, -150) > 0) in mpfr_get_flt()
75 (rnd_mode == MPFR_RNDZ || rnd_mode == MPFR_RNDU ? in mpfr_get_flt()
77 (rnd_mode == MPFR_RNDZ || rnd_mode == MPFR_RNDD ? in mpfr_get_flt()
[all …]
HDget_d.c36 mpfr_get_d (mpfr_srcptr src, mpfr_rnd_t rnd_mode) in mpfr_get_d() argument
59 if (MPFR_UNLIKELY(rnd_mode == MPFR_RNDA)) in mpfr_get_d()
60 rnd_mode = negative ? MPFR_RNDD : MPFR_RNDU; in mpfr_get_d()
70 (rnd_mode == MPFR_RNDD || in mpfr_get_d()
71 (rnd_mode == MPFR_RNDN && mpfr_cmp_si_2exp(src, -1, -1075) < 0) in mpfr_get_d()
73 (rnd_mode == MPFR_RNDU || in mpfr_get_d()
74 (rnd_mode == MPFR_RNDN && mpfr_cmp_si_2exp(src, 1, -1075) > 0) in mpfr_get_d()
84 (rnd_mode == MPFR_RNDZ || rnd_mode == MPFR_RNDU ? in mpfr_get_d()
86 (rnd_mode == MPFR_RNDZ || rnd_mode == MPFR_RNDD ? in mpfr_get_d()
106 nbits, rnd_mode); in mpfr_get_d()
[all …]
HDsub.c26 mpfr_sub (mpfr_ptr a, mpfr_srcptr b, mpfr_srcptr c, mpfr_rnd_t rnd_mode) in mpfr_sub() argument
31 mpfr_get_prec (c), mpfr_log_prec, c, rnd_mode), in mpfr_sub()
65 int sign = rnd_mode != MPFR_RNDD in mpfr_sub()
73 return mpfr_neg (a, c, rnd_mode); in mpfr_sub()
78 return mpfr_set (a, b, rnd_mode); in mpfr_sub()
89 return mpfr_sub1sp (a, b, c, rnd_mode); in mpfr_sub()
91 return mpfr_sub1 (a, b, c, rnd_mode); in mpfr_sub()
98 rnd_mode = MPFR_INVERT_RND (rnd_mode); in mpfr_sub()
101 inexact = mpfr_add1sp (a, c, b, rnd_mode); in mpfr_sub()
103 inexact = mpfr_add1 (a, c, b, rnd_mode); in mpfr_sub()
[all …]
HDget_d64.c289 mpfr_get_decimal64 (mpfr_srcptr src, mpfr_rnd_t rnd_mode) in mpfr_get_decimal64() argument
312 if (MPFR_UNLIKELY(rnd_mode == MPFR_RNDA)) in mpfr_get_decimal64()
313 rnd_mode = negative ? MPFR_RNDD : MPFR_RNDU; in mpfr_get_decimal64()
319 if (rnd_mode == MPFR_RNDZ || rnd_mode == MPFR_RNDN in mpfr_get_decimal64()
320 || (rnd_mode == MPFR_RNDD && negative == 0) in mpfr_get_decimal64()
321 || (rnd_mode == MPFR_RNDU && negative != 0)) in mpfr_get_decimal64()
329 if (rnd_mode == MPFR_RNDZ in mpfr_get_decimal64()
330 || (rnd_mode == MPFR_RNDU && negative != 0) in mpfr_get_decimal64()
331 || (rnd_mode == MPFR_RNDD && negative == 0)) in mpfr_get_decimal64()
341 mpfr_get_str (s, &e, 10, 16, src, rnd_mode); in mpfr_get_decimal64()
[all …]
HDpow.c34 mpfr_rnd_t rnd_mode, int *inexact) in mpfr_pow_is_exact() argument
100 *inexact = mpfr_pow_z (z, tmp, c, rnd_mode); in mpfr_pow_is_exact()
161 mpfr_rnd_t rnd_mode, int y_is_integer, mpfr_save_expo_t *expo) in mpfr_pow_general() argument
178 mpfr_get_prec (y), mpfr_log_prec, y, rnd_mode), in mpfr_pow_general()
191 rnd_mode = MPFR_INVERT_RND (rnd_mode); in mpfr_pow_general()
251 inexact = mpfr_underflow (z, rnd_mode == MPFR_RNDN ? MPFR_RNDZ in mpfr_pow_general()
252 : rnd_mode, MPFR_SIGN_POS); in mpfr_pow_general()
272 inexact = mpfr_overflow (z, rnd_mode, MPFR_SIGN_POS); in mpfr_pow_general()
296 if (MPFR_LIKELY (MPFR_CAN_ROUND (t, Nt - err, Nz, rnd_mode))) in mpfr_pow_general()
298 inexact = mpfr_set (z, t, rnd_mode); in mpfr_pow_general()
[all …]
HDhypot.c30 mpfr_hypot (mpfr_ptr z, mpfr_srcptr x, mpfr_srcptr y, mpfr_rnd_t rnd_mode) in mpfr_hypot() argument
47 mpfr_get_prec (y), mpfr_log_prec, y, rnd_mode), in mpfr_hypot()
67 return mpfr_abs (z, y, rnd_mode); in mpfr_hypot()
69 return mpfr_abs (z, x, rnd_mode); in mpfr_hypot()
87 threshold = (MAX (N, Nz) + (rnd_mode == MPFR_RNDN ? 1 : 0)) << 1; in mpfr_hypot()
88 if (rnd_mode == MPFR_RNDA) in mpfr_hypot()
89 rnd_mode = MPFR_RNDU; /* since the result is positive, RNDA = RNDU */ in mpfr_hypot()
96 if (MPFR_UNLIKELY (rnd_mode == MPFR_RNDU)) in mpfr_hypot()
100 if (mpfr_abs (z, x, rnd_mode) == 0) in mpfr_hypot()
108 mpfr_abs (z, x, rnd_mode); /* exact */ in mpfr_hypot()
[all …]
HDsi_op.c28 mpfr_add_si (mpfr_ptr y, mpfr_srcptr x, long int u, mpfr_rnd_t rnd_mode) in mpfr_add_si() argument
31 return mpfr_add_ui (y, x, u, rnd_mode); in mpfr_add_si()
33 return mpfr_sub_ui (y, x, -u, rnd_mode); in mpfr_add_si()
37 mpfr_sub_si (mpfr_ptr y, mpfr_srcptr x, long int u, mpfr_rnd_t rnd_mode) in mpfr_sub_si() argument
40 return mpfr_sub_ui (y, x, u, rnd_mode); in mpfr_sub_si()
42 return mpfr_add_ui (y, x, -u, rnd_mode); in mpfr_sub_si()
46 mpfr_si_sub (mpfr_ptr y, long int u, mpfr_srcptr x, mpfr_rnd_t rnd_mode) in mpfr_si_sub() argument
49 return mpfr_ui_sub (y, u, x, rnd_mode); in mpfr_si_sub()
52 int res = -mpfr_add_ui (y, x, -u, MPFR_INVERT_RND (rnd_mode)); in mpfr_si_sub()
HDexp.c28 mpfr_exp (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode) in mpfr_exp() argument
36 (("x[%Pu]=%.*Rg rnd=%d", mpfr_get_prec (x), mpfr_log_prec, x, rnd_mode), in mpfr_exp()
59 return mpfr_set_ui (y, 1, rnd_mode); in mpfr_exp()
82 return mpfr_overflow (y, rnd_mode, 1); in mpfr_exp()
96 return mpfr_underflow (y, rnd_mode == MPFR_RNDN ? MPFR_RNDZ : rnd_mode, in mpfr_exp()
117 if (MPFR_IS_NEG_SIGN (signx) && (rnd_mode == MPFR_RNDD || in mpfr_exp()
118 rnd_mode == MPFR_RNDZ)) in mpfr_exp()
130 if (MPFR_IS_POS_SIGN (signx) && (rnd_mode == MPFR_RNDU || in mpfr_exp()
131 rnd_mode == MPFR_RNDA)) in mpfr_exp()
153 inexact = mpfr_exp_3 (y, x, rnd_mode); /* O(M(n) log(n)^2) */ in mpfr_exp()
[all …]
HDadd.c26 mpfr_add (mpfr_ptr a, mpfr_srcptr b, mpfr_srcptr c, mpfr_rnd_t rnd_mode) in mpfr_add() argument
31 mpfr_get_prec (c), mpfr_log_prec, c, rnd_mode), in mpfr_add()
71 (rnd_mode != MPFR_RNDD ? in mpfr_add()
77 return mpfr_set (a, c, rnd_mode); in mpfr_add()
82 return mpfr_set (a, b, rnd_mode); in mpfr_add()
93 return mpfr_sub1sp(a, b, c, rnd_mode); in mpfr_add()
95 return mpfr_sub1(a, b, c, rnd_mode); in mpfr_add()
102 return mpfr_add1sp(a, c, b, rnd_mode); in mpfr_add()
104 return mpfr_add1sp(a, b, c, rnd_mode); in mpfr_add()
107 return mpfr_add1(a, c, b, rnd_mode); in mpfr_add()
[all …]
HDdiv_ui.c28 mpfr_div_ui (mpfr_ptr y, mpfr_srcptr x, unsigned long int u, mpfr_rnd_t rnd_mode) in mpfr_div_ui() argument
40 mpfr_get_prec(x), mpfr_log_prec, x, u, rnd_mode), in mpfr_div_ui()
84 return mpfr_set (y, x, rnd_mode); in mpfr_div_ui()
87 return mpfr_div_2si (y, x, MPFR_INT_CEIL_LOG2 (u), rnd_mode); in mpfr_div_ui()
120 if (MPFR_LIKELY (rnd_mode == MPFR_RNDN)) in mpfr_div_ui()
193 return mpfr_underflow (y, rnd_mode == MPFR_RNDN ? MPFR_RNDZ : rnd_mode, in mpfr_div_ui()
200 MPFR_UPDATE2_RND_MODE(rnd_mode, MPFR_SIGN (y)); in mpfr_div_ui()
201 switch (rnd_mode) in mpfr_div_ui()
214 MPFR_ASSERTD (rnd_mode == MPFR_RNDN); in mpfr_div_ui()
259 return mpfr_check_range (y, inexact, rnd_mode); in mpfr_div_ui()
[all …]
HDexpm1.c31 mpfr_expm1 (mpfr_ptr y, mpfr_srcptr x , mpfr_rnd_t rnd_mode) in mpfr_expm1() argument
38 (("x[%Pu]=%.*Rg rnd=%d", mpfr_get_prec (x), mpfr_log_prec, x, rnd_mode), in mpfr_expm1()
59 return mpfr_set_si (y, -1, rnd_mode); in mpfr_expm1()
76 MPFR_FAST_COMPUTE_IF_SMALL_INPUT (y, x, - ex, 0, 1, rnd_mode, {}); in mpfr_expm1()
78 MPFR_FAST_COMPUTE_IF_SMALL_INPUT (y, x, - ex, 1, 0, rnd_mode, {}); in mpfr_expm1()
101 MPFR_SMALL_INPUT_AFTER_SAVE_EXPO (y, minus_one, err, 0, 0, rnd_mode, in mpfr_expm1()
138 inexact = mpfr_overflow (y, rnd_mode, MPFR_SIGN_POS); in mpfr_expm1()
144 inexact = mpfr_set_si (y, -1, rnd_mode); in mpfr_expm1()
147 if (MPFR_IS_LIKE_RNDZ (rnd_mode, 1)) in mpfr_expm1()
162 if (MPFR_LIKELY (MPFR_CAN_ROUND (t, err, Ny, rnd_mode))) in mpfr_expm1()
[all …]
HDsub1sp.c31 int mpfr_sub1sp2 (mpfr_ptr a, mpfr_srcptr b, mpfr_srcptr c, mpfr_rnd_t rnd_mode);
32 int mpfr_sub1sp (mpfr_ptr a, mpfr_srcptr b, mpfr_srcptr c, mpfr_rnd_t rnd_mode) in mpfr_sub1sp() argument
47 inexact2 = mpfr_sub1 (tmpa, tmpb, tmpc, rnd_mode); in mpfr_sub1sp()
48 inexact = mpfr_sub1sp2(a, b, c, rnd_mode); in mpfr_sub1sp()
54 mpfr_print_rnd_mode (rnd_mode), (unsigned long) MPFR_PREC (a), in mpfr_sub1sp()
135 mpfr_sub1sp (mpfr_ptr a, mpfr_srcptr b, mpfr_srcptr c, mpfr_rnd_t rnd_mode) in mpfr_sub1sp() argument
176 if (rnd_mode == MPFR_RNDD) in mpfr_sub1sp()
274 if (rnd_mode == MPFR_RNDN && in mpfr_sub1sp()
277 rnd_mode = MPFR_RNDZ; in mpfr_sub1sp()
278 return mpfr_underflow (a, rnd_mode, MPFR_SIGN(a)); in mpfr_sub1sp()
[all …]
HDroot.c42 mpfr_root (mpfr_ptr y, mpfr_srcptr x, unsigned long k, mpfr_rnd_t rnd_mode) in mpfr_root() argument
52 mpfr_get_prec (x), mpfr_log_prec, x, k, rnd_mode), in mpfr_root()
73 return mpfr_set_ui (y, 1, rnd_mode); /* 1 may be Out of Range */ in mpfr_root()
75 return mpfr_set_ui (y, 0, rnd_mode); /* 0+ */ in mpfr_root()
88 return mpfr_set (y, x, rnd_mode); in mpfr_root()
143 n = MPFR_PREC (y) + (rnd_mode == MPFR_RNDN); in mpfr_root()
182 rnd_mode = MPFR_INVERT_RND (rnd_mode); in mpfr_root()
183 if (rnd_mode == MPFR_RNDU || rnd_mode == MPFR_RNDA in mpfr_root()
184 || (rnd_mode == MPFR_RNDN && mpz_tstbit (m, 0))) in mpfr_root()
204 return mpfr_check_range (y, inexact, rnd_mode); in mpfr_root()
HDli2.c34 li2_series (mpfr_t sum, mpfr_srcptr z, mpfr_rnd_t rnd_mode) in li2_series() argument
98 if (MPFR_CAN_ROUND (s, (mpfr_exp_t) p - err, sump, rnd_mode)) in li2_series()
108 mpfr_set (sum, s, rnd_mode); in li2_series()
135 mpfr_li2_asympt_pos (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode) in mpfr_li2_asympt_pos() argument
163 MPFR_CAN_ROUND (g, w - 4, MPFR_PREC (y), rnd_mode)) in mpfr_li2_asympt_pos()
164 inex = mpfr_set (y, g, rnd_mode); in mpfr_li2_asympt_pos()
181 mpfr_li2_asympt_neg (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode) in mpfr_li2_asympt_neg() argument
206 MPFR_CAN_ROUND (g, w - 4, MPFR_PREC (y), rnd_mode)) in mpfr_li2_asympt_neg()
207 inex = mpfr_neg (y, g, rnd_mode); in mpfr_li2_asympt_neg()
218 mpfr_li2 (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode) in mpfr_li2() argument
[all …]
HDdigamma.c121 mpfr_digamma_reflection (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode) in mpfr_digamma_reflection() argument
148 inex = mpfr_digamma (y, u, rnd_mode); in mpfr_digamma_reflection()
182 if (MPFR_CAN_ROUND (v, p - e1, MPFR_PREC(y), rnd_mode)) in mpfr_digamma_reflection()
190 inex = mpfr_set (y, v, rnd_mode); in mpfr_digamma_reflection()
202 mpfr_digamma_positive (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode) in mpfr_digamma_positive() argument
264 if (MPFR_CAN_ROUND (t, p - errt, MPFR_PREC(y), rnd_mode)) in mpfr_digamma_positive()
271 inex = mpfr_set (y, t, rnd_mode); in mpfr_digamma_positive()
279 mpfr_digamma (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode) in mpfr_digamma() argument
285 (("x[%Pu]=%.*Rg rnd=%d", mpfr_get_prec(x), mpfr_log_prec, x, rnd_mode), in mpfr_digamma()
345 inex = mpfr_si_div (y, -1, x, rnd_mode); in mpfr_digamma()
[all …]
HDsub1.c32 mpfr_sub1 (mpfr_ptr a, mpfr_srcptr b, mpfr_srcptr c, mpfr_rnd_t rnd_mode) in mpfr_sub1() argument
53 if (rnd_mode == MPFR_RNDD) in mpfr_sub1()
99 rnd_mode, MPFR_SIGN (a), in mpfr_sub1()
101 inexact = mpfr_overflow (a, rnd_mode, MPFR_SIGN (a))); in mpfr_sub1()
109 if (MPFR_IS_LIKE_RNDZ (rnd_mode, MPFR_IS_NEG (a))) in mpfr_sub1()
202 mpfr_print_rnd_mode (rnd_mode), shift_b, shift_c, in mpfr_sub1()
329 if (MPFR_LIKELY(rnd_mode == MPFR_RNDN)) in mpfr_sub1()
349 if (MPFR_IS_RNDUTEST_OR_RNDDNOTTEST(rnd_mode, MPFR_IS_NEG(a))) in mpfr_sub1()
350 rnd_mode = MPFR_RNDZ; in mpfr_sub1()
354 if (rnd_mode == MPFR_RNDZ) in mpfr_sub1()
[all …]
HDdiv_2si.c26 mpfr_div_2si (mpfr_ptr y, mpfr_srcptr x, long int n, mpfr_rnd_t rnd_mode) in mpfr_div_2si() argument
32 mpfr_get_prec(x), mpfr_log_prec, x, n, rnd_mode), in mpfr_div_2si()
37 return mpfr_set (y, x, rnd_mode); in mpfr_div_2si()
41 MPFR_SETRAW (inexact, y, x, exp, rnd_mode); in mpfr_div_2si()
45 if (rnd_mode == MPFR_RNDN && in mpfr_div_2si()
49 rnd_mode = MPFR_RNDZ; in mpfr_div_2si()
50 return mpfr_underflow (y, rnd_mode, MPFR_SIGN(y)); in mpfr_div_2si()
54 return mpfr_overflow (y, rnd_mode, MPFR_SIGN(y)); in mpfr_div_2si()
HDmul_2si.c26 mpfr_mul_2si (mpfr_ptr y, mpfr_srcptr x, long int n, mpfr_rnd_t rnd_mode) in mpfr_mul_2si() argument
32 mpfr_get_prec (x), mpfr_log_prec, x, n, rnd_mode), in mpfr_mul_2si()
37 return mpfr_set (y, x, rnd_mode); in mpfr_mul_2si()
41 MPFR_SETRAW (inexact, y, x, exp, rnd_mode); in mpfr_mul_2si()
44 return mpfr_overflow (y, rnd_mode, MPFR_SIGN(y)); in mpfr_mul_2si()
48 if (rnd_mode == MPFR_RNDN && in mpfr_mul_2si()
52 rnd_mode = MPFR_RNDZ; in mpfr_mul_2si()
53 return mpfr_underflow (y, rnd_mode, MPFR_SIGN(y)); in mpfr_mul_2si()
HDsinh.c30 mpfr_sinh (mpfr_ptr y, mpfr_srcptr xt, mpfr_rnd_t rnd_mode) in mpfr_sinh() argument
36 (("x[%Pu]=%.*Rg rnd=%d", mpfr_get_prec (xt), mpfr_log_prec, xt, rnd_mode), in mpfr_sinh()
64 rnd_mode, {}); in mpfr_sinh()
110 inexact = mpfr_overflow (y, rnd_mode, MPFR_SIGN (xt)); in mpfr_sinh()
123 inexact = mpfr_overflow (y, rnd_mode, MPFR_SIGN (xt)); in mpfr_sinh()
132 inexact = mpfr_overflow (y, rnd_mode, MPFR_SIGN (xt)); in mpfr_sinh()
140 rnd_mode))) in mpfr_sinh()
142 inexact = mpfr_set4 (y, t, rnd_mode, MPFR_SIGN (xt)); in mpfr_sinh()
165 rnd_mode))) in mpfr_sinh()
167 inexact = mpfr_set4 (y, t, rnd_mode, MPFR_SIGN (xt)); in mpfr_sinh()
[all …]
HDsinh_cosh.c35 mpfr_sinh_cosh (mpfr_ptr sh, mpfr_ptr ch, mpfr_srcptr xt, mpfr_rnd_t rnd_mode) in mpfr_sinh_cosh() argument
44 mpfr_get_prec (xt), mpfr_log_prec, xt, rnd_mode), in mpfr_sinh_cosh()
71 inexact_ch = mpfr_set_ui (ch, 1, rnd_mode); /* cosh(0) = 1 */ in mpfr_sinh_cosh()
113 inexact_ch = mpfr_overflow (ch, rnd_mode, MPFR_SIGN_POS); in mpfr_sinh_cosh()
115 inexact_sh = mpfr_sinh (sh, xt, rnd_mode); in mpfr_sinh_cosh()
137 rnd_mode) && \ in mpfr_sinh_cosh()
139 rnd_mode))) in mpfr_sinh_cosh()
141 inexact_sh = mpfr_set4 (sh, s, rnd_mode, MPFR_SIGN (xt)); in mpfr_sinh_cosh()
142 inexact_ch = mpfr_set (ch, c, rnd_mode); in mpfr_sinh_cosh()
157 inexact_sh = mpfr_check_range (sh, inexact_sh, rnd_mode); in mpfr_sinh_cosh()
[all …]
HDexp2.c30 mpfr_exp2 (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode) in mpfr_exp2() argument
38 (("x[%Pu]=%.*Rg rnd=%d", mpfr_get_prec(x), mpfr_log_prec, x, rnd_mode), in mpfr_exp2()
61 return mpfr_set_ui (y, 1, rnd_mode); in mpfr_exp2()
70 mpfr_rnd_t rnd2 = rnd_mode; in mpfr_exp2()
72 if (rnd_mode == MPFR_RNDN && in mpfr_exp2()
80 return mpfr_overflow (y, rnd_mode, 1); in mpfr_exp2()
90 MPFR_SIGN(x) > 0, rnd_mode, expo, {}); in mpfr_exp2()
129 if (MPFR_LIKELY (MPFR_CAN_ROUND (t, err, Ny, rnd_mode))) in mpfr_exp2()
138 inexact = mpfr_set (y, t, rnd_mode); in mpfr_exp2()
150 return mpfr_check_range (y, inexact, rnd_mode); in mpfr_exp2()

12345