Method get_marginal

runtime/common/SampleResult.cpp:399–424

}

sample_result
sample_result::get_marginal(const std::vector<std::size_t> &marginalIndices,
                            const std::string_view registerName) const {
  const auto &counts = retrieve_result(registerName.data()).counts;
  auto mutableIndices = marginalIndices;

  std::sort(mutableIndices.begin(), mutableIndices.end());

  ExecutionResult sr;
  for (auto &[bits, count] : counts) {
    std::string newBits;
    for ([[maybe_unused]] auto &m : mutableIndices)
      newBits += "0";
    for (int counter = 0; auto &index : mutableIndices) {
      if (index > bits.size())
        throw std::runtime_error("Invalid marginal index (" +
                                 std::to_string(index) +
                                 ", size=" + std::to_string(bits.size()));

      newBits[counter++] = bits[index];
    }
    sr.appendResult(newBits, count);
  }

  return sample_result(sr);
}

void sample_result::clear() {
  sampleResults.clear();