Esempio n. 1
0
 def make_output(self):
     '''
         Dump Yaff parameters, Yaff system, plot energy contributions along
         perturbation trajectories and dump perturbation trajectories to XYZ
         files.
     '''
     if self.settings.fn_yaff is not None:
         dump_yaff(self.valence, self.settings.fn_yaff)
     if self.settings.fn_charmm22_prm is not None:
         dump_charmm22_prm(self.valence, self.settings.fn_charmm22_prm)
     if self.settings.fn_charmm22_psf is not None:
         dump_charmm22_psf(self.system, self.valence, self.settings.fn_charmm22_psf)
     if self.settings.fn_sys is not None:
         self.system.to_file(self.settings.fn_sys)
     if self.settings.plot_traj is not None and self.settings.plot_traj.lower() in ['Ehc3', 'final', 'all']:
         self.plot_trajectories(do_valence=True, suffix='_Ehc3')
     if self.settings.xyz_traj:
         self.write_trajectories()
Esempio n. 2
0
 def make_output(self):
     '''
         Dump Yaff parameters, Yaff system, plot energy contributions along
         perturbation trajectories and dump perturbation trajectories to XYZ
         files.
     '''
     fn_yaff = self.kwargs.get('fn_yaff', None)
     if fn_yaff is None:
         fn_yaff = 'pars_cov%s.txt' % (self.kwargs.get('suffix', ''))
     self.valence.dump_yaff(fn_yaff)
     fn_charmm22_prm = self.kwargs.get('fn_charmm22_prm')
     if fn_charmm22_prm is not None:
         dump_charmm22_prm(self.valence, fn_charmm22_prm)
     fn_charmm22_psf = self.kwargs.get('fn_charmm22_psf')
     if fn_charmm22_psf is not None:
         dump_charmm22_psf(self.system, self.valence, fn_charmm22_psf)
     fn_sys = self.kwargs.get('fn_sys', None)
     if fn_sys is None:
         fn_sys = 'system%s.chk' % (self.kwargs.get('suffix', ''))
     self.system.to_file(fn_sys)
     self.plot_trajectories(do_valence=True)