Minor updates to the IDA documentation and test cases.

diff --git a/include/deal.II/sundials/ida.h b/include/deal.II/sundials/ida.h
index 0a392906bbb0264ee1f0fd38a67871449403a9ac..8dc48791aae46753708c521b9d6271e31c4920e5 100644 (file)
--- a/include/deal.II/sundials/ida.h
+++ b/include/deal.II/sundials/ida.h
@@ -151,10 +151,10 @@ namespace SUNDIALS
    * \end{matrix}
    * \f]
    *
-   * That is $F(y', y, t) = y' + A y = 0 $
+   * That is, $F(y', y, t) = y' + A y = 0 $
    * where
-   * A =
    * \f[
+   * A =
    * \begin{pmatrix}
    * 0 & -1 \\
    * k^2 &0
@@ -206,8 +206,8 @@ namespace SUNDIALS
    * {
    *   J = A;
    *
-   *   J(0,0) = alpha;
-   *   J(1,1) = alpha;
+   *   J(0,0) += alpha;
+   *   J(1,1) += alpha;
    *
    *   Jinv.invert(J);
    * };

diff --git a/tests/sundials/ida_01.cc b/tests/sundials/ida_01.cc
index afbf86750843c469fb4e18fb46ff0e5d8f1b90f4..1df9c52d708a84f5f696da2f90f99ef714ffa502 100644 (file)
--- a/tests/sundials/ida_01.cc
+++ b/tests/sundials/ida_01.cc
@@ -30,7 +30,7 @@
  * u (0) = 0
  * u'(0) = k
  *
- * write in terms of a first order ode:
+ * written in terms of a first order ode:
  *
  * y[0]' -     y[1]  = 0
  * y[1]' + k^2 y[0]  = 0
@@ -89,8 +89,8 @@ public:
 
       J = A;
 
-      J(0, 0) = alpha;
-      J(1, 1) = alpha;
+      J(0, 0) += alpha;
+      J(1, 1) += alpha;
 
       Jinv.invert(J);
     };

diff --git a/tests/sundials/ida_02.cc b/tests/sundials/ida_02.cc
index 0af53390ab1a293964d1f8cb74f05b389d47107a..ecc7882dea24fb494578b73e4df4bc955edff18a 100644 (file)
--- a/tests/sundials/ida_02.cc
+++ b/tests/sundials/ida_02.cc
@@ -33,7 +33,7 @@
  * u (0) = 0
  * u'(0) = k
  *
- * write in terms of a first order ode:
+ * written in terms of a first order ode:
  *
  * y[0]' -     y[1]  = 0
  * y[1]' + k^2 y[0]  = 0
@@ -92,8 +92,8 @@ public:
 
       J = A;
 
-      J(0, 0) = alpha;
-      J(1, 1) = alpha;
+      J(0, 0) += alpha;
+      J(1, 1) += alpha;
     };
 
     time_stepper.solve_with_jacobian =