From 2d24a24c4c7d70a2244da8aed6b18443cbe2a879 Mon Sep 17 00:00:00 2001 From: Daniel Peter Date: Wed, 30 Jun 2021 15:49:24 +0200 Subject: [PATCH] avoids zero division for zero traces in corrcoef --- utils/compare_seismogram_correlations.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/utils/compare_seismogram_correlations.py b/utils/compare_seismogram_correlations.py index e93525413..6e79d81d9 100755 --- a/utils/compare_seismogram_correlations.py +++ b/utils/compare_seismogram_correlations.py @@ -247,11 +247,11 @@ def plot_correlations(out_dir,ref_dir): err = norm(ref-syn) #debug - #print('norm syn = %e norm ref = %e' % (norm(syn),fac_norm)) + #print('syn_file: %s\nref_file: %s\nnorm syn = %e\nnorm ref = %e' % (syn_file,ref_file,norm(syn),norm(ref))) # correlation test # total length - if fac_norm > 0.0: + if norm(ref) > 0.0 and norm(syn) > 0.0: corr_mat = np.corrcoef(ref, syn) else: if norm(ref-syn) > 0.0: