LevelWaveOperator.H

Go to the documentation of this file.

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

#ifndef LEVEL_WAVEOPERATOR_H
#define LEVEL_WAVEOPERATOR_H

// Original design:                     Phil Colella
// Original implementation:
// Implementation date:

#include "FArrayBox.H"
#include "DisjointBoxLayout.H"
#include "LevelData.H"
#include "QuadCFInterp.H"
#include "ExtrapFillPatch.H"
#include "LevelFluxRegister.H"
#include "ProblemDomain.H"
#include "WaveOpF_F.H"
#include "UsingNamespace.H"


class LevelWaveOperator
{
public:

  LevelWaveOperator();


  ~LevelWaveOperator();


  void define(const DisjointBoxLayout&  a_thisDisjointBoxLayout,
              const DisjointBoxLayout&  a_coarserDisjointBoxLayout,
              const ProblemDomain&      a_domain,
              const int&                a_refineCoarse,
              const int&                a_numStates,
              const Real&               a_dx,
              const bool&               a_hasCoarser,
              const bool&               a_hasFiner);


  void eval(LevelData<FArrayBox>&       a_Phi,
            LevelData<FArrayBox>&       a_LOfPhi,
            LevelFluxRegister&          a_finerFluxRegister,
            LevelFluxRegister&          a_coarserFluxRegister,
            const LevelData<FArrayBox>& a_phiCoarseOld,
            const Real&                 a_TCoarseOld,
            const LevelData<FArrayBox>& a_phiCoarseNew,
            const Real&                 a_TCoarseNew,
            Real                        a_time,
            Real                        a_fluxweight);

  void avgdown(
            LevelData<FArrayBox>&       a_Phi,
            LevelData<FArrayBox>&       a_phiCoarse);

protected:
  // define() has been called
  bool m_defined;

  // Box layout for this level
  DisjointBoxLayout m_grids;

  // Interpolator for filling in ghost cells from the next coarser level
  QuadCFInterp m_patcher;

  // Grid spacing
  Real m_dx;

  // Problem domain - index space for this level
  ProblemDomain m_domain;

  // Refinement ratio between this level and the next coarser
  int m_refineCoarse;

  // Number of conservative fluxes and other face centered quantities
  int m_numFluxes;

  int m_numGhost;
  int m_numFields;

  // Flags for the existence of coarser and finer levels, respectively
  bool m_hasCoarser;
  bool m_hasFiner;

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

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

#endif

---