BaseEBCellFAB.H

Go to the documentation of this file.

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

//  ANAG, LBNL

#ifndef _BASEEBCELLFAB_H_
#define _BASEEBCELLFAB_H_

#include <cmath>
#include <cstdlib>
#include "Vector.H"
#include "EBISBox.H"
#include "BaseIVFAB.H"
#include "IntVectSet.H"
#include "BaseFab.H"
#include "MiniIVFAB.H"
#include "NamespaceHeader.H"


template <class T>
class BaseEBCellFAB
{
public:
  BaseEBCellFAB();

  BaseEBCellFAB(const EBISBox& a_ebisBox, const Box& a_region, int a_nVar);

  void  define(const EBISBox& a_ebisBox, const Box& a_region, int a_nVar);


  BaseEBCellFAB(const Interval&   a_comps,
                BaseEBCellFAB<T>& a_original);

  virtual ~BaseEBCellFAB();

  void clear();

  void  setVal(const T& value);

  void  setVal(int ivar,const T& value);

  void  setVal(const T& value,
               const Box& a_box,
               int        a_nstart,
               int        a_numcomp);

  bool isDefined() const;

  int nComp() const ;

  const Box& getRegion() const ;

  const Box& box() const ;

  const IntVectSet& getMultiCells() const ;

  const BaseIVFAB<T>& getMultiValuedFAB() const;

  BaseIVFAB<T>& getMultiValuedFAB() ;

  const BaseFab<T>& getSingleValuedFAB() const;

  BaseFab<T>& getSingleValuedFAB();

  const T& operator() (const VolIndex& a_ndin, int a_nVarLoc,
                       int a_isKnownMultiValued=0) const;

  T& operator() (const VolIndex& a_ndin,int  a_nVarLoc,
                 int a_isKnownMultiValued=0);

  void fill(T* array, const VolIndex& a_ndin, const Interval& a_comps) const;

  void assign(const T* array, const VolIndex& a_ndin, const Interval& a_comps);
  void copy(const Box& RegionFrom,
            const Interval& destInt,
            const Box& RegionTo,
            const BaseEBCellFAB<T>& source,
            const Interval& srcInt);

  BaseEBCellFAB<T>& copy(const BaseEBCellFAB<T>& source);


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

  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);

  BaseEBCellFAB(const Box& a_region, int a_nVar)
  {
    MayDay::Error("invalid constructor called for baseebcellfab");
  }

  const EBISBox& getEBISBox() const;

  void
  setCoveredCellVal(const T&    a_val,
                    const int&  a_comp,
                    const bool& a_doMulti=true);

protected:

  //BaseIVFAB<T> m_irrFAB;
  MiniIVFAB<T> m_irrFAB;

  BaseFab<T> m_regFAB;

  EBISBox m_ebisBox;
  //number of data values at each VolIndex
  int m_nComp;

  Box m_region;

  IntVectSet m_multiCells;

  bool m_hasMultiCells;

  bool m_isDefined;

protected:
  void
  setDefaultValues();

private:
  //disallowed for perfomance reasons
  void operator= (const BaseEBCellFAB<T>& ebcin)
  {
    MayDay::Error("operator= for baseebcellfab undefined");
  }
  BaseEBCellFAB (const BaseEBCellFAB<T>& ebcin)
  {
    MayDay::Error("copy constructor for baseebcellfab undefined");
  }
};

#include "NamespaceFooter.H"

#ifndef CH_EXPLICIT_TEMPLATES
#include "BaseEBCellFABI.H"
#endif // CH_EXPLICIT_TEMPLATES

#endif

---