Teaching numerical methods : a first experience. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai.

Transcription

1 Teaching numerical methods : a first experience Ronojoy Adhikari The Institute of Mathematical Sciences Chennai.

2 Numerical Methods : the first class 14 students, spread across second and third years of undergraduate study. Physics, Mathematics and Computer Science. Physics : knew why numerical methods are necessary, moderate programming skills. Mathematics : very little idea about numerical methods, very poor programming skills. Computer Science : poor idea about numerical methods, very good programming skills.

3 Purpose The purpose of computing is insight, not numbers. Numerical methods can give insight into problems we cannot otherwise solve. R. W. Hamming

4 No purpose! The purpose of computing numbers is not yet in sight.

5 Motivation : physical models Computational Fluid Dynamics Molecular Dynamics Nonlinear Dynamics CFD MD NLD

6 Motivation : mathematics of algorithms A higher order polynomial interpolation does not always guarantee a smaller error. Indeed, it may so happen that the error increases with the degree of the interpolating polynomial. This is called the Runge phenomenon. Example: interpolating the Lorentzian f(x) = x 2

7 Motivation : mathematics of algorithms x k = cos( kπ N ) 1.5 equispaced points 1.5 Chebyshev points !0.5 max error = !0.5 max error = !1!1! !1!1! Black magic?

8 Motivation : algorithmic complexity Toeplitz (constant entries along diagonals) and circulant (reflecting pbc). f 1.. f N = h f 1.. f N pbc: f N+1 = f 1, f 0 = f N N 1 2 Other boundary conditions also possible.

9 Motivation : algorithmic complexity f 1... f N = h f 1.. f N With quartic interpolation the matrix is pentadiagonal: The bandwidth increases with increasing order of interpolation. What is the best way to compute a spectral differentiation matrix?

10 Mathematical methods What do we learn in a mathematical methods course? Calculus Linear Algebra Orthogonal Functions Differential Equations Partial Differential Equations...

11 Numerical methods What do we learn in a numerical methods course? Calculus Linear Algebra Orthogonal Functions Differential Equations Partial Differential Equations...

12 Interpolation : the central idea data (x k, y k ) k = 1... N interpolant I(x; {x k }, {y k }) derivative quadrature di(x; {x k }, {y k }) dx x = x i b a I(x; {x k }, {y k })dx

13 Numerical analysis and computer programming What does numerical analysis have to do with C, C++..? Nothing really... Numerical analysis can be done in almost any computer language. Historically, Fortran was the language of choice. Modern codes are written in C, C++, Java, Python, Ruby... The course was language agnostic!

14 Open resources Wikipedia is an excellent starting point for topics.

15 Open resource : programs

16 Open resource : libraries software/gsl/

17 Open resource : environments

18 Open resource : books and websites

19 Course website for out-of-class interaction. Undergraduate course given to the students of the Chennai Mathematical Institute. Mix of physics, mathematics and computer science students, spread across different years. Standard topics were covered. Choice was given to code in C/C++ or Matlab/ Octave, but teaching was language agnostic. GSL was recommended to C/C++ users. Format was 1.5 hours of classroom lecture per week and 1.5 hours of lab. Students knew at least one computer language. Weekly assignments and a final project were used for evaluation.

20 Outlook Undergraduate courses in numerical methods extremely necessary. Must focus on basics first, the building blocks of numerical analysis, before proceeding towards packages or molecular dynamics and Monte Carlo. Do not reinvent the wheel : use libraries! GSL is a very good general purpose library, suited to work in physics. (eg : has a Monte Carlo module). For moderate work (roughly 10,0000 degrees of freedom) there is an enormous saving in human time in using environments like Matlab/Octave, while the saving in computer time is not more than half. Strongly encourage use of environments. Try to make numerical methods fun, not a drudgery of using the VI editor, messing with compiler flags, linking...

21 Thanks to the organisers for inviting me to speak. Thank you for your attention.