Merge pull request #2440 from chebfun/fix-chebpoly

Fix chebpoly
chebfun · Apr 17, 2023 · d9c5e7b · d9c5e7b
2 parents b839628 + e20e785
commit d9c5e7b
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 0 deletions.
diff --git a/chebpoly.m b/chebpoly.m
@@ -52,6 +52,8 @@
 
 % Construct the Chebyshev coefficients:
 N = max(n) + 1;
+% c = eye(N); % c = c(:,n+1); % Replaced by the below. See Issue #2438.
+% Code below does the same as above without forming the NxN identity matrix.
 M = numel(n);
 c = zeros(N, M);
 c(n+1,:) = eye(M);

diff --git a/tests/chebfun/test_chebpoly.m b/tests/chebfun/test_chebpoly.m
@@ -72,4 +72,13 @@
 
 warning(warnState);
 
+% Test on very high degree polynomials (#2438)
+try 
+    chebpoly(1e5);
+    pass(9) = true;
+catch
+    pass(9) = false;
+end
+
+
 end