qflux 0.0.2__py3-none-any.whl → 0.0.3__py3-none-any.whl

qflux/closed_systems/VarQTE.py

@@ -3,15 +3,14 @@ import numpy.typing as npt
 from typing import List, Optional, Tuple
 
 from qiskit import QuantumCircuit
-from qiskit.primitives import Estimator
+from qiskit_aer.primitives import EstimatorV2 as Estimator
 from qiskit.quantum_info import SparsePauliOp
-from qiskit.providers.aer.noise import NoiseModel
-from qiskit.providers.fake_provider import FakeSherbrooke
-
+from qiskit_aer.noise import NoiseModel
+from qiskit_ibm_runtime.fake_provider import FakeSherbrooke
 
 # To change the ansatz, apply_param and measure_der must both be modified.
 def apply_param(
-    params: npt.NDArray[np.float_], i: int, qc: QuantumCircuit, N: int
+    params: npt.NDArray[np.float64], i: int, qc: QuantumCircuit, N: int
 ) -> None:
     """Apply parameter i to the quantum circuit currently constructing the ansatz.
         The ansatz must be built in a peicewise manner to allow for hadamard tests
@@ -24,7 +23,7 @@ def apply_param(
         qc (QuantumCircuit): The qiskit ansatz quantum circuit currently being constructed.
         N (int): Number of qubits
     """
-    qc.rx(params[i], i % N))
+    qc.rx(params[i], i % N)
     if i % N == N - 1 and i != len(params) - 1:
         for i in range(N - 1):
             qc.cz(i, i + 1)
@@ -43,7 +42,7 @@ def measure_der(i: int, qc: QuantumCircuit, N: int) -> None:
     qc.cx(N, i % N)
 
 
-def A_Circuit(params: npt.NDArray[np.float_], i: int, j: int, N: int) -> QuantumCircuit:
+def A_Circuit(params: npt.NDArray[np.float64], i: int, j: int, N: int) -> QuantumCircuit:
     """Constructs the qiskit quantum circuits used to measure each element of the A_ij matrix.
 
     Args:
@@ -73,11 +72,11 @@ def A_Circuit(params: npt.NDArray[np.float_], i: int, j: int, N: int) -> Quantum
 
 def Measure_A(
     init_circ: QuantumCircuit,
-    params: npt.NDArray[np.float_],
+    params: npt.NDArray[np.float64],
     N: int,
     shots: int = 2**10,
     noisy: bool = False,
-) -> npt.NDArray[np.float_]:
+) -> npt.NDArray[np.float64]:
     """Create the A_ij matrix through measuring quantum circuits corresponding to each element.
 
     Args:
@@ -104,23 +103,19 @@ def Measure_A(
                 coupling_map = device_backend.coupling_map
                 noise_model = NoiseModel.from_backend(device_backend)
                 basis_gates = noise_model.basis_gates
-                estimator = Estimator(
-                    options={
-                        "shots": shots,
-                        "noise_model": noise_model,
-                        "coupling_map": coupling_map,
-                        "basis_gates": basis_gates,
-                    }
+                estimator = Estimator(options={
+                                    "backend_options":{"noise_model": noise_model},
+                                    "run_options":{"shots": shots}}
                 )
             else:
-                estimator = Estimator(options={"shots": shots})
-            result = estimator.run(qc, observable).result()
-            A[i][i + j] = result.values[0]
+                estimator = Estimator(options={"run_options":{"shots": shots}})
+            result = estimator.run([(qc, observable)]).result()
+            A[i][i + j] = result[0].data.evs
     return np.array(A)
 
 
 def C_Circuit(
-    params: npt.NDArray[np.float_],
+    params: npt.NDArray[np.float64],
     i: int,
     pauli_string: str,
     N: int,
@@ -162,13 +157,13 @@ def C_Circuit(
 
 def Measure_C(
     init_circ: QuantumCircuit,
-    params: npt.NDArray[np.float_],
+    params: npt.NDArray[np.float64],
     H: SparsePauliOp,
     N: int,
     shots: int = 2**10,
     evolution_type: str = "real",
     noisy: bool = False,
-) -> npt.NDArray[np.float_]:
+) -> npt.NDArray[np.float64]:
     """Create the C_i vector through measuring quantum circuits corresponding to each element.
 
     Args:
@@ -199,19 +194,15 @@ def Measure_C(
                 coupling_map = device_backend.coupling_map
                 noise_model = NoiseModel.from_backend(device_backend)
                 basis_gates = noise_model.basis_gates
-                estimator = Estimator(
-                    options={
-                        "shots": shots,
-                        "noise_model": noise_model,
-                        "coupling_map": coupling_map,
-                        "basis_gates": basis_gates,
-                    }
-                )
+                estimator = Estimator(options={
+                                        "backend_options":{"noise_model": noise_model},
+                                        "run_options":{"shots": shots}}
+                    )
             else:
-                estimator = Estimator(options={"shots": shots})
-            result = estimator.run(qc, observable).result()
+                estimator = Estimator(options={"run_options":{"shots": shots}})
+            result = estimator.run([(qc, observable)]).result()
 
-            C[i] -= 1 / 2 * H.coeffs[pauli_string].real * result.values[0]
+            C[i] -= 1 / 2 * H.coeffs[pauli_string].real * result[0].data.evs
     return np.array(C)
 
 
@@ -233,7 +224,7 @@ def pauli_measure(qc: QuantumCircuit, pauli_string: str) -> None:
 
 
 def Construct_Ansatz(
-    init_circ: QuantumCircuit, params: npt.NDArray[np.float_], N: int
+    init_circ: QuantumCircuit, params: npt.NDArray[np.float64], N: int
 ) -> QuantumCircuit:
     """Construct the full ansatz for use in measuring observables.
 
@@ -258,7 +249,7 @@ def Construct_Ansatz(
 
 def ansatz_energy(
     init_circ: QuantumCircuit,
-    params: npt.NDArray[np.float_],
+    params: npt.NDArray[np.float64],
     H: SparsePauliOp,
     shots: int = 2**14,
     noisy: bool = False,
@@ -282,19 +273,15 @@ def ansatz_energy(
         coupling_map = device_backend.coupling_map
         noise_model = NoiseModel.from_backend(device_backend)
         basis_gates = noise_model.basis_gates
-        estimator = Estimator(
-            options={
-                "shots": shots,
-                "noise_model": noise_model,
-                "coupling_map": coupling_map,
-                "basis_gates": basis_gates,
-            }
-        )
+        estimator = Estimator(options={
+                                "backend_options":{"noise_model": noise_model},
+                                "run_options":{"shots": shots}}
+                                )
     else:
-        estimator = Estimator(options={"shots": shots})
+        estimator = Estimator(options={"run_options":{"shots": shots}})
     qc = Construct_Ansatz(init_circ, params, N)
-    result = estimator.run(qc, H).result()
-    return result.values[0], result.metadata[0]["variance"]
+    result = estimator.run([(qc, H)]).result()
+    return result[0].data.evs, result[0].data.stds
 
 
 def VarQRTE(
@@ -305,7 +292,7 @@ def VarQRTE(
     init_circ: Optional[QuantumCircuit] = None,
     shots: int = 2**10,
     noisy: bool = False,
-) -> List[npt.NDArray[np.float_]]:
+) -> List[npt.NDArray[np.float64]]:
     """The Variational Quantum Real Time Evolution (VarQRTE) algorithm.  This uses quantum circuits to measure
         the elements of two objects, the A_ij matrix and the C_i vector.
 
@@ -368,7 +355,7 @@ def VarQITE(
     init_circ: Optional[QuantumCircuit] = None,
     shots: int = 2**10,
     noisy: bool = False,
-) -> List[npt.NDArray[np.float_]]:
+) -> List[npt.NDArray[np.float64]]:
     """The Variational Quantum Imaginary Time Evolution (VarQITE) algorithm.  This uses quantum circuits to measure
         the elements of two objects, the A_ij matrix and the C_i vector.
 

qflux/closed_systems/qubit_methods.py

@@ -220,7 +220,7 @@ class QubitDynamicsCS(DynamicsCS):
 
 
     def propagate_qmatvec(self, backend=None, n_shots: int = 1024, hamiltonian_matrix=None, initial_state=None):
-        """
+        r"""
             Function to propagate dynamics object with the qubit matvec method.
 
             Args:

{qflux-0.0.2.dist-info → qflux-0.0.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: qflux
-Version: 0.0.2
+Version: 0.0.3
 Summary: qflux is a package for running quantum dynamics calculations on quantum devices.
 Author-email: Brandon Allen <brandon.allen@yale.edu>, Delmar Cabral <delmar.azevedocabral@yale.edu>, Alexander Soudackov <alexander.soudackov@yale.edu>, Anton Morgunov <anton@ischemist.com>
 License: MIT

{qflux-0.0.2.dist-info → qflux-0.0.3.dist-info}/RECORD

@@ -5,12 +5,12 @@ qflux/GQME/params.py,sha256=u7mHRdbbf88ccTIvw9Hb4AunlfNBrzara0OZTL80vwk,1822
 qflux/GQME/readwrite.py,sha256=JUUEsDL8Ep705P8iwlQ4YjAm7MswQ-eNpH1wa8_P-_I,3583
 qflux/GQME/tdvp.py,sha256=32Z_cLJornzM4D63nmqekhwtc-OosrQcZhIeEZAqq9s,7136
 qflux/GQME/tt_tfd.py,sha256=-RPH9aiCinLzG0VdZo43UBVKWcuoplkFmi77oOtv5zU,15204
-qflux/closed_systems/VarQTE.py,sha256=8Znw6R0jBsHnV_e5eGGMumcSG9HlUfRltuqevlBka9A,17956
+qflux/closed_systems/VarQTE.py,sha256=rNDFnRmF_jDv-wSSaBnFp0ZwrgjOJ3NdqyoT4WUlX4E,17840
 qflux/closed_systems/__init__.py,sha256=CBnnqCBp0MD7jCGq1P0nzwCyR3okp4jkC8tvqWUTBCY,472
 qflux/closed_systems/classical_methods.py,sha256=Etqjd9OkqCgPzhC2-FgQWlIiTu66V2_7LRftPDwBrb8,15882
 qflux/closed_systems/custom_execute.py,sha256=Mk24DEvEbQlnu52MqhnSOGBpplHoAMz5__AcBKo1E_4,656
 qflux/closed_systems/hamiltonians.py,sha256=ArdljgB7ASUq5QIKdN8CRdhf7b_2nFyb8OPacOGGvEs,3156
-qflux/closed_systems/qubit_methods.py,sha256=NPC0pLM8Jok4bpkTBcDvDNO3YqsCpw9S6MVMxRHwcuU,10882
+qflux/closed_systems/qubit_methods.py,sha256=2TFeR7lxAqAu3CYUPQqUJtN9bMPpKTukHgGohiJJASk,10883
 qflux/closed_systems/spin_dynamics_oo.py,sha256=cn4aRM3wYkP6CBbL-LhlUMR4dqQSQK7bpSV-twSzFbU,13925
 qflux/closed_systems/spin_propagators.py,sha256=pMg4ttThkXw2kzBlY2yDpdomTn-7alRtpKUAxD-vSl8,11496
 qflux/closed_systems/utils.py,sha256=eiB5V0xvi43-jC35nKl2kMY_vr9r_-2njNn5qPLh494,6975
@@ -32,8 +32,8 @@ qflux/variational_methods/qmad/ansatz.py,sha256=2ttVOrueLYycyXaBy1uxPkaTNEtkctL9
 qflux/variational_methods/qmad/ansatzVect.py,sha256=OJdRDhKBYou93fwetri5KXFoAj9N4cdjj799Nqu-qmE,2934
 qflux/variational_methods/qmad/effh.py,sha256=l74femDGDDVCiHerfDDj2k30phyrYIlI_ssn4YU426g,3722
 qflux/variational_methods/qmad/solver.py,sha256=1B9DpcuT7bpKcIWt1kHkxT3-gwQECdEzKFcBtJlFYDg,12751
-qflux-0.0.2.dist-info/licenses/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
-qflux-0.0.2.dist-info/METADATA,sha256=1rRhjqS2949C8qCYtHxjAnemMNy7vEr-N6p_ZQ-B3qg,6871
-qflux-0.0.2.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-qflux-0.0.2.dist-info/top_level.txt,sha256=GZgw1Z1DZDPg78NCF1dXCHw_f4usRh5HP5yLtgECIpo,6
-qflux-0.0.2.dist-info/RECORD,,
+qflux-0.0.3.dist-info/licenses/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
+qflux-0.0.3.dist-info/METADATA,sha256=sZdUGSirgVQ6W6yT5E-jZp4BTVCltcDlme9AJ6ILGx4,6871
+qflux-0.0.3.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+qflux-0.0.3.dist-info/top_level.txt,sha256=GZgw1Z1DZDPg78NCF1dXCHw_f4usRh5HP5yLtgECIpo,6
+qflux-0.0.3.dist-info/RECORD,,