Updated eval implementation and tests (design-docs #19)

Author: bbbales2

stan/math/prim/fun/eval.hpp

@@ -10,26 +10,30 @@ namespace stan {
 namespace math {
 
 /**
- * Inputs that are not Eigen types are forwarded unmodified
+ * Inputs which have a plain_type equal to the own time are forwarded
+ * unmodified (for Eigen expressions these types are different)
  *
  * @tparam T Input type
  * @param[in] arg Input argument
  * @return Forwarded input argument
  **/
-template <typename T, require_not_eigen_t<T>* = nullptr>
+template <typename T,
+	  require_same_t<std::decay_t<T>, plain_type_t<T>>* = nullptr>
 inline decltype(auto) eval(T&& arg) {
   return std::forward<T>(arg);
 }
 
 /**
- * Inputs that are Eigen types are eval'd and returned
+ * Inputs which have a plain_type different from their own type are
+ * Eval'd (this catches Eigen expressions)
  *
  * @tparam T Input type
  * @param[in] arg Input argument
  * @return Eval'd argument
  **/
-template <typename T, require_eigen_t<T>* = nullptr>
-inline decltype(auto) eval(T arg) {
+template <typename T,
+	  require_not_same_t<std::decay_t<T>, plain_type_t<T>>* = nullptr>
+inline decltype(auto) eval(const T& arg) {
   return arg.eval();
 }
 

test/unit/math/prim/fun/eval_test.cpp

@@ -149,10 +149,9 @@ TEST(MathFunctions, eval_return_type_short_circuit_std_vector) {
 TEST(MathFunctions, eval_return_type_short_circuit_vector_xd) {
   Eigen::Matrix<double, Eigen::Dynamic, 1> a(5);
   const Eigen::Matrix<double, Eigen::Dynamic, 1> b(5);
-  stan::math::eval(a);
   EXPECT_TRUE(
       (std::is_same<decltype(stan::math::eval(a)),
-                    const Eigen::Matrix<double, Eigen::Dynamic, 1>&>::value));
+                    Eigen::Matrix<double, Eigen::Dynamic, 1>&>::value));
   EXPECT_TRUE(
       (std::is_same<decltype(stan::math::eval(b)),
                     const Eigen::Matrix<double, Eigen::Dynamic, 1>&>::value));
@@ -163,7 +162,7 @@ TEST(MathFunctions, eval_return_type_short_circuit_row_vector_xd) {
   const Eigen::Matrix<double, 1, Eigen::Dynamic> b(5);
   EXPECT_TRUE(
       (std::is_same<decltype(stan::math::eval(a)),
-                    const Eigen::Matrix<double, 1, Eigen::Dynamic>&>::value));
+                    Eigen::Matrix<double, 1, Eigen::Dynamic>&>::value));
   EXPECT_TRUE(
       (std::is_same<decltype(stan::math::eval(b)),
                     const Eigen::Matrix<double, 1, Eigen::Dynamic>&>::value));
@@ -173,7 +172,7 @@ TEST(MathFunctions, eval_return_type_short_circuit_matrix_xd) {
   Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic> a(5, 4);
   const Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic> b(5, 4);
   EXPECT_TRUE((std::is_same<decltype(stan::math::eval(a)),
-                            const Eigen::Matrix<double, Eigen::Dynamic,
+	                    Eigen::Matrix<double, Eigen::Dynamic,
                                                 Eigen::Dynamic>&>::value));
   EXPECT_TRUE((std::is_same<decltype(stan::math::eval(b)),
                             const Eigen::Matrix<double, Eigen::Dynamic,
@@ -201,7 +200,7 @@ TEST(MathFunctions, eval_return_type_short_circuit_static_sized_matrix) {
   Eigen::Matrix<double, 5, 4> a;
   const Eigen::Matrix<double, 5, 4> b;
   EXPECT_TRUE((std::is_same<decltype(stan::math::eval(a)),
-                            const Eigen::Matrix<double, 5, 4>&>::value));
+                            Eigen::Matrix<double, 5, 4>&>::value));
   EXPECT_TRUE((std::is_same<decltype(stan::math::eval(b)),
                             const Eigen::Matrix<double, 5, 4>&>::value));
 }