SimulatorsEstimatorInterface.h

Go to the documentation of this file.

 
#pragma once
 
#ifndef __SIMULATORS_ESTIMATOR_INTERFACE_H_
#define __SIMULATORS_ESTIMATOR_INTERFACE_H_
 
#include "../Circuit/Circuit.h"
#include "../Simulators/Simulator.h"
 
namespace Estimators {

template <typename Time = Types::time_type>
class SimulatorsEstimatorInterface {
 public:
  virtual ~SimulatorsEstimatorInterface() = default;
 
  virtual bool IsInitialized() const = 0;
 
  virtual std::shared_ptr<Simulators::ISimulator> ChooseBestSimulator(
      const std::vector<std::pair<Simulators::SimulatorType,
                                  Simulators::SimulationType>> &simulatorTypes,
      std::shared_ptr<Circuits::Circuit<Time>> &dcirc, size_t &counts,
      size_t nrQubits, size_t nrCbits, size_t nrResultCbits,
      Simulators::SimulatorType &simType, Simulators::SimulationType &method,
      std::vector<bool> &executed, const long long int maxBondDim,
      const double singularValueThreshold, const std::string &truncationMode,
      const std::string &mpsSample,
      size_t maxSimulators, const std::vector<std::string> *paulis, bool multithreading = false) const = 0;
 
  static void ExecuteUpToMeasurements(
      const std::shared_ptr<Circuits::Circuit<Time>> &dcirc, size_t nrQubits,
      size_t nrCbits, size_t nrResultCbits,
      const std::shared_ptr<Simulators::ISimulator> &sim,
      std::vector<bool> &executed, size_t *curMaxBondDim = nullptr) {
    Circuits::OperationState state;
    state.AllocateBits(nrCbits);
 
    const bool hasMeasurementsOnlyAtEnd = !dcirc->HasOpsAfterMeasurements();
    const bool specialOptimizationForStatevector =
        sim->GetSimulationType() == Simulators::SimulationType::kStatevector &&
        hasMeasurementsOnlyAtEnd;
    const bool specialOptimizationForMPS =
        sim->GetSimulationType() ==
            Simulators::SimulationType::kMatrixProductState &&
        hasMeasurementsOnlyAtEnd;
 
    executed = dcirc->ExecuteNonMeasurements(sim, state, curMaxBondDim);
 
    if (!specialOptimizationForStatevector && !specialOptimizationForMPS)
      sim->SaveState();
  }
};
 
}  // namespace Estimators
 
#endif  // !__SIMULATORS_ESTIMATOR_INTERFACE_H_