I have some very big TLS point clouds (up to 400 million points) and want to assign the eigenvalues to each point in order to later calculate geometric features like planarity or sphericity. I use R and the package `lidR`. I realized that the function `point_metrics()` in combination with `fast_eigen_values()` is too slow. The function `segment_shape()` is way faster and also computes eigenvalues in the background. However, looking at the C++ Code behind `lidR` - although I can't write C++ - it seems to be that it is written in a way that it only able to return booleans, not doubles. Does someone know how to solve this problem? Was there already a version somewhere implemented, which enables me to use the spatial indexing to return eigenvalues without needing to learn C++ just for this? My next resort would be to try using `voxel_metrics()`. However, I would prefer to have one set of eigenvalues for each point.