| 901 | #endif |
| 902 | } |
| 903 | bool seissol::initializer::isElasticSideOfElasticAcousticInterface(CellMaterialData &material, |
| 904 | unsigned int face) { |
| 905 | #ifdef USE_ANISOTROPIC |
| 906 | return false; |
| 907 | #else |
| 908 | constexpr auto eps = std::numeric_limits<real>::epsilon(); |
| 909 | return material.local.mu > eps && material.neighbor[face].mu < eps; |
| 910 | #endif |
| 911 | } |
| 912 | |
| 913 | bool seissol::initializer::isAtElasticAcousticInterface(CellMaterialData &material, unsigned int face) { |
| 914 | // We define the interface cells as all cells that are in the elastic domain but have a |