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 IsValid()
const {
61 if (obj)
return lib->MPSIsValid(obj);
66 bool IsCreated()
const {
67 if (obj)
return lib->MPSIsCreated(obj);
72 bool SetDataType(
int useDoublePrecision) {
73 if (obj)
return lib->MPSSetDataType(obj, useDoublePrecision);
78 bool IsDoublePrecision()
const {
79 if (obj)
return lib->MPSIsDoublePrecision(obj);
84 bool SetCutoff(
double val) {
85 if (obj)
return lib->MPSSetCutoff(obj, val);
90 double GetCutoff()
const {
91 if (obj)
return lib->MPSGetCutoff(obj);
96 bool SetGesvdJ(
int val) {
97 if (obj)
return lib->MPSSetGesvdJ(obj, val);
102 bool GetGesvdJ()
const {
103 if (obj)
return lib->MPSGetGesvdJ(obj);
108 bool SetMaxExtent(
long int val) {
109 if (obj)
return lib->MPSSetMaxExtent(obj, val);
114 long int GetMaxExtent()
const {
115 if (obj)
return lib->MPSGetMaxExtent(obj);
120 int GetNrQubits()
const {
121 if (obj)
return lib->MPSGetNrQubits(obj);
126 bool Amplitude(
long int numFixedValues,
long int *fixedValues,
double *real,
127 double *imaginary)
const {
129 return lib->MPSAmplitude(obj, numFixedValues, fixedValues, real,
135 double Probability0(
unsigned int qubit)
const {
136 if (obj)
return lib->MPSProbability0(obj, qubit);
141 bool Measure(
unsigned int qubit) {
142 if (obj)
return lib->MPSMeasure(obj, qubit);
147 bool MeasureQubits(
long int numQubits,
unsigned int *qubits,
int *result) {
148 if (obj)
return lib->MPSMeasureQubits(obj, numQubits, qubits, result);
153 std::unordered_map<std::vector<bool>, int64_t> *GetMapForSample()
const {
154 if (lib)
return lib->MPSGetMapForSample();
159 bool FreeMapForSample(
160 std::unordered_map<std::vector<bool>, int64_t> *map)
const {
161 if (lib)
return lib->MPSFreeMapForSample(map);
166 bool Sample(
long int numShots,
long int numQubits,
unsigned int *qubits,
168 if (obj)
return lib->MPSSample(obj, numShots, numQubits, qubits, resultMap);
174 if (obj)
return lib->MPSSaveState(obj);
180 if (obj)
return lib->MPSRestoreState(obj);
185 bool CleanSavedState() {
186 if (obj)
return lib->MPSCleanSavedState(obj);
191 std::unique_ptr<GpuLibMPSSim> Clone()
const {
192 if (obj)
return std::make_unique<GpuLibMPSSim>(lib, lib->MPSClone(obj));
197 double ExpectationValue(
const std::string &pauliString)
const {
199 return lib->MPSExpectationValue(obj, pauliString.c_str(),
200 pauliString.length());
205 bool ApplyX(
unsigned int siteA) {
206 if (obj)
return lib->MPSApplyX(obj, siteA);
211 bool ApplyY(
unsigned int siteA) {
212 if (obj)
return lib->MPSApplyY(obj, siteA);
217 bool ApplyZ(
unsigned int siteA) {
218 if (obj)
return lib->MPSApplyZ(obj, siteA);
223 bool ApplyH(
unsigned int siteA) {
224 if (obj)
return lib->MPSApplyH(obj, siteA);
229 bool ApplyS(
unsigned int siteA) {
230 if (obj)
return lib->MPSApplyS(obj, siteA);
236 if (obj)
return lib->MPSApplySDG(obj, siteA);
241 bool ApplyT(
unsigned int siteA) {
242 if (obj)
return lib->MPSApplyT(obj, siteA);
248 if (obj)
return lib->MPSApplyTDG(obj, siteA);
253 bool ApplySX(
unsigned int siteA) {
254 if (obj)
return lib->MPSApplySX(obj, siteA);
260 if (obj)
return lib->MPSApplySXDG(obj, siteA);
265 bool ApplyK(
unsigned int siteA) {
266 if (obj)
return lib->MPSApplyK(obj, siteA);
271 bool ApplyP(
unsigned int siteA,
double theta) {
272 if (obj)
return lib->MPSApplyP(obj, siteA, theta);
277 bool ApplyRx(
unsigned int siteA,
double theta) {
278 if (obj)
return lib->MPSApplyRx(obj, siteA, theta);
283 bool ApplyRy(
unsigned int siteA,
double theta) {
284 if (obj)
return lib->MPSApplyRy(obj, siteA, theta);
289 bool ApplyRz(
unsigned int siteA,
double theta) {
290 if (obj)
return lib->MPSApplyRz(obj, siteA, theta);
295 bool ApplyU(
unsigned int siteA,
double theta,
double phi,
double lambda,
297 if (obj)
return lib->MPSApplyU(obj, siteA, theta, phi, lambda, gamma);
302 bool ApplySwap(
unsigned int controlQubit,
unsigned int targetQubit) {
303 if (obj)
return lib->MPSApplySwap(obj, controlQubit, targetQubit);
308 bool ApplyCX(
unsigned int controlQubit,
unsigned int targetQubit) {
309 if (obj)
return lib->MPSApplyCX(obj, controlQubit, targetQubit);
314 bool ApplyCY(
unsigned int controlQubit,
unsigned int targetQubit) {
315 if (obj)
return lib->MPSApplyCY(obj, controlQubit, targetQubit);
320 bool ApplyCZ(
unsigned int controlQubit,
unsigned int targetQubit) {
321 if (obj)
return lib->MPSApplyCZ(obj, controlQubit, targetQubit);
326 bool ApplyCH(
unsigned int controlQubit,
unsigned int targetQubit) {
327 if (obj)
return lib->MPSApplyCH(obj, controlQubit, targetQubit);
332 bool ApplyCSX(
unsigned int controlQubit,
unsigned int targetQubit) {
333 if (obj)
return lib->MPSApplyCSX(obj, controlQubit, targetQubit);
338 bool ApplyCSXDG(
unsigned int controlQubit,
unsigned int targetQubit) {
339 if (obj)
return lib->MPSApplyCSXDG(obj, controlQubit, targetQubit);
344 bool ApplyCP(
unsigned int controlQubit,
unsigned int targetQubit,
346 if (obj)
return lib->MPSApplyCP(obj, controlQubit, targetQubit, theta);
351 bool ApplyCRx(
unsigned int controlQubit,
unsigned int targetQubit,
353 if (obj)
return lib->MPSApplyCRx(obj, controlQubit, targetQubit, theta);
358 bool ApplyCRy(
unsigned int controlQubit,
unsigned int targetQubit,
360 if (obj)
return lib->MPSApplyCRy(obj, controlQubit, targetQubit, theta);
365 bool ApplyCRz(
unsigned int controlQubit,
unsigned int targetQubit,
367 if (obj)
return lib->MPSApplyCRz(obj, controlQubit, targetQubit, theta);
372 bool ApplyCU(
unsigned int controlQubit,
unsigned int targetQubit,
373 double theta,
double phi,
double lambda,
double gamma) {
375 return lib->MPSApplyCU(obj, controlQubit, targetQubit, theta, phi, lambda,
382 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)