R中Rcpp导出的C++函数循环调用结果逐次异常累加

该异常是Armadillo容器未初始化内存导致的典型累积错误，完全匹配观测到的逐次倍增行为：

• 问题根源在函数内部的向量初始化代码：arma::colvec N(K);
  按照Armadillo的内存分配规则，仅传入维度参数构造向量/矩阵时，系统仅申请对应大小的内存空间，不会自动将内存中的残留值清零。
• 倍增现象的形成链路：
  1. 首次调用Demo函数时，为N分配的内存区域此前未被其他程序写入，残留值恰好为0，因此累加计算后得到正确结果。函数返回时，这块大小为K的colvec对应的内存被标记为可复用，但内存中存储的计算结果不会被主动擦除。
  2. 第二次、第三次循环调用时，系统为新的同长度colvec申请内存，会优先复用刚被释放、大小完全匹配的那块内存。此时新的N向量读取到的就是上一次函数计算得到的结果值。
  3. 代码内的累加逻辑为N[j] += accu(...)，相当于在残留的旧值基础上叠加本次的计算结果。因为输入矩阵n未发生变化，每次叠加的增量是固定的正确结果值，因此第二次输出为正确值的2倍，第三次输出为正确值的3倍。
     整个过程中R环境下的输入矩阵n、参数K因为始终被R引用、不会被释放复用，所以全程保持初始值不变，与观测结果完全一致。

• 修改向量初始化逻辑，构造N时显式指定用0值填充内存，将原代码：

  arma::colvec N(K);
  

  替换为：

  arma::colvec N(K, arma::fill::zeros);
  

  修改后N的所有元素默认值为0，从根源上避免读取内存残留的旧值，多次循环调用时会稳定返回相同的正确结果。
• 通用编写注意事项：
  使用Armadillo容器时，只要涉及逐元素累加、增量计算的逻辑，初始化时必须显式指定填充规则，不要依赖未初始化内存的默认值。这类bug通常在单次调用时因为内存恰好为0被掩盖，只有在多次重复调用、触发内存复用的场景下才会暴露。

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