Why isn't `ts` in `integrate_ode` using `ordered` type?


#1

Is there any reason we are not mandating it be ordered type?


#2

ordered type is a constrained type for block-level variables, but you can pass local variables to this function as well.


#3

Also, the times in the integrate_ode functions are real arrays (real[]). The ordered type maps to vector.

In C++, the current implementation takes std::vector and not Eigen::Matrix.


#4

I see. Thanks.

Since we are checking ts to be ordered anyway, I’d think using ordered type is more straightforward.


#5

It’s not something that can be enforced. Bob’s written it up in Section 3.1, “Overview of Data Types.” (in the v2.17.0 version of the manual). Here’s the relevant sentence:

Arguments for built-in and user-defined functions and local variables are required to be basic data types, meaning an unconstrained primitive, vector, or matrix type or an array of such.

In this context, constrained vectors are things like ordered or simplex, and unconstrained vector is vector. There’s a little more information in Chapter 3 that might be useful to understand. The manual wasn’t designed to describe how the Stan types map to C++ types, but hopefully understanding how it maps makes this section a little clearer.


#6

Yes. Your quote really makes it clear. Thanks.