Anisotropy

[lsawade]

Reading the materials and adding a new material

In the IO/mesh/impl/read_materials.cpp a section needs to be added to
add the correct values to the new material this is done in the local function
read_materials(), which returns the specification. This would use the
implemented material class from above to populate the materials

• Add Anisotropic Material and Read it

  • Create Tests to Read Anisotropic Mesh and Run Anisotropic case for the Newmark Timescheme #239
  • include/enumerations/medium.hpp Add anisotropic medium
  • Add conditional to read new material properties read_materials #252
    IO/mesh/impl/read_materials.cpp::read_materials() from line 73ish
  • Add test to check whether the properties are correctly assigned. #264
  • include/mesh/materials/materials.hpp Add new material to the struct/class #258

• Move and implement anisotropic specfem::compute::impl::properties_container to specfem::medium::compute_properties #251

  • specfem::medium::compute_properties< ... > should store least possible implementation of the data required for computation of forces as Kokkos views.
  • Add tests for compute_properties. #277

• Move and implement specfem::point::properties to specfem::medium::point_properties #250

  • specfem::medium::point_properties< ... > should hold all the data required for computations
  • specfem::material::properties (-> specfem::medium::properties) needs specfem::medium::point_properties get_properties(...) method.

• specfem::compute::compute_mesh - has to be updated as well due to the new tags. Refactor below #281

• Update specfem::domain::impl::sources - Add anisotropic sources case. We need to refactor this - see below #279

• Implement specfem::domain::impl::receivers to implement anisotropic receivers. Refactor to use compute_wavefield routine. See below #280

Implement Wave propagation physics

• include/compute/properties/impl/properties_container.hpp Done above
• Create include/medium/elastic_anisotropic2d/elastic_anisotropic2d.hpp #270
• Implement Stacey boundary condition for anisotropic media. refactor see below #287
• Implement anisotropic mass matrix computations - Need refactor, see below #285
• Add explicit template instantiations for various anisotropic boundaries. Needs refactor see below #289
• include/domain/impl/kernels.hpp Need refactor. See #278  #276

Implement kernels

• include/compute/kernels/kernels.hpp #274
• include/compute/kernels/kernels/impl/kernels_container.hpp #275
• include/writer/kernel.tpp #299

Implement frechet derivatives

• Create frechet_derivatives/impl/element_kernel/elastic_anisotropic.hpp #295
• include/point/kernels.hpp #297
• Update naming in the kernel computation. It sounds like adjoint and backward wavefield are named incorrectly switched, and it just works because stiffness matrix is symmetric.
  rho_kl = dot(adjoint_field.acceleration, backward_field.displacement);
  Compare:
  
  So the naming is confusing.
• frechet derivative driver function #303

Examples

• Implement anisotropic crystal #306
• Implement isotropic/anisotropic example #324

Refactors

Refactor domain class

• Refactor kernels within domain class. #278
  • We currently assume that every material systems have implmenetation for all property types. Which is not the case i.e. (acoustic, anisotropic) is not a valid combination.
• Add anisotropic kernel tests

Refactor sources class

• Sources are independent of medium they lie in. i.e. acceleration(i) = stf(i) * lagrange_interpolant(i) for all components #339

Refactor receivers class

• Receiver implementation need to updated to use specfem::medium::compute_wavefield

Refactor medium

• Move and implement anisotropy to specfem::medium::properties #259

  • specfem::medium::properties should store exactly the values that are read from the mesh database

• Move specfem::domain::impl::element::impl_mass_matrix to specfem::medium::mass_matrix #319

• Move propeties_container kernels_container to medium namespace;

Refactor frechet kernels

• Refactor frechet kernels calling routine frechet_kernels.hpp \ .cpp #318
• Move specfem::frechet_derivatives::impl_compute_frechet_derivatives to specfem::medium::compute_frechet_derivatives #320

Refactor properties container and kernels container

• Move specfem::compute::impl::properties_container to specfem::medium::compute_properties
• Move specfem::compute::impl::kernels_container to specfem::medium::compute_kernels
• Make specfem::point::properties public parent of specfem::medium::point_properties
• Make specfem::point::kernels public parent of specfem::medium::point_kernels

Refactor compute_mesh

• Refactor compute mesh so that new tags are added automatically #292

Refactor src/domain/impl/boundary_conditions

• Refactor boundary condition cpp files so that explicit template instantiations are auto generated #313