OSDN Git Service

50f4bad3b93e1821c2f17a7e7146a6770115dc5d
[pf3gnuchains/gcc-fork.git] / 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 }