CppAD与Eigen

时间:2018-04-19 10:38:25

标签: c++ eigen eigen3

我似乎无法找到关于如何在CppAD中使用Eigen的优秀文档。我希望看到的具体问题是如何将双精度矩阵与CppAD::AD<double>矩阵相乘。这是一个不编译的MWE:

#include "cppad/example/cppad_eigen.hpp"
#include "Eigen/Dense"

int main() {

    using ADdouble= CppAD::AD<double>;

    const int n = 3;
    Eigen::Matrix<double, n, n> A = Eigen::Matrix<double, n, n>::Zero();
    Eigen::Matrix<ADdouble, n, n> B = Eigen::Matrix<ADdouble, n, n>::Zero();
    Eigen::Matrix<ADdouble, n, n> C = A * B;
    return 0;
}

我的印象是我只需要包含

#include "cppad/example/cppad_eigen.hpp"

那会解决这些问题。

修改

进一步检查后,documentation提供了两个例子:

  1. https://www.coin-or.org/CppAD/Doc/eigen_array.cpp.htm
  2. https://www.coin-or.org/CppAD/Doc/eigen_det.cpp.htm
  3. 但两者都没有显示如何混合类型。

    修改

    以下是完整的错误消息

    C:\Users\matth\.CLion2018.1\system\cygwin_cmake\bin\cmake.exe --build /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/src/main/cpp/cmake-build-debug --target tester -- -j 4
    Scanning dependencies of target tester
    [ 50%] Building CXX object CMakeFiles/tester.dir/test.cpp.o
    In file included from /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/Core:463:0,
                     from /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/cppad/example/cppad_eigen.hpp:73,
                     from /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/src/main/cpp/test.cpp:1:
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/Product.h: In instantiation of 'struct Eigen::internal::traits<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >':
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/Product.h:115:7:   required from 'class Eigen::internal::dense_product_base<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0, 3>'
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/Product.h:147:7:   required from 'class Eigen::ProductImpl<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0, Eigen::Dense>'
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/Product.h:71:7:   required from 'class Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>'
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/src/main/cpp/test.cpp:11:43:   required from here
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/Product.h:29:127: error: no type named 'ReturnType' in 'struct Eigen::ScalarBinaryOpTraits<double, CppAD::AD<double>, Eigen::internal::scalar_product_op<double, CppAD::AD<double> > >'
       typedef typename ScalarBinaryOpTraits<typename traits<LhsCleaned>::Scalar, typename traits<RhsCleaned>::Scalar>::ReturnType Scalar;
                                                                                                                                   ^~~~~~
    In file included from /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/Core:459:0,
                     from /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/cppad/example/cppad_eigen.hpp:73,
                     from /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/src/main/cpp/test.cpp:1:
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/DenseBase.h: In instantiation of 'class Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >':
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/MatrixBase.h:48:34:   required from 'class Eigen::MatrixBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >'
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/Product.h:115:7:   required from 'class Eigen::internal::dense_product_base<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0, 3>'
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/Product.h:147:7:   required from 'class Eigen::ProductImpl<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0, Eigen::Dense>'
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/Product.h:71:7:   required from 'class Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>'
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/src/main/cpp/test.cpp:11:43:   required from here
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/DenseBase.h:83:17: error: no members matching 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}::coeff' in 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}'
         using Base::coeff;
                     ^~~~~
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/DenseBase.h:84:17: error: no members matching 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}::coeffByOuterInner' in 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}'
         using Base::coeffByOuterInner;
                     ^~~~~~~~~~~~~~~~~
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/DenseBase.h:85:26: error: no members matching 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}::operator()' in 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}'
         using Base::operator();
                              ^
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/DenseBase.h:86:26: error: no members matching 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}::operator[]' in 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}'
         using Base::operator[];
                              ^
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/DenseBase.h:87:17: error: no members matching 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}::x' in 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}'
         using Base::x;
                     ^
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/DenseBase.h:88:17: error: no members matching 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}::y' in 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}'
         using Base::y;
                     ^
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/DenseBase.h:89:17: error: no members matching 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}::z' in 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}'
         using Base::z;
                     ^
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/DenseBase.h:90:17: error: no members matching 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}::w' in 'Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseCoeffsBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>, 0>}'
         using Base::w;
                     ^
    In file included from /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/Core:460:0,
                     from /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/cppad/example/cppad_eigen.hpp:73,
                     from /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/src/main/cpp/test.cpp:1:
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/MatrixBase.h: In instantiation of 'class Eigen::MatrixBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >':
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/Product.h:115:7:   required from 'class Eigen::internal::dense_product_base<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0, 3>'
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/Product.h:147:7:   required from 'class Eigen::ProductImpl<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0, Eigen::Dense>'
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/Product.h:71:7:   required from 'class Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0>'
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/src/main/cpp/test.cpp:11:43:   required from here
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/MatrixBase.h:75:17: error: no members matching 'Eigen::MatrixBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >}::coeff' in 'Eigen::MatrixBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >}'
         using Base::coeff;
                     ^~~~~
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/MatrixBase.h:78:17: error: no members matching 'Eigen::MatrixBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >}::eval' in 'Eigen::MatrixBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >}'
         using Base::eval;
                     ^~~~
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/MatrixBase.h:79:25: error: no members matching 'Eigen::MatrixBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >}::operator-' in 'Eigen::MatrixBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >}'
         using Base::operator-;
                             ^
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/MatrixBase.h:82:25: error: no members matching 'Eigen::MatrixBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >}::operator*=' in 'Eigen::MatrixBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >}'
         using Base::operator*=;
                             ^~
    /cygdrive/c/Dropbox/android/TrajectoryOptimization/app/libs/include/Eigen/src/Core/MatrixBase.h:83:25: error: no members matching 'Eigen::MatrixBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >}::operator/=' in 'Eigen::MatrixBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >::Base {aka class Eigen::DenseBase<Eigen::Product<Eigen::Matrix<double, 3, 3>, Eigen::Matrix<CppAD::AD<double>, 3, 3>, 0> >}'
         using Base::operator/=;
                             ^~
    make[3]: *** [CMakeFiles/tester.dir/build.make:63: CMakeFiles/tester.dir/test.cpp.o] Error 1
    make[2]: *** [CMakeFiles/Makefile2:68: CMakeFiles/tester.dir/all] Error 2
    make[1]: *** [CMakeFiles/Makefile2:80: CMakeFiles/tester.dir/rule] Error 2
    make: *** [Makefile:118: tester] Error 2
    

    修改

    根据建议的答案,我创建了标题mixed_types.h

    #ifndef MIXED_TYPES_HEADER
    #define MIXED_TYPES_HEADER
    
    #include <Eigen/Core>
    #include <cppad/cppad.hpp>
    
    namespace Eigen {
    
    template<typename S, typename BinOp>
    struct ScalarBinaryOpTraits<CppAD::AD<S>, S, BinOp> {
        typedef CppAD::AD<S> ReturnType;
    };
    
    template<typename S, typename BinOp>
    struct ScalarBinaryOpTraits<S, CppAD::AD<S>, BinOp> {
        typedef CppAD::AD<S> ReturnType;
    };
    
    }
    
    #endif /* MIXED_TYPES_HEADER */
    

    我尝试在

    形式的产品中使用它
    A = B * C
    

    其中Scalar类型AB类型为CppAD::AD<double>Scalar类型C类型为double 。但是,这会生成以下格式的非常长的编译器错误:

    In file included from ../src\main\cpp/mixed_types.h:4:
    In file included from ../../../../libs/include\Eigen/Core:516:
    ../../../../libs/include\Eigen/src/Core/products/GeneralBlockPanelKernel.h:999:13: error: no matching member function for call to 'madd'
                EIGEN_GEBP_ONESTEP(0);
                ^~~~~~~~~~~~~~~~~~~~~
    ../../../../libs/include\Eigen/src/Core/products/GeneralBlockPanelKernel.h:980:22: note: expanded from macro 'EIGEN_GEBP_ONESTEP'
                  traits.madd(A0, B_0, C0, T0); \
                  ~~~~~~~^~~~
    ../../../../libs/include\Eigen/src/Core/products/GeneralMatrixMatrix.h:194:11: note: in instantiation of member function 'Eigen::internal::gebp_kernel<CppAD::AD<double>, double, long, Eigen::internal::blas_data_mapper<CppAD::AD<double>, long, 0, 0>, 2, 4, false, false>::operator()' requested here
              gebp(res.getSubMapper(i2, j2), blockA, blockB, actual_mc, actual_kc, actual_nc, alpha);
              ^
    ../../../../libs/include\Eigen/src/Core/products/GeneralMatrixMatrix.h:226:11: note: in instantiation of member function 'Eigen::internal::general_matrix_matrix_product<long, CppAD::AD<double>, 0, false, double, 0, false, 0>::run' requested here
        Gemm::run(rows, cols, m_lhs.cols(),
    

3 个答案:

答案 0 :(得分:0)

目前无法验证,但似乎您需要的代码位于 eigen_mat_mul.hpp

查看文件here

答案 1 :(得分:0)

嗯。令人惊讶的是,从原始类型doubleCppAD::AD<double>的强制转换似乎可以解决问题。例如,此代码

#include <iostream>
#include <cppad/cppad.hpp>
#include <cppad/example/cppad_eigen.hpp>
#include <Eigen/Dense>

template<typename T>
using Vector = Eigen::Matrix<T, Eigen::Dynamic, 1>;

int main() {

    /* Types and sizes */
    using ADdouble = CppAD::AD<double>;
    const int rows = 5;
    const int cols = 3;

    /* Define the independent variable that we will be differentiating with respect to */
    Vector<ADdouble> X = Vector<ADdouble>::Ones(cols);
    CppAD::Independent(X);

    /* Now define the function. The function is simply a matrix multiplication */
    Eigen::Matrix<double, rows, cols> A = Eigen::Matrix<double, rows, cols>::Random();
    Vector<ADdouble> Y = A.cast<ADdouble>() * X;
    CppAD::ADFun<double> f(X, Y);

    /* Now define the value of X that we want to evaluate the jacobian at */
    Vector<double> x = Vector<double>::Random(cols);

    /* Compute the jacobian and compare to the know value */
    Vector<double> jac = f.Jacobian(x);

    /* We are going to map the jacobian back to a matrix so that it is easier to visualize.
     * Note that CppAD seems to use row major format, while Eigen defaults to
     * column major */
    Eigen::Map<Eigen::Matrix<double, rows, cols, Eigen::RowMajor>> Jac(jac.data());
    std::cout << "CppAD Jacobian: " << std::endl << Jac << std::endl << std::endl;

    /* Now compare to the know jacobian, which is simply "A" */
    std::cout << "Known Jacobian: " << std::endl << A << std::endl << std::endl;
    return 0;
}

产生

CppAD Jacobian: 
 0.380002 -0.484536 -0.862444
0.0108368 -0.585236 -0.180186
 0.182981  0.252524  0.759988
  0.10957 -0.319746 -0.361039
-0.243142  0.687704  0.961136

Known Jacobian: 
 0.380002 -0.484536 -0.862444
0.0108368 -0.585236 -0.180186
 0.182981  0.252524  0.759988
  0.10957 -0.319746 -0.361039
-0.243142  0.687704  0.961136

但是,我不确定这是否会一直有效。我似乎对此表示怀疑,因为案例必须初始化double类型的常量矩阵的某些属性。我猜这里有一些技术上未定义的行为。如果一个Eigen / CppAD开发者可以发表评论,那就太好了。

答案 2 :(得分:0)

正确的解决方案是专门化ScalarBinaryOpTraits

namespace Eigen{
  template<typename S, typename BinOp>
  struct ScalarBinaryOpTraits<CppAD::AD<S>,S,BinOp>
  { typedef CppAD::AD<S> ReturnType; };
  template<typename S, typename BinOp>
  struct ScalarBinaryOpTraits<S,CppAD::AD<S>,BinOp>
  { typedef CppAD::AD<S> ReturnType; };
}

当然,这仅适用于operator*(double,CppAD::AD<double>)operator*(CppAD::AD<double>,double)的实施。

double矩阵投放到CppAD::AD<double>也可以,但每次都需要制作更昂贵的CppAD::AD<double>*CppAD::AD<double>产品而不是double*CppAD::AD<double>产品。