EBFluxFAB.H

Go to the documentation of this file.

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

// DTG  7-31-2001

#ifndef _EBFLUXFAB_H_
#define _EBFLUXFAB_H_

#include "Box.H"
#include "Vector.H"
#include "EBFaceFAB.H"
#include "NamespaceHeader.H"


class EBFluxFAB
{

public:


  static int preAllocatable()
  {
    return 1; // symmetric messaging
  }

  EBFluxFAB();

  ~EBFluxFAB();

  void clear();

  EBFluxFAB& operator+=(const EBFluxFAB& a_ebfab)
  {
    for(int idir = 0; idir < SpaceDim; idir++)
      {
        (*m_fluxes[idir]) += (*a_ebfab.m_fluxes[idir]);
      }
    return *this;
  }


  EBFluxFAB& operator-=(const EBFluxFAB& a_ebfab)
  {
    for(int idir = 0; idir < SpaceDim; idir++)
      {
        (*m_fluxes[idir]) -= (*a_ebfab.m_fluxes[idir]);
      }
    return *this;
  }


  EBFluxFAB& operator*=(const EBFluxFAB& a_ebfab)
  {
    for(int idir = 0; idir < SpaceDim; idir++)
      {
        (*m_fluxes[idir]) *= (*a_ebfab.m_fluxes[idir]);
      }
    return *this;
  }


  EBFluxFAB& operator/=(const EBFluxFAB& a_ebfab)
  {
    for(int idir = 0; idir < SpaceDim; idir++)
      {
        (*m_fluxes[idir]) /= (*a_ebfab.m_fluxes[idir]);
      }
    return *this;
  }


  int nComp() const;

  const Box& getRegion() const;

  EBFaceFAB& operator[] (const int dir);

  const EBFaceFAB& operator[] (const int dir) const;

  void setVal(const Real& val);


  void copy(const Box& Rfrom,
            const Interval& Cdest,
            const Box& Rto,
            const EBFluxFAB& src,
            const Interval& Csrc);

  bool isDefined() const;

  EBFluxFAB(const EBISBox& a_ebisBox, const Box& bx, int n);

  void define(const EBISBox& a_ebisBox, const Box& bx, int n);

  int size(const Box& R, const Interval& comps) const;

  void linearOut(void* buf, const Box& R, const Interval& comps) const;

  void linearIn(void* buf, const Box& R, const Interval& comps);


  void alias(Vector<EBFaceFAB*> a_inputFAB);

  EBFluxFAB(const Box& bx, int n)
  {
    MayDay::Error("invalid constructor called for EBFluxFAB");
  }
  EBISBox getEBISBox() const
  {
    return m_ebisBox;
  }

protected:
  void setDefaultValues();

  bool m_aliased;
  Box m_region;

  EBISBox m_ebisBox;

  int m_nComp;

  bool m_isDefined;

  Tuple<EBFaceFAB*, CH_SPACEDIM> m_fluxes;


private:

  EBFluxFAB (const EBFluxFAB&)
  {
    MayDay::Error("bogus constructor");
  }

  void operator = (const EBFluxFAB&)
  {
    MayDay::Error("bogus constructor");
  }

};

#include "NamespaceFooter.H"
#endif

---