在OS X上剖析Rcpp代码

Question

我有兴趣分析OS X（Mountain Lion 10.8.2）下的一些Rcpp代码，但剖析器在运行时会崩溃。

玩具示例，使用inline，只是设计为花费足够的时间让探查器注意。

library(Rcpp) 
library(inline) 

src.cpp <- " 
    RNGScope scope; 
    int n = as<int>(n_); 
    double x = 0.0; 
    for (int i = 0; i < n; i++) 
    x += (unif_rand()-.5); 
    return wrap(x);" 

src.c <- " 
    int i, n = INTEGER(n_)[0]; 
    double x = 0.0; 
    GetRNGstate(); 
    for (i = 0; i < n; i++) 
    x += (unif_rand()-.5); 
    PutRNGstate(); 
    return ScalarReal(x);" 

f.cpp <- cxxfunction(signature(n_="integer"), src.cpp, plugin="Rcpp") 
f.c <- cfunction(signature(n_="integer"), src.c) 

如果我使用GUI仪器（在Xcode，版本4.5（4523））或命令行sample，两者碰撞：器械马上崩溃，而样品崩溃之前完成处理样品：

# (in R) 
set.seed(1) 
f.cpp(200000000L) 

# (in a separate terminal window) 
~ » sample R # this invokes the profiler 
Sampling process 81337 for 10 seconds with 1 millisecond of run time between samples 
Sampling completed, processing symbols... 
[1] 81654 segmentation fault sample 81337 

如果我做同样的过程，但与C版本（即f.c(200000000L)）这两项文书和sample做工精细，并产生像

Call graph: 
1832 Thread_6890779 DispatchQueue_1: com.apple.main-thread (serial) 
    1832 start (in R) + 52 [0x100000e74] 
    1832 main (in R) + 27 [0x100000eeb] 
     1832 run_Rmainloop (in libR.dylib) + 80 [0x1000e4020] 
     1832 R_ReplConsole (in libR.dylib) + 161 [0x1000e3b11] 
      1832 Rf_ReplIteration (in libR.dylib) + 514 [0x1000e3822] 
      1832 Rf_eval (in libR.dylib) + 1010 [0x1000aa402] 
       1832 Rf_applyClosure (in libR.dylib) + 849 [0x1000af5d1] 
       1832 Rf_eval (in libR.dylib) + 1672 [0x1000aa698] 
        1832 do_dotcall (in libR.dylib) + 16315 [0x10007af3b] 
        1382 file1412f6e212474 (in file1412f6e212474.so) + 53 [0x1007fded5] file1412f6e212474.cpp:16 
        + 862 unif_rand (in libR.dylib) + 1127,1099,... [0x10000b057,0x10000b03b,...] 
        + 520 fixup (in libR.dylib) + 39,67,... [0x10000aab7,0x10000aad3,...] 
        356 file1412f6e212474 (in file1412f6e212474.so) + 70,61,... [0x1007fdee6,0x1007fdedd,...] file1412f6e212474.cpp:16 
        56 unif_rand (in libR.dylib) + 1133 [0x10000b05d] 
        38 DYLD-STUB$$unif_rand (in file1412f6e212474.so) + 0 [0x1007fdf1c] 

输出

我真的很感激一些建议，如果有什么我做错了，如果有一些其他的首选方式，或者如果这是不可能的。鉴于Rcpp的主要用途之一似乎是加速R代码，我很惊讶没有找到更多的信息，但也许我看错了地方。

这是在OS X 10.8.2上，R 2.15.1（x86_64-apple-darwin9.8.0），Rcpp 0.9.15和g ++ --version报告“i686-apple-darwin11-llvm-g ++ - 4.2 （GCC）4.2.1（基于Apple Inc. build 5658）（LLVM build 2336.11.00）“。

溶液

得益于以下德克的回答，和他的谈话在这里http://dirk.eddelbuettel.com/papers/ismNov2009introHPCwithR.pdf，我至少使用谷歌perftools的部分解决方案。首先，从这里安装http://code.google.com/p/gperftools/，并在编译C++代码时将-lprofiler添加到PKG_LIBS。然后要么

（a）运行R作为CPUPROFILE=samples.log R，运行所有代码和退出（或使用RSCRIPT）

（b）利用两个小的效用函数，以打开/关闭分析：

RcppExport SEXP start_profiler(SEXP str) { 
    ProfilerStart(as<const char*>(str)); 
    return R_NilValue; 
} 

RcppExport SEXP stop_profiler() { 
    ProfilerStop(); 
    return R_NilValue; 
} 

然后，R内，你可以做

.Call("start_profiler", "samples.log") 
# code that calls C++ code to be profiled 
.Call("stop_profiler") 

无论哪种方式，文件samples.log将包含分析信息。这可以在与

pprof --text /Library/Frameworks/R.framework/Resources/bin/exec/x86_64/R samples.log 

产生像

Using local file /Library/Frameworks/R.framework/Resources/bin/exec/x86_64/R. 
Using local file samples.log. 
Removing __sigtramp from all stack traces. 
Total: 112 samples 
    64 57.1% 57.1%  64 57.1% _unif_rand 
    30 26.8% 83.9%  30 26.8% _process_system_Renviron 
    14 12.5% 96.4%  101 90.2% _for_profile 
    3 2.7% 99.1%  3 2.7% Rcpp::internal::expr_eval_methods 
    1 0.9% 100.0%  1 0.9% _Rf_PrintValueRec 
    0 0.0% 100.0%  1 0.9% 0x0000000102bbc1ff 
    0 0.0% 100.0%  15 13.4% 0x00007fff5fbfe06f 
    0 0.0% 100.0%  1 0.9% _Rf_InitFunctionHashing 
    0 0.0% 100.0%  1 0.9% _Rf_PrintValueEnv 
    0 0.0% 100.0%  112 100.0% _Rf_ReplIteration 

输出，很可能是在一个真实的例子更多的信息来看。

Answer 1

我很困惑，你的例子是不完整的：

• 你不显示的cfunction()的（简单）调用和cxxfunction()

• 你不告诉你如何调用分析器

• 您没有概要分析的C或C++代码（!!）

你可以编辑问题并使其更清楚吗？

另外，当我运行这个时，这两个例子确实给出了相同的速度结果，因为它们本质上是相同的。 [Rcpp会让你用糖随机数函数调用。 ]

R> library(Rcpp) 
R> library(inline) 
R> 
R> src.cpp <- " 
+ RNGScope scope; 
+ int n = as<int>(n_); 
+ double x = 0.0; 
+ for (int i = 0; i < n; i++) 
+  x += (unif_rand()-.5); 
+ return wrap(x);" 
R> 
R> src.c <- " 
+ int i, n = INTEGER(n_)[0]; 
+ double x = 0.0; 
+ GetRNGstate(); 
+ for (i = 0; i < n; i++) 
+  x += (unif_rand()-.5); 
+ PutRNGstate(); 
+ return Rf_ScalarReal(x);" 
R> 
R> fc <- cfunction(signature(n_="int"), body=src.c) 
R> fcpp <- cxxfunction(signature(n_="int"), body=src.c, plugin="Rcpp") 
R> 
R> library(rbenchmark) 
R> 
R> print(benchmark(fc(10000L), fcpp(10000L))) 
     test replications elapsed relative user.self sys.self user.child sys.child 
1 fc(10000)   100 0.013  1  0.012  0   0   0 
2 fcpp(10000)   100 0.013  1  0.012  0   0   0 
R>