mailr17276 - /branches/frame_order_testing/specific_fns/frame_order.py


Others Months | Index by Date | Thread Index
>>   [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Header


Content

Posted by edward on July 17, 2012 - 17:02:
Author: bugman
Date: Tue Jul 17 17:02:23 2012
New Revision: 17276

URL: http://svn.gna.org/viewcvs/relax?rev=17276&view=rev
Log:
Updates to the frame order module for the interatomic data design.

This includes import cleanups.


Modified:
    branches/frame_order_testing/specific_fns/frame_order.py

Modified: branches/frame_order_testing/specific_fns/frame_order.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/frame_order_testing/specific_fns/frame_order.py?rev=17276&r1=17275&r2=17276&view=diff
==============================================================================
--- branches/frame_order_testing/specific_fns/frame_order.py (original)
+++ branches/frame_order_testing/specific_fns/frame_order.py Tue Jul 17 
17:02:23 2012
@@ -30,26 +30,27 @@
 from numpy import arccos, array, dot, eye, float64, identity, ones, 
transpose, zeros
 from numpy.linalg import inv
 from re import search
-from string import upper
 from warnings import warn
 
 # relax module imports.
-from api_base import API_base
-from api_common import API_common
+import arg_check
 from float import isNaN, isInf
-from generic_fns import align_tensor, pipes
+from generic_fns import pipes
 from generic_fns.angles import wrap_angles
+from generic_fns.interatomic import interatomic_loop
 from generic_fns.mol_res_spin import return_spin, spin_loop
+from generic_fns.structure import geometric
 from generic_fns.structure.cones import Iso_cone, Pseudo_elliptic
-from generic_fns.structure.geometric import create_cone_pdb, 
generate_vector_dist, generate_vector_residues
 from generic_fns.structure.internal import Internal
 from maths_fns import frame_order, order_parameters
 from maths_fns.coord_transform import spherical_to_cartesian
 from maths_fns.rotation_matrix import euler_to_R_zyz, two_vect_to_R
 from physical_constants import dipolar_constant, g1H, 
return_gyromagnetic_ratio
-from relax_errors import RelaxError, RelaxInfError, RelaxModelError, 
RelaxNaNError, RelaxNoModelError, RelaxNoValueError, RelaxSpinTypeError
+from relax_errors import RelaxError, RelaxInfError, RelaxNaNError, 
RelaxNoModelError, RelaxNoValueError, RelaxSpinTypeError
 from relax_io import open_write_file
-from relax_warnings import RelaxWarning, RelaxDeselectWarning
+from relax_warnings import RelaxWarning
+from specific_fns.api_base import API_base
+from specific_fns.api_common import API_common
 
 
 class Frame_order(API_base, API_common):
@@ -238,8 +239,8 @@
         list = []
 
         # RDC search.
-        for spin in spin_loop():
-            if hasattr(spin, 'rdc'):
+        for interatom in interatomic_loop():
+            if hasattr(interatom, 'rdc'):
                 list.append('rdc')
                 break
 
@@ -409,11 +410,11 @@
 
             # Generate the axis vectors.
             print("\nGenerating the axis vectors.")
-            res_num = generate_vector_residues(mol=mol, vector=axis_pos, 
atom_name='z-ax', res_name_vect='AXE', sim_vectors=axis_sim_pos, res_num=2, 
origin=cdp.pivot, scale=size)
+            res_num = geometric.generate_vector_residues(mol=mol, 
vector=axis_pos, atom_name='z-ax', res_name_vect='AXE', 
sim_vectors=axis_sim_pos, res_num=2, origin=cdp.pivot, scale=size)
 
             # The negative.
             if neg_cone:
-                res_num = generate_vector_residues(mol=mol_neg, 
vector=axis_neg, atom_name='z-ax', res_name_vect='AXE', 
sim_vectors=axis_sim_neg, res_num=2, origin=cdp.pivot, scale=size)
+                res_num = geometric.generate_vector_residues(mol=mol_neg, 
vector=axis_neg, atom_name='z-ax', res_name_vect='AXE', 
sim_vectors=axis_sim_neg, res_num=2, origin=cdp.pivot, scale=size)
 
         # The full axis system.
         else:
@@ -458,9 +459,9 @@
                     axis_sim_neg = axes_sim_neg[:, :, j]
 
                 # The vectors.
-                res_num = generate_vector_residues(mol=mol, 
vector=axes_pos[:, j], atom_name='%s-ax'%label[j], res_name_vect='AXE', 
sim_vectors=axis_sim_pos, res_num=2, origin=cdp.pivot, scale=size)
+                res_num = geometric.generate_vector_residues(mol=mol, 
vector=axes_pos[:, j], atom_name='%s-ax'%label[j], res_name_vect='AXE', 
sim_vectors=axis_sim_pos, res_num=2, origin=cdp.pivot, scale=size)
                 if neg_cone:
-                    res_num = generate_vector_residues(mol=mol_neg, 
vector=axes_neg[:, j], atom_name='%s-ax'%label[j], res_name_vect='AXE', 
sim_vectors=axis_sim_neg, res_num=2, origin=cdp.pivot, scale=size)
+                    res_num = 
geometric.generate_vector_residues(mol=mol_neg, vector=axes_neg[:, j], 
atom_name='%s-ax'%label[j], res_name_vect='AXE', sim_vectors=axis_sim_neg, 
res_num=2, origin=cdp.pivot, scale=size)
 
 
         # The cone object.
@@ -495,11 +496,11 @@
                 cone = Iso_cone(cone_theta)
 
             # Create the positive and negative cones.
-            create_cone_pdb(mol=mol, cone=cone, start_res=mol.res_num[-1]+1, 
apex=cdp.pivot, R=R_pos, inc=inc, distribution='regular')
+            geometric.create_cone_pdb(mol=mol, cone=cone, 
start_res=mol.res_num[-1]+1, apex=cdp.pivot, R=R_pos, inc=inc, 
distribution='regular')
 
             # The negative.
             if neg_cone:
-                create_cone_pdb(mol=mol_neg, cone=cone, 
start_res=mol_neg.res_num[-1]+1, apex=cdp.pivot, R=R_neg, inc=inc, 
distribution='regular')
+                geometric.create_cone_pdb(mol=mol_neg, cone=cone, 
start_res=mol_neg.res_num[-1]+1, apex=cdp.pivot, R=R_neg, inc=inc, 
distribution='regular')
 
 
         # Create the PDB file.




Related Messages


Powered by MHonArc, Updated Wed Jul 18 11:00:02 2012