summaryrefslogtreecommitdiff
path: root/Build/source/libs/mpfr/mpfr-src/src/add_ui.c
diff options
context:
space:
mode:
Diffstat (limited to 'Build/source/libs/mpfr/mpfr-src/src/add_ui.c')
-rw-r--r--Build/source/libs/mpfr/mpfr-src/src/add_ui.c70
1 files changed, 46 insertions, 24 deletions
diff --git a/Build/source/libs/mpfr/mpfr-src/src/add_ui.c b/Build/source/libs/mpfr/mpfr-src/src/add_ui.c
index ef3b5ec923f..e6ad5db7831 100644
--- a/Build/source/libs/mpfr/mpfr-src/src/add_ui.c
+++ b/Build/source/libs/mpfr/mpfr-src/src/add_ui.c
@@ -23,7 +23,7 @@ http://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc.,
#define MPFR_NEED_LONGLONG_H
#include "mpfr-impl.h"
-int
+MPFR_HOT_FUNCTION_ATTR int
mpfr_add_ui (mpfr_ptr y, mpfr_srcptr x, unsigned long int u, mpfr_rnd_t rnd_mode)
{
MPFR_LOG_FUNC
@@ -31,29 +31,51 @@ mpfr_add_ui (mpfr_ptr y, mpfr_srcptr x, unsigned long int u, mpfr_rnd_t rnd_mode
mpfr_get_prec(x), mpfr_log_prec, x, u, rnd_mode),
("y[%Pu]=%.*Rg", mpfr_get_prec (y), mpfr_log_prec, y));
- if (MPFR_LIKELY(u != 0) ) /* if u=0, do nothing */
+ /* (unsigned long) 0 is assumed to be a real 0 (unsigned) */
+ if (MPFR_UNLIKELY (u == 0))
+ return mpfr_set (y, x, rnd_mode);
+
+ /* This block is actually useless, but this is a minor optimization. */
+ if (MPFR_UNLIKELY (MPFR_IS_SINGULAR (x)))
{
- mpfr_t uu;
- mp_limb_t up[1];
- unsigned long cnt;
- int inex;
- MPFR_SAVE_EXPO_DECL (expo);
-
- MPFR_TMP_INIT1 (up, uu, GMP_NUMB_BITS);
- MPFR_ASSERTD (u == (mp_limb_t) u);
- count_leading_zeros(cnt, (mp_limb_t) u);
- up[0] = (mp_limb_t) u << cnt;
-
- /* Optimization note: Exponent save/restore operations may be
- removed if mpfr_add works even when uu is out-of-range. */
- MPFR_SAVE_EXPO_MARK (expo);
- MPFR_SET_EXP (uu, GMP_NUMB_BITS - cnt);
- inex = mpfr_add(y, x, uu, rnd_mode);
- MPFR_SAVE_EXPO_UPDATE_FLAGS (expo, __gmpfr_flags);
- MPFR_SAVE_EXPO_FREE (expo);
- return mpfr_check_range(y, inex, rnd_mode);
+ if (MPFR_IS_NAN (x))
+ {
+ MPFR_SET_NAN (y);
+ MPFR_RET_NAN;
+ }
+ if (MPFR_IS_INF (x))
+ {
+ MPFR_SET_INF (y);
+ MPFR_SET_SAME_SIGN (y, x);
+ MPFR_RET (0); /* +/-infinity is exact */
+ }
+ MPFR_ASSERTD (MPFR_IS_ZERO (x) && u != 0);
+ /* Note: the fact that u != 0 is important due to signed zeros. */
+ return mpfr_set_ui (y, u, rnd_mode);
}
- else
- /* (unsigned long) 0 is assumed to be a real 0 (unsigned) */
- return mpfr_set (y, x, rnd_mode);
+
+ /* Main code */
+ {
+ mpfr_t uu;
+ mp_limb_t up[1];
+ int cnt;
+ int inex;
+ MPFR_SAVE_EXPO_DECL (expo);
+
+ MPFR_TMP_INIT1 (up, uu, GMP_NUMB_BITS);
+ MPFR_STAT_STATIC_ASSERT (MPFR_LIMB_MAX >= ULONG_MAX);
+ /* So, u fits in a mp_limb_t, which justifies the casts below. */
+ MPFR_ASSERTD (u != 0);
+ count_leading_zeros (cnt, (mp_limb_t) u);
+ up[0] = (mp_limb_t) u << cnt;
+
+ /* Optimization note: Exponent save/restore operations may be
+ removed if mpfr_add works even when uu is out-of-range. */
+ MPFR_SAVE_EXPO_MARK (expo);
+ MPFR_SET_EXP (uu, GMP_NUMB_BITS - cnt);
+ inex = mpfr_add (y, x, uu, rnd_mode);
+ MPFR_SAVE_EXPO_UPDATE_FLAGS (expo, __gmpfr_flags);
+ MPFR_SAVE_EXPO_FREE (expo);
+ return mpfr_check_range (y, inex, rnd_mode);
+ }
}