compiler: add missing Slice_info_expression::do_traverse
[official-gcc.git] / gcc / fold-const-call.cc
blob43819c1f984391e6f18b912548e970b4c04ddacc
1 /* Constant folding for calls to built-in and internal functions.
2 Copyright (C) 1988-2023 Free Software Foundation, Inc.
4 This file is part of GCC.
6 GCC is free software; you can redistribute it and/or modify it under
7 the terms of the GNU General Public License as published by the Free
8 Software Foundation; either version 3, or (at your option) any later
9 version.
11 GCC is distributed in the hope that it will be useful, but WITHOUT ANY
12 WARRANTY; without even the implied warranty of MERCHANTABILITY or
13 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
14 for more details.
16 You should have received a copy of the GNU General Public License
17 along with GCC; see the file COPYING3. If not see
18 <http://www.gnu.org/licenses/>. */
20 #include "config.h"
21 #include "system.h"
22 #include "coretypes.h"
23 #include "realmpfr.h"
24 #include "tree.h"
25 #include "stor-layout.h"
26 #include "options.h"
27 #include "fold-const.h"
28 #include "fold-const-call.h"
29 #include "case-cfn-macros.h"
30 #include "tm.h" /* For C[LT]Z_DEFINED_AT_ZERO. */
31 #include "builtins.h"
32 #include "gimple-expr.h"
33 #include "tree-vector-builder.h"
35 /* Functions that test for certain constant types, abstracting away the
36 decision about whether to check for overflow. */
38 static inline bool
39 integer_cst_p (tree t)
41 return TREE_CODE (t) == INTEGER_CST && !TREE_OVERFLOW (t);
44 static inline bool
45 real_cst_p (tree t)
47 return TREE_CODE (t) == REAL_CST && !TREE_OVERFLOW (t);
50 static inline bool
51 complex_cst_p (tree t)
53 return TREE_CODE (t) == COMPLEX_CST;
56 /* Return true if ARG is a size_type_node constant.
57 Store it in *SIZE_OUT if so. */
59 static inline bool
60 size_t_cst_p (tree t, unsigned HOST_WIDE_INT *size_out)
62 if (types_compatible_p (size_type_node, TREE_TYPE (t))
63 && integer_cst_p (t)
64 && tree_fits_uhwi_p (t))
66 *size_out = tree_to_uhwi (t);
67 return true;
69 return false;
72 /* RES is the result of a comparison in which < 0 means "less", 0 means
73 "equal" and > 0 means "more". Canonicalize it to -1, 0 or 1 and
74 return it in type TYPE. */
76 tree
77 build_cmp_result (tree type, int res)
79 return build_int_cst (type, res < 0 ? -1 : res > 0 ? 1 : 0);
82 /* M is the result of trying to constant-fold an expression (starting
83 with clear MPFR flags) and INEXACT says whether the result in M is
84 exact or inexact. Return true if M can be used as a constant-folded
85 result in format FORMAT, storing the value in *RESULT if so. */
87 static bool
88 do_mpfr_ckconv (real_value *result, mpfr_srcptr m, bool inexact,
89 const real_format *format)
91 /* Proceed iff we get a normal number, i.e. not NaN or Inf and no
92 overflow/underflow occurred. If -frounding-math, proceed iff the
93 result of calling FUNC was exact. */
94 if (!mpfr_number_p (m)
95 || mpfr_overflow_p ()
96 || mpfr_underflow_p ()
97 || (flag_rounding_math && inexact))
98 return false;
100 REAL_VALUE_TYPE tmp;
101 real_from_mpfr (&tmp, m, format, MPFR_RNDN);
103 /* Proceed iff GCC's REAL_VALUE_TYPE can hold the MPFR values.
104 If the REAL_VALUE_TYPE is zero but the mpft_t is not, then we
105 underflowed in the conversion. */
106 if (!real_isfinite (&tmp)
107 || ((tmp.cl == rvc_zero) != (mpfr_zero_p (m) != 0)))
108 return false;
110 real_convert (result, format, &tmp);
111 return real_identical (result, &tmp);
114 /* Try to evaluate:
116 *RESULT = f (*ARG)
118 in format FORMAT, given that FUNC is the MPFR implementation of f.
119 Return true on success. */
121 static bool
122 do_mpfr_arg1 (real_value *result,
123 int (*func) (mpfr_ptr, mpfr_srcptr, mpfr_rnd_t),
124 const real_value *arg, const real_format *format)
126 /* To proceed, MPFR must exactly represent the target floating point
127 format, which only happens when the target base equals two. */
128 if (format->b != 2 || !real_isfinite (arg))
129 return false;
131 int prec = format->p;
132 mpfr_rnd_t rnd = format->round_towards_zero ? MPFR_RNDZ : MPFR_RNDN;
133 mpfr_t m;
135 mpfr_init2 (m, prec);
136 mpfr_from_real (m, arg, MPFR_RNDN);
137 mpfr_clear_flags ();
138 bool inexact = func (m, m, rnd);
139 bool ok = do_mpfr_ckconv (result, m, inexact, format);
140 mpfr_clear (m);
142 return ok;
145 /* Try to evaluate:
147 *RESULT_SIN = sin (*ARG);
148 *RESULT_COS = cos (*ARG);
150 for format FORMAT. Return true on success. */
152 static bool
153 do_mpfr_sincos (real_value *result_sin, real_value *result_cos,
154 const real_value *arg, const real_format *format)
156 /* To proceed, MPFR must exactly represent the target floating point
157 format, which only happens when the target base equals two. */
158 if (format->b != 2 || !real_isfinite (arg))
159 return false;
161 int prec = format->p;
162 mpfr_rnd_t rnd = format->round_towards_zero ? MPFR_RNDZ : MPFR_RNDN;
163 mpfr_t m, ms, mc;
165 mpfr_inits2 (prec, m, ms, mc, NULL);
166 mpfr_from_real (m, arg, MPFR_RNDN);
167 mpfr_clear_flags ();
168 bool inexact = mpfr_sin_cos (ms, mc, m, rnd);
169 bool ok = (do_mpfr_ckconv (result_sin, ms, inexact, format)
170 && do_mpfr_ckconv (result_cos, mc, inexact, format));
171 mpfr_clears (m, ms, mc, NULL);
173 return ok;
176 /* Try to evaluate:
178 *RESULT = f (*ARG0, *ARG1)
180 in format FORMAT, given that FUNC is the MPFR implementation of f.
181 Return true on success. */
183 static bool
184 do_mpfr_arg2 (real_value *result,
185 int (*func) (mpfr_ptr, mpfr_srcptr, mpfr_srcptr, mpfr_rnd_t),
186 const real_value *arg0, const real_value *arg1,
187 const real_format *format)
189 /* To proceed, MPFR must exactly represent the target floating point
190 format, which only happens when the target base equals two. */
191 if (format->b != 2 || !real_isfinite (arg0) || !real_isfinite (arg1))
192 return false;
194 int prec = format->p;
195 mpfr_rnd_t rnd = format->round_towards_zero ? MPFR_RNDZ : MPFR_RNDN;
196 mpfr_t m0, m1;
198 mpfr_inits2 (prec, m0, m1, NULL);
199 mpfr_from_real (m0, arg0, MPFR_RNDN);
200 mpfr_from_real (m1, arg1, MPFR_RNDN);
201 mpfr_clear_flags ();
202 bool inexact = func (m0, m0, m1, rnd);
203 bool ok = do_mpfr_ckconv (result, m0, inexact, format);
204 mpfr_clears (m0, m1, NULL);
206 return ok;
209 /* Try to evaluate:
211 *RESULT = f (ARG0, *ARG1)
213 in format FORMAT, given that FUNC is the MPFR implementation of f.
214 Return true on success. */
216 static bool
217 do_mpfr_arg2 (real_value *result,
218 int (*func) (mpfr_ptr, long, mpfr_srcptr, mpfr_rnd_t),
219 const wide_int_ref &arg0, const real_value *arg1,
220 const real_format *format)
222 if (format->b != 2 || !real_isfinite (arg1))
223 return false;
225 int prec = format->p;
226 mpfr_rnd_t rnd = format->round_towards_zero ? MPFR_RNDZ : MPFR_RNDN;
227 mpfr_t m;
229 mpfr_init2 (m, prec);
230 mpfr_from_real (m, arg1, MPFR_RNDN);
231 mpfr_clear_flags ();
232 bool inexact = func (m, arg0.to_shwi (), m, rnd);
233 bool ok = do_mpfr_ckconv (result, m, inexact, format);
234 mpfr_clear (m);
236 return ok;
239 /* Try to evaluate:
241 *RESULT = f (*ARG0, *ARG1, *ARG2)
243 in format FORMAT, given that FUNC is the MPFR implementation of f.
244 Return true on success. */
246 static bool
247 do_mpfr_arg3 (real_value *result,
248 int (*func) (mpfr_ptr, mpfr_srcptr, mpfr_srcptr,
249 mpfr_srcptr, mpfr_rnd_t),
250 const real_value *arg0, const real_value *arg1,
251 const real_value *arg2, const real_format *format)
253 /* To proceed, MPFR must exactly represent the target floating point
254 format, which only happens when the target base equals two. */
255 if (format->b != 2
256 || !real_isfinite (arg0)
257 || !real_isfinite (arg1)
258 || !real_isfinite (arg2))
259 return false;
261 int prec = format->p;
262 mpfr_rnd_t rnd = format->round_towards_zero ? MPFR_RNDZ : MPFR_RNDN;
263 mpfr_t m0, m1, m2;
265 mpfr_inits2 (prec, m0, m1, m2, NULL);
266 mpfr_from_real (m0, arg0, MPFR_RNDN);
267 mpfr_from_real (m1, arg1, MPFR_RNDN);
268 mpfr_from_real (m2, arg2, MPFR_RNDN);
269 mpfr_clear_flags ();
270 bool inexact = func (m0, m0, m1, m2, rnd);
271 bool ok = do_mpfr_ckconv (result, m0, inexact, format);
272 mpfr_clears (m0, m1, m2, NULL);
274 return ok;
277 /* M is the result of trying to constant-fold an expression (starting
278 with clear MPFR flags) and INEXACT says whether the result in M is
279 exact or inexact. Return true if M can be used as a constant-folded
280 result in which the real and imaginary parts have format FORMAT.
281 Store those parts in *RESULT_REAL and *RESULT_IMAG if so. */
283 static bool
284 do_mpc_ckconv (real_value *result_real, real_value *result_imag,
285 mpc_srcptr m, bool inexact, const real_format *format)
287 /* Proceed iff we get a normal number, i.e. not NaN or Inf and no
288 overflow/underflow occurred. If -frounding-math, proceed iff the
289 result of calling FUNC was exact. */
290 if (!mpfr_number_p (mpc_realref (m))
291 || !mpfr_number_p (mpc_imagref (m))
292 || mpfr_overflow_p ()
293 || mpfr_underflow_p ()
294 || (flag_rounding_math && inexact))
295 return false;
297 REAL_VALUE_TYPE tmp_real, tmp_imag;
298 real_from_mpfr (&tmp_real, mpc_realref (m), format, MPFR_RNDN);
299 real_from_mpfr (&tmp_imag, mpc_imagref (m), format, MPFR_RNDN);
301 /* Proceed iff GCC's REAL_VALUE_TYPE can hold the MPFR values.
302 If the REAL_VALUE_TYPE is zero but the mpft_t is not, then we
303 underflowed in the conversion. */
304 if (!real_isfinite (&tmp_real)
305 || !real_isfinite (&tmp_imag)
306 || (tmp_real.cl == rvc_zero) != (mpfr_zero_p (mpc_realref (m)) != 0)
307 || (tmp_imag.cl == rvc_zero) != (mpfr_zero_p (mpc_imagref (m)) != 0))
308 return false;
310 real_convert (result_real, format, &tmp_real);
311 real_convert (result_imag, format, &tmp_imag);
313 return (real_identical (result_real, &tmp_real)
314 && real_identical (result_imag, &tmp_imag));
317 /* Try to evaluate:
319 RESULT = f (ARG)
321 in format FORMAT, given that FUNC is the mpc implementation of f.
322 Return true on success. Both RESULT and ARG are represented as
323 real and imaginary pairs. */
325 static bool
326 do_mpc_arg1 (real_value *result_real, real_value *result_imag,
327 int (*func) (mpc_ptr, mpc_srcptr, mpc_rnd_t),
328 const real_value *arg_real, const real_value *arg_imag,
329 const real_format *format)
331 /* To proceed, MPFR must exactly represent the target floating point
332 format, which only happens when the target base equals two. */
333 if (format->b != 2
334 || !real_isfinite (arg_real)
335 || !real_isfinite (arg_imag))
336 return false;
338 int prec = format->p;
339 mpc_rnd_t crnd = format->round_towards_zero ? MPC_RNDZZ : MPC_RNDNN;
340 mpc_t m;
342 mpc_init2 (m, prec);
343 mpfr_from_real (mpc_realref (m), arg_real, MPFR_RNDN);
344 mpfr_from_real (mpc_imagref (m), arg_imag, MPFR_RNDN);
345 mpfr_clear_flags ();
346 bool inexact = func (m, m, crnd);
347 bool ok = do_mpc_ckconv (result_real, result_imag, m, inexact, format);
348 mpc_clear (m);
350 return ok;
353 /* Try to evaluate:
355 RESULT = f (ARG0, ARG1)
357 in format FORMAT, given that FUNC is the mpc implementation of f.
358 Return true on success. RESULT, ARG0 and ARG1 are represented as
359 real and imaginary pairs. */
361 static bool
362 do_mpc_arg2 (real_value *result_real, real_value *result_imag,
363 int (*func)(mpc_ptr, mpc_srcptr, mpc_srcptr, mpc_rnd_t),
364 const real_value *arg0_real, const real_value *arg0_imag,
365 const real_value *arg1_real, const real_value *arg1_imag,
366 const real_format *format)
368 if (!real_isfinite (arg0_real)
369 || !real_isfinite (arg0_imag)
370 || !real_isfinite (arg1_real)
371 || !real_isfinite (arg1_imag))
372 return false;
374 int prec = format->p;
375 mpc_rnd_t crnd = format->round_towards_zero ? MPC_RNDZZ : MPC_RNDNN;
376 mpc_t m0, m1;
378 mpc_init2 (m0, prec);
379 mpc_init2 (m1, prec);
380 mpfr_from_real (mpc_realref (m0), arg0_real, MPFR_RNDN);
381 mpfr_from_real (mpc_imagref (m0), arg0_imag, MPFR_RNDN);
382 mpfr_from_real (mpc_realref (m1), arg1_real, MPFR_RNDN);
383 mpfr_from_real (mpc_imagref (m1), arg1_imag, MPFR_RNDN);
384 mpfr_clear_flags ();
385 bool inexact = func (m0, m0, m1, crnd);
386 bool ok = do_mpc_ckconv (result_real, result_imag, m0, inexact, format);
387 mpc_clear (m0);
388 mpc_clear (m1);
390 return ok;
393 /* Try to evaluate:
395 *RESULT = logb (*ARG)
397 in format FORMAT. Return true on success. */
399 static bool
400 fold_const_logb (real_value *result, const real_value *arg,
401 const real_format *format)
403 switch (arg->cl)
405 case rvc_nan:
406 /* If arg is +-NaN, then return it. */
407 *result = *arg;
408 return true;
410 case rvc_inf:
411 /* If arg is +-Inf, then return +Inf. */
412 *result = *arg;
413 result->sign = 0;
414 return true;
416 case rvc_zero:
417 /* Zero may set errno and/or raise an exception. */
418 return false;
420 case rvc_normal:
421 /* For normal numbers, proceed iff radix == 2. In GCC,
422 normalized significands are in the range [0.5, 1.0). We
423 want the exponent as if they were [1.0, 2.0) so get the
424 exponent and subtract 1. */
425 if (format->b == 2)
427 real_from_integer (result, format, REAL_EXP (arg) - 1, SIGNED);
428 return true;
430 return false;
434 /* Try to evaluate:
436 *RESULT = significand (*ARG)
438 in format FORMAT. Return true on success. */
440 static bool
441 fold_const_significand (real_value *result, const real_value *arg,
442 const real_format *format)
444 switch (arg->cl)
446 case rvc_zero:
447 case rvc_nan:
448 case rvc_inf:
449 /* If arg is +-0, +-Inf or +-NaN, then return it. */
450 *result = *arg;
451 return true;
453 case rvc_normal:
454 /* For normal numbers, proceed iff radix == 2. */
455 if (format->b == 2)
457 *result = *arg;
458 /* In GCC, normalized significands are in the range [0.5, 1.0).
459 We want them to be [1.0, 2.0) so set the exponent to 1. */
460 SET_REAL_EXP (result, 1);
461 return true;
463 return false;
467 /* Try to evaluate:
469 *RESULT = f (*ARG)
471 where FORMAT is the format of *ARG and PRECISION is the number of
472 significant bits in the result. Return true on success. */
474 static bool
475 fold_const_conversion (wide_int *result,
476 void (*fn) (real_value *, format_helper,
477 const real_value *),
478 const real_value *arg, unsigned int precision,
479 const real_format *format)
481 if (!real_isfinite (arg))
482 return false;
484 real_value rounded;
485 fn (&rounded, format, arg);
487 bool fail = false;
488 *result = real_to_integer (&rounded, &fail, precision);
489 return !fail;
492 /* Try to evaluate:
494 *RESULT = pow (*ARG0, *ARG1)
496 in format FORMAT. Return true on success. */
498 static bool
499 fold_const_pow (real_value *result, const real_value *arg0,
500 const real_value *arg1, const real_format *format)
502 if (do_mpfr_arg2 (result, mpfr_pow, arg0, arg1, format))
503 return true;
505 /* Check for an integer exponent. */
506 REAL_VALUE_TYPE cint1;
507 HOST_WIDE_INT n1 = real_to_integer (arg1);
508 real_from_integer (&cint1, VOIDmode, n1, SIGNED);
509 /* Attempt to evaluate pow at compile-time, unless this should
510 raise an exception. */
511 if (real_identical (arg1, &cint1)
512 && (n1 > 0
513 || (!flag_trapping_math && !flag_errno_math)
514 || !real_equal (arg0, &dconst0)))
516 bool inexact = real_powi (result, format, arg0, n1);
517 /* Avoid the folding if flag_signaling_nans is on. */
518 if (flag_unsafe_math_optimizations
519 || (!inexact
520 && !(flag_signaling_nans
521 && REAL_VALUE_ISSIGNALING_NAN (*arg0))))
522 return true;
525 return false;
528 /* Try to evaluate:
530 *RESULT = nextafter (*ARG0, *ARG1)
534 *RESULT = nexttoward (*ARG0, *ARG1)
536 in format FORMAT. Return true on success. */
538 static bool
539 fold_const_nextafter (real_value *result, const real_value *arg0,
540 const real_value *arg1, const real_format *format)
542 if (REAL_VALUE_ISSIGNALING_NAN (*arg0)
543 || REAL_VALUE_ISSIGNALING_NAN (*arg1))
544 return false;
546 /* Don't handle composite modes, nor decimal, nor modes without
547 inf or denorm at least for now. */
548 if (format->pnan < format->p
549 || format->b == 10
550 || !format->has_inf
551 || !format->has_denorm)
552 return false;
554 if (real_nextafter (result, format, arg0, arg1)
555 /* If raising underflow or overflow and setting errno to ERANGE,
556 fail if we care about those side-effects. */
557 && (flag_trapping_math || flag_errno_math))
558 return false;
559 /* Similarly for nextafter (0, 1) raising underflow. */
560 else if (flag_trapping_math
561 && arg0->cl == rvc_zero
562 && result->cl != rvc_zero)
563 return false;
565 real_convert (result, format, result);
567 return true;
570 /* Try to evaluate:
572 *RESULT = ldexp (*ARG0, ARG1)
574 in format FORMAT. Return true on success. */
576 static bool
577 fold_const_builtin_load_exponent (real_value *result, const real_value *arg0,
578 const wide_int_ref &arg1,
579 const real_format *format)
581 /* Bound the maximum adjustment to twice the range of the
582 mode's valid exponents. Use abs to ensure the range is
583 positive as a sanity check. */
584 int max_exp_adj = 2 * labs (format->emax - format->emin);
586 /* The requested adjustment must be inside this range. This
587 is a preliminary cap to avoid things like overflow, we
588 may still fail to compute the result for other reasons. */
589 if (wi::les_p (arg1, -max_exp_adj) || wi::ges_p (arg1, max_exp_adj))
590 return false;
592 /* Don't perform operation if we honor signaling NaNs and
593 operand is a signaling NaN. */
594 if (!flag_unsafe_math_optimizations
595 && flag_signaling_nans
596 && REAL_VALUE_ISSIGNALING_NAN (*arg0))
597 return false;
599 REAL_VALUE_TYPE initial_result;
600 real_ldexp (&initial_result, arg0, arg1.to_shwi ());
602 /* Ensure we didn't overflow. */
603 if (real_isinf (&initial_result))
604 return false;
606 /* Only proceed if the target mode can hold the
607 resulting value. */
608 *result = real_value_truncate (format, initial_result);
609 return real_equal (&initial_result, result);
612 /* Fold a call to __builtin_nan or __builtin_nans with argument ARG and
613 return type TYPE. QUIET is true if a quiet rather than signalling
614 NaN is required. */
616 static tree
617 fold_const_builtin_nan (tree type, tree arg, bool quiet)
619 REAL_VALUE_TYPE real;
620 const char *str = c_getstr (arg);
621 if (str && real_nan (&real, str, quiet, TYPE_MODE (type)))
622 return build_real (type, real);
623 return NULL_TREE;
626 /* Fold a call to IFN_REDUC_<CODE> (ARG), returning a value of type TYPE. */
628 static tree
629 fold_const_reduction (tree type, tree arg, tree_code code)
631 unsigned HOST_WIDE_INT nelts;
632 if (TREE_CODE (arg) != VECTOR_CST
633 || !VECTOR_CST_NELTS (arg).is_constant (&nelts))
634 return NULL_TREE;
636 tree res = VECTOR_CST_ELT (arg, 0);
637 for (unsigned HOST_WIDE_INT i = 1; i < nelts; i++)
639 res = const_binop (code, type, res, VECTOR_CST_ELT (arg, i));
640 if (res == NULL_TREE || !CONSTANT_CLASS_P (res))
641 return NULL_TREE;
643 return res;
646 /* Fold a call to IFN_VEC_CONVERT (ARG) returning TYPE. */
648 static tree
649 fold_const_vec_convert (tree ret_type, tree arg)
651 enum tree_code code = NOP_EXPR;
652 tree arg_type = TREE_TYPE (arg);
653 if (TREE_CODE (arg) != VECTOR_CST)
654 return NULL_TREE;
656 gcc_checking_assert (VECTOR_TYPE_P (ret_type) && VECTOR_TYPE_P (arg_type));
658 if (INTEGRAL_TYPE_P (TREE_TYPE (ret_type))
659 && SCALAR_FLOAT_TYPE_P (TREE_TYPE (arg_type)))
660 code = FIX_TRUNC_EXPR;
661 else if (INTEGRAL_TYPE_P (TREE_TYPE (arg_type))
662 && SCALAR_FLOAT_TYPE_P (TREE_TYPE (ret_type)))
663 code = FLOAT_EXPR;
665 /* We can't handle steps directly when extending, since the
666 values need to wrap at the original precision first. */
667 bool step_ok_p
668 = (INTEGRAL_TYPE_P (TREE_TYPE (ret_type))
669 && INTEGRAL_TYPE_P (TREE_TYPE (arg_type))
670 && (TYPE_PRECISION (TREE_TYPE (ret_type))
671 <= TYPE_PRECISION (TREE_TYPE (arg_type))));
672 tree_vector_builder elts;
673 if (!elts.new_unary_operation (ret_type, arg, step_ok_p))
674 return NULL_TREE;
676 unsigned int count = elts.encoded_nelts ();
677 for (unsigned int i = 0; i < count; ++i)
679 tree elt = fold_unary (code, TREE_TYPE (ret_type),
680 VECTOR_CST_ELT (arg, i));
681 if (elt == NULL_TREE || !CONSTANT_CLASS_P (elt))
682 return NULL_TREE;
683 elts.quick_push (elt);
686 return elts.build ();
689 /* Try to evaluate:
691 IFN_WHILE_ULT (ARG0, ARG1, (TYPE) { ... })
693 Return the value on success and null on failure. */
695 static tree
696 fold_while_ult (tree type, poly_uint64 arg0, poly_uint64 arg1)
698 if (known_ge (arg0, arg1))
699 return build_zero_cst (type);
701 if (maybe_ge (arg0, arg1))
702 return NULL_TREE;
704 poly_uint64 diff = arg1 - arg0;
705 poly_uint64 nelts = TYPE_VECTOR_SUBPARTS (type);
706 if (known_ge (diff, nelts))
707 return build_all_ones_cst (type);
709 unsigned HOST_WIDE_INT const_diff;
710 if (known_le (diff, nelts) && diff.is_constant (&const_diff))
712 tree minus_one = build_minus_one_cst (TREE_TYPE (type));
713 tree zero = build_zero_cst (TREE_TYPE (type));
714 return build_vector_a_then_b (type, const_diff, minus_one, zero);
716 return NULL_TREE;
719 /* Try to evaluate:
721 *RESULT = FN (*ARG)
723 in format FORMAT. Return true on success. */
725 static bool
726 fold_const_call_ss (real_value *result, combined_fn fn,
727 const real_value *arg, const real_format *format)
729 switch (fn)
731 CASE_CFN_SQRT:
732 CASE_CFN_SQRT_FN:
733 return (real_compare (GE_EXPR, arg, &dconst0)
734 && do_mpfr_arg1 (result, mpfr_sqrt, arg, format));
736 CASE_CFN_CBRT:
737 CASE_CFN_CBRT_FN:
738 return do_mpfr_arg1 (result, mpfr_cbrt, arg, format);
740 CASE_CFN_ASIN:
741 CASE_CFN_ASIN_FN:
742 return (real_compare (GE_EXPR, arg, &dconstm1)
743 && real_compare (LE_EXPR, arg, &dconst1)
744 && do_mpfr_arg1 (result, mpfr_asin, arg, format));
746 CASE_CFN_ACOS:
747 CASE_CFN_ACOS_FN:
748 return (real_compare (GE_EXPR, arg, &dconstm1)
749 && real_compare (LE_EXPR, arg, &dconst1)
750 && do_mpfr_arg1 (result, mpfr_acos, arg, format));
752 CASE_CFN_ATAN:
753 CASE_CFN_ATAN_FN:
754 return do_mpfr_arg1 (result, mpfr_atan, arg, format);
756 CASE_CFN_ASINH:
757 CASE_CFN_ASINH_FN:
758 return do_mpfr_arg1 (result, mpfr_asinh, arg, format);
760 CASE_CFN_ACOSH:
761 CASE_CFN_ACOSH_FN:
762 return (real_compare (GE_EXPR, arg, &dconst1)
763 && do_mpfr_arg1 (result, mpfr_acosh, arg, format));
765 CASE_CFN_ATANH:
766 CASE_CFN_ATANH_FN:
767 return (real_compare (GE_EXPR, arg, &dconstm1)
768 && real_compare (LE_EXPR, arg, &dconst1)
769 && do_mpfr_arg1 (result, mpfr_atanh, arg, format));
771 CASE_CFN_SIN:
772 CASE_CFN_SIN_FN:
773 return do_mpfr_arg1 (result, mpfr_sin, arg, format);
775 CASE_CFN_COS:
776 CASE_CFN_COS_FN:
777 return do_mpfr_arg1 (result, mpfr_cos, arg, format);
779 CASE_CFN_TAN:
780 CASE_CFN_TAN_FN:
781 return do_mpfr_arg1 (result, mpfr_tan, arg, format);
783 CASE_CFN_SINH:
784 CASE_CFN_SINH_FN:
785 return do_mpfr_arg1 (result, mpfr_sinh, arg, format);
787 CASE_CFN_COSH:
788 CASE_CFN_COSH_FN:
789 return do_mpfr_arg1 (result, mpfr_cosh, arg, format);
791 CASE_CFN_TANH:
792 CASE_CFN_TANH_FN:
793 return do_mpfr_arg1 (result, mpfr_tanh, arg, format);
795 CASE_CFN_ERF:
796 CASE_CFN_ERF_FN:
797 return do_mpfr_arg1 (result, mpfr_erf, arg, format);
799 CASE_CFN_ERFC:
800 CASE_CFN_ERFC_FN:
801 return do_mpfr_arg1 (result, mpfr_erfc, arg, format);
803 CASE_CFN_TGAMMA:
804 CASE_CFN_TGAMMA_FN:
805 return do_mpfr_arg1 (result, mpfr_gamma, arg, format);
807 CASE_CFN_EXP:
808 CASE_CFN_EXP_FN:
809 return do_mpfr_arg1 (result, mpfr_exp, arg, format);
811 CASE_CFN_EXP2:
812 CASE_CFN_EXP2_FN:
813 return do_mpfr_arg1 (result, mpfr_exp2, arg, format);
815 CASE_CFN_EXP10:
816 CASE_CFN_POW10:
817 return do_mpfr_arg1 (result, mpfr_exp10, arg, format);
819 CASE_CFN_EXPM1:
820 CASE_CFN_EXPM1_FN:
821 return do_mpfr_arg1 (result, mpfr_expm1, arg, format);
823 CASE_CFN_LOG:
824 CASE_CFN_LOG_FN:
825 return (real_compare (GT_EXPR, arg, &dconst0)
826 && do_mpfr_arg1 (result, mpfr_log, arg, format));
828 CASE_CFN_LOG2:
829 CASE_CFN_LOG2_FN:
830 return (real_compare (GT_EXPR, arg, &dconst0)
831 && do_mpfr_arg1 (result, mpfr_log2, arg, format));
833 CASE_CFN_LOG10:
834 CASE_CFN_LOG10_FN:
835 return (real_compare (GT_EXPR, arg, &dconst0)
836 && do_mpfr_arg1 (result, mpfr_log10, arg, format));
838 CASE_CFN_LOG1P:
839 CASE_CFN_LOG1P_FN:
840 return (real_compare (GT_EXPR, arg, &dconstm1)
841 && do_mpfr_arg1 (result, mpfr_log1p, arg, format));
843 CASE_CFN_J0:
844 return do_mpfr_arg1 (result, mpfr_j0, arg, format);
846 CASE_CFN_J1:
847 return do_mpfr_arg1 (result, mpfr_j1, arg, format);
849 CASE_CFN_Y0:
850 return (real_compare (GT_EXPR, arg, &dconst0)
851 && do_mpfr_arg1 (result, mpfr_y0, arg, format));
853 CASE_CFN_Y1:
854 return (real_compare (GT_EXPR, arg, &dconst0)
855 && do_mpfr_arg1 (result, mpfr_y1, arg, format));
857 CASE_CFN_FLOOR:
858 CASE_CFN_FLOOR_FN:
859 if (!REAL_VALUE_ISSIGNALING_NAN (*arg))
861 real_floor (result, format, arg);
862 return true;
864 return false;
866 CASE_CFN_CEIL:
867 CASE_CFN_CEIL_FN:
868 if (!REAL_VALUE_ISSIGNALING_NAN (*arg))
870 real_ceil (result, format, arg);
871 return true;
873 return false;
875 CASE_CFN_TRUNC:
876 CASE_CFN_TRUNC_FN:
877 if (!REAL_VALUE_ISSIGNALING_NAN (*arg))
879 real_trunc (result, format, arg);
880 return true;
882 return false;
884 CASE_CFN_ROUND:
885 CASE_CFN_ROUND_FN:
886 if (!REAL_VALUE_ISSIGNALING_NAN (*arg))
888 real_round (result, format, arg);
889 return true;
891 return false;
893 CASE_CFN_ROUNDEVEN:
894 CASE_CFN_ROUNDEVEN_FN:
895 if (!REAL_VALUE_ISSIGNALING_NAN (*arg))
897 real_roundeven (result, format, arg);
898 return true;
900 return false;
902 CASE_CFN_LOGB:
903 CASE_CFN_LOGB_FN:
904 return fold_const_logb (result, arg, format);
906 CASE_CFN_SIGNIFICAND:
907 return fold_const_significand (result, arg, format);
909 default:
910 return false;
914 /* Try to evaluate:
916 *RESULT = FN (*ARG)
918 where FORMAT is the format of ARG and PRECISION is the number of
919 significant bits in the result. Return true on success. */
921 static bool
922 fold_const_call_ss (wide_int *result, combined_fn fn,
923 const real_value *arg, unsigned int precision,
924 const real_format *format)
926 switch (fn)
928 CASE_CFN_SIGNBIT:
929 if (real_isneg (arg))
930 *result = wi::one (precision);
931 else
932 *result = wi::zero (precision);
933 return true;
935 CASE_CFN_ILOGB:
936 CASE_CFN_ILOGB_FN:
937 /* For ilogb we don't know FP_ILOGB0, so only handle normal values.
938 Proceed iff radix == 2. In GCC, normalized significands are in
939 the range [0.5, 1.0). We want the exponent as if they were
940 [1.0, 2.0) so get the exponent and subtract 1. */
941 if (arg->cl == rvc_normal && format->b == 2)
943 *result = wi::shwi (REAL_EXP (arg) - 1, precision);
944 return true;
946 return false;
948 CASE_CFN_ICEIL:
949 CASE_CFN_LCEIL:
950 CASE_CFN_LLCEIL:
951 return fold_const_conversion (result, real_ceil, arg,
952 precision, format);
954 CASE_CFN_LFLOOR:
955 CASE_CFN_IFLOOR:
956 CASE_CFN_LLFLOOR:
957 return fold_const_conversion (result, real_floor, arg,
958 precision, format);
960 CASE_CFN_IROUND:
961 CASE_CFN_LROUND:
962 CASE_CFN_LROUND_FN:
963 CASE_CFN_LLROUND:
964 CASE_CFN_LLROUND_FN:
965 return fold_const_conversion (result, real_round, arg,
966 precision, format);
968 CASE_CFN_IRINT:
969 CASE_CFN_LRINT:
970 CASE_CFN_LRINT_FN:
971 CASE_CFN_LLRINT:
972 CASE_CFN_LLRINT_FN:
973 /* Not yet folded to a constant. */
974 return false;
976 CASE_CFN_FINITE:
977 case CFN_BUILT_IN_FINITED32:
978 case CFN_BUILT_IN_FINITED64:
979 case CFN_BUILT_IN_FINITED128:
980 case CFN_BUILT_IN_ISFINITE:
981 *result = wi::shwi (real_isfinite (arg) ? 1 : 0, precision);
982 return true;
984 case CFN_BUILT_IN_ISSIGNALING:
985 *result = wi::shwi (real_issignaling_nan (arg) ? 1 : 0, precision);
986 return true;
988 CASE_CFN_ISINF:
989 case CFN_BUILT_IN_ISINFD32:
990 case CFN_BUILT_IN_ISINFD64:
991 case CFN_BUILT_IN_ISINFD128:
992 if (real_isinf (arg))
993 *result = wi::shwi (arg->sign ? -1 : 1, precision);
994 else
995 *result = wi::shwi (0, precision);
996 return true;
998 CASE_CFN_ISNAN:
999 case CFN_BUILT_IN_ISNAND32:
1000 case CFN_BUILT_IN_ISNAND64:
1001 case CFN_BUILT_IN_ISNAND128:
1002 *result = wi::shwi (real_isnan (arg) ? 1 : 0, precision);
1003 return true;
1005 default:
1006 return false;
1010 /* Try to evaluate:
1012 *RESULT = FN (ARG)
1014 where ARG_TYPE is the type of ARG and PRECISION is the number of bits
1015 in the result. Return true on success. */
1017 static bool
1018 fold_const_call_ss (wide_int *result, combined_fn fn, const wide_int_ref &arg,
1019 unsigned int precision, tree arg_type)
1021 switch (fn)
1023 CASE_CFN_FFS:
1024 *result = wi::shwi (wi::ffs (arg), precision);
1025 return true;
1027 CASE_CFN_CLZ:
1029 int tmp;
1030 if (wi::ne_p (arg, 0))
1031 tmp = wi::clz (arg);
1032 else if (!CLZ_DEFINED_VALUE_AT_ZERO (SCALAR_INT_TYPE_MODE (arg_type),
1033 tmp))
1034 tmp = TYPE_PRECISION (arg_type);
1035 *result = wi::shwi (tmp, precision);
1036 return true;
1039 CASE_CFN_CTZ:
1041 int tmp;
1042 if (wi::ne_p (arg, 0))
1043 tmp = wi::ctz (arg);
1044 else if (!CTZ_DEFINED_VALUE_AT_ZERO (SCALAR_INT_TYPE_MODE (arg_type),
1045 tmp))
1046 tmp = TYPE_PRECISION (arg_type);
1047 *result = wi::shwi (tmp, precision);
1048 return true;
1051 CASE_CFN_CLRSB:
1052 *result = wi::shwi (wi::clrsb (arg), precision);
1053 return true;
1055 CASE_CFN_POPCOUNT:
1056 *result = wi::shwi (wi::popcount (arg), precision);
1057 return true;
1059 CASE_CFN_PARITY:
1060 *result = wi::shwi (wi::parity (arg), precision);
1061 return true;
1063 case CFN_BUILT_IN_BSWAP16:
1064 case CFN_BUILT_IN_BSWAP32:
1065 case CFN_BUILT_IN_BSWAP64:
1066 case CFN_BUILT_IN_BSWAP128:
1067 *result = wide_int::from (arg, precision, TYPE_SIGN (arg_type)).bswap ();
1068 return true;
1070 default:
1071 return false;
1075 /* Try to evaluate:
1077 RESULT = FN (*ARG)
1079 where FORMAT is the format of ARG and of the real and imaginary parts
1080 of RESULT, passed as RESULT_REAL and RESULT_IMAG respectively. Return
1081 true on success. */
1083 static bool
1084 fold_const_call_cs (real_value *result_real, real_value *result_imag,
1085 combined_fn fn, const real_value *arg,
1086 const real_format *format)
1088 switch (fn)
1090 CASE_CFN_CEXPI:
1091 /* cexpi(x+yi) = cos(x)+sin(y)*i. */
1092 return do_mpfr_sincos (result_imag, result_real, arg, format);
1094 default:
1095 return false;
1099 /* Try to evaluate:
1101 *RESULT = fn (ARG)
1103 where FORMAT is the format of RESULT and of the real and imaginary parts
1104 of ARG, passed as ARG_REAL and ARG_IMAG respectively. Return true on
1105 success. */
1107 static bool
1108 fold_const_call_sc (real_value *result, combined_fn fn,
1109 const real_value *arg_real, const real_value *arg_imag,
1110 const real_format *format)
1112 switch (fn)
1114 CASE_CFN_CABS:
1115 CASE_CFN_CABS_FN:
1116 return do_mpfr_arg2 (result, mpfr_hypot, arg_real, arg_imag, format);
1118 default:
1119 return false;
1123 /* Try to evaluate:
1125 RESULT = fn (ARG)
1127 where FORMAT is the format of the real and imaginary parts of RESULT
1128 (RESULT_REAL and RESULT_IMAG) and of ARG (ARG_REAL and ARG_IMAG).
1129 Return true on success. */
1131 static bool
1132 fold_const_call_cc (real_value *result_real, real_value *result_imag,
1133 combined_fn fn, const real_value *arg_real,
1134 const real_value *arg_imag, const real_format *format)
1136 switch (fn)
1138 CASE_CFN_CCOS:
1139 CASE_CFN_CCOS_FN:
1140 return do_mpc_arg1 (result_real, result_imag, mpc_cos,
1141 arg_real, arg_imag, format);
1143 CASE_CFN_CCOSH:
1144 CASE_CFN_CCOSH_FN:
1145 return do_mpc_arg1 (result_real, result_imag, mpc_cosh,
1146 arg_real, arg_imag, format);
1148 CASE_CFN_CPROJ:
1149 CASE_CFN_CPROJ_FN:
1150 if (real_isinf (arg_real) || real_isinf (arg_imag))
1152 *result_real = dconstinf;
1153 *result_imag = dconst0;
1154 result_imag->sign = arg_imag->sign;
1156 else
1158 *result_real = *arg_real;
1159 *result_imag = *arg_imag;
1161 return true;
1163 CASE_CFN_CSIN:
1164 CASE_CFN_CSIN_FN:
1165 return do_mpc_arg1 (result_real, result_imag, mpc_sin,
1166 arg_real, arg_imag, format);
1168 CASE_CFN_CSINH:
1169 CASE_CFN_CSINH_FN:
1170 return do_mpc_arg1 (result_real, result_imag, mpc_sinh,
1171 arg_real, arg_imag, format);
1173 CASE_CFN_CTAN:
1174 CASE_CFN_CTAN_FN:
1175 return do_mpc_arg1 (result_real, result_imag, mpc_tan,
1176 arg_real, arg_imag, format);
1178 CASE_CFN_CTANH:
1179 CASE_CFN_CTANH_FN:
1180 return do_mpc_arg1 (result_real, result_imag, mpc_tanh,
1181 arg_real, arg_imag, format);
1183 CASE_CFN_CLOG:
1184 CASE_CFN_CLOG_FN:
1185 return do_mpc_arg1 (result_real, result_imag, mpc_log,
1186 arg_real, arg_imag, format);
1188 CASE_CFN_CSQRT:
1189 CASE_CFN_CSQRT_FN:
1190 return do_mpc_arg1 (result_real, result_imag, mpc_sqrt,
1191 arg_real, arg_imag, format);
1193 CASE_CFN_CASIN:
1194 CASE_CFN_CASIN_FN:
1195 return do_mpc_arg1 (result_real, result_imag, mpc_asin,
1196 arg_real, arg_imag, format);
1198 CASE_CFN_CACOS:
1199 CASE_CFN_CACOS_FN:
1200 return do_mpc_arg1 (result_real, result_imag, mpc_acos,
1201 arg_real, arg_imag, format);
1203 CASE_CFN_CATAN:
1204 CASE_CFN_CATAN_FN:
1205 return do_mpc_arg1 (result_real, result_imag, mpc_atan,
1206 arg_real, arg_imag, format);
1208 CASE_CFN_CASINH:
1209 CASE_CFN_CASINH_FN:
1210 return do_mpc_arg1 (result_real, result_imag, mpc_asinh,
1211 arg_real, arg_imag, format);
1213 CASE_CFN_CACOSH:
1214 CASE_CFN_CACOSH_FN:
1215 return do_mpc_arg1 (result_real, result_imag, mpc_acosh,
1216 arg_real, arg_imag, format);
1218 CASE_CFN_CATANH:
1219 CASE_CFN_CATANH_FN:
1220 return do_mpc_arg1 (result_real, result_imag, mpc_atanh,
1221 arg_real, arg_imag, format);
1223 CASE_CFN_CEXP:
1224 CASE_CFN_CEXP_FN:
1225 return do_mpc_arg1 (result_real, result_imag, mpc_exp,
1226 arg_real, arg_imag, format);
1228 default:
1229 return false;
1233 /* Subroutine of fold_const_call, with the same interface. Handle cases
1234 where the arguments and result are numerical. */
1236 static tree
1237 fold_const_call_1 (combined_fn fn, tree type, tree arg)
1239 machine_mode mode = TYPE_MODE (type);
1240 machine_mode arg_mode = TYPE_MODE (TREE_TYPE (arg));
1242 if (integer_cst_p (arg))
1244 if (SCALAR_INT_MODE_P (mode))
1246 wide_int result;
1247 if (fold_const_call_ss (&result, fn, wi::to_wide (arg),
1248 TYPE_PRECISION (type), TREE_TYPE (arg)))
1249 return wide_int_to_tree (type, result);
1251 return NULL_TREE;
1254 if (real_cst_p (arg))
1256 gcc_checking_assert (SCALAR_FLOAT_MODE_P (arg_mode));
1257 if (mode == arg_mode)
1259 /* real -> real. */
1260 REAL_VALUE_TYPE result;
1261 if (fold_const_call_ss (&result, fn, TREE_REAL_CST_PTR (arg),
1262 REAL_MODE_FORMAT (mode)))
1263 return build_real (type, result);
1265 else if (COMPLEX_MODE_P (mode)
1266 && GET_MODE_INNER (mode) == arg_mode)
1268 /* real -> complex real. */
1269 REAL_VALUE_TYPE result_real, result_imag;
1270 if (fold_const_call_cs (&result_real, &result_imag, fn,
1271 TREE_REAL_CST_PTR (arg),
1272 REAL_MODE_FORMAT (arg_mode)))
1273 return build_complex (type,
1274 build_real (TREE_TYPE (type), result_real),
1275 build_real (TREE_TYPE (type), result_imag));
1277 else if (INTEGRAL_TYPE_P (type))
1279 /* real -> int. */
1280 wide_int result;
1281 if (fold_const_call_ss (&result, fn,
1282 TREE_REAL_CST_PTR (arg),
1283 TYPE_PRECISION (type),
1284 REAL_MODE_FORMAT (arg_mode)))
1285 return wide_int_to_tree (type, result);
1287 return NULL_TREE;
1290 if (complex_cst_p (arg))
1292 gcc_checking_assert (COMPLEX_MODE_P (arg_mode));
1293 machine_mode inner_mode = GET_MODE_INNER (arg_mode);
1294 tree argr = TREE_REALPART (arg);
1295 tree argi = TREE_IMAGPART (arg);
1296 if (mode == arg_mode
1297 && real_cst_p (argr)
1298 && real_cst_p (argi))
1300 /* complex real -> complex real. */
1301 REAL_VALUE_TYPE result_real, result_imag;
1302 if (fold_const_call_cc (&result_real, &result_imag, fn,
1303 TREE_REAL_CST_PTR (argr),
1304 TREE_REAL_CST_PTR (argi),
1305 REAL_MODE_FORMAT (inner_mode)))
1306 return build_complex (type,
1307 build_real (TREE_TYPE (type), result_real),
1308 build_real (TREE_TYPE (type), result_imag));
1310 if (mode == inner_mode
1311 && real_cst_p (argr)
1312 && real_cst_p (argi))
1314 /* complex real -> real. */
1315 REAL_VALUE_TYPE result;
1316 if (fold_const_call_sc (&result, fn,
1317 TREE_REAL_CST_PTR (argr),
1318 TREE_REAL_CST_PTR (argi),
1319 REAL_MODE_FORMAT (inner_mode)))
1320 return build_real (type, result);
1322 return NULL_TREE;
1325 return NULL_TREE;
1328 /* Try to fold FN (ARG) to a constant. Return the constant on success,
1329 otherwise return null. TYPE is the type of the return value. */
1331 tree
1332 fold_const_call (combined_fn fn, tree type, tree arg)
1334 switch (fn)
1336 case CFN_BUILT_IN_STRLEN:
1337 if (const char *str = c_getstr (arg))
1338 return build_int_cst (type, strlen (str));
1339 return NULL_TREE;
1341 CASE_CFN_NAN:
1342 CASE_FLT_FN_FLOATN_NX (CFN_BUILT_IN_NAN):
1343 case CFN_BUILT_IN_NAND32:
1344 case CFN_BUILT_IN_NAND64:
1345 case CFN_BUILT_IN_NAND128:
1346 return fold_const_builtin_nan (type, arg, true);
1348 CASE_CFN_NANS:
1349 CASE_FLT_FN_FLOATN_NX (CFN_BUILT_IN_NANS):
1350 case CFN_BUILT_IN_NANSF16B:
1351 case CFN_BUILT_IN_NANSD32:
1352 case CFN_BUILT_IN_NANSD64:
1353 case CFN_BUILT_IN_NANSD128:
1354 return fold_const_builtin_nan (type, arg, false);
1356 case CFN_REDUC_PLUS:
1357 return fold_const_reduction (type, arg, PLUS_EXPR);
1359 case CFN_REDUC_MAX:
1360 return fold_const_reduction (type, arg, MAX_EXPR);
1362 case CFN_REDUC_MIN:
1363 return fold_const_reduction (type, arg, MIN_EXPR);
1365 case CFN_REDUC_AND:
1366 return fold_const_reduction (type, arg, BIT_AND_EXPR);
1368 case CFN_REDUC_IOR:
1369 return fold_const_reduction (type, arg, BIT_IOR_EXPR);
1371 case CFN_REDUC_XOR:
1372 return fold_const_reduction (type, arg, BIT_XOR_EXPR);
1374 case CFN_VEC_CONVERT:
1375 return fold_const_vec_convert (type, arg);
1377 default:
1378 return fold_const_call_1 (fn, type, arg);
1382 /* Fold a call to IFN_FOLD_LEFT_<CODE> (ARG0, ARG1), returning a value
1383 of type TYPE. */
1385 static tree
1386 fold_const_fold_left (tree type, tree arg0, tree arg1, tree_code code)
1388 if (TREE_CODE (arg1) != VECTOR_CST)
1389 return NULL_TREE;
1391 unsigned HOST_WIDE_INT nelts;
1392 if (!VECTOR_CST_NELTS (arg1).is_constant (&nelts))
1393 return NULL_TREE;
1395 for (unsigned HOST_WIDE_INT i = 0; i < nelts; i++)
1397 arg0 = const_binop (code, type, arg0, VECTOR_CST_ELT (arg1, i));
1398 if (arg0 == NULL_TREE || !CONSTANT_CLASS_P (arg0))
1399 return NULL_TREE;
1401 return arg0;
1404 /* Try to evaluate:
1406 *RESULT = FN (*ARG0, *ARG1)
1408 in format FORMAT. Return true on success. */
1410 static bool
1411 fold_const_call_sss (real_value *result, combined_fn fn,
1412 const real_value *arg0, const real_value *arg1,
1413 const real_format *format)
1415 switch (fn)
1417 CASE_CFN_DREM:
1418 CASE_CFN_REMAINDER:
1419 CASE_CFN_REMAINDER_FN:
1420 return do_mpfr_arg2 (result, mpfr_remainder, arg0, arg1, format);
1422 CASE_CFN_ATAN2:
1423 CASE_CFN_ATAN2_FN:
1424 return do_mpfr_arg2 (result, mpfr_atan2, arg0, arg1, format);
1426 CASE_CFN_FDIM:
1427 CASE_CFN_FDIM_FN:
1428 return do_mpfr_arg2 (result, mpfr_dim, arg0, arg1, format);
1430 CASE_CFN_FMOD:
1431 CASE_CFN_FMOD_FN:
1432 return do_mpfr_arg2 (result, mpfr_fmod, arg0, arg1, format);
1434 CASE_CFN_HYPOT:
1435 CASE_CFN_HYPOT_FN:
1436 return do_mpfr_arg2 (result, mpfr_hypot, arg0, arg1, format);
1438 CASE_CFN_COPYSIGN:
1439 CASE_CFN_COPYSIGN_FN:
1440 *result = *arg0;
1441 real_copysign (result, arg1);
1442 return true;
1444 CASE_CFN_FMIN:
1445 CASE_CFN_FMIN_FN:
1446 return do_mpfr_arg2 (result, mpfr_min, arg0, arg1, format);
1448 CASE_CFN_FMAX:
1449 CASE_CFN_FMAX_FN:
1450 return do_mpfr_arg2 (result, mpfr_max, arg0, arg1, format);
1452 CASE_CFN_POW:
1453 CASE_CFN_POW_FN:
1454 return fold_const_pow (result, arg0, arg1, format);
1456 CASE_CFN_NEXTAFTER:
1457 CASE_CFN_NEXTAFTER_FN:
1458 case CFN_BUILT_IN_NEXTAFTERF16B:
1459 CASE_CFN_NEXTTOWARD:
1460 return fold_const_nextafter (result, arg0, arg1, format);
1462 default:
1463 return false;
1467 /* Try to evaluate:
1469 *RESULT = FN (*ARG0, ARG1)
1471 where FORMAT is the format of *RESULT and *ARG0. Return true on
1472 success. */
1474 static bool
1475 fold_const_call_sss (real_value *result, combined_fn fn,
1476 const real_value *arg0, const wide_int_ref &arg1,
1477 const real_format *format)
1479 switch (fn)
1481 CASE_CFN_LDEXP:
1482 CASE_CFN_LDEXP_FN:
1483 return fold_const_builtin_load_exponent (result, arg0, arg1, format);
1485 CASE_CFN_SCALBN:
1486 CASE_CFN_SCALBN_FN:
1487 CASE_CFN_SCALBLN:
1488 CASE_CFN_SCALBLN_FN:
1489 return (format->b == 2
1490 && fold_const_builtin_load_exponent (result, arg0, arg1,
1491 format));
1493 CASE_CFN_POWI:
1494 /* Avoid the folding if flag_signaling_nans is on and
1495 operand is a signaling NaN. */
1496 if (!flag_unsafe_math_optimizations
1497 && flag_signaling_nans
1498 && REAL_VALUE_ISSIGNALING_NAN (*arg0))
1499 return false;
1501 real_powi (result, format, arg0, arg1.to_shwi ());
1502 return true;
1504 default:
1505 return false;
1509 /* Try to evaluate:
1511 *RESULT = FN (ARG0, *ARG1)
1513 where FORMAT is the format of *RESULT and *ARG1. Return true on
1514 success. */
1516 static bool
1517 fold_const_call_sss (real_value *result, combined_fn fn,
1518 const wide_int_ref &arg0, const real_value *arg1,
1519 const real_format *format)
1521 switch (fn)
1523 CASE_CFN_JN:
1524 return do_mpfr_arg2 (result, mpfr_jn, arg0, arg1, format);
1526 CASE_CFN_YN:
1527 return (real_compare (GT_EXPR, arg1, &dconst0)
1528 && do_mpfr_arg2 (result, mpfr_yn, arg0, arg1, format));
1530 default:
1531 return false;
1535 /* Try to evaluate:
1537 RESULT = fn (ARG0, ARG1)
1539 where FORMAT is the format of the real and imaginary parts of RESULT
1540 (RESULT_REAL and RESULT_IMAG), of ARG0 (ARG0_REAL and ARG0_IMAG)
1541 and of ARG1 (ARG1_REAL and ARG1_IMAG). Return true on success. */
1543 static bool
1544 fold_const_call_ccc (real_value *result_real, real_value *result_imag,
1545 combined_fn fn, const real_value *arg0_real,
1546 const real_value *arg0_imag, const real_value *arg1_real,
1547 const real_value *arg1_imag, const real_format *format)
1549 switch (fn)
1551 CASE_CFN_CPOW:
1552 CASE_CFN_CPOW_FN:
1553 return do_mpc_arg2 (result_real, result_imag, mpc_pow,
1554 arg0_real, arg0_imag, arg1_real, arg1_imag, format);
1556 default:
1557 return false;
1561 /* Subroutine of fold_const_call, with the same interface. Handle cases
1562 where the arguments and result are numerical. */
1564 static tree
1565 fold_const_call_1 (combined_fn fn, tree type, tree arg0, tree arg1)
1567 machine_mode mode = TYPE_MODE (type);
1568 machine_mode arg0_mode = TYPE_MODE (TREE_TYPE (arg0));
1569 machine_mode arg1_mode = TYPE_MODE (TREE_TYPE (arg1));
1571 if (mode == arg0_mode
1572 && real_cst_p (arg0)
1573 && real_cst_p (arg1))
1575 gcc_checking_assert (SCALAR_FLOAT_MODE_P (arg0_mode));
1576 REAL_VALUE_TYPE result;
1577 if (arg0_mode == arg1_mode)
1579 /* real, real -> real. */
1580 if (fold_const_call_sss (&result, fn, TREE_REAL_CST_PTR (arg0),
1581 TREE_REAL_CST_PTR (arg1),
1582 REAL_MODE_FORMAT (mode)))
1583 return build_real (type, result);
1585 else if (arg1_mode == TYPE_MODE (long_double_type_node))
1586 switch (fn)
1588 CASE_CFN_NEXTTOWARD:
1589 /* real, long double -> real. */
1590 if (fold_const_call_sss (&result, fn, TREE_REAL_CST_PTR (arg0),
1591 TREE_REAL_CST_PTR (arg1),
1592 REAL_MODE_FORMAT (mode)))
1593 return build_real (type, result);
1594 break;
1595 default:
1596 break;
1598 return NULL_TREE;
1601 if (real_cst_p (arg0)
1602 && integer_cst_p (arg1))
1604 gcc_checking_assert (SCALAR_FLOAT_MODE_P (arg0_mode));
1605 if (mode == arg0_mode)
1607 /* real, int -> real. */
1608 REAL_VALUE_TYPE result;
1609 if (fold_const_call_sss (&result, fn, TREE_REAL_CST_PTR (arg0),
1610 wi::to_wide (arg1),
1611 REAL_MODE_FORMAT (mode)))
1612 return build_real (type, result);
1614 return NULL_TREE;
1617 if (integer_cst_p (arg0)
1618 && real_cst_p (arg1))
1620 gcc_checking_assert (SCALAR_FLOAT_MODE_P (arg1_mode));
1621 if (mode == arg1_mode)
1623 /* int, real -> real. */
1624 REAL_VALUE_TYPE result;
1625 if (fold_const_call_sss (&result, fn, wi::to_wide (arg0),
1626 TREE_REAL_CST_PTR (arg1),
1627 REAL_MODE_FORMAT (mode)))
1628 return build_real (type, result);
1630 return NULL_TREE;
1633 if (arg0_mode == arg1_mode
1634 && complex_cst_p (arg0)
1635 && complex_cst_p (arg1))
1637 gcc_checking_assert (COMPLEX_MODE_P (arg0_mode));
1638 machine_mode inner_mode = GET_MODE_INNER (arg0_mode);
1639 tree arg0r = TREE_REALPART (arg0);
1640 tree arg0i = TREE_IMAGPART (arg0);
1641 tree arg1r = TREE_REALPART (arg1);
1642 tree arg1i = TREE_IMAGPART (arg1);
1643 if (mode == arg0_mode
1644 && real_cst_p (arg0r)
1645 && real_cst_p (arg0i)
1646 && real_cst_p (arg1r)
1647 && real_cst_p (arg1i))
1649 /* complex real, complex real -> complex real. */
1650 REAL_VALUE_TYPE result_real, result_imag;
1651 if (fold_const_call_ccc (&result_real, &result_imag, fn,
1652 TREE_REAL_CST_PTR (arg0r),
1653 TREE_REAL_CST_PTR (arg0i),
1654 TREE_REAL_CST_PTR (arg1r),
1655 TREE_REAL_CST_PTR (arg1i),
1656 REAL_MODE_FORMAT (inner_mode)))
1657 return build_complex (type,
1658 build_real (TREE_TYPE (type), result_real),
1659 build_real (TREE_TYPE (type), result_imag));
1661 return NULL_TREE;
1664 return NULL_TREE;
1667 /* Try to fold FN (ARG0, ARG1) to a constant. Return the constant on success,
1668 otherwise return null. TYPE is the type of the return value. */
1670 tree
1671 fold_const_call (combined_fn fn, tree type, tree arg0, tree arg1)
1673 const char *p0, *p1;
1674 char c;
1675 switch (fn)
1677 case CFN_BUILT_IN_STRSPN:
1678 if ((p0 = c_getstr (arg0)) && (p1 = c_getstr (arg1)))
1679 return build_int_cst (type, strspn (p0, p1));
1680 return NULL_TREE;
1682 case CFN_BUILT_IN_STRCSPN:
1683 if ((p0 = c_getstr (arg0)) && (p1 = c_getstr (arg1)))
1684 return build_int_cst (type, strcspn (p0, p1));
1685 return NULL_TREE;
1687 case CFN_BUILT_IN_STRCMP:
1688 if ((p0 = c_getstr (arg0)) && (p1 = c_getstr (arg1)))
1689 return build_cmp_result (type, strcmp (p0, p1));
1690 return NULL_TREE;
1692 case CFN_BUILT_IN_STRCASECMP:
1693 if ((p0 = c_getstr (arg0)) && (p1 = c_getstr (arg1)))
1695 int r = strcmp (p0, p1);
1696 if (r == 0)
1697 return build_cmp_result (type, r);
1699 return NULL_TREE;
1701 case CFN_BUILT_IN_INDEX:
1702 case CFN_BUILT_IN_STRCHR:
1703 if ((p0 = c_getstr (arg0)) && target_char_cst_p (arg1, &c))
1705 const char *r = strchr (p0, c);
1706 if (r == NULL)
1707 return build_int_cst (type, 0);
1708 return fold_convert (type,
1709 fold_build_pointer_plus_hwi (arg0, r - p0));
1711 return NULL_TREE;
1713 case CFN_BUILT_IN_RINDEX:
1714 case CFN_BUILT_IN_STRRCHR:
1715 if ((p0 = c_getstr (arg0)) && target_char_cst_p (arg1, &c))
1717 const char *r = strrchr (p0, c);
1718 if (r == NULL)
1719 return build_int_cst (type, 0);
1720 return fold_convert (type,
1721 fold_build_pointer_plus_hwi (arg0, r - p0));
1723 return NULL_TREE;
1725 case CFN_BUILT_IN_STRSTR:
1726 if ((p1 = c_getstr (arg1)))
1728 if ((p0 = c_getstr (arg0)))
1730 const char *r = strstr (p0, p1);
1731 if (r == NULL)
1732 return build_int_cst (type, 0);
1733 return fold_convert (type,
1734 fold_build_pointer_plus_hwi (arg0, r - p0));
1736 if (*p1 == '\0')
1737 return fold_convert (type, arg0);
1739 return NULL_TREE;
1741 case CFN_FOLD_LEFT_PLUS:
1742 return fold_const_fold_left (type, arg0, arg1, PLUS_EXPR);
1744 default:
1745 return fold_const_call_1 (fn, type, arg0, arg1);
1749 /* Try to evaluate:
1751 *RESULT = FN (*ARG0, *ARG1, *ARG2)
1753 in format FORMAT. Return true on success. */
1755 static bool
1756 fold_const_call_ssss (real_value *result, combined_fn fn,
1757 const real_value *arg0, const real_value *arg1,
1758 const real_value *arg2, const real_format *format)
1760 switch (fn)
1762 CASE_CFN_FMA:
1763 CASE_CFN_FMA_FN:
1764 return do_mpfr_arg3 (result, mpfr_fma, arg0, arg1, arg2, format);
1766 case CFN_FMS:
1768 real_value new_arg2 = real_value_negate (arg2);
1769 return do_mpfr_arg3 (result, mpfr_fma, arg0, arg1, &new_arg2, format);
1772 case CFN_FNMA:
1774 real_value new_arg0 = real_value_negate (arg0);
1775 return do_mpfr_arg3 (result, mpfr_fma, &new_arg0, arg1, arg2, format);
1778 case CFN_FNMS:
1780 real_value new_arg0 = real_value_negate (arg0);
1781 real_value new_arg2 = real_value_negate (arg2);
1782 return do_mpfr_arg3 (result, mpfr_fma, &new_arg0, arg1,
1783 &new_arg2, format);
1786 default:
1787 return false;
1791 /* Subroutine of fold_const_call, with the same interface. Handle cases
1792 where the arguments and result are numerical. */
1794 static tree
1795 fold_const_call_1 (combined_fn fn, tree type, tree arg0, tree arg1, tree arg2)
1797 machine_mode mode = TYPE_MODE (type);
1798 machine_mode arg0_mode = TYPE_MODE (TREE_TYPE (arg0));
1799 machine_mode arg1_mode = TYPE_MODE (TREE_TYPE (arg1));
1800 machine_mode arg2_mode = TYPE_MODE (TREE_TYPE (arg2));
1802 if (arg0_mode == arg1_mode
1803 && arg0_mode == arg2_mode
1804 && real_cst_p (arg0)
1805 && real_cst_p (arg1)
1806 && real_cst_p (arg2))
1808 gcc_checking_assert (SCALAR_FLOAT_MODE_P (arg0_mode));
1809 if (mode == arg0_mode)
1811 /* real, real, real -> real. */
1812 REAL_VALUE_TYPE result;
1813 if (fold_const_call_ssss (&result, fn, TREE_REAL_CST_PTR (arg0),
1814 TREE_REAL_CST_PTR (arg1),
1815 TREE_REAL_CST_PTR (arg2),
1816 REAL_MODE_FORMAT (mode)))
1817 return build_real (type, result);
1819 return NULL_TREE;
1822 return NULL_TREE;
1825 /* Try to fold FN (ARG0, ARG1, ARG2) to a constant. Return the constant on
1826 success, otherwise return null. TYPE is the type of the return value. */
1828 tree
1829 fold_const_call (combined_fn fn, tree type, tree arg0, tree arg1, tree arg2)
1831 const char *p0, *p1;
1832 char c;
1833 unsigned HOST_WIDE_INT s0, s1, s2 = 0;
1834 switch (fn)
1836 case CFN_BUILT_IN_STRNCMP:
1837 if (!size_t_cst_p (arg2, &s2))
1838 return NULL_TREE;
1839 if (s2 == 0
1840 && !TREE_SIDE_EFFECTS (arg0)
1841 && !TREE_SIDE_EFFECTS (arg1))
1842 return build_int_cst (type, 0);
1843 else if ((p0 = c_getstr (arg0)) && (p1 = c_getstr (arg1)))
1844 return build_int_cst (type, strncmp (p0, p1, MIN (s2, SIZE_MAX)));
1845 return NULL_TREE;
1847 case CFN_BUILT_IN_STRNCASECMP:
1848 if (!size_t_cst_p (arg2, &s2))
1849 return NULL_TREE;
1850 if (s2 == 0
1851 && !TREE_SIDE_EFFECTS (arg0)
1852 && !TREE_SIDE_EFFECTS (arg1))
1853 return build_int_cst (type, 0);
1854 else if ((p0 = c_getstr (arg0))
1855 && (p1 = c_getstr (arg1))
1856 && strncmp (p0, p1, MIN (s2, SIZE_MAX)) == 0)
1857 return build_int_cst (type, 0);
1858 return NULL_TREE;
1860 case CFN_BUILT_IN_BCMP:
1861 case CFN_BUILT_IN_MEMCMP:
1862 if (!size_t_cst_p (arg2, &s2))
1863 return NULL_TREE;
1864 if (s2 == 0
1865 && !TREE_SIDE_EFFECTS (arg0)
1866 && !TREE_SIDE_EFFECTS (arg1))
1867 return build_int_cst (type, 0);
1868 if ((p0 = getbyterep (arg0, &s0))
1869 && (p1 = getbyterep (arg1, &s1))
1870 && s2 <= s0
1871 && s2 <= s1)
1872 return build_cmp_result (type, memcmp (p0, p1, s2));
1873 return NULL_TREE;
1875 case CFN_BUILT_IN_MEMCHR:
1876 if (!size_t_cst_p (arg2, &s2))
1877 return NULL_TREE;
1878 if (s2 == 0
1879 && !TREE_SIDE_EFFECTS (arg0)
1880 && !TREE_SIDE_EFFECTS (arg1))
1881 return build_int_cst (type, 0);
1882 if ((p0 = getbyterep (arg0, &s0))
1883 && s2 <= s0
1884 && target_char_cst_p (arg1, &c))
1886 const char *r = (const char *) memchr (p0, c, s2);
1887 if (r == NULL)
1888 return build_int_cst (type, 0);
1889 return fold_convert (type,
1890 fold_build_pointer_plus_hwi (arg0, r - p0));
1892 return NULL_TREE;
1894 case CFN_WHILE_ULT:
1896 poly_uint64 parg0, parg1;
1897 if (poly_int_tree_p (arg0, &parg0) && poly_int_tree_p (arg1, &parg1))
1898 return fold_while_ult (type, parg0, parg1);
1899 return NULL_TREE;
1902 default:
1903 return fold_const_call_1 (fn, type, arg0, arg1, arg2);