ElasticSolver.h

Go to the documentation of this file.

//
// ExaHyPE2 solver file
// Generated by Peano's Python API
// www.peano-framework.org
//
// This is generated. If you change fundamental properties, you will have to
// generate this file. Backup your manual changes before you do so.
//
#pragma once
 
 
 
#include "AbstractElasticSolver.h"
#include "VariableShortcuts.h"
#include "tarch/logging/Log.h"
 
#define _CUSTOM_COORDINATES
#define ASAGI_NOMPI
#define ISOTROPIC
#define _TOP 1
 
#include "../ExaSeis_core/Context/ContextCurvilinearPML.h"
#include "../ExaSeis_core/Numerics/routines_curvilinear.h"
#include "../ExaSeis_core/Numerics/pml.h"
#include "../ExaSeis_core/Refinement/refinement.h"
 
namespace exahype2 {
  namespace elastic {
 
    class ElasticSolver;
    class VariableShortcutsExpanded;
 
  } // namespace elastic
} // namespace exahype2
 
class exahype2::elastic::VariableShortcutsExpanded: public ::exahype2::elastic::VariableShortcuts {
public:
  static constexpr int NumberOfUnknowns = AbstractElasticSolver::NumberOfUnknowns;
  static constexpr int NumberOfAuxiliaryVariables = AbstractElasticSolver::NumberOfAuxiliaryVariables;
  static constexpr int Size = NumberOfUnknowns + NumberOfAuxiliaryVariables;
};
 
class exahype2::elastic::ElasticSolver: public AbstractElasticSolver {
private:
  static tarch::logging::Log _log;
 
  VariableShortcuts s;
 
  ContextCurvilinear<VariableShortcutsExpanded, Order + 1>*                context;
  std::vector<Refinement::RefinementCriterion<VariableShortcutsExpanded>*> refinementCriteria;
  int                                                              pml_cell_width;
  double                                                           pml_alpha_const;
  double                                                           pml_alpha_scalar;
  double                                                           pml_rel_error;
  int                                                              pml_power;
 
public:
  ElasticSolver();
  double* QuadraturePoints1d = kernels::ElasticSolver::Quadrature<double>::nodes;
  ::exahype2::RefinementCommand refinementCriterion(
    const double* __restrict__ Q, // Q[36+19]
    const tarch::la::Vector<Dimensions, double>& x,
    const tarch::la::Vector<Dimensions, double>& h,
    double                                       t
  ) override;
 
  void initialCondition(
    double* __restrict__ Q, // Q[36+19]
    const tarch::la::Vector<Dimensions, double>& x,
    const tarch::la::Vector<Dimensions, double>& h,
    bool                                         gridIsConstructed
  ) override;
 
  void boundaryConditions(
    const double* __restrict__ Qinside, // Qinside[36+19]
    double* __restrict__ Qoutside,      // Qoutside[36+19]
    const tarch::la::Vector<Dimensions, double>& faceCentre,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t,
    int                                          normal
  ) final override;
 
  void riemannSolver(
    double* const                                FL, // FL[36
    double* const                                FR, // FR[36
    const double* const                          QL, // QL[36+19]
    const double* const                          QR, // QR[36+19]
    const double                                 t,
    const double                                 dt,
    const tarch::la::Vector<Dimensions, double>& x,
    const tarch::la::Vector<Dimensions, double>& h,
    const int                                    direction,
    bool                                         isBoundaryFace,
    int                                          faceIndex
  ) final override;
 
  double maxEigenvalue(
    const double* __restrict__ Q,
    const tarch::la::Vector<Dimensions, double>& volumeX,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t,
    double                                       dt,
    int                                          normal
  ) override;
 
  void flux(
    const double* __restrict__ Q, // Q[36+19]
    const tarch::la::Vector<Dimensions, double>& volumeX,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t,
    double                                       dt,
    int                                          normal,
    double* __restrict__ F // F[36]
  ) override;
 
  void nonconservativeProduct(
    const double* __restrict__ Q,
    const double* __restrict__ deltaQ,
    const tarch::la::Vector<Dimensions, double>& faceCentre,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t,
    double                                       dt,
    int                                          normal,
    double* __restrict__ BgradQ
  ) override;
 
  void algebraicSource(
    const double* const Q,
    const tarch::la::Vector<Dimensions,double>&  x,
    const tarch::la::Vector<Dimensions,double>&  h,
    double                                       t,
    double                                       dt,
    double* S
  ) override;
 
  void multiplyMaterialParameterMatrix(
    const double* const                          Q,
    const tarch::la::Vector<Dimensions, double>& volumeX,
    const tarch::la::Vector<Dimensions, double>& volumeH,
    double                                       t,
    double                                       dt,
    int                                          normal,
    double* const                                rhs
  ) override;
 
  void pointSource(
    const double* const Q, const double* const x, const double t, const double dt, double* const forceVector, int n
  ) override;
 
  void initPointSourceLocations(double sourceLocation[NumberOfPointSources][Dimensions]) override;
};