fnrizzi/ElasticShearWaves

Scope

This repository contains code for simulating elastic seismic shear waves in an axisymmetric domain. The implementation uses Kokkos and Kokkos-kernels to enable performance portability, and provides capabilities to run both the full order model (FOM) as well as the Galerkin reduced order model (ROM).

Specifically, this code has been developed mainly for this paper.

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

This repo is not being updated anymore. The code is being moved to https://github.com/Pressio/SHAW.

For all the latest developments, refer to that repo.

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Building

To build the code, you can use one of the following two alternatives:

• If you are an "expert" and want to use CMake directly, you can refer to this page.

• If you prefer a more guided procedure where we show you step-by-step how to install the TPLs, and provide you a way to build the code using helper bash scripts that we prepared, you can
  follow this step-by-step guide.

Step-by-step for running a full-order model simulation

Follow this guide for an example to run a FOM run.

Step-by-step for a rank-1 Galerkin ROM simulation

Follow this guide for an example to run a rank-1 Galerkin run.

Step-by-step for a rank-2 Galerkin ROM simulation

Coming soon.

High-level content

• C++ source code
• C++ tests
• meshing
• Python scripts for processing and workflows

C++ code references in detail

• FOM

  • main file
  • rank-1 problem class
  • rank-2 problem class
  • kernel call for velocity
  • kernel call for stresses
  • time loop driver

• ROM

  • main file
  • rank-1 problem class
  • rank-2 problem class
  • rank-1 solution loop
  • rank-2 solution loop

License

File is here.