ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

do now write out per-type pair cutoff with kspace enabled

Browse Source
This commit is contained in:
Axel Kohlmeyer
2023-01-06 12:12:13 -05:00
parent 909bbcfdbd
commit ce1e997de0
1 changed files with 16 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
src/LEPTON/pair_lepton.cpp
View File

@ -388,8 +388,13 @@ void PairLepton::read_restart_settings(FILE *fp)
void PairLepton::write_data(FILE *fp)
{
for (int i = 1; i <= atom->ntypes; i++)
fprintf(fp, "%d %s %g\n", i, expressions[type2expression[i][i]].c_str(), cut[i][i]);
if (pppmflag || ewaldflag || msmflag || dispersionflag || tip4pflag) {
for (int i = 1; i <= atom->ntypes; i++)
fprintf(fp, "%d %s\n", i, expressions[type2expression[i][i]].c_str());
} else {
for (int i = 1; i <= atom->ntypes; i++)
fprintf(fp, "%d %s %g\n", i, expressions[type2expression[i][i]].c_str(), cut[i][i]);
}
}
/* ----------------------------------------------------------------------
@ -398,9 +403,15 @@ void PairLepton::write_data(FILE *fp)
void PairLepton::write_data_all(FILE *fp)
{
for (int i = 1; i <= atom->ntypes; i++)
for (int j = i; j <= atom->ntypes; j++)
fprintf(fp, "%d %d %s %g\n", i, j, expressions[type2expression[i][j]].c_str(), cut[i][j]);
if (pppmflag || ewaldflag || msmflag || dispersionflag || tip4pflag) {
for (int i = 1; i <= atom->ntypes; i++)
for (int j = i; j <= atom->ntypes; j++)
fprintf(fp, "%d %d %s\n", i, j, expressions[type2expression[i][j]].c_str());
} else {
for (int i = 1; i <= atom->ntypes; i++)
for (int j = i; j <= atom->ntypes; j++)
fprintf(fp, "%d %d %s %g\n", i, j, expressions[type2expression[i][j]].c_str(), cut[i][j]);
}
}
/* ---------------------------------------------------------------------- */