PC5215, ``Numerical Recipes with Applications'' Semester I, Aug-Nov 2022

INSTRUCTOR: Prof. Wang Jian-Sheng

CLASS SCHEDULE: Tue/Thu 10:00-12:00 at S11-02-05. Final exam Thursday 1pm, 24 Nov 2022.

Lab briefing will be conducted at Computer Lab 2 (S13-M-08) on ??day ?-?pm and only on the weeks of new lab assignments (i.e. weeks 2, 4, 9, 11), as posted below.

TEXTBOOK: ``Numerical Recipes in C (1992), the Art of Scientific Computing,'' W. H. Press, etc.

MODULE DESCRIPTION: Covers computational techniques for the solution of problems arising in physics and engineering, with an emphasis on molecular simulation and modelling. Topics will be from the text, "Numerical Recipes", Press et al, supplemented with example problems in materials and condensed matter physics. The textbook is in C, but this course will use Python for programming. This course insures that graduate students intend to do research in computational physics will have sufficient background in computational methods and programming experience. Assessment: 4 labs + assignments 35%, midterm 15%, final 50%.

COURSE OUTLINE (lecture slides)

week 0, 1-5 Aug, no class

week 1, 11 Aug, introduction, Python (9 Aug is National day)

week 2 , 16,17 Aug, Chap 1, Python continued, error, accuracy, etc

week 3, 23,25 Aug, Chap 2, LU decomposition lab 1

week 4, 30 Aug, 1 Sep, Chap 3, interpolation and extrapolation

week 5, 6,8 Sep, Chap 4, integration, gaussian quadrature, lab 2

week 6, 13,15 Sep, Chap 7, Monte Carlo integration

week 7, recess week, no class

week 8, 27,29 sep, introduction to learning algorithm/neural network (midterm test at Thursday lecture time)

week 9, 4,6 Oct, Chap 10, optimization, conjugate gradient, lab 3

week 10, 11,13 Oct, Chap 12, fast fourier transform

week 11, 18,20 Oct, Chap 15, modeling of data (least square), the final lab, lab 4

week 12, 25,27 Oct, Chap 16, ordinary differential equations, molecular dynamics

week 13, 1,3 Nov, ODE, continued

week 14, 8,10 Nov, boundary value problem, revision (last week)


Review Problems (with partial answers)

2008, 2010, 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018, 2019, 2020. final exam and answers.

Some sample codes: stability.py, ludcmp.py, polint.py, trapzd.py, mnist.py, conjugategradient.py, nfit.dat, testieee.c, main-eigen.c, eigen.c, sort.c, fft.c.

Some interesting sites related to the course: IEEE 754 (PDF) floating point arithmetic, LAPACK, an easy to read article on conjugate gradient method, fast Fourier transform FFTW, symplectic integrator.

To install a python interpreter on Windows, go to https://www.python.org, and follow the download link for Windows or your preferred OS.