From cbca54da080972adb5b275b68ec7089ddbdb8e41 Mon Sep 17 00:00:00 2001 From: kanschat Date: Thu, 7 Mar 2013 10:30:56 +0000 Subject: [PATCH] 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 --- deal.II/examples/step-21/step-21.cc | 13 ++++++++++--- 1 file changed, 10 insertions(+), 3 deletions(-) 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 src.size() CG iterations before the solver // in the vmult() 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 void TwoPhaseFlowProblem::run () { @@ -1246,7 +1253,7 @@ namespace Step21 << std::endl << std::endl; } - while (time <= 250); + while (time <= 1.); } } -- 2.39.5