debug eigen tests

Author: wgledbetter

include/boost/math/differentiation/autodiff_eigen.hpp

@@ -20,24 +20,34 @@ struct NumTraits<boost::math::differentiation::autodiff_v1::detail::
   };
 };
 
-template <typename RealType, size_t Order, typename BinaryOp, typename A>
-struct ScalarBinaryOpTraits<boost::math::differentiation::autodiff_v1::detail::
-                                template fvar<RealType, Order>,
-                            A, BinaryOp> {
-  typedef boost::math::differentiation::autodiff_v1::detail::template fvar<
-      RealType, Order>
-      ReturnType;
-};
+#define BOOST_AUTODIFF_EIGEN_SCALAR_TRAITS(A)                                 \
+  template <class Scalar, size_t Order, typename BinaryOp>                    \
+  struct ScalarBinaryOpTraits<boost::math::differentiation::autodiff_v1::     \
+                                  detail::template fvar<Scalar, Order>,       \
+                              A, BinaryOp> {                                  \
+    typedef boost::math::differentiation::autodiff_v1::detail::template fvar< \
+        Scalar, Order>                                                        \
+        ReturnType;                                                           \
+  };                                                                          \
+  template <class Scalar, size_t Order, typename BinaryOp>                    \
+  struct ScalarBinaryOpTraits<A,                                              \
+                              boost::math::differentiation::autodiff_v1::     \
+                                  detail::template fvar<Scalar, Order>,       \
+                              BinaryOp> {                                     \
+    typedef boost::math::differentiation::autodiff_v1::detail::template fvar< \
+        Scalar, Order>                                                        \
+        ReturnType;                                                           \
+  };
 
-template <typename RealType, size_t Order, typename BinaryOp, typename A>
-struct ScalarBinaryOpTraits<A,
-                            boost::math::differentiation::autodiff_v1::detail::
-                                template fvar<RealType, Order>,
-                            BinaryOp> {
-  typedef boost::math::differentiation::autodiff_v1::detail::template fvar<
-      RealType, Order>
-      ReturnType;
-};
+BOOST_AUTODIFF_EIGEN_SCALAR_TRAITS(float);
+BOOST_AUTODIFF_EIGEN_SCALAR_TRAITS(double);
+BOOST_AUTODIFF_EIGEN_SCALAR_TRAITS(long double);
+BOOST_AUTODIFF_EIGEN_SCALAR_TRAITS(short);
+BOOST_AUTODIFF_EIGEN_SCALAR_TRAITS(unsigned short);
+BOOST_AUTODIFF_EIGEN_SCALAR_TRAITS(int);
+BOOST_AUTODIFF_EIGEN_SCALAR_TRAITS(unsigned int);
+BOOST_AUTODIFF_EIGEN_SCALAR_TRAITS(long);
+BOOST_AUTODIFF_EIGEN_SCALAR_TRAITS(unsigned long);
 
 template <typename RealType, size_t Order, typename RealType2, size_t Order2,
           typename BinaryOp>
@@ -53,7 +63,7 @@ struct ScalarBinaryOpTraits<
   typedef boost::math::differentiation::autodiff_v1::detail::template fvar<
       RealReturn, ReturnOrder>
       ReturnType;
-}
+};
 }  // namespace Eigen
 
 #endif  // BOOST_MATH_DIFFERENTIATION_AUTODIFF_EIGEN_HPP

test/test_autodiff_9.cpp

@@ -224,6 +224,32 @@ BOOST_AUTO_TEST_CASE_TEMPLATE(eigen_vector, T, all_float_types) {
           xD0 * sin(xD0) - 4 * cos(xD0));
 }
 
-// BOOST_AUTO_TEST_CASE_TEMPLATE(eigen_determinant, T, all_float_types) {
-//   constexpr int dim = 4;
-// }
+BOOST_AUTO_TEST_CASE_TEMPLATE(eigen_determinant, T, all_float_types) {
+  using std::cos;
+  using std::sin;
+  constexpr int dim = 4;
+  constexpr size_t n = 1;
+
+  typedef typename fvar<T, n> fTn;
+  fTn x = 3;
+  T xD0 = x.derivative(0);
+  Eigen::Matrix<fTn, dim, dim> M = 10 * Eigen::Matrix<fTn, dim, dim>::Random();
+  M(0, 3) = x;
+  M(1, 0) = 3 * x;
+  M(1, 2) = 1 / x;
+  M(2, 2) = x * x;
+  M(3, 1) = sin(x);
+  fTn y = M.determinant();
+
+  Eigen::Matrix<T, dim, dim> dMdx = Eigen::Matrix<T, dim, dim>::Zero();
+  dMdx(0, 3) = 1;
+  dMdx(1, 0) = 3;
+  dMdx(1, 2) = -1 / (xD0 * xD0);
+  dMdx(2, 2) = 2 * xD0;
+  dMdx(3, 1) = cos(xD0);
+
+  T ans = y.derivative(0) * (M.inverse() * dMdx).trace();
+  BOOST_CHECK_EQUAL(y.derivative(1), ans);
+}
+
+BOOST_AUTO_TEST_SUITE_END()