From 554fd8c5195424bdbcabf5de30fdc183aba391bd Mon Sep 17 00:00:00 2001 From: upstream source tree Date: Sun, 15 Mar 2015 20:14:05 -0400 Subject: obtained gcc-4.6.4.tar.bz2 from upstream website; verified gcc-4.6.4.tar.bz2.sig; imported gcc-4.6.4 source tree from verified upstream tarball. downloading a git-generated archive based on the 'upstream' tag should provide you with a source tree that is binary identical to the one extracted from the above tarball. if you have obtained the source via the command 'git clone', however, do note that line-endings of files in your working directory might differ from line-endings of the respective files in the upstream repository. --- .../special_functions/03_beta/check_value.cc | 465 +++++++++++++++++++++ 1 file changed, 465 insertions(+) create mode 100644 libstdc++-v3/testsuite/tr1/5_numerical_facilities/special_functions/03_beta/check_value.cc (limited to 'libstdc++-v3/testsuite/tr1/5_numerical_facilities/special_functions/03_beta/check_value.cc') diff --git a/libstdc++-v3/testsuite/tr1/5_numerical_facilities/special_functions/03_beta/check_value.cc b/libstdc++-v3/testsuite/tr1/5_numerical_facilities/special_functions/03_beta/check_value.cc new file mode 100644 index 000000000..dfadaa67c --- /dev/null +++ b/libstdc++-v3/testsuite/tr1/5_numerical_facilities/special_functions/03_beta/check_value.cc @@ -0,0 +1,465 @@ +// 2007-02-04 Edward Smith-Rowland <3dw4rd@verizon.net> +// +// Copyright (C) 2007, 2009 Free Software Foundation, Inc. +// +// This file is part of the GNU ISO C++ Library. This library is free +// software; you can redistribute it and/or modify it under the +// terms of the GNU General Public License as published by the +// Free Software Foundation; either version 3, or (at your option) +// any later version. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +// GNU General Public License for more details. +// +// You should have received a copy of the GNU General Public License along +// with this library; see the file COPYING3. If not see +// . + +// beta + + +// Compare against values generated by the GNU Scientific Library. +// The GSL can be found on the web: http://www.gnu.org/software/gsl/ + +#include +#if defined(__TEST_DEBUG) +#include +#define VERIFY(A) \ +if (!(A)) \ + { \ + std::cout << "line " << __LINE__ \ + << " max_abs_frac = " << max_abs_frac \ + << std::endl; \ + } +#else +#include +#endif +#include "../testcase.h" + + +// Test data for x=10.000000000000000. +testcase_beta data001[] = { + { 1.0825088224469029e-06, 10.000000000000000, 10.000000000000000 }, + { 4.9925087406346778e-09, 10.000000000000000, 20.000000000000000 }, + { 1.5729567312509485e-10, 10.000000000000000, 30.000000000000000 }, + { 1.2168673582561288e-11, 10.000000000000000, 40.000000000000000 }, + { 1.5916380099863291e-12, 10.000000000000000, 50.000000000000000 }, + { 2.9408957938463963e-13, 10.000000000000000, 60.000000000000000 }, + { 6.9411637980691676e-14, 10.000000000000000, 70.000000000000000 }, + { 1.9665612972502651e-14, 10.000000000000000, 80.000000000000000 }, + { 6.4187824828154399e-15, 10.000000000000000, 90.000000000000000 }, + { 2.3455339739604842e-15, 10.000000000000000, 100.00000000000000 }, +}; + +// Test function for x=10.000000000000000. +template +void test001() +{ + const Tp eps = std::numeric_limits::epsilon(); + Tp max_abs_diff = -Tp(1); + Tp max_abs_frac = -Tp(1); + unsigned int num_datum = sizeof(data001) + / sizeof(testcase_beta); + for (unsigned int i = 0; i < num_datum; ++i) + { + const Tp f = std::tr1::beta(Tp(data001[i].x), Tp(data001[i].y)); + const Tp f0 = data001[i].f0; + const Tp diff = f - f0; + if (std::abs(diff) > max_abs_diff) + max_abs_diff = std::abs(diff); + if (std::abs(f0) > Tp(10) * eps + && std::abs(f) > Tp(10) * eps) + { + const Tp frac = diff / f0; + if (std::abs(frac) > max_abs_frac) + max_abs_frac = std::abs(frac); + } + } + VERIFY(max_abs_frac < Tp(5.0000000000000029e-12)); +} + +// Test data for x=20.000000000000000. +testcase_beta data002[] = { + { 4.9925087406346778e-09, 20.000000000000000, 10.000000000000000 }, + { 7.2544445519248436e-13, 20.000000000000000, 20.000000000000000 }, + { 1.7681885473062028e-15, 20.000000000000000, 30.000000000000000 }, + { 1.7891885039182335e-17, 20.000000000000000, 40.000000000000000 }, + { 4.3240677875623635e-19, 20.000000000000000, 50.000000000000000 }, + { 1.8857342309689050e-20, 20.000000000000000, 60.000000000000000 }, + { 1.2609804003539998e-21, 20.000000000000000, 70.000000000000000 }, + { 1.1660809542079041e-22, 20.000000000000000, 80.000000000000000 }, + { 1.3907944279729071e-23, 20.000000000000000, 90.000000000000000 }, + { 2.0365059099917614e-24, 20.000000000000000, 100.00000000000000 }, +}; + +// Test function for x=20.000000000000000. +template +void test002() +{ + const Tp eps = std::numeric_limits::epsilon(); + Tp max_abs_diff = -Tp(1); + Tp max_abs_frac = -Tp(1); + unsigned int num_datum = sizeof(data002) + / sizeof(testcase_beta); + for (unsigned int i = 0; i < num_datum; ++i) + { + const Tp f = std::tr1::beta(Tp(data002[i].x), Tp(data002[i].y)); + const Tp f0 = data002[i].f0; + const Tp diff = f - f0; + if (std::abs(diff) > max_abs_diff) + max_abs_diff = std::abs(diff); + if (std::abs(f0) > Tp(10) * eps + && std::abs(f) > Tp(10) * eps) + { + const Tp frac = diff / f0; + if (std::abs(frac) > max_abs_frac) + max_abs_frac = std::abs(frac); + } + } + VERIFY(max_abs_frac < Tp(2.5000000000000020e-13)); +} + +// Test data for x=30.000000000000000. +testcase_beta data003[] = { + { 1.5729567312509485e-10, 30.000000000000000, 10.000000000000000 }, + { 1.7681885473062028e-15, 30.000000000000000, 20.000000000000000 }, + { 5.6370779640482451e-19, 30.000000000000000, 30.000000000000000 }, + { 1.0539424603796547e-21, 30.000000000000000, 40.000000000000000 }, + { 6.0118197777273836e-24, 30.000000000000000, 50.000000000000000 }, + { 7.4279528553260165e-26, 30.000000000000000, 60.000000000000000 }, + { 1.6212207780604767e-27, 30.000000000000000, 70.000000000000000 }, + { 5.4783729715317616e-29, 30.000000000000000, 80.000000000000000 }, + { 2.6183005659681346e-30, 30.000000000000000, 90.000000000000000 }, + { 1.6587948222122229e-31, 30.000000000000000, 100.00000000000000 }, +}; + +// Test function for x=30.000000000000000. +template +void test003() +{ + const Tp eps = std::numeric_limits::epsilon(); + Tp max_abs_diff = -Tp(1); + Tp max_abs_frac = -Tp(1); + unsigned int num_datum = sizeof(data003) + / sizeof(testcase_beta); + for (unsigned int i = 0; i < num_datum; ++i) + { + const Tp f = std::tr1::beta(Tp(data003[i].x), Tp(data003[i].y)); + const Tp f0 = data003[i].f0; + const Tp diff = f - f0; + if (std::abs(diff) > max_abs_diff) + max_abs_diff = std::abs(diff); + if (std::abs(f0) > Tp(10) * eps + && std::abs(f) > Tp(10) * eps) + { + const Tp frac = diff / f0; + if (std::abs(frac) > max_abs_frac) + max_abs_frac = std::abs(frac); + } + } + VERIFY(max_abs_frac < Tp(5.0000000000000039e-13)); +} + +// Test data for x=40.000000000000000. +testcase_beta data004[] = { + { 1.2168673582561288e-11, 40.000000000000000, 10.000000000000000 }, + { 1.7891885039182335e-17, 40.000000000000000, 20.000000000000000 }, + { 1.0539424603796547e-21, 40.000000000000000, 30.000000000000000 }, + { 4.6508509140090659e-25, 40.000000000000000, 40.000000000000000 }, + { 7.5161712118557719e-28, 40.000000000000000, 50.000000000000000 }, + { 3.0311331979886071e-30, 40.000000000000000, 60.000000000000000 }, + { 2.4175035070466313e-32, 40.000000000000000, 70.000000000000000 }, + { 3.2734839142758369e-34, 40.000000000000000, 80.000000000000000 }, + { 6.7690629601315579e-36, 40.000000000000000, 90.000000000000000 }, + { 1.9797337118812366e-37, 40.000000000000000, 100.00000000000000 }, +}; + +// Test function for x=40.000000000000000. +template +void test004() +{ + const Tp eps = std::numeric_limits::epsilon(); + Tp max_abs_diff = -Tp(1); + Tp max_abs_frac = -Tp(1); + unsigned int num_datum = sizeof(data004) + / sizeof(testcase_beta); + for (unsigned int i = 0; i < num_datum; ++i) + { + const Tp f = std::tr1::beta(Tp(data004[i].x), Tp(data004[i].y)); + const Tp f0 = data004[i].f0; + const Tp diff = f - f0; + if (std::abs(diff) > max_abs_diff) + max_abs_diff = std::abs(diff); + if (std::abs(f0) > Tp(10) * eps + && std::abs(f) > Tp(10) * eps) + { + const Tp frac = diff / f0; + if (std::abs(frac) > max_abs_frac) + max_abs_frac = std::abs(frac); + } + } + VERIFY(max_abs_frac < Tp(2.5000000000000015e-12)); +} + +// Test data for x=50.000000000000000. +testcase_beta data005[] = { + { 1.5916380099863291e-12, 50.000000000000000, 10.000000000000000 }, + { 4.3240677875623635e-19, 50.000000000000000, 20.000000000000000 }, + { 6.0118197777273836e-24, 50.000000000000000, 30.000000000000000 }, + { 7.5161712118557719e-28, 50.000000000000000, 40.000000000000000 }, + { 3.9646612085674138e-31, 50.000000000000000, 50.000000000000000 }, + { 5.8425643906418403e-34, 50.000000000000000, 60.000000000000000 }, + { 1.8672362180783552e-36, 50.000000000000000, 70.000000000000000 }, + { 1.0939382296458962e-38, 50.000000000000000, 80.000000000000000 }, + { 1.0442781609881063e-40, 50.000000000000000, 90.000000000000000 }, + { 1.4904121110954370e-42, 50.000000000000000, 100.00000000000000 }, +}; + +// Test function for x=50.000000000000000. +template +void test005() +{ + const Tp eps = std::numeric_limits::epsilon(); + Tp max_abs_diff = -Tp(1); + Tp max_abs_frac = -Tp(1); + unsigned int num_datum = sizeof(data005) + / sizeof(testcase_beta); + for (unsigned int i = 0; i < num_datum; ++i) + { + const Tp f = std::tr1::beta(Tp(data005[i].x), Tp(data005[i].y)); + const Tp f0 = data005[i].f0; + const Tp diff = f - f0; + if (std::abs(diff) > max_abs_diff) + max_abs_diff = std::abs(diff); + if (std::abs(f0) > Tp(10) * eps + && std::abs(f) > Tp(10) * eps) + { + const Tp frac = diff / f0; + if (std::abs(frac) > max_abs_frac) + max_abs_frac = std::abs(frac); + } + } + VERIFY(max_abs_frac < Tp(2.5000000000000020e-13)); +} + +// Test data for x=60.000000000000000. +testcase_beta data006[] = { + { 2.9408957938463963e-13, 60.000000000000000, 10.000000000000000 }, + { 1.8857342309689050e-20, 60.000000000000000, 20.000000000000000 }, + { 7.4279528553260165e-26, 60.000000000000000, 30.000000000000000 }, + { 3.0311331979886071e-30, 60.000000000000000, 40.000000000000000 }, + { 5.8425643906418403e-34, 60.000000000000000, 50.000000000000000 }, + { 3.4501231469782229e-37, 60.000000000000000, 60.000000000000000 }, + { 4.7706855386086599e-40, 60.000000000000000, 70.000000000000000 }, + { 1.2902663809722593e-42, 60.000000000000000, 80.000000000000000 }, + { 6.0105571058570508e-45, 60.000000000000000, 90.000000000000000 }, + { 4.3922898898347209e-47, 60.000000000000000, 100.00000000000000 }, +}; + +// Test function for x=60.000000000000000. +template +void test006() +{ + const Tp eps = std::numeric_limits::epsilon(); + Tp max_abs_diff = -Tp(1); + Tp max_abs_frac = -Tp(1); + unsigned int num_datum = sizeof(data006) + / sizeof(testcase_beta); + for (unsigned int i = 0; i < num_datum; ++i) + { + const Tp f = std::tr1::beta(Tp(data006[i].x), Tp(data006[i].y)); + const Tp f0 = data006[i].f0; + const Tp diff = f - f0; + if (std::abs(diff) > max_abs_diff) + max_abs_diff = std::abs(diff); + if (std::abs(f0) > Tp(10) * eps + && std::abs(f) > Tp(10) * eps) + { + const Tp frac = diff / f0; + if (std::abs(frac) > max_abs_frac) + max_abs_frac = std::abs(frac); + } + } + VERIFY(max_abs_frac < Tp(1.0000000000000008e-12)); +} + +// Test data for x=70.000000000000000. +testcase_beta data007[] = { + { 6.9411637980691676e-14, 70.000000000000000, 10.000000000000000 }, + { 1.2609804003539998e-21, 70.000000000000000, 20.000000000000000 }, + { 1.6212207780604767e-27, 70.000000000000000, 30.000000000000000 }, + { 2.4175035070466313e-32, 70.000000000000000, 40.000000000000000 }, + { 1.8672362180783552e-36, 70.000000000000000, 50.000000000000000 }, + { 4.7706855386086599e-40, 70.000000000000000, 60.000000000000000 }, + { 3.0453137143486369e-43, 70.000000000000000, 70.000000000000000 }, + { 4.0192274082013779e-46, 70.000000000000000, 80.000000000000000 }, + { 9.5865870063501807e-49, 70.000000000000000, 90.000000000000000 }, + { 3.7409127305819802e-51, 70.000000000000000, 100.00000000000000 }, +}; + +// Test function for x=70.000000000000000. +template +void test007() +{ + const Tp eps = std::numeric_limits::epsilon(); + Tp max_abs_diff = -Tp(1); + Tp max_abs_frac = -Tp(1); + unsigned int num_datum = sizeof(data007) + / sizeof(testcase_beta); + for (unsigned int i = 0; i < num_datum; ++i) + { + const Tp f = std::tr1::beta(Tp(data007[i].x), Tp(data007[i].y)); + const Tp f0 = data007[i].f0; + const Tp diff = f - f0; + if (std::abs(diff) > max_abs_diff) + max_abs_diff = std::abs(diff); + if (std::abs(f0) > Tp(10) * eps + && std::abs(f) > Tp(10) * eps) + { + const Tp frac = diff / f0; + if (std::abs(frac) > max_abs_frac) + max_abs_frac = std::abs(frac); + } + } + VERIFY(max_abs_frac < Tp(1.0000000000000008e-12)); +} + +// Test data for x=80.000000000000000. +testcase_beta data008[] = { + { 1.9665612972502651e-14, 80.000000000000000, 10.000000000000000 }, + { 1.1660809542079041e-22, 80.000000000000000, 20.000000000000000 }, + { 5.4783729715317616e-29, 80.000000000000000, 30.000000000000000 }, + { 3.2734839142758369e-34, 80.000000000000000, 40.000000000000000 }, + { 1.0939382296458962e-38, 80.000000000000000, 50.000000000000000 }, + { 1.2902663809722593e-42, 80.000000000000000, 60.000000000000000 }, + { 4.0192274082013779e-46, 80.000000000000000, 70.000000000000000 }, + { 2.7160590828669411e-49, 80.000000000000000, 80.000000000000000 }, + { 3.4593773902125368e-52, 80.000000000000000, 90.000000000000000 }, + { 7.4807039968503468e-55, 80.000000000000000, 100.00000000000000 }, +}; + +// Test function for x=80.000000000000000. +template +void test008() +{ + const Tp eps = std::numeric_limits::epsilon(); + Tp max_abs_diff = -Tp(1); + Tp max_abs_frac = -Tp(1); + unsigned int num_datum = sizeof(data008) + / sizeof(testcase_beta); + for (unsigned int i = 0; i < num_datum; ++i) + { + const Tp f = std::tr1::beta(Tp(data008[i].x), Tp(data008[i].y)); + const Tp f0 = data008[i].f0; + const Tp diff = f - f0; + if (std::abs(diff) > max_abs_diff) + max_abs_diff = std::abs(diff); + if (std::abs(f0) > Tp(10) * eps + && std::abs(f) > Tp(10) * eps) + { + const Tp frac = diff / f0; + if (std::abs(frac) > max_abs_frac) + max_abs_frac = std::abs(frac); + } + } + VERIFY(max_abs_frac < Tp(5.0000000000000029e-12)); +} + +// Test data for x=90.000000000000000. +testcase_beta data009[] = { + { 6.4187824828154399e-15, 90.000000000000000, 10.000000000000000 }, + { 1.3907944279729071e-23, 90.000000000000000, 20.000000000000000 }, + { 2.6183005659681346e-30, 90.000000000000000, 30.000000000000000 }, + { 6.7690629601315579e-36, 90.000000000000000, 40.000000000000000 }, + { 1.0442781609881063e-40, 90.000000000000000, 50.000000000000000 }, + { 6.0105571058570508e-45, 90.000000000000000, 60.000000000000000 }, + { 9.5865870063501807e-49, 90.000000000000000, 70.000000000000000 }, + { 3.4593773902125368e-52, 90.000000000000000, 80.000000000000000 }, + { 2.4416737907558032e-55, 90.000000000000000, 90.000000000000000 }, + { 3.0238531916564246e-58, 90.000000000000000, 100.00000000000000 }, +}; + +// Test function for x=90.000000000000000. +template +void test009() +{ + const Tp eps = std::numeric_limits::epsilon(); + Tp max_abs_diff = -Tp(1); + Tp max_abs_frac = -Tp(1); + unsigned int num_datum = sizeof(data009) + / sizeof(testcase_beta); + for (unsigned int i = 0; i < num_datum; ++i) + { + const Tp f = std::tr1::beta(Tp(data009[i].x), Tp(data009[i].y)); + const Tp f0 = data009[i].f0; + const Tp diff = f - f0; + if (std::abs(diff) > max_abs_diff) + max_abs_diff = std::abs(diff); + if (std::abs(f0) > Tp(10) * eps + && std::abs(f) > Tp(10) * eps) + { + const Tp frac = diff / f0; + if (std::abs(frac) > max_abs_frac) + max_abs_frac = std::abs(frac); + } + } + VERIFY(max_abs_frac < Tp(1.0000000000000008e-12)); +} + +// Test data for x=100.00000000000000. +testcase_beta data010[] = { + { 2.3455339739604842e-15, 100.00000000000000, 10.000000000000000 }, + { 2.0365059099917614e-24, 100.00000000000000, 20.000000000000000 }, + { 1.6587948222122229e-31, 100.00000000000000, 30.000000000000000 }, + { 1.9797337118812366e-37, 100.00000000000000, 40.000000000000000 }, + { 1.4904121110954370e-42, 100.00000000000000, 50.000000000000000 }, + { 4.3922898898347209e-47, 100.00000000000000, 60.000000000000000 }, + { 3.7409127305819802e-51, 100.00000000000000, 70.000000000000000 }, + { 7.4807039968503468e-55, 100.00000000000000, 80.000000000000000 }, + { 3.0238531916564246e-58, 100.00000000000000, 90.000000000000000 }, + { 2.2087606931991853e-61, 100.00000000000000, 100.00000000000000 }, +}; + +// Test function for x=100.00000000000000. +template +void test010() +{ + const Tp eps = std::numeric_limits::epsilon(); + Tp max_abs_diff = -Tp(1); + Tp max_abs_frac = -Tp(1); + unsigned int num_datum = sizeof(data010) + / sizeof(testcase_beta); + for (unsigned int i = 0; i < num_datum; ++i) + { + const Tp f = std::tr1::beta(Tp(data010[i].x), Tp(data010[i].y)); + const Tp f0 = data010[i].f0; + const Tp diff = f - f0; + if (std::abs(diff) > max_abs_diff) + max_abs_diff = std::abs(diff); + if (std::abs(f0) > Tp(10) * eps + && std::abs(f) > Tp(10) * eps) + { + const Tp frac = diff / f0; + if (std::abs(frac) > max_abs_frac) + max_abs_frac = std::abs(frac); + } + } + VERIFY(max_abs_frac < Tp(5.0000000000000029e-12)); +} + +int main(int, char**) +{ + test001(); + test002(); + test003(); + test004(); + test005(); + test006(); + test007(); + test008(); + test009(); + test010(); + return 0; +} -- cgit v1.2.3