S(hierarchy=all, molecule="Seh1")) add_distance_restraint(S(hierarchy=all, molecule="Nup145C", residue_indexes=[(0,423)]), S(hierarchy=all, molecule="Sec13")) # now do the actual work (m,all)= create_representation() create_restraints(m, all) # we can get the full dependency graph for the whole model with all the restraints # but it is pretty complex dg= IMP.get_dependency_graph(m) IMP.show_graphviz(dg); # better thing to do is to get the "pruned" graph pdg= IMP.get_pruned_dependency_graph(m) try: # these all open new windows which must be closed to continue # also, the graph is no where near as nice as displayed by # IMP.show_graphviz below #import matplotlib # the engine to be used must be selected before pyplot is imported #matplotlib.use("macosx") #import matplotlib.pyplot as plt # the method below requires the altgraph python package #xg=IMP.get_networkx_graph(pdg) #import networkx #networkx.draw_spectral(xg)