@@ -35,7 +35,8 @@ def select_levels(self, data_container):
3535 that will be used in the entropy calculations.
3636
3737 Args:
38- arg_DataContainer: MDAnalysis universe object containing the system of interest
38+ arg_DataContainer: MDAnalysis universe object containing the system of
39+ interest
3940
4041 Returns:
4142 number_molecules (int): Number of molecules in the system.
@@ -84,7 +85,8 @@ def get_matrices(
8485 level (str): 'polymer', 'residue', or 'united_atom'.
8586 number_frames (int): Number of frames being processed.
8687 highest_level (bool): Whether this is the top (largest bead size) level.
87- force_matrix, torque_matrix (np.ndarray or None): Accumulated matrices to add to.
88+ force_matrix, torque_matrix (np.ndarray or None): Accumulated matrices to add
89+ to.
8890
8991 Returns:
9092 force_matrix (np.ndarray): Accumulated force covariance matrix.
@@ -126,12 +128,8 @@ def get_matrices(
126128
127129 for i in range (number_beads ):
128130 for j in range (i , number_beads ):
129- f_sub = self .create_submatrix (
130- weighted_forces [i ], weighted_forces [j ]
131- )
132- t_sub = self .create_submatrix (
133- weighted_torques [i ], weighted_torques [j ]
134- )
131+ f_sub = self .create_submatrix (weighted_forces [i ], weighted_forces [j ])
132+ t_sub = self .create_submatrix (weighted_torques [i ], weighted_torques [j ])
135133 force_submatrix [i ][j ] = f_sub
136134 force_submatrix [j ][i ] = f_sub .T
137135 torque_submatrix [i ][j ] = t_sub
@@ -330,7 +328,7 @@ def get_beads(self, data_container, level):
330328 atom_group = "resindex " + str (residue )
331329 list_of_beads .append (data_container .select_atoms (atom_group ))
332330
333- #NOTE this could cause problems for hydrogen or helium molecules
331+ # NOTE this could cause problems for hydrogen or helium molecules
334332 if level == "united_atom" :
335333 list_of_beads = []
336334 heavy_atoms = data_container .select_atoms ("not name H*" )
@@ -570,11 +568,12 @@ def get_weighted_forces(
570568 bead : The part of the molecule to be considered.
571569 trans_axes (np.ndarray): The axes relative to which the forces are located.
572570 highest_level (bool): Is this the largest level of the length scale hierarchy
573- force_partitioning (float): Factor to adjust force contributions to avoid over
574- counting correlated forces, default is 0.5.
571+ force_partitioning (float): Factor to adjust force contributions to avoid
572+ over counting correlated forces, default is 0.5.
575573
576574 Returns:
577- weighted_force (np.ndarray): The mass-weighted sum of the forces in the bead.
575+ weighted_force (np.ndarray): The mass-weighted sum of the forces in the
576+ bead.
578577 """
579578
580579 forces_trans = np .zeros ((3 ,))
0 commit comments