Problem with Chebyshev basis?
See original GitHub issueHi Mikael,
I have started to implement my Rayleigh-Benard solver and wanted to play a bit with the Poisson and Biharmonic solvers first. To do so, I downloaded the scripts available in the demo
folder and I have some weird things going on:
- For both the Poisson and Biharmonic scipts (2D ones), if I set the basis used in the non-homogeneous direction to be Chebyshev, the right-hand-side vector
f_hat = inner(v, fj, output_array=f_hat)
all its entries turn out be to equal to0 + 0j
. - If instead I set the basis to be Legendre polynomials, then everything works out correctly.
For information, I have installed shenfun using conda.
Issue Analytics
- State:
- Created 5 years ago
- Comments:9 (4 by maintainers)
Top Results From Across the Web
Chebyshev polynomials - Wikipedia
The Chebyshev polynomials Tn are polynomials with the largest possible leading coefficient whose absolute value on the interval [−1, 1] is bounded by...
Read more >Asymptotic coefficients and errors for ... - Springer Link
Despite having nearly identical error norms, the error in the Chebyshev basis is concentrated in boundary layers near both endpoints, ...
Read more >Numerical evaluation of polynomials in Chebyshev basis
Evaluating polynomials of arbitrarily large degree in a Chebyshev basis is practical, and provably numerically stable, using a barycentric ...
Read more >Chebyshev Polynomials
They are easy to compute and to convert to and from a power series form. These properties together produce an approximating polynomial which...
Read more >Chapter 8. Chebyshev spectral methods - People
The answer to this question tells us something about the ability of various kinds of spectral methods to resolve data. POINTS PER WAVELENGTH....
Read more >Top Related Medium Post
No results found
Top Related StackOverflow Question
No results found
Troubleshoot Live Code
Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start FreeTop Related Reddit Thread
No results found
Top Related Hackernoon Post
No results found
Top Related Tweet
No results found
Top Related Dev.to Post
No results found
Top Related Hashnode Post
No results found
Top GitHub Comments
Awesome! It now works like a charm. Thanks a lot.
I will get back to you whenever I am done with my 2D Rayleigh-Benard solver. Once it works (and has been validated), I can eventually send a pull request to add it to the demos if you want.
Should be ready now. Build version 7. Just realised, though, that python=2.7 version has not been uploaded to Anaconda cloud for linux after the move to Circle-Ci 2.0 some months ago.