Update assignment 6 for Fall 2010.

diff --git a/assignments/archive/poisson/index.shtml.itex2MML b/assignments/archive/poisson/index.2009.11.shtml.itex2MML
similarity index 100%
rename from assignments/archive/poisson/index.shtml.itex2MML
rename to assignments/archive/poisson/index.2009.11.shtml.itex2MML

diff --git a/assignments/archive/poisson/index.shtml b/assignments/archive/poisson/index.shtml
new file mode 100644 (file)
index 0000000..1f0b249
--- /dev/null
+++ b/assignments/archive/poisson/index.shtml
@@ -0,0 +1,30 @@
+<!--#set var="root_directory" value="../../.." --><!--#include virtual="$root_directory/shared/header.shtml"-->
+
+<h1>Assignment #6</h1>
+<p><em>Due Friday, November 19</em></p>
+
+<h2>Purpose</h2>
+
+<p>Learn how to implement Gauss-Siedel in parallel.</p>
+
+<p>Note: Please identify all your work.</p>
+
+<p>This assignment consists in building your own parallel code to
+generate the solution of Poisson equation.
+The <a href="../../../content/poisson/#Gauss-Siedel">Gauss-Siedel</a>
+and <a href="../../../content/poisson/#SOR">SOR</a> algorithms are
+explained in the course web pages. The steps toward a parallel code
+are also explained.</p>
+
+<p>Start from the code <a href="../../../src/possion/p.c">p.c</a> and
+implement the parallel version. Pay particular attention to minimize
+communications.</p>
+
+<ol>
+  <li>Produce a correct, parallel code.</li>
+  <li>Instrument both the serial and parallel codes for timing.</li>
+  <li>Find the time dependence of both codes on the <em>grid size</em>.</li>
+  <li>How efficient is your parallel version?</li>
+</ol>
+
+<!--#include virtual="$root_directory/shared/footer.shtml"-->

diff --git a/assignments/current/6 b/assignments/current/6
index 986997746d4d4cce64da818dead3d7a2d479d02f..80c8f41c1b2a844418da739564bad372d96e8f0e 120000 (symlink)
--- a/assignments/current/6
+++ b/assignments/current/6
@@ -1 +1 @@
-../archive/mandelbrot/
\ No newline at end of file
+../archive/poisson
\ No newline at end of file

diff --git a/content/2D_Domain_Decomposition/index.shtml.itex2MML b/content/2D_Domain_Decomposition/index.shtml.itex2MML
index 5094a33303727226742df071e525c6b67f8860e9..415c1857a9f5cba5c974916a89c1a2220f2278df 100644 (file)
--- a/content/2D_Domain_Decomposition/index.shtml.itex2MML
+++ b/content/2D_Domain_Decomposition/index.shtml.itex2MML
@@ -169,13 +169,13 @@ identical. This stems from the handling of the <code>GHOST</code>
 lines in the parallel implementation.</p>
 
 <p>Consider the update of the field at the location marked 0 below.
-The Gauss-Siedel solution calls for an <em>in-place</em> update of the
-field via the finite difference formula in a systematic sweep of the
-numerical lattice.  This is done via the <a href="#iter_2d">2D</a>
-finite difference formula. The values of the field at the points 2 and
-3 have already been updated if the sweep originates from the
-upper-left corner of the lattice, while the field at the point 1 and 4
-are yet to be updated.</p>
+The <a href="../poisson/#Gauss-Siedel">Gauss-Siedel</a> solution calls
+for an <em>in-place</em> update of the field via the finite difference
+formula in a systematic sweep of the numerical lattice.  This is done
+via the <a href="#iter_2d">2D</a> finite difference formula. The
+values of the field at the points 2 and 3 have already been updated if
+the sweep originates from the upper-left corner of the lattice, while
+the field at the point 1 and 4 are yet to be updated.</p>
 
 <img border="0" src="img/Red_Black_1.gif" width="221" height="223" /> <!-- TODO: convert to SVG -->
 

diff --git a/content/poisson/index.shtml.itex2MML b/content/poisson/index.shtml.itex2MML
index e39fbcce05c450bd78dda027ed612b1964c279a8..9e84cb8c8d306a39c481fd9e9074383556030f41 100644 (file)
--- a/content/poisson/index.shtml.itex2MML
+++ b/content/poisson/index.shtml.itex2MML
@@ -169,9 +169,9 @@ field $u_{i,j}$ at the interior grid points.</p>
   <embed type="image/svg+xml" src="img/jacobi.svg" width="190" height="174" />
 </object>
 
-<p>The Gauss-Siedel method consists in replacing the new value of the
-field within the same "array". This corresponds to the following
-form:</p>
+<p id="Gauss-Siedel">The Gauss-Siedel method consists in replacing the
+new value of the field within the same "array". This corresponds to
+the following form:</p>
 
 <p class="equation">\[
   u_{i,j}^{k+1} =
@@ -203,7 +203,7 @@ the regions of high curvature of the field and on the ability of the
 numerical grid to cover the domain.</p>
 
 <p id="SOR">The Gauss-Siedel method converges faster than the Jacobi
-method. Yet it is a notoriously slow method, often requiring very
+method. However, it is a notoriously slow method, often requiring very
 large number of iterations to achieve convergence. A better rate of
 convergence is achieved by the Succesive Over Relaxation (SOR)
 method. In this approach, the old and new fields are mixed via a