From a1de92582d1a6ef728bf73fd92836d03ab4f2acc Mon Sep 17 00:00:00 2001
From: Rene Gassmoeller <rene.gassmoeller@mailbox.org>
Date: Wed, 29 May 2024 16:44:46 -0600
Subject: [PATCH] Check that composition type is valid

---
 source/simulator/introspection.cc | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/source/simulator/introspection.cc b/source/simulator/introspection.cc
index e7c0e0d220d..4296c4896bc 100644
--- a/source/simulator/introspection.cc
+++ b/source/simulator/introspection.cc
@@ -439,6 +439,7 @@ namespace aspect
   const std::vector<unsigned int> &
   Introspection<dim>::get_indices_for_fields_of_type (const CompositionalFieldDescription::Type &type) const
   {
+    Assert(type < CompositionalFieldDescription::n_types, ExcInternalError());
     return composition_indices_for_type[type];
   }
 
@@ -448,6 +449,7 @@ namespace aspect
   const std::vector<std::string> &
   Introspection<dim>::get_names_for_fields_of_type (const CompositionalFieldDescription::Type &type) const
   {
+    Assert(type < CompositionalFieldDescription::n_types, ExcInternalError());
     return composition_names_for_type[type];
   }
 
@@ -457,6 +459,7 @@ namespace aspect
   unsigned int
   Introspection<dim>::get_number_of_fields_of_type (const CompositionalFieldDescription::Type &type) const
   {
+    Assert(type < CompositionalFieldDescription::n_types, ExcInternalError());
     return composition_indices_for_type[type].size();
   }