Author: bugman Date: Mon Aug 9 11:52:09 2010 New Revision: 11440 URL: http://svn.gna.org/viewcvs/relax?rev=11440&view=rev Log: Added a second unit test for perfect order in compile_2nd_matrix_iso_cone(). This checks values at and close to zero, as the sinc() function is differently defined at this position. Modified: 1.3/test_suite/unit_tests/_maths_fns/test_frame_order_matrix_ops.py Modified: 1.3/test_suite/unit_tests/_maths_fns/test_frame_order_matrix_ops.py URL: http://svn.gna.org/viewcvs/relax/1.3/test_suite/unit_tests/_maths_fns/test_frame_order_matrix_ops.py?rev=11440&r1=11439&r2=11440&view=diff ============================================================================== --- 1.3/test_suite/unit_tests/_maths_fns/test_frame_order_matrix_ops.py (original) +++ 1.3/test_suite/unit_tests/_maths_fns/test_frame_order_matrix_ops.py Mon Aug 9 11:52:09 2010 @@ -149,6 +149,28 @@ self.assertAlmostEqual(f2[i, j], self.I_order[i, j]) + def test_compile_2nd_matrix_iso_cone_order2(self): + """2nd check if compile_2nd_matrix_iso_cone() can return the identity matrix for order.""" + + # Init. + R = eye(3) + z_axis = array([0, 0, 1], float64) + cone_axis = zeros(3, float64) + + # Calculate the frame order matrix. + f2 = compile_2nd_matrix_iso_cone(self.f2_temp, R, 0.0, 0.0, 0.0, 0.0, 0.0) + + # Print outs. + print_frame_order_2nd_degree(self.I_order, "Identity for order") + print_frame_order_2nd_degree(f2, "Compiled frame order") + + # Check the values. + for i in range(9): + for j in range(9): + print "Element %s, %s." % (i, j) + self.assertAlmostEqual(f2[i, j], self.I_order[i, j]) + + def test_compile_2nd_matrix_iso_cone_free_rotor_disorder(self): """Check if compile_2nd_matrix_iso_cone_free_rotor() can return the identity matrix for disorder."""