]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fix a bug with the computation of the right hand side contributions.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 12 Sep 2008 15:52:10 +0000 (15:52 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 12 Sep 2008 15:52:10 +0000 (15:52 +0000)
git-svn-id: https://svn.dealii.org/trunk@16820 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/doc/news/changes.h
deal.II/examples/step-23/step-23.cc

index 393600768981c7ecab0efc920492403f4911a24d..ab79646a640c55d72187ac76c482e72648669ba6 100644 (file)
@@ -112,6 +112,17 @@ inconvenience this causes.
 <h3>General</h3>
 
 <ol>
+  <li>
+  <p>
+  Fixed: In the @ref step_23 "step-23" tutorial program the terms corresponding
+  to the external force were computed incorrectly. The error wasn't visible
+  right away since in the program the right hand side is assumed to be zero.
+  This has now been fixed.
+  <br>
+  (Johan Lorentzon, WB 2008/09/12)
+  </p>
+  </li>
+
   <li>
   <p>
   Fixed: In the new @ref step_33 "step-33" tutorial program there was
index b0e5499299223e48fd7257389789205b5686910b..0c1552440225a6cf293301e81a1ebbe8b8081636 100644 (file)
@@ -4,7 +4,7 @@
 /*    $Id$       */
 /*    Version: $Name:  $                                          */
 /*                                                                */
-/*    Copyright (C) 2006, 2007 by the deal.II authors */
+/*    Copyright (C) 2006, 2007, 2008 by the deal.II authors */
 /*                                                                */
 /*    This file is subject to QPL and may not be  distributed     */
 /*    without copyright and license information. Please refer     */
@@ -647,13 +647,12 @@ void WaveEquation<dim>::run ()
       forcing_terms = tmp;
       forcing_terms *= theta * time_step;
       
-      system_rhs.add (theta * time_step, forcing_terms);
-
       rhs_function.set_time (time-time_step);
       VectorTools::create_right_hand_side (dof_handler, QGauss<dim>(2),
                                           rhs_function, tmp);
 
       forcing_terms.add ((1-theta) * time_step, tmp);
+
       system_rhs.add (theta * time_step, forcing_terms);
 
                                       // After so constructing the right hand

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.