+2015-10-27 Richard Sandiford <richard.sandiford@arm.com>
+
+ * builtins.c (fold_builtin_hypot): Delete.
+ (fold_builtin_2): Handle constant hypot arguments here.
+ * match.pd: Fold hypot(x, 0) and hypot(0, x) to x. Canonicalize
+ hypot(x, x) to fabs(x)*sqrt(2).
+
2015-10-27 Richard Sandiford <richard.sandiford@arm.com>
* gimple-match-head.c (maybe_push_res_to_seq): Use create_tmp_reg
return NULL_TREE;
}
-/* Fold a builtin function call to hypot, hypotf, or hypotl. Return
- NULL_TREE if no simplification can be made. */
-
-static tree
-fold_builtin_hypot (location_t loc, tree arg0, tree arg1, tree type)
-{
- tree res;
-
- if (!validate_arg (arg0, REAL_TYPE)
- || !validate_arg (arg1, REAL_TYPE))
- return NULL_TREE;
-
- /* Calculate the result when the argument is a constant. */
- if ((res = do_mpfr_arg2 (arg0, arg1, type, mpfr_hypot)))
- return res;
-
- /* If either argument is zero, hypot is fabs of the other. */
- if (real_zerop (arg0))
- return fold_build1_loc (loc, ABS_EXPR, type, arg1);
- else if (real_zerop (arg1))
- return fold_build1_loc (loc, ABS_EXPR, type, arg0);
-
- /* hypot(x,x) -> fabs(x)*sqrt(2). */
- if (flag_unsafe_math_optimizations
- && operand_equal_p (arg0, arg1, OEP_PURE_SAME))
- return fold_build2_loc (loc, MULT_EXPR, type,
- fold_build1_loc (loc, ABS_EXPR, type, arg0),
- build_real_truncate (type, dconst_sqrt2 ()));
-
- return NULL_TREE;
-}
-
-
/* Fold a builtin function call to pow, powf, or powl. Return
NULL_TREE if no simplification can be made. */
static tree
break;
CASE_FLT_FN (BUILT_IN_HYPOT):
- return fold_builtin_hypot (loc, arg0, arg1, type);
+ if (validate_arg (arg0, REAL_TYPE)
+ && validate_arg (arg1, REAL_TYPE))
+ return do_mpfr_arg2 (arg0, arg1, type, mpfr_hypot);
+ break;
CASE_FLT_FN (BUILT_IN_CPOW):
if (validate_arg (arg0, COMPLEX_TYPE)
(RINT integer_valued_real_p@0)
@0))
+/* hypot(x,0) and hypot(0,x) -> abs(x). */
+(simplify
+ (hypot:c @0 real_zerop@1)
+ (abs @0))
+
/* Canonicalization of sequences of math builtins. These rules represent
IL simplifications but are not necessarily optimizations.
(CABS (complex @0 @0))
(mult (abs @0) { build_real_truncate (type, dconst_sqrt2 ()); }))
+ /* hypot(x,x) -> fabs(x)*sqrt(2). */
+ (simplify
+ (HYPOT @0 @0)
+ (mult (abs @0) { build_real_truncate (type, dconst_sqrt2 ()); }))
+
/* cexp(x+yi) -> exp(x)*cexpi(y). */
(for cexps (CEXP)
exps (EXP)