Use xsimd::sincos when possible in xsimd::pow for complex number

serge-sans-paille · serge-sans-paille · commit 3b05fd63c808 · 2024-05-21T05:50:28.000Z
diff --git a/include/xsimd/arch/generic/xsimd_generic_math.hpp b/include/xsimd/arch/generic/xsimd_generic_math.hpp
@@ -2026,7 +2026,8 @@ namespace xsimd
             auto cond = (y == ze);
             r = select(cond, r, r * exp(-y * arga));
             theta = select(cond, theta, theta + y * log(absa));
-            return select(absa == ze, cplx_batch(ze), cplx_batch(r * cos(theta), r * sin(theta)));
+            auto sincosTheta = xsimd::sincos(theta);
+            return select(absa == ze, cplx_batch(ze), cplx_batch(r * sincosTheta.second, r * sincosTheta.first));
         }
 
         // reciprocal

Original file line number	Diff line number	Diff line change
`@@ -2026,7 +2026,8 @@ namespace xsimd`
`2026`	`2026`	`auto cond = (y == ze);`
`2027`	`2027`	`r = select(cond, r, r * exp(-y * arga));`
`2028`	`2028`	`theta = select(cond, theta, theta + y * log(absa));`
`2029`		`- return select(absa == ze, cplx_batch(ze), cplx_batch(r * cos(theta), r * sin(theta)));`
	`2029`	`+ auto sincosTheta = xsimd::sincos(theta);`
	`2030`	`+ return select(absa == ze, cplx_batch(ze), cplx_batch(r * sincosTheta.second, r * sincosTheta.first));`
`2030`	`2031`	`}`
`2031`	`2032`
`2032`	`2033`	`// reciprocal`