Added symmetrization of cutoff in init_one() and fixed a print

This commit is contained in:
maitanemuba
2024-04-04 18:18:10 +02:00
committed by Axel Kohlmeyer
parent e8294aa207
commit c114938867

View File

@ -300,6 +300,7 @@ double PairYLZ::init_one(int i, int j)
zeta[j][i] = zeta[i][j];
mu[j][i] = mu[i][j];
beta[j][i] = beta[i][j];
cut[j][i] = cut[i][j];
return cut[i][j];
}
@ -409,7 +410,7 @@ void PairYLZ::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 %g %g %g %g %g %g\n", i, j, epsilon[i][i], sigma[i][i], cut[i][j],
fprintf(fp, "%d %d %g %g %g %g %g %g\n", i, j, epsilon[i][j], sigma[i][j], cut[i][j],
zeta[i][j], mu[i][j], beta[i][j]);
}