From 6710fb13c2f9cc4cf72e60f353597bdbb141578e Mon Sep 17 00:00:00 2001 From: L. Tibaldo Date: Fri, 15 Jun 2018 18:13:23 +0200 Subject: [PATCH] if model = 0 error on residual for SUBDIV and SUBDIVSQRT set to 0 --- examples/show_residuals.py | 10 ++++++++-- 1 file changed, 8 insertions(+), 2 deletions(-) diff --git a/examples/show_residuals.py b/examples/show_residuals.py index 414ebaa..395acae 100755 --- a/examples/show_residuals.py +++ b/examples/show_residuals.py @@ -114,9 +114,15 @@ def fill_cmr(row, counts, model, resid, e_counts, e_resid, c_counts, c_model, if algorithm == 'SUB': e_resid.append(err) elif algorithm == 'SUBDIV': - e_resid.append(err / c_model.real(row)) + if c_model.real(row) > 0.: + e_resid.append(err / c_model.real(row)) + else: + e_resid.append(0.) elif algorithm == 'SUBDIVSQRT': - e_resid.append(err / math.sqrt(c_model.real(row))) + if c_model.real(row) > 0.: + e_resid.append(err / math.sqrt(c_model.real(row))) + else: + e_resid.append(0.) elif algorithm == 'SIGNIFICANCE': e_resid.append(1.0) else: -- libgit2 0.21.4