diff --git a/qiskit/synthesis/two_qubit/two_qubit_decompose.py b/qiskit/synthesis/two_qubit/two_qubit_decompose.py
index 84989dc78d8d..1ac5b35fc9de 100644
--- a/qiskit/synthesis/two_qubit/two_qubit_decompose.py
+++ b/qiskit/synthesis/two_qubit/two_qubit_decompose.py
@@ -501,7 +501,7 @@ def _closest_partial_swap(a, b, c) -> float:
     am, bm, cm = a - m, b - m, c - m
     ab, bc, ca = a - b, b - c, c - a
 
-    return m + am * bm * cm * (6 + ab * ab + bc * bc * ca * ca) / 18
+    return m + am * bm * cm * (6 + ab * ab + bc * bc + ca * ca) / 18
 
 
 class TwoQubitWeylPartialSWAPEquiv(TwoQubitWeylDecomposition):