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, 2009 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
  28 void
  29 test01()
  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
  58 void
  59 test02()
  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
  88 void
  89 test03()
  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
 118 int
 119 main()
 120 {
 121   test01();
 122   test02();
 123   test03();
 124   return 0;
 125 }
 126