栈上声明大数组会崩溃，必须转堆分配；因栈空间有限（Windows约1MB，Linux约8MB），超限导致运行时SIGSEGV；std::vector是默认推荐方案，兼顾安全与易用性。
栈上直接声明大数组会崩溃，必须转到堆上分配——这是最核心的判断，别犹豫，也别先想着调栈大小。
为什么 double phi[N][Nx][Ny] 一变大就崩
因为这是栈分配：编译器在函数进入时就预留连续内存，大小固定、不可扩展。Nx 和 Ny 到几千时，phi 很容易超过几 MB，默认栈（Windows 常为 1MB，Linux 通常 8MB）直接溢出，程序启动即 crash 或 SIGSEGV。这不是编译错误，是运行时栈溢出，gdb 里常看到 Segmentation fault (core dumped) 却找不到具体行——它发生在函数 prologue 阶段。
- 栈分配无异常机制，失败不抛
std::bad_alloc，而是直接终止 - 全局/静态声明虽能绕过栈限制，但污染命名空间、无法按需构造、难以管理生命周期
- VC 里改 Linker 的
/STACK选项只是临时止血，不同环境行为不一致，不可移植
std::vector 是最稳妥的默认选择
它把数据放堆上，自动管理容量增长和异常安全，且接口与原生数组接近。对三维数组这种场景，用嵌套 vector 虽稍慢，但清晰可靠：
#include <vector> std::vector<std::vector<std::vector<double>>> phi(N, std::vector<std::vector<double>>(Nx, std::vector<double>(Ny)));
- 优点：自动初始化、边界检查（
.at()）、移动语义支持、RAII 安全释放 - 缺点：三层指针跳转，缓存局部性差；每层 vector 有额外 24 字节开销（指针+size+capacity）
- 若性能敏感，改用一维展平：
std::vector<double> phi(N * Nx * Ny);，访问用phi[i * Nx * Ny + j * Ny + k]
手动 new[] / delete[] 要小心三件事
它比 vector 轻量，但错误成本高，尤其涉及类类型或异常路径时：
立即学习“C++免费学习笔记（深入）”；
- 必须配对使用
new[]和delete[]—— 用delete会只析构首元素，其余内存泄漏 - 内置类型（如
double）可直接new double[Nx*Ny*N]，但自定义类型必须显式调用构造函数（placement new），否则对象未构造 - 若分配途中抛异常（如内存不足），
new已抛std::bad_alloc，但你没机会清理前面已分配的块 —— 必须用 try/catch 或 RAII 封装
超大规模（GB 级）或高频小对象，考虑内存池
当数组不是单一大块，而是成千上万个中等大小数组（比如每个 64KB 的粒子缓冲区），new 的锁竞争和元数据开销会成为瓶颈。此时标准分配器失效，需定制：
- 定长内存池：预分配大块内存，切分为固定尺寸 slot，分配/释放 O(1)，零碎片，适合同构数组池
- 避免用
std::allocator直接替换 vector —— 它只管内存，不处理对象构造；必须配合construct/destroy手动调用 - 生产环境慎用自研池：先 profile 确认
malloc真是瓶颈（用perf record -e syscalls:sys_enter_brk查系统调用频次）
真正麻烦的从来不是“怎么分配”，而是“谁负责释放”和“中间出异常怎么办”。vector 把这两点封死了；裸指针把它们全摊开在你面前。选哪个，取决于你愿不愿意为那点性能多写二十行防御代码。
