PhysIBC.H

Go to the documentation of this file.

#ifdef CH_LANG_CC
/*
*      _______              __
*     / ___/ /  ___  __ _  / /  ___
*    / /__/ _ \/ _ \/  V \/ _ \/ _ \
*    \___/_//_/\___/_/_/_/_.__/\___/
*    Please refer to Copyright.txt, in Chombo's root directory.
*/
#endif

//
//  PhysIBC.H
// ============
//
// Virtual base class through which a user specifies the initial and boundary
// conditions for a hyperbolic system of PDEs.
//

#ifndef _PHYSIBC_H_
#define _PHYSIBC_H_

#include "FArrayBox.H"
#include "REAL.H"
#include "LevelData.H"
#include "ProblemDomain.H"
#include "NamespaceHeader.H"


class PhysIBC
{
public:

  PhysIBC();


  virtual ~PhysIBC();


  virtual void define(const ProblemDomain& a_domain,
                      const Real&          a_dx);


  virtual PhysIBC* new_physIBC() = 0;


  virtual void initialize(LevelData<FArrayBox>& a_U) = 0;


  virtual void primBC(FArrayBox&            a_WGdnv,
                      const FArrayBox&      a_Wextrap,
                      const FArrayBox&      a_W,
                      const int&            a_dir,
                      const Side::LoHiSide& a_side,
                      const Real&           a_time) = 0;


  virtual
  void setBdrySlopes(FArrayBox&       a_dW,
                     const FArrayBox& a_W,
                     const int&       a_dir,
                     const Real&      a_time) = 0;


  virtual
  void artViscBC(FArrayBox&       a_F,
                 const FArrayBox& a_U,
                 const FArrayBox& a_divVel,
                 const int&       a_dir,
                 const Real&      a_time) = 0;


protected:

  virtual
  void getBoundaryFaces(
                        Box&                   a_boundaryBox,
                        const Box&             a_dataFaceBox,
                        const int&             a_dir,
                        const Side::LoHiSide&  a_side);

  // Has this object been defined
  bool          m_isDefined;

  // The current level's problem domain
  ProblemDomain m_domain;

  // The current level's grid spacing
  Real          m_dx;


private:
  // Disallowed for all the usual reasons
  void operator=(const PhysIBC& a_input)
  {
    MayDay::Error("invalid operator");
  }

  // Disallowed for all the usual reasons
  PhysIBC(const PhysIBC& a_input)
  {
    MayDay::Error("invalid operator");
  }
};

#include "NamespaceFooter.H"
#endif

---