qflux 0.0.2__tar.gz → 0.0.3__tar.gz

{qflux-0.0.2/src/qflux.egg-info → qflux-0.0.3}/PKG-INFO

@@ -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 → qflux-0.0.3}/pyproject.toml

@@ -11,7 +11,7 @@ where = ["src"]
 
 [project]
 name = "qflux"
-version = "0.0.2"
+version = "0.0.3"
 description = "qflux is a package for running quantum dynamics calculations on quantum devices."
 readme = "README.md"
 requires-python = ">=3.10"

{qflux-0.0.2 → qflux-0.0.3}/src/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-0.0.2 → qflux-0.0.3}/src/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 → qflux-0.0.3/src/qflux.egg-info}/PKG-INFO

@@ -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