qcor.cpp 7.53 KB
Newer Older
Mccaskey, Alex's avatar
Mccaskey, Alex committed
1
#include "qcor.hpp"
2
#include "/home/project/xacc/xacc/XACC.hpp"
3
#include "qpu_handler.hpp"
Mccaskey, Alex's avatar
Mccaskey, Alex committed
4

5
#include "AcceleratorBuffer.hpp"
Mccaskey, Alex's avatar
Mccaskey, Alex committed
6
7
#include "IRProvider.hpp"
#include "XACC.hpp"
8
9
#include "xacc_service.hpp"

10
#include "FermionOperator.hpp"
11
#include "PauliOperator.hpp"
12

13
#include <regex>
Mccaskey, Alex's avatar
Mccaskey, Alex committed
14
15
16
17

using namespace xacc;

namespace qcor {
18
std::map<std::string, InstructionParameter> runtimeMap = {};
Mccaskey, Alex's avatar
Mccaskey, Alex committed
19

20
21
22
23
24
25
void Initialize(int argc, char **argv) {
  std::vector<const char *> tmp(argv, argv + argc);
  std::vector<std::string> newargv;
  for (auto &t : tmp)
    newargv.push_back(std::string(t));
  Initialize(newargv);
26
27
}
void Initialize(std::vector<std::string> argv) {
28
29
30
  argv.push_back("--logger-name");
  argv.push_back("qcor");
  xacc::Initialize(argv);
31
32
}

33
34
35
36
const std::string persistCompiledCircuit(std::shared_ptr<Function> function,
                                         std::shared_ptr<Accelerator> acc) {
  srand(time(NULL));
  std::function<char()> randChar = []() -> char {
37
38
39
40
41
42
43
44
45
46
47
48
49
    const char charset[] = "0123456789"
                           "ABCDEFGHIJKLMNOPQRSTUVWXYZ"
                           "abcdefghijklmnopqrstuvwxyz";
    const size_t max_index = (sizeof(charset) - 1);
    return charset[rand() % max_index];
  };

  auto generateRandomString = [&](const int length = 10) -> const std::string {
    std::string str(length, 0);
    std::generate_n(str.begin(), length, randChar);
    return str;
  };

50
51
52
53
54
55
56
57
58
  std::string file_name;
  if (xacc::optionExists("qcor-compiled-filename")) {
    file_name = xacc::getOption("qcor-compiled-filename");
  } else {

    file_name = generateRandomString();
    // std::cout << "Generating random string " << file_name << "\n";
  }

59
60
61
  std::stringstream ss;
  function->persist(ss);
  auto persistedFunction = ss.str();
62
63
64
65
//   std::cout << "PERSISTED\n" << persistedFunction << "\n";
  //   xacc::getCompiler("xacc-py")->translate("", function);
  //   persistedFunction = persistedFunction.substr(7,
  //   persistedFunction.length());
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
  xacc::appendCache(file_name, "compiled",
                    InstructionParameter(persistedFunction), ".qcor_cache");

  if (acc) {
    xacc::appendCache(file_name, "accelerator",
                      InstructionParameter(acc->name()), ".qcor_cache");
  } else {
    xacc::appendCache(file_name, "accelerator",
                      InstructionParameter("default-sim"), ".qcor_cache");
  }

  if (function->hasIRGenerators()) {
    xacc::appendCache(file_name, "requires-jit", InstructionParameter("true"),
                      ".qcor_cache");
  } else {
    xacc::appendCache(file_name, "requires-jit", InstructionParameter("false"),
                      ".qcor_cache");
  }

85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
  bool isDw = false;
  for (auto& inst : function->getInstructions()) {
      if (inst->name() == "dw-qmi" || inst->name() == "anneal") {
          isDw = true;
          break;
      }
  }

  if (isDw) {
    xacc::appendCache(file_name, "ir-type", InstructionParameter("anneal"),
                      ".qcor_cache");
  } else {
    xacc::appendCache(file_name, "ir-type", InstructionParameter("gate"),
                      ".qcor_cache");
  }
100
  return file_name;
Mccaskey, Alex's avatar
Mccaskey, Alex committed
101
102
}

103
std::shared_ptr<Function> loadCompiledCircuit(const std::string &fileName) {
104
  //   std::cout << "Loading Circuit " << fileName << "\n";
105
106
107
108
  auto cache = xacc::getCache(fileName, ".qcor_cache");
  if (!cache.count("compiled")) {
    xacc::error("Invalid quantum compilation cache.");
  }
109

110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
  std::shared_ptr<Accelerator> targetAccelerator;
  if (cache["accelerator"] == "default-sim") {
    // First, if this is compiled for simulation,
    // let users set the simulator at the command line
    // if they didnt, check for tnqvm, then local-ibm
    if (xacc::optionExists("accelerator")) {
      targetAccelerator = xacc::getAccelerator();
    } else if (xacc::hasAccelerator("tnqvm")) {
      targetAccelerator = xacc::getAccelerator("tnqvm");
    } else if (xacc::hasAccelerator("local-ibm")) {
      targetAccelerator = xacc::getAccelerator("local-ibm");
    }
  } else {
    auto accStr = cache["accelerator"].as<std::string>();
    if ((accStr == "tnqvm" || accStr == "local-ibm") &&
        xacc::optionExists("accelerator")) {
      targetAccelerator = xacc::getAccelerator();
    } else {
      targetAccelerator =
          xacc::getAccelerator(cache["accelerator"].as<std::string>());
    }
  }

  // If for some reason we still dont have an
  // accelerator, force them to specify at command line
  if (!targetAccelerator) {
    targetAccelerator = xacc::getAccelerator();
  }

  xacc::setAccelerator(targetAccelerator->name());

141
  auto type = cache["ir-type"].as<std::string>();
142
  auto compiled = cache["compiled"].as<std::string>();
143
144
  auto loaded =
      xacc::getService<IRProvider>("quantum")->createFunction("loaded", {}, {InstructionParameter(type)});
145
146
  std::istringstream iss(compiled);
  loaded->load(iss);
147

148
  //   std::cout << "Lodaded: " << loaded->toString() << "\n";
149
150
151
  if (cache["requires-jit"].as<std::string>() == "true") {
    auto runtimeMap = getRuntimeMap();

152
    // for (auto& kv : runtimeMap) {
153
154
    // std::cout << "Runtime: " << kv.first << ", " << kv.second.toString() <<
    // "\n";
155
    // }
156
157
158
159
160
    loaded->expandIRGenerators(runtimeMap);

    // Kick off quantum compilation
    auto qcor = xacc::getCompiler("qcor");
    loaded = qcor->compile(loaded, targetAccelerator);
161
    // std::cout << "NPARAMS: " << loaded->nParameters() << "\n";
162
163
  }

164
  //   std::cout<< "Loaded IR:\n" << loaded->toString() <<"\n";
165
  return loaded;
Mccaskey, Alex's avatar
Mccaskey, Alex committed
166
167
}

168
169
170
void storeRuntimeVariable(const std::string name, InstructionParameter param) {
  //   std::cout << "Storing Runtime Variable " << name << ", " <<
  //   param.toString() << "\n";
171
172
173
174
175
176
177
178
179
  runtimeMap.insert({name, param});
}

std::map<std::string, InstructionParameter> getRuntimeMap() {
  return runtimeMap;
}

std::future<std::shared_ptr<AcceleratorBuffer>>
submit(HandlerLambda &&totalJob) {
Mccaskey, Alex's avatar
Mccaskey, Alex committed
180
181
  // Create the QPU Handler to pass to the given
  // Handler HandlerLambda
182
183
  return std::async(std::launch::async, [=]() { // bug must be by value...
    qpu_handler handler;
Mccaskey, Alex's avatar
Mccaskey, Alex committed
184
    totalJob(handler);
185
    return handler.getResults();
Mccaskey, Alex's avatar
Mccaskey, Alex committed
186
187
188
  });
}

189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
std::shared_ptr<Optimizer> getOptimizer(const std::string &name) {
  return xacc::getService<qcor::Optimizer>(name);
}
std::shared_ptr<Optimizer>
getOptimizer(const std::string &name,
             std::map<std::string, InstructionParameter> &&options) {
  auto opt = getOptimizer(name);
  opt->setOptions(options);
  return opt;
}

std::shared_ptr<Observable> getObservable(const std::string &type,
                                          const std::string &representation) {
  using namespace xacc::quantum;
  if (type == "pauli") {
    return representation.empty()
               ? std::make_shared<PauliOperator>()
               : std::make_shared<PauliOperator>(representation);
207
  } else if (type == "fermion") {
208
    return representation.empty()
209
210
               ? std::make_shared<FermionOperator>()
               : std::make_shared<FermionOperator>(representation);
211
  } else {
212
213
214
215
216
217
218
219
220
221
222
223
    xacc::error("Invalid observable type");
    return std::make_shared<PauliOperator>();
  }
}

std::shared_ptr<Observable> getObservable() {
  return getObservable("pauli", "");
}
std::shared_ptr<Observable> getObservable(const std::string &representation) {
  return getObservable("pauli", representation);
}

224
225
226
227
228
229
std::shared_ptr<Observable>
getObservable(const std::string &type,
              std::map<std::string, InstructionParameter> &&options) {
  auto observable = xacc::getService<Observable>(type);
  observable->fromOptions(options);
  return observable;
230
}
231
232
std::shared_ptr<algorithm::Algorithm> getAlgorithm(const std::string name) {
  return xacc::getService<qcor::algorithm::Algorithm>(name);
233
234
}

Mccaskey, Alex's avatar
Mccaskey, Alex committed
235
} // namespace qcor