18 #define QINTERFACE_TO_QALU(qReg) std::dynamic_pointer_cast<QAlu>(qReg)
19 #define QINTERFACE_TO_QPARITY(qReg) std::dynamic_pointer_cast<QParity>(qReg)
35 class QStabilizerHybrid;
78 bool TrimControls(
const std::vector<bitLenInt>& lControls, std::vector<bitLenInt>& output,
bool anti =
false);
83 for (
size_t i = 0; i <
shards.size(); ++i) {
90 for (
size_t i = 0
U; i < maxLcv; ++i) {
104 for (
size_t i = 0
U; i < maxLcv; ++i) {
109 if (shard->IsHPhase() || shard->IsHInvert()) {
112 if (shard->IsInvert()) {
115 if (!shard->IsPhase()) {
137 std::unique_ptr<complex[]> dMtrx(
new complex[4]);
144 std::conj(
shards[qubit]->gate[1]), std::conj(
shards[qubit]->gate[3]) };
146 mul2x2(dMtrx.get(), adj, out);
153 template <
typename F>
155 std::vector<real1_f>& rng, F fn)
157 for (int64_t shot = rng.size() - 1U; shot >= 0; --shot) {
158 if (rng[shot] >= partProb) {
163 for (
size_t i = 0
U; i < qPowers.size(); ++i) {
164 if (m & qPowers[i]) {
170 rng.erase(rng.begin() + shot);
179 std::vector<real1_f> rng;
181 for (
unsigned shot = 0
U; shot < shots; ++shot) {
182 rng.push_back(
Rand());
184 std::sort(rng.begin(), rng.end());
194 while (angle >= Period) {
202 if (!isGateSuppressed) {
219 real1_f correctionAngle = angle - (sector * sectorAngle);
220 if (correctionAngle >
PI_R1) {
221 correctionAngle -= Period;
223 if (correctionAngle <= -
PI_R1) {
224 correctionAngle += Period;
227 return correctionAngle;
238 if (shard->IsHPhase() || shard->IsHInvert()) {
241 if (shard->IsInvert()) {
244 if (!shard->IsPhase()) {
255 shard->gate[0
U] =
complex(angleCos, -angleSin);
256 shard->gate[3U] =
complex(angleCos, angleSin);
268 clone->RdmCloneFlush(
ONE_R1 / 2);
275 const std::vector<real1_f>& weights,
bitCapInt offset,
bool roundRz)
278 return isFloat ?
engine->ExpectationFloatsFactorizedRdm(roundRz, bits, weights)
279 :
engine->ExpectationBitsFactorizedRdm(roundRz, bits, perms, offset);
286 return isFloat ?
stabilizer->ExpectationFloatsFactorizedRdm(roundRz, bits, weights)
287 :
stabilizer->ExpectationBitsFactorizedRdm(roundRz, bits, perms, offset);
290 return isFloat ?
RdmCloneHelper()->stabilizer->ExpectationFloatsFactorizedRdm(roundRz, bits, weights)
291 :
RdmCloneHelper()->stabilizer->ExpectationBitsFactorizedRdm(roundRz, bits, perms, offset);
302 if (i != deadIndex) {
322 bool randomGlobalPhase =
true,
bool useHostMem =
false, int64_t deviceId = -1,
bool useHardwareRNG =
true,
323 bool useSparseStateVec =
false,
real1_f norm_thresh =
REAL1_EPSILON, std::vector<int64_t> devList = {},
328 bool useHostMem =
false, int64_t deviceId = -1,
bool useHardwareRNG =
true,
bool useSparseStateVec =
false,
332 useHostMem, deviceId, useHardwareRNG, useSparseStateVec, norm_thresh, devList, qubitThreshold,
378 mul2x2(dMtrx.get(), pauliZ, pMtrx);
379 return (
ONE_R1 - std::real(pMtrx[0] + pMtrx[1])) / 2;
393 CNOT(control, target);
395 CNOT(control, target);
419 return Compose(std::dynamic_pointer_cast<QStabilizerHybrid>(toCopy), start);
424 return ComposeNoClone(std::dynamic_pointer_cast<QStabilizerHybrid>(toCopy));
429 Decompose(start, std::dynamic_pointer_cast<QStabilizerHybrid>(dest));
446 engine->SetAmplitude(perm, amp);
465 bool ForceM(
bitLenInt qubit,
bool result,
bool doForce =
true,
bool doApply =
true);
479 const std::vector<bitLenInt>& controls,
bitLenInt qubitIndex,
const complex* mtrxs);
481 std::map<bitCapInt, int>
MultiShotMeasureMask(
const std::vector<bitCapInt>& qPowers,
unsigned shots);
482 void MultiShotMeasureMask(
const std::vector<bitCapInt>& qPowers,
unsigned shots,
unsigned long long* shotsArray);
515 engine->INC(toAdd, start, length);
524 engine->DEC(toSub, start, length);
533 engine->DECS(toSub, start, length, overflowIndex);
542 engine->CINC(toAdd, inOutStart, length, controls);
551 engine->INCS(toAdd, start, length, overflowIndex);
560 engine->INCDECC(toAdd, start, length, carryIndex);
610 const std::vector<bitLenInt>& controls)
616 const std::vector<bitLenInt>& controls)
622 const std::vector<bitLenInt>& controls)
628 const std::vector<bitLenInt>& controls)
634 const std::vector<bitLenInt>& controls)
641 const unsigned char* values,
bool resetValue =
true)
645 indexStart, indexLength, valueStart, valueLength, values, resetValue);
648 bitLenInt carryIndex,
const unsigned char* values)
652 indexStart, indexLength, valueStart, valueLength, carryIndex, values);
655 bitLenInt carryIndex,
const unsigned char* values)
659 indexStart, indexLength, valueStart, valueLength, carryIndex, values);
679 engine->ZeroPhaseFlip(start, length);
690 engine->SqrtSwap(qubitIndex1, qubitIndex2);
700 engine->ISqrtSwap(qubitIndex1, qubitIndex2);
704 const std::vector<bitLenInt> controls{ qubit1 };
714 ISwap(qubit1, qubit2);
727 engine->FSim(theta, phi, qubit1, qubit2);
733 return engine->ProbMask(mask, permutation);
738 return ApproxCompareHelper(std::dynamic_pointer_cast<QStabilizerHybrid>(toCompare),
false);
743 ApproxCompareHelper(std::dynamic_pointer_cast<QStabilizerHybrid>(toCompare),
true, error_tol);
749 engine->UpdateRunningNorm(norm_thresh);
764 return engine->ExpectationBitsAll(bits, offset);
769 return engine->ExpectationBitsAllRdm(roundRz, bits, offset);
775 return stabilizer->ExpectationBitsAll(bits, offset);
778 return RdmCloneHelper()->stabilizer->ExpectationBitsAll(bits, offset);
781 const std::vector<bitLenInt>& bits,
const std::vector<bitCapInt>& perms,
bitCapInt offset = 0
U)
787 return engine->ExpectationBitsFactorized(bits, perms, offset);
790 bool roundRz,
const std::vector<bitLenInt>& bits,
const std::vector<bitCapInt>& perms,
bitCapInt offset = 0
U)
800 return engine->ExpectationFloatsFactorized(bits, weights);
803 bool roundRz,
const std::vector<bitLenInt>& bits,
const std::vector<real1_f>& weights)
830 return engine->GetMaxSize();
unsigned GetConcurrencyLevel()
Definition: parallel_for.hpp:38
A "Qrack::QInterface" is an abstract interface exposing qubit permutation state vector with methods t...
Definition: qinterface.hpp:146
virtual void SetConcurrency(uint32_t threadsPerEngine)
Set the number of threads in parallel for loops, per component QEngine.
Definition: qinterface.hpp:249
virtual bitLenInt Allocate(bitLenInt length)
Allocate new "length" count of |0> state qubits at end of qubit index position.
Definition: qinterface.hpp:434
virtual bitLenInt Compose(QInterfacePtr toCopy)
Combine another QInterface with this one, after the last bit index of this one.
Definition: qinterface.hpp:338
bitLenInt qubitCount
Definition: qinterface.hpp:151
real1_f Rand()
Generate a random real number between 0 and 1.
Definition: qinterface.hpp:260
Definition: qparity.hpp:22
A "Qrack::QStabilizerHybrid" internally switched between Qrack::QStabilizer and Qrack::QEngine to max...
Definition: qstabilizerhybrid.hpp:43
bool TrimControls(const std::vector< bitLenInt > &lControls, std::vector< bitLenInt > &output, bool anti=false)
Definition: qstabilizerhybrid.cpp:270
complex GetAmplitude(bitCapInt perm)
Get the representational amplitude of a full permutation.
Definition: qstabilizerhybrid.hpp:440
void MACPhase(const std::vector< bitLenInt > &controls, complex topLeft, complex bottomRight, bitLenInt target)
Apply a single bit transformation that only effects phase, with arbitrary (anti-)control bits.
Definition: qstabilizerhybrid.cpp:1266
void MUL(bitCapInt toMul, bitLenInt inOutStart, bitLenInt carryStart, bitLenInt length)
Multiply by integer.
Definition: qstabilizerhybrid.hpp:584
bitLenInt ComposeNoClone(QInterfacePtr toCopy)
Definition: qstabilizerhybrid.hpp:422
bitLenInt maxEngineQubitCount
Definition: qstabilizerhybrid.hpp:56
void MACMtrx(const std::vector< bitLenInt > &controls, const complex *mtrx, bitLenInt target)
Apply an arbitrary single bit unitary transformation, with arbitrary (anti-)control bits.
Definition: qstabilizerhybrid.cpp:1240
void SetAmplitude(bitCapInt perm, complex amp)
Sets the representational amplitude of a full permutation.
Definition: qstabilizerhybrid.hpp:443
real1_f ProbAllRdm(bool roundRz, bitCapInt fullRegister)
Direct measure of full permutation probability, treating all ancillary qubits as post-selected T gate...
Definition: qstabilizerhybrid.cpp:378
void CISqrtSwap(const std::vector< bitLenInt > &lControls, bitLenInt qubit1, bitLenInt qubit2)
Apply an inverse square root of swap with arbitrary control bits.
Definition: qstabilizerhybrid.cpp:959
bitLenInt ComposeEither(QStabilizerHybridPtr toCopy, bool willDestroy)
Definition: qstabilizerhybrid.cpp:483
void IISwap(bitLenInt qubit1, bitLenInt qubit2)
Inverse ISwap - Swap values of two bits in register, and apply phase factor of -i if bits are differe...
Definition: qstabilizerhybrid.hpp:452
void SetTInjection(bool useGadget)
Set the option to use T-injection gadgets (off by default)
Definition: qstabilizerhybrid.hpp:337
real1_f FractionalRzAngleWithFlush(bitLenInt i, real1_f angle, bool isGateSuppressed=false)
Definition: qstabilizerhybrid.hpp:190
int64_t devID
Definition: qstabilizerhybrid.hpp:60
real1_f ProbMaskRdm(bool roundRz, bitCapInt mask, bitCapInt permutation)
Direct measure of masked permutation probability, treating all ancillary qubits as post-selected T ga...
Definition: qstabilizerhybrid.cpp:391
bitLenInt maxStateMapCacheQubitCount
Definition: qstabilizerhybrid.hpp:58
void INC(bitCapInt toAdd, bitLenInt start, bitLenInt length)
Add integer (without sign)
Definition: qstabilizerhybrid.hpp:508
bool isFinished()
Returns "false" if asynchronous work is still running, and "true" if all previously dispatched asynch...
Definition: qstabilizerhybrid.hpp:349
QStabilizerHybridPtr RdmCloneHelper()
Definition: qstabilizerhybrid.hpp:264
void FlushBuffers()
Definition: qstabilizerhybrid.cpp:251
void XMask(bitCapInt mask)
Masked X gate.
Definition: qstabilizerhybrid.cpp:992
std::vector< int64_t > deviceIDs
Definition: qstabilizerhybrid.hpp:64
bool doNormalize
Definition: qstabilizerhybrid.hpp:49
friend std::ostream & operator<<(std::ostream &os, const QStabilizerHybridPtr s)
Definition: qstabilizerhybrid.cpp:2182
void SqrtSwap(bitLenInt qubitIndex1, bitLenInt qubitIndex2)
Square root of Swap gate.
Definition: qstabilizerhybrid.hpp:682
void MULModNOut(bitCapInt toMul, bitCapInt modN, bitLenInt inStart, bitLenInt outStart, bitLenInt length)
Multiplication modulo N by integer, (out of place)
Definition: qstabilizerhybrid.hpp:594
virtual bitLenInt Allocate(bitLenInt length)
Allocate new "length" count of |0> state qubits at end of qubit index position.
Definition: qinterface.hpp:434
void AntiCSqrtSwap(const std::vector< bitLenInt > &lControls, bitLenInt qubit1, bitLenInt qubit2)
Apply a square root of swap with arbitrary (anti) control bits.
Definition: qstabilizerhybrid.cpp:943
real1_f ProbRdm(bitLenInt qubit)
Direct measure of bit probability to be in |1> state, treating all ancillary qubits as post-selected ...
Definition: qstabilizerhybrid.hpp:368
void Swap(bitLenInt qubit1, bitLenInt qubit2)
Swap values of two bits in register.
Definition: qstabilizerhybrid.cpp:897
void INCDECSC(bitCapInt toAdd, bitLenInt start, bitLenInt length, bitLenInt carryIndex)
Common driver method behind INCSC and DECSC (without overflow flag)
Definition: qstabilizerhybrid.hpp:567
void GetProbs(real1 *outputProbs)
Get the pure quantum state representation.
Definition: qstabilizerhybrid.cpp:670
void RdmCloneFlush(real1_f threshold=FP_NORM_EPSILON)
Flush non-Clifford phase gate gadgets with angle below a threshold.
Definition: qstabilizerhybrid.cpp:1936
void CSqrtSwap(const std::vector< bitLenInt > &lControls, bitLenInt qubit1, bitLenInt qubit2)
Apply a square root of swap with arbitrary control bits.
Definition: qstabilizerhybrid.cpp:927
void PhaseFlip()
Phase flip always - equivalent to Z X Z X on any bit in the QInterface.
Definition: qstabilizerhybrid.hpp:668
void FlushIfBlocked(bitLenInt control, bitLenInt target, bool isPhase=false)
Definition: qstabilizerhybrid.cpp:139
bitCapInt MAll()
Measure permutation state of all coherent bits.
Definition: qstabilizerhybrid.cpp:1528
virtual bitLenInt Compose(QInterfacePtr toCopy)
Combine another QInterface with this one, after the last bit index of this one.
Definition: qinterface.hpp:338
real1_f ExpectationBitsFactorized(const std::vector< bitLenInt > &bits, const std::vector< bitCapInt > &perms, bitCapInt offset=0U)
Get expectation value of bits, given an array of qubit weights.
Definition: qstabilizerhybrid.hpp:780
real1_f ProbAll(bitCapInt perm)
Direct measure of full permutation probability.
Definition: qstabilizerhybrid.hpp:441
bool isClifford()
Returns "true" if current state is identifiably within the Clifford set, or "false" if it is not or c...
Definition: qstabilizerhybrid.hpp:407
void FlushCliffordFromBuffers()
Definition: qstabilizerhybrid.hpp:230
void ClearAncilla(bitLenInt i)
Definition: qstabilizerhybrid.hpp:294
bool useTGadget
Definition: qstabilizerhybrid.hpp:51
QUnitCliffordPtr MakeStabilizer(bitCapInt perm=0U)
Definition: qstabilizerhybrid.cpp:98
bitCapInt IndexedADC(bitLenInt indexStart, bitLenInt indexLength, bitLenInt valueStart, bitLenInt valueLength, bitLenInt carryIndex, const unsigned char *values)
Add to entangled 8 bit register state with a superposed index-offset-based read from classical memory...
Definition: qstabilizerhybrid.hpp:647
void CDIV(bitCapInt toDiv, bitLenInt inOutStart, bitLenInt carryStart, bitLenInt length, const std::vector< bitLenInt > &controls)
Controlled division by power of integer.
Definition: qstabilizerhybrid.hpp:615
real1_f ExpectationFloatsFactorizedRdm(bool roundRz, const std::vector< bitLenInt > &bits, const std::vector< real1_f > &weights)
Get (reduced density matrix) expectation value of bits, given a (floating-point) array of qubit weigh...
Definition: qstabilizerhybrid.hpp:802
void DECS(bitCapInt toSub, bitLenInt start, bitLenInt length, bitLenInt overflowIndex)
Add a classical integer to the register, with sign and without carry.
Definition: qstabilizerhybrid.hpp:526
void MACInvert(const std::vector< bitLenInt > &controls, complex topRight, complex bottomLeft, bitLenInt target)
Apply a single bit transformation that reverses bit probability and might effect phase,...
Definition: qstabilizerhybrid.cpp:1315
void CacheEigenstate(bitLenInt target)
Definition: qstabilizerhybrid.cpp:310
void InvertBuffer(bitLenInt qubit)
Definition: qstabilizerhybrid.cpp:120
bitLenInt ancillaCount
Definition: qstabilizerhybrid.hpp:54
void DumpBuffers()
Definition: qstabilizerhybrid.hpp:81
real1_f SumSqrDiff(QInterfacePtr toCompare)
Definition: qstabilizerhybrid.hpp:736
std::vector< MpsShardPtr > shards
Definition: qstabilizerhybrid.hpp:67
bitCapInt IndexedSBC(bitLenInt indexStart, bitLenInt indexLength, bitLenInt valueStart, bitLenInt valueLength, bitLenInt carryIndex, const unsigned char *values)
Subtract from an entangled 8 bit register state with a superposed index-offset-based read from classi...
Definition: qstabilizerhybrid.hpp:654
real1_f ACProbRdm(bitLenInt control, bitLenInt target)
Definition: qstabilizerhybrid.hpp:391
bool IsBuffered()
Definition: qstabilizerhybrid.hpp:99
complex phaseFactor
Definition: qstabilizerhybrid.hpp:61
void MCMtrx(const std::vector< bitLenInt > &controls, const complex *mtrx, bitLenInt target)
Apply an arbitrary single bit unitary transformation, with arbitrary control bits.
Definition: qstabilizerhybrid.cpp:1113
QInterfacePtr MakeEngine(bitCapInt perm=0U)
Definition: qstabilizerhybrid.cpp:103
void UpdateRunningNorm(real1_f norm_thresh=REAL1_DEFAULT_ARG)
Force a calculation of the norm of the state vector, in order to make it unit length before the next ...
Definition: qstabilizerhybrid.hpp:746
bitLenInt maxAncillaCount
Definition: qstabilizerhybrid.hpp:57
virtual void UniformlyControlledSingleBit(const std::vector< bitLenInt > &controls, bitLenInt qubitIndex, const complex *mtrxs)
Apply a "uniformly controlled" arbitrary single bit unitary transformation.
Definition: qinterface.hpp:590
bool isSparse
Definition: qstabilizerhybrid.hpp:50
real1_f ProbParity(bitCapInt mask)
Overall probability of any odd permutation of the masked set of bits.
Definition: qstabilizerhybrid.cpp:1802
bitLenInt Compose(QInterfacePtr toCopy, bitLenInt start)
Definition: qstabilizerhybrid.hpp:417
bitLenInt ComposeNoClone(QStabilizerHybridPtr toCopy)
Definition: qstabilizerhybrid.hpp:421
bool TrySeparate(bitLenInt qubit)
Single-qubit TrySeparate()
Definition: qstabilizerhybrid.cpp:2141
QInterfacePtr engine
Definition: qstabilizerhybrid.hpp:62
real1_f ExpectationBitsAllRdm(bool roundRz, const std::vector< bitLenInt > &bits, bitCapInt offset=0U)
Get permutation expectation value of bits, treating all ancillary qubits as post-selected T gate gadg...
Definition: qstabilizerhybrid.hpp:766
void INCDECC(bitCapInt toAdd, bitLenInt start, bitLenInt length, bitLenInt carryIndex)
Common driver method behind INCC and DECC (without sign, with carry)
Definition: qstabilizerhybrid.hpp:553
void ISwap(bitLenInt qubit1, bitLenInt qubit2)
Swap values of two bits in register, and apply phase factor of i if bits are different.
Definition: qstabilizerhybrid.hpp:451
void X(bitLenInt q)
Definition: qstabilizerhybrid.hpp:496
bitCapInt IndexedLDA(bitLenInt indexStart, bitLenInt indexLength, bitLenInt valueStart, bitLenInt valueLength, const unsigned char *values, bool resetValue=true)
Set 8 bit register bits by a superposed index-offset-based read from classical memory.
Definition: qstabilizerhybrid.hpp:640
std::vector< QInterfaceEngine > cloneEngineTypes
Definition: qstabilizerhybrid.hpp:66
void ISwapHelper(bitLenInt qubit1, bitLenInt qubit2, bool inverse)
Definition: qstabilizerhybrid.cpp:2085
std::map< bitCapInt, int > MultiShotMeasureMask(const std::vector< bitCapInt > &qPowers, unsigned shots)
Statistical measure of masked permutation probability.
Definition: qstabilizerhybrid.cpp:1636
bitLenInt Compose(QStabilizerHybridPtr toCopy)
Definition: qstabilizerhybrid.hpp:414
void Mtrx(const complex *mtrx, bitLenInt target)
Apply an arbitrary single bit unitary transformation.
Definition: qstabilizerhybrid.cpp:1035
void MCInvert(const std::vector< bitLenInt > &controls, complex topRight, complex bottomLeft, bitLenInt target)
Apply a single bit transformation that reverses bit probability and might effect phase,...
Definition: qstabilizerhybrid.cpp:1192
bool EitherIsBuffered(bool logical)
Definition: qstabilizerhybrid.hpp:87
real1_f ExpectationFloatsFactorized(const std::vector< bitLenInt > &bits, const std::vector< real1_f > &weights)
Get expectation value of bits, given a (floating-point) array of qubit weights.
Definition: qstabilizerhybrid.hpp:794
bool M(bitLenInt q)
Definition: qstabilizerhybrid.hpp:494
real1_f ExpectationBitsFactorizedRdm(bool roundRz, const std::vector< bitLenInt > &bits, const std::vector< bitCapInt > &perms, bitCapInt offset=0U)
Get (reduced density matrix) expectation value of bits, given an array of qubit weights.
Definition: qstabilizerhybrid.hpp:789
bitCapIntOcl GetMaxSize()
Definition: qstabilizerhybrid.hpp:824
bool useHostRam
Definition: qstabilizerhybrid.hpp:48
void CUniformParityRZ(const std::vector< bitLenInt > &controls, bitCapInt mask, real1_f angle)
If the controls are set and the target qubit set parity is odd, this applies a phase factor of .
Definition: qstabilizerhybrid.hpp:486
bool GetTInjection()
Get the option to use T-injection gadgets.
Definition: qstabilizerhybrid.hpp:338
void AntiCISqrtSwap(const std::vector< bitLenInt > &lControls, bitLenInt qubit1, bitLenInt qubit2)
Apply an inverse square root of swap with arbitrary (anti) control bits.
Definition: qstabilizerhybrid.cpp:975
void ZeroPhaseFlip(bitLenInt start, bitLenInt length)
Reverse the phase of the state where the register equals zero.
Definition: qstabilizerhybrid.hpp:676
QInterfacePtr Clone()
Clone this QInterface.
Definition: qstabilizerhybrid.cpp:350
void MCPhase(const std::vector< bitLenInt > &controls, complex topLeft, complex bottomRight, bitLenInt target)
Apply a single bit transformation that only effects phase, with arbitrary control bits.
Definition: qstabilizerhybrid.cpp:1139
void INCS(bitCapInt toAdd, bitLenInt start, bitLenInt length, bitLenInt overflowIndex)
Add a classical integer to the register, with sign and without carry.
Definition: qstabilizerhybrid.hpp:544
bool IsLogicalBuffered()
Definition: qstabilizerhybrid.hpp:100
void CSwap(const std::vector< bitLenInt > &lControls, bitLenInt qubit1, bitLenInt qubit2)
Apply a swap with arbitrary control bits.
Definition: qstabilizerhybrid.cpp:911
void DIV(bitCapInt toDiv, bitLenInt inOutStart, bitLenInt carryStart, bitLenInt length)
Divide by integer.
Definition: qstabilizerhybrid.hpp:589
void SetDevice(int64_t dID)
Set the device index, if more than one device is available.
Definition: qstabilizerhybrid.hpp:814
void CPOWModNOut(bitCapInt base, bitCapInt modN, bitLenInt inStart, bitLenInt outStart, bitLenInt length, const std::vector< bitLenInt > &controls)
Controlled, raise a classical base to a quantum power, modulo N, (out of place)
Definition: qstabilizerhybrid.hpp:633
real1_f ApproxCompareHelper(QStabilizerHybridPtr toCompare, bool isDiscreteBool, real1_f error_tol=TRYDECOMPOSE_EPSILON)
Definition: qstabilizerhybrid.cpp:2002
void ISqrtSwap(bitLenInt qubitIndex1, bitLenInt qubitIndex2)
Inverse square root of Swap gate.
Definition: qstabilizerhybrid.hpp:692
void CMULModNOut(bitCapInt toMul, bitCapInt modN, bitLenInt inStart, bitLenInt outStart, bitLenInt length, const std::vector< bitLenInt > &controls)
Controlled multiplication modulo N by integer, (out of place)
Definition: qstabilizerhybrid.hpp:621
bool isBinaryDecisionTree()
Returns "true" if current state representation is definitely a binary decision tree,...
Definition: qstabilizerhybrid.hpp:411
bool EitherIsProbBuffered(bool logical)
Definition: qstabilizerhybrid.hpp:101
void CPhaseFlipIfLess(bitCapInt greaterPerm, bitLenInt start, bitLenInt length, bitLenInt flagIndex)
The 6502 uses its carry flag also as a greater-than/less-than flag, for the CMP operation.
Definition: qstabilizerhybrid.hpp:497
void Dispose(bitLenInt start, bitLenInt length)
Minimally decompose a set of contiguous bits from the separably composed unit, and discard the separa...
Definition: qstabilizerhybrid.cpp:613
bool isRoundingFlushed
Definition: qstabilizerhybrid.hpp:52
real1_f ProbMask(bitCapInt mask, bitCapInt permutation)
Direct measure of masked permutation probability.
Definition: qstabilizerhybrid.hpp:730
real1_f separabilityThreshold
Definition: qstabilizerhybrid.hpp:59
complex GetAmplitudeOrProb(bitCapInt perm, bool isProb=false)
Definition: qstabilizerhybrid.cpp:686
void Decompose(bitLenInt start, QInterfacePtr dest)
Minimally decompose a set of contiguous bits from the separably composed unit, into "destination".
Definition: qstabilizerhybrid.hpp:427
void YMask(bitCapInt mask)
Masked Y gate.
Definition: qstabilizerhybrid.cpp:1006
real1_f ExpectationBitsAll(const std::vector< bitLenInt > &bits, bitCapInt offset=0)
Get permutation expectation value of bits.
Definition: qstabilizerhybrid.hpp:758
QUnitCliffordPtr stabilizer
Definition: qstabilizerhybrid.hpp:63
void SetQuantumState(const complex *inputState)
Set an arbitrary pure quantum state representation.
Definition: qstabilizerhybrid.cpp:841
void IMULModNOut(bitCapInt toMul, bitCapInt modN, bitLenInt inStart, bitLenInt outStart, bitLenInt length)
Inverse of multiplication modulo N by integer, (out of place)
Definition: qstabilizerhybrid.hpp:599
bitLenInt deadAncillaCount
Definition: qstabilizerhybrid.hpp:55
void CMUL(bitCapInt toMul, bitLenInt inOutStart, bitLenInt carryStart, bitLenInt length, const std::vector< bitLenInt > &controls)
Controlled multiplication by integer.
Definition: qstabilizerhybrid.hpp:609
bool ForceMParity(bitCapInt mask, bool result, bool doForce=true)
Act as if is a measurement of parity of the masked set of qubits was applied, except force the (usual...
Definition: qstabilizerhybrid.cpp:1815
bitLenInt thresholdQubits
Definition: qstabilizerhybrid.hpp:53
QStabilizerHybrid(std::vector< QInterfaceEngine > eng, bitLenInt qBitCount, bitCapInt initState=0U, qrack_rand_gen_ptr rgp=nullptr, complex phaseFac=CMPLX_DEFAULT_ARG, bool doNorm=false, bool randomGlobalPhase=true, bool useHostMem=false, int64_t deviceId=-1, bool useHardwareRNG=true, bool useSparseStateVec=false, real1_f norm_thresh=REAL1_EPSILON, std::vector< int64_t > devList={}, bitLenInt qubitThreshold=0U, real1_f separation_thresh=FP_NORM_EPSILON_F)
Definition: qstabilizerhybrid.cpp:35
std::map< bitCapInt, complex > stateMapCache
Definition: qstabilizerhybrid.hpp:68
void SwitchToEngine()
Switches between CPU and GPU used modes.
Definition: qstabilizerhybrid.cpp:408
std::vector< QInterfaceEngine > engineTypes
Definition: qstabilizerhybrid.hpp:65
void CheckShots(unsigned shots, bitCapInt m, real1_f partProb, const std::vector< bitCapInt > &qPowers, std::vector< real1_f > &rng, F fn)
Definition: qstabilizerhybrid.hpp:154
void NormalizeState(real1_f nrm=REAL1_DEFAULT_ARG, real1_f norm_thresh=REAL1_DEFAULT_ARG, real1_f phaseArg=ZERO_R1_F)
Apply the normalization factor found by UpdateRunningNorm() or on the fly by a single bit gate.
Definition: qstabilizerhybrid.cpp:2128
bool isClifford(bitLenInt qubit)
Returns "true" if current qubit state is identifiably within the Clifford set, or "false" if it is no...
Definition: qstabilizerhybrid.hpp:409
void Finish()
If asynchronous work is still running, block until it finishes.
Definition: qstabilizerhybrid.hpp:340
bool CollapseSeparableShard(bitLenInt qubit)
Definition: qstabilizerhybrid.cpp:227
bool IsLogicalProbBuffered()
Definition: qstabilizerhybrid.hpp:124
friend std::istream & operator>>(std::istream &is, const QStabilizerHybridPtr s)
Definition: qstabilizerhybrid.cpp:2210
void INCDECBCDC(bitCapInt toAdd, bitLenInt start, bitLenInt length, bitLenInt carryIndex)
Common driver method behind INCSC and DECSC (without overflow flag)
Definition: qstabilizerhybrid.hpp:578
void SetPermutation(bitCapInt perm, complex phaseFac=CMPLX_DEFAULT_ARG)
Set to a specific permutation of all qubits.
Definition: qstabilizerhybrid.cpp:880
real1_f Prob(bitLenInt qubit)
Direct measure of bit probability to be in |1> state.
Definition: qstabilizerhybrid.cpp:1363
std::unique_ptr< complex[]> GetQubitReducedDensityMatrix(bitLenInt qubit)
Definition: qstabilizerhybrid.hpp:126
void DEC(bitCapInt toSub, bitLenInt start, bitLenInt length)
Add integer (without sign)
Definition: qstabilizerhybrid.hpp:517
QStabilizerHybrid(bitLenInt qBitCount, bitCapInt initState=0U, qrack_rand_gen_ptr rgp=nullptr, complex phaseFac=CMPLX_DEFAULT_ARG, bool doNorm=false, bool randomGlobalPhase=true, bool useHostMem=false, int64_t deviceId=-1, bool useHardwareRNG=true, bool useSparseStateVec=false, real1_f norm_thresh=REAL1_EPSILON, std::vector< int64_t > devList={}, bitLenInt qubitThreshold=0U, real1_f separation_thresh=FP_NORM_EPSILON_F)
Definition: qstabilizerhybrid.hpp:326
real1_f ExpectationFactorized(bool isFloat, const std::vector< bitLenInt > &bits, const std::vector< bitCapInt > &perms, const std::vector< real1_f > &weights, bitCapInt offset, bool roundRz)
Definition: qstabilizerhybrid.hpp:274
bool ForceM(bitLenInt qubit, bool result, bool doForce=true, bool doApply=true)
Act as if is a measurement was applied, except force the (usually random) result.
Definition: qstabilizerhybrid.cpp:1447
void GetQuantumState(complex *outputState)
Get the pure quantum state representation.
Definition: qstabilizerhybrid.cpp:653
void INCBCD(bitCapInt toAdd, bitLenInt start, bitLenInt length)
Add classical BCD integer (without sign)
Definition: qstabilizerhybrid.hpp:573
void Hash(bitLenInt start, bitLenInt length, const unsigned char *values)
Transform a length of qubit register via lookup through a hash table.
Definition: qstabilizerhybrid.hpp:661
void FlushH(bitLenInt qubit)
Definition: qstabilizerhybrid.cpp:129
void CombineAncillae()
Definition: qstabilizerhybrid.cpp:1831
void ZMask(bitCapInt mask)
Masked Z gate.
Definition: qstabilizerhybrid.cpp:1020
bool ApproxCompare(QInterfacePtr toCompare, real1_f error_tol=TRYDECOMPOSE_EPSILON)
Compare state vectors approximately, component by component, to determine whether this state vector i...
Definition: qstabilizerhybrid.hpp:740
void POWModNOut(bitCapInt base, bitCapInt modN, bitLenInt inStart, bitLenInt outStart, bitLenInt length)
Raise a classical base to a quantum power, modulo N, (out of place)
Definition: qstabilizerhybrid.hpp:604
void FSim(real1_f theta, real1_f phi, bitLenInt qubit1, bitLenInt qubit2)
The 2-qubit "fSim" gate, (useful in the simulation of particles with fermionic statistics)
Definition: qstabilizerhybrid.hpp:702
void SetConcurrency(uint32_t threadCount)
Set the number of threads in parallel for loops, per component QEngine.
Definition: qstabilizerhybrid.hpp:360
void CIMULModNOut(bitCapInt toMul, bitCapInt modN, bitLenInt inStart, bitLenInt outStart, bitLenInt length, const std::vector< bitLenInt > &controls)
Inverse of controlled multiplication modulo N by integer, (out of place)
Definition: qstabilizerhybrid.hpp:627
bitLenInt Compose(QInterfacePtr toCopy)
Combine another QInterface with this one, after the last bit index of this one.
Definition: qstabilizerhybrid.hpp:415
real1_f CProbRdm(bitLenInt control, bitLenInt target)
Definition: qstabilizerhybrid.hpp:382
void CINC(bitCapInt toAdd, bitLenInt inOutStart, bitLenInt length, const std::vector< bitLenInt > &controls)
Add integer (without sign, with controls)
Definition: qstabilizerhybrid.hpp:535
void PhaseFlipIfLess(bitCapInt greaterPerm, bitLenInt start, bitLenInt length)
This is an expedient for an adaptive Grover's search for a function's global minimum.
Definition: qstabilizerhybrid.hpp:502
std::vector< real1_f > GenerateShotProbs(unsigned shots)
Definition: qstabilizerhybrid.hpp:177
void Dump()
If asynchronous work is still running, let the simulator know that it can be aborted.
Definition: qstabilizerhybrid.hpp:351
void INCDECSC(bitCapInt toAdd, bitLenInt start, bitLenInt length, bitLenInt overflowIndex, bitLenInt carryIndex)
Common driver method behind INCSC and DECSC (with overflow flag)
Definition: qstabilizerhybrid.hpp:562
int64_t GetDeviceID()
Definition: qstabilizerhybrid.hpp:822
bool IsProbBuffered()
Definition: qstabilizerhybrid.hpp:123
Half-precision floating-point type.
Definition: half.hpp:2222
virtual void INCDECC(bitCapInt toAdd, bitLenInt start, bitLenInt length, bitLenInt carryIndex)
Common driver method behind INCC and DECC.
Definition: arithmetic.cpp:63
virtual void CINC(bitCapInt toAdd, bitLenInt inOutStart, bitLenInt length, const std::vector< bitLenInt > &controls)
Add integer (without sign, with controls)
Definition: arithmetic.cpp:114
virtual void INCS(bitCapInt toAdd, bitLenInt start, bitLenInt length, bitLenInt overflowIndex)
Add a classical integer to the register, with sign and without carry.
Definition: arithmetic.cpp:167
virtual void DEC(bitCapInt toSub, bitLenInt start, bitLenInt length)
Subtract classical integer (without sign)
Definition: arithmetic.cpp:57
virtual void INC(bitCapInt toAdd, bitLenInt start, bitLenInt length)
Add integer (without sign)
Definition: arithmetic.cpp:24
virtual void DECS(bitCapInt toSub, bitLenInt start, bitLenInt length, bitLenInt overflowIndex)
Subtract a classical integer from the register, with sign and without carry.
Definition: arithmetic.cpp:182
virtual void CNOT(bitLenInt control, bitLenInt target)
Controlled NOT gate.
Definition: qinterface.hpp:672
virtual bool M(bitLenInt qubitIndex)
Measurement gate.
Definition: qinterface.hpp:976
virtual void UniformlyControlledSingleBit(const std::vector< bitLenInt > &controls, bitLenInt qubitIndex, const complex *mtrxs)
Apply a "uniformly controlled" arbitrary single bit unitary transformation.
Definition: qinterface.hpp:590
virtual void X(bitLenInt qubit)
X gate.
Definition: qinterface.hpp:1054
virtual void AntiCNOT(bitLenInt control, bitLenInt target)
Anti controlled NOT gate.
Definition: qinterface.hpp:683
virtual void U(bitLenInt target, real1_f theta, real1_f phi, real1_f lambda)
General unitary gate.
Definition: rotational.cpp:18
virtual real1_f ExpectationBitsAll(const std::vector< bitLenInt > &bits, bitCapInt offset=0U)
Get permutation expectation value of bits.
Definition: qinterface.hpp:2394
bitCapIntOcl GetMaxSize()
Get maximum number of amplitudes that can be allocated on current device.
Definition: qinterface.hpp:2704
virtual real1_f ExpectationFloatsFactorized(const std::vector< bitLenInt > &bits, const std::vector< real1_f > &weights)
Get expectation value of bits, given a (floating-point) array of qubit weights.
Definition: qinterface.cpp:510
virtual real1_f ExpectationBitsFactorized(const std::vector< bitLenInt > &bits, const std::vector< bitCapInt > &perms, bitCapInt offset=0U)
Get expectation value of bits, given an array of qubit weights.
Definition: qinterface.cpp:469
Definition: complex16x2simd.hpp:25
std::complex< half_float::half > complex
Definition: qrack_types.hpp:62
@ QINTERFACE_OPTIMAL_BASE
Definition: qinterface.hpp:129
std::shared_ptr< QInterface > QInterfacePtr
Definition: qinterface.hpp:28
QRACK_CONST real1_f TRYDECOMPOSE_EPSILON
Definition: qrack_types.hpp:244
std::shared_ptr< QStabilizerHybrid > QStabilizerHybridPtr
Definition: qstabilizerhybrid.hpp:35
constexpr real1_f ZERO_R1_F
Definition: qrack_types.hpp:152
constexpr real1_f FP_NORM_EPSILON_F
Definition: qrack_types.hpp:245
const real1 ONE_R1
Definition: qrack_types.hpp:153
void mul2x2(complex const *left, complex const *right, complex *out)
Definition: functions.cpp:97
constexpr real1_f ONE_R1_F
Definition: qrack_types.hpp:154
half_float::half real1
Definition: qrack_types.hpp:63
QRACK_CONST real1 FP_NORM_EPSILON
Definition: qrack_types.hpp:243
bitCapInt pow2(const bitLenInt &p)
Definition: qrack_functions.hpp:22
void reverse(BidirectionalIterator first, BidirectionalIterator last, bitCapInt stride)
std::shared_ptr< QUnitClifford > QUnitCliffordPtr
Definition: qunitclifford.hpp:19
double norm(const complex2 &c)
Definition: complex16x2simd.hpp:101
const real1 REAL1_DEFAULT_ARG
Definition: qrack_types.hpp:155
const real1 PI_R1
Definition: qrack_types.hpp:158
QRACK_CONST complex ONE_CMPLX
Definition: qrack_types.hpp:239
const real1 ZERO_R1
Definition: qrack_types.hpp:151
float real1_f
Definition: qrack_types.hpp:64
float real1_s
Definition: qrack_types.hpp:65
QRACK_CONST complex CMPLX_DEFAULT_ARG
Definition: qrack_types.hpp:242
std::shared_ptr< MpsShard > MpsShardPtr
Definition: mpsshard.hpp:18
QRACK_CONST complex I_CMPLX
Definition: qrack_types.hpp:241
const real1 REAL1_EPSILON
Definition: qrack_types.hpp:157
QRACK_CONST complex ZERO_CMPLX
Definition: qrack_types.hpp:240
HALF_CONSTEXPR half abs(half arg)
Absolute value.
Definition: half.hpp:2975
half sin(half arg)
Sine function.
Definition: half.hpp:3885
half cos(half arg)
Cosine function.
Definition: half.hpp:3922
half round(half arg)
Nearest integer.
Definition: half.hpp:4495
half exp(half arg)
Exponential function.
Definition: half.hpp:3201
#define QRACK_CONST
Definition: qrack_types.hpp:150
#define bitLenInt
Definition: qrack_types.hpp:44
#define qrack_rand_gen_ptr
Definition: qrack_types.hpp:146
#define bitCapInt
Definition: qrack_types.hpp:105
#define bitCapIntOcl
Definition: qrack_types.hpp:91
#define QINTERFACE_TO_QALU(qReg)
Definition: qstabilizerhybrid.hpp:18
#define QINTERFACE_TO_QPARITY(qReg)
Definition: qstabilizerhybrid.hpp:19
Definition: qstabilizerhybrid.hpp:23
QUnitCliffordAmp(complex a, QUnitCliffordPtr s)
Definition: qstabilizerhybrid.hpp:27
QUnitCliffordPtr stabilizer
Definition: qstabilizerhybrid.hpp:25
complex amp
Definition: qstabilizerhybrid.hpp:24