]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix some matrix coefficients in the step-25 intro.
authorDavid Wells <wellsd2@rpi.edu>
Sat, 25 Jun 2016 23:27:34 +0000 (19:27 -0400)
committerDavid Wells <wellsd2@rpi.edu>
Sat, 25 Jun 2016 23:42:17 +0000 (19:42 -0400)
The above formula has the correct sign on N. Additionally, the Laplace
matrix here is missing the coefficients.

examples/step-25/doc/intro.dox

index 03ae995ab2953fdd57f42743f56f64a51201eb7a..f1a5761816e6e47b872aa14e1c9b2c4d0904a034 100644 (file)
@@ -176,9 +176,10 @@ i.e. $N_{ij}(f,g) = \left( \cos\left[ \theta f + (1-\theta) g\right]\varphi_i,
 What solvers can we use for the first equation? Let's look at the matrix we
 have to invert:
 @f[
-  (M+k^2\theta^2(A-N))_{ij} = 
-  \int_\Omega (1+k^2\theta^2 \cos \alpha) 
-  \varphi_i\varphi_j \; dx+\int_\Omega \nabla\varphi_i\nabla\varphi_j \; dx,
+  (M+k^2\theta^2(A + N))_{ij} =
+  \int_\Omega (1+k^2\theta^2 \cos \alpha)
+  \varphi_i\varphi_j \; dx
+  + k^2 \theta^2 \int_\Omega \nabla\varphi_i\nabla\varphi_j \; dx,
 @f]
 for some $\alpha$ that depends on the present and previous solution. First,
 note that the matrix is symmetric. In addition, if the time step $k$ is small

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.