1 // Copyright John Maddock 2006.
2 // Copyright Paul A. Bristow 2007.
3
4 // Use, modification and distribution are subject to the
5 // Boost Software License, Version 1.0. (See accompanying file
6 // LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
7
8 #include <pch.hpp>
9
10 #include <boost/math/tools/minima.hpp>
11 #define BOOST_TEST_MAIN
12 #include <boost/test/unit_test.hpp>
13 #include <boost/test/floating_point_comparison.hpp>
14 #include <boost/math/special_functions/gamma.hpp>
15 #include <iostream>
16 #include <iomanip>
17
18 template <class T>
19 struct poly_test
20 {
21 // minima is at (3,4):
operator ()poly_test22 T operator()(const T& v)
23 {
24 T a = v - 3;
25 return 3 * a * a + 4;
26 }
27 };
28
29 template <class T>
test_minima(T,const char *)30 void test_minima(T, const char* /* name */)
31 {
32 std::pair<T, T> m = boost::math::tools::brent_find_minima(poly_test<T>(), T(-10), T(10), 50);
33 BOOST_CHECK_CLOSE(m.first, T(3), T(0.001));
34 BOOST_CHECK_CLOSE(m.second, T(4), T(0.001));
35
36 T (*fp)(T);
37 #if defined(BOOST_MATH_NO_DEDUCED_FUNCTION_POINTERS)
38 fp = boost::math::lgamma<T>;
39 #else
40 fp = boost::math::lgamma;
41 #endif
42
43 m = boost::math::tools::brent_find_minima(fp, T(0.5), T(10), 50);
44 BOOST_CHECK_CLOSE(m.first, T(1.461632), T(0.1));
45 #if defined(BOOST_MATH_NO_DEDUCED_FUNCTION_POINTERS)
46 fp = boost::math::tgamma<T>;
47 #else
48 fp = boost::math::tgamma;
49 #endif
50 m = boost::math::tools::brent_find_minima(fp, T(0.5), T(10), 50);
51 BOOST_CHECK_CLOSE(m.first, T(1.461632), T(0.1));
52 }
53
BOOST_AUTO_TEST_CASE(test_main)54 BOOST_AUTO_TEST_CASE( test_main )
55 {
56 test_minima(0.1f, "float");
57 test_minima(0.1, "double");
58 #ifndef BOOST_MATH_NO_LONG_DOUBLE_MATH_FUNCTIONS
59 test_minima(0.1L, "long double");
60 #endif
61
62 }
63
64
65