mx/Jungfraujoch
1
1
Fork 0
Code Issues Pull Requests 1 Actions Packages Projects Releases 53 Wiki Activity

Regression: Better guard for not-finite values

Browse Source
This commit is contained in:
leonarski_f
2026-04-14 15:10:33 +02:00
parent 4a852b4d6b
commit f6a0b784d5
1 changed files with 6 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files
image_analysis/bragg_integration/Regression.h
+6 -4
View File
@@ -33,7 +33,7 @@ RegressionResult regression(std::vector<float> &x, std::vector<float> &y, const
int64_t nelem = 0;
for (int i = 0; i < N; i++) {
if (std::isfinite(y[i])) {
if (std::isfinite(y[i]) && std::isfinite(x[i])) {
x_sum += x[i];
y_sum += y[i];
x2_sum += x[i] * x[i];
@@ -58,9 +58,11 @@ RegressionResult regression(std::vector<float> &x, std::vector<float> &y, const
fptype ss_reg = 0.0;
for (int i = 0; i < nelem; i++) {
fptype pred = slope * x[i] + intercept;
ss_tot += (y[i] - y_mean) * (y[i] - y_mean);
ss_reg += (pred - y_mean) * (pred - y_mean);
if (std::isfinite(y[i]) && std::isfinite(x[i])) {
fptype pred = slope * x[i] + intercept;
ss_tot += (y[i] - y_mean) * (y[i] - y_mean);
ss_reg += (pred - y_mean) * (pred - y_mean);
}
}
fptype r_square = 0.0;