gcc/testsuite/gcc.dg/torture/builtin-power-1.c

   1 /* Copyright (C) 2004  Free Software Foundation.
   2
   3    Verify that built-in folding of various math "power" functions is
   4    correctly performed by the compiler.
   5
   6    Written by Kaveh Ghazi, 2004-03-11.  */
   7
   8 /* { dg-do link } */
   9 /* { dg-options "-ffast-math" } */
  10 /* { dg-add-options c99_runtime } */
  11 /* { dg-skip-if "PR44214" { *-*-* } { "-O0" } { "" } } */
  12
  13 #include "../builtins-config.h"
  14
  15 #ifdef HAVE_C99_RUNTIME
  16 #define C99CODE(CODE) CODE
  17 #else
  18 #define C99CODE(CODE) 0
  19 #endif
  20
  21 #define PROTOTYPE(FN) extern double FN(double); extern float FN##f(float); \
  22   extern long double FN##l(long double);
  23 #define PROTOTYPE2(FN) extern double FN(double, double); \
  24   extern float FN##f(float, float); \
  25   extern long double FN##l(long double, long double);
  26
  27 PROTOTYPE(fabs)
  28 PROTOTYPE(sqrt)
  29 PROTOTYPE(cbrt)
  30 PROTOTYPE2(pow)
  31
  32 void test(double d1, double d2, double d3,
  33           float f1, float f2, float f3,
  34           long double ld1, long double ld2, long double ld3)
  35 {
  36   /* Test N1root(N2root(x)) -> pow(x,1/(N1*N2)).  */
  37   /* E.g. sqrt(cbrt(x)) -> pow(x,1/6).  */
  38   /* The `ABS' argument is `fabs' when the transformation only works
  39      for nonnegative arguments.  Otherwise it's blank.  */
  40 #define ROOT_ROOT(FN1,N1,FN2,N2,ABS) \
  41  extern void link_failure_##FN1##_##FN2(void); \
  42  if (FN1(FN2(ABS(d1))) != pow(ABS(d1),1.0/(N1*N2)) \
  43      || C99CODE (FN1##f(FN2##f(ABS(f1))) != powf(ABS(f1),1.0F/(N1*N2))) \
  44      || C99CODE (FN1##l(FN2##l(ABS(ld1))) != powl(ABS(ld1),1.0L/(N1*N2)))) \
  45     link_failure_##FN1##_##FN2()
  46
  47   ROOT_ROOT(sqrt,2,sqrt,2,);
  48   ROOT_ROOT(sqrt,2,cbrt,3,);
  49   ROOT_ROOT(cbrt,3,sqrt,2,);
  50   ROOT_ROOT(cbrt,3,cbrt,3,fabs);
  51
  52   /* Test pow(Nroot(x),y) -> pow(x,y/N).  */
  53   /* The `ABS' argument is `fabs' when the transformation only works
  54      for nonnegative arguments.  Otherwise it's blank.  */
  55 #define POW_ROOT(FN,N,ABS) \
  56  extern void link_failure_pow_##FN(void); \
  57  if (pow(FN(ABS(d1)), d2) != pow(ABS(d1),d2/N) \
  58      || powf(FN##f(ABS(f1)),f2) != powf(ABS(f1),f2/N) \
  59      || powl(FN##l(ABS(ld1)),ld2) != powl(ABS(ld1),ld2/N)) \
  60     link_failure_pow_##FN()
  61
  62   POW_ROOT(sqrt,2,);
  63   POW_ROOT(cbrt,3,fabs);
  64
  65   /* Test Nroot(pow(x,y)) -> pow(x,y/N).  */
  66   /* The `ABS' argument is `fabs' when the transformation only works
  67      for nonnegative arguments.  Otherwise it's blank.  */
  68 #define ROOT_POW(FN,N,ABS) \
  69  extern void link_failure_##FN##_pow(void); \
  70  if (FN(pow(ABS(d1), d2)) != pow(ABS(d1),d2/N) \
  71      || FN##f(powf(ABS(f1),f2)) != powf(ABS(f1),f2/N) \
  72      || FN##l(powl(ABS(ld1),ld2)) != powl(ABS(ld1),ld2/N)) \
  73     link_failure_##FN##_pow()
  74
  75   ROOT_POW(sqrt,2,fabs);
  76   ROOT_POW(cbrt,3,fabs);
  77
  78   /* Test pow(pow(x,y),z) -> pow(x,y*z).  */
  79 #define POW_POW \
  80  extern void link_failure_pow_pow(void); \
  81  if (pow(pow(fabs(d1), d2), d3) != pow(fabs(d1),d2*d3) \
  82      || powf(powf(fabs(f1),f2),f3) != powf(fabs(f1),f2*f3) \
  83      || powl(powl(fabs(ld1),ld2),ld3) != powl(fabs(ld1),ld2*ld3)) \
  84     link_failure_pow_pow()
  85
  86   POW_POW;
  87
  88   /* Test Nroot(x)*Nroot(y) -> Nroot(x*y).  */
  89 #define ROOT_X_ROOT(FN) \
  90  extern void link_failure_root_x_root(void); \
  91  if (FN(d1)*FN(d2) != FN(d1*d2) \
  92      || FN##f(f1)*FN##f(f2) != FN##f(f1*f2) \
  93      || FN##l(ld1)*FN##l(ld2) != FN##l(ld1*ld2)) \
  94     link_failure_root_x_root()
  95
  96   ROOT_X_ROOT(sqrt);
  97   ROOT_X_ROOT(cbrt);
  98
  99   /* Test pow(x,y)*pow(x,z) -> pow(x,y+z).  */
 100 #define POW_X_POW \
 101  extern void link_failure_pow_x_pow(void); \
 102  if (pow(d1,d2)*pow(d1,d3) != pow(d1,d2+d3) \
 103      || powf(f1,f2)*powf(f1,f3) != powf(f1,f2+f3) \
 104      || powl(ld1,ld2)*powl(ld1,ld3) != powl(ld1,ld2+ld3)) \
 105     link_failure_pow_x_pow()
 106
 107   POW_X_POW;
 108
 109 }
 110
 111 int main (void)
 112 {
 113   return 0;
 114 }