1 #ifndef STAN_MATH_PRIM_MAT_FUN_CSR_MATRIX_TIMES_VECTOR_HPP
2 #define STAN_MATH_PRIM_MAT_FUN_CSR_MATRIX_TIMES_VECTOR_HPP
10 #include <boost/math/tools/promotion.hpp>
75 template <
typename T1,
typename T2>
77 Eigen::Matrix<typename boost::math::tools::promote_args<T1, T2>::type,
81 const Eigen::Matrix<T1, Eigen::Dynamic, 1>& w,
82 const std::vector<int>& v,
83 const std::vector<int>& u,
84 const Eigen::Matrix<T2, Eigen::Dynamic, 1>& b) {
85 typedef typename boost::math::tools::promote_args<T1, T2>::type
94 u[m - 1] +
csr_u_to_z(u, m - 1) - 1,
"v", v.size());
95 for (
unsigned int i = 0; i < v.size(); ++i)
96 check_range(
"csr_matrix_times_vector",
"v[]", n, v[i]);
98 Eigen::Matrix<result_t, Eigen::Dynamic, 1> result(m);
104 Eigen::Matrix<result_t, Eigen::Dynamic, 1> b_sub(idx);
107 nze < row_end_in_w; ++nze, ++i) {
108 check_range(
"csr_matrix_times_vector",
"j", n, v[nze]);
111 Eigen::Matrix<T1, Eigen::Dynamic, 1>
bool check_range(const char *function, const char *name, const int max, const int index, const int nested_level, const char *error_msg)
Return true if specified index is within range.
Eigen::Matrix< T, 1, Eigen::Dynamic > row(const Eigen::Matrix< T, Eigen::Dynamic, Eigen::Dynamic > &m, size_t i)
Return the specified row of the specified matrix, using start-at-1 indexing.
bool check_positive(const char *function, const char *name, const T_y &y)
Return true if y is positive.
bool check_size_match(const char *function, const char *name_i, T_size1 i, const char *name_j, T_size2 j)
Return true if the provided sizes match.
fvar< T > dot_product(const Eigen::Matrix< fvar< T >, R1, C1 > &v1, const Eigen::Matrix< fvar< T >, R2, C2 > &v2)
Eigen::Matrix< typename boost::math::tools::promote_args< T1, T2 >::type, Eigen::Dynamic, 1 > csr_matrix_times_vector(const int &m, const int &n, const Eigen::Matrix< T1, Eigen::Dynamic, 1 > &w, const std::vector< int > &v, const std::vector< int > &u, const Eigen::Matrix< T2, Eigen::Dynamic, 1 > &b)