OSDN Git Service

* testsuite/tr1/5_numerical_facilities/special_functions/
[pf3gnuchains/gcc-fork.git] / libstdc++-v3 / testsuite / tr1 / 5_numerical_facilities / special_functions / 14_ellint_3 / check_nan.cc
1 // { dg-require-c-std "" }
2 // { dg-options "-mieee" { target sh*-*-* } }
3
4 // 2007-01-10  Edward Smith-Rowland <3dw4rd@verizon.net>
5 //
6 // Copyright (C) 2007 Free Software Foundation, Inc.
7 //
8 // This file is part of the GNU ISO C++ Library.  This library is free
9 // software; you can redistribute it and/or modify it under the
10 // terms of the GNU General Public License as published by the
11 // Free Software Foundation; either version 2, or (at your option)
12 // any later version.
13 //
14 // This library is distributed in the hope that it will be useful,
15 // but WITHOUT ANY WARRANTY; without even the implied warranty of
16 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
17 // GNU General Public License for more details.
18 //
19 // You should have received a copy of the GNU General Public License along
20 // with this library; see the file COPYING.  If not, write to the Free
21 // Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307,
22 // USA.
23
24 // 5.2.1.13 ellint_3
25
26 #include <tr1/cmath>
27 #include <testsuite_hooks.h>
28
29 void
30 test01()
31 {
32   float kf = std::numeric_limits<float>::quiet_NaN();
33   double kd = std::numeric_limits<double>::quiet_NaN();
34   long double kl = std::numeric_limits<long double>::quiet_NaN();
35
36   float nuf = 0.2F;
37   double nud = 0.2;
38   long double nul = 0.2L;
39
40   float phif = std::atan2(1.0F, 1.0F);
41   double phid = std::atan2(1.0, 1.0);
42   long double phil = std::atan2(1.0L, 1.0L);
43
44   float a = std::tr1::ellint_3(kf, nuf, phif);
45   float b = std::tr1::ellint_3f(kf, nuf, phif);
46   double c = std::tr1::ellint_3(kd, nud, phid);
47   long double d = std::tr1::ellint_3(kl, nul, phil);
48   long double e = std::tr1::ellint_3l(kl, nul, phil);
49
50   VERIFY(std::tr1::isnan<float>(a));
51   VERIFY(std::tr1::isnan<float>(b));
52   VERIFY(std::tr1::isnan<double>(c));
53   VERIFY(std::tr1::isnan<long double>(d));
54   VERIFY(std::tr1::isnan<long double>(e));
55
56   return;
57 }
58
59 void
60 test02()
61 {
62   float kf = 0.5F;
63   double kd = 0.5;
64   long double kl = 0.5L;
65
66   float nuf = std::numeric_limits<float>::quiet_NaN();
67   double nud = std::numeric_limits<double>::quiet_NaN();
68   long double nul = std::numeric_limits<long double>::quiet_NaN();
69
70   float phif = std::atan2(1.0F, 1.0F);
71   double phid = std::atan2(1.0, 1.0);
72   long double phil = std::atan2(1.0L, 1.0L);
73
74   float a = std::tr1::ellint_3(kf, nuf, phif);
75   float b = std::tr1::ellint_3f(kf, nuf, phif);
76   double c = std::tr1::ellint_3(kd, nud, phid);
77   long double d = std::tr1::ellint_3(kl, nul, phil);
78   long double e = std::tr1::ellint_3l(kl, nul, phil);
79
80   VERIFY(std::tr1::isnan<float>(a));
81   VERIFY(std::tr1::isnan<float>(b));
82   VERIFY(std::tr1::isnan<double>(c));
83   VERIFY(std::tr1::isnan<long double>(d));
84   VERIFY(std::tr1::isnan<long double>(e));
85
86   return;
87 }
88
89 void
90 test03()
91 {
92   float kf = 0.5F;
93   double kd = 0.5;
94   long double kl = 0.5L;
95
96   float nuf = 0.2F;
97   double nud = 0.2;
98   long double nul = 0.2L;
99
100   float phif = std::numeric_limits<float>::quiet_NaN();
101   double phid = std::numeric_limits<double>::quiet_NaN();
102   long double phil = std::numeric_limits<long double>::quiet_NaN();
103
104   float a = std::tr1::ellint_3(kf, nuf, phif);
105   float b = std::tr1::ellint_3f(kf, nuf, phif);
106   double c = std::tr1::ellint_3(kd, nud, phid);
107   long double d = std::tr1::ellint_3(kl, nul, phil);
108   long double e = std::tr1::ellint_3l(kl, nul, phil);
109
110   VERIFY(std::tr1::isnan<float>(a));
111   VERIFY(std::tr1::isnan<float>(b));
112   VERIFY(std::tr1::isnan<double>(c));
113   VERIFY(std::tr1::isnan<long double>(d));
114   VERIFY(std::tr1::isnan<long double>(e));
115
116   return;
117 }
118
119 int
120 main()
121 {
122   test01();
123   test02();
124   test03();
125   return 0;
126 }
127