使用boost::odeint搭配std::array出现EXC_BAD_ACCESS栈溢出问题

问题根因 #

错误类型为栈溢出（stack-overflow），核心差异来自std::array和std::vector的内存分配机制不同：

• std::vector的元素存储在堆内存，只有3个指针大小的控制结构存储在栈上，即使存储1万个double，栈占用也只有几十字节，不会有栈溢出风险。
• std::array的所有元素直接存储在栈上，单个std::array<double, 10000>的大小为 10000 * 8字节 = 80KB，而boost::odeint的Runge-Kutta dopri5步进器内部会持有至少5_{7个中间状态数组（用于存储RK算法的k1}k6中间量、误差估计值等），叠加后总栈占用会达到数百KB，再加上AddressSanitizer的插桩逻辑会额外占用大量栈空间、压缩可用栈上限，最终触发栈溢出。

当n_variables降到1000时，单个std::array的大小降到8KB，总栈占用远小于可用上限，因此问题消失。

解决方案 #

方案1（推荐）：继续使用std::vector作为状态类型 #

这是兼容性最好、无额外开发成本的方案，已验证可以正常运行，堆内存分配足以支撑大维度的状态存储需求。

方案2：强制大对象堆分配 #

如果一定要使用std::array，可以将步进器、ODE实例等包含大std::array成员的对象改为堆分配，避免占用栈空间：

// 示例：堆分配ODE2实例和stepper
auto o = std::make_unique<ODE2>(-60.0);
auto stepper = std::make_unique<decltype(make_controlled<error_stepper_type>(abs_err, rel_err))>(abs_err, rel_err);
integrate_adaptive(*stepper, *o, o->GetVoltage(), t_start, t_stop, dt);

方案3：手动增大栈空间（不推荐） #

macOS下可以在编译时添加链接参数调整栈大小，以下参数将栈上限设为16MB：

c++ -O0 -g -fsanitize=address -fno-omit-frame-pointer -Wall -Wextra -Wpedantic -std=c++17 ode.cpp -o ode -Wl,-stack_size,0x1000000

该方案可移植性极差，不同操作系统、不同运行环境的默认栈大小规则不同，不建议在生产环境使用。

内容的提问来源于stack exchange，提问作者bhache