Updated elastic tag to elastic_sv tag

PrincetonUniversity · Feb 19, 2025 · e14cf16 · e14cf16
1 parent c5d4792
commit e14cf16
Show file tree

Hide file tree

Showing 72 changed files with 440 additions and 410 deletions.
diff --git a/include/IO/kernel/writer.tpp b/include/IO/kernel/writer.tpp
@@ -37,7 +37,7 @@ void specfem::IO::kernel_writer<OutputLibrary>::write(specfem::compute::assembly
     typename OutputLibrary::Group elastic = file.createGroup("/ElasticIsotropic");
 
     const auto element_indices = element_types.get_elements_on_host(
-        specfem::element::medium_tag::elastic,
+        specfem::element::medium_tag::elastic_sv,
         specfem::element::property_tag::isotropic);
     n_elastic_isotropic = element_indices.size();
 
@@ -60,7 +60,7 @@ void specfem::IO::kernel_writer<OutputLibrary>::write(specfem::compute::assembly
           const specfem::point::index<specfem::dimension::type::dim2> index(
               ispec, iz, ix);
           specfem::point::kernels<specfem::dimension::type::dim2,
-                                  specfem::element::medium_tag::elastic,
+                                  specfem::element::medium_tag::elastic_sv,
                                   specfem::element::property_tag::isotropic,
                                   false>
               point_kernels;
@@ -91,7 +91,7 @@ void specfem::IO::kernel_writer<OutputLibrary>::write(specfem::compute::assembly
     typename OutputLibrary::Group elastic = file.createGroup("/ElasticAnisotropic");
 
     const auto element_indices = element_types.get_elements_on_host(
-        specfem::element::medium_tag::elastic,
+        specfem::element::medium_tag::elastic_sv,
         specfem::element::property_tag::anisotropic);
     n_elastic_anisotropic = element_indices.size();
 
@@ -115,7 +115,7 @@ void specfem::IO::kernel_writer<OutputLibrary>::write(specfem::compute::assembly
           const specfem::point::index<specfem::dimension::type::dim2> index(
               ispec, iz, ix);
           specfem::point::kernels<specfem::dimension::type::dim2,
-                                  specfem::element::medium_tag::elastic,
+                                  specfem::element::medium_tag::elastic_sv,
                                   specfem::element::property_tag::anisotropic,
                                   false>
               point_kernels;

diff --git a/include/IO/property/writer.tpp b/include/IO/property/writer.tpp
@@ -37,7 +37,7 @@ void specfem::IO::property_writer<OutputLibrary>::write(specfem::compute::assemb
     typename OutputLibrary::Group elastic = file.createGroup("/ElasticIsotropic");
 
     const auto element_indices = element_types.get_elements_on_host(
-        specfem::element::medium_tag::elastic,
+        specfem::element::medium_tag::elastic_sv,
         specfem::element::property_tag::isotropic);
     n_elastic_isotropic = element_indices.size();
 
@@ -66,7 +66,7 @@ void specfem::IO::property_writer<OutputLibrary>::write(specfem::compute::assemb
     typename OutputLibrary::Group elastic = file.createGroup("/ElasticAnisotropic");
 
     const auto element_indices = element_types.get_elements_on_host(
-        specfem::element::medium_tag::elastic,
+        specfem::element::medium_tag::elastic_sv,
         specfem::element::property_tag::anisotropic);
     n_elastic_anisotropic = element_indices.size();
 

diff --git a/include/boundary_conditions/stacey/stacey.tpp b/include/boundary_conditions/stacey/stacey.tpp
@@ -9,9 +9,9 @@
 #include <type_traits>
 
 namespace {
-using elastic_type =
+using elastic_sv_type =
     std::integral_constant<specfem::element::medium_tag,
-                           specfem::element::medium_tag::elastic>;
+                           specfem::element::medium_tag::elastic_sv>;
 
 using acoustic_type =
     std::integral_constant<specfem::element::medium_tag,
@@ -112,7 +112,7 @@ template <
     typename PointFieldType, typename ViewType,
     typename std::enable_if_t<!PointBoundaryType::simd::using_simd, int> = 0>
 KOKKOS_FUNCTION void
-impl_enforce_traction(const elastic_type &, const isotropic_type &,
+impl_enforce_traction(const elastic_sv_type &, const isotropic_type &,
                       const PointBoundaryType &boundary,
                       const PointPropertyType &property,
                       const PointFieldType &field, ViewType &traction) {
@@ -122,7 +122,7 @@ impl_enforce_traction(const elastic_type &, const isotropic_type &,
                 "Boundary tag must be stacey");
 
   static_assert(PointPropertyType::medium_tag ==
-                    specfem::element::medium_tag::elastic,
+                    specfem::element::medium_tag::elastic_sv,
                 "Medium tag must be elastic");
 
   static_assert(PointPropertyType::property_tag ==
@@ -164,7 +164,7 @@ template <
     typename PointFieldType, typename ViewType,
     typename std::enable_if_t<PointBoundaryType::simd::using_simd, int> = 0>
 KOKKOS_FUNCTION void
-impl_enforce_traction(const elastic_type &, const isotropic_type &,
+impl_enforce_traction(const elastic_sv_type &, const isotropic_type &,
                       const PointBoundaryType &boundary,
                       const PointPropertyType &property,
                       const PointFieldType &field, ViewType &traction) {
@@ -174,7 +174,7 @@ impl_enforce_traction(const elastic_type &, const isotropic_type &,
                 "Boundary tag must be stacey");
 
   static_assert(PointPropertyType::medium_tag ==
-                    specfem::element::medium_tag::elastic,
+                    specfem::element::medium_tag::elastic_sv,
                 "Medium tag must be elastic");
 
   static_assert(PointPropertyType::property_tag ==
@@ -224,7 +224,7 @@ template <
     typename PointFieldType, typename ViewType,
     typename std::enable_if_t<!PointBoundaryType::simd::using_simd, int> = 0>
 KOKKOS_FUNCTION void
-impl_enforce_traction(const elastic_type &, const anisotropic_type &,
+impl_enforce_traction(const elastic_sv_type &, const anisotropic_type &,
                       const PointBoundaryType &boundary,
                       const PointPropertyType &property,
                       const PointFieldType &field, ViewType &traction) {
@@ -234,7 +234,7 @@ impl_enforce_traction(const elastic_type &, const anisotropic_type &,
                 "Boundary tag must be stacey");
 
   static_assert(PointPropertyType::medium_tag ==
-                    specfem::element::medium_tag::elastic,
+                    specfem::element::medium_tag::elastic_sv,
                 "Medium tag must be elastic");
 
   static_assert(PointPropertyType::property_tag ==
@@ -277,7 +277,7 @@ template <
     typename PointFieldType, typename ViewType,
     typename std::enable_if_t<PointBoundaryType::simd::using_simd, int> = 0>
 KOKKOS_FUNCTION void
-impl_enforce_traction(const elastic_type &, const anisotropic_type &,
+impl_enforce_traction(const elastic_sv_type &, const anisotropic_type &,
                       const PointBoundaryType &boundary,
                       const PointPropertyType &property,
                       const PointFieldType &field, ViewType &traction) {
@@ -287,7 +287,7 @@ impl_enforce_traction(const elastic_type &, const anisotropic_type &,
                 "Boundary tag must be stacey");
 
   static_assert(PointPropertyType::medium_tag ==
-                    specfem::element::medium_tag::elastic,
+                    specfem::element::medium_tag::elastic_sv,
                 "Medium tag must be elastic");
 
   static_assert(PointPropertyType::property_tag ==

diff --git a/include/compute/boundary_values/boundary_values_container.hpp b/include/compute/boundary_values/boundary_values_container.hpp
@@ -27,7 +27,7 @@ class boundary_value_container {
       DimensionType, specfem::element::medium_tag::acoustic, BoundaryTag>
       acoustic;
   specfem::compute::impl::boundary_medium_container<
-      DimensionType, specfem::element::medium_tag::elastic, BoundaryTag>
+      DimensionType, specfem::element::medium_tag::elastic_sv, BoundaryTag>
       elastic;
 
   boundary_value_container() = default;
@@ -89,7 +89,7 @@ store_on_device(const int istep, const IndexType index,
   if constexpr (MediumTag == specfem::element::medium_tag::acoustic) {
     boundary_value_container.acoustic.store_on_device(istep, l_index,
                                                       acceleration);
-  } else if constexpr (MediumTag == specfem::element::medium_tag::elastic) {
+  } else if constexpr (MediumTag == specfem::element::medium_tag::elastic_sv) {
     boundary_value_container.elastic.store_on_device(istep, l_index,
                                                      acceleration);
   }
@@ -123,7 +123,7 @@ load_on_device(const int istep, const IndexType index,
   if constexpr (MediumType == specfem::element::medium_tag::acoustic) {
     boundary_value_container.acoustic.load_on_device(istep, l_index,
                                                      acceleration);
-  } else if constexpr (MediumType == specfem::element::medium_tag::elastic) {
+  } else if constexpr (MediumType == specfem::element::medium_tag::elastic_sv) {
     boundary_value_container.elastic.load_on_device(istep, l_index,
                                                     acceleration);
   }

diff --git a/include/compute/boundary_values/boundary_values_container.tpp b/include/compute/boundary_values/boundary_values_container.tpp
@@ -27,7 +27,7 @@ specfem::compute::boundary_value_container<DimensionType, BoundaryTag>::
       nstep, mesh, element_types, boundaries, h_property_index_mapping);
 
   elastic = specfem::compute::impl::boundary_medium_container<
-      DimensionType, specfem::element::medium_tag::elastic, BoundaryTag>(
+      DimensionType, specfem::element::medium_tag::elastic_sv, BoundaryTag>(
       nstep, mesh, element_types, boundaries, h_property_index_mapping);
 
   Kokkos::deep_copy(property_index_mapping, h_property_index_mapping);

diff --git a/include/compute/coupled_interfaces/coupled_interfaces.hpp b/include/compute/coupled_interfaces/coupled_interfaces.hpp
@@ -61,8 +61,9 @@ struct coupled_interfaces {
   specfem::compute::interface_container<medium1, medium2>
   get_interface_container() const;
 
-  specfem::compute::interface_container<specfem::element::medium_tag::elastic,
-                                        specfem::element::medium_tag::acoustic>
+  specfem::compute::interface_container<
+      specfem::element::medium_tag::elastic_sv,
+      specfem::element::medium_tag::acoustic>
       elastic_acoustic; ///< Elastic-acoustic interface
 
   specfem::compute::interface_container<
@@ -71,7 +72,7 @@ struct coupled_interfaces {
       acoustic_poroelastic; ///< Acoustic-poroelastic interface
 
   specfem::compute::interface_container<
-      specfem::element::medium_tag::elastic,
+      specfem::element::medium_tag::elastic_sv,
       specfem::element::medium_tag::poroelastic>
       elastic_poroelastic; ///< Elastic-poroelastic interface
 };

diff --git a/include/compute/coupled_interfaces/coupled_interfaces.tpp b/include/compute/coupled_interfaces/coupled_interfaces.tpp
@@ -10,11 +10,11 @@ template <specfem::element::medium_tag medium1,
           specfem::element::medium_tag medium2>
 specfem::compute::interface_container<medium1, medium2>
 specfem::compute::coupled_interfaces::get_interface_container() const {
-  if constexpr (medium1 == specfem::element::medium_tag::elastic &&
+  if constexpr (medium1 == specfem::element::medium_tag::elastic_sv &&
                 medium2 == specfem::element::medium_tag::acoustic) {
     return elastic_acoustic;
   } else if constexpr (medium1 == specfem::element::medium_tag::acoustic &&
-                       medium2 == specfem::element::medium_tag::elastic) {
+                       medium2 == specfem::element::medium_tag::elastic_sv) {
     return specfem::compute::interface_container<medium1, medium2>(
         elastic_acoustic);
   } else if constexpr (medium1 == specfem::element::medium_tag::acoustic &&
@@ -24,11 +24,11 @@ specfem::compute::coupled_interfaces::get_interface_container() const {
                        medium2 == specfem::element::medium_tag::acoustic) {
     return specfem::compute::interface_container<medium1, medium2>(
         acoustic_poroelastic);
-  } else if constexpr (medium1 == specfem::element::medium_tag::elastic &&
+  } else if constexpr (medium1 == specfem::element::medium_tag::elastic_sv &&
                        medium2 == specfem::element::medium_tag::poroelastic) {
     return elastic_poroelastic;
   } else if constexpr (medium1 == specfem::element::medium_tag::poroelastic &&
-                       medium2 == specfem::element::medium_tag::elastic) {
+                       medium2 == specfem::element::medium_tag::elastic_sv) {
     return specfem::compute::interface_container<medium1, medium2>(
         elastic_poroelastic);
   }