diff --git a/src/scalar/complex.rs b/src/scalar/complex.rs
index 4c58ea3..df2fb22 100644
--- a/src/scalar/complex.rs
+++ b/src/scalar/complex.rs
@@ -1455,5 +1455,6 @@ impl<N: RealField + PartialOrd> ComplexField for num_complex::Complex<N> {
 
 #[inline]
 fn complex_from_polar<N: RealField>(r: N, theta: N) -> num_complex::Complex<N> {
-    num_complex::Complex::new(r.clone() * theta.clone().cos(), r * theta.sin())
+    let (sin, cos) = theta.sin_cos();
+    num_complex::Complex::new(r.clone() * cos, r * sin)
 }
diff --git a/src/simd/auto_simd_impl.rs b/src/simd/auto_simd_impl.rs
index 9526d08..108d58f 100644
--- a/src/simd/auto_simd_impl.rs
+++ b/src/simd/auto_simd_impl.rs
@@ -1479,7 +1479,8 @@ macro_rules! impl_float_simd (
 
 #[inline]
 fn simd_complex_from_polar<N: SimdRealField>(r: N, theta: N) -> num_complex::Complex<N> {
-    num_complex::Complex::new(r.clone() * theta.clone().simd_cos(), r * theta.simd_sin())
+    let (sin, cos) = theta.simd_sin_cos();
+    num_complex::Complex::new(r.clone() * cos, r * sin)
 }
 
 impl_float_simd!(
diff --git a/src/simd/portable_simd_impl.rs b/src/simd/portable_simd_impl.rs
index 68ad1ac..c64e508 100644
--- a/src/simd/portable_simd_impl.rs
+++ b/src/simd/portable_simd_impl.rs
@@ -1501,7 +1501,8 @@ macro_rules! impl_float_simd (
 
 #[inline]
 fn simd_complex_from_polar<N: SimdRealField>(r: N, theta: N) -> num_complex::Complex<N> {
-    num_complex::Complex::new(r.clone() * theta.clone().simd_cos(), r * theta.simd_sin())
+    let (sin, cos) = theta.simd_sin_cos();
+    num_complex::Complex::new(r.clone() * cos, r * sin)
 }
 
 impl_float_simd!(
diff --git a/src/simd/wide_simd_impl.rs b/src/simd/wide_simd_impl.rs
index 70f042d..5c430a3 100644
--- a/src/simd/wide_simd_impl.rs
+++ b/src/simd/wide_simd_impl.rs
@@ -1584,5 +1584,6 @@ impl_wide_f32!(f32, f32x8, WideF32x8, WideBoolF32x8, 8; 1, 2, 3, 4, 5, 6, 7);
 
 #[inline]
 fn simd_complex_from_polar<N: SimdRealField>(r: N, theta: N) -> num_complex::Complex<N> {
-    num_complex::Complex::new(r.clone() * theta.clone().simd_cos(), r * theta.simd_sin())
+    let (sin, cos) = theta.simd_sin_cos();
+    num_complex::Complex::new(r.clone() * cos, r * sin)
 }