Switch to qiskit_algorithms

README.md

@@ -33,8 +33,8 @@ evaluate a fixed income asset with uncertain interest rates.
 
 ```python
 import numpy as np
-from qiskit.algorithms import AmplitudeEstimation
 from qiskit.primitives import Sampler
+from qiskit_algorithms import AmplitudeEstimation
 from qiskit_finance.circuit.library import NormalDistribution
 from qiskit_finance.applications import FixedIncomePricing
 

docs/conf.py

@@ -191,6 +191,7 @@
     "python": ("https://docs.python.org/3", None),
     "numpy": ("https://numpy.org/doc/stable", None),
     "scipy": ("https://docs.scipy.org/doc/scipy", None),
+    "qiskit_algorithms": ("https://qiskit.org/ecosystem/algorithms", None),
     "qiskit_optimization": ("https://qiskit.org/ecosystem/optimization", None),
     "qiskit": ("https://qiskit.org/documentation", None),
 }

qiskit_finance/applications/estimation/estimation_application.py

@@ -13,7 +13,7 @@
 """An abstract class for estimation application classes."""
 from abc import ABC, abstractmethod
 
-from qiskit.algorithms.amplitude_estimators import (
+from qiskit_algorithms import (
     EstimationProblem,
     AmplitudeEstimatorResult,
 )
@@ -27,14 +27,14 @@ class EstimationApplication(ABC):
     @abstractmethod
     def to_estimation_problem(self) -> EstimationProblem:
         """Convert a problem instance into a
-        :class:`qiskit.algorithms.EstimationProblem`
+        :class:`qiskit_algorithms.EstimationProblem`
         """
         pass
 
     @abstractmethod
     def interpret(self, result: AmplitudeEstimatorResult) -> float:
         """Convert the calculation result of the problem
-        (:class:`qiskit.algorithms.AmplitudeEstimatorResult`)
+        (:class:`qiskit_algorithms.AmplitudeEstimatorResult`)
         to the answer of the problem.
 
         Args:

qiskit_finance/applications/estimation/european_call_delta.py

@@ -14,7 +14,7 @@
 from typing import Tuple
 
 from qiskit.circuit import QuantumCircuit
-from qiskit.algorithms.amplitude_estimators import (
+from qiskit_algorithms import (
     EstimationProblem,
     AmplitudeEstimatorResult,
 )
@@ -60,10 +60,10 @@ def __init__(
 
     def to_estimation_problem(self) -> EstimationProblem:
         """Convert a problem instance into a
-        :class:`qiskit.algorithms.EstimationProblem`
+        :class:`qiskit_algorithms.EstimationProblem`
 
         Returns:
-            The :class:`qiskit.algorithms.EstimationProblem` created
+            The :class:`qiskit_algorithms.EstimationProblem` created
             from the European call delta problem instance.
         """
         problem = EstimationProblem(
@@ -75,7 +75,7 @@ def to_estimation_problem(self) -> EstimationProblem:
 
     def interpret(self, result: AmplitudeEstimatorResult) -> float:
         """Convert the calculation result of the problem
-        (:class:`qiskit.algorithms.AmplitudeEstimatorResult`)
+        (:class:`qiskit_algorithms.AmplitudeEstimatorResult`)
         to the answer of the problem.
 
         Args:

qiskit_finance/applications/estimation/european_call_pricing.py

@@ -14,7 +14,7 @@
 from typing import Tuple
 
 from qiskit.circuit import QuantumCircuit
-from qiskit.algorithms.amplitude_estimators import (
+from qiskit_algorithms import (
     EstimationProblem,
     AmplitudeEstimatorResult,
 )
@@ -65,10 +65,10 @@ def __init__(
 
     def to_estimation_problem(self) -> EstimationProblem:
         """Convert a problem instance into a
-        :class:`qiskit.algorithms.EstimationProblem`
+        :class:`qiskit_algorithms.EstimationProblem`
 
         Returns:
-            The :class:`qiskit.algorithms.EstimationProblem` created
+            The :class:`qiskit_algorithms.EstimationProblem` created
             from the European call pricing problem instance.
         """
         problem = EstimationProblem(
@@ -80,7 +80,7 @@ def to_estimation_problem(self) -> EstimationProblem:
 
     def interpret(self, result: AmplitudeEstimatorResult) -> float:
         """Convert the calculation result of the problem
-        (:class:`qiskit.algorithms.AmplitudeEstimatorResult`)
+        (:class:`qiskit_algorithms.AmplitudeEstimatorResult`)
         to the answer of the problem.
 
         Args:

qiskit_finance/applications/estimation/fixed_income_pricing.py

@@ -16,7 +16,7 @@
 import numpy as np
 
 from qiskit.circuit import QuantumCircuit
-from qiskit.algorithms.amplitude_estimators import (
+from qiskit_algorithms import (
     EstimationProblem,
     AmplitudeEstimatorResult,
 )
@@ -83,10 +83,10 @@ def __init__(
 
     def to_estimation_problem(self) -> EstimationProblem:
         """Convert a problem instance into a
-        :class:`qiskit.algorithms.EstimationProblem`
+        :class:`qiskit_algorithms.EstimationProblem`
 
         Returns:
-            The :class:`qiskit.algorithms.EstimationProblem` created
+            The :class:`qiskit_algorithms.EstimationProblem` created
             from the Fixed problem instance.
         """
         problem = EstimationProblem(
@@ -98,7 +98,7 @@ def to_estimation_problem(self) -> EstimationProblem:
 
     def interpret(self, result: AmplitudeEstimatorResult) -> float:
         """Convert the calculation result of the problem
-        (:class:`qiskit.algorithms.AmplitudeEstimatorResult`)
+        (:class:`qiskit_algorithms.AmplitudeEstimatorResult`)
         to the answer of the problem.
 
         Args:

qiskit_finance/data_providers/_base_data_provider.py

@@ -20,7 +20,7 @@
 import numpy as np
 import fastdtw
 
-from qiskit.utils import algorithm_globals
+from qiskit_algorithms.utils import algorithm_globals
 from ..exceptions import QiskitFinanceError
 
 logger = logging.getLogger(__name__)

requirements.txt

@@ -1,4 +1,5 @@
 qiskit-terra>=0.20.0
+qiskit-algorithms>=0.2.0
 qiskit-optimization>=0.2.0
 scipy>=1.4
 numpy>=1.17

test/applications/test_european_call_delta.py

@@ -15,7 +15,6 @@
 import unittest
 from test import QiskitFinanceTestCase
 
-from qiskit.utils import algorithm_globals
 from qiskit.circuit.library import IntegerComparator
 from qiskit.quantum_info import Operator
 from qiskit_finance.applications.estimation import EuropeanCallDelta
@@ -25,11 +24,6 @@
 class TestEuropeanCallDelta(QiskitFinanceTestCase):
     """Tests the EuropeanCallDelta application"""
 
-    def setUp(self):
-        super().setUp()
-        self.seed = 457
-        algorithm_globals.random_seed = self.seed
-
     def test_to_estimation_problem(self):
         """Test the expected circuit."""
         num_qubits = 3

test/applications/test_european_call_pricing.py

@@ -17,19 +17,13 @@
 
 from qiskit.circuit.library import LinearAmplitudeFunction
 from qiskit.quantum_info import Operator
-from qiskit.utils import algorithm_globals
 from qiskit_finance.applications.estimation import EuropeanCallPricing
 from qiskit_finance.circuit.library.probability_distributions import UniformDistribution
 
 
 class TestEuropeanCallPricing(QiskitFinanceTestCase):
     """Tests the EuropeanCallPricing application"""
 
-    def setUp(self):
-        super().setUp()
-        self.seed = 457
-        algorithm_globals.random_seed = self.seed
-
     def test_to_estimation_problem(self):
         """Test the expected circuit."""
         num_qubits = 3

test/applications/test_fixed_income_pricing.py

@@ -19,7 +19,6 @@
 
 from qiskit import QuantumCircuit
 from qiskit.quantum_info import Operator
-from qiskit.utils import algorithm_globals
 from qiskit_finance.applications.estimation import FixedIncomePricing
 from qiskit_finance.circuit.library.probability_distributions import UniformDistribution
 from qiskit_finance.circuit.library.payoff_functions import FixedIncomePricingObjective
@@ -28,11 +27,6 @@
 class TestFixedIncomePricing(QiskitFinanceTestCase):
     """Tests the FixedIncomePricing application"""
 
-    def setUp(self):
-        super().setUp()
-        self.seed = 457
-        algorithm_globals.random_seed = self.seed
-
     def test_to_estimation_problem(self):
         """Test the expected circuit."""
         num_qubits = [2, 2]

test/applications/test_portfolio_diversification.py

@@ -17,7 +17,6 @@
 from test import QiskitFinanceTestCase
 
 import numpy as np
-from qiskit.utils import algorithm_globals
 from qiskit_optimization.problems import QuadraticProgram, VarType
 from qiskit_finance.applications.optimization import PortfolioDiversification
 
@@ -30,7 +29,6 @@ class TestPortfolioDiversification(QiskitFinanceTestCase):
     def setUp(self):
         """Set up for the tests"""
         super().setUp()
-        algorithm_globals.random_seed = 100
         self.n = 2
         self.q = 1
         self.similarity_matrix = np.ones((self.n, self.n))

test/applications/test_portfolio_optimization.py

@@ -18,7 +18,6 @@
 from ddt import ddt, data, unpack
 import numpy as np
 
-from qiskit.utils import algorithm_globals
 from qiskit_optimization.problems import QuadraticProgram
 from qiskit_finance.applications.optimization import PortfolioOptimization
 from qiskit_finance.exceptions import QiskitFinanceError
@@ -33,7 +32,6 @@ class TestPortfolioDiversification(QiskitFinanceTestCase):
     def setUp(self):
         """Set up for the tests"""
         super().setUp()
-        algorithm_globals.random_seed = 100
         self.num_assets = 4
         self.expected_returns = [0.01528439, -0.00078095, 0.00051792, 0.00087001]
         self.covariances = [

test/circuit/test_european_call_delta_objective.py

@@ -18,10 +18,9 @@
 import numpy as np
 from qiskit.circuit.library import IntegerComparator
 from qiskit.quantum_info import Operator
-from qiskit.utils import QuantumInstance, optionals
-from qiskit.algorithms import IterativeAmplitudeEstimation, EstimationProblem
+from qiskit.primitives import Sampler
+from qiskit_algorithms import IterativeAmplitudeEstimation, EstimationProblem
 from qiskit_finance.circuit.library import LogNormalDistribution
-
 from qiskit_finance.circuit.library.payoff_functions import EuropeanCallDeltaObjective
 
 
@@ -46,7 +45,6 @@ def test_circuit(self):
 
         self.assertTrue(Operator(ecd).equiv(comparator))
 
-    @unittest.skipUnless(optionals.HAS_AER, "qiskit-aer is required to run this test")
     def test_application(self):
         """Test an end-to-end application."""
 
@@ -93,17 +91,10 @@ def test_application(self):
             post_processing=european_call_delta.post_processing,
         )
 
-        # run amplitude estimation
-        from qiskit_aer import Aer
-
-        q_i = QuantumInstance(
-            Aer.get_backend("aer_simulator"),
-            seed_simulator=125,
-            seed_transpiler=80,
-        )
-        iae = IterativeAmplitudeEstimation(epsilon_target=0.01, alpha=0.05, quantum_instance=q_i)
+        sampler = Sampler(options={"shots": 1024, "seed": 12})
+        iae = IterativeAmplitudeEstimation(epsilon_target=0.01, alpha=0.05, sampler=sampler)
         result = iae.estimate(problem)
-        self.assertAlmostEqual(result.estimation_processed, 0.8088790606143996)
+        self.assertAlmostEqual(result.estimation_processed, 0.8088153392162598)
 
 
 if __name__ == "__main__":

test/circuit/test_european_call_pricing_objective.py

@@ -16,21 +16,16 @@
 from test import QiskitFinanceTestCase
 
 import numpy as np
-from qiskit.utils import algorithm_globals, QuantumInstance, optionals
-from qiskit.algorithms import IterativeAmplitudeEstimation, EstimationProblem
 from qiskit.circuit.library import LinearAmplitudeFunction, TwoLocal
 from qiskit.quantum_info import Operator
+from qiskit.primitives import Sampler
+from qiskit_algorithms import IterativeAmplitudeEstimation, EstimationProblem
 from qiskit_finance.circuit.library import EuropeanCallPricingObjective, NormalDistribution
 
 
 class TestEuropeanCallExpectedValue(QiskitFinanceTestCase):
     """Tests EuropeanCallPricingObjective."""
 
-    def setUp(self):
-        super().setUp()
-        self.seed = 457
-        algorithm_globals.random_seed = self.seed
-
     def test_ecev_circuit(self):
         """Test the expected circuit.
 
@@ -59,7 +54,6 @@ def test_ecev_circuit(self):
 
         self.assertTrue(Operator(ecev).equiv(linear_function))
 
-    @unittest.skipUnless(optionals.HAS_AER, "qiskit-aer is required to run this test")
     def test_application(self):
         """Test an end-to-end application."""
 
@@ -101,16 +95,10 @@ def test_application(self):
             post_processing=european_call.post_processing,
         )
 
-        from qiskit_aer import Aer
-
-        q_i = QuantumInstance(
-            Aer.get_backend("aer_simulator"),
-            seed_simulator=125,
-            seed_transpiler=80,
-        )
-        iae = IterativeAmplitudeEstimation(epsilon_target=0.01, alpha=0.05, quantum_instance=q_i)
+        sampler = Sampler(options={"shots": 1024, "seed": 12})
+        iae = IterativeAmplitudeEstimation(epsilon_target=0.01, alpha=0.05, sampler=sampler)
         result = iae.estimate(problem)
-        self.assertAlmostEqual(result.estimation_processed, 1.0364776997977694)
+        self.assertAlmostEqual(result.estimation_processed, 1.0341976859652098)
 
 
 if __name__ == "__main__":

test/circuit/test_fixed_income_pricing_objective.py

@@ -18,9 +18,9 @@
 import numpy as np
 
 from qiskit import QuantumCircuit
-from qiskit.utils import QuantumInstance, optionals
-from qiskit.algorithms import IterativeAmplitudeEstimation, EstimationProblem
 from qiskit.quantum_info import Operator
+from qiskit.primitives import Sampler
+from qiskit_algorithms import IterativeAmplitudeEstimation, EstimationProblem
 from qiskit_finance.circuit.library import NormalDistribution
 from qiskit_finance.circuit.library.payoff_functions import FixedIncomePricingObjective
 
@@ -50,7 +50,6 @@ def test_circuit(self):
 
         self.assertTrue(Operator(circuit).equiv(expected))
 
-    @unittest.skipUnless(optionals.HAS_AER, "qiskit-aer is required to run this test")
     def test_application(self):
         """Test an end-to-end application."""
 
@@ -87,19 +86,12 @@ def test_application(self):
             post_processing=fixed_income.post_processing,
         )
 
-        # run simulation
-        from qiskit_aer import Aer
-
-        q_i = QuantumInstance(
-            Aer.get_backend("aer_simulator"),
-            seed_simulator=2,
-            seed_transpiler=2,
-        )
-        iae = IterativeAmplitudeEstimation(epsilon_target=0.01, alpha=0.05, quantum_instance=q_i)
+        sampler = Sampler(options={"shots": 1024, "seed": 12})
+        iae = IterativeAmplitudeEstimation(epsilon_target=0.01, alpha=0.05, sampler=sampler)
         result = iae.estimate(problem)
 
         # compare to precomputed solution
-        self.assertAlmostEqual(result.estimation_processed, 2.3389012822103044)
+        self.assertAlmostEqual(result.estimation_processed, 2.329154511815111)
 
 
 if __name__ == "__main__":