check_lower_triangular.hpp

Go to the documentation of this file.

#ifndef STAN_MATH_PRIM_ERR_CHECK_LOWER_TRIANGULAR_HPP
#define STAN_MATH_PRIM_ERR_CHECK_LOWER_TRIANGULAR_HPP
 
#include <stan/math/prim/meta.hpp>
#include <stan/math/prim/fun/Eigen.hpp>
#include <stan/math/prim/err/throw_domain_error.hpp>
#include <sstream>
#include <string>
 
namespace stan {
namespace math {
 
template <typename T_y, require_eigen_t<T_y>* = nullptr>
inline void check_lower_triangular(const char* function, const char* name,
                                   const T_y& y) {
  const auto& y_ref = to_ref(y);
  for (int n = 1; n < y.cols(); ++n) {
    for (int m = 0; m < n && m < y.rows(); ++m) {
      if (y_ref(m, n) != 0) {
        std::stringstream msg;
        msg << "is not lower triangular;"
            << " " << name << "[" << stan::error_index::value + m << ","
            << stan::error_index::value + n << "]=";
        std::string msg_str(msg.str());
        throw_domain_error(function, name, y_ref(m, n), msg_str.c_str());
      }
    }
  }
}
 
}  // namespace math
}  // namespace stan
#endif