Merge commit '5c6a4dd7a814a60bd3ea766ab8701ea2234724b3' into HEAD

Author: yashikno

.github/FUNDING.yml

@@ -0,0 +1,2 @@
+github: stan-dev
+custom: https://mc-stan.org/support/

stan/math/fwd/functor.hpp

@@ -1,6 +1,7 @@
 #ifndef STAN_MATH_FWD_FUNCTOR_HPP
 #define STAN_MATH_FWD_FUNCTOR_HPP
 
+#include <stan/math/fwd/functor/apply_scalar_unary.hpp>
 #include <stan/math/fwd/functor/gradient.hpp>
 #include <stan/math/fwd/functor/hessian.hpp>
 #include <stan/math/fwd/functor/jacobian.hpp>

stan/math/fwd/meta/apply_scalar_unary.hpp …/math/fwd/functor/apply_scalar_unary.hppstan/math/fwd/meta/apply_scalar_unary.hpp renamed to stan/math/fwd/functor/apply_scalar_unary.hpp stan/math/fwd/meta/apply_scalar_unary.hpp renamed to stan/math/fwd/functor/apply_scalar_unary.hpp

@@ -1,7 +1,7 @@
-#ifndef STAN_MATH_FWD_META_APPLY_SCALAR_UNARY_HPP
-#define STAN_MATH_FWD_META_APPLY_SCALAR_UNARY_HPP
+#ifndef STAN_MATH_FWD_FUNCTOR_APPLY_SCALAR_UNARY_HPP
+#define STAN_MATH_FWD_FUNCTOR_APPLY_SCALAR_UNARY_HPP
 
-#include <stan/math/prim/meta/apply_scalar_unary.hpp>
+#include <stan/math/prim/functor/apply_scalar_unary.hpp>
 #include <stan/math/fwd/core/fvar.hpp>
 
 namespace stan {

stan/math/fwd/meta.hpp

@@ -1,7 +1,6 @@
 #ifndef STAN_MATH_FWD_META_HPP
 #define STAN_MATH_FWD_META_HPP
 
-#include <stan/math/fwd/meta/apply_scalar_unary.hpp>
 #include <stan/math/fwd/meta/is_fvar.hpp>
 #include <stan/math/fwd/meta/partials_type.hpp>
 #include <stan/math/fwd/meta/operands_and_partials.hpp>

stan/math/prim/eigen_plugins.h

@@ -53,14 +53,14 @@ struct val_Op{
   //Returns value from a vari*
   template<typename T = Scalar>
   EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
-    std::enable_if_t<std::is_pointer<T>::value, reverse_return_t<T>>
+    std::enable_if_t<std::is_pointer<T>::value, const double&>
       operator()(T &v) const { return v->val_; }
 
   //Returns value from a var
   template<typename T = Scalar>
   EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE
     std::enable_if_t<(!std::is_pointer<T>::value && !is_fvar<T>::value
-                      && !std::is_arithmetic<T>::value), reverse_return_t<T>>
+                      && !std::is_arithmetic<T>::value), const double&>
       operator()(T &v) const { return v.vi_->val_; }
 
   //Returns value from an fvar

stan/math/prim/fun/Phi.hpp

@@ -7,6 +7,7 @@
 #include <stan/math/prim/fun/erf.hpp>
 #include <stan/math/prim/fun/erfc.hpp>
 #include <stan/math/prim/fun/Phi.hpp>
+#include <stan/math/prim/functor/apply_scalar_unary.hpp>
 
 namespace stan {
 namespace math {

stan/math/prim/fun/Phi_approx.hpp

@@ -3,6 +3,7 @@
 
 #include <stan/math/prim/meta.hpp>
 #include <stan/math/prim/fun/inv_logit.hpp>
+#include <stan/math/prim/functor/apply_scalar_unary.hpp>
 #include <cmath>
 
 namespace stan {

stan/math/prim/fun/acos.hpp

@@ -10,6 +10,8 @@
 #include <stan/math/prim/fun/isfinite.hpp>
 #include <stan/math/prim/fun/isnan.hpp>
 #include <stan/math/prim/fun/polar.hpp>
+#include <stan/math/prim/functor/apply_scalar_unary.hpp>
+#include <stan/math/prim/functor/apply_vector_unary.hpp>
 #include <cmath>
 #include <complex>
 

stan/math/prim/fun/acosh.hpp

@@ -8,6 +8,7 @@
 #include <stan/math/prim/fun/is_nan.hpp>
 #include <stan/math/prim/fun/log.hpp>
 #include <stan/math/prim/fun/sqrt.hpp>
+#include <stan/math/prim/functor/apply_scalar_unary.hpp>
 #include <cmath>
 #include <complex>
 

stan/math/prim/fun/asin.hpp

@@ -8,6 +8,8 @@
 #include <stan/math/prim/fun/Eigen.hpp>
 #include <stan/math/prim/fun/i_times.hpp>
 #include <stan/math/prim/fun/value_of_rec.hpp>
+#include <stan/math/prim/functor/apply_scalar_unary.hpp>
+#include <stan/math/prim/functor/apply_vector_unary.hpp>
 #include <cmath>
 #include <complex>