Research Projects
Here is a list of possible projects. You are welcome to suggest a
project yourself, possibly related to an application arising in
your own research. The projects listed can usually be done in
various ways, so multiple people can work on different variants of
the same project. Most of the projects listed are an indication of
what might be done; the exact project will be decided after some
discussion with the instructor. More projects may be added, but
don't wait too long to decide, as you will need to start in time
to complete your project by the end of the term.
Linear Equations
- Adaptive GMRES method
- Hybrid BiCGSTAB/GMRES method
- Experiment with loss of orthogonality in CG
- Incomplete LU-type (ILU) preconditioners to improve convergence of
iterative methods
- Incomplete Choleski (IC) preconditioners to improve convergence
of iterative methods
- Block-ILU or -IC preconditioners
- Multi-grid method (various possibilities)
- Test MATLAB's sparse direct solvers
Eigenvalue Problems
- Test MATLAB's sparse eigenvalue solver
- Implement and test variant(s) of Jacobi-Davidson method
Nonlinear Equations
- Robust Newton methods (various possibilities)
- Newton-Krylov methods (hybrid Newton-GMRES type method)
- Broyden's method or other secant updating methods
Optimization
- Simplex method for linear programming
- Interior point method for linear programming
- Broyden's method or other secant updating methods
- Robust Newton methods (various possibilities)
Discretization of ODEs or PDEs
- One-dimensional, adaptive discretization of wave equation
- Discretization of one-dimensional channel flow
- Two-dimensional Galerkin finite element discretization of PDE on
quadrilaterals (various possibilities)
- Two-dimensional Galerkin finite element discretization of PDE on
triangles (various possibilities)
- Test various MATLAB ODE methods
- Test MATLAB's PDE toolbox
Applications
- Medical image reconstruction (various possibilities)
- Experiments with pattern formation in biology (various possibilities)
- Image processing using SVD, FFT, or wavelets