]> https://gitweb.dealii.org/ - dealii.git/commitdiff
minor documentation cleanup
authorDenis Davydov <davydden@gmail.com>
Wed, 27 Sep 2017 08:34:53 +0000 (10:34 +0200)
committerDenis Davydov <davydden@gmail.com>
Wed, 27 Sep 2017 14:50:24 +0000 (16:50 +0200)
include/deal.II/lac/lapack_full_matrix.h
include/deal.II/lac/lapack_support.h

index 6bf20a3a8f52c4f09857c867373336b8b10a09c5..a1a2c413a58694bf9abeae396472aa041e4dd5ec 100644 (file)
@@ -413,19 +413,29 @@ public:
   /**
    * Compute the Cholesky factorization of the matrix using LAPACK function Xpotrf.
    *
-   * @note The factorization is stored in the lower diagonal.
+   * @note The factorization is stored in the lower-triangular part of the matrix.
    */
   void compute_cholesky_factorization ();
 
   /**
-   * Estimate reciprocal of the condition number (in 1-norm) of a SPD matrix
-   * using Cholesky factorization. This function can only be called if the matrix
-   * is already factorized.
+   * Estimate the reciprocal of the condition number $1/k(A)$ in $L_1$ norm ($1/(||A||_1 ||A^{-1}||_1)$)
+   * of a symmetric positive definite matrix using Cholesky factorization. This function can only
+   * be called if the matrix is already factorized.
    *
-   * @param a_norm Is the 1-norm of the matrix before calling Cholesky
-   * factorization.
+   * @note The condition number $k(A)$ can be used to estimate the numerical
+   * error related to the matrix inversion or the solution of the
+   * system of linear algebraic equations as
+   * <code>error = std::numeric_limits<Number>::epsilon * k</code>.
+   * Alternatively one can get the number of accurate digits
+   * <code>std::floor(std::log10(k))</code>.
+   *
+   * @note The function computes reciprocal of the condition number to
+   * avoid possible overflow if the matrix is nearly singular.
+   *
+   * @param[in] l1_norm Is the $L_1$ norm of the matrix before calling Cholesky
+   * factorization. It can be obtained by calling l1_norm().
    */
-  number reciprocal_condition_number(const number a_norm) const;
+  number reciprocal_condition_number(const number l1_norm) const;
 
   /**
    * Compute the determinant of a matrix. As it requires the LU factorization of
index 6b4032999a48c55e74539eb4af7faec004077d91..aa9dd7c3c16aabcfdc9f52d78c39f362a464f9ba 100644 (file)
@@ -176,7 +176,7 @@ namespace LAPACKSupport
    * operation.
    */
   DeclException1(ExcProperty, Property,
-                 << "The function cannot be called while the "
+                 << "The function cannot be called with a "
                  << property_name(arg1)
                  << " matrix.");
 

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.