libstdc++-v3/testsuite/26_numerics/random/poisson_distribution/operators/values.cc

   1 // { dg-do run { target c++11 } }
   2 // { dg-require-cstdint "" }
   3 // { dg-require-cmath "" }
   4 //
   5 // Copyright (C) 2011-2018 Free Software Foundation, Inc.
   6 //
   7 // This file is part of the GNU ISO C++ Library.  This library is free
   8 // software; you can redistribute it and/or modify it under the
   9 // terms of the GNU General Public License as published by the
  10 // Free Software Foundation; either version 3, or (at your option)
  11 // any later version.
  12 //
  13 // This library is distributed in the hope that it will be useful,
  14 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  15 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  16 // GNU General Public License for more details.
  17 //
  18 // You should have received a copy of the GNU General Public License along
  19 // with this library; see the file COPYING3.  If not see
  20 // <http://www.gnu.org/licenses/>.
  21
  22 // 26.5.8.4.1 Class template poisson_distribution [rand.dist.pois.poisson]
  23
  24 #include <random>
  25 #include <functional>
  26 #include <testsuite_random.h>
  27
  28 void test01()
  29 {
  30   using namespace __gnu_test;
  31
  32   std::mt19937 eng;
  33
  34   std::poisson_distribution<> pd1(3.0);
  35   auto bpd1 = std::bind(pd1, eng);
  36   testDiscreteDist(bpd1, [](int n) { return poisson_pdf(n, 3.0); } );
  37
  38   std::poisson_distribution<> pd2(15.0);
  39   auto bpd2 = std::bind(pd2, eng);
  40   testDiscreteDist(bpd2, [](int n) { return poisson_pdf(n, 15.0); } );
  41
  42   std::poisson_distribution<> pd3(30.0);
  43   auto bpd3 = std::bind(pd3, eng);
  44   testDiscreteDist(bpd3, [](int n) { return poisson_pdf(n, 30.0); } );
  45
  46   // libstdc++/83237
  47   std::poisson_distribution<> pd4(37.17);
  48   auto bpd4 = std::bind(pd4, eng);
  49   testDiscreteDist<100, 2000000>(bpd4, [](int n)
  50                                  { return poisson_pdf(n, 37.17); } );
  51 }
  52
  53 int main()
  54 {
  55   test01();
  56   return 0;
  57 }