gnu/dist/gcc4/libstdc++-v3/testsuite/26_numerics/complex/13450.cc

   1 // { dg-do run { xfail broken_cplxf_arg } }
   2
   3 // Copyright (C) 2004 Free Software Foundation
   4 //
   5 // This file is part of the GNU ISO C++ Library.  This library is free
   6 // software; you can redistribute it and/or modify it under the
   7 // terms of the GNU General Public License as published by the
   8 // Free Software Foundation; either version 2, or (at your option)
   9 // any later version.
  10
  11 // This library is distributed in the hope that it will be useful,
  12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14 // GNU General Public License for more details.
  15
  16 // You should have received a copy of the GNU General Public License along
  17 // with this library; see the file COPYING.  If not, write to the Free
  18 // Software Foundation, 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301,
  19 // USA.
  20
  21 // 26.2.8 complex transcendentals
  22
  23 #include <complex>
  24 #include <limits>
  25 #include <testsuite_hooks.h>
  26
  27 template<typename T>
  28   void test01_do(T a, T b)
  29   {
  30     using namespace std;
  31     bool test __attribute__((unused)) = true;
  32     typedef complex<T> cplx;
  33
  34     T eps = numeric_limits<T>::epsilon() * 100;
  35
  36     cplx ref = pow(cplx(a, T()), cplx(b, T()));
  37     cplx res1 = pow(a, cplx(b, T()));
  38     cplx res2 = pow(cplx(a, T()), b);
  39
  40     VERIFY( abs(ref - res1) < eps );
  41     VERIFY( abs(ref - res2) < eps );
  42     VERIFY( abs(res1 - res2) < eps );
  43   }
  44
  45 // libstdc++/13450
  46 void test01()
  47 {
  48   float f1 = -1.0f;
  49   float f2 = 0.5f;
  50   test01_do(f1, f2);
  51
  52   f1 = -3.2f;
  53   f2 = 1.4f;
  54   test01_do(f1, f2);
  55
  56   double d1 = -1.0;
  57   double d2 = 0.5;
  58   test01_do(d1, d2);
  59
  60   d1 = -3.2;
  61   d2 = 1.4;
  62   test01_do(d1, d2);
  63
  64 #if __LDBL_MANT_DIG__ != 106
  65   /* For IBM long double, epsilon is too small (since 1.0 plus any
  66      double is representable) to be able to expect results within
  67      epsilon * 100 (which may be much less than 1ulp for a particular
  68      long double value).  */
  69   long double ld1 = -1.0l;
  70   long double ld2 = 0.5l;
  71   test01_do(ld1, ld2);
  72
  73   ld1 = -3.2l;
  74   ld2 = 1.4l;
  75   test01_do(ld1, ld2);
  76 #endif
  77 }
  78
  79 int main()
  80 {
  81   test01();
  82   return 0;
  83 }