Allow set_mesh_smoothing() on non-empty triangulations.

diff --git a/doc/news/changes/minor/20230822Fehling b/doc/news/changes/minor/20230822Fehling
new file mode 100644 (file)
index 0000000..68fe48e
--- /dev/null
+++ b/doc/news/changes/minor/20230822Fehling
@@ -0,0 +1,4 @@
+Changed: You can now set MeshSmoothing flags in non-empty triangulations
+with Triangulation::set_mesh_smoothing().
+<br>
+(Marc Fehling, 2023/08/22)

diff --git a/include/deal.II/grid/tria.h b/include/deal.II/grid/tria.h
index 5dac93b8a203f974105eadfb45c0b09e10e36704..548a120d85d24bf081bd64bcef123d55722d265c 100644 (file)
--- a/include/deal.II/grid/tria.h
+++ b/include/deal.II/grid/tria.h
@@ -1849,8 +1849,7 @@ public:
 
   /**
    * Set the mesh smoothing to @p mesh_smoothing. This overrides the
-   * MeshSmoothing given to the constructor. It is allowed to call this
-   * function only if the triangulation is empty.
+   * MeshSmoothing given to the constructor.
    */
   virtual void
   set_mesh_smoothing(const MeshSmoothing mesh_smoothing);

diff --git a/source/grid/tria.cc b/source/grid/tria.cc
index 83218e36319b694bf04e5a47995b200f58d68e67..4e1d7c0595bf85809efc48b0d3b2ec36aa51ccff 100644 (file)
--- a/source/grid/tria.cc
+++ b/source/grid/tria.cc
@@ -12259,8 +12259,6 @@ DEAL_II_CXX20_REQUIRES((concepts::is_valid_dim_spacedim<dim, spacedim>))
 void Triangulation<dim, spacedim>::set_mesh_smoothing(
   const MeshSmoothing mesh_smoothing)
 {
-  Assert(n_levels() == 0,
-         ExcTriangulationNotEmpty(vertices.size(), levels.size()));
   smooth_grid = mesh_smoothing;
 }