Ejemplo n.º 1
0
 def test_noreference(self):
     from pytraj.datafiles import load_cpptraj_output, tz2_ortho_trajin
     traj = pt.iterload("./data/tz2.ortho.nc", "./data/tz2.ortho.parm7")
     cout = load_cpptraj_output(tz2_ortho_trajin + """
     rmsd first @CA perres range 2-7""")
     d = pt.rmsd_perres(traj,
                        ref=0,
                        mask='@CA',
                        resrange='2-7',
                        dtype='ndarray')
     aa_eq(cout[1:].values, d)
Ejemplo n.º 2
0
 def test_noreference(self):
     from pytraj.datafiles import load_cpptraj_output, tz2_ortho_trajin
     traj = pt.iterload(fn("tz2.ortho.nc"), fn("tz2.ortho.parm7"))
     cout = load_cpptraj_output(tz2_ortho_trajin + """
     rmsd first @CA perres range 2-7""")
     d = pt.rmsd_perres(traj,
                        ref=0,
                        mask='@CA',
                        resrange='2-7',
                        dtype='ndarray')
     aa_eq(cout[1:].values, d)