def testFeedback(self): A, B, C, D, A1, A2 = get_symbols(3, 2, 1, 1, 1, 1) circuit_identity(1) # self.assertRaises(Exception, Feedback, ()) # self.assertRaises(Exception, Feedback, (C,)) # self.assertRaises(Exception, Feedback, (C + D,)) # self.assertRaises(Exception, Feedback, (C << D,)) # self.assertRaises(Exception, Feedback, (circuit_identity(n),)) # self.assertRaises(Exception, Feedback.create, (circuit_identity(0))) # self.assertEquals(Feedback.create(circuit_identity(n)), circuit_identity(n-1)) self.assertEquals(FB(A+B), A + FB(B)) smq = map_signals_circuit({2:1}, 3) # == 'cid(1) + X' self.assertEquals(smq, smq.series_inverse()) # import metapost as mp # mp.display_circuit(Feedback.apply_with_rules(smq.series_inverse() << (B + C) << smq)) # mp.display_circuit(B.feedback() + C) self.assertEquals(( smq << (B + C)).feedback(out_index = 2, in_index = 1), B.feedback() + C) print((smq << (B + C) << smq).feedback()) self.assertEquals(( smq << (B + C) << smq).feedback(), B.feedback() + C) self.assertEquals((B + C).feedback(1,1), B.feedback() + C) #check that feedback is resolved into series when possible self.assertEquals(B.feedback(1,0).substitute({B:(C+D)}), C << D) self.assertEquals((A << (B + cid(1))).feedback(), A.feedback() << B) self.assertEquals((A << (B + cid(1)) << (cid(1) + P_sigma(1,0))).feedback(2,1), A.feedback() << B) self.assertEquals((A << (cid(1) + P_sigma(1,0)) << (B + cid(1)) << (cid(1) + P_sigma(1,0))).feedback(1,1), A.feedback(1,1) << B) self.assertEquals((B << (cid(1) + C)).feedback(0,1).substitute({B: (A1 + A2)}), A2 << C << A1) self.assertEquals(((cid(1) + C)<< P_sigma(1,0) << B).feedback(1,1).substitute({B: (A1 + A2)}), A2 << C << A1) self.assertEquals(((cid(1) + C)<< P_sigma(1,0) << B << (cid(1) + D)).feedback(1,1).substitute({B: (A1 + A2)}), A2 << D<< C << A1)
def testCPermutation(self): test_perm = (0,1,2,5,6,3,4) qtp = CPermutation.create(test_perm) self.assertRaises(Exception, CPermutation.create,(list(test_perm))) self.assertEqual(qtp.series_inverse(), CPermutation.create(invert_permutation(test_perm))) self.assertEqual(qtp.block_structure, (1,1,1,4)) id1 = circuit_identity(1) self.assertEqual(qtp.get_blocks(), (id1, id1, id1, CPermutation.create((2,3,0,1)))) self.assertEqual(CPermutation((1,0,3,2)).get_blocks(), (CPermutation((1,0)), CPermutation((1,0)))) nt = len(test_perm) self.assertEqual(qtp << qtp.series_inverse(), circuit_identity(nt)) self.assertEqual(permute(list(invert_permutation(test_perm)), test_perm), list(range(nt)))
def testPermutation(self): n = 5 self.assertEqual(CPermutation.create(()), circuit_identity(0)) invalid_permutation = (1,1) self.assertRaises(Exception, CPermutation.create, (invalid_permutation,)) p_id = tuple(range(n)) self.assertEqual(CPermutation.create(p_id), circuit_identity(n)) self.assertEqual(map_signals({0:1,1:0}, 2), (1,0)) self.assertEqual(map_signals({0:5,1:0}, 6), (5,0,1,2,3,4)) self.assertEqual(map_signals({0:5,1:0, 3:2}, 6), invert_permutation(map_signals({5:0,0:1, 2:3}, 6)))
def testSeriesFilterIdentities(self): for n in (1,2,3, 10): A, B = get_symbol(n), get_symbol(n) idn = circuit_identity(n) self.assertEquals(A << idn, A) self.assertEquals(idn << A, A) self.assertEquals(SeriesProduct.create(idn, idn, A, idn, idn, B, idn, idn), A << B)
def testSeriesFilterIdentities(self): for n in (1, 2, 3, 10): A, B = get_symbol(n), get_symbol(n) idn = circuit_identity(n) self.assertEquals(A << idn, A) self.assertEquals(idn << A, A) self.assertEquals( SeriesProduct.create(idn, idn, A, idn, idn, B, idn, idn), A << B)
def testFeedback(self): A, B, C, D, A1, A2 = get_symbols(3, 2, 1, 1, 1, 1) circuit_identity(1) # self.assertRaises(Exception, Feedback, ()) # self.assertRaises(Exception, Feedback, (C,)) # self.assertRaises(Exception, Feedback, (C + D,)) # self.assertRaises(Exception, Feedback, (C << D,)) # self.assertRaises(Exception, Feedback, (circuit_identity(n),)) # self.assertRaises(Exception, Feedback.create, (circuit_identity(0))) # self.assertEquals(Feedback.create(circuit_identity(n)), circuit_identity(n-1)) self.assertEquals(FB(A + B), A + FB(B)) smq = map_signals_circuit({2: 1}, 3) # == 'cid(1) + X' self.assertEquals(smq, smq.series_inverse()) # import metapost as mp # mp.display_circuit(Feedback.apply_with_rules(smq.series_inverse() << (B + C) << smq)) # mp.display_circuit(B.feedback() + C) self.assertEquals((smq << (B + C)).feedback(out_index=2, in_index=1), B.feedback() + C) print((smq << (B + C) << smq).feedback()) self.assertEquals((smq << (B + C) << smq).feedback(), B.feedback() + C) self.assertEquals((B + C).feedback(1, 1), B.feedback() + C) #check that feedback is resolved into series when possible self.assertEquals(B.feedback(1, 0).substitute({B: (C + D)}), C << D) self.assertEquals((A << (B + cid(1))).feedback(), A.feedback() << B) self.assertEquals((A << (B + cid(1)) << (cid(1) + P_sigma(1, 0))).feedback(2, 1), A.feedback() << B) self.assertEquals((A << (cid(1) + P_sigma(1, 0)) << (B + cid(1)) << (cid(1) + P_sigma(1, 0))).feedback(1, 1), A.feedback(1, 1) << B) self.assertEquals( (B << (cid(1) + C)).feedback(0, 1).substitute({B: (A1 + A2)}), A2 << C << A1) self.assertEquals(((cid(1) + C) << P_sigma(1, 0) << B).feedback( 1, 1).substitute({B: (A1 + A2)}), A2 << C << A1) self.assertEquals( ((cid(1) + C) << P_sigma(1, 0) << B << (cid(1) + D)).feedback(1, 1).substitute({B: (A1 + A2)}), A2 << D << C << A1)
def testCPermutation(self): test_perm = (0, 1, 2, 5, 6, 3, 4) qtp = CPermutation.create(test_perm) self.assertRaises(Exception, CPermutation.create, (list(test_perm))) self.assertEqual(qtp.series_inverse(), CPermutation.create(invert_permutation(test_perm))) self.assertEqual(qtp.block_structure, (1, 1, 1, 4)) id1 = circuit_identity(1) self.assertEqual(qtp.get_blocks(), (id1, id1, id1, CPermutation.create((2, 3, 0, 1)))) self.assertEqual( CPermutation((1, 0, 3, 2)).get_blocks(), (CPermutation( (1, 0)), CPermutation((1, 0)))) nt = len(test_perm) self.assertEqual(qtp << qtp.series_inverse(), circuit_identity(nt)) self.assertEqual( permute(list(invert_permutation(test_perm)), test_perm), list(range(nt)))
def testPermutation(self): n = 5 self.assertEqual(CPermutation.create(()), circuit_identity(0)) invalid_permutation = (1, 1) self.assertRaises(Exception, CPermutation.create, (invalid_permutation, )) p_id = tuple(range(n)) self.assertEqual(CPermutation.create(p_id), circuit_identity(n)) self.assertEqual(map_signals({0: 1, 1: 0}, 2), (1, 0)) self.assertEqual(map_signals({0: 5, 1: 0}, 6), (5, 0, 1, 2, 3, 4)) self.assertEqual( map_signals({ 0: 5, 1: 0, 3: 2 }, 6), invert_permutation(map_signals({ 5: 0, 0: 1, 2: 3 }, 6)))
def testConcatenation(self): n = 4 A, B = get_symbol(n), get_symbol(n) id0 = circuit_identity(0) self.assertEquals(A+B, Concatenation(A,B)) self.assertEquals(A+B, Concatenation.create(A,B)) self.assertEquals(id0 + id0 + A + id0 + id0 + B + id0 + id0, A + B) # self.assertRaises(Exception, Concatenation, ()) # self.assertRaises(Exception, Concatenation, (A,)) self.assertEquals((A+B).block_structure, (n,n)) self.assertEquals((A+B).get_blocks((n,n)), (A,B)) #test index_in_block() self.assertEquals((A+B).index_in_block(0), (0,0)) self.assertEquals((A+B).index_in_block(1), (1,0)) self.assertEquals((A+B).index_in_block(2), (2,0)) self.assertEquals((A+B).index_in_block(3), (3,0)) self.assertEquals((A+B).index_in_block(4), (0,1)) self.assertEquals((A+B).index_in_block(5), (1,1)) self.assertEquals((A+B).index_in_block(7), (3,1))
def testConcatenation(self): n = 4 A, B = get_symbol(n), get_symbol(n) id0 = circuit_identity(0) self.assertEquals(A + B, Concatenation(A, B)) self.assertEquals(A + B, Concatenation.create(A, B)) self.assertEquals(id0 + id0 + A + id0 + id0 + B + id0 + id0, A + B) # self.assertRaises(Exception, Concatenation, ()) # self.assertRaises(Exception, Concatenation, (A,)) self.assertEquals((A + B).block_structure, (n, n)) self.assertEquals((A + B).get_blocks((n, n)), (A, B)) #test index_in_block() self.assertEquals((A + B).index_in_block(0), (0, 0)) self.assertEquals((A + B).index_in_block(1), (1, 0)) self.assertEquals((A + B).index_in_block(2), (2, 0)) self.assertEquals((A + B).index_in_block(3), (3, 0)) self.assertEquals((A + B).index_in_block(4), (0, 1)) self.assertEquals((A + B).index_in_block(5), (1, 1)) self.assertEquals((A + B).index_in_block(7), (3, 1))