]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Comments to the source code added.
authorschrage <schrage@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 7 Jun 1999 14:31:52 +0000 (14:31 +0000)
committerschrage <schrage@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 7 Jun 1999 14:31:52 +0000 (14:31 +0000)
git-svn-id: https://svn.dealii.org/trunk@1380 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/doc/tutorial/chapter-1.elements/rhs.html

index 9490c77f1c3e18fc7dc5ad23e3a241601f1214f2..b0e8da5431b6b4e707a96bda7aa65fb55c77a8e8 100644 (file)
@@ -119,24 +119,27 @@ Vector&lt;double&gt; elvec(fe.total_dofs);
 FullMatrix&lt;double&gt; elmat(fe.total_dofs);
   
 for (DoFHandler&lt;2&gt;::active_cell_iterator c = dof.begin_active()
-                                               ; c != dof.end() ; ++c)
+       ; c != dof.end() ; ++c)
 {
   fevalues.reinit(c, stb);
   elmat.clear();
   elvec.clear();
   c-&gt;get_dof_indices(indices);
-    
+
+  // The loop over the quadrature points starts here. In effect, the 
+  // integration is performed in this loop.
   for (unsigned k=0;k&lt;qc.n_quadrature_points;++k)
   {
     for (unsigned i=0;i&lt;fe.total_dofs;++i)
     {
-      const Point&lt;2&gt; dv = fevalues.shape_grad(i,k);
+       const Point&lt;2&gt; dv = fevalues.shape_grad(i,k);
 
        
       for (unsigned j=0;j&lt;fe.total_dofs;++j)
       {
        const Point&lt;2&gt; du = fevalues.shape_grad(j,k);
-         
+       
+       // Perform the integration using the discretized Laplace operator.
        elmat(i,j) += fevalues.JxW(k)
                       * du * dv
                      ;

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.