Skip to content
Snippets Groups Projects
Commit b19c5324 authored by James Graham's avatar James Graham
Browse files

style: resolve a couple of linting complaints

parent 0dbce7d5
No related branches found
No related tags found
No related merge requests found
...@@ -45,7 +45,7 @@ class PyCGTOOL: ...@@ -45,7 +45,7 @@ class PyCGTOOL:
frame_number=0) frame_number=0)
if self.out_frame.n_frames > 1: if self.out_frame.n_frames > 1:
self.train_backmapper(self.in_frame, self.out_frame) self.train_backmapper()
self.bondset = None self.bondset = None
if self.config.bondset: if self.config.bondset:
...@@ -93,7 +93,8 @@ class PyCGTOOL: ...@@ -93,7 +93,8 @@ class PyCGTOOL:
self.config.out_dir) self.config.out_dir)
logger.info('Finished writing bond measurements to file') logger.info('Finished writing bond measurements to file')
def get_coords(self, frame: Frame, resname: str) -> np.ndarray: @staticmethod
def get_coords(frame: Frame, resname: str) -> np.ndarray:
return np.concatenate([ return np.concatenate([
frame._trajectory.atom_slice( frame._trajectory.atom_slice(
[atom.index for atom in residue.atoms]).xyz [atom.index for atom in residue.atoms]).xyz
...@@ -101,15 +102,16 @@ class PyCGTOOL: ...@@ -101,15 +102,16 @@ class PyCGTOOL:
if residue.name == resname if residue.name == resname
]) ])
def train_backmapper(self, aa_frame: Frame, cg_frame: Frame): def train_backmapper(self):
# resname = 'POPC' # resname = 'POPC'
# aa_coords = get_coords(aa_frame, resname) # aa_coords = get_coords(aa_frame, resname)
# cg_coords = get_coords(cg_frame, resname) # cg_coords = get_coords(cg_frame, resname)
sel = 'resid 0'
cg_subset_traj = cg_frame._trajectory.atom_slice( aa_subset_traj = self.in_frame._trajectory.atom_slice(
cg_frame._trajectory.topology.select('resid 0')) self.in_frame._trajectory.topology.select(sel))
aa_subset_traj = aa_frame._trajectory.atom_slice( cg_subset_traj = self.out_frame._trajectory.atom_slice(
aa_frame._trajectory.topology.select('resid 0')) self.out_frame._trajectory.topology.select(sel))
cg_subset_traj.save('cg_test.gro') cg_subset_traj.save('cg_test.gro')
aa_subset_traj.save('aa_test.gro') aa_subset_traj.save('aa_test.gro')
...@@ -143,13 +145,12 @@ class BooleanAction(argparse.Action): ...@@ -143,13 +145,12 @@ class BooleanAction(argparse.Action):
def parse_arguments(arg_list): def parse_arguments(arg_list):
# yapf: disable
parser = argparse.ArgumentParser( parser = argparse.ArgumentParser(
formatter_class=argparse.ArgumentDefaultsHelpFormatter, formatter_class=argparse.ArgumentDefaultsHelpFormatter,
description= description="Generate coarse-grained molecular dynamics models from atomistic trajectories."
"Generate coarse-grained molecular dynamics models from atomistic trajectories."
) )
# yapf: disable
# Input files # Input files
input_files = parser.add_argument_group("input files") input_files = parser.add_argument_group("input files")
...@@ -298,9 +299,8 @@ def main(): ...@@ -298,9 +299,8 @@ def main():
pycgtool = PyCGTOOL(args) pycgtool = PyCGTOOL(args)
elapsed_time = time.time() - start_time elapsed_time = time.time() - start_time
logger.info( logger.info('Processed %d frames in %.2f s',
f'Processed {pycgtool.out_frame.n_frames} frames in {elapsed_time:.2f}s' pycgtool.out_frame.n_frames, elapsed_time)
)
logger.info('Finished processing - goodbye!') logger.info('Finished processing - goodbye!')
except Exception as exc: except Exception as exc:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment