1 #define chkder_log10e 0.43429448190325182765 2 #define chkder_factor 100. 8 template<
typename Scalar>
23 typedef DenseIndex
Index;
25 const Scalar eps =
sqrt(NumTraits<Scalar>::epsilon());
26 const Scalar epsf = chkder_factor * NumTraits<Scalar>::epsilon();
27 const Scalar epslog = chkder_log10e *
log(eps);
30 const Index m = fvec.size(), n = x.size();
35 for (Index j = 0; j < n; ++j) {
36 temp = eps *
abs(x[j]);
45 for (Index j = 0; j < n; ++j) {
49 err += temp * fjac.col(j);
51 for (Index i = 0; i < m; ++i) {
53 if (fvec[i] != 0. && fvecp[i] != 0. &&
abs(fvecp[i] - fvec[i]) >= epsf *
abs(fvec[i]))
54 temp = eps *
abs((fvecp[i] - fvec[i]) / eps - err[i]) / (
abs(fvec[i]) +
abs(fvecp[i]));
56 if (temp > NumTraits<Scalar>::epsilon() && temp < eps)
57 err[i] = (chkder_log10e *
log(temp) - epslog) / epslog;
const Eigen::CwiseUnaryOp< Eigen::internal::scalar_sqrt_op< typename Derived::Scalar >, const Derived > sqrt(const Eigen::ArrayBase< Derived > &x)
Namespace containing all symbols from the Eigen library.
EIGEN_DEFAULT_DENSE_INDEX_TYPE Index
const Eigen::CwiseUnaryOp< Eigen::internal::scalar_abs_op< typename Derived::Scalar >, const Derived > abs(const Eigen::ArrayBase< Derived > &x)
const Eigen::CwiseUnaryOp< Eigen::internal::scalar_log_op< typename Derived::Scalar >, const Derived > log(const Eigen::ArrayBase< Derived > &x)