Enable TriaAccessor::diameter() for wedges

diff --git a/include/deal.II/grid/tria_accessor.templates.h b/include/deal.II/grid/tria_accessor.templates.h
index aeb62ddc6f20ffb39b78d65ce3c6ade4fce2360f..1653dd49b010bf280637f4f98ea47f40a1fccb93 100644 (file)
--- a/include/deal.II/grid/tria_accessor.templates.h
+++ b/include/deal.II/grid/tria_accessor.templates.h
@@ -1904,6 +1904,14 @@ TriaAccessor<structdim, dim, spacedim>::diameter() const
                             (this->vertex(3) - this->vertex(0)).norm())),
           std::max((this->vertex(3) - this->vertex(1)).norm(),
                    (this->vertex(3) - this->vertex(2)).norm()));
+      case ReferenceCell::Type::Wedge:
+        return std::max(
+          std::max(std::max((this->vertex(4) - this->vertex(0)).norm(),
+                            (this->vertex(3) - this->vertex(1)).norm()),
+                   std::max((this->vertex(5) - this->vertex(1)).norm(),
+                            (this->vertex(4) - this->vertex(2)).norm())),
+          std::max((this->vertex(5) - this->vertex(0)).norm(),
+                   (this->vertex(3) - this->vertex(2)).norm()));
       case ReferenceCell::Type::Hex:
         return std::max(std::max((this->vertex(7) - this->vertex(0)).norm(),
                                  (this->vertex(6) - this->vertex(1)).norm()),