如何确定我的R安装使用的矩阵库?

时间:2015-03-31 15:45:39

标签: r build lapack blas atlas

在我的大学使用我无法在本地计算机上重现的计算机群集时出现矩阵错误。我认为这可能是由于矩阵库(BLAS,LAPACK,ATLAS等)的不同。除了我读过的here之外,我对这些库了解不多,但我认为应该有办法让R告诉我它安装了哪些矩阵库(即它是哪些库)使用),类似于sessionInfo()告诉我它正在使用哪个版本的R包。

3 个答案:

答案 0 :(得分:5)

虽然可能没有明确的R函数,但也许你可以利用shell命令(例如fileldd)获得一些线索,而不需要重建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 "/lib""/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

答案 1 :(得分:4)

据我所知,没有特定的功能(特别是不在基础R中)会给你BLAS版本。 R很可能会使用系统附带的标准库(除非有人使用特定的BLAS版本从源代码构建它)。

您必须从终端进行操作,输入

locate libblas.so

它可能会在/usr/lib/libblas.so/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

答案 2 :(得分:-3)