2.096/6.336/16.910 Introduction to Numerical Simulation
Fall 2009
Instructor: Luca Daniel
TAs: Zohaib Mahmood, Tarek Ali El Moselhy, Brad Bond
Lecture: TR1-2.30 (35-225)
Information:
is an introduction to computational techniques for the simulation of a large variety of engineering and physical systems. Applications are drawn from aerospace, mechanical, electrical, and chemical engineering, biology, and materials science. Topics include mathematical formulations (techniques for automatic assembly of mathematical problems from physics' principles); sparse, direct and iterative solution techniques for linear systems; Newton methods for nonlinear problems; discretization methods for ordinary, time-periodic and partial differential equations; accelerated methods for integral equations; techniques for dynamical system model reduction.
Announcements
Problem Set 5
- For problem 1, you need a second initial condition (such as x[-1] or x[1]) to implement the method. There are many ways you can select a value for this. Describe and justify your choice in your writeup. There is no "right" answer here.
Announced on 15 November 2009 9:06 a.m. by Brad Bond
Guidline elaboration (submitting code)
Please include in your writeup (in the hard copy you turn in to us) important parts of matlab code you write or modify.
So what defines "important"? Changes you make to code we give you are important. Segments of code used to generate plots are not important. Use your best judgement. Anything that shows us how you solved the problem or will help you explain your algorithm is important.
For example, in Problem 5 of PS4, you do not need to include the entire function newton.m, you simply need to paste in the small portion that you add in order to solve the linear systems.
It's ok if you exceed the 10 page limit because of code printouts (but try to keep it under 15 pages including code).
Don't worry if you didn't include code in your writeup for PS3, we will not penalize you for that. But please include it in the future because it makes it easier for us to understand your explanations and to see where things go wrong if there are problems.
Also, please continue to submit all code via stellar in addition.
Announced on 29 October 2009 8:36 p.m. by Brad Bond
Problem Set Guidelines
- Be sure to answer all parts of the questions in your writeup. Include all plots that are asked for, and make comparisons to Matlab if asked for. Many people are not including things that are explicitly asked for.
- Include important code in your writeup. This means anything you write from scratch or any modifications you make to code that we give to you.
- List on the front page of your writeup the names of anybody you collaborated with on that assignment. This is very important!
- Even if you collaborate with others, you still must write your own code and perform your own analysis for every question.
- List any external references you used to help solve the problems. For example, if you found an algorithm on a website somewhere that was used to solve one of the problems, please mention this and include a citation.
- If you are stuck on a problem, ask the TAs for help. Do not leave the question blank, you will just have to redo it later.
Announced on 29 October 2009 11:52 a.m. by Brad Bond
Problem Set 4
- When using damping, you do not need to pick alpha to minimize ||F(xk+alfa*dx)||, you simply need to find an alpha such that ||F(xk+alfa*dx)|| < ||F(xk)||
- For 1c, you can consider functions other than polynomials and exponentials
Announced on 21 October 2009 12:38 p.m. by Brad Bond
Problem Set 3
- Oct 08, Modified problem 2 to eliminate a redundant question.
Announced on 07 October 2009 12:06 p.m. by Brad Bond