1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
|
// Copyright (C) 2009 by Thomas Moulard, AIST, CNRS, INRIA.
//
// This file is part of the roboptim.
//
// roboptim is free software: you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// roboptim 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 Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with roboptim. If not, see <http://www.gnu.org/licenses/>.
#include "shared-tests/fixture.hh"
#include <iostream>
#include <roboptim/core/io.hh>
#include <roboptim/core/numeric-quadratic-function.hh>
using namespace roboptim;
typedef DummySolver solver_t;
BOOST_FIXTURE_TEST_SUITE (core, TestSuiteConfiguration)
BOOST_AUTO_TEST_CASE (numeric_quadratic_function)
{
boost::shared_ptr<boost::test_tools::output_test_stream>
output = retrievePattern ("numeric-quadratic-function");
NumericQuadraticFunction::matrix_t a (5, 5);
NumericQuadraticFunction::vector_t b (5);
NumericQuadraticFunction::vector_t x (5);
a.setZero ();
b.setZero ();
x.setZero ();
a(0, 0) = 1.1, a(0, 1) = 1.2, a(0, 2) = 1.3, a(0, 3) = 1.4, a(0, 4) = 1.5;
a(1, 0) = 1.2, a(1, 1) = 2.2, a(1, 2) = 2.3, a(1, 3) = 2.4, a(1, 4) = 2.5;
a(2, 0) = 1.3, a(2, 1) = 2.3, a(2, 2) = 3.3, a(2, 3) = 3.4, a(2, 4) = 3.5;
a(3, 0) = 1.4, a(3, 1) = 2.4, a(3, 2) = 3.4, a(3, 3) = 4.4, a(3, 4) = 4.5;
a(4, 0) = 1.5, a(4, 1) = 2.5, a(4, 2) = 3.5, a(4, 3) = 4.5, a(4, 4) = 5.5;
b[0] = 2.1;
b[1] = 4.3;
b[2] = 6.5;
b[3] = 8.7;
NumericQuadraticFunction f (a, b);
(*output) << f << std::endl;
x[0] = 0.1;
x[1] = 1.2;
x[2] = 2.3;
x[3] = 3.4;
x[4] = 4.5;
(*output) << "f(x) = " << f (x) << std::endl;
(*output) << "J(x) = " << f.jacobian (x) << std::endl;
(*output) << "G(x) = " << f.gradient (x, 0) << std::endl;
(*output) << "H(x) = " << f.hessian (x, 0) << std::endl;
std::cout << output->str () << std::endl;
//BOOST_CHECK (output->match_pattern ());
}
BOOST_AUTO_TEST_SUITE_END ()
|
