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