test_suite.unit_tests._lib._frame_order.test_matrix

46 """Unit tests for the lib.frame_order_matrix_ops relax module.""" 47

48 - def __init__(self, methodName='runTest'):

49 """Skip the tests if scipy is not installed. 50 51 @keyword methodName: The name of the test. 52 @type methodName: str 53 """ 54 55 # Missing module. 56 if not dep_check.scipy_module: 57 # Store in the status object. 58 status.skipped_tests.append([methodName, 'Scipy', 'unit']) 59 60 # Execute the base class method. 61 super(Test_matrix_ops, self).__init__(methodName)

62 63

64 - def setUp(self):

65 """Initialise a few data structures for the tests.""" 66 67 # Temp storage. 68 self.f2_temp = zeros((9, 9), float64) 69 self.R_temp = zeros((3, 3), float64) 70 self.z_axis = array([0, 0, 1], float64) 71 self.cone_axis = zeros(3, float64) 72 73 # Set up the identity matrices. 74 self.setup_identity() 75 76 # Set up the identity matrices for free rotors. 77 self.setup_identity_free_rotor() 78 79 # Out of frame Euler angles. 80 self.out_of_frame_alpha = 1.10714871779 81 self.out_of_frame_beta = 0.841068670568 82 self.out_of_frame_gamma = 5.81953769818

83 84

85 - def calc_Rx2_eigen_axis(self, axis_theta, axis_phi):

86 """Calculate the Kronecker product of the eigenframe rotation for the z-axis based frame.""" 87 88 # Generate the cone axis from the spherical angles. 89 spherical_to_cartesian([1.0, axis_theta, axis_phi], self.cone_axis) 90 91 # Pre-calculate the eigenframe rotation matrix. 92 two_vect_to_R(self.z_axis, self.cone_axis, self.R_temp) 93 94 # The Kronecker product of the eigenframe rotation. 95 return kron_prod(self.R_temp, self.R_temp)

96 97

98 - def calc_Rx2_eigen_full(self, eigen_alpha, eigen_beta, eigen_gamma):

99 """Calculate the Kronecker product of the eigenframe rotation for the full eigenframe.""" 100 101 # Average position rotation. 102 euler_to_R_zyz(eigen_alpha, eigen_beta, eigen_gamma, self.R_temp) 103 104 # The Kronecker product of the eigenframe rotation. 105 return kron_prod(self.R_temp, self.R_temp)

106 107

108 - def setup_identity(self):

109 """Set up a few identity matrices.""" 110 111 # The order identity matrix. 112 self.I_order = zeros((9, 9), float64) 113 for i in range(9): 114 self.I_order[i, i] = 1.0 115 116 # The disorder identity matrix. 117 data = [[0, 0, 1.0/3.0], 118 [4, 4, 1.0/3.0], 119 [8, 8, 1.0/3.0], 120 [0, 4, 1.0/3.0], 121 [4, 0, 1.0/3.0], 122 [0, 8, 1.0/3.0], 123 [8, 0, 1.0/3.0], 124 [4, 8, 1.0/3.0], 125 [8, 4, 1.0/3.0]] 126 self.I_disorder = zeros((9, 9), float64) 127 for i, j, val in data: 128 self.I_disorder[i, j] = val 129 130 # The half cone matrix. 131 data = [[0, 0, 1.0/3.0], 132 [4, 4, 1.0/3.0], 133 [8, 8, 1.0/3.0], 134 [0, 4, 1.0/3.0], 135 [4, 0, 1.0/3.0], 136 [0, 8, 1.0/3.0], 137 [8, 0, 1.0/3.0], 138 [4, 8, 1.0/3.0], 139 [8, 4, 1.0/3.0], 140 [1, 3, -0.25], 141 [3, 1, -0.25], 142 [1, 1, 0.25], 143 [3, 3, 0.25]] 144 self.f2_half_cone = zeros((9, 9), float64) 145 for i, j, val in data: 146 self.f2_half_cone[i, j] = val 147 148 # The half cone matrix rotated 90 degrees about y. 149 data = [[0, 0, 1.0/3.0], 150 [4, 4, 1.0/3.0], 151 [8, 8, 1.0/3.0], 152 [0, 4, 1.0/3.0], 153 [4, 0, 1.0/3.0], 154 [0, 8, 1.0/3.0], 155 [8, 0, 1.0/3.0], 156 [4, 8, 1.0/3.0], 157 [8, 4, 1.0/3.0], 158 [5, 7, -0.25], 159 [7, 5, -0.25], 160 [7, 7, 0.25], 161 [5, 5, 0.25]] 162 self.f2_half_cone_90_y = zeros((9, 9), float64) 163 for i, j, val in data: 164 self.f2_half_cone_90_y[i, j] = val

165 166

167 - def setup_identity_free_rotor(self):

168 """Set up a few identity matrices.""" 169 170 # The order identity matrix for the free rotors. 171 data = [[0, 0, 0.5], 172 [1, 1, 0.5], 173 [3, 3, 0.5], 174 [4, 4, 0.5], 175 [0, 4, 0.5], 176 [4, 0, 0.5], 177 [1, 3, -0.5], 178 [3, 1, -0.5], 179 [8, 8, 1.0]] 180 self.I_order_free_rotor = zeros((9, 9), float64) 181 for i, j, val in data: 182 self.I_order_free_rotor[i, j] = val 183 184 # The disorder identity matrix for the free rotors. 185 data = [[0, 0, 0.25], 186 [1, 1, 0.125], 187 [3, 3, 0.125], 188 [4, 4, 0.25], 189 [0, 4, 0.25], 190 [4, 0, 0.25], 191 [1, 3, -0.125], 192 [3, 1, -0.125], 193 [0, 8, 0.5], 194 [8, 0, 0.5], 195 [4, 8, 0.5], 196 [8, 4, 0.5]] 197 self.I_disorder_free_rotor = zeros((9, 9), float64) 198 for i, j, val in data: 199 self.I_disorder_free_rotor[i, j] = val

200 201

202 - def test_compile_2nd_matrix_free_rotor_point1(self):

203 """Check the operation of the compile_2nd_matrix_free_rotor() function.""" 204 205 # The simulated in frame free rotor 2nd degree frame order matrix (1e6 ensembles). 206 real = array( 207 [[ 0.5001, 0.0001, 0, 0.0001, 0.4999, 0, 0, 0, 0], 208 [ -0.0001, 0.5001, 0, -0.4999, 0.0001, 0, 0, 0, 0], 209 [ 0, 0, 0.0006, 0, 0, -0.0005, 0, 0, 0], 210 [ -0.0001, -0.4999, 0, 0.5001, 0.0001, 0, 0, 0, 0], 211 [ 0.4999, -0.0001, 0, -0.0001, 0.5001, 0, 0, 0, 0], 212 [ 0, 0, 0.0005, 0, 0, 0.0006, 0, 0, 0], 213 [ 0, 0, 0, 0, 0, 0, 0.0006, -0.0005, 0], 214 [ 0, 0, 0, 0, 0, 0, 0.0005, 0.0006, 0], 215 [ 0, 0, 0, 0, 0, 0, 0, 0, 1.0000]]) 216 217 # The Kronecker product of the eigenframe rotation. 218 Rx2_eigen = self.calc_Rx2_eigen_axis(0.0, 0.0) 219 220 # Calculate the matrix. 221 f2 = compile_2nd_matrix_free_rotor(self.f2_temp, Rx2_eigen) 222 223 # Print out. 224 print_frame_order_2nd_degree(real, "real") 225 print_frame_order_2nd_degree(f2, "calculated") 226 print_frame_order_2nd_degree(real-f2, "difference") 227 228 # Check the values. 229 for i in range(9): 230 for j in range(9): 231 print("Element %s, %s; diff %s." % (i, j, f2[i, j] - real[i, j])) 232 self.assert_(abs(f2[i, j] - real[i, j]) < 1e-3)

233 234

235 - def test_compile_2nd_matrix_free_rotor_point2(self):

236 """Check the operation of the compile_2nd_matrix_free_rotor() function.""" 237 238 # The simulated free rotor 2nd degree frame order matrix (1e6 ensembles, axis=[2,1,3]). 239 real = array( 240 [[ 0.3367, -0.0100, -0.0307, -0.0100, 0.3521, -0.0152, -0.0307, -0.0152, 0.3112], 241 [ -0.0104, 0.3520, -0.0152, -0.2908, -0.0559, 0.2602, 0.1989, -0.1685, 0.0664], 242 [ -0.0306, -0.0155, 0.3112, 0.1991, -0.1683, 0.0666, 0.2399, 0.2092, 0.1989], 243 [ -0.0104, -0.2908, 0.1989, 0.3520, -0.0559, -0.1685, -0.0152, 0.2602, 0.0664], 244 [ 0.3520, -0.0563, -0.1684, -0.0563, 0.4362, -0.0841, -0.1684, -0.0841, 0.2118], 245 [ -0.0153, 0.2602, 0.0661, -0.1684, -0.0844, 0.2117, 0.2093, -0.0740, 0.0997], 246 [ -0.0306, 0.1991, 0.2399, -0.0155, -0.1683, 0.2092, 0.3112, 0.0666, 0.1989], 247 [ -0.0153, -0.1684, 0.2093, 0.2602, -0.0844, -0.0740, 0.0661, 0.2117, 0.0997], 248 [ 0.3113, 0.0663, 0.1991, 0.0663, 0.2117, 0.0993, 0.1991, 0.0993, 0.4770]]) 249 250 # The cone axis. 251 r, theta, phi = cartesian_to_spherical([2, 1, 3]) 252 253 # The Kronecker product of the eigenframe rotation. 254 Rx2_eigen = self.calc_Rx2_eigen_axis(theta, phi) 255 256 # Calculate the matrix. 257 f2 = compile_2nd_matrix_free_rotor(self.f2_temp, Rx2_eigen) 258 259 # Print out. 260 print_frame_order_2nd_degree(real, "real") 261 print_frame_order_2nd_degree(f2, "calculated") 262 print_frame_order_2nd_degree(real-f2, "difference") 263 264 # Check the values. 265 for i in range(9): 266 for j in range(9): 267 print("Element %s, %s; diff %s." % (i, j, f2[i, j] - real[i, j])) 268 self.assert_(abs(f2[i, j] - real[i, j]) < 1e-3)

269 270

271 - def test_compile_2nd_matrix_iso_cone_disorder(self):

272 """Check if compile_2nd_matrix_iso_cone() can return the identity matrix for disorder.""" 273 274 # The Kronecker product of the eigenframe rotation. 275 Rx2_eigen = self.calc_Rx2_eigen_axis(0.0, 0.0) 276 277 # Calculate the frame order matrix. 278 f2 = compile_2nd_matrix_iso_cone(self.f2_temp, Rx2_eigen, pi, pi) 279 280 # Print outs. 281 print_frame_order_2nd_degree(self.I_disorder, "Identity for disorder") 282 print_frame_order_2nd_degree(f2, "Compiled frame order") 283 284 # Check the values. 285 for i in range(9): 286 for j in range(9): 287 print("Element %s, %s." % (i, j)) 288 self.assertAlmostEqual(f2[i, j], self.I_disorder[i, j])

289 290

291 - def test_compile_2nd_matrix_iso_cone_half_cone(self):

292 """Check if compile_2nd_matrix_iso_cone() can return the matrix for a half cone.""" 293 294 # The Kronecker product of the eigenframe rotation. 295 Rx2_eigen = self.calc_Rx2_eigen_axis(0.0, 0.0) 296 297 # Calculate the frame order matrix. 298 f2 = compile_2nd_matrix_iso_cone(self.f2_temp, Rx2_eigen, pi/2.0, pi) 299 300 # Print outs. 301 print_frame_order_2nd_degree(self.f2_half_cone, "The half cone frame order matrix") 302 print_frame_order_2nd_degree(f2, "Compiled frame order") 303 304 # Check the values. 305 for i in range(9): 306 for j in range(9): 307 print("Element %s, %s." % (i, j)) 308 self.assertAlmostEqual(f2[i, j], self.f2_half_cone[i, j])

309 310

311 - def test_compile_2nd_matrix_iso_cone_half_cone_90_y(self):

312 """Check if compile_2nd_matrix_iso_cone() can return the matrix for a half cone rotated 90 degrees about y.""" 313 314 # The Kronecker product of the eigenframe rotation. 315 Rx2_eigen = self.calc_Rx2_eigen_axis(pi/2.0, 0.0) 316 317 # Calculate the frame order matrix. 318 f2 = compile_2nd_matrix_iso_cone(self.f2_temp, Rx2_eigen, pi/2.0, pi) 319 320 # Print outs. 321 print_frame_order_2nd_degree(self.f2_half_cone_90_y, "The half cone frame order matrix") 322 print_frame_order_2nd_degree(f2, "Compiled frame order") 323 324 # Check the values. 325 for i in range(9): 326 for j in range(9): 327 print("Element %s, %s." % (i, j)) 328 self.assertAlmostEqual(f2[i, j], self.f2_half_cone_90_y[i, j])

329 330

331 - def test_compile_2nd_matrix_iso_cone_order(self):

332 """Check if compile_2nd_matrix_iso_cone() can return the identity matrix for order.""" 333 334 # The Kronecker product of the eigenframe rotation. 335 Rx2_eigen = self.calc_Rx2_eigen_axis(0.0, 0.0) 336 337 # Calculate the frame order matrix. 338 f2 = compile_2nd_matrix_iso_cone(self.f2_temp, Rx2_eigen, 1e-5, 1e-10) 339 340 # Print outs. 341 print_frame_order_2nd_degree(self.I_order, "Identity for order") 342 print_frame_order_2nd_degree(f2, "Compiled frame order") 343 344 # Check the values. 345 for i in range(9): 346 for j in range(9): 347 print("Element %s, %s." % (i, j)) 348 self.assertAlmostEqual(f2[i, j], self.I_order[i, j])

349 350

351 - def test_compile_2nd_matrix_iso_cone_order2(self):

352 """2nd check if compile_2nd_matrix_iso_cone() can return the identity matrix for order.""" 353 354 # The Kronecker product of the eigenframe rotation. 355 Rx2_eigen = self.calc_Rx2_eigen_axis(0.0, 0.0) 356 357 # Calculate the frame order matrix. 358 f2 = compile_2nd_matrix_iso_cone(self.f2_temp, Rx2_eigen, 0.0, 0.0) 359 360 # Print outs. 361 print_frame_order_2nd_degree(self.I_order, "Identity for order") 362 print_frame_order_2nd_degree(f2, "Compiled frame order") 363 364 # Check the values. 365 for i in range(9): 366 for j in range(9): 367 print("Element %s, %s." % (i, j)) 368 self.assertAlmostEqual(f2[i, j], self.I_order[i, j])

369 370

371 - def test_compile_2nd_matrix_iso_cone_restriction_test(self):

372 """Check if compile_2nd_matrix_iso_cone() can approximate compile_2nd_matrix_iso_cone_free_rotor().""" 373 374 # The Kronecker product of the eigenframe rotation. 375 Rx2_eigen_a = self.calc_Rx2_eigen_axis(0.0, 0.0) 376 Rx2_eigen_b = self.calc_Rx2_eigen_axis(0.0, 1.0) 377 378 # Calculate the frame order matrix. 379 f2a = compile_2nd_matrix_iso_cone(self.f2_temp, Rx2_eigen_a, pi/4.6, pi) 380 f2b = compile_2nd_matrix_iso_cone_free_rotor(self.f2_temp, Rx2_eigen_b, pi/4.6) 381 382 # Print outs. 383 print_frame_order_2nd_degree(f2a, "Isotropic cone frame order") 384 print_frame_order_2nd_degree(f2b, "Free rotor isotropic cone frame order") 385 print_frame_order_2nd_degree(f2b-f2a, "difference") 386 387 # Check the values. 388 for i in range(9): 389 for j in range(9): 390 print("Element %s, %s." % (i, j)) 391 self.assertAlmostEqual(f2a[i, j], f2b[i, j])

392 393

394 - def test_compile_2nd_matrix_iso_cone_restriction_test2(self):

395 """Check if compile_2nd_matrix_iso_cone() can approximate compile_2nd_matrix_iso_cone_torsionless().""" 396 397 # The Kronecker product of the eigenframe rotation. 398 Rx2_eigen_a = self.calc_Rx2_eigen_axis(0.0, 0.0) 399 Rx2_eigen_b = self.calc_Rx2_eigen_axis(0.0, 0.0) 400 401 # Calculate the frame order matrix. 402 f2a = compile_2nd_matrix_iso_cone(self.f2_temp, Rx2_eigen_a, pi/4.6, 0) 403 f2b = compile_2nd_matrix_iso_cone_torsionless(self.f2_temp, Rx2_eigen_b, pi/4.6) 404 405 # Print outs. 406 print_frame_order_2nd_degree(f2a, "Isotropic cone frame order") 407 print_frame_order_2nd_degree(f2b, "Torsionless isotropic cone frame order") 408 print_frame_order_2nd_degree(f2b-f2a, "difference") 409 410 # Check the values. 411 for i in range(9): 412 for j in range(9): 413 print("Element %s, %s." % (i, j)) 414 self.assertAlmostEqual(f2a[i, j], f2b[i, j])

415 416

417 - def test_compile_2nd_matrix_iso_cone_free_rotor_disorder(self):

418 """Check if compile_2nd_matrix_iso_cone_free_rotor() can return the identity matrix for disorder.""" 419 420 # The Kronecker product of the eigenframe rotation. 421 Rx2_eigen = self.calc_Rx2_eigen_axis(0.0, 1.0) 422 423 # Calculate the frame order matrix. 424 f2 = compile_2nd_matrix_iso_cone_free_rotor(self.f2_temp, Rx2_eigen, pi) 425 426 # Cannot differentiate full disorder from the half cone in this model! 427 f2_ref = self.f2_half_cone 428 429 # Print outs. 430 print_frame_order_2nd_degree(self.I_disorder, "Identity for disorder") 431 print_frame_order_2nd_degree(f2, "Compiled frame order") 432 433 # Check the values. 434 for i in range(9): 435 for j in range(9): 436 print("Element %s, %s." % (i, j)) 437 self.assertAlmostEqual(f2[i, j], self.I_disorder[i, j])

438 439

440 - def test_compile_2nd_matrix_iso_cone_free_rotor_half_cone(self):

441 """Check if compile_2nd_matrix_iso_cone_free_rotor() can return the matrix for a half cone.""" 442 443 # The Kronecker product of the eigenframe rotation. 444 Rx2_eigen = self.calc_Rx2_eigen_axis(0.0, 1.0) 445 446 # Calculate the frame order matrix. 447 f2 = compile_2nd_matrix_iso_cone_free_rotor(self.f2_temp, Rx2_eigen, pi/2.0) 448 449 # Print outs. 450 print_frame_order_2nd_degree(self.f2_half_cone, "The half cone frame order matrix") 451 print_frame_order_2nd_degree(f2, "Compiled frame order") 452 453 # Check the values. 454 for i in range(9): 455 for j in range(9): 456 print("Element %s, %s." % (i, j)) 457 self.assertAlmostEqual(f2[i, j], self.f2_half_cone[i, j])

458 459

460 - def test_compile_2nd_matrix_iso_cone_free_rotor_half_cone_90_y(self):

461 """Check if compile_2nd_matrix_iso_cone_free_rotor() can return the matrix for a half cone rotated 90 degrees about y.""" 462 463 # Init. 464 self.z_axis = array([1, 0, 0], float64) 465 466 # The Kronecker product of the eigenframe rotation. 467 Rx2_eigen = self.calc_Rx2_eigen_axis(0.0, 1.0) 468 469 # Calculate the frame order matrix. 470 f2 = compile_2nd_matrix_iso_cone_free_rotor(self.f2_temp, Rx2_eigen, pi/2.0) 471 472 # Print outs. 473 print_frame_order_2nd_degree(self.f2_half_cone_90_y, "The half cone frame order matrix") 474 print_frame_order_2nd_degree(f2, "Compiled frame order") 475 476 # Check the values. 477 for i in range(9): 478 for j in range(9): 479 print("Element %s, %s." % (i, j)) 480 self.assertAlmostEqual(f2[i, j], self.f2_half_cone_90_y[i, j])

481 482

483 - def test_compile_2nd_matrix_iso_cone_free_rotor_order(self):

484 """Check if compile_2nd_matrix_iso_cone_free_rotor() can return the identity matrix for order.""" 485 486 # The Kronecker product of the eigenframe rotation. 487 Rx2_eigen = self.calc_Rx2_eigen_axis(0.0, 1.0) 488 489 # Calculate the frame order matrix. 490 f2 = compile_2nd_matrix_iso_cone_free_rotor(self.f2_temp, Rx2_eigen, 0.0) 491 492 # Print outs. 493 print_frame_order_2nd_degree(self.I_order_free_rotor, "Free rotor identity for order") 494 print_frame_order_2nd_degree(f2, "Compiled frame order") 495 496 # Check the values. 497 for i in range(9): 498 for j in range(9): 499 print("Element %s, %s." % (i, j)) 500 self.assertAlmostEqual(f2[i, j], self.I_order_free_rotor[i, j])

501 502

503 - def test_compile_2nd_matrix_pseudo_ellipse_point1(self):

504 """Check the operation of the compile_2nd_matrix_pseudo_ellipse() function.""" 505 506 # The simulated in frame pseudo-ellipse 2nd degree frame order matrix. 507 real = array( 508 [[ 0.7901, 0, 0, 0, 0.7118, 0, 0, 0, 0.6851], 509 [ 0, 0.0816, 0, -0.0606, 0, 0, 0, 0, 0], 510 [ 0, 0, 0.1282, 0, 0, 0, -0.1224, 0, 0], 511 [ 0, -0.0606, 0, 0.0708, 0, 0, 0, 0, 0], 512 [ 0.7118, 0, 0, 0, 0.6756, 0, 0, 0, 0.6429], 513 [ 0, 0, 0, 0, 0, 0.2536, 0, -0.2421, 0], 514 [ 0, 0, -0.1224, 0, 0, 0, 0.1391, 0, 0], 515 [ 0, 0, 0, 0, 0, -0.2421, 0, 0.2427, 0], 516 [ 0.6851, 0, 0, 0, 0.6429, 0, 0, 0, 0.6182]], float64) 517 transpose_23(real) 518 519 # Init. 520 x = pi/4.0 521 y = 3.0*pi/8.0 522 z = pi/6.0 523 524 # The Kronecker product of the eigenframe rotation. 525 Rx2_eigen = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 526 527 # Calculate the matrix. 528 f2 = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen, x, y, z) 529 530 # Print out. 531 print_frame_order_2nd_degree(real, "real") 532 print_frame_order_2nd_degree(f2, "calculated") 533 print_frame_order_2nd_degree(real-f2, "difference") 534 535 # Check the values. 536 for i in range(9): 537 for j in range(9): 538 print("Element %s, %s; diff %s." % (i, j, f2[i, j] - real[i, j])) 539 self.assert_(abs(f2[i, j] - real[i, j]) < 1e-4)

540 541

542 - def test_compile_2nd_matrix_pseudo_ellipse_point2(self):

543 """Check the operation of the compile_2nd_matrix_pseudo_ellipse() function.""" 544 545 # The simulated in frame pseudo-ellipse 2nd degree frame order matrix (1e6 ensembles). 546 real = array( 547 [[ 0.7379, 0, 0, 0, 0.1338, 0, 0, 0, 0.1284], 548 [ 0, 0.6637, 0, -0.1085, 0, 0, 0, 0, 0], 549 [ 0, 0, 0.6603, 0, 0, 0, -0.1181, 0, 0], 550 [ 0, -0.1085, 0, 0.6637, 0, 0, 0, 0, 0], 551 [ 0.1154, 0, 0, 0, 0.6309, 0, 0, 0, 0.2536], 552 [ 0, 0, 0, 0, 0, 0.6196, 0, -0.2336, 0], 553 [ 0, 0, -0.1181, 0, 0, 0, 0.6603, 0, 0], 554 [ 0, 0, 0, 0, 0, -0.2336, 0, 0.6196, 0], 555 [ 0.1467, 0, 0, 0, 0.2353, 0, 0, 0, 0.6180]], float64) 556 557 # Init. 558 x = pi/4.0 559 y = 3.0*pi/8.0 560 z = 40.0 / 360.0 * 2.0 * pi 561 562 # The Kronecker product of the eigenframe rotation. 563 Rx2_eigen = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 564 565 # Calculate the matrix. 566 f2 = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen, x, y, z) 567 568 # Print out. 569 print_frame_order_2nd_degree(real, "real") 570 print_frame_order_2nd_degree(f2, "calculated") 571 print_frame_order_2nd_degree(real-f2, "difference") 572 573 # Check the values. 574 for i in range(9): 575 for j in range(9): 576 print("Element %s, %s; diff %s." % (i, j, f2[i, j] - real[i, j])) 577 self.assert_(abs(f2[i, j] - real[i, j]) < 1e-3)

578 579

580 - def test_compile_2nd_matrix_pseudo_ellipse_point3(self):

581 """Check the operation of the compile_2nd_matrix_pseudo_ellipse() function.""" 582 583 # The simulated out of frame pseudo-ellipse 2nd degree frame order matrix (1e6 ensembles). 584 real = array( 585 [[ 0.6314, 0.0232, -0.0344, 0.0232, 0.1558, -0.0222, -0.0344, -0.0222, 0.2128], 586 [ 0.0220, 0.6366, 0.0069, -0.1352, 0.0243, -0.0722, 0.0206, -0.0277, -0.0464], 587 [ -0.0332, 0.0097, 0.6137, 0.0222, 0.0668, 0.0173, -0.1967, 0.0489, -0.0336], 588 [ 0.0220, -0.1352, 0.0206, 0.6366, 0.0243, -0.0277, 0.0069, -0.0722, -0.0464], 589 [ 0.1554, 0.0233, 0.0669, 0.0233, 0.6775, 0.0113, 0.0669, 0.0113, 0.1671], 590 [ -0.0222, -0.0738, 0.0188, -0.0286, 0.0113, 0.6310, 0.0507, -0.1502, 0.0109], 591 [ -0.0332, 0.0222, -0.1967, 0.0097, 0.0668, 0.0489, 0.6137, 0.0173, -0.0336], 592 [ -0.0222, -0.0286, 0.0507, -0.0738, 0.0113, -0.1502, 0.0188, 0.6310, 0.0109], 593 [ 0.2132, -0.0465, -0.0324, -0.0465, 0.1667, 0.0110, -0.0324, 0.0110, 0.6201]], float64) 594 595 # Init. 596 x = 60.0 / 360.0 * 2.0 * pi 597 y = 3.0 * pi / 8.0 598 z = pi / 6.0 599 600 # The Kronecker product of the eigenframe rotation. 601 Rx2_eigen = self.calc_Rx2_eigen_full(self.out_of_frame_alpha, self.out_of_frame_beta, self.out_of_frame_gamma) 602 603 # Calculate the matrix. 604 f2 = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen, x, y, z) 605 606 # Print out. 607 print_frame_order_2nd_degree(real, "real") 608 print_frame_order_2nd_degree(f2, "calculated") 609 print_frame_order_2nd_degree(real-f2, "difference") 610 611 # Check the values. 612 for i in range(9): 613 for j in range(9): 614 print("Element %s, %s; diff %s." % (i, j, f2[i, j] - real[i, j])) 615 self.assert_(abs(f2[i, j] - real[i, j]) < 1e-3)

616 617

618 - def test_compile_2nd_matrix_pseudo_ellipse_disorder(self):

619 """Check if compile_2nd_matrix_pseudo_ellipse() can return the identity matrix for disorder.""" 620 621 # The Kronecker product of the eigenframe rotation. 622 Rx2_eigen = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 623 624 # Calculate the frame order matrix. 625 f2 = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen, pi, pi, pi) 626 627 # Print outs. 628 print_frame_order_2nd_degree(self.I_disorder, "Identity for disorder") 629 print_frame_order_2nd_degree(f2, "Compiled frame order") 630 631 # Check the values. 632 for i in range(9): 633 for j in range(9): 634 print("Element %s, %s." % (i, j)) 635 self.assertAlmostEqual(f2[i, j], self.I_disorder[i, j])

636 637

638 - def test_compile_2nd_matrix_pseudo_ellipse_half_cone(self):

639 """Check if compile_2nd_matrix_pseudo_ellipse() can return the matrix for a half cone.""" 640 641 # The Kronecker product of the eigenframe rotation. 642 Rx2_eigen = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 643 644 # Calculate the frame order matrix. 645 f2 = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen, pi/2.0, pi/2.0, pi) 646 647 # Print outs. 648 print_frame_order_2nd_degree(self.f2_half_cone, "The half cone frame order matrix") 649 print_frame_order_2nd_degree(f2, "Compiled frame order") 650 651 # Check the values. 652 for i in range(9): 653 for j in range(9): 654 print("Element %s, %s." % (i, j)) 655 self.assertAlmostEqual(f2[i, j], self.f2_half_cone[i, j])

656 657

658 - def test_compile_2nd_matrix_pseudo_ellipse_half_cone_90_y(self):

659 """Check if compile_2nd_matrix_pseudo_ellipse() can return the matrix for a half cone rotated 90 degrees about y.""" 660 661 # The Kronecker product of the eigenframe rotation. 662 Rx2_eigen = self.calc_Rx2_eigen_full(0.0, pi/2.0, 0.0) 663 664 # Calculate the frame order matrix. 665 f2 = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen, pi/2.0, pi/2.0, pi) 666 667 # Print outs. 668 print_frame_order_2nd_degree(self.f2_half_cone_90_y, "The half cone frame order matrix") 669 print_frame_order_2nd_degree(f2, "Compiled frame order") 670 671 # Check the values. 672 for i in range(9): 673 for j in range(9): 674 print("Element %s, %s." % (i, j)) 675 self.assertAlmostEqual(f2[i, j], self.f2_half_cone_90_y[i, j])

676 677

678 - def test_compile_2nd_matrix_pseudo_ellipse_order(self):

679 """Check if compile_2nd_matrix_pseudo_ellipse() can return the identity matrix for order.""" 680 681 # The Kronecker product of the eigenframe rotation. 682 Rx2_eigen = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 683 684 # Calculate the frame order matrix. 685 f2 = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen, 1e-2, 1e-2, 1e-10) 686 687 # Print outs. 688 print_frame_order_2nd_degree(self.I_order, "Identity for order") 689 print_frame_order_2nd_degree(f2, "Compiled frame order") 690 print_frame_order_2nd_degree(f2-self.I_order, "difference") 691 692 # Check the values. 693 for i in range(9): 694 for j in range(9): 695 print("Element %s, %s." % (i, j)) 696 self.assertAlmostEqual(f2[i, j], self.I_order[i, j], 4)

697 698

699 - def test_compile_2nd_matrix_pseudo_ellipse_restriction_test(self):

700 """Check if compile_2nd_matrix_pseudo_ellipse() can approximate compile_2nd_matrix_pseudo_ellipse_free_rotor().""" 701 702 # The Kronecker product of the eigenframe rotation. 703 Rx2_eigen_a = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 704 Rx2_eigen_b = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 705 706 # Calculate the frame order matrix. 707 f2a = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen_a, pi/1.6, pi/5.8, pi) 708 f2b = compile_2nd_matrix_pseudo_ellipse_free_rotor(self.f2_temp, Rx2_eigen_b, pi/1.6, pi/5.8) 709 710 # Print outs. 711 print_frame_order_2nd_degree(f2a, "Pseudo-ellipse frame order") 712 print_frame_order_2nd_degree(f2b, "Free rotor pseudo-ellipse frame order") 713 print_frame_order_2nd_degree(f2b-f2a, "difference") 714 715 # Check the values. 716 for i in range(9): 717 for j in range(9): 718 print("Element %s, %s." % (i, j)) 719 self.assertAlmostEqual(f2a[i, j], f2b[i, j])

720 721

722 - def test_compile_2nd_matrix_pseudo_ellipse_restriction_test2(self):

723 """Check if compile_2nd_matrix_pseudo_ellipse() can approximate a pi/2 rotated compile_2nd_matrix_pseudo_ellipse_free_rotor().""" 724 725 # The Kronecker product of the eigenframe rotation. 726 Rx2_eigen_a = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 727 Rx2_eigen_b = self.calc_Rx2_eigen_full(pi/2, 0.0, 0.0) 728 729 # Calculate the frame order matrix. 730 f2a = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen_a, pi/1.6, pi/5.8, pi) 731 f2b = compile_2nd_matrix_pseudo_ellipse_free_rotor(self.f2_temp, Rx2_eigen_b, pi/5.8, pi/1.6) 732 733 # Print outs. 734 print_frame_order_2nd_degree(f2a, "Pseudo-ellipse frame order") 735 print_frame_order_2nd_degree(f2b, "pi/2 rotated free rotor pseudo-ellipse frame order") 736 print_frame_order_2nd_degree(f2b-f2a, "difference") 737 738 # Check the values. 739 for i in range(9): 740 for j in range(9): 741 print("Element %s, %s." % (i, j)) 742 self.assertAlmostEqual(f2a[i, j], f2b[i, j])

743 744

745 - def test_compile_2nd_matrix_pseudo_ellipse_restriction_test3(self):

746 """Check if compile_2nd_matrix_pseudo_ellipse() can approximate compile_2nd_matrix_pseudo_ellipse_torsionless().""" 747 748 # The Kronecker product of the eigenframe rotation. 749 Rx2_eigen_a = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 750 Rx2_eigen_b = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 751 752 # Calculate the frame order matrix. 753 f2a = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen_a, pi/2.1, pi/4.6, 0) 754 f2b = compile_2nd_matrix_pseudo_ellipse_torsionless(self.f2_temp, Rx2_eigen_b, pi/2.1, pi/4.6) 755 756 # Print outs. 757 print_frame_order_2nd_degree(f2a, "Pseudo-ellipse frame order") 758 print_frame_order_2nd_degree(f2b, "Torsionless pseudo-ellipse frame order") 759 print_frame_order_2nd_degree(f2b-f2a, "difference") 760 761 # Check the values. 762 for i in range(9): 763 for j in range(9): 764 print("Element %s, %s; diff %s." % (i, j, f2b[i, j] - f2a[i, j])) 765 self.assertAlmostEqual(f2a[i, j], f2b[i, j])

766 767

768 - def test_compile_2nd_matrix_pseudo_ellipse_restriction_test4(self):

769 """Check if compile_2nd_matrix_pseudo_ellipse() can approximate compile_2nd_matrix_iso_cone().""" 770 771 # The Kronecker product of the eigenframe rotation. 772 Rx2_eigen_a = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 773 Rx2_eigen_b = self.calc_Rx2_eigen_axis(0.0, 0.0) 774 775 # Calculate the frame order matrix. 776 f2a = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen_a, pi/4.6, pi/4.6, 0.2) 777 f2b = compile_2nd_matrix_iso_cone(self.f2_temp, Rx2_eigen_b, pi/4.6, 0.2) 778 779 # Print outs. 780 print_frame_order_2nd_degree(f2a, "Pseudo-ellipse frame order") 781 print_frame_order_2nd_degree(f2b, "Isotropic cone frame order") 782 print_frame_order_2nd_degree(f2b-f2a, "difference") 783 784 # Check the values. 785 for i in range(9): 786 for j in range(9): 787 print("Element %s, %s." % (i, j)) 788 self.assertAlmostEqual(f2a[i, j], f2b[i, j])

789 790

791 - def test_compile_2nd_matrix_pseudo_ellipse_restriction_test5(self):

792 """Check if compile_2nd_matrix_pseudo_ellipse() can approximate compile_2nd_matrix_iso_cone_free_rotor().""" 793 794 # The Kronecker product of the eigenframe rotation. 795 Rx2_eigen_a = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 796 Rx2_eigen_b = self.calc_Rx2_eigen_axis(0.0, 1.0) 797 798 # Calculate the frame order matrix. 799 f2a = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen_a, pi/4.6, pi/4.6, pi) 800 f2b = compile_2nd_matrix_iso_cone_free_rotor(self.f2_temp, Rx2_eigen_b, pi/4.6) 801 802 # Print outs. 803 print_frame_order_2nd_degree(f2a, "Pseudo-ellipse frame order") 804 print_frame_order_2nd_degree(f2b, "Free rotor isotropic cone frame order") 805 print_frame_order_2nd_degree(f2b-f2a, "difference") 806 807 # Check the values. 808 for i in range(9): 809 for j in range(9): 810 print("Element %s, %s." % (i, j)) 811 self.assertAlmostEqual(f2a[i, j], f2b[i, j])

812 813

814 - def test_compile_2nd_matrix_pseudo_ellipse_restriction_test6(self):

815 """Check if compile_2nd_matrix_pseudo_ellipse() can approximate compile_2nd_matrix_iso_cone_torsionless().""" 816 817 # The Kronecker product of the eigenframe rotation. 818 Rx2_eigen_a = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 819 Rx2_eigen_b = self.calc_Rx2_eigen_axis(0.0, 0.0) 820 821 # Calculate the frame order matrix. 822 f2a = compile_2nd_matrix_pseudo_ellipse(self.f2_temp, Rx2_eigen_a, pi/8.6, pi/8.6, 0) 823 f2b = compile_2nd_matrix_iso_cone_torsionless(self.f2_temp, Rx2_eigen_b, pi/8.6) 824 825 # Print outs. 826 print_frame_order_2nd_degree(f2a, "Pseudo-ellipse frame order") 827 print_frame_order_2nd_degree(f2b, "Torsionless isotropic cone frame order") 828 print_frame_order_2nd_degree(f2b-f2a, "difference") 829 830 # Check the values. 831 for i in range(9): 832 for j in range(9): 833 print("Element %s, %s." % (i, j)) 834 self.assertAlmostEqual(f2a[i, j], f2b[i, j])

835 836

837 - def test_compile_2nd_matrix_pseudo_ellipse_free_rotor_disorder(self):

838 """Check if compile_2nd_matrix_pseudo_ellipse_free_rotor() can return the identity matrix for disorder.""" 839 840 # The Kronecker product of the eigenframe rotation. 841 Rx2_eigen = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 842 843 # Calculate the frame order matrix. 844 f2 = compile_2nd_matrix_pseudo_ellipse_free_rotor(self.f2_temp, Rx2_eigen, pi, pi) 845 846 # Print outs. 847 print_frame_order_2nd_degree(self.I_disorder, "Identity for disorder") 848 print_frame_order_2nd_degree(f2, "Compiled frame order") 849 850 # Check the values. 851 for i in range(9): 852 for j in range(9): 853 print("Element %s, %s." % (i, j)) 854 self.assertAlmostEqual(f2[i, j], self.I_disorder[i, j])

855 856

857 - def test_compile_2nd_matrix_pseudo_ellipse_free_rotor_half_cone(self):

858 """Check if compile_2nd_matrix_pseudo_ellipse() can return the matrix for a half cone.""" 859 860 # The Kronecker product of the eigenframe rotation. 861 Rx2_eigen = self.calc_Rx2_eigen_full(pi, 0.0, pi) 862 863 # Calculate the frame order matrix (rotated about z by 2pi). 864 f2 = compile_2nd_matrix_pseudo_ellipse_free_rotor(self.f2_temp, Rx2_eigen, pi/2.0, pi/2.0) 865 866 # Print outs. 867 print_frame_order_2nd_degree(self.f2_half_cone, "The half cone frame order matrix") 868 print_frame_order_2nd_degree(f2, "Compiled frame order") 869 870 # Check the values. 871 for i in range(9): 872 for j in range(9): 873 print("Element %s, %s." % (i, j)) 874 self.assertAlmostEqual(f2[i, j], self.f2_half_cone[i, j])

875 876

877 - def test_compile_2nd_matrix_pseudo_ellipse_free_rotor_half_cone_90_y(self):

878 """Check if compile_2nd_matrix_pseudo_ellipse() can return the matrix for a half cone rotated 90 degrees about y.""" 879 880 # The Kronecker product of the eigenframe rotation. 881 Rx2_eigen = self.calc_Rx2_eigen_full(pi, pi/2.0, pi) 882 883 # Calculate the frame order matrix (rotated about z by 2pi). 884 f2 = compile_2nd_matrix_pseudo_ellipse_free_rotor(self.f2_temp, Rx2_eigen, pi/2.0, pi/2.0) 885 886 # Print outs. 887 print_frame_order_2nd_degree(self.f2_half_cone_90_y, "The half cone frame order matrix") 888 print_frame_order_2nd_degree(f2, "Compiled frame order") 889 890 # Check the values. 891 for i in range(9): 892 for j in range(9): 893 print("Element %s, %s." % (i, j)) 894 self.assertAlmostEqual(f2[i, j], self.f2_half_cone_90_y[i, j])

895 896

897 - def fixme_test_compile_2nd_matrix_pseudo_ellipse_free_rotor_order(self):

898 """Check if compile_2nd_matrix_pseudo_ellipse_free_rotor() can return the identity matrix for order.""" 899 900 # The Kronecker product of the eigenframe rotation. 901 Rx2_eigen = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 902 903 # Calculate the frame order matrix. 904 f2 = compile_2nd_matrix_pseudo_ellipse_free_rotor(self.f2_temp, Rx2_eigen, 1e-10, 1e-10) 905 906 # Print outs. 907 print_frame_order_2nd_degree(self.I_order_free_rotor, "Free rotor identity for order") 908 print_frame_order_2nd_degree(f2, "Compiled frame order") 909 910 # Check the values. 911 for i in range(9): 912 for j in range(9): 913 print("Element %s, %s." % (i, j)) 914 self.assertAlmostEqual(f2[i, j], self.I_order_free_rotor[i, j])

915 916

917 - def test_compile_2nd_matrix_pseudo_ellipse_free_rotor_point1(self):

918 """Check the operation of the compile_2nd_matrix_pseudo_ellipse_free_rotor() function.""" 919 920 # The simulated out of frame free rotor pseudo-ellipse 2nd degree frame order matrix (1e6 ensembles). 921 real = array( 922 [[ 0.3428, -0.0193, 0.0389, -0.0193, 0.3137, -0.0194, 0.0389, -0.0194, 0.3435], 923 [ -0.0225, 0.2313, 0.0034, -0.1413, 0.0449, 0.2309, -0.1830, -0.1412, -0.0224], 924 [ 0.0417, 0.0091, 0.2142, -0.1767, -0.0838, 0.0092, 0.1211, -0.1770, 0.0421], 925 [ -0.0225, -0.1413, -0.1830, 0.2313, 0.0449, -0.1412, 0.0034, 0.2309, -0.0224], 926 [ 0.3124, 0.0418, -0.0840, 0.0418, 0.3758, 0.0418, -0.0840, 0.0418, 0.3118], 927 [ -0.0193, 0.2251, 0.0151, -0.1476, 0.0389, 0.2251, -0.1706, -0.1468, -0.0196], 928 [ 0.0417, -0.1767, 0.1211, 0.0091, -0.0838, -0.1770, 0.2142, 0.0092, 0.0421], 929 [ -0.0193, -0.1476, -0.1706, 0.2251, 0.0389, -0.1468, 0.0151, 0.2251, -0.0196], 930 [ 0.3448, -0.0225, 0.0450, -0.0225, 0.3104, -0.0224, 0.0450, -0.0224, 0.3447]], float64) 931 932 # Init. 933 x = pi/4.0 934 y = 50.0 / 360.0 * 2.0 * pi 935 936 # The Kronecker product of the eigenframe rotation. 937 Rx2_eigen = self.calc_Rx2_eigen_full(self.out_of_frame_alpha, self.out_of_frame_beta, self.out_of_frame_gamma) 938 939 # Calculate the matrix. 940 f2 = compile_2nd_matrix_pseudo_ellipse_free_rotor(self.f2_temp, Rx2_eigen, x, y) 941 942 # Print out. 943 print_frame_order_2nd_degree(real, "real") 944 print_frame_order_2nd_degree(f2, "calculated") 945 print_frame_order_2nd_degree(real-f2, "difference") 946 947 # Check the values. 948 for i in range(9): 949 for j in range(9): 950 print("Element %s, %s." % (i, j)) 951 self.assert_(abs(f2[i, j] - real[i, j]) < 1e-3)

952 953

954 - def test_compile_2nd_matrix_pseudo_ellipse_free_rotor_point2(self):

955 """Check the operation of the compile_2nd_matrix_pseudo_ellipse_free_rotor() function.""" 956 957 # The simulated out of frame free rotor pseudo-ellipse 2nd degree frame order matrix (1e6 ensembles). 958 real = array( 959 [[ 0.3251, 0.0163, -0.0324, 0.0163, 0.3493, 0.0164, -0.0324, 0.0164, 0.3256], 960 [ -0.0248, 0.1481, -0.0480, -0.0500, 0.0492, 0.1475, -0.1472, -0.0500, -0.0244], 961 [ 0.0079, 0.0328, 0.0572, -0.0661, -0.0163, 0.0331, 0.0074, -0.0662, 0.0084], 962 [ -0.0248, -0.0500, -0.1472, 0.1481, 0.0492, -0.0500, -0.0480, 0.1475, -0.0244], 963 [ 0.3289, 0.0081, -0.0167, 0.0081, 0.3426, 0.0080, -0.0167, 0.0080, 0.3285], 964 [ 0.0163, 0.0669, 0.1139, -0.1322, -0.0324, 0.0662, 0.0157, -0.1307, 0.0161], 965 [ 0.0079, -0.0661, 0.0074, 0.0328, -0.0163, -0.0662, 0.0572, 0.0331, 0.0084], 966 [ 0.0163, -0.1322, 0.0157, 0.0669, -0.0324, -0.1307, 0.1139, 0.0662, 0.0161], 967 [ 0.3459, -0.0245, 0.0491, -0.0245, 0.3081, -0.0244, 0.0491, -0.0244, 0.3460]], float64) 968 969 # Init. 970 x = pi / 4.0 971 y = 150.0 / 360.0 * 2.0 * pi 972 973 # The Kronecker product of the eigenframe rotation. 974 Rx2_eigen = self.calc_Rx2_eigen_full(self.out_of_frame_alpha, self.out_of_frame_beta, self.out_of_frame_gamma) 975 976 # Calculate the matrix. 977 f2 = compile_2nd_matrix_pseudo_ellipse_free_rotor(self.f2_temp, Rx2_eigen, x, y) 978 979 # Print out. 980 print_frame_order_2nd_degree(real, "real") 981 print_frame_order_2nd_degree(f2, "calculated") 982 print_frame_order_2nd_degree(real-f2, "difference") 983 984 # Check the values. 985 for i in range(9): 986 for j in range(9): 987 print("Element %s, %s." % (i, j)) 988 self.assert_(abs(f2[i, j] - real[i, j]) < 1e-2)

989 990

991 - def test_compile_2nd_matrix_pseudo_ellipse_free_rotor_restriction_test(self):

992 """Check if compile_2nd_matrix_pseudo_ellipse_free_rotor() can approximate compile_2nd_matrix_iso_cone_free_rotor().""" 993 994 # The Kronecker product of the eigenframe rotation. 995 Rx2_eigen_a = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 996 Rx2_eigen_b = self.calc_Rx2_eigen_axis(0.0, 1.0) 997 998 # Calculate the frame order matrix. 999 f2a = compile_2nd_matrix_pseudo_ellipse_free_rotor(self.f2_temp, Rx2_eigen_a, pi/4.6, pi/4.6) 1000 f2b = compile_2nd_matrix_iso_cone_free_rotor(self.f2_temp, Rx2_eigen_b, pi/4.6) 1001 1002 # Print outs. 1003 print_frame_order_2nd_degree(f2a, "Free rotor pseudo-ellipse frame order") 1004 print_frame_order_2nd_degree(f2b, "Free rotor isotropic cone frame order") 1005 print_frame_order_2nd_degree(f2b-f2a, "difference") 1006 1007 # Check the values. 1008 for i in range(9): 1009 for j in range(9): 1010 print("Element %s, %s." % (i, j)) 1011 self.assertAlmostEqual(f2a[i, j], f2b[i, j])

1012 1013

1014 - def test_compile_2nd_matrix_pseudo_ellipse_torsionless_restriction_test(self):

1015 """Check if compile_2nd_matrix_pseudo_ellipse_torsionless() can approximate compile_2nd_matrix_iso_cone_torsionless().""" 1016 1017 # The Kronecker product of the eigenframe rotation. 1018 Rx2_eigen_a = self.calc_Rx2_eigen_full(0.0, 0.0, 0.0) 1019 Rx2_eigen_b = self.calc_Rx2_eigen_axis(0.0, 0.0) 1020 1021 # Calculate the frame order matrix. 1022 f2a = compile_2nd_matrix_pseudo_ellipse_torsionless(self.f2_temp, Rx2_eigen_a, pi/4.6, pi/4.6) 1023 f2b = compile_2nd_matrix_iso_cone_torsionless(self.f2_temp, Rx2_eigen_b, pi/4.6) 1024 1025 # Print outs. 1026 print_frame_order_2nd_degree(f2a, "Torsionless pseudo-ellipse frame order") 1027 print_frame_order_2nd_degree(f2b, "Torsionless isotropic cone frame order") 1028 print_frame_order_2nd_degree(f2b-f2a, "difference") 1029 1030 # Check the values. 1031 for i in range(9): 1032 for j in range(9): 1033 print("Element %s, %s." % (i, j)) 1034 self.assertAlmostEqual(f2a[i, j], f2b[i, j])

1035 1036

1037 - def test_compile_2nd_matrix_rotor_point1(self):

1038 """Check the operation of the compile_2nd_matrix_rotor() function.""" 1039 1040 # The simulated in frame rotor 2nd degree frame order matrix (1e7 ensembles). 1041 real = array( 1042 [[ 7.06775425e-01, 1.36710179e-04, 0.00000000e+00, 1.36710179e-04, 2.93224575e-01, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00], 1043 [ -1.36710179e-04, 7.06775425e-01, 0.00000000e+00, -2.93224575e-01, 1.36710179e-04, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00], 1044 [ 0.00000000e+00, 0.00000000e+00, 8.27014112e-01, 0.00000000e+00, 0.00000000e+00, 2.19539417e-04, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00], 1045 [ -1.36710179e-04, -2.93224575e-01, 0.00000000e+00, 7.06775425e-01, 1.36710179e-04, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00], 1046 [ 2.93224575e-01, -1.36710179e-04, 0.00000000e+00, -1.36710179e-04, 7.06775425e-01, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00], 1047 [ 0.00000000e+00, 0.00000000e+00, -2.19539417e-04, 0.00000000e+00, 0.00000000e+00, 8.27014112e-01, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00], 1048 [ 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 8.27014112e-01, 2.19539417e-04, 0.00000000e+00], 1049 [ 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, -2.19539417e-04, 8.27014112e-01, 0.00000000e+00], 1050 [ 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 0.00000000e+00, 1.00000000e+00]]) 1051 1052 # Init. 1053 sigma_max = 60.0 / 180.0 * pi 1054 1055 # The Kronecker product of the eigenframe rotation. 1056 Rx2_eigen = self.calc_Rx2_eigen_axis(0.0, 0.0) 1057 1058 # Calculate the matrix. 1059 f2 = compile_2nd_matrix_rotor(self.f2_temp, Rx2_eigen, sigma_max) 1060 1061 # Print out. 1062 print_frame_order_2nd_degree(real, "real") 1063 print_frame_order_2nd_degree(f2, "calculated") 1064 print_frame_order_2nd_degree(real-f2, "difference") 1065 1066 # Check the values. 1067 for i in range(9): 1068 for j in range(9): 1069 print("Element %s, %s." % (i, j)) 1070 self.assert_(abs(f2[i, j] - real[i, j]) < 1e-3)

1071 1072

1073 - def test_reduce_alignment_tensor_order(self):

1074 """Test the alignment tensor reduction for the order identity matrix.""" 1075 1076 # The tensors. 1077 A = array([1, 2, 3, 4, 5], float64) 1078 red = zeros(5, float64) 1079 1080 # Reduce. 1081 reduce_alignment_tensor(self.I_order, A, red) 1082 1083 # Check. 1084 for i in range(5): 1085 self.assertEqual(A[i], red[i])

1086 1087

1088 - def test_reduce_alignment_tensor_disorder(self):

1089 """Test the alignment tensor reduction for the order identity matrix.""" 1090 1091 # The tensors. 1092 A = array([1, 2, 3, 4, 5], float64) 1093 red = zeros(5, float64) 1094 1095 # Reduce. 1096 reduce_alignment_tensor(self.I_disorder, A, red) 1097 1098 # Check. 1099 for i in range(5): 1100 self.assertEqual(red[i], 0.0)

1101 1102

1103 - def test_reduce_alignment_tensor_half_cone(self):

1104 """Test the alignment tensor reduction for the order identity matrix.""" 1105 1106 # The tensors. 1107 A = array([1, 2, 3, 4, 5], float64) 1108 red = zeros(5, float64) 1109 1110 # Reduce. 1111 reduce_alignment_tensor(self.f2_half_cone, A, red) 1112 1113 # Check. 1114 for i in range(5): 1115 self.assertEqual(red[i], 0.0)

Source Code for Module test_suite.unit_tests._lib._frame_order.test_matrix_ops