.. _oceanstack: ==================== Ocean Software Stack ==================== The Ocean software stack provides a chain of tools that implements the steps needed to solve your problem on a CPU/GPU or a D-Wave system. As described in the :ref:`solving_problems` section, these steps include formulating the problem in a way the quantum computer understands (as a :term:`binary quadratic model`) and solving the formulated problem by submitting it to a D-Wave system or classical :term:`sampler` (the component used to minimize a BQM and therefore solve the original problem). It's helpful to visualize the tool chain as layers of abstraction, each of which handles one part of the solution procedure. Abstraction Layers ================== .. _fig_stack: .. figure:: ../_images/ocean_stack.png :name: stack :scale: 100 % :alt: Overview of the software stack. :height: 400 pt :width: 400 pt Ocean Software Stack The :ref:`fig_stack` graphic above divides Ocean software and its context into the following layers of functionality: * Compute Resources The hardware on which the problem is solved. This might be a D-Wave quantum processor but it can also be the CPU of your laptop computer. * Samplers Abstraction layer of the :term:`sampler` functionality. Ocean tools implement several samplers that use the D-Wave system and classical compute resources. You can use the Ocean tools to customize a D-Wave sampler, create your own sampler, or use existing (classical) samplers to your code as you develop it. * Sampler API Abstraction layer that represents the problem in a form that can access the selected sampler; for example, a :doc:`dimod ` binary quadratic model (BQM) class representing your problem wrapped in a :term:`minor-embedding` composite that handles the mapping between your problem's variables and the sampler's graph. * Methods Tools that help formulate a problem as binary quadratic models; for example :doc:`dwave_networkx ` (`repo `_\ ) for graph-related problems. * Application Original problem in its context ("problem space"); for example, circuit fault diagnosis attempts to identify failed logic gates during chip manufacturing. Problem-to-Solution Tool Chain ============================== As described in the :ref:`solving_problems` section, problems can be posed in a variety of formulations; the D-Wave system solves binary quadratic models. Ocean tools assist you in converting the problem from its original form to a form native to the D-Wave system and sending the compatible problem for solving. This section will familiarize you with the different tools and how you can fit them together to solve your problem. Bottom-Up Approach ------------------ One approach to envisioning how you can map your problem-solving process to Ocean software is to start from the bottom---the hardware doing the computations---and work your way up the Ocean stack to see the complete picture. This section shows how you might map each stage of the process to a layer of the Ocean stack. 1. **Compute resource** You will likely use some combination of both local classical resources and a D-Wave system in your work with Ocean software. When would you use which? * CPU/GPU: for offline testing, small problems that can be solved exactly or heuristically in a reasonable amount of time. * QPU: hard problems or for learning how to use quantum resources to solve such problems. * Hybrid of both QPU and CPU/GPU: large, complex problems that need to run classically but may benefit from having some parts allocated to a quantum computer for solution. 2. **Sampler** Your sampler provides access to the compute resource that solves your problem. The table below shows some Ocean samplers and considerations for selecting one or another. .. list-table:: Ocean Samplers :widths: 10 20 50 40 :header-rows: 1 * - Computation - Tool & Sampler - Usage - Notes * - Classical - :doc:`dimod ` :class:`~dimod.reference.samplers.ExactSolver` - Find all states for small (<20 variables) problems. - For code-development testing. * - Classical - :doc:`dimod ` :class:`~dimod.reference.samplers.random_sampler.RandomSampler` - Random sampler for testing. - For code-development testing. * - Classical - :doc:`dimod ` :class:`~dimod.reference.samplers.simulated_annealing.SimulatedAnnealingSampler` - Simulated annealing sampler for testing. - For code-development testing. * - Classical - :doc:`dwave-greedy ` :class:`~greedy.sampler.SteepestDescentSolver`. - A steepest-descent solver for binary quadratic models. - For post-processing and convex problems. * - Classical - :doc:`dwave-neal ` :class:`~neal.sampler.SimulatedAnnealingSampler` - Simulated annealing sampler. - * - Quantum - :doc:`dwave-system ` :class:`~dwave.system.samplers.DWaveSampler` - Quick incorporation of the D-Wave system as a sampler. - Typically part of a composite that handles :term:`minor-embedding`. * - Quantum - :doc:`dwave-system ` :class:`~dwave.system.samplers.DWaveCliqueSampler` - Quick incorporation of the D-Wave system as a sampler. - Handles :term:`minor-embedding` for clique (:term:`complete graph`) problems. * - Quantum - :doc:`dwave-cloud-client ` :code:`Solver()` - D-Wave system as a sampler.\ [#]_ - For low-level control of problem submission. * - Hybrid - :doc:`dwave-hybrid ` :class:`~hybrid.reference.kerberos.KerberosSampler` - *dimod*-compatible hybrid asynchronous decomposition sampler. - For problems of arbitrary structure and size. * - Hybrid - `Leap `_\ 's :class:`~dwave.system.samplers.LeapHybridSampler` - Cloud-based quantum-classical hybrid solver. - For problems of arbitrary structure and size, especially large problems. * - Hybrid - `Leap `_\ 's :class:`~dwave.system.samplers.LeapHybridCQMSampler` - Cloud-based quantum-classical hybrid solver. - For constrained quadratic models (:term:`CQM`) of arbitrary structure and size. * - - :doc:`dimod ` custom - Write a custom sampler for special cases. - See examples in :doc:`dimod `. .. [#] This sampler is for low-level work on communicating with SAPI and is not a dimod sampler. 3. **Pre- and Post-Processing** Samplers can be composed of `composite patterns `_ that layer pre- and post-processing to binary quadratic programs without changing the underlying sampler. The table below shows some Ocean composites and considerations for selecting one or another. .. list-table:: Ocean Composites :widths: 10 50 50 :header-rows: 1 * - Tool & Composite - Usage - Notes * - :doc:`dwave-system ` :class:`~dwave.system.composites.EmbeddingComposite` - Maps unstructured problems to a structured sampler. - Enables quick incorporation of the D-Wave system as a sampler by handling the :term:`minor-embedding` to the QPU's :term:`Chimera` topology of qubits. * - :doc:`dwave-system ` :class:`~dwave.system.composites.FixedEmbeddingComposite` - Maps unstructured problems to a structured sampler. - Uses a pre-calculated minor-embedding for improved performance. * - :doc:`dwave-system ` :class:`~dwave.system.composites.TilingComposite` - Tiles small problems multiple times to a Chimera-structured sampler. - Enables parallel sampling for small problems. * - :doc:`dwave-system ` :class:`~dwave.system.composites.VirtualGraphComposite` - Uses the D-Wave virtual graph feature for improved minor-embedding. - Calibrates qubits in chains to compensate for the effects of biases and enables easy creation, optimization, use, and reuse of an embedding for a given working graph. * - :std:doc:`dwave-preprocessing ` :class:`~dwave.preprocessing.composites.SpinReversalTransformComposite` - Applies spin reversal transform preprocessing. - Improves QPU results by reducing the impact of possible analog and systematic errors. * - :doc:`dimod ` :class:`~dimod.reference.composites.structure.StructureComposite` - Creates a structured composed sampler from an unstructured sampler. - Maps from a problem graph (e.g., a square graph) to a sampler's graph. In addition to composites that provide pre- and post-processing, Ocean also provides stand-alone tools to handle complex or large problems. For example: * :doc:`minorminer ` for :term:`minor-embedding` might be used to improve solutions by fine tuning parameters or incorporating problem knowledge into the embedding. * :doc:`dwave-greedy ` provides a steepest-descent solver for binary quadratic models that can be run on the samples returned from solvers such as :class:`~dwave.system.samplers.DWaveSampler` to find local minima in the neighbourhoods * :doc:`dwave-hybrid ` splits problems too large for the QPU into pieces solved either via a D-Wave system or a classical solver. 4. **Map to a Supported Format** Typically, you formulate your problem as a binary quadratic model (BQM), which you solve by submitting to the sampler (with its pre- and post-processing composite layers) you select based on the considerations listed above. Ocean provides tools for formulating the BQM: * :doc:`dwavebinarycsp ` for constraint satisfaction problems with small constraints over binary variables. For example, many problems can be posed as satisfiability problems or with Boolean logic. * :doc:`dwave_networkx ` for implementing graph-theory algorithms of the D-Wave system. Many problems can be posed in a form of graphs---this tool handles the construction of BQMs for several standard graph algorithms such as maximum cut, cover, and coloring. You might formulate a BQM mathematically; see :ref:`not` for a mathematical formulation for a two-variable problem. See the :std:doc:`system documents ` for more information on techniques for formulating problems as BQMs. 5. **Formulate** The first step in solving a problem is to express it in a mathematical formulation. For example, the :ref:`map_coloring` problem is to assign a color to each region of a map such that any two regions sharing a border have different colors. To begin solving this problem on any computer, classical or quantum, it must be concretely defined; an intuitive approach, for the map problem, is to think of the regions as variables representing the possible set of colors, the values of which must be selected from some numerical scheme, such as natural numbers. The selection function must express the problem’s constraints: * Each region is assigned one color only, of C possible colors. * The color assigned to one region cannot be assigned to adjacent regions. Now solving the problem means finding a permissible value for each of the variables. When formulating a problem for the D-Wave system, bear in mind a few considerations: * Mathematical formulations must use binary variables because the solution is implemented physically with qubits, and so must translate to spins :math:`s_i \in {−1, +1}` or equivalent binary values :math:`x_i \in {0, 1}`. * Relationships between variables must be reducible to quadratic (e.g., a QUBO) because the problem’s parameters are represented by qubits’ weights and couplers’ strengths on a QPU. * Formulations should be sparing in its number of variables because a QPU has a limited number of qubits and couplers. * Alternative formulations may have different implications for performance. Ocean demo applications, which formulate known problems, include: * `Structural Imbalance `_. * `Circuit-Fault Diagnosis `_. Top-Down Approach ----------------- Another approach to envisioning how you can map your problem-solving process to Ocean software is to start from the top---your (possibly abstractly defined) problem---and work your way down the Ocean stack. .. list-table:: Ocean Software :widths: 10 120 :header-rows: 1 * - Step - Description * - State the Problem - Define your problem concretely/mathematically; for example, as a constraint satisfaction problem or a graph problem. * - Formulate as a BQM - Reformulate an integer problem to use binary variables, for example, or convert a nonquadratic (high-order) polynomial to a QUBO. Ocean's :doc:`dwavebinarycsp ` and :doc:`dwave_networkx ` can be helpful for some problems. * - Decompose - Allocate large problems to classical and quantum resources. Ocean's :doc:`dwave-hybrid ` provides a framework and building blocks to help you create hybrid workflows. * - Embed - Consider whether your problem has repeated elements, such as logic gates, when deciding what tool to use to :term:`minor-embed` your BQM on the QPU. You might start with fully automated embedding (using :class:`~dwave.system.composites.EmbeddingComposite` for example) and then seek performance improvements through :doc:`minorminer `. * - Configure the QPU - Use spin-reversal transforms to reduce errors, for example, or examine the annealing with reverse anneal. See the :std:doc:`system documents ` for more information of features that improve performance.