Undergrad Question about weights using Chebyshev polynomials as quadrature

Click For Summary
SUMMARY

The discussion centers on the application of the Clenshaw–Curtis quadrature scheme using Chebyshev polynomials for optimizing integral calculations. It specifically addresses the calculation of weights for collocation points when the number of points, N, is even, highlighting that for N=6, the weight w3 is defined only once despite appearing to be calculated twice. Participants clarify that this is a deliberate simplification in the article to avoid redundancy in notation. The resolution confirms that the weights are correctly defined and that the perceived duplication is a misunderstanding of the notation.

PREREQUISITES
  • Understanding of Clenshaw–Curtis quadrature scheme
  • Familiarity with Chebyshev polynomials
  • Knowledge of numerical integration techniques
  • Basic comprehension of mathematical notation and weight calculations
NEXT STEPS
  • Study the Clenshaw–Curtis quadrature method in detail
  • Explore the properties and applications of Chebyshev polynomials
  • Learn about numerical integration techniques and their optimizations
  • Investigate the implications of odd and even collocation points in quadrature methods
USEFUL FOR

Mathematicians, numerical analysts, and optimization researchers interested in advanced numerical integration methods and their applications in computational mathematics.

confused_engineer
Messages
34
Reaction score
1
Hello everyone.

I am studying this article since I am interested in optimization. The article makes use of Clenshaw–Curtis quadrature scheme to discretize the integral part of the cost function to a finite sum using Chebyshev polynomials.

The article differentiates between the case of odd and even number of collocation points. In equation 27 and 28 (fourth page), the case of N even is discussed. If N is even, then the N+1 collocation points, including 0, form a vector of odd length.

Then weights are calculated as ws=wN-s=... for s=1, 2, ..., N/2. Meanwhile, w0 and wN are calculated on a different way.

This in turn means that one of the elements is calculated twice, as shown in the following example:

N=6; N+1=7; N/2=3; s=1, 2, 3; N-s=6-1=5, 6-2=4, 6-3=3.

As you can see, the fourth element of the vector, number 3, appears two times. I find this weird. Can someone please tell me if I am understanding the article wrong or if this is intended to happen?

Thanks for reading.
Regards.
 
Physics news on Phys.org
I do not see where the problem is. For ##N=6## we have ##w_0=w_6=1/35## and then formulas for ##w_1=w_5##, ##w_2=w_4##, and ##w_3##. Since ##s=3= N-s=6-3## and the formula for ##w_3## only depends on the given number ##N## and ##s=3##, where is it defined twice?
 
  • Like
Likes confused_engineer
fresh_42 said:
I do not see where the problem is. For ##N=6## we have ##w_0=w_6=1/35## and then formulas for ##w_1=w_5##, ##w_2=w_4##, and ##w_3##. Since ##s=3= N-s=6-3## and the formula for ##w_3## only depends on the given number ##N## and ##s=3##, where is it defined twice?
First of all, thanks for answering my question.

I am confused because the article tallks about ωs and ωN-s and I find extrange that following the paper one arrives to ω33=4/N*... for this particular example.

Therefore, I understand from your answer that ω33=4/N*... is the weight of the fourth element of the vector. I thought that since I had ω33, the middle element might have twice the weight.

Sorry if I haven't expressed myself clearly.
Thanks again for your answer.
 
No, it was only a way to avoid writing an extra line for ##\omega_{N/2}##. ##N/2 = N-N/2## for even ##N## so it is only defined once.
 
  • Like
Likes confused_engineer
fresh_42 said:
No, it was only a way to avoid writing an extra line for ##\omega_{N/2}##. ##N/2 = N-N/2## for even ##N## so it is only defined once.
Then is all clear. Thank you very much.
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 18 ·
Replies
18
Views
3K
Replies
2
Views
2K
  • · Replies 4 ·
Replies
4
Views
1K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
6
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 4 ·
Replies
4
Views
3K