| 62 | } |
| 63 | |
| 64 | QCircuit QAdder( |
| 65 | QVec &adder1, |
| 66 | QVec &adder2, |
| 67 | Qubit* c, |
| 68 | Qubit* is_carry) |
| 69 | { |
| 70 | if ((adder1.size() == 0) || adder1.size() != adder2.size()) |
| 71 | { |
| 72 | QCERR("adder1 and adder2 must be equal, but not equal to 0!"); |
| 73 | throw run_fail("adder1 and adder2 must be equal, but not equal to 0!"); |
| 74 | } |
| 75 | |
| 76 | int nbit = adder1.size(); |
| 77 | |
| 78 | QCircuit circuit; |
| 79 | circuit << MAJ(c, adder1[0], adder2[0]); |
| 80 | |
| 81 | for (auto i = 1; i < nbit; i++) |
| 82 | { |
| 83 | circuit << MAJ(adder2[i - 1], adder1[i], adder2[i]); |
| 84 | } |
| 85 | |
| 86 | circuit << CNOT(adder2[adder2.size() - 1], is_carry); |
| 87 | |
| 88 | for (auto i = nbit - 1; i > 0; i = i - 1) |
| 89 | { |
| 90 | circuit << UMA(adder2[i - 1], adder1[i], adder2[i]); |
| 91 | } |
| 92 | |
| 93 | circuit << UMA(c, adder1[0], adder2[0]); |
| 94 | |
| 95 | return circuit; |
| 96 | |
| 97 | } |
| 98 | |
| 99 | QCircuit QAdder( |
| 100 | QVec &adder1, |