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

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