在我的大学使用计算机集群时出现矩阵错误,我无法在本地计算机上重现。我认为这可能是由于矩阵库(BLAS、LAPACK、ATLAS 等)的不同。除了我在这里阅读的内容之外,我对这些库了解不多,但我认为应该有一种方法让 R 告诉我它安装了哪些矩阵库(即它正在使用哪些矩阵库),类似到 sessionInfo() 告诉我它正在使用哪个版本的 R 包。
3 回答
尽管可能没有明确的 R 函数,但也许您可以利用 shell 命令(例如file
和ldd
)来获得一些线索,而无需重建 R 或集群上的 root 访问权限:
(rpath <- Sys.which("R"))
# R
# "/usr/bin/R"
为确保它ldd
能够正常工作,请查看它是什么类型的文件:
system2("file", rpath)
# /usr/bin/R: Bourne-Again shell script, ASCII text executable
如果您的显示一个实际的可执行文件(例如ELF 64-bit LSB executable, x86-64, ...
),则跳过这一步。
script <- readLines(rpath)
script[grepl("/bin/", script)]
# [1] "#!/bin/bash"
# [2] " if [ -x \"/usr/${libnn}/R/bin/exec/R\" ]; then"
# [3] " elif [ -x \"/usr/${libnn_fallback}/R/bin/exec/R\" ]; then"
# [4] "## some systems have a more portable sed, e.g. /usr/xpg4/bin/sed on Solaris,"
# [5] "SED=/bin/sed"
# [6] " exec sh \"${R_HOME}/bin/Rcmd\" \"${@}\" ;;"
# [7] "R_binary=\"${R_HOME}/bin/exec${R_ARCH}/R\""
这告诉我实际的可执行文件是/usr/lib/R/bin/exec/R
(或/usr/lib64/...
)。它正在采取一些推论,但这是一个步骤。这对我在 ubuntu 上使用 R-3.3.2 有效,所以我只能假设它在不同的操作系统上是相似的。如果这没有提供信息,您还可以 grep for "/lib"
or "/exec"
(或者只是检查整个脚本文件以获取其他线索)。
一旦你知道脚本调用的路径(或者如果它不是开始的脚本),找到共享库依赖项:
system2("ldd", "/usr/lib/R/bin/exec/R")
# linux-vdso.so.1 => (0x00007ffcfadcd000)
# libR.so => /usr/lib/R/lib/libR.so (0x00007f746f615000)
# libpthread.so.0 => /lib/x86_64-linux-gnu/libpthread.so.0 (0x00007f746f3eb000)
# libc.so.6 => /lib/x86_64-linux-gnu/libc.so.6 (0x00007f746f025000)
# libblas.so.3 => /usr/lib/libblas.so.3 (0x00007f746eda8000)
# libm.so.6 => /lib/x86_64-linux-gnu/libm.so.6 (0x00007f746eaa2000)
# libreadline.so.6 => /lib/x86_64-linux-gnu/libreadline.so.6 (0x00007f746e85b000)
# libpcre.so.3 => /lib/x86_64-linux-gnu/libpcre.so.3 (0x00007f746e61d000)
# liblzma.so.5 => /lib/x86_64-linux-gnu/liblzma.so.5 (0x00007f746e3fb000)
# libbz2.so.1.0 => /lib/x86_64-linux-gnu/libbz2.so.1.0 (0x00007f746e1ea000)
# libz.so.1 => /lib/x86_64-linux-gnu/libz.so.1 (0x00007f746dfd1000)
# librt.so.1 => /lib/x86_64-linux-gnu/librt.so.1 (0x00007f746ddc9000)
# libdl.so.2 => /lib/x86_64-linux-gnu/libdl.so.2 (0x00007f746dbc4000)
# libgomp.so.1 => /usr/lib/x86_64-linux-gnu/libgomp.so.1 (0x00007f746d9b5000)
# /lib64/ld-linux-x86-64.so.2 (0x0000560abd5fa000)
# libtinfo.so.5 => /lib/x86_64-linux-gnu/libtinfo.so.5 (0x00007f746d78b000)
此输出的第 5 行表明 R 正在使用 BLAS 库,特别是libblas3
.
据我所知,没有一个特定的函数(尤其是在 base R 中没有)可以为您提供 BLAS 版本。R 很可能会使用系统附带的标准库(除非有人使用特定 BLAS 版本从源代码构建它)。
您必须从终端执行此操作,输入
locate libblas.so
它可能会在/usr/lib/libblas.so
or/usr/lib64/libblas.so
中,并且该位置可能是一个符号链接。所以按照链接
ls -l /usr/lib/libblas.so
如果要更改符号链接以指向不同的 BLAS 库,则可以执行
update-alternatives --config libblas.so.3gf
并选择要使用的版本。这篇文章非常详细地解释了它:https ://www.r-bloggers.com/for-faster-r-use-openblas-instead-better-than-atlas-trivial-to-switch-to-on- ubuntu/
另一种确定您正在使用的版本的方法是安装您想要的任何 BLAS 实现——例如 OpenBLAS,然后用指向新库的符号链接替换标准 R BLAS 库。确保使用--with-shared-blas
选项安装 R,如下所述:https ://cran.r-project.org/doc/manuals/r-release/R-admin.html#BLAS
然后你只需要:
# navigate to your R directory
cd …/path-to-R-installation/lib/R/lib
# backup the original
mv libRblas.so libRblas.backup
# create simlink
ln –s …/path-to-new-0blas-library/new-lib.so libRblas.so
谷歌搜索就够了吗?: http: //www.stat.cmu.edu/~nmv/2013/07/09/for-faster-r-use-openblas-instead-better-than-atlas-trivial-to -切换到 ubuntu/