def test5(): ''' Null test ''' a = Glycan() s, o = overlap_score(a, a) test(5, 0, s)
def test3(): ''' Identity test, size 3 ''' VIM = "alpha-D-NeupAc-(2->3)-beta-D-Galp-(1->4)-beta-D-GlcpNAc(1->" b = Glycan(VIM) s, o = overlap_score(b, b) test(3, 3, s)
def test8(): ''' Long test glycan ''' # https://glytoucan.org/Structures/Glycans/G27293OK iupac = '''alpha-D-Manp-(1->3)beta-D-GlcpNAc-(1->4)[beta-D-Galp-(1->4)[alpha-L-Fucp-(1->3)]-beta-D-GlcpNAc-(1->2)-alpha-D-Manp-(1->6)]-beta-D-Manp-(1->4)-beta-D-GlcpNAc-(1->4)[alpha-L-Fucp-(1->6)]-beta-D-GlcpNAc(1->''' a = Glycan(iupac) s, o = overlap_score(a, a) test(8, 10, s)
def test6(): ''' Null test ''' a = Glycan() VIM = "alpha-D-NeupAc-(2->3)-beta-D-Galp-(1->4)-beta-D-GlcpNAc(1->" b = Glycan(VIM) s, o = overlap_score(a, b) test(6, 0, s) s, o = overlap_score(b, a) test(6.1, 0, s)
def test7(): ''' Singleton test ''' a = Glycan('beta-D-GlcpNAc(1->') VIM = "alpha-D-NeupAc-(2->3)-beta-D-Galp-(1->4)-beta-D-GlcpNAc(1->" b = Glycan(VIM) s, o = overlap_score(a, b) test(7, 1, s) s, o = overlap_score(b, a) test(7.1, 1, s)
def test4(): ''' Small glycans with extra crap ''' lactosamine = "beta-D-Galp-(1->4)-beta-D-GlcpNAc(1->" VIM = "alpha-D-NeupAc-(2->3)-beta-D-Galp-(1->4)-beta-D-GlcpNAc(1->" test21 = "alpha-D-NeupAc-(2->3)-" + lactosamine test22 = "beta-bargledook-(2->5)-" + VIM a = Glycan(test21) b = Glycan(test22) s, o = overlap_score(a, b) test(4, 3, s)
def test10(): ''' Long against short not subtree ''' iupac = '''alpha-D-Manp-(1->3)beta-D-GlcpNAc-(1->4)[beta-D-Galp-(1->4)[alpha-L-Fucp-(1->3)]-beta-D-GlcpNAc-(1->2)-alpha-D-Manp-(1->6)]-beta-D-Manp-(1->4)-beta-D-GlcpNAc-(1->4)[alpha-L-Fucp-(1->6)]-beta-D-GlcpNAc(1->''' a = Glycan(iupac) VIM = "alpha-D-NeupAc-(2->3)-beta-D-Galp-(1->4)-beta-D-GlcpNAc(1->" b = Glycan(VIM) s, o = overlap_score(a, b) test(10, 2, s) s, o = overlap_score(b, a) test(10.1, 2, s)
def test9(): ''' Long against short subtree ''' iupac = '''alpha-D-Manp-(1->3)beta-D-GlcpNAc-(1->4)[beta-D-Galp-(1->4)[alpha-L-Fucp-(1->3)]-beta-D-GlcpNAc-(1->2)-alpha-D-Manp-(1->6)]-beta-D-Manp-(1->4)-beta-D-GlcpNAc-(1->4)[alpha-L-Fucp-(1->6)]-beta-D-GlcpNAc(1->''' a = Glycan(iupac) lactosamine = "beta-D-Galp-(1->4)-beta-D-GlcpNAc(1->" b = Glycan(lactosamine) s, o = overlap_score(a, b) test(9, 2, s) s, o = overlap_score(b, a) test(9.1, 2, s)
def test1(): ''' Two small glycans ''' #https://glytoucan.org/Structures/Glycans/G00055MO #https://glytoucan.org/Structures/Glycans/G00065MO lactosamine = "beta-D-Galp-(1->4)-beta-D-GlcpNAc(1->" a = Glycan(lactosamine) VIM = "alpha-D-NeupAc-(2->3)-beta-D-Galp-(1->4)-beta-D-GlcpNAc(1->" b = Glycan(VIM) s, o = overlap_score(a, b) test(1, 2, s)
def test11(): ''' Non-Overlapping Fragments ''' a = Glycan(['L-Fucp', 'D-Galp'], [{ 1: ('alpha', 1, 2) }, { 0: ('alpha', 2, 1) }]) b = Glycan(['L-Fucp', 'D-GlcpNAc'], [{ 1: ('alpha', 1, 3) }, { 0: ('alpha', 3, 1) }]) s, o = overlap_score(a, b) test(11, 0, s)
def test3_5(verbose=False): ''' Identity test, switching the indices ''' lactosamine = "beta-D-Galp-(1->4)-beta-D-GlcpNAc(1->" a = Glycan(lactosamine) b = Glycan(lactosamine) b.names.append(b.names.pop(0)) b.bonds.append(b.bonds.pop(0)) b.bonds[0][1] = b.bonds[0].pop(0) b.bonds[1][0] = b.bonds[1].pop(1) if verbose: print('Same glycan, different ordering:') print(a) print(b) s, o = overlap_score(a, b) test(3.5, 2, s)
def test2(): ''' Combining two non-subtree glycans https://glytoucan.org/Structures/Glycans/G00053MO https://glytoucan.org/Structures/Glycans/G00047MO https://glytoucan.org/Structures/Glycans/G00065MO ''' sialyl_lewis_a = "alpha-D-NeupAc-(2->3)-beta-D-Galp-(1->3)[alpha-L-Fucp-(1->4)]-beta-D-GlcpNAc(1->" lewis_a = "beta-D-Galp-(1->3)[alpha-L-Fucp-(1->4)]-?-D-GlcpNAc(1->" notVIM = "alpha-D-NeupAc-(2->3)-beta-D-Galp-(1->3)-beta-D-GlcpNAc(1->" a = Glycan(sialyl_lewis_a) b = Glycan(lewis_a) c = Glycan(notVIM) s, o = overlap_score(b, c) #print(len(b),len(c),s) d = overlap_combine(b, c, o.pop()) #s,o = overlap_score(a,d) #print(len(a), len(d), s) test(2, a, d)
def test12(): ''' Non-Overlapping Fragments with nuances ''' a = Glycan(['L-Fucp', 'D-Galp', 'D-Manp'], [{ 1: ('alpha', 1, 2) }, { 0: ('alpha', 2, 1), 2: ('beta', 1, 4) }, { 1: ('beta', 4, 1) }]) b = Glycan(['L-Fucp', 'D-Galp', 'D-GlcpNAc'], [{ 1: ('alpha', 1, 2) }, { 0: ('alpha', 2, 1), 2: ('beta', 1, 4) }, { 1: ('beta', 4, 1) }]) s, o = overlap_score(a, b) test(12, 0, s)
#!/usr/bin/python import simple_test simple_test.test("test23", [ "blah", "--blah", "-s=bill", "-i=9", "-i=8", "-B", "homer", "marge", "bart", ])
#!/usr/bin/python import simple_test simple_test.test("test5", [ "-b", "asdf", "-c", "fdas", "-g", "asdf", "-j", "homer", ])
#!/usr/bin/python import simple_test simple_test.test("test38", [])
#!/usr/bin/python import simple_test simple_test.test("test7", ["-n", "mike", "2", "1", ], expect_fail=True)
#!/usr/bin/python import simple_test simple_test.test("test8", [ "-s=one", "homer", "-B", ], expect_fail=True)
#!/usr/bin/python import simple_test simple_test.test("test12", [ "-v", "a 1 0.3", ], expect_fail=True)
#!/usr/bin/python import simple_test simple_test.test("test1", [], expect_fail=True)
#!/usr/bin/python import simple_test simple_test.test("test5", [ "-a", "asdf", "-c", "fdas", "--eee", "blah", "-i", "sss", "-i", "fdsf", ])
#!/usr/bin/python import simple_test simple_test.test("test42", [ "-h", ], head=6)
#!/usr/bin/python import simple_test simple_test.test("test20", [ "-b", "-a", ], expect_fail=True)
#!/usr/bin/python import simple_test simple_test.test("test31", [])
#!/usr/bin/python import simple_test simple_test.test("test8", ["-s", "bbb", "homer", "marge", "bart", "--", "-hv", "two", ], expect_fail=True)
#!/usr/bin/python import simple_test simple_test.test("test47", ['-p', '1 2.3'])
#!/usr/bin/python import simple_test simple_test.test("test4", [ "--help", ])
#!/usr/bin/python import simple_test simple_test.test("test1", [ "-n", "mike", "-r", ])
#!/usr/bin/python import simple_test simple_test.test("test33", [])
#!/usr/bin/python import simple_test simple_test.test("test8", [ "-f=9", "-f=1.0.0", "-s=asdf", "asdf", "asdf", ], expect_fail=True)