devtools/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2026-04-10 11:34:33 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fix ambiguous static_cast in JacobiSVD blocking threshold computation

Browse Source 
The L2 cache size threshold computation used numext::sqrt with a
static_cast<RealScalar>, which fails to compile when RealScalar is
AnnoyingScalar (a test-only type with multiple conversion constructors).
Since this is a pure cache-size computation unrelated to the matrix
scalar type, use std::sqrt(double) instead.

Co-Authored-By: Claude Opus 4.6 (1M context) <noreply@anthropic.com>
This commit is contained in:
Rasmus Munk Larsen
2026-02-25 19:03:01 -08:00
parent a31de4778d
commit e567151ce3
1 changed files with 1 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
Eigen/src/SVD/JacobiSVD.h
View File

@@ -774,7 +774,7 @@ JacobiSVD<MatrixType, Options>& JacobiSVD<MatrixType, Options>::compute_impl(con
const Index blockingThreshold = EIGEN_JACOBI_SVD_BLOCKING_THRESHOLD; const Index blockingThreshold = EIGEN_JACOBI_SVD_BLOCKING_THRESHOLD;
#else #else
const Index blockingThreshold = const Index blockingThreshold =
static_cast<Index>(numext::sqrt(static_cast<RealScalar>(l2CacheSize() / sizeof(float)))); static_cast<Index>(std::sqrt(static_cast<double>(l2CacheSize() / sizeof(float))));
#endif #endif
if (n >= blockingThreshold) { if (n >= blockingThreshold) {