We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent d12454d commit 814a225Copy full SHA for 814a225
source/lib/src/NNPInter.cc
@@ -1700,11 +1700,11 @@ compute_relative_std_f (vector<VALUETYPE> &std,
1700
for (unsigned ii = 0; ii < nloc; ++ii){
1701
const VALUETYPE * tmp_avg = &(avg[ii*3]);
1702
VALUETYPE vdiff[3];
1703
- vdiff[0] = tmp_avg[0];
+ vdiff[0] = tmp_avg[0];
1704
vdiff[1] = tmp_avg[1];
1705
vdiff[2] = tmp_avg[2];
1706
- VALUETYPE f_norm = MathUtilities::dot(vdiff, vdiff);
1707
- // relative std = std/(abs(f)+eps)
+ VALUETYPE f_norm = sqrt(MathUtilities::dot(vdiff, vdiff));
+ // relative std = std/(abs(f)+eps)
1708
std[ii] /= f_norm + eps;
1709
}
1710
0 commit comments