| 244 | } |
| 245 | |
| 246 | void Encode::amplitude_encode(const QVec &q, const std::vector<double>& data) |
| 247 | { |
| 248 | |
| 249 | vector<double>data_temp(data); |
| 250 | |
| 251 | if (!_check_normalized(data_temp)) { |
| 252 | throw run_fail("Data is not normalized"); |
| 253 | } |
| 254 | |
| 255 | if (data_temp.size() > (1 << q.size())) |
| 256 | { |
| 257 | throw run_fail("Amplitude_encode parameter error."); |
| 258 | } |
| 259 | |
| 260 | QVec qubits; |
| 261 | int k = 0; |
| 262 | |
| 263 | for (auto i : q) { |
| 264 | if (k >= std::ceil(std::log2(data.size()))) break; |
| 265 | qubits.push_back(i); |
| 266 | ++k; |
| 267 | } |
| 268 | |
| 269 | while (data_temp.size() < (1 << qubits.size())) |
| 270 | { |
| 271 | data_temp.push_back(0); |
| 272 | } |
| 273 | |
| 274 | std::vector<std::vector<double>>betas(qubits.size()); |
| 275 | std::vector<double>input(data_temp); |
| 276 | _recursive_compute_beta(data_temp, betas, (int)(qubits.size() - 1)); |
| 277 | _generate_circuit(betas, qubits); |
| 278 | |
| 279 | for (int i = 0; i < std::ceil(std::log2(data.size())); ++i) |
| 280 | { |
| 281 | m_out_qubits.push_back(q[i]); |
| 282 | } |
| 283 | |
| 284 | return; |
| 285 | } |
| 286 | void Encode::amplitude_encode(const QVec &q, const std::vector<complex<double>>& data) |
| 287 | { |
| 288 | vector<complex<double>> data_temp(data); |