TEST(MathFunctions, ones_vector) {
  for (int K = 0; K < 5; K++) {
    Eigen::VectorXd y = Eigen::VectorXd::Constant(K, 1);
    EXPECT_MATRIX_FLOAT_EQ(y, stan::math::ones_vector(K));
  }
}

TEST(MathFunctions, ones_vector_throw) {
  EXPECT_THROW(stan::math::ones_vector(-1), std::domain_error);
}
