1#ifndef STAN_MATH_PRIM_PROB_NORMAL_LPDF_HPP
2#define STAN_MATH_PRIM_PROB_NORMAL_LPDF_HPP
41template <
bool propto,
typename T_y,
typename T_loc,
typename T_scale,
43 T_y, T_loc, T_scale>* =
nullptr>
46 const T_scale& sigma) {
51 static constexpr const char* function =
"normal_lpdf";
53 mu,
"Scale parameter", sigma);
56 T_sigma_ref sigma_ref = sigma;
76 = to_ref_if<!is_constant_all<T_y, T_scale, T_loc>::value>(
inv(sigma_val));
77 const auto& y_scaled =
to_ref((y_val - mu_val) * inv_sigma);
78 const auto& y_scaled_sq
79 = to_ref_if<!is_constant_all<T_scale>::value>(y_scaled * y_scaled);
82 T_partials_return logp = -0.5 *
sum(y_scaled_sq);
91 auto scaled_diff = to_ref_if<!is_constant_all<T_y>::value
94 >= 2>(inv_sigma * y_scaled);
96 partials<0>(ops_partials) = -scaled_diff;
99 partials<2>(ops_partials) = inv_sigma * y_scaled_sq - inv_sigma;
102 partials<1>(ops_partials) = std::move(scaled_diff);
105 return ops_partials.build(logp);
108template <
typename T_y,
typename T_loc,
typename T_scale>
111 const T_scale& sigma) {
112 return normal_lpdf<false>(y, mu, sigma);
require_all_not_t< is_nonscalar_prim_or_rev_kernel_expression< std::decay_t< Types > >... > require_all_not_nonscalar_prim_or_rev_kernel_expression_t
Require none of the types satisfy is_nonscalar_prim_or_rev_kernel_expression.
return_type_t< T_y_cl, T_loc_cl, T_scale_cl > normal_lpdf(const T_y_cl &y, const T_loc_cl &mu, const T_scale_cl &sigma)
The log of the normal density for the specified scalar(s) given the specified mean(s) and deviation(s...
size_t size(const T &m)
Returns the size (number of the elements) of a matrix_cl or var_value<matrix_cl<T>>.
typename return_type< Ts... >::type return_type_t
Convenience type for the return type of the specified template parameters.
size_t max_size(const T1 &x1, const Ts &... xs)
Calculate the size of the largest input.
bool size_zero(const T &x)
Returns 1 if input is of length 0, returns 0 otherwise.
fvar< T > log(const fvar< T > &x)
const double NEG_LOG_SQRT_TWO_PI
The value of minus the natural logarithm of the square root of , .
auto as_value_column_array_or_scalar(T &&a)
Extract the value from an object and for eigen vectors and std::vectors convert to an eigen column ar...
void check_consistent_sizes(const char *)
Trivial no input case, this function is a no-op.
fvar< T > sum(const std::vector< fvar< T > > &m)
Return the sum of the entries of the specified standard vector.
ref_type_t< T && > to_ref(T &&a)
This evaluates expensive Eigen expressions.
void check_finite(const char *function, const char *name, const T_y &y)
Return true if all values in y are finite.
void check_not_nan(const char *function, const char *name, const T_y &y)
Check if y is not NaN.
void check_positive(const char *function, const char *name, const T_y &y)
Check if y is positive.
fvar< T > inv(const fvar< T > &x)
auto make_partials_propagator(Ops &&... ops)
Construct an partials_propagator.
typename ref_type_if<!is_constant< T >::value, T >::type ref_type_if_not_constant_t
typename partials_return_type< Args... >::type partials_return_t
The lgamma implementation in stan-math is based on either the reentrant safe lgamma_r implementation ...
Extends std::true_type when instantiated with zero or more template parameters, all of which extend t...
Template metaprogram to calculate whether a summand needs to be included in a proportional (log) prob...