/* Sort eigenvalues in ascending order */
#define SWAPPER(i) \
- if (fabs(dd[i+1]) < fabs(dd[i])) { \
+ if (std::abs(dd[i+1]) < std::abs(dd[i])) { \
temp = dd[i]; \
for (j = 0; (j < NDIM); j++) { tvec[j] = ev[j][i]; } \
dd[i] = dd[i+1]; \
{
if (bQ)
{
- m0 = fabs(atom[ii].q);
+ m0 = std::abs(atom[ii].q);
}
else
{