Anisotropy Implementation

.pre-commit-config.yaml

@@ -11,7 +11,7 @@ repos:
       - id: requirements-txt-fixer
       - id: trailing-whitespace
   - repo: https://github.com/pre-commit/mirrors-clang-format
-    rev: "v14.0.5"
+    rev: "v15.0.6"
     hooks:
       - id: clang-format
         types_or: [c++]

CMakeLists.txt

@@ -10,7 +10,7 @@ option(BUILD_TESTS "Tests included" OFF)
 option(BUILD_EXAMPLES "Examples included" ON)
 option(ENABLE_SIMD "Enable SIMD" OFF)
 option(ENABLE_PROFILING "Enable profiling" OFF)
-option(SPECFEMPP_ENABLE_PYTHON "Enable Python binding" OFF)
+option(SPECFEMPP_BINDING_PYTHON "Enable Python binding" OFF)
 # set(CMAKE_BUILD_TYPE Release)
 set(CHUNK_SIZE 32)
 set(NUM_CHUNKS 1)
@@ -25,7 +25,7 @@ set(CMAKE_RUNTIME_OUTPUT_DIRECTORY ${CMAKE_BINARY_DIR}/bin)
 # For external modules that need to be built suppress build output
 set(FETCHCONTENT_QUIET TRUE)
 
-if (SPECFEMPP_ENABLE_PYTHON)
+if (SPECFEMPP_BINDING_PYTHON)
     message("-- Adding -fPIC flag for Python binding.")
     if(CMAKE_CXX_COMPILER_ID STREQUAL "IntelLLVM")
             set(CMAKE_CXX_FLAGS "-fp-model=precise -fPIC")
@@ -36,15 +36,15 @@ if (SPECFEMPP_ENABLE_PYTHON)
     else()
             set(CMAKE_CXX_FLAGS "-fPIC")
     endif()
-else (SPECFEMPP_ENABLE_PYTHON)
+else (SPECFEMPP_BINDING_PYTHON)
     if(CMAKE_CXX_COMPILER_ID STREQUAL "IntelLLVM")
             set(CMAKE_CXX_FLAGS "-fp-model=precise")
     elseif(CMAKE_CXX_COMPILER_ID STREQUAL "Intel")
             message("-- Detected Intel classic compiler which will be deprecated soon.")
             message("-- It is recommended you use IntelLLVM compiler.")
             set(CMAKE_CXX_FLAGS "-diag-disable=10441 -fp-model=precise")
     endif()
-endif (SPECFEMPP_ENABLE_PYTHON)
+endif (SPECFEMPP_BINDING_PYTHON)
 
 # Check if MacOS
 if(${CMAKE_SYSTEM_NAME} MATCHES "Darwin")
@@ -110,8 +110,9 @@ find_package(Boost 1.73.0 COMPONENTS program_options filesystem system)
 
 if (NOT ${Boost_FOUND})
     # Add boost lib sources
-    set(BOOST_INCLUDE_LIBRARIES program_options filesystem system algorithm tokenizer)
-    set(BOOST_LIBS Boost::program_options Boost::filesystem Boost::system Boost::algorithm Boost::tokenizer)
+    set(BOOST_INCLUDE_LIBRARIES program_options filesystem system algorithm tokenizer preprocessor)
+    set(BOOST_LIBS Boost::program_options Boost::filesystem Boost::system
+                   Boost::algorithm Boost::tokenizer Boost::preprocessor)
     set(BOOST_ENABLE_CMAKE ON)
     set(BOOST_ENABLE_MPI OFF CACHE INTERNAL "Boost MPI Switch") # Assume outer variable
     set(BOOST_ENABLE_PYTHON OFF CACHE INTERNAL "Boost Python Switch") # Assume outer variable
@@ -193,8 +194,6 @@ target_link_libraries(
 add_library(
         IO
         src/IO/fortranio/fortran_io.cpp
-        src/IO/HDF5/native_type.cpp
-        # src/IO/ASCII/native_type.cpp
         src/IO/sources.cpp
         src/IO/receivers.cpp
         src/IO/mesh.cpp
@@ -248,6 +247,12 @@ add_library(
         src/enumerations/medium.cpp
 )
 
+target_link_libraries(
+        enumerations
+        ${BOOST_LIBS}
+)
+
+
 add_library(
         edge
         src/edge/interface.cpp
@@ -334,7 +339,7 @@ target_link_libraries(
 
 add_library(
         read_seismogram
-        src/reader/seismogram.cpp
+        src/IO/seismogram/reader.cpp
 )
 
 target_link_libraries(
@@ -344,7 +349,8 @@ target_link_libraries(
 
 add_library(
         reader
-        src/reader/wavefield.cpp
+        src/IO/wavefield/reader.cpp
+        src/IO/property/reader.cpp
 )
 
 target_link_libraries(
@@ -420,6 +426,8 @@ target_link_libraries(
 add_library(
         compute
         src/compute/compute_mesh.cpp
+        src/compute/impl/element_types.cpp
+        src/compute/impl/value_containers.cpp
         src/compute/compute_partial_derivatives.cpp
         src/compute/compute_properties.cpp
         src/compute/compute_kernels.cpp
@@ -450,20 +458,7 @@ target_link_libraries(
 
 add_library(
         domain
-        src/domain/impl/boundary_conditions/none/none.cpp
-        src/domain/impl/boundary_conditions/none/print.cpp
-        src/domain/impl/boundary_conditions/dirichlet/dirichlet.cpp
-        src/domain/impl/boundary_conditions/dirichlet/print.cpp
-        src/domain/impl/boundary_conditions/stacey/stacey.cpp
-        src/domain/impl/boundary_conditions/stacey/mass_matrix.cpp
-        src/domain/impl/boundary_conditions/stacey/print.cpp
-        src/domain/impl/boundary_conditions/composite_stacey_dirichlet/composite_stacey_dirichlet.cpp
-        src/domain/impl/boundary_conditions/composite_stacey_dirichlet/mass_matrix.cpp
-        src/domain/impl/boundary_conditions/composite_stacey_dirichlet/print.cpp
         src/domain/impl/boundary_conditions/boundary_conditions.cpp
-        src/domain/impl/elements/acoustic/acoustic2d.cpp
-        src/domain/impl/elements/elastic/elastic2d.cpp
-        src/domain/impl/elements/element.cpp
         src/domain/impl/elements/kernel.cpp
         src/domain/impl/kernels.cpp
         src/domain/domain.cpp
@@ -475,40 +470,28 @@ target_link_libraries(
         enumerations
 )
 
-add_library(coupled_interface
-        src/coupled_interface/coupled_interface.cpp
+add_library(
+        medium
+        src/medium/compute_mass_matrix.cpp
+        src/medium/compute_stress_integrand.cpp
 )
 
 target_link_libraries(
-        coupled_interface
+        medium
         Kokkos::kokkos
         compute
+        domain
+        enumerations
 )
 
-# add_library(
-#         kernels
-#         src/kernels/kernels.cpp
-# )
-
-# target_link_libraries(
-#         kernels
-#         Kokkos::kokkos
-#         compute
-#         coupled_interface
-# )
-
-add_library(
-        frechet_derivatives
-        src/frechet_derivatives/impl/element_kernel/elastic_isotropic.cpp
-        src/frechet_derivatives/impl/element_kernel/acoustic_isotropic.cpp
-        src/frechet_derivatives/impl/frechet_element.cpp
-        src/frechet_derivatives/frechet_derivatives.cpp
+add_library(coupled_interface
+        src/coupled_interface/coupled_interface.cpp
 )
 
 target_link_libraries(
-        frechet_derivatives
+        coupled_interface
         Kokkos::kokkos
-        point
+        compute
 )
 
 add_library(
@@ -520,7 +503,6 @@ target_link_libraries(
         kernels
         Kokkos::kokkos
         compute
-        frechet_derivatives
 )
 
 add_library(
@@ -546,13 +528,16 @@ target_link_libraries(
         Kokkos::kokkos
         timescheme
         domain
+        kernels
+        medium
 )
 
 add_library(
         writer
-        src/writer/seismogram.cpp
-        src/writer/wavefield.cpp
-        src/writer/kernel.cpp
+        src/IO/seismogram/writer.cpp
+        src/IO/wavefield/writer.cpp
+        src/IO/kernel/writer.cpp
+        src/IO/property/writer.cpp
 )
 
 target_link_libraries(
@@ -604,6 +589,7 @@ add_library(
         src/parameter_parser/writer/wavefield.cpp
         src/parameter_parser/writer/plot_wavefield.cpp
         src/parameter_parser/writer/kernel.cpp
+        src/parameter_parser/writer/property.cpp
 )
 
 target_link_libraries(
@@ -616,6 +602,7 @@ target_link_libraries(
         reader
         kernels
         domain
+        medium
         solver
         ${BOOST_LIBS}
 )
@@ -641,6 +628,7 @@ target_link_libraries(
         plotter
         reader
         domain
+        medium
         coupled_interface
         kernels
         solver
@@ -684,10 +672,10 @@ endif (DOXYGEN_FOUND)
 
 if (SPECFEMPP_USE_SKBUILD AND EXISTS ${SKBUILD_SCRIPTS_DIR})
     install(TARGETS specfem2d DESTINATION ${SKBUILD_SCRIPTS_DIR})
-    install(FILES ${CMAKE_BINARY_DIR}/xmeshfem2D DESTINATION ${SKBUILD_SCRIPTS_DIR} PERMISSIONS OWNER_READ OWNER_WRITE OWNER_EXECUTE)
+    install(FILES ${CMAKE_BINARY_DIR}/bin/xmeshfem2D DESTINATION ${SKBUILD_SCRIPTS_DIR} PERMISSIONS OWNER_READ OWNER_WRITE OWNER_EXECUTE)
 endif ()
 
-if (SPECFEMPP_ENABLE_PYTHON)
+if (SPECFEMPP_BINDING_PYTHON)
     find_package(Python REQUIRED COMPONENTS Interpreter Development.Module)
     find_package(pybind11 CONFIG REQUIRED)
     python_add_library(_core MODULE src/specfem2d.cpp WITH_SOABI)
@@ -702,15 +690,17 @@ if (SPECFEMPP_ENABLE_PYTHON)
         parameter_reader
         receiver_class
         writer
+        plotter
         reader
         domain
+        medium
         coupled_interface
         kernels
         solver
         ${BOOST_LIBS}
         pybind11::headers
     )
-    target_compile_definitions(_core PRIVATE SPECFEMPP_ENABLE_PYTHON)
+    target_compile_definitions(_core PRIVATE SPECFEMPP_BINDING_PYTHON)
     target_compile_definitions(_core PRIVATE VERSION_INFO=${PROJECT_VERSION})
     install(TARGETS _core DESTINATION specfempp_core)
-endif (SPECFEMPP_ENABLE_PYTHON)
+endif (SPECFEMPP_BINDING_PYTHON)

docs/api/IO/index.rst

@@ -12,9 +12,9 @@ formats for the reading of the mesh is binary, and for sources and receivers it
 is yaml.
 
 In addition to these basic read functions, there are also the two
-reader and writer classes, :cpp:class:`specfem::reader::wavefield` and
-:cpp:class:`specfem::writer::wavefield`, which support both HDF5 and ASCII I/O.
-And, to write seismograms, we can use :cpp:class:`specfem::writer::seismogram`.
+reader and writer classes, :cpp:class:`specfem::IO::wavefield_reader` and
+:cpp:class:`specfem::IO::wavefield_writer`, which support both HDF5 and ASCII I/O.
+And, to write seismograms, we can use :cpp:class:`specfem::IO::seismogram_writer`.
 Seismogram I/O is only supported in ASCII format thus far.
 
 The slightly-lower level functionality to read and write data to and from disk

docs/api/IO/reader/index.rst

@@ -4,7 +4,7 @@ Reader
 
 The ``reader`` class provides interfaces to read simulation data as output files.
 
-.. doxygenclass:: specfem::writer::writer
+.. doxygenclass:: specfem::IO::writer
 
 Types of readers
 ----------------

docs/api/IO/reader/wavefield.rst

@@ -4,5 +4,5 @@
 Wavefield Reader
 ================
 
-.. doxygenclass:: specfem::reader::wavefield
+.. doxygenclass:: specfem::IO::wavefield_reader
     :members:

docs/api/IO/writer/index.rst

@@ -4,7 +4,7 @@ Writer
 
 The ``writer`` class provides interfaces to write simulation data as output files.
 
-.. doxygenclass:: specfem::writer::writer
+.. doxygenclass:: specfem::IO::writer
    :members:
 
 Types of writers

docs/api/IO/writer/seismogram_writer.rst

@@ -4,5 +4,5 @@ Seismogram Writer
 
 The ``seismogram`` provides methods to write seismograms to a file.
 
-.. doxygenclass:: specfem::writer::seismogram
+.. doxygenclass:: specfem::IO::seismogram_writer
     :members:

docs/api/IO/writer/wavefield.rst

@@ -4,5 +4,5 @@
 Wavefield Writer
 ================
 
-.. doxygenclass:: specfem::writer::wavefield
+.. doxygenclass:: specfem::IO::wavefield_writer
     :members:

docs/api/material/index.rst

@@ -4,20 +4,20 @@
 Material Properties
 ===================
 
-.. doxygenclass:: specfem::material::material
+.. doxygenclass:: specfem::medium::material
    :members:
 
 
 Implementation Details
 ----------------------
 
-.. doxygenclass:: specfem::material::properties
+.. doxygenclass:: specfem::medium::properties
    :members:
 
-.. doxygenclass:: specfem::material::properties< specfem::element::medium_tag::elastic, specfem::element::property_tag::isotropic >
+.. doxygenclass:: specfem::medium::properties< specfem::element::medium_tag::elastic, specfem::element::property_tag::isotropic >
    :members:
    :private-members:
 
-.. doxygenclass:: specfem::material::properties< specfem::element::medium_tag::acoustic, specfem::element::property_tag::isotropic >
+.. doxygenclass:: specfem::medium::properties< specfem::element::medium_tag::acoustic, specfem::element::property_tag::isotropic >
    :members:
    :private-members: