sourcestats: add upper bound for skew
This commit is contained in:
parent
11f7cc0507
commit
2e311d1766
1 changed files with 3 additions and 4 deletions
|
@ -47,8 +47,9 @@
|
||||||
2000ppm, which would be pretty bad */
|
2000ppm, which would be pretty bad */
|
||||||
#define WORST_CASE_FREQ_BOUND (2000.0/1.0e6)
|
#define WORST_CASE_FREQ_BOUND (2000.0/1.0e6)
|
||||||
|
|
||||||
/* The minimum allowed skew */
|
/* The minimum and maximum assumed skew */
|
||||||
#define MIN_SKEW 1.0e-12
|
#define MIN_SKEW 1.0e-12
|
||||||
|
#define MAX_SKEW 1.0e+02
|
||||||
|
|
||||||
/* The minimum assumed std dev for weighting */
|
/* The minimum assumed std dev for weighting */
|
||||||
#define MIN_WEIGHT_SD 1.0e-9
|
#define MIN_WEIGHT_SD 1.0e-9
|
||||||
|
@ -541,9 +542,7 @@ SST_DoNewRegression(SST_Stats inst)
|
||||||
inst->std_dev = sqrt(est_var);
|
inst->std_dev = sqrt(est_var);
|
||||||
inst->nruns = nruns;
|
inst->nruns = nruns;
|
||||||
|
|
||||||
if (inst->skew < MIN_SKEW)
|
inst->skew = CLAMP(MIN_SKEW, inst->skew, MAX_SKEW);
|
||||||
inst->skew = MIN_SKEW;
|
|
||||||
|
|
||||||
stress = fabs(old_freq - inst->estimated_frequency) / old_skew;
|
stress = fabs(old_freq - inst->estimated_frequency) / old_skew;
|
||||||
|
|
||||||
DEBUG_LOG(LOGF_SourceStats, "off=%e freq=%e skew=%e n=%d bs=%d runs=%d asym=%f arun=%d",
|
DEBUG_LOG(LOGF_SourceStats, "off=%e freq=%e skew=%e n=%d bs=%d runs=%d asym=%f arun=%d",
|
||||||
|
|
Loading…
Reference in a new issue