Merge branch 'feature/parameter-pack-odes' of https://github.com/stan-dev/math into feature/parameter-pack-odes

Author: bbbales2

stan/math/prim/fun/eval.hpp

@@ -16,8 +16,7 @@ namespace math {
  * @param[in] arg Input argument
  * @return Forwarded input argument
  **/
-template <typename T,
-	  require_not_eigen_t<T>* = nullptr>
+template <typename T, require_not_eigen_t<T>* = nullptr>
 inline decltype(auto) eval(T&& arg) {
   return std::forward<T>(arg);
 }
@@ -29,8 +28,7 @@ inline decltype(auto) eval(T&& arg) {
  * @param[in] arg Input argument
  * @return Eval'd argument
  **/
-template <typename T,
-	  require_eigen_t<T>* = nullptr>
+template <typename T, require_eigen_t<T>* = nullptr>
 inline decltype(auto) eval(T arg) {
   return arg.eval();
 }

stan/math/rev/functor/ode_store_sensitivities.hpp

@@ -55,7 +55,8 @@ Eigen::Matrix<var, Eigen::Dynamic, 1> ode_store_sensitivities(
 
   Eigen::VectorXd f_y0_t0;
   if (is_var<T_t0>::value)
-    f_y0_t0 = f(value_of(t0), eval(value_of(y0)), msgs, eval(value_of(args))...);
+    f_y0_t0
+        = f(value_of(t0), eval(value_of(y0)), msgs, eval(value_of(args))...);
 
   const size_t total_vars
       = num_y0_vars + num_args_vars + num_t0_vars + num_t_vars;

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

@@ -32,8 +32,8 @@ TEST(MathMatrixPrimArr, eval) {
 TEST(MathFunctions, eval_int_return_type_short_circuit) {
   std::vector<int> a(5, 0);
   const std::vector<int> b(5, 0);
-  EXPECT_TRUE((std::is_same<decltype(stan::math::eval(a)),
-                            std::vector<int>&>::value));
+  EXPECT_TRUE(
+      (std::is_same<decltype(stan::math::eval(a)), std::vector<int>&>::value));
   EXPECT_TRUE((std::is_same<decltype(stan::math::eval(b)),
                             const std::vector<int>&>::value));
 }
@@ -150,8 +150,9 @@ 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));
+  EXPECT_TRUE(
+      (std::is_same<decltype(stan::math::eval(a)),
+                    const 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));
@@ -160,8 +161,9 @@ TEST(MathFunctions, eval_return_type_short_circuit_vector_xd) {
 TEST(MathFunctions, eval_return_type_short_circuit_row_vector_xd) {
   Eigen::Matrix<double, 1, Eigen::Dynamic> a(5);
   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));
+  EXPECT_TRUE(
+      (std::is_same<decltype(stan::math::eval(a)),
+                    const 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));
@@ -170,9 +172,9 @@ TEST(MathFunctions, eval_return_type_short_circuit_row_vector_xd) {
 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::Dynamic>&>::value));
+  EXPECT_TRUE((std::is_same<decltype(stan::math::eval(a)),
+                            const Eigen::Matrix<double, Eigen::Dynamic,
+                                                Eigen::Dynamic>&>::value));
   EXPECT_TRUE((std::is_same<decltype(stan::math::eval(b)),
                             const Eigen::Matrix<double, Eigen::Dynamic,
                                                 Eigen::Dynamic>&>::value));
@@ -185,18 +187,21 @@ TEST(MathFunctions, eval_return_type_expression) {
   const auto& expr_a = 3 * a;
   auto expr_b = b * b;
 
-  EXPECT_TRUE((std::is_same<decltype(stan::math::eval(expr_a)),
-	       const Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic>>::value));
-  EXPECT_TRUE((std::is_same<decltype(stan::math::eval(expr_b)),
-	       const Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic>>::value));
+  EXPECT_TRUE(
+      (std::is_same<
+          decltype(stan::math::eval(expr_a)),
+          const Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic>>::value));
+  EXPECT_TRUE(
+      (std::is_same<
+          decltype(stan::math::eval(expr_b)),
+          const Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic>>::value));
 }
 
 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));
+                            const Eigen::Matrix<double, 5, 4>&>::value));
   EXPECT_TRUE((std::is_same<decltype(stan::math::eval(b)),
-	       const Eigen::Matrix<double, 5, 4>&>::value));
+                            const Eigen::Matrix<double, 5, 4>&>::value));
 }
-

test/unit/math/rev/functor/ode_store_sensitivities_test.cpp

@@ -41,10 +41,9 @@ struct aytm {
   template <typename T0, typename T_y, typename T_Arg, int R, int C>
   inline auto operator()(const T0& t, const T_y& y, std::ostream* msgs,
                          const Eigen::Matrix<T_Arg, R, C>& args) const {
-    std::vector<typename stan::return_type<T_Arg>::type> vec
-        = {sum_(args)};
-    Eigen::Matrix<stan::return_type_t<T0, T_y, T_Arg>, Eigen::Dynamic, 1>
-        out(2);
+    std::vector<typename stan::return_type<T_Arg>::type> vec = {sum_(args)};
+    Eigen::Matrix<stan::return_type_t<T0, T_y, T_Arg>, Eigen::Dynamic, 1> out(
+        2);
     out(0) = -sum_(vec) * y(0) * t;
     out(1) = -1.7 * sum_(vec) * y(1) * t;
     return out;