use eigen internal traits for getting the scalar type of a eigen type

Author: SteveBronder

stan/math/prim/meta/is_eigen.hpp

@@ -41,7 +41,7 @@ struct scalar_type<T, std::enable_if_t<is_eigen<T>::value>> {
  */
 template <typename T>
 struct value_type<T, std::enable_if_t<is_eigen<T>::value>> {
-  using type = typename std::decay_t<T>::Scalar;
+ using type = typename Eigen::internal::traits<std::decay_t<T>>::Scalar;
 };
 
 /*! \ingroup require_eigens_types */

test/unit/math/prim/meta/value_type_test.cpp

@@ -1,5 +1,6 @@
-#include <stan/math/prim/meta.hpp>
 #include <test/unit/util.hpp>
+#include <stan/math/prim/meta.hpp>
+#include <unsupported/Eigen/KroneckerProduct>
 #include <gtest/gtest.h>
 #include <vector>
 
@@ -35,3 +36,11 @@ TEST(MathMetaPrim, value_type_matrix) {
   EXPECT_SAME_TYPE(Eigen::RowVectorXd,
                    value_type<std::vector<Eigen::RowVectorXd> >::type);
 }
+
+TEST(MathMetaPrim, value_type_kronecker) {
+  Eigen::Matrix<double,2,2> A;
+  const auto B = Eigen::kroneckerProduct(A, Eigen::Matrix<double,2,2>::Identity());
+  Eigen::Matrix<double,4,1> C = Eigen::Matrix<double,4,1>::Random(4, 1);
+  EXPECT_TRUE((std::is_same<double, stan::value_type_t<decltype(B)>>::value));
+  Eigen::MatrixXd D = B * C;
+}