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

use math.h function without std:: prefix

Browse Source
This commit is contained in:
Axel Kohlmeyer
2025-06-28 02:22:25 -04:00
parent c28869124b
commit 11c4a8b49c
2 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/EXTRA-PAIR/pair_coul_ctip.cpp
View File

@ -309,7 +309,7 @@ void PairCoulCTIP::init_style()
for (int elt2 = 0; elt2 < nelements; elt2++) {
shield[elt1][elt2] = sqrt(params[elt1].gamma * params[elt2].gamma);
shieldcu[elt1][elt2] = shield[elt1][elt2] * shield[elt1][elt2] * shield[elt1][elt2];
reffc[elt1][elt2] = std::cbrt(cut_coulcu + 1.0 / shieldcu[elt1][elt2]);
reffc[elt1][elt2] = cbrt(cut_coulcu + 1.0 / shieldcu[elt1][elt2]);
reffcsq[elt1][elt2] = reffc[elt1][elt2] * reffc[elt1][elt2];
reffc4[elt1][elt2] = reffcsq[elt1][elt2] * reffcsq[elt1][elt2];
reffc7[elt1][elt2] = reffc4[elt1][elt2] * reffcsq[elt1][elt2] * reffc[elt1][elt2];

2
src/QEQ/fix_qeq_ctip.cpp
View File

@ -122,7 +122,7 @@ void FixQEqCTIP::init()
(t_cut * (A1 + t_cut * (A2 + t_cut * (A3 + t_cut * (A4 + t_cut * A5)))) * erfcd_cut) / r;
for (int elt1 = 0; elt1 < ntypes; elt1++) {
reff[elt1] = std::cbrt(rsq * r + 1.0 / (gamma[elt1 + 1] * gamma[elt1 + 1] * gamma[elt1 + 1]));
reff[elt1] = cbrt(rsq * r + 1.0 / (gamma[elt1 + 1] * gamma[elt1 + 1] * gamma[elt1 + 1]));
reffsq[elt1] = reff[elt1] * reff[elt1];
reff4[elt1] = reffsq[elt1] * reffsq[elt1];
reff7[elt1] = reff4[elt1] * reffsq[elt1] * reff[elt1];