From 1821b980df34f92e32dabab2e52d9f8dd1d1461d Mon Sep 17 00:00:00 2001
From: Jake Lishman <jake.lishman@ibm.com>
Date: Fri, 22 Sep 2023 20:31:31 +0100
Subject: [PATCH] Replace use of `numpy.lib.scimath` with preferred
 `numpy.emath` (#10892)

`numpy.emath` is the preferred alias to `numpy.lib.scimath`, and direct
access to the latter is scheduled for removal in Numpy 2.

(cherry picked from commit af242271f8efc2af1d6d2ac1b6eaf82cbdcf7fe5)
---
 qiskit/quantum_info/synthesis/qsd.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/qiskit/quantum_info/synthesis/qsd.py b/qiskit/quantum_info/synthesis/qsd.py
index b2f00a56045d..90278dad62ed 100644
--- a/qiskit/quantum_info/synthesis/qsd.py
+++ b/qiskit/quantum_info/synthesis/qsd.py
@@ -175,7 +175,7 @@ def _demultiplex(um0, um1, opt_a1=False, opt_a2=False, *, _depth=0):
     else:
         evals, vmat = scipy.linalg.schur(um0um1, output="complex")
         eigvals = evals.diagonal()
-    dvals = np.lib.scimath.sqrt(eigvals)
+    dvals = np.emath.sqrt(eigvals)
     dmat = np.diag(dvals)
     wmat = dmat @ vmat.T.conjugate() @ um1