On Thu, 19 Aug 2021 at 04:53, Dirk Eddelbuettel <e...@debian.org> wrote: > > > Naeem, > > I would simplify, simplify, simplify -- as 'Rcpp FAQ 7.31' reminds us all, > testing _equality_ of doubles is challenging anyway. > > Besides, it may make sense to would ascertain first you get what you want in > _purely serial modes_ and then move to OpenMP.
Exactly. Serial execution should be fine. I.e., if you set the number of threads to 1, then all platforms will return the same result. However, you have defined a number of variables outside the parallel region, and then you modify them inside the parallel region. OpenMP takes those variables as shared by default, which leads to the unexpected results you are getting. You need to tell OpenMP that those variables are threadprivate. Or you could just define them inside the parallel region, so that OpenMP knows that they are private without additional hints. -- Iñaki Úcar _______________________________________________ Rcpp-devel mailing list Rcpp-devel@lists.r-forge.r-project.org https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel