Hhs1_Transposed.h

Go to the documentation of this file.

#ifndef EXASEIS_SCENARIOHHSONE_TRANSPOSED_HEADER
#define EXASEIS_SCENARIOHHSONE_TRANSPOSED_HEADER
#include "Scenario.h"
#include "kernels/GaussLegendreBasis.h"
 
template <class Shortcuts, int basisSize>
class Hhs1_Transposed: public Scenario<Shortcuts, basisSize> {
private:
  double m_pointSourceLocation[3];
 
public:
  Hhs1_Transposed(DomainInformation* info):
    Scenario<Shortcuts, basisSize>(info) {
    // place point source in center of topography 2km below surface
    m_pointSourceLocation[0] = 0.693;
    m_pointSourceLocation[1] = 0.0;
    m_pointSourceLocation[2] = 0.0;
  };
 
  void initUnknownsPointwise(
    const double* const                          x,
    const tarch::la::Vector<Dimensions, double>& center,
    const double                                 t,
    const double                                 dt,
    double*                                      Q
  ) {
    Shortcuts s;
 
#ifdef ISOTROPIC
    Q[s.rho] = 2.7;
    Q[s.cp]  = 6.0;
    Q[s.cs]  = 3.464;
#elif defined ANISOTROPIC
    assertion1(false, "Hhs1_transposed Not implemented in the Anisotropic case");
#else
#error Whole space problem only defined for Isotropic and Anisotropic Scenario
#endif
  }
 
  void initPointSourceLocation(double pointSourceLocation[][3]) {
    pointSourceLocation[0][0] = m_pointSourceLocation[0];
    pointSourceLocation[0][1] = m_pointSourceLocation[1];
    pointSourceLocation[0][2] = m_pointSourceLocation[2];
  }
 
  void setPointSourceVector(
    const double* const Q, const double* const x, const double t, const double dt, double* forceVector, int n
  ) {
 
    assertion2(n == 0, "Only a single pointSource for Hhs1_Transposed_transposed", n);
 
    constexpr double t0 = 0.1;
    constexpr double M0 = 1000.0;
 
    double f = M0 * t / (t0 * t0) * std::exp(-t / t0);
 
    Shortcuts s;
    forceVector[s.sigma + 5] = f;
  }
 
#if defined _EXASEIS_DYNAMIC
  void refinementCriteria(std::vector<Refinement::RefinementCriterion<Shortcuts>*>& criteria
  ) override {
    Shortcuts s;
 
    criteria.push_back(new Refinement::RefinePosition<Shortcuts>(m_pointSourceLocation));
    criteria.push_back(new Refinement::trackVelocity<Shortcuts>(basisSize, 0.5, 0.2));
 
    // criteria.push_back(new Refinement::CoarseBoundaryLayer<Shortcuts>(
    //   solver->getMaximumAdaptiveMeshLevel(),
    //   solver->getCoarsestMeshLevel(),
    //   solver->getCoarsestMeshSize(),
    //   0,
    //   &(solver->getDomainSize()[0]),
    //   &(solver->getDomainOffset()[0])
    // ));
 
    // criteria.push_back(new Refinement::RefineDownToPosition<Shortcuts>(
    //   basisSize, solver->getMaximumAdaptiveMeshLevel(), m_pointSourceLocation, 0
    // ));
  }
#else
  void refinementCriteria(std::vector<Refinement::RefinementCriterion<Shortcuts>*>& criteria) override {
    Shortcuts s;
 
    criteria.push_back(new Refinement::StaticAMR<Shortcuts>);
    // criteria.push_back(new Refinement::CoarseBoundaryLayer<Shortcuts>(
    //   solver->getMaximumAdaptiveMeshLevel(),
    //   solver->getCoarsestMeshLevel(),
    //   solver->getCoarsestMeshSize(),
    //   0,
    //   &(solver->getDomainSize()[0]),
    //   &(solver->getDomainOffset()[0])
    // ));
#ifdef _CUSTOM_COORDINATES
    // criteria.push_back(new Refinement::RefineDownToPositionCustomCoordinates<Shortcuts>(
    //   basisSize,
    //   solver->getMaximumAdaptiveMeshLevel(),
    //   m_pointSourceLocation,
    //   0,
    //   static_cast<exahype::solvers::ADERDGSolver*>(solver)->FCoeff[basisSize - 1][0]
    // ));
#else
    // criteria.push_back(new Refinement::RefineDownToPosition<Shortcuts>(
    //   basisSize, solver->getMaximumAdaptiveMeshLevel(), m_pointSourceLocation, 0
    // ));
#endif
  };
#endif
};
#endif