Add probabilities to MeasurementOutcomes

doc/source/api-reference/qibo.rst

@@ -1133,10 +1133,11 @@ Therefore, one of the three objects :class:`qibo.result.QuantumState`,
 is going to be returned by the circuit execution. The first gives acces to the final
 state and probabilities via the :meth:`qibo.result.QuantumState.state` and
 :meth:`qibo.result.QuantumState.probabilities` methods, whereas the second
-allows to retrieve the final samples and frequencies with the
-:meth:`qibo.result.MeasurementOutcomes.samples` and
-:meth:`qibo.result.MeasurementOutcomes.frequencies` methods. The
-third includes all the above instead.
+allows to retrieve the final samples, the frequencies and the probabilities (calculated
+as ``frequencies/nshots``) with the :meth:`qibo.result.MeasurementOutcomes.samples`,
+:meth:`qibo.result.MeasurementOutcomes.frequencies` and
+:meth:`qibo.result.MeasurementOutcomes.probabilities` methods respectively. The
+:class:`qibo.result.CircuitResult` object includes all the above instead.
 
 Every time some measurement is performed at the end of the execution, the result
 will be a ``CircuitResult`` unless the final state could not be represented with the

src/qibo/backends/numpy.py

@@ -426,21 +426,21 @@ def execute_circuit(self, circuit, initial_state=None, nshots=1000):
                 # execute_circuit_repeated would have been called
                 if circuit.measurements:
                     circuit._final_state = CircuitResult(
-                        state, circuit.measurements, self, nshots=nshots
+                        state, circuit.measurements, backend=self, nshots=nshots
                     )
                     return circuit._final_state
                 else:
-                    circuit._final_state = QuantumState(state, self)
+                    circuit._final_state = QuantumState(state, backend=self)
                     return circuit._final_state
 
             else:
                 if circuit.measurements:
                     circuit._final_state = CircuitResult(
-                        state, circuit.measurements, self, nshots=nshots
+                        state, circuit.measurements, backend=self, nshots=nshots
                     )
                     return circuit._final_state
                 else:
-                    circuit._final_state = QuantumState(state, self)
+                    circuit._final_state = QuantumState(state, backend=self)
                     return circuit._final_state
 
         except self.oom_error:
@@ -516,7 +516,9 @@ def execute_circuit_repeated(self, circuit, nshots, initial_state=None):
             if circuit.density_matrix:
                 final_states.append(state)
             if circuit.measurements:
-                result = CircuitResult(state, circuit.measurements, self, nshots=1)
+                result = CircuitResult(
+                    state, circuit.measurements, backend=self, nshots=1
+                )
                 sample = result.samples()[0]
                 results.append(sample)
                 if not circuit.density_matrix:
@@ -532,18 +534,18 @@ def execute_circuit_repeated(self, circuit, nshots, initial_state=None):
                 final_result = CircuitResult(
                     final_state,
                     circuit.measurements,
-                    self,
+                    backend=self,
                     samples=self.aggregate_shots(results),
                     nshots=nshots,
                 )
             else:
-                final_result = QuantumState(final_state, self)
+                final_result = QuantumState(final_state, backend=self)
             circuit._final_state = final_result
             return final_result
         else:
             final_result = MeasurementOutcomes(
                 circuit.measurements,
-                self,
+                backend=self,
                 samples=self.aggregate_shots(results),
                 nshots=nshots,
             )