︠4671c85c-d778-427b-8096-b5e7eec13e19︠
def rcd(g, ancilla): 
        
        # Maybe some Check Up here if g is of correct format        

        # Creating Graph
     G = Graph(g, boundary = ancilla)

        # Required Matrix and Vector Space, p-list and Field
     F = GF(2)
     p = G.vertices()    
     MS = MatrixSpace(F, G.order())
     V = VectorSpace(F, G.order())

     for x in G.get_boundary(): p.remove(x)   
     
        # Creating subspace
     SemiSubSpace = V.basis_matrix().matrix_from_rows(p).row_space()
     SubSpace = SemiSubSpace.direct_sum(SemiSubSpace)     

        # Creating Generator Matrix and Stabilizer Space
     adj = MS(G.adjacency_matrix())
     Mgenerator = block_matrix([MS.identity_matrix(), adj], ncols = 2)
     Stabilizer = Mgenerator.row_space()        

        # Projection onto spanned subspace
     TraceOut = Stabilizer.intersection(SubSpace)


     return TraceOut
︡138a2c76-df1b-4909-82c8-39bd59a705eb︡︡
︠a8519c62-1e0b-4dfe-a3f2-3a1d552163c7︠
@interact
def _(graph = input_box(default = {0:[1], 1:[2], 2:[3], 3:[4], 5:[0]}, label = 'Insert Graph here:',type = dict), ancilla = input_box(default = [0], label = 'Ancilla Vertices', type = list)):
     p = rcd(graph, ancilla)    
     Graph(graph).plot
     print p
︡bc4c05a4-9f78-4395-9a65-0e9e60322275︡{"html": "<!--notruncate--><div padding=6 id=\"div-interact-2\"> <table width=800px height=20px bgcolor=\"#c5c5c5\"\n                 cellpadding=15><tr><td bgcolor=\"#f9f9f9\" valign=top align=left><table><tr><td align=right><font color=\"black\">Insert Graph here:&nbsp;</font></td><td><input type=\"text\" value=\"{0: [1], 1: [2], 2: [3], 3: [4], 5: [0]}\" size=80 onchange=\"interact(2, '_interact_.update(\\'2\\', \\'graph\\', 22, _interact_.standard_b64decode(\\''+encode64(this.value)+'\\'), globals()); _interact_.recompute(\\'2\\');')\"></input></td></tr>\n<tr><td align=right><font color=\"black\">Ancilla Vertices&nbsp;</font></td><td><input type=\"text\" value=\"[0]\" size=80 onchange=\"interact(2, '_interact_.update(\\'2\\', \\'ancilla\\', 23, _interact_.standard_b64decode(\\''+encode64(this.value)+'\\'), globals()); _interact_.recompute(\\'2\\');')\"></input></td></tr>\n</table><div id=\"cell-interact-2\"><?__SAGE__START>\n        <table border=0 bgcolor=\"white\" width=100% height=100%>\n        <tr><td bgcolor=\"white\" align=left valign=top><pre><?__SAGE__TEXT></pre></td></tr>\n        <tr><td  align=left valign=top><?__SAGE__HTML></td></tr>\n        </table><?__SAGE__END></div></td>\n                 </tr></table></div>\n                 "}︡