Merge pull request #13505 from peterrum/ac_is_closed

Add AffineConstraints::is_closed()
dealii · Mar 8, 2022 · a5d36df · a5d36df
2 parents 21aaeb8 + 23a8fd7
commit a5d36df
Show file tree

Hide file tree

Showing 2 changed files with 33 additions and 0 deletions.
diff --git a/include/deal.II/lac/affine_constraints.h b/include/deal.II/lac/affine_constraints.h
@@ -770,6 +770,22 @@ class AffineConstraints : public Subscriptor
   void
   close();
 
+  /**
+   * Check if the function close() was called or there are no
+   * constraints locally, which is normally the case if a dummy
+   * AffineConstraints was created for the DG case.
+   */
+  bool
+  is_closed() const;
+
+  /**
+   * Check if the function close() was called or there are no
+   * constraints globally, which is normally the case if a dummy
+   * AffineConstraints was created for the DG case.
+   */
+  bool
+  is_closed(const MPI_Comm &comm) const;
+
   /**
    * Merge the constraints represented by the object given as argument into
    * the constraints represented by this object. Both objects may or may not

diff --git a/include/deal.II/lac/affine_constraints.templates.h b/include/deal.II/lac/affine_constraints.templates.h
@@ -938,6 +938,23 @@ AffineConstraints<number>::close()
 
 
 
+template <typename number>
+bool
+AffineConstraints<number>::is_closed() const
+{
+  return sorted || (n_constraints() == 0);
+}
+
+
+
+template <typename number>
+bool
+AffineConstraints<number>::is_closed(const MPI_Comm &comm) const
+{
+  return Utilities::MPI::min(static_cast<unsigned int>(is_closed()), comm) == 1;
+}
+
+
 template <typename number>
 void
 AffineConstraints<number>::merge(