Can anyone spot the problem with calling this exported function auto nonstiff(Eigen::MatrixXd A, Eigen::VectorXd y0) { using stan::math::integrate_ode_rk45; using stan::math::to_vector; using stan::math::to_array_1d; std::vector<double> theta; auto y = integrate_ode_rk45([&A](auto t, auto…

It doesn’t take ode or matrix_exp to crash: auto nonstiff(Eigen::MatrixXd A, Eigen::VectorXd y0) { std::vector<std::vector<double>> y(1, {2, 2}); return stan::math::to_vector(y[0]) - A * y0; } TEST(ode, debug) { Eigen::MatrixXd A(2,2); A << 1, 1, 1, 1; Eigen::VectorXd y0(2); y0 << 1, …

[image] yizhang: the run passes when the return of nonstiff is specified as VectorXd . So my guess is it has to do with auto type deduction of Eigen objects. I’m always amazed at what people are able to debug. I wonder what it thought the return type was (not that it matters at this point).

Heap use after free

Developers

bgoodri July 19, 2019, 11:12pm 3

Thanks

Topic		Replies	Views
Undefined behavior in ODE solver Developers rstan	22	1088	January 23, 2020
Am I Leaking? [Eigen autodiff] Developers	12	2012	April 29, 2019
Refactoring math with more general signatures Developers math	7	988	October 25, 2019
ODEs: Cost of checking Developers	28	1931	May 12, 2017
The Math Library is not passing tests on Windows. This is the next thing I'm working on Developers math	22	1441	November 30, 2018

Heap use after free

Related topics