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

first test for a basic wrapper (and only wrapper!) for cminpack functions

Browse Source
This commit is contained in:
Thomas Capricelli 2009-08-09 03:07:34 +02:00
parent 2b9f110639
commit f19eda7cf6
2 changed files with 33 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
unsupported/Eigen/src/NonLinear/MathFunctions.h
View File

@ -27,4 +27,22 @@
#include <cminpack.h> #include <cminpack.h>
template<typename Functor, typename VectorType>
// TODO : fixe Scalar here
int ei_hybrd1(
VectorType &x,
VectorType &fvec,
// ei_traits<VectorType>::Scalar tol
double tol
// = ei::sqrt(machine_epsilon<VectorType::Scalar>())
)
{
typedef typename VectorType::Scalar Scalar;
int lwa = (x.size()*(3*x.size()+13))/2;
VectorType wa(lwa);
fvec.resize(x.size());
return hybrd1(Functor::f, 0, x.size(), x.data(), fvec.data(), tol, wa.data(), lwa);
}
#endif // EIGEN_NONLINEAR_MATHFUNCTIONS_H #endif // EIGEN_NONLINEAR_MATHFUNCTIONS_H

26
unsupported/test/NonLinear.cpp
View File

@ -483,32 +483,36 @@ int fcn_hybrd1(void * /*p*/, int n, const double *x, double *fvec, int /*iflag*/
return 0; return 0;
} }
struct myfunctor {
static int f(void *p, int n, const double *x, double *fvec, int iflag )
{ return fcn_hybrd1(p,n,x,fvec,iflag) ; }
};
void testHybrd1() void testHybrd1()
{ {
int j, n, info, lwa; int j, n=9, info;
double tol, fnorm; double fnorm;
double x[9], fvec[9], wa[180]; Eigen::VectorXd x(9), fvec(9);
n = 9;
/* the following starting values provide a rough solution. */ /* the following starting values provide a rough solution. */
for (j=1; j<=9; j++) for (j=1; j<=n; j++)
{ {
x[j-1] = -1.; x[j-1] = -1.;
} }
lwa = 180;
/* set tol to the square root of the machine precision. */ /* set tol to the square root of the machine precision. */
/* unless high solutions are required, */ /* unless high solutions are required, */
/* this is the recommended setting. */ /* this is the recommended setting. */
tol = sqrt(dpmpar(1)); info = ei_hybrd1<myfunctor,VectorXd>(x, fvec, sqrt(dpmpar(1)));
info = hybrd1(fcn_hybrd1, 0, n, x, fvec, tol, wa, lwa);
fnorm = enorm(n, fvec); fnorm = enorm(fvec.size(), fvec.data());
VERIFY_IS_APPROX(fvec.norm(), 1.192636e-08);
VERIFY_IS_APPROX(fnorm, 1.192636e-08);
VERIFY(info==1); VERIFY(info==1);
double x_ref[] = { double x_ref[] = {