R在Fedora下矩阵乘法异常占用全部CPU核心的求助

问题根源 #

Fedora默认打包的R会链接多线程线性代数库（如OpenBLAS或MKL），这类库会自动将矩阵运算（如%*%）分配到所有CPU核心执行。而Windows/Ubuntu Server默认使用的是单线程BLAS实现，因此不会出现此问题。

解决方案 #

1. 临时设置单线程（单次R会话生效） #

在R启动后，根据当前使用的BLAS库设置线程数：

• 如果是OpenBLAS：

Sys.setenv(OPENBLAS_NUM_THREADS = 1)

• 如果是MKL：

Sys.setenv(MKL_NUM_THREADS = 1)

验证当前BLAS库：在终端执行以下命令查看R依赖的线性代数库：

ldd $(which R) | grep -E '(openblas|mkl)'

2. 永久设置（所有R会话生效） #

方式一：修改R启动文件

编辑用户目录下的.Rprofile文件（不存在则创建），添加对应环境变量设置：

# 针对OpenBLAS
Sys.setenv(OPENBLAS_NUM_THREADS = 1)
# 若使用MKL则替换为：
# Sys.setenv(MKL_NUM_THREADS = 1)

方式二：系统层面设置环境变量

创建/etc/profile.d/r-blas.sh文件，添加：

# 针对OpenBLAS
export OPENBLAS_NUM_THREADS=1
# 针对MKL则替换为：
# export MKL_NUM_THREADS=1

保存后重启终端或执行source /etc/profile.d/r-blas.sh生效。

3. 适配手动并行场景 #

使用foreach等手动并行框架时，需确保每个工作线程都限制BLAS线程数，避免嵌套并行导致CPU过载：

library(doParallel)
library(foreach)
N <- 10000
A <- matrix(0.1, nrow = N, ncol = N)
B <- matrix(0.1, nrow = N, ncol = N)

cl <- makeCluster(4)
# 给所有工作线程设置单线程BLAS
clusterEvalQ(cl, Sys.setenv(OPENBLAS_NUM_THREADS = 1))
registerDoParallel(cl)

res <- foreach(i = 1:10) %dopar% {
  AB <- A %*% B
  return(list(i = i, AB = AB))
}

stopCluster(cl)

4. 重新编译R（备选方案） #

若上述方法无效，可重新编译R并指定链接单线程BLAS库。编译时添加参数：

./configure --with-blas=-lblas --with-lapack=-llapack

后续按常规编译流程完成安装。

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