def test_create_one_from_pdb(self):
"""Testing destroy rigid body"""
m= IMP.kernel.Model()
ps=[]
ds=[]
for i in range(0,4):
p= IMP.kernel.Particle(m)
ps.append(p)
d= IMP.core.XYZ.setup_particle(p, IMP.algebra.Vector3D(i, 0, 0))
ds.append(d)
p0= IMP.kernel.Particle(m)
p1= IMP.kernel.Particle(m)
p01= IMP.kernel.Particle(m)
rb0= IMP.core.RigidBody.setup_particle(p0, [ps[0], ps[1]])
rb1= IMP.core.RigidBody.setup_particle(p1, [ps[2], ps[3]])
rb01= IMP.core.RigidBody.setup_particle(p01, [rb0, rb1])
print rb0, rb1, rb01
print IMP.core.RigidMember(rb0).get_internal_transformation(), IMP.core.RigidMember(rb1).get_internal_transformation()
dg= IMP.get_dependency_graph(m)
dg.show()
#IMP.show_graphviz(dg)
ordered= IMP.get_update_order(m.get_score_states())
print [o.get_name() for o in ordered]
