LMU/musrfit
5
0
Fork 0
Code Issues 6 Pull Requests Actions Packages Projects Releases 4 Activity

add err==0.0 handling for rounding.

Browse Source
This commit is contained in:
suter_a
2025-07-03 10:22:11 +02:00
parent 781a2e7dbd
commit 25bcdd25e5
1 changed files with 7 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
src/classes/PFitter.cpp
View File

@ -2943,9 +2943,13 @@ PDoubleVector PFitter::ParamRound(const PDoubleVector &par, const PDoubleVector
int exp;
double dval;
for (unsigned int i=0; i<par.size(); i++) {
exp = round(log10(fabs(err[i])))-2;
dval = round(par[i]*pow(10.0, -exp))/pow(10.0, -exp);
par_r[i] = dval;
if (err[i] != 0.0) {
exp = round(log10(fabs(err[i])))-2;
dval = round(par[i]*pow(10.0, -exp))/pow(10.0, -exp);
par_r[i] = dval;
} else {
par_r[i] = par[i];
}
}
return par_r;