GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-12 19:59:05 +08:00
Code Issues Packages Projects Releases Wiki Activity

BoxBOD : oops.. shame on me, i did a mistake in the derivative.... now we need 16

Browse Source 
iterations instead of 7693 ;-)
the first test still fails though.
This commit is contained in:
Thomas Capricelli 2009-08-12 02:34:22 +02:00
parent fd307b8f3f
commit 7b922eb634
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
unsupported/test/NonLinear.cpp
View File

@ -1273,7 +1273,7 @@ struct BoxBOD_functor {
for(i=0; i<6; i++) { for(i=0; i<6; i++) {
double e = exp(-b[1]*x[i]); double e = exp(-b[1]*x[i]);
fjac[i+ldfjac*0] = 1.-e; fjac[i+ldfjac*0] = 1.-e;
fjac[i+ldfjac*1] = b[1]*x[i]*e; fjac[i+ldfjac*1] = b[0]*x[i]*e;
} }
} }
return 0; return 0;
@ -1296,7 +1296,7 @@ void testNistBoxBOD(void)
*/ */
x<< 1., 1.; x<< 1., 1.;
// do the computation // do the computation
info = ei_lmder<BoxBOD_functor, double>(x, fvec, nfev, njev, fjac, ipvt, wa1, diag, info = ei_lmder<BoxBOD_functor, double>(x, fvec, nfev, njev, fjac, ipvt, wa1, diag);
1, 100., 14000, Eigen::machine_epsilon<double>(), Eigen::machine_epsilon<double>()); 1, 100., 14000, Eigen::machine_epsilon<double>(), Eigen::machine_epsilon<double>());
// check return value // check return value
@ -1304,8 +1304,8 @@ void testNistBoxBOD(void)
printf("norm2 = %.50g\n", fvec.squaredNorm()); printf("norm2 = %.50g\n", fvec.squaredNorm());
std::cout << x << std::endl; std::cout << x << std::endl;
VERIFY( 1 == info); VERIFY( 1 == info);
VERIFY( 55 == nfev); VERIFY( 10 == nfev);
VERIFY( 11 == njev); VERIFY( 6 == njev);
// check norm^2 // check norm^2
VERIFY_IS_APPROX(fvec.squaredNorm(), 1.1680088766E+03); VERIFY_IS_APPROX(fvec.squaredNorm(), 1.1680088766E+03);
// check x // check x
@ -1322,9 +1322,9 @@ void testNistBoxBOD(void)
1, 100., 14000, Eigen::machine_epsilon<double>(), Eigen::machine_epsilon<double>()); 1, 100., 14000, Eigen::machine_epsilon<double>(), Eigen::machine_epsilon<double>());
// check return value // check return value
VERIFY( 2 == info); VERIFY( 1 == info);
VERIFY( 7693 == nfev); VERIFY( 16 == nfev);
VERIFY( 5871 == njev); VERIFY( 15 == njev);
// check norm^2 // check norm^2
VERIFY_IS_APPROX(fvec.squaredNorm(), 1.1680088766E+03); VERIFY_IS_APPROX(fvec.squaredNorm(), 1.1680088766E+03);
// check x // check x