Update

vutuanhai237 · Mar 1, 2024 · 98bba9f · 98bba9f
1 parent c91d873
commit 98bba9f
Show file tree

Hide file tree

Showing 14 changed files with 13,220 additions and 645 deletions.
diff --git a/codes/Wchain_zxz.eps b/codes/Wchain_zxz.eps
diff --git a/codes/cnot.eps b/codes/cnot.eps
diff --git a/codes/compare.ipynb b/codes/compare.ipynb
diff --git a/codes/custom.eps b/codes/custom.eps
diff --git a/codes/evolution_circuit_XY.eps b/codes/evolution_circuit_XY.eps
diff --git a/codes/magnetization.eps b/codes/magnetization.eps
diff --git a/codes/qsp.ipynb b/codes/qsp.ipynb
diff --git a/codes/zxz_WchainCNOT_2_2.eps b/codes/zxz_WchainCNOT_2_2.eps
diff --git a/codes/zxz_walltoallcnot.eps b/codes/zxz_walltoallcnot.eps
diff --git a/codes/zxz_walltoallcnot4.eps b/codes/zxz_walltoallcnot4.eps
diff --git a/codes/zxz_walternatingcnot.eps b/codes/zxz_walternatingcnot.eps
diff --git a/codes/zxz_wchainCNOT_eq2.eps.eps b/codes/zxz_wchainCNOT_eq2.eps.eps
diff --git a/codes/zxz_wchaincnot.eps b/codes/zxz_wchaincnot.eps
diff --git a/qsp.py b/qsp.py
@@ -0,0 +1,40 @@
+from qsee.compilation.qsp import QuantumStatePreparation
+from qsee.core import ansatz, state
+from qsee.backend import constant, utilities
+import matplotlib.pyplot as plt
+import numpy as np, qiskit
+from scipy.linalg import expm
+from qiskit.quantum_info import Pauli, Statevector
+from qsee.compilation.qsp import QuantumStatePreparation
+
+# Define Pauli matrices
+XX = Pauli("XX").to_matrix()
+YY = Pauli("YY").to_matrix()
+ZZ = Pauli("ZZ").to_matrix()
+ZI = Pauli("ZI").to_matrix()
+IZ = Pauli("IZ").to_matrix()
+# Set coefficients
+J = 1.0
+D = 0.0
+h = 0.0
+# Define H
+H = J * (XX + YY) + D * ZZ + h * (ZI - IZ)
+# define V(t), calculate U and local magnetization
+time = np.linspace(0, 10, 10)
+magnetization = []
+for t in time:
+    V = expm(-1j * t * H)
+    compiler = QuantumStatePreparation(u=ansatz.zxz_WchainCNOT(2), target_state=V).fit(num_steps = 30)
+    statevector = Statevector.from_instruction(
+        compiler.u.assign_parameters(compiler.thetas)
+    )
+    # print(V)
+    # print(V)
+    print(compiler.fidelity)
+    # print(compiler.thetas)
+    psi_t = np.array(statevector.data).reshape(-1, 1)
+    psi_t_dag = np.conj(psi_t.T)
+    mag = (psi_t_dag @ ZI @ psi_t)[0][0] - (psi_t_dag @ IZ @ psi_t)[0][0]
+    magnetization.append(0.5 * mag)
+
+plt.plot(time, magnetization)