MA 5490, Fall 2010
Professor Craig C. Douglas
Tuesday - Thursday 1:20-2:35, Ross Hall 247
Notes   Homework   Syllabus

Homework and Projects

All homework should be emailed to me. I require that you send me your files as a single zip or gzipped tar file using the naming convention, or lastname.tgz unless it is a single file, which should lastname.extension, where extension makes sense based on the type of file. Absolutely no .rar files!

I will send you a reply when I get your email. If you do not get a reply, I probably did not get your email. You are responsible for having an email account that can communicate successfully with me. I am flexible on which of my email accounts you use up to a degree. I will send you email from GMail, however. Please send email to my GMail account since that is where mail to my university address will forward to eventually.

Description Category Due Worth
hw1 homework 8/31/2010 0
hw2 homework 9/2/2010 0
hw3 homework 10/5/2010 0
hw4 homework 11/2/2010 0
hw5 homework 12/?/2010 0


This will be detailed later in the course.


Fill out the course questionaire and email it back to me as an attachment. Note that it is a .doc Word document. Please return it as such (not pdf, docx, txt, rtf, xls, or anything else).


Read the two papers emailed and mark them up with what you do not understand, which will probably a lot. Do not worry that you do not understand the whole paper now.


Work out on paper the finite difference problems from 9/30/2010's class.


This is a continuation of Hw3's first problem only. Do the first problem over again from Hw3 with variable mesh spacing only (case H2) and use finite elements. Do not try to do this assignment in a finite difference approach by using multiple elements simultaneously (which is a complete mistake in finite element philosophy will result in having to do the assignment all over again). Do the following steps in order (similar to what is done for the CHIP problem in Chapter 4 of my textbook):

  1. Choose piecewise linear basis functions.
  2. Define a global element.
  3. Define mappings in both directions between a local and the global element.
  4. Define all of the shape functions on the global element and their drivatives.
  5. Define a complete set of element stiffness matrices on the global element.
  6. Define a complete set of element load vectors on the global element.
  7. Generate a sample mesh with at least 5 elements (note the word elements rather than nodal points). Create the local/global numbering table for the nodes, too.
  8. Show all work for generating the global stiffness matrix.
  9. Show all work for generating the global load vector.
  10. Correct the global stiffness matrix for each of the boundary conditions in Hw3's problem 1.
  11. Correct the global load vector for each of the boundary conditions in Hw3's problem 1.

If you really want to understand what is going on, additionally do the problem with piecewise constant basis functions (this is really quick to do since the basis functions are utterly trivial to work with).


Implement Algorithm 2.4 from chapter two of my dissertation. Do this in steps:

  1. Choose an elliptic boundary value problem in 1D including boundary conditions.
  2. Choose finite elements or differences as the discretization method.
  3. Implement generating a stiffness matrix and load vector (finite elements) or finite difference matrix and right hand side vector for a given mesh.
  4. Generate a single level problem (the coarse level per se) so that there is one unknown.
  5. Implement a smoother and show that it works on the single level problem in no more than one iteration (thus, it is a direct solver when only one unknown exists).
  6. Refine the mesh by diving each mesh element in two equal parts (corresponding to sigma = 2 in the algorithm).
  7. Implement a linear interpolation algorithm for Ej.
  8. Generate a second level problem (without losing the coarse level problem).
  9. Implement a residual calculation algorithm for a problem on level j.
  10. Implement a weighted averaging (or L2 projection) method for the right hand side of the next coarser level.
  11. Finish implementing the multigrid algorithm based only on the pieces implemented above.

Do steps 1-5 by Tuesday, November 23. Be prepared to discuss your implementation in class and have questions. If necessary, we will have an evening class to clarify the homework.


Craig C. Douglas

Last modified: