GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-14 04:35:57 +08:00
Code Issues Packages Projects Releases Wiki Activity

iwa is not really an argument, but just an old fashioned 'work array' :

Browse Source 
remove it from the eigen API
This commit is contained in:
Thomas Capricelli 2009-08-20 22:46:38 +02:00
parent b423e640a6
commit 8d2f6ad7e1
2 changed files with 2 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -250,17 +250,16 @@ template<typename Functor, typename Scalar>
int ei_lmdif1( int ei_lmdif1(
Matrix< Scalar, Dynamic, 1 > &x, Matrix< Scalar, Dynamic, 1 > &x,
Matrix< Scalar, Dynamic, 1 > &fvec, Matrix< Scalar, Dynamic, 1 > &fvec,
VectorXi &iwa,
Scalar tol = ei_sqrt(epsilon<Scalar>()) Scalar tol = ei_sqrt(epsilon<Scalar>())
) )
{ {
int n = x.size(); int n = x.size();
int ldfjac = fvec.size(); int ldfjac = fvec.size();
int lwa = ldfjac*n+5*n+ldfjac; int lwa = ldfjac*n+5*n+ldfjac;
VectorXi iwa(n);
Matrix< Scalar, Dynamic, 1 > wa(lwa); Matrix< Scalar, Dynamic, 1 > wa(lwa);
Matrix< Scalar, Dynamic, Dynamic > fjac(ldfjac, n); Matrix< Scalar, Dynamic, Dynamic > fjac(ldfjac, n);
iwa.resize(n);
wa.resize(lwa); wa.resize(lwa);
return lmdif1_template<Scalar> ( return lmdif1_template<Scalar> (
Functor::f, 0, Functor::f, 0,

3
unsupported/test/NonLinear.cpp
View File

@ -681,13 +681,12 @@ void testLmdif1()
int m=15, n=3, info; int m=15, n=3, info;
VectorXd x(n), fvec(m); VectorXd x(n), fvec(m);
VectorXi iwa;
/* the following starting values provide a rough fit. */ /* the following starting values provide a rough fit. */
x.setConstant(n, 1.); x.setConstant(n, 1.);
// do the computation // do the computation
info = ei_lmdif1<lmdif1_functor,double>(x, fvec, iwa); info = ei_lmdif1<lmdif1_functor,double>(x, fvec);
// check return value // check return value
VERIFY( 1 == info); VERIFY( 1 == info);