def lammps_to_lammps(name, lmppath='', lmpbin='lmp_openmpi', energy=True, verbose=False): """Test lammps to lammps conversion """ lmp_in = os.path.join('Inputs/Lammps/', name, 'data.lmp') if not os.path.isfile(lmp_in): raise Exception("File not found: {0}!".format(lmp_in)) lmp_out = os.path.join('Outputs/LammpsToLammps/', name, 'data.lmp') # calc input energies if energy: e_in = lammps_energies(name, 'in', lmppath, lmpbin, verbose) # where the magic happens Driver.initSystem(name) Driver.load(lmp_in) Driver.write(lmp_out) # calc output energies if energy: e_out = lammps_energies(name, 'LtoL', lmppath, lmpbin) if energy: return combine_energy_results(e_in, e_out)
def lammps_to_gromacs(name, gropath='', grosuff='', lmppath='', lmpbin='lmp_openmpi', energy=True, clean=True): """Test lammps to gromacs conversion """ lmp_in = os.path.join('Inputs/Lammps/', name, 'data.lmp') if not os.path.isfile(lmp_in): raise Exception("File not found: {0}!".format(lmp_in)) top_out = os.path.join('Outputs/LammpsToGromacs/', name, 'topol.top') gro_out = os.path.join('Outputs/LammpsToGromacs/', name, 'conf.gro') # calc input energies if energy: e_in = lammps_energies(name, 'in', lmppath, lmpbin) # where the magic happens Driver.initSystem(name) Driver.load(lmp_in) Driver.write(top_out, gro_out) # calc output energies if energy: e_out = gromacs_energies(name, top_out, gro_out, 'LtoG', gropath, grosuff) # delete gromacs backup files if clean: import glob filelist = glob.glob("Inputs/Gromacs/{name}/#*#".format(name=name)) filelist += glob.glob("Outputs/GromacsToGromacs/{name}/#*#".format(name=name)) for f in filelist: os.remove(f) if energy: return combine_energy_results(e_in, e_out)