当Valgrind未指向自定义代码时，如何修复R语言内存泄漏？

问题场景 #

我正在开发fixest2包，在标准误计算环节，执行以下代码的第二个循环时出现内存问题：

library(fixest2)

lhs = 'y1'
rhs = 'x3'
all_rhs = c('', 'x2', 'x3')
n_rhs = 1
k = 1

base = readRDS("dev/base.rds")

# 重复执行以触发内存泄漏（无问题）
for (i in 1:10) {
  est_multi <- feols(c(y1, y2) ~ x1 + csw0(x2, x3) + x4 | species + fe2, base, fsplit = ~species)
}

# 重复执行触发内存问题的代码块
for (i in 1:10) {
  mat <- coeftable(object = est_multi[[1]],
      vcov = NULL, ssc = NULL, cluster = "fe3",
      keep = NULL, drop = NULL, order = NULL)

  mat[, 2]
}

但Valgrind的输出仅指向R内部代码，无法定位到我自己的代码：

...
==52939==  Uninitialised value was created by a heap allocation
==52939==    at 0x4848899: malloc (in /usr/libexec/valgrind/vgpreload_memcheck-amd64-linux.so)
==52939==    by 0x49F0E14: ??? (in /usr/lib/R/lib/libR.so)
...

使用的调试启动命令：

R --vanilla -d 'valgrind -s --track-origins=yes'

排查方案 #

1. 编译包时开启调试符号与无优化编译 #

在包的Makevars（Linux/macOS）或Makevars.win（Windows）文件中添加以下编译参数，确保Valgrind能关联到你的代码行号：

CXXFLAGS += -g -O0 -Wall -Wextra
CXX11FLAGS += -g -O0 -Wall -Wextra

• -g：生成调试符号，让Valgrind能映射到源代码行
• -O0：关闭编译器优化，避免代码行号被打乱
• -Wall/-Wextra：开启编译警告，提前发现未初始化变量、内存越界等潜在问题

2. 逐步缩小代码范围，锁定触发点 #

• 简化第二个循环的代码：先单独调用coeftable并去掉cluster参数，看问题是否消失；再逐步添加参数，定位到具体触发问题的参数或逻辑
• 检查est_multi[[1]]的生命周期：确认该对象在循环中是否被正确覆盖或回收，是否存在未释放的内部资源
• 替换cluster = "fe3"为其他聚类变量或值，验证是否是聚类标准误计算逻辑导致的问题

3. 用R内置工具辅助排查 #

• 在循环内添加gc()手动触发垃圾回收，观察内存变化，判断是否是对象未被正确回收
• 用tracemem(est_multi[[1]])跟踪对象的内存引用链，查看是否有意外引用导致对象无法被GC回收
• 使用profmem包记录内存分配细节，定位哪段代码分配了大量内存且未释放

4. 调整Valgrind参数获取更详细信息 #

更新启动命令，添加更全面的泄漏检查参数：

R --vanilla -d 'valgrind -s --track-origins=yes --leak-check=full --show-leak-kinds=all'

如果仍指向R内部，大概率是你的代码传递了未初始化的变量给R内置函数，导致R触发内存问题。此时重点检查coeftable中处理cluster参数的代码，确保所有传递给R API的变量都已正确初始化。

5. 检查C/C++后端的内存管理 #

如果coeftable包含C/C++实现：

• 核对所有malloc/free、new/delete的配对，确保没有遗漏释放
• 检查是否存在访问已释放内存的野指针问题
• 确认所有变量在使用前都已完成初始化

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