UW AMath High Performance Scientific Computing
AMath 483/583 Class Notes
Spring Quarter, 2013

Table Of Contents

Previous topic

Lab 18: Thursday May 29, 2014

Next topic

Lab 20: Thursday June 5, 2014

This Page

Lab 19: Tuesday June 3, 2014

Some things to try

  • You might want to make a copy of this notebook before you start playing with it.
  • Experiment with different initial conditions for the heat equation.
  • Create an animation (in the notebook) of the numerical solution to the heat equation along with the true solution.
  • The Crank-Nicolson method is second order accurate: the error should go to zero like \({\cal O}(\Delta t^2 + \Delta x^2)\) as the grid is refined. So increasing both n (the number of spatial points) and nsteps (the number of time steps) by a factor of 2 should reduce the error by a factor of 4. Test this out.
  • Compute or look up the Fourier sine series for some interesting function and try this as initial conditions for the heat equation. Compare the true solution with the numerical solution (where the “true solution” might be estimated by adding up a large but finite number of terms in the Fourier series).
  • Try using Sympy to compute the coefficients in the Fourier sine series.

There is quiz for Lab 19