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

reduce outlier margin from 5 sigma to 4 sigma

Browse Source
This commit is contained in:
Axel Kohlmeyer
2024-08-24 11:40:15 -04:00
parent 1b3652583b
commit d98a3d61da
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
tools/lammps-gui/chartviewer.cpp
View File

@ -424,7 +424,7 @@ void ChartViewer::add_data(int step, double data)
if (last_step < step) { if (last_step < step) {
last_step = step; last_step = step;
// do not add data that deviates by more than 5 sigma from the average // do not add data that deviates by more than 4 sigma from the average
// over the last 5 to 20 data items. this is a hack to work around // over the last 5 to 20 data items. this is a hack to work around
// getting corrupted data from lammps_get_last_thermo() // getting corrupted data from lammps_get_last_thermo()
const auto &points = series->points(); const auto &points = series->points();
@ -442,7 +442,7 @@ void ChartViewer::add_data(int step, double data)
const double num = count - first; const double num = count - first;
const double avg = ysum / num; const double avg = ysum / num;
const double avgsq = ysumsq / num; const double avgsq = ysumsq / num;
if (fabs(data - avg) > (5.0 * sqrt(avgsq - (avg * avg)))) return; if (fabs(data - avg) > (4.0 * sqrt(avgsq - (avg * avg)))) return;
} }
series->append(step, data); series->append(step, data);