5 template<
typename FunctorType,
typename Scalar>
7 const FunctorType &Functor,
8 Matrix< Scalar, Dynamic, 1 > &x,
9 Matrix< Scalar, Dynamic, 1 > &fvec,
10 Matrix< Scalar, Dynamic, Dynamic > &fjac,
11 DenseIndex ml, DenseIndex mu,
14 typedef DenseIndex Index;
25 const Scalar epsmch = NumTraits<Scalar>::epsilon();
26 const Index n = x.size();
27 assert(fvec.size()==n);
28 Matrix< Scalar, Dynamic, 1 > wa1(n);
29 Matrix< Scalar, Dynamic, 1 > wa2(n);
31 eps = sqrt((std::max)(epsfcn,epsmch));
35 for (j = 0; j < n; ++j) {
41 iflag = Functor(x, wa1);
45 fjac.col(j) = (wa1-fvec)/h;
50 for (k = 0; k < msum; ++k) {
51 for (j = k; (msum<0) ? (j>n): (j<n); j += msum) {
53 h = eps * abs(wa2[j]);
57 iflag = Functor(x, wa1);
60 for (j = k; (msum<0) ? (j>n): (j<n); j += msum) {
62 h = eps * abs(wa2[j]);
64 fjac.col(j).setZero();
65 start = std::max<Index>(0,j-mu);
66 length = (std::min)(n-1, j+ml) - start + 1;
67 fjac.col(j).segment(start, length) = ( wa1.segment(start, length)-fvec.segment(start, length))/h;