Line data    Source code

       1              : //
       2              : // Class of BC operators that work with :ref:`Operator::Elastic`.
       3              : //
       4              : 
       5              : #ifndef BC_OPERATOR_ELASTIC_H
       6              : #define BC_OPERATOR_ELASTIC_H
       7              : 
       8              : // #include "Operator/Elastic.H"
       9              : #include "IO/ParmParse.H"
      10              : #include "Numeric/Interpolator/Linear.H"
      11              : 
      12              : namespace BC
      13              : {
      14              : namespace Operator
      15              : {
      16              : namespace Elastic
      17              : {
      18              : class Elastic
      19              : {
      20              : public:
      21              :     static const int min = 1, max = 2;
      22              : 
      23           19 :     virtual ~Elastic() = default;
      24              : 
      25              : public:
      26              :     enum Type {Displacement, Traction, Periodic, Neumann, None}; 
      27              : 
      28              :     #if AMREX_SPACEDIM==2
      29              :     enum Face{
      30              :         XLO, YLO, XHI, YHI, 
      31              :         XLO_YLO, XLO_YHI, XHI_YLO, XHI_YHI,
      32              :         INT
      33              :     };
      34              :     #elif AMREX_SPACEDIM==3
      35              :     enum Face{
      36              :         XLO, YLO, ZLO, XHI, YHI, ZHI, // 6
      37              :         YLO_ZLO, YLO_ZHI, YHI_ZLO, YHI_ZHI, // 12
      38              :         ZLO_XLO, ZLO_XHI, ZHI_XLO, ZHI_XHI,
      39              :         XLO_YLO, XLO_YHI, XHI_YLO, XHI_YHI,
      40              :         XLO_YLO_ZLO, XLO_YLO_ZHI, XLO_YHI_ZLO, XLO_YHI_ZHI, // 8
      41              :         XHI_YLO_ZLO, XHI_YLO_ZHI, XHI_YHI_ZLO, XHI_YHI_ZHI,
      42              :         INT
      43              :     };
      44              :     #endif
      45              : 
      46              :     enum Direction {AMREX_D_DECL(X=0,Y=1,Z=2)}; 
      47              : 
      48              :     void 
      49          426 :     SetTime(const Set::Scalar a_time) {m_time = a_time;}
      50              : 
      51              :     virtual void
      52              :     Init(amrex::MultiFab * a_rhs,
      53              :         const amrex::Geometry &a_geom,
      54              :         bool a_homogeneous = false) const = 0;
      55              : 
      56              :     virtual void
      57              :     Init(amrex::FabArray<amrex::BaseFab<Set::Vector>> * a_rhs,
      58              :         const amrex::Geometry &a_geom,
      59              :         bool a_homogeneous = false) const = 0;
      60              : 
      61              :     void
      62              :     Init(Set::Field<Set::Scalar> &a_rhs,
      63              :         const amrex::Vector<amrex::Geometry> &a_geom,
      64              :         bool a_homogeneous = false) const
      65              :     {
      66              :         for (int ilev = 0; ilev <= a_rhs.finest_level; ilev++)
      67              :             Init(a_rhs[ilev].get(),a_geom[ilev],a_homogeneous);
      68              :     }
      69              : 
      70              :     void
      71          426 :     Init(Set::Field<Set::Vector> &a_rhs,
      72              :         const amrex::Vector<amrex::Geometry> &a_geom,
      73              :         bool a_homogeneous = false) const
      74              :     {
      75          898 :         for (int ilev = 0; ilev <= a_rhs.finest_level; ilev++)
      76          472 :             Init(a_rhs[ilev].get(),a_geom[ilev],a_homogeneous);
      77          426 :     }
      78              : 
      79              : #define SQRT3INV 0.57735026919
      80              : #define SQRT2INV 0.70710678118
      81              : 
      82              :     virtual
      83              :     std::array<Type,AMREX_SPACEDIM> getType (
      84              :                 const int &i, const int &j, const int &k,
      85              :                 const amrex::Box &domain) = 0;
      86              : 
      87              :     virtual
      88              :     Set::Vector operator () (const Set::Vector &u,
      89              :                 const Set::Matrix &gradu,
      90              :                 const Set::Matrix &sigma,
      91              :                 const int &i, const int &j, const int &k,
      92              :                 const amrex::Box &domain) = 0;
      93              : 
      94              : protected:
      95              :     Set::Scalar m_time = 0.0;
      96              : };
      97              : }
      98              : }
      99              : }
     100              : #endif