From: kanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Thu, 7 Mar 2013 10:30:56 +0000 (+0000)
Subject: do less time steps for a regular run
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=cbca54da080972adb5b275b68ec7089ddbdb8e41;p=dealii-svn.git

do less time steps for a regular run
  and add comments about how to change back to the original

git-svn-id: https://svn.dealii.org/trunk@28779 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/examples/step-21/step-21.cc b/deal.II/examples/step-21/step-21.cc
index 861cbbf64a..806e83a01a 100644
--- a/deal.II/examples/step-21/step-21.cc
+++ b/deal.II/examples/step-21/step-21.cc
@@ -2,7 +2,7 @@
 
 /*    $Id$       */
 /*                                                                */
-/*    Copyright (C) 2006-2012 by the deal.II authors */
+/*    Copyright (C) 2006-2013 by the deal.II authors */
 /*                                                                */
 /*    This file is subject to QPL and may not be  distributed     */
 /*    without copyright and license information. Please refer     */
@@ -441,7 +441,9 @@ namespace Step21
 
   // The linear solvers we use are also completely analogous to the ones used
   // in step-20. The following classes are therefore copied verbatim from
-  // there. There is a single change: if the size of a linear system is small,
+  // there. Note that the classes here are not only copied from
+  // step-20, but also duplicate classes in deal.II. In a future version of this example, they should be
+  // replaced by an efficient method, though. There is a single change: if the size of a linear system is small,
   // i.e. when the mesh is very coarse, then it is sometimes not sufficient to
   // set a maximum of <code>src.size()</code> CG iterations before the solver
   // in the <code>vmult()</code> function converges. (This is, of course, a
@@ -1209,6 +1211,11 @@ namespace Step21
   // the present time step in the middle of solving the linear system
   // corresponding to each time step. We can therefore output the present end
   // time of a time step only at the end of the time step.
+  //
+  // The function as it is here does actually not compute the results
+  // found on the web page. The reason is, that even on a decent
+  // computer it runs more than a day. If you want to reproduce these
+  // results, set the final time at the end of the do loop to 250.
   template <int dim>
   void TwoPhaseFlowProblem<dim>::run ()
   {
@@ -1246,7 +1253,7 @@ namespace Step21
                   << std::endl
                   << std::endl;
       }
-    while (time <= 250);
+    while (time <= 1.);
   }
 }