15#ifndef _GPU_LIB_MPS_SIM_H_
16#define _GPU_LIB_MPS_SIM_H_
28 explicit GpuLibMPSSim(
const std::shared_ptr<GpuLibrary> &lib) : lib(lib) {
30 obj = lib->CreateMPS();
35 GpuLibMPSSim(
const std::shared_ptr<GpuLibrary> &lib,
void *obj)
36 : lib(lib), obj(obj) {}
38 GpuLibMPSSim() =
delete;
39 GpuLibMPSSim(
const GpuLibMPSSim &) =
delete;
40 GpuLibMPSSim &operator=(
const GpuLibMPSSim &) =
delete;
41 GpuLibMPSSim(GpuLibMPSSim &&) =
default;
42 GpuLibMPSSim &operator=(GpuLibMPSSim &&) =
default;
45 if (lib && obj) lib->DestroyMPS(obj);
48 bool Create(
unsigned int nrQubits) {
49 if (obj)
return lib->MPSCreate(obj, nrQubits);
55 if (obj)
return lib->MPSReset(obj);
60 bool SetInitialQubitsMap(
const std::vector<long long int> &initialMap) {
61 if (obj)
return lib->MPSSetInitialQubitsMap(obj, initialMap);
66 bool SetUseOptimalMeetingPosition(
bool useOptimalMeetingPosition) {
67 if (obj)
return lib->MPSSetUseOptimalMeetingPosition(obj, useOptimalMeetingPosition);
71 bool IsValid()
const {
72 if (obj)
return lib->MPSIsValid(obj);
77 bool IsCreated()
const {
78 if (obj)
return lib->MPSIsCreated(obj);
83 bool SetDataType(
int useDoublePrecision) {
84 if (obj)
return lib->MPSSetDataType(obj, useDoublePrecision);
89 bool IsDoublePrecision()
const {
90 if (obj)
return lib->MPSIsDoublePrecision(obj);
95 bool SetCutoff(
double val) {
96 if (obj)
return lib->MPSSetCutoff(obj, val);
101 double GetCutoff()
const {
102 if (obj)
return lib->MPSGetCutoff(obj);
107 bool SetGesvdJ(
int val) {
108 if (obj)
return lib->MPSSetGesvdJ(obj, val);
113 bool GetGesvdJ()
const {
114 if (obj)
return lib->MPSGetGesvdJ(obj);
119 bool SetMaxExtent(
long int val) {
120 if (obj)
return lib->MPSSetMaxExtent(obj, val);
125 long int GetMaxExtent()
const {
126 if (obj)
return lib->MPSGetMaxExtent(obj);
131 int GetNrQubits()
const {
132 if (obj)
return lib->MPSGetNrQubits(obj);
137 bool Amplitude(
long int numFixedValues,
long int *fixedValues,
double *real,
138 double *imaginary)
const {
140 return lib->MPSAmplitude(obj, numFixedValues, fixedValues, real,
146 double Probability0(
unsigned int qubit)
const {
147 if (obj)
return lib->MPSProbability0(obj, qubit);
152 bool Measure(
unsigned int qubit) {
153 if (obj)
return lib->MPSMeasure(obj, qubit);
158 bool MeasureQubits(
long int numQubits,
unsigned int *qubits,
int *result) {
159 if (obj)
return lib->MPSMeasureQubits(obj, numQubits, qubits, result);
164 std::unordered_map<std::vector<bool>, int64_t> *GetMapForSample()
const {
165 if (lib)
return lib->MPSGetMapForSample();
170 bool FreeMapForSample(
171 std::unordered_map<std::vector<bool>, int64_t> *map)
const {
172 if (lib)
return lib->MPSFreeMapForSample(map);
177 bool Sample(
long int numShots,
long int numQubits,
unsigned int *qubits,
179 if (obj)
return lib->MPSSample(obj, numShots, numQubits, qubits, resultMap);
185 if (obj)
return lib->MPSSaveState(obj);
191 if (obj)
return lib->MPSRestoreState(obj);
196 bool CleanSavedState() {
197 if (obj)
return lib->MPSCleanSavedState(obj);
202 std::unique_ptr<GpuLibMPSSim> Clone()
const {
203 if (obj)
return std::make_unique<GpuLibMPSSim>(lib, lib->MPSClone(obj));
208 double ExpectationValue(
const std::string &pauliString) {
210 return lib->MPSExpectationValue(obj, pauliString.c_str(),
211 pauliString.length());
216 std::complex<double> ProjectOnZero() {
217 if (obj)
return lib->MPSProjectOnZero(obj);
219 return std::complex<double>(0., 0.);
222 bool ApplyX(
unsigned int siteA) {
223 if (obj)
return lib->MPSApplyX(obj, siteA);
228 bool ApplyY(
unsigned int siteA) {
229 if (obj)
return lib->MPSApplyY(obj, siteA);
234 bool ApplyZ(
unsigned int siteA) {
235 if (obj)
return lib->MPSApplyZ(obj, siteA);
240 bool ApplyH(
unsigned int siteA) {
241 if (obj)
return lib->MPSApplyH(obj, siteA);
246 bool ApplyS(
unsigned int siteA) {
247 if (obj)
return lib->MPSApplyS(obj, siteA);
253 if (obj)
return lib->MPSApplySDG(obj, siteA);
258 bool ApplyT(
unsigned int siteA) {
259 if (obj)
return lib->MPSApplyT(obj, siteA);
265 if (obj)
return lib->MPSApplyTDG(obj, siteA);
270 bool ApplySX(
unsigned int siteA) {
271 if (obj)
return lib->MPSApplySX(obj, siteA);
277 if (obj)
return lib->MPSApplySXDG(obj, siteA);
282 bool ApplyK(
unsigned int siteA) {
283 if (obj)
return lib->MPSApplyK(obj, siteA);
288 bool ApplyP(
unsigned int siteA,
double theta) {
289 if (obj)
return lib->MPSApplyP(obj, siteA, theta);
294 bool ApplyRx(
unsigned int siteA,
double theta) {
295 if (obj)
return lib->MPSApplyRx(obj, siteA, theta);
300 bool ApplyRy(
unsigned int siteA,
double theta) {
301 if (obj)
return lib->MPSApplyRy(obj, siteA, theta);
306 bool ApplyRz(
unsigned int siteA,
double theta) {
307 if (obj)
return lib->MPSApplyRz(obj, siteA, theta);
312 bool ApplyU(
unsigned int siteA,
double theta,
double phi,
double lambda,
314 if (obj)
return lib->MPSApplyU(obj, siteA, theta, phi, lambda, gamma);
319 bool ApplySwap(
unsigned int controlQubit,
unsigned int targetQubit) {
320 if (obj)
return lib->MPSApplySwap(obj, controlQubit, targetQubit);
325 bool ApplyCX(
unsigned int controlQubit,
unsigned int targetQubit) {
326 if (obj)
return lib->MPSApplyCX(obj, controlQubit, targetQubit);
331 bool ApplyCY(
unsigned int controlQubit,
unsigned int targetQubit) {
332 if (obj)
return lib->MPSApplyCY(obj, controlQubit, targetQubit);
337 bool ApplyCZ(
unsigned int controlQubit,
unsigned int targetQubit) {
338 if (obj)
return lib->MPSApplyCZ(obj, controlQubit, targetQubit);
343 bool ApplyCH(
unsigned int controlQubit,
unsigned int targetQubit) {
344 if (obj)
return lib->MPSApplyCH(obj, controlQubit, targetQubit);
349 bool ApplyCSX(
unsigned int controlQubit,
unsigned int targetQubit) {
350 if (obj)
return lib->MPSApplyCSX(obj, controlQubit, targetQubit);
355 bool ApplyCSXDG(
unsigned int controlQubit,
unsigned int targetQubit) {
356 if (obj)
return lib->MPSApplyCSXDG(obj, controlQubit, targetQubit);
361 bool ApplyCP(
unsigned int controlQubit,
unsigned int targetQubit,
363 if (obj)
return lib->MPSApplyCP(obj, controlQubit, targetQubit, theta);
368 bool ApplyCRx(
unsigned int controlQubit,
unsigned int targetQubit,
370 if (obj)
return lib->MPSApplyCRx(obj, controlQubit, targetQubit, theta);
375 bool ApplyCRy(
unsigned int controlQubit,
unsigned int targetQubit,
377 if (obj)
return lib->MPSApplyCRy(obj, controlQubit, targetQubit, theta);
382 bool ApplyCRz(
unsigned int controlQubit,
unsigned int targetQubit,
384 if (obj)
return lib->MPSApplyCRz(obj, controlQubit, targetQubit, theta);
389 bool ApplyCU(
unsigned int controlQubit,
unsigned int targetQubit,
390 double theta,
double phi,
double lambda,
double gamma) {
392 return lib->MPSApplyCU(obj, controlQubit, targetQubit, theta, phi, lambda,
399 std::shared_ptr<GpuLibrary> lib;
int ApplyK(void *sim, int qubit)
int RestoreState(void *sim)
int ApplyRx(void *sim, int qubit, double theta)
int ApplyX(void *sim, int qubit)
int ApplyU(void *sim, int qubit, double theta, double phi, double lambda, double gamma)
int ApplyCRy(void *sim, int controlQubit, int targetQubit, double theta)
int ApplyTDG(void *sim, int qubit)
int ApplyCSXDG(void *sim, int controlQubit, int targetQubit)
int ApplyS(void *sim, int qubit)
int ApplyCX(void *sim, int controlQubit, int targetQubit)
int ApplyCRz(void *sim, int controlQubit, int targetQubit, double theta)
int ApplyCP(void *sim, int controlQubit, int targetQubit, double theta)
int ApplySXDG(void *sim, int qubit)
int ApplySDG(void *sim, int qubit)
unsigned long long int Measure(void *sim, const unsigned long int *qubits, unsigned long int nrQubits)
int ApplyY(void *sim, int qubit)
double * Amplitude(void *sim, unsigned long long int outcome)
int ApplyZ(void *sim, int qubit)
int ApplyH(void *sim, int qubit)
int ApplyCY(void *sim, int controlQubit, int targetQubit)
int ApplyCU(void *sim, int controlQubit, int targetQubit, double theta, double phi, double lambda, double gamma)
int ApplySwap(void *sim, int qubit1, int qubit2)
int ApplyRy(void *sim, int qubit, double theta)
int ApplyP(void *sim, int qubit, double theta)
int ApplyCH(void *sim, int controlQubit, int targetQubit)
int ApplySX(void *sim, int qubit)
int ApplyCZ(void *sim, int controlQubit, int targetQubit)
int ApplyRz(void *sim, int qubit, double theta)
int ApplyT(void *sim, int qubit)
int ApplyCRx(void *sim, int controlQubit, int targetQubit, double theta)
int ApplyCSX(void *sim, int controlQubit, int targetQubit)