Re: [Rcpp-devel] Numerical precision in rotations with Eigen

2024-01-10 Thread Rafael Ayala Hernandez
Dear Dirk, Thanks a lot. This has indeed been quite insightful, and the link to the classic you provided makes me realize how little I knew about machine precision issues. The link to the R FAQ you provided made me wonder if the numerical error was arising at the level of R or already from the

Re: [Rcpp-devel] Numerical precision in rotations with Eigen

2024-01-10 Thread Dirk Eddelbuettel
On 10 January 2024 at 09:52, Rafael Ayala Hernandez wrote: | Hi, | | I have implemented a function to rotate a 3D vector a given angle around a given axis (basically wrapping the functionality provided by Eigen::AngleAxis) as an Rcpp function. | Below is an extract from the source file: | | #

[Rcpp-devel] Numerical precision in rotations with Eigen

2024-01-10 Thread Rafael Ayala Hernandez
Hi, I have implemented a function to rotate a 3D vector a given angle around a given axis (basically wrapping the functionality provided by Eigen::AngleAxis) as an Rcpp function. Below is an extract from the source file: #include #include #include #include using namespace Rcpp; using Rcpp::