Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 9 April 15, 2011
Today:
Computer architecture
Cache considerations
Monday:
Optimizing Fortran codes
Debugging Fortran
Read: Class notes and references
R.J. LeVeque, University of Washington
AMath 483/583, Lecture 9, April 15, 201
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 23 May 18, 2011
Today:
MPI in subroutines
Comments on Homework 6
Python plotting
Friday:
Grady on GPUs
Read: Class notes and references
New MPI examples.
Send me your info if you want totalview
Make sure Python plotting works
R.J
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 25 May 23, 2011
Today:
Python concepts and objects
Data types, lists, tuples
Modules
NumPy arrays
Next:
More about NumPy and SciPy
Python exception handling
Python functions
Iterators, booleans
Read: Class notes and reference
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 26 May 25, 2011
Today:
NumPy arrays
Python functions
Iterators, booleans
Next:
Python exception handling
Python plus Fortran
Read: Class notes and references
R.J. LeVeque, University of Washington
AMath 483/583, Lecture 26, May
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 28 June 1, 2011
Today:
Python plus Fortran
Comments on quadtests.py for project
Linear vs. loglog plots
Visualization
Friday:
Animation: plots to movies
Binary I/O
Parallel IPython
Reproducible research
Course evaluations (
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 27 May 27, 2011
Today:
Python exception handling
Python plus Fortran: f2py
Next week:
More Python plus Fortran
Visualization
Parallel IPython
Read: Class notes and references
R.J. LeVeque, University of Washington
AMath 483/583,
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 22 May 16, 2011
Today:
GaussSeidel and SOR iterative methods
Totalview debugging
Next:
GPUs, Python
Read: Class notes and references
R.J. LeVeque, University of Washington
AMath 483/583, Lecture 22, May 16, 2011
Jacobi iteration
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 21 May 13, 2011
Today:
OpenMP and MPI versions of Jacobi iteration
GaussSeidel and SOR iterative methods
Next week:
More MPI
Debugging and totalview
GPU computing
Read: Class notes and references
$CLASSHG/codes/openmp/jacobi1.f
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 18 May 6, 2011
Today:
MPI concepts
Communicators, broadcast, reduce
Next week:
MPI send and receive
Iterative methods
Read: Class notes and references
$CLASSHG/codes/mpi
R.J. LeVeque, University of Washington
AMath 483/583, Lectu
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 16 May 2, 2011
Today:
Fine grain vs. coarse grain parallelism
Manually splitting do loops among threads
Wednesday:
Adaptive quadrature, recursive functions
Start MPI?
Read: Class notes and references
R.J. LeVeque, University of W
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 17 May 4, 2011
Today:
Adaptive quadrature, recursive functions
Load balancing with OpenMP
nested forking
Friday:
MPI
Read: Class notes and references
$CLASSHG/codes/adaptive_quadrature
R.J. LeVeque, University of Washington
AMath
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 19 May 9, 2011
Today:
MPI send and receive
Heat equation and discretization
Wednesday:
Iterative methods
Read: Class notes and references
R.J. LeVeque, University of Washington
AMath 483/583, Lecture 19, May 9, 2011
MPI Simple exa
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 24 May 20, 2011
Today:
The Graphical Processing Unit (GPU)
GPU Programming
Todays lecture developed and presented by Grady Lemoine
References:
Andreas Kloeckners High Performance Scientic Computing
course at NYU:
http:/cs.nyu.edu/c
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 15 April 29, 2011
Today:
OpenMP
Fine grain vs. coarse grain parallelism
Next week:
Iterative methods for linear systems
Adaptive quadrature
Start MPI
Read: Class notes and references
R.J. LeVeque, University of Washington
AMath
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 13 April 25, 2011
Today:
Amdahls law
Speed up, strong and weak scaling
OpenMP
Monday:
More OpenMP
Read: Class notes and references
Chapter 2 and Section 6.3 of T. Rauber and G. Rnger,
Parallel Programming For Multicore and Cluste
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 8 April 13, 2011
Today:
Makeles
Friday:
Computer architecture
Cache considerations
Optimizing Fortran codes
Read: Class notes and references
R.J. LeVeque, University of Washington
AMath 483/583, Lecture 8, April 13, 2011
Splittin
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 10 April 18, 2011
Today:
Cache considerations
Optimizing Fortran codes
Debugging Fortran
Wednesday:
Software packages
LAPACK and BLAS
Read: Class notes and references
There are several new sections!
R.J. LeVeque, University of W
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 11 April 20, 2011
Today:
Debugging Fortran
Software packages
zeroin for nding zeros of a function
LAPACK and BLAS
Friday:
Parallel computing concepts
Read: Class notes and references
There are several new sections!
R.J. LeVeque,
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 14 April 27, 2011
Today:
OpenMP
Friday:
More OpenMP
Read: Class notes and references
Homework 4 is posted, due next Thursday.
R.J. LeVeque, University of Washington
AMath 483/583, Lecture 14, April 27, 2011
Memory stack
Note: Paral
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 12 April 22, 2011
Today:
LAPACK and BLAS
Parallel computing concepts
Monday:
OpenMP
Read: Class notes and references
R.J. LeVeque, University of Washington
AMath 483/583, Lecture 12, April 22, 2011
The BLAS
Basic Linear Algebra Su
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 7 April 11, 2011
Today:
Array storage in Fortran
Passing arrays to subroutines
Fortran modules
Multile Fortran codes
Wednesday:
Makele
Read: Class notes and references
There are several new Fortran sections.
R.J. LeVeque, Unive
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 6 April 8, 2011
Today:
Computer arithmetic
Fortran subroutines and functions
Monday:
Fortran array storage
Fortran modules
Multile Fortran codes
Read: Class notes and references.
R.J. LeVeque, University of Washington
AMath 483
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 2 March 30, 2011
Today:
Unix
Version control main ideas
Clientserver version control, e.g., CVS, svn
Distributed version control, e.g., Mercurial, GIT
Friday:
Mercurial examples
R.J. LeVeque, University of Washington
AMath 483/
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 3 April 1, 2011
Today:
Unix
Mercurial examples
Monday:
Compiled languages, Fortran 90
R.J. LeVeque, University of Washington
AMath 483/583, Lecture 3, April 1, 2011
Bash shell
The shell is the program that prints a prompt in a ter
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 4 April 4, 2011
Today:
Interpreted vs. compiled languages
Fortran
Wednesday:
More Fortran
Computer storage of numbers
Read: Class notes and references.
R.J. LeVeque, University of Washington
AMath 483/583, Lecture 4, April 4, 201
Applied Linear Algebra and Introductory Numerical Analysis
AMATH 584

Summer 2011
AMath 483/583 Lecture 29 June 3, 2011
Today:
Reproducible research
Binary I/O
Animation: plots to movies
Sage
Parallel IPython
Course evaluations
Some new examples:
$CLASSHG/codes/io
$CLASSHG/codes/graphics
$CLASSHG/codes/python/mectest.py
R.J. LeVe