Trees | Indices | Help |
|
---|
|
1 ############################################################################### 2 # # 3 # Copyright (C) 2009-2013 Edward d'Auvergne # 4 # # 5 # This file is part of the program relax (http://www.nmr-relax.com). # 6 # # 7 # This program is free software: you can redistribute it and/or modify # 8 # it under the terms of the GNU General Public License as published by # 9 # the Free Software Foundation, either version 3 of the License, or # 10 # (at your option) any later version. # 11 # # 12 # This program is distributed in the hope that it will be useful, # 13 # but WITHOUT ANY WARRANTY; without even the implied warranty of # 14 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the # 15 # GNU General Public License for more details. # 16 # # 17 # You should have received a copy of the GNU General Public License # 18 # along with this program. If not, see <http://www.gnu.org/licenses/>. # 19 # # 20 ############################################################################### 21 22 # Module docstring. 23 """Module for the handling of Frame Order.""" 24 25 # relax module imports. 26 from lib.frame_order.matrix_ops import rotate_daeg 27 2830 """Generate the rotated 2nd degree Frame Order matrix for the free rotor model. 31 32 The rotor axis is assumed to be parallel to the z-axis in the eigenframe. 33 34 35 @param matrix: The Frame Order matrix, 2nd degree to be populated. 36 @type matrix: numpy 9D, rank-2 array 37 @param Rx2_eigen: The Kronecker product of the eigenframe rotation matrix with itself. 38 @type Rx2_eigen: numpy 9D, rank-2 array 39 """ 40 41 # Zeros. 42 for i in range(9): 43 for j in range(9): 44 matrix[i, j] = 0.0 45 46 # Diagonal. 47 matrix[0, 0] = 0.5 48 matrix[1, 1] = 0.5 49 matrix[3, 3] = 0.5 50 matrix[4, 4] = 0.5 51 matrix[8, 8] = 1 52 53 # Off diagonal set 1. 54 matrix[0, 4] = matrix[4, 0] = 0.5 55 56 # Off diagonal set 2. 57 matrix[1, 3] = matrix[3, 1] = -0.5 58 59 # Rotate and return the frame order matrix. 60 return rotate_daeg(matrix, Rx2_eigen)61
Trees | Indices | Help |
|
---|
Generated by Epydoc 3.0.1 on Thu Jul 3 13:39:02 2014 | http://epydoc.sourceforge.net |