Debug classical shadows

unitaryfund · Nov 10, 2024 · a2af0af · a2af0af
1 parent 538c692
commit a2af0af
Showing 1 changed file with 6 additions and 10 deletions.
diff --git a/src/qunit.cpp b/src/qunit.cpp
@@ -3901,18 +3901,16 @@ void QUnit::ApplyBuffer(PhaseShardPtr phaseShard, bitLenInt control, bitLenInt t
 
         const complex& polarBottom = isAnti ? polarDiff : polarSame;
         const complex& polarTop = isAnti ? polarSame : polarDiff;
-        complex mtrx[4U];
-        TransformPhase(polarTop, polarBottom, mtrx);
         const real1_f cProb = isAnti ? (ONE_R1_F - Prob(control)) : Prob(control);
         const real1_f tProb = Prob(target);
 
         if (tProb > cProb) {
             if (tProb > (ONE_R1_F / 2)) {
-                Mtrx(mtrx, control);
+                Phase(ONE_CMPLX, polarBottom, control);
             }
         } else {
             if (cProb > (ONE_R1_F / 2)) {
-                Mtrx(mtrx, target);
+                Phase(ONE_CMPLX, polarBottom, target);
             }
         }
 
@@ -3921,18 +3919,16 @@ void QUnit::ApplyBuffer(PhaseShardPtr phaseShard, bitLenInt control, bitLenInt t
         const real1_f antiCProb = ONE_R1_F - cProb;
         if (tProb > antiCProb) {
             if (tProb > (ONE_R1_F / 2)) {
-                Mtrx(mtrx, control);
+                Phase(ONE_CMPLX, polarTop, control);
             }
         } else {
             if (antiCProb > (ONE_R1_F / 2)) {
-                Mtrx(mtrx, target);
+                Phase(ONE_CMPLX, polarTop, target);
             }
         }
 
-        X(target);
-
-        if (phaseShard->isInvert) {
-            H(target);
+        if (!phaseShard->isInvert) {
+            X(target);
         }
     }
     freezeBasis2Qb = false;