libstdc++-v3/testsuite/tr1/5_numerical_facilities/special_functions/11_cyl_neumann/check_nan.cc

   1 // { dg-require-c-std "" }
   2 // { dg-add-options ieee }
   3 // { dg-prune-output ".*warning: tgammal has lower than advertised precision" }
   4
   5 // 2007-01-10  Edward Smith-Rowland <3dw4rd@verizon.net>
   6 //
   7 // Copyright (C) 2007-2023 Free Software Foundation, Inc.
   8 //
   9 // This file is part of the GNU ISO C++ Library.  This library is free
  10 // software; you can redistribute it and/or modify it under the
  11 // terms of the GNU General Public License as published by the
  12 // Free Software Foundation; either version 3, or (at your option)
  13 // any later version.
  14 //
  15 // This library is distributed in the hope that it will be useful,
  16 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  18 // GNU General Public License for more details.
  19 //
  20 // You should have received a copy of the GNU General Public License along
  21 // with this library; see the file COPYING3.  If not see
  22 // <http://www.gnu.org/licenses/>.
  23
  24 // 5.2.1.11 cyl_neumann
  25
  26 #include <tr1/cmath>
  27 #include <testsuite_hooks.h>
  28
  29 void
  30 test01()
  31 {
  32   float xf = std::numeric_limits<float>::quiet_NaN();
  33   double xd = std::numeric_limits<double>::quiet_NaN();
  34   long double xl = std::numeric_limits<long double>::quiet_NaN();
  35
  36   float nuf = 0.0F;
  37   double nud = 0.0;
  38   long double nul = 0.0L;
  39
  40   float a = std::tr1::cyl_neumann(nuf, xf);
  41   float b = std::tr1::cyl_neumannf(nuf, xf);
  42   double c = std::tr1::cyl_neumann(nud, xd);
  43   long double d = std::tr1::cyl_neumann(nul, xl);
  44   long double e = std::tr1::cyl_neumannl(nul, xl);
  45
  46   VERIFY(std::tr1::isnan<float>(a));
  47   VERIFY(std::tr1::isnan<float>(b));
  48   VERIFY(std::tr1::isnan<double>(c));
  49   VERIFY(std::tr1::isnan<long double>(d));
  50   VERIFY(std::tr1::isnan<long double>(e));
  51
  52   return;
  53 }
  54
  55 void
  56 test02()
  57 {
  58   float xf = 1.0F;
  59   double xd = 1.0;
  60   long double xl = 1.0L;
  61
  62   float nuf = std::numeric_limits<float>::quiet_NaN();
  63   double nud = std::numeric_limits<double>::quiet_NaN();
  64   long double nul = std::numeric_limits<long double>::quiet_NaN();
  65
  66   float a = std::tr1::cyl_neumann(nuf, xf);
  67   float b = std::tr1::cyl_neumannf(nuf, xf);
  68   double c = std::tr1::cyl_neumann(nud, xd);
  69   long double d = std::tr1::cyl_neumann(nul, xl);
  70   long double e = std::tr1::cyl_neumannl(nul, xl);
  71
  72   VERIFY(std::tr1::isnan<float>(a));
  73   VERIFY(std::tr1::isnan<float>(b));
  74   VERIFY(std::tr1::isnan<double>(c));
  75   VERIFY(std::tr1::isnan<long double>(d));
  76   VERIFY(std::tr1::isnan<long double>(e));
  77
  78   return;
  79 }
  80
  81 int
  82 main()
  83 {
  84   test01();
  85   test02();
  86   return 0;
  87 }
  88