def test_inner_approx(self, backend): psi = peps.computational_zeros(2, 3, backend=backend) psi.apply_circuit([ Gate('H', [], [0]), Gate('CX', [], [0,3]), Gate('H', [], [3]), ], update_option=peps.DirectUpdate(ImplicitRandomizedSVD(rank=2))) phi = peps.computational_zeros(2, 3, backend=backend) contract_option = peps.BMPS(ReducedSVD(rank=2)) self.assertTrue(backend.isclose(psi.inner(phi, contract_option), 0.5))
def test_amplitude_local_gram_qr_svd_update(self, backend): qstate = peps.computational_zeros(2, 3, backend=backend) qstate.apply_circuit([ Gate('X', [], [0]), Gate('H', [], [1]), Gate('CX', [], [0,3]), Gate('CX', [], [1,4]), Gate('S', [], [1]), ], update_option=peps.LocalGramQRSVDUpdate(rank=2)) contract_option = peps.BMPS(ReducedSVD(rank=2)) self.assertTrue(backend.isclose(qstate.amplitude([1,0,0,1,0,0], contract_option), 1/np.sqrt(2))) self.assertTrue(backend.isclose(qstate.amplitude([1,1,0,1,1,0], contract_option), 1j/np.sqrt(2)))
def test_amplitude_approx(self, backend): qstate = peps.computational_zeros(2, 3, backend=backend) qstate.apply_circuit([ Gate('X', [], [0]), Gate('H', [], [1]), Gate('CX', [], [0,3]), Gate('CX', [], [1,4]), Gate('S', [], [1]), ], update_option=peps.DirectUpdate(ImplicitRandomizedSVD(rank=2))) contract_option = peps.BMPS(ReducedSVD(rank=2)) self.assertTrue(backend.isclose(qstate.amplitude([1,0,0,1,0,0], contract_option), 1/np.sqrt(2))) self.assertTrue(backend.isclose(qstate.amplitude([1,1,0,1,1,0], contract_option), 1j/np.sqrt(2)))
def test_expectation_use_cache_approx(self, backend): qstate = peps.computational_zeros(2, 3, backend=backend) qstate.apply_circuit([ Gate('X', [], [0]), Gate('CX', [], [0,3]), Gate('H', [], [2]), ], update_option=peps.DirectUpdate(ImplicitRandomizedSVD(rank=2))) observable = 1.5 * Observable.sum([ Observable.Z(0) * 2, Observable.Z(1), Observable.Z(2) * 2, Observable.Z(3), ]) contract_option = peps.BMPS(ReducedSVD(rank=2)) self.assertTrue(backend.isclose(qstate.expectation(observable, use_cache=True, contract_option=contract_option), -3))