Debug phase-aware QUnitClifford

unitaryfund · Aug 7, 2023 · 5e537e6 · 5e537e6
1 parent 16fc6e7
commit 5e537e6
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/include/qunitclifford.hpp b/include/qunitclifford.hpp
@@ -469,6 +469,10 @@ class QUnitClifford : public QInterface {
             return false;
         }
 
+        if (this == toCompare.get()) {
+            return true;
+        }
+
         return EntangleAll()->ApproxCompare(toCompare->EntangleAll(), error_tol);
     }
 
@@ -484,18 +488,21 @@ class QUnitClifford : public QInterface {
         ThrowIfQubitInvalid(t, std::string("QUnitClifford::Mtrx"));
         CliffordShard& shard = shards[t];
         shard.unit->Mtrx(mtrx, shard.mapped);
+        CombinePhaseOffsets(shard.unit);
     }
     void Phase(complex topLeft, complex bottomRight, bitLenInt t)
     {
         ThrowIfQubitInvalid(t, std::string("QUnitClifford::Phase"));
         CliffordShard& shard = shards[t];
         shard.unit->Phase(topLeft, bottomRight, shard.mapped);
+        CombinePhaseOffsets(shard.unit);
     }
     void Invert(complex topRight, complex bottomLeft, bitLenInt t)
     {
         ThrowIfQubitInvalid(t, std::string("QUnitClifford::Invert"));
         CliffordShard& shard = shards[t];
         shard.unit->Invert(topRight, bottomLeft, shard.mapped);
+        CombinePhaseOffsets(shard.unit);
     }
     void MCPhase(const std::vector<bitLenInt>& controls, complex topLeft, complex bottomRight, bitLenInt t)
     {