Introduce LAPACKFullMatrix::get_state()

diff --git a/doc/news/changes/minor/20241015Munch-1 b/doc/news/changes/minor/20241015Munch-1
new file mode 100644 (file)
index 0000000..b62e30c
--- /dev/null
+++ b/doc/news/changes/minor/20241015Munch-1
@@ -0,0 +1,6 @@
+Improved: The new function LAPACKFullMatrix::get_state() allows to
+query the current state of LAPACKFullMatrix, which allows to decide
+which method (e.g., LAPACKFullMatrix::solve() vs.
+LAPACKFullMatrix::vmult()) to use.
+<br>
+(Peter Munch, 2024/10/15)

diff --git a/include/deal.II/lac/lapack_full_matrix.h b/include/deal.II/lac/lapack_full_matrix.h
index aa15249f94027cd59f610951182fcc2e15c40b54..7d9596077a2db0c7a6e0ccacb4e15c007a5d6da8 100644 (file)
--- a/include/deal.II/lac/lapack_full_matrix.h
+++ b/include/deal.II/lac/lapack_full_matrix.h
@@ -917,6 +917,12 @@ public:
                   const double       threshold   = 0.,
                   const char        *separator   = " ") const;
 
+  /**
+   * Return current state after the last operation here.
+   */
+  LAPACKSupport::State
+  get_state() const;
+
 private:
   /**
    * Internal function to compute various norms.

diff --git a/source/lac/lapack_full_matrix.cc b/source/lac/lapack_full_matrix.cc
index 4f09243e36bf2c087424b44b8d4938b03a0e8aa8..38d52117bfe436dc287aa210df1cba8f38315353 100644 (file)
--- a/source/lac/lapack_full_matrix.cc
+++ b/source/lac/lapack_full_matrix.cc
@@ -2574,6 +2574,15 @@ LAPACKFullMatrix<number>::print_formatted(std::ostream      &out,
 }
 
 
+
+template <typename number>
+LAPACKSupport::State
+LAPACKFullMatrix<number>::get_state() const
+{
+  return this->state;
+}
+
+
 //----------------------------------------------------------------------//
 
 template <typename number>