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 std::shared_ptr<Simulators::ISimulator> ChooseBestSimulator(
      const std::vector<std::pair<Simulators::SimulatorType,
                                  Simulators::SimulationType>> &simulatorTypes,
      const 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 std::string &maxBondDim,
      const std::string &singularValueThreshold, const std::string &mpsSample,
      size_t maxSimulators, const std::vector<std::string> *paulis,
      bool multithreading = false, bool dontRunCircuitStart = 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, bool multithreading) {
    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;
 
    sim->AllocateQubits(nrQubits);
    sim->Initialize();
 
    executed = dcirc->ExecuteNonMeasurements(sim, state);
 
    if (!specialOptimizationForStatevector && !specialOptimizationForMPS)
      sim->SaveState();
  }
};
 
} // namespace Estimators
 
#endif // !__SIMULATORS_ESTIMATOR_INTERFACE_H_