libstdc++-v3/testsuite/tr1/5_numerical_facilities/special_functions/09_cyl_bessel_j/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.9 cyl_bessel_j
  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_bessel_j(nuf, xf);
  41   float b = std::tr1::cyl_bessel_jf(nuf, xf);
  42   double c = std::tr1::cyl_bessel_j(nud, xd);
  43   long double d = std::tr1::cyl_bessel_j(nul, xl);
  44   long double e = std::tr1::cyl_bessel_jl(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_bessel_j(nuf, xf);
  67   float b = std::tr1::cyl_bessel_jf(nuf, xf);
  68   double c = std::tr1::cyl_bessel_j(nud, xd);
  69   long double d = std::tr1::cyl_bessel_j(nul, xl);
  70   long double e = std::tr1::cyl_bessel_jl(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