From 378aafdb4beb73747f8b9deca6eec0744c913e6e Mon Sep 17 00:00:00 2001 From: bangerth Date: Tue, 20 May 2008 16:47:02 +0000 Subject: [PATCH] Document the new compute_forcing_vector function. git-svn-id: https://svn.dealii.org/trunk@16138 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-33/step-33.cc | 16 ++++++++++++++++ 1 file changed, 16 insertions(+) diff --git a/deal.II/examples/step-33/step-33.cc b/deal.II/examples/step-33/step-33.cc index 44da0a71fd..d539a1232d 100644 --- a/deal.II/examples/step-33/step-33.cc +++ b/deal.II/examples/step-33/step-33.cc @@ -370,6 +370,22 @@ struct EulerEquations } + // In the same way as describing the flux + // function $\mathbf F(\mathbf w), we + // also need to have a way to describe + // the right hand side forcing term. As + // mentioned in the introduction, we + // consider only gravity here, which + // leads to the specific form $\mathbf + // G(\mathbf w) = \left( \begin{array}{c} + // g_1\rho \\ g_2\rho \\ g_3\rho \\ 0 + // \\ \rho \mathbf g \cdot \mathbf v + // \end{array} \right)$, shown here for + // the 3d case. More specifically, we + // will consider only $\mathbf + // g=(0,0,-1)^T$ in 3d, or $\mathbf + // g=(0,-1)^T$ in 2d. This naturally + // leads to the following function: template static void compute_forcing_vector (const InputVector &W, -- 2.39.5