我在使用MKL库触发R的线程时遇到问题。我目前正在使用Ubuntu 18.04.2 LTS。
Linux pedro-HP-EliteOne-800-G1-AiO 4.18.0-15-generic #16~18.04.1-Ubuntu SMP Thu Feb 7 14:06:04 UTC 2019 x86_64 x86_64 x86_64 GNU/Linux
r-open$ screenfetch
./+o+- pedro@pedro-HP-EliteOne-800-G1-AiO
yyyyy- -yyyyyy+ OS: Ubuntu 18.04 bionic
://+//////-yyyyyyo Kernel: x86_64 Linux 4.18.0-15-generic
.++ .:/++++++/-.+sss/` Uptime: 19h 27m
.:++o: /++++++++/:--:/- Packages: 1834
o:+o+:++.`..```.-/oo+++++/ Shell: bash
.:+o:+o/. `+sssoo+/ Resolution: 1920x1080
.++/+:+oo+o:` /sssooo. DE: GNOME
/+++//+:`oo+o /::--:. WM: GNOME Shell
\+/+o+++`o++o ++////. WM Theme: Adwaita
.++.o+++oo+:` /dddhhh. GTK Theme: Ambiance [GTK2/3]
.+.o+oo:. `oddhhhh+ Icon Theme: ubuntu-mono-dark
\+.++o+o``-````.:ohdhhhhh+ Font: Ubuntu 11
`:o+++ `ohhhhhhhhyo++os: CPU: Intel Core i5-4590S @ 4x 3.7GHz [27.8°C]
.o:`.syhhhhhhh/.oo++o` GPU: intel
/osyyyyyyo++ooo+++/ RAM: 2659MiB / 7867MiB
````` +oo+++o\:
`oo++.
注意:我的机器有4个核心。
最初我在intel网站上下载了文件l_mkl_2019.2.187.tgz,并为/home/pedro/
进行了安装。后来我下载了R-3.5.2.tar.gz
文件来编译R,支持多线程MKL。我通过执行以下操作来执行编译:
source /home/pedro/intel/mkl/bin/mklvars.sh intel64
cd /home/pedro/Downloads/R-3.5.2 && ./configure --enable-R-shlib --enable-threads=posix --with-lapack --with-blas="-fopenmp -m64 -I$MKLROOT/include -L$MKLROOT/lib/intel64 -lmkl_gf_lp64 -lmkl_gnu_thread -lmkl_core -lpthread -lm"
make
make install
显然编译没有任何问题。使用sessionInfo()
命令我得到以下结果:
> sessionInfo()
R version 3.5.2 (2018-12-20)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 18.04.2 LTS
Matrix products: default
BLAS/LAPACK: /home/pedro/intel/compilers_and_libraries_2019.2.187/linux/mkl/lib/intel64_lin/libmkl_gf_lp64.so
locale:
[1] LC_CTYPE=pt_BR.UTF-8 LC_NUMERIC=C
[3] LC_TIME=pt_BR.UTF-8 LC_COLLATE=pt_BR.UTF-8
[5] LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=pt_BR.UTF-8
[7] LC_PAPER=pt_BR.UTF-8 LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
loaded via a namespace (and not attached):
[1] compiler_3.5.2
基本上我正在尝试使用parallel :: mclapply
函数并行化蒙特卡罗模拟。事实上,串行(1核)我使用MKL有更好的结果,但在尝试并行化时,我意识到线程不会在4核中被触发。
我使用parallel::mclapply
函数做的是:
parallel::mclapply(1:8, FUN = function(x) func_metodos(), mc.cores = 4))
,其中func_metodos()
是我在R中实现的并且想要重复8次的函数。
注意:除了使用Microsoft R Open之外,我想要一个解决方案。我不喜欢使用微软所关注的东西。
最好的祝福。