]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Reformat a comment in LAPACKFullMatrix.
authorDavid Wells <wellsd2@rpi.edu>
Sat, 21 Nov 2015 23:04:01 +0000 (18:04 -0500)
committerDavid Wells <wellsd2@rpi.edu>
Sat, 21 Nov 2015 23:51:29 +0000 (18:51 -0500)
The information is the same but this makes it easier to read.

source/lac/lapack_full_matrix.cc

index 45b51c0293701026d229468532eedfda74b4e1a6..bb414c30cffdbfb2353374d2c6db9e6a66767e84 100644 (file)
@@ -635,11 +635,15 @@ LAPACKFullMatrix<number>::compute_eigenvalues(const bool right,
   const char *const jobvr = (right) ? (&V) : (&N);
   const char *const jobvl = (left)  ? (&V) : (&N);
 
-  // The LAPACK routine DGEEV requires a sufficient large workspace variable,
-  // minimum requirement is
-  //    work.size>=4*nn.
-  // However, to improve performance, a somewhat larger workspace may be
-  // needed.
+  /*
+   * The LAPACK routine xGEEV requires a sufficiently large work array; the
+   * minimum requirement is
+   *
+   * work.size >= 4*nn.
+   *
+   * However, for better performance, a larger work array may be needed. The
+   * first call determines the optimal work size and the second does the work.
+   */
   lwork = -1;
   work.resize(1);
 
@@ -701,11 +705,15 @@ LAPACKFullMatrix<number>::compute_eigenvalues_symmetric(const number        lowe
   std::vector<int> ifail(static_cast<size_type> (nn));
 
 
-  // The LAPACK routine ?SYEVX requires a sufficient large workspace variable,
-  // minimum requirement is
-  //    work.size>=3*nn-1.
-  // However, to improve performance, a somewhat larger workspace may be
-  // needed.
+  /*
+   * The LAPACK routine xSYEVX requires a sufficiently large work array; the
+   * minimum requirement is
+   *
+   * work.size >= 8*nn.
+   *
+   * However, for better performance, a larger work array may be needed. The
+   * first call determines the optimal work size and the second does the work.
+   */
   work.resize(1);
 
   syevx (jobz, range,
@@ -789,11 +797,15 @@ LAPACKFullMatrix<number>::compute_generalized_eigenvalues_symmetric(
   std::vector<int> ifail(static_cast<size_type> (nn));
 
 
-  // The LAPACK routine ?SYGVX requires a sufficient large workspace variable,
-  // minimum requirement is
-  //    work.size>=3*nn-1.
-  // However, to improve performance, a somewhat larger workspace may be
-  // needed.
+  /*
+   * The LAPACK routine xSYGVX requires a sufficiently large work array; the
+   * minimum requirement is
+   *
+   * work.size >= 8*nn.
+   *
+   * However, for better performance, a larger work array may be needed. The
+   * first call determines the optimal work size and the second does the work.
+   */
   work.resize(1);
 
   sygvx (&itype, jobz, range, uplo, &nn, values_A, &nn,
@@ -868,11 +880,15 @@ LAPACKFullMatrix<number>::compute_generalized_eigenvalues_symmetric (
   const char *const jobz = (eigenvectors.size() > 0) ? (&V) : (&N);
   const char *const uplo = (&U);
 
-  // The LAPACK routine DSYGV requires a sufficient large workspace variable,
-  // minimum requirement is
-  //    work.size>=3*nn-1.
-  // However, to improve performance, a somewhat larger workspace may be
-  // needed.
+  /*
+   * The LAPACK routine xSYGV requires a sufficiently large work array; the
+   * minimum requirement is
+   *
+   * work.size >= 3*nn - 1.
+   *
+   * However, for better performance, a larger work array may be needed. The
+   * first call determines the optimal work size and the second does the work.
+   */
   work.resize(1);
 
   sygv (&itype, jobz, uplo, &nn, values_A, &nn,

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.