1// { dg-require-c-std "" }
2// { dg-add-options ieee }
3
4// 2007-01-10  Edward Smith-Rowland <3dw4rd@verizon.net>
5//
6// Copyright (C) 2007-2015 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 3, 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 COPYING3.  If not see
21// <http://www.gnu.org/licenses/>.
22
23// 5.2.1.13 ellint_3
24
25#include <tr1/cmath>
26#include <testsuite_hooks.h>
27
28void
29test01()
30{
31  float kf = std::numeric_limits<float>::quiet_NaN();
32  double kd = std::numeric_limits<double>::quiet_NaN();
33  long double kl = std::numeric_limits<long double>::quiet_NaN();
34
35  float nuf = 0.2F;
36  double nud = 0.2;
37  long double nul = 0.2L;
38
39  float phif = std::atan2(1.0F, 1.0F);
40  double phid = std::atan2(1.0, 1.0);
41  long double phil = std::atan2(1.0L, 1.0L);
42
43  float a = std::tr1::ellint_3(kf, nuf, phif);
44  float b = std::tr1::ellint_3f(kf, nuf, phif);
45  double c = std::tr1::ellint_3(kd, nud, phid);
46  long double d = std::tr1::ellint_3(kl, nul, phil);
47  long double e = std::tr1::ellint_3l(kl, nul, phil);
48
49  VERIFY(std::tr1::isnan<float>(a));
50  VERIFY(std::tr1::isnan<float>(b));
51  VERIFY(std::tr1::isnan<double>(c));
52  VERIFY(std::tr1::isnan<long double>(d));
53  VERIFY(std::tr1::isnan<long double>(e));
54
55  return;
56}
57
58void
59test02()
60{
61  float kf = 0.5F;
62  double kd = 0.5;
63  long double kl = 0.5L;
64
65  float nuf = std::numeric_limits<float>::quiet_NaN();
66  double nud = std::numeric_limits<double>::quiet_NaN();
67  long double nul = std::numeric_limits<long double>::quiet_NaN();
68
69  float phif = std::atan2(1.0F, 1.0F);
70  double phid = std::atan2(1.0, 1.0);
71  long double phil = std::atan2(1.0L, 1.0L);
72
73  float a = std::tr1::ellint_3(kf, nuf, phif);
74  float b = std::tr1::ellint_3f(kf, nuf, phif);
75  double c = std::tr1::ellint_3(kd, nud, phid);
76  long double d = std::tr1::ellint_3(kl, nul, phil);
77  long double e = std::tr1::ellint_3l(kl, nul, phil);
78
79  VERIFY(std::tr1::isnan<float>(a));
80  VERIFY(std::tr1::isnan<float>(b));
81  VERIFY(std::tr1::isnan<double>(c));
82  VERIFY(std::tr1::isnan<long double>(d));
83  VERIFY(std::tr1::isnan<long double>(e));
84
85  return;
86}
87
88void
89test03()
90{
91  float kf = 0.5F;
92  double kd = 0.5;
93  long double kl = 0.5L;
94
95  float nuf = 0.2F;
96  double nud = 0.2;
97  long double nul = 0.2L;
98
99  float phif = std::numeric_limits<float>::quiet_NaN();
100  double phid = std::numeric_limits<double>::quiet_NaN();
101  long double phil = std::numeric_limits<long double>::quiet_NaN();
102
103  float a = std::tr1::ellint_3(kf, nuf, phif);
104  float b = std::tr1::ellint_3f(kf, nuf, phif);
105  double c = std::tr1::ellint_3(kd, nud, phid);
106  long double d = std::tr1::ellint_3(kl, nul, phil);
107  long double e = std::tr1::ellint_3l(kl, nul, phil);
108
109  VERIFY(std::tr1::isnan<float>(a));
110  VERIFY(std::tr1::isnan<float>(b));
111  VERIFY(std::tr1::isnan<double>(c));
112  VERIFY(std::tr1::isnan<long double>(d));
113  VERIFY(std::tr1::isnan<long double>(e));
114
115  return;
116}
117
118int
119main()
120{
121  test01();
122  test02();
123  test03();
124  return 0;
125}
126
127