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

use template parameter Scalar instead of VectorType, fix a segfault.

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

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

@ -27,22 +27,19 @@
#include <cminpack.h>
template<typename Functor, typename VectorType>
template<typename Functor, typename Scalar>
// TODO : fixe Scalar here
int ei_hybrd1(
VectorType &x,
VectorType &fvec,
// ei_traits<VectorType>::Scalar tol
double tol
// = ei_sqrt(machine_epsilon<VectorType::Scalar>())
Eigen::Matrix< double, Eigen::Dynamic, 1 > &x,
Eigen::Matrix< double, Eigen::Dynamic, 1 > &fvec,
Scalar tol = Eigen::ei_sqrt(Eigen::machine_epsilon<Scalar>())
)
{
typedef typename VectorType::Scalar Scalar;
int lwa = (x.size()*(3*x.size()+13))/2;
VectorType wa(lwa);
Eigen::Matrix< double, Eigen::Dynamic, 1 > 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

4
unsupported/test/NonLinear.cpp
View File

@ -501,7 +501,7 @@ void testHybrd1()
/* unless high solutions are required, */
/* this is the recommended setting. */
info = ei_hybrd1<myfunctor,VectorXd>(x, fvec, ei_sqrt(machine_epsilon<double>()));
info = ei_hybrd1<myfunctor,double>(x, fvec);
// check return value
VERIFY( 1 == info);
@ -510,7 +510,7 @@ void testHybrd1()
VERIFY_IS_APPROX(fvec.norm(), 1.192636e-08);
// check x
VectorXd x_ref;
VectorXd x_ref(n);
x_ref << -0.5706545, -0.6816283, -0.7017325, -0.7042129, -0.701369, -0.6918656, -0.665792, -0.5960342, -0.4164121;
VERIFY_IS_APPROX(x, x_ref);
}