libstdc++-v3/testsuite/tr1/5_numerical_facilities/special_functions/08_cyl_bessel_i/pr56216.cc

   1 // 2013-02-08  Edward Smith-Rowland <3dw4rd@verizon.net>
   2 //
   3 // Copyright (C) 2013-2024 Free Software Foundation, Inc.
   4 //
   5 // This file is part of the GNU ISO C++ Library.  This library is free
   6 // software; you can redistribute it and/or modify it under the
   7 // terms of the GNU General Public License as published by the
   8 // Free Software Foundation; either version 3, or (at your option)
   9 // any later version.
  10 //
  11 // This library is distributed in the hope that it will be useful,
  12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14 // GNU General Public License for more details.
  15 //
  16 // You should have received a copy of the GNU General Public License along
  17 // with this library; see the file COPYING3.  If not see
  18 // <http://www.gnu.org/licenses/>.
  19
  20 // PR libstdc++/56216 - Crash of Bessel functions at x==0!
  21
  22 #include <testsuite_hooks.h>
  23 #include <tr1/cmath>
  24
  25 void
  26 test01()
  27 {
  28   double j0 = std::tr1::cyl_bessel_j(0.0, 0.0);
  29   double i0 = std::tr1::cyl_bessel_i(0.0, 0.0);
  30   double j1 = std::tr1::cyl_bessel_j(1.0, 0.0);
  31   double i1 = std::tr1::cyl_bessel_i(1.0, 0.0);
  32
  33   VERIFY(j0 == 1.0);
  34   VERIFY(i0 == 1.0);
  35   VERIFY(j1 == 0.0);
  36   VERIFY(i1 == 0.0);
  37 }
  38
  39 int
  40 main()
  41 {
  42   test01();
  43
  44   return 0;
  45 }