Stan Math Library  2.12.0
reverse mode automatic differentiation
check_ordered.hpp
Go to the documentation of this file.
1 #ifndef STAN_MATH_PRIM_MAT_ERR_CHECK_ORDERED_HPP
2 #define STAN_MATH_PRIM_MAT_ERR_CHECK_ORDERED_HPP
3 
8 #include <sstream>
9 #include <string>
10 #include <vector>
11 
12 namespace stan {
13  namespace math {
14 
30  template <typename T_y>
31  bool check_ordered(const char* function,
32  const char* name,
33  const Eigen::Matrix<T_y, Eigen::Dynamic, 1>& y) {
34  using Eigen::Dynamic;
35  using Eigen::Matrix;
36 
37  typedef typename index_type<Matrix<T_y, Dynamic, 1> >::type size_t;
38 
39  if (y.size() == 0)
40  return true;
41 
42  for (size_t n = 1; n < y.size(); n++) {
43  if (!(y[n] > y[n-1])) {
44  std::ostringstream msg1;
45  msg1 << "is not a valid ordered vector."
46  << " The element at " << stan::error_index::value + n
47  << " is ";
48  std::string msg1_str(msg1.str());
49  std::ostringstream msg2;
50  msg2 << ", but should be greater than the previous element, "
51  << y[n-1];
52  std::string msg2_str(msg2.str());
53  domain_error(function, name, y[n],
54  msg1_str.c_str(), msg2_str.c_str());
55  return false;
56  }
57  }
58  return true;
59  }
60 
61  }
62 }
63 #endif
Primary template class for the metaprogram to compute the index type of a container.
Definition: index_type.hpp:18
bool check_ordered(const char *function, const char *name, const std::vector< T_y > &y)
Return true if the specified vector is sorted into strictly increasing order.
void domain_error(const char *function, const char *name, const T &y, const char *msg1, const char *msg2)
Throw a domain error with a consistently formatted message.

     [ Stan Home Page ] © 2011–2016, Stan Development Team.