42 if (
os_) *
os_ <<
"testDouble()\n";
45 ublas::matrix<double>
A(2, 2);
46 A(0,0) = 1;
A(0,1) = 2;
47 A(1,0) = 3;
A(1,1) = 4;
49 ublas::vector<double>
y(2);
53 ublas::vector<double>
x = lls.solve(A, y);
55 if (
os_) *
os_ <<
"A: " << A << endl;
56 if (
os_) *
os_ <<
"y: " << y << endl;
57 if (
os_) *
os_ <<
"x: " << x << endl;
65 if (
os_) *
os_ <<
"testDoubleQR()\n";
68 ublas::matrix<double>
A(2, 2);
69 A(0,0) = 1;
A(0,1) = 2;
70 A(1,0) = 3;
A(1,1) = 4;
72 ublas::vector<double>
y(2);
76 ublas::vector<double>
x = lls.
solve(A, y);
78 if (
os_) *
os_ <<
"A: " << A << endl;
79 if (
os_) *
os_ <<
"y: " << y << endl;
80 if (
os_) *
os_ <<
"x: " << x << endl;
82 if (
os_) *
os_ <<
"x(0) = " <<
x(0) - 1. <<
", x(1) = " <<
x(1) - 2. << endl;
89 if (
os_) *
os_ <<
"***************************\n";
90 if (
os_) *
os_ <<
"testExactFit()\n";
109 if (
os_) *
os_ <<
"testExactFit(): success\n";
115 if (
os_) *
os_ <<
"***************************\n";
116 if (
os_) *
os_ <<
"testSimpleRectangleQR()\n";
143 if (
os_) *
os_ <<
"testSimpleRectangleQR(): success\n";
148 if (
os_) *
os_ <<
"***************************\n";
149 if (
os_) *
os_ <<
"testLeastSquaresQR()\n";
179 if (
os_) *
os_ <<
"testLeastSquaresQR(): success\n";
182 int main(
int argc,
char* argv[])
188 if (argc>1 && !strcmp(argv[1],
"-v"))
os_ = &cout;
189 if (
os_) *
os_ <<
"LinearLeastSquaresTest\n";
#define unit_assert_equal(x, y, epsilon)
boost::numeric::ublas::vector< double > dvector
void testSimpleRectangleQR()
void testLeastSquaresQR()
int main(int argc, char *argv[])
boost::numeric::ublas::vector< T > solve(const boost::numeric::ublas::matrix< T > &A, const boost::numeric::ublas::vector< T > &x)
#define TEST_PROLOG(argc, argv)
KernelTraitsBase< Kernel >::space_type::abscissa_type x
boost::numeric::ublas::matrix< double > dmatrix
KernelTraitsBase< Kernel >::space_type::ordinate_type y