Example shows numeric quadratic function use.
#include "common.hh"
#include <iostream>
using namespace roboptim;
int run_test ()
{
a.clear ();
b.clear ();
x.clear ();
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;
std::cout << f << std::endl;
x[0] = 0.1;
x[1] = 1.2;
x[2] = 2.3;
x[3] = 3.4;
x[4] = 4.5;
std::cout << "f(x) = " << f (x) << std::endl;
std::cout <<
"J(x) = " << f.
jacobian (x) << std::endl;
std::cout <<
"G(x) = " << f.
gradient (x, 0) << std::endl;
std::cout <<
"H(x) = " << f.
hessian (x, 0) << std::endl;
return 0;
}
GENERATE_TEST ()