change prz to remove zero weights to avoid spurious poles

chebfun · Apr 15, 2023 · 9ed8763 · 9ed8763
1 parent 37d1c85
commit 9ed8763
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/prz.m b/prz.m
@@ -10,6 +10,8 @@
 % Copyright 2018 by The University of Oxford and The Chebfun Developers.
 % See http://www.chebfun.org/ for Chebfun information.
 
+ix = find(wj==0); % indices w/ zero weights; to be removed to avoid spurious poles
+wj(ix) = []; fj(ix) = []; zj(ix) = [];
 m = length(wj);
 
 % Compute poles via generalized eigenvalue problem: