Add branchless w() variants sadly we had worse performance

2026-05-20 22:23:07 +02:00
parent 2cd9031913
commit c3f6fa3b27
4 changed files with 951 additions and 0 deletions
@@ -96,4 +96,7 @@ const double EPS = 1e-9;
 // Squared numerical tolerance
 const double EPS_SQ = EPS * EPS;

+
+const double cos_phi_sq = cos_phi * cos_phi;
+
 const size_t CACHE_BLOCK_SIZE = std::hardware_destructive_interference_size;