use log2l

SwayamInSync · SwayamInSync · commit d4d5b4e3029d · 2025-10-16T12:44:40.000+05:30
diff --git a/quaddtype/numpy_quaddtype/src/ops.hpp b/quaddtype/numpy_quaddtype/src/ops.hpp
@@ -832,13 +832,12 @@ ld_logaddexp2(const long double *x, const long double *y)
         return *x;
     }
     
-    // Numerically stable computation
     // log2(2^x + 2^y) = max(x, y) + log2(1 + 2^(-abs(x - y)))
     long double diff = *x - *y;
     long double abs_diff = fabsl(diff);
     long double max_val = (*x > *y) ? *x : *y;
-    // log2(1 + z) = log(1 + z) / log(2)
-    return max_val + log1pl(exp2l(-abs_diff)) / M_LN2;
+    // Use native log2l function for base-2 logarithm
+    return max_val + log2l(1.0L + exp2l(-abs_diff));
 }
 
 // comparison quad functions

Original file line number	Diff line number	Diff line change
`@@ -832,13 +832,12 @@ ld_logaddexp2(const long double x, const long double y)`
`832`	`832`	`return *x;`
`833`	`833`	`}`
`834`	`834`
`835`		`- // Numerically stable computation`
`836`	`835`	`// log2(2^x + 2^y) = max(x, y) + log2(1 + 2^(-abs(x - y)))`
`837`	`836`	`long double diff = x - y;`
`838`	`837`	`long double abs_diff = fabsl(diff);`
`839`	`838`	`long double max_val = (x > y) ? x : y;`
`840`		`- // log2(1 + z) = log(1 + z) / log(2)`
`841`		`- return max_val + log1pl(exp2l(-abs_diff)) / M_LN2;`
	`839`	`+ // Use native log2l function for base-2 logarithm`
	`840`	`+ return max_val + log2l(1.0L + exp2l(-abs_diff));`
`842`	`841`	`}`
`843`	`842`
`844`	`843`	`// comparison quad functions`