From 9623f0699dc57153db109d4e34150b92d1222640 Mon Sep 17 00:00:00 2001 From: bangerth Date: Fri, 12 Sep 2008 15:52:10 +0000 Subject: [PATCH] Fix a bug with the computation of the right hand side contributions. git-svn-id: https://svn.dealii.org/trunk@16820 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/doc/news/changes.h | 11 +++++++++++ deal.II/examples/step-23/step-23.cc | 5 ++--- 2 files changed, 13 insertions(+), 3 deletions(-) diff --git a/deal.II/doc/news/changes.h b/deal.II/doc/news/changes.h index 3936007689..ab79646a64 100644 --- a/deal.II/doc/news/changes.h +++ b/deal.II/doc/news/changes.h @@ -112,6 +112,17 @@ inconvenience this causes.

General

    +
  1. +

    + 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. +
    + (Johan Lorentzon, WB 2008/09/12) +

    +
  2. +
  3. Fixed: In the new @ref step_33 "step-33" tutorial program there was diff --git a/deal.II/examples/step-23/step-23.cc b/deal.II/examples/step-23/step-23.cc index b0e5499299..0c15524402 100644 --- a/deal.II/examples/step-23/step-23.cc +++ b/deal.II/examples/step-23/step-23.cc @@ -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::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(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 -- 2.39.5