Stack imbalance with RcppParallel

Question

I wrote the following code to train myself to use RcppParallel. It is just a toy example.

// [[Rcpp::depends(RcppParallel)]]
#include 
#include         


        

Answer 1

Look more closely at an example such as the parallel distance from the Rcpp Gallery.

It does not use NumericMatrix but rather RMatrix<double>. I would do the same here, and have generally advocated not to rely on contact with R types while running parallel segments.