def test_one_atom_pdb(self):
self.o.read_pdb(DataPath+'1ATM.pdb')
axis = 'x'
frame = 0
theta=numpy.pi/2.0
#
self.o.rotate(frame,axis,theta)
result_coor = self.o.coor()
result_com = self.o.calccom(0)
print '\nresult_coor:\n'; util.printfl([result_coor]); print '\nresult_com:\n',util.printfl([result_com])
#
expected_coor = numpy.array([[[73.944, -41.652, 41.799]]], floattype)
expected_com = numpy.array([73.944, -41.652, 41.799], floattype)
self.assert_list_almost_equal(expected_coor, result_coor,3)
self.assert_list_almost_equal(expected_com, result_com,3)
def test_two_aa_pdb(self):
self.o.read_pdb(DataPath+'2AAD.pdb')
axis = 'y'
frame = 0
theta=numpy.pi/2.0
#
self.o.rotate(frame,axis,theta)
result_coor = self.o.coor()
result_com = self.o.calccom(0)
print '\nresult_coor:\n'; util.printfl([result_coor]); print '\nresult_com:\n',util.printfl([result_com])
#
expected_coor = numpy.array([[[41.652, 41.799, -73.944], [40.456, 42.563, -74.229], [40.463, 43.093, -75.667], [39.401, 43.279, -76.264], [40.336, 43.734, -73.210], [39.926, 43.168, -71.856], [39.354, 44.782, -73.67], [39.946, 44.177, -70.721], [41.647, 43.330, -76.231], [41.730, 43.852, -77.592], [42.184, 42.820, -78.617], [42.656, 43.169, -79.712], [42.648, 45.097, -77.671], [43.910, 44.816, -77.054], [42.000, 46.273, -76.970]]], floattype)
expected_com = numpy.array([41.276, 43.708, -75.680], floattype)
self.assert_list_almost_equal(expected_coor, result_coor,1)
self.assert_list_almost_equal(expected_com, result_com,2)
def test_one_atom_pdb_norotate(self):
self.o.read_pdb(DataPath+'1ATM.pdb')
frame = 0
phi=0.0
theta=0.0
psi=0.0
#
self.o.euler_rotate(frame,phi,theta,psi)
result_coor = self.o.coor()
result_com = self.o.calccom(0)
print '\nresult_coor:\n'; util.printfl([result_coor]); print '\nresult_com:\n',util.printfl([result_com])
#
expected_coor = numpy.array([[[73.944, 41.799, 41.652]]], floattype)
expected_com = numpy.array([73.944, 41.799, 41.652], floattype)
self.assert_list_almost_equal(expected_coor, result_coor,3)
self.assert_list_almost_equal(expected_com, result_com,3)
def test_one_atom_pdb(self):
self.o.read_pdb(DataPath + "1ATM.pdb")
axis = "x"
frame = 0
theta = numpy.pi / 2.0
#
self.o.general_axis_rotate(frame, theta, 0.2, 1.3, -3.5)
result_coor = self.o.coor()
result_com = self.o.calccom(0)
print "\nresult_coor:\n"
util.printfl([result_coor])
print "\nresult_com:\n", util.printfl([result_com])
#
expected_coor = numpy.array([[[-215.775, 167.484, 356.058]]], floattype)
expected_com = numpy.array([-215.775, 167.484, 356.058], floattype)
self.assert_list_almost_equal(expected_coor, result_coor, 3)
self.assert_list_almost_equal(expected_com, result_com, 3)
def test_one_atom_pdb_norotate(self):
self.o.read_pdb(DataPath + "1ATM.pdb")
axis = "x"
frame = 0
theta = 0.0
#
self.o.general_axis_rotate(frame, theta, 0.2, 1.3, -3.5)
result_coor = self.o.coor()
result_com = self.o.calccom(0)
print "\nresult_coor:\n"
util.printfl([result_coor])
print "\nresult_com:\n", util.printfl([result_com])
#
expected_coor = numpy.array([[[73.944, 41.799, 41.652]]], floattype)
expected_com = numpy.array([73.944, 41.799, 41.652], floattype)
self.assert_list_almost_equal(expected_coor, result_coor, 3)
self.assert_list_almost_equal(expected_com, result_com, 3)
def test_one_atom_pdb_norotate(self):
self.o.read_pdb(DataPath + '1ATM.pdb')
axis = 'x'
frame = 0
theta = 0.0
#
self.o.general_axis_rotate(frame, theta, 0.2, 1.3, -3.5)
result_coor = self.o.coor()
result_com = self.o.calccom(0)
print '\nresult_coor:\n'
util.printfl([result_coor])
print '\nresult_com:\n', util.printfl([result_com])
#
expected_coor = numpy.array([[[73.944, 41.799, 41.652]]], floattype)
expected_com = numpy.array([73.944, 41.799, 41.652], floattype)
self.assert_list_almost_equal(expected_coor, result_coor, 3)
self.assert_list_almost_equal(expected_com, result_com, 3)
def test_one_atom_pdb(self):
self.o.read_pdb(DataPath + '1ATM.pdb')
axis = 'x'
frame = 0
theta = numpy.pi / 2.0
#
self.o.general_axis_rotate(frame, theta, 0.2, 1.3, -3.5)
result_coor = self.o.coor()
result_com = self.o.calccom(0)
print '\nresult_coor:\n'
util.printfl([result_coor])
print '\nresult_com:\n', util.printfl([result_com])
#
expected_coor = numpy.array([[[-215.775, 167.484, 356.058]]],
floattype)
expected_com = numpy.array([-215.775, 167.484, 356.058], floattype)
self.assert_list_almost_equal(expected_coor, result_coor, 3)
self.assert_list_almost_equal(expected_com, result_com, 3)
def test_1KP8_pdb(self):
self.o.read_pdb(DataPath + '1KP8.pdb')
frame = 0
theta = 12.0
#
self.o.general_axis_rotate(frame, theta, 0, 1, 0)
result_com = self.o.calccom(0)
print '\nresult_com:\n', util.printfl([result_com])
#
expected_com = numpy.array([84.358, 0.251, -22.552], floattype)
self.assert_list_almost_equal(expected_com, result_com, 2)
def test_rna_pdb(self):
self.o.read_pdb(DataPath + 'rna.pdb')
frame = 0
theta = numpy.pi / 3.0
#
self.o.general_axis_rotate(frame, theta, 0.2, 1.3, -3.5)
result_com = self.o.calccom(0)
print '\nresult_com:\n', util.printfl([result_com])
#
expected_com = numpy.array([-30.425, -38.942, 44.267], floattype)
self.assert_list_almost_equal(expected_com, result_com, 2)
def test_1CRN_pdb(self):
self.o.read_pdb(DataPath + "1CRN.pdb")
frame = 0
theta = numpy.pi / 2.0
#
self.o.general_axis_rotate(frame, theta, 0, 1, 0)
result_com = self.o.calccom(0)
print "\nresult_com:\n", util.printfl([result_com])
#
expected_com = numpy.array([-6.978, 9.775, 9.300], floattype)
self.assert_list_almost_equal(expected_com, result_com, 2)
def test_two_aa_pdb(self):
self.o.read_pdb(DataPath + "2AAD.pdb")
frame = 0
theta = numpy.pi / 2.0
#
self.o.general_axis_rotate(frame, theta, 0.2, 1.3, -3.5)
result_com = self.o.calccom(0)
print "\nresult_com:\n", util.printfl([result_com])
#
expected_com = numpy.array([-221.139, 178.873, 343.429], floattype)
self.assert_list_almost_equal(expected_com, result_com, 2)
def test_1KP8_pdb(self):
self.o.read_pdb(DataPath + "1KP8.pdb")
frame = 0
theta = 12.0
#
self.o.general_axis_rotate(frame, theta, 0, 1, 0)
result_com = self.o.calccom(0)
print "\nresult_com:\n", util.printfl([result_com])
#
expected_com = numpy.array([84.358, 0.251, -22.552], floattype)
self.assert_list_almost_equal(expected_com, result_com, 2)
def test_rna_pdb(self):
self.o.read_pdb(DataPath + "rna.pdb")
frame = 0
theta = numpy.pi / 3.0
#
self.o.general_axis_rotate(frame, theta, 0.2, 1.3, -3.5)
result_com = self.o.calccom(0)
print "\nresult_com:\n", util.printfl([result_com])
#
expected_com = numpy.array([-30.425, -38.942, 44.267], floattype)
self.assert_list_almost_equal(expected_com, result_com, 2)
def test_two_aa_pdb(self):
self.o.read_pdb(DataPath + '2AAD.pdb')
frame = 0
theta = numpy.pi / 2.0
#
self.o.general_axis_rotate(frame, theta, 0.2, 1.3, -3.5)
result_com = self.o.calccom(0)
print '\nresult_com:\n', util.printfl([result_com])
#
expected_com = numpy.array([-221.139, 178.873, 343.429], floattype)
self.assert_list_almost_equal(expected_com, result_com, 2)
def test_1CRN_pdb(self):
self.o.read_pdb(DataPath+'1CRN.pdb')
frame = 0
theta=numpy.pi/2.0
#
self.o.euler_rotate(frame,theta,theta,theta)
result_com = self.o.calccom(0)
print '\nresult_com:\n',util.printfl([result_com])
#
expected_com = numpy.array([6.978, -9.775, 9.300], floattype)
self.assert_list_almost_equal(expected_com, result_com,2)
def test_rna_pdb(self):
self.o.read_pdb(DataPath+'rna.pdb')
axis = 'z'
frame = 0
theta=numpy.pi/2.0
#
self.o.rotate(frame,axis,theta)
result_com = self.o.calccom(0)
print '\nresult_com:\n',util.printfl([result_com])
#
expected_com = numpy.array([-4.352, -8.033, 9.231], floattype)
self.assert_list_almost_equal(expected_com, result_com,2)
def test_1KP8_pdb(self):
self.o.read_pdb(DataPath+'1KP8.pdb')
axis = 'x'
frame = 0
theta=12.0
#
self.o.rotate(frame,axis,theta)
result_com = self.o.calccom(0)
print '\nresult_com:\n',util.printfl([result_com])
#
expected_com = numpy.array([83.286, 14.288, 22.003], floattype)
self.assert_list_almost_equal(expected_com, result_com,2)
def test_1KP8_pdb(self):
self.o.read_pdb(DataPath+'1KP8.pdb')
frame = 0
phi=0.3
theta=2.0
psi=2.3
#
self.o.euler_rotate(frame,phi,theta,psi)
result_com = self.o.calccom(0)
print '\nresult_com:\n',util.printfl([result_com])
#
expected_com = numpy.array([47.025, 47.438, 56.242], floattype)
self.assert_list_almost_equal(expected_com, result_com,2)
