MiniIFFAB.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, ANAG, LBNL, 2009

#ifndef _MINIIFFAB_H_
#define _MINIIFFAB_H_

#include <cmath>
#include <cstdlib>
#include "SPACE.H"
#include "Vector.H"
#include "IntVectSet.H"
#include "FaceIndex.H"
#include "BaseFab.H"
#include "EBGraph.H"
#include "Arena.H"
#include "NamespaceHeader.H"



template <class T>
class MiniIFFAB
{
public:

  MiniIFFAB();


  MiniIFFAB(const IntVectSet& a_region,
            const EBGraph& a_ebgraph,
            const int& a_direction,
            const int& a_nvarin);

  ~MiniIFFAB();


  void define(const IntVectSet& a_region,
              const EBGraph& a_ebgraph,
              const int& a_direction,
              const int& a_nvarin);


  void setVal(const T& value);


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


  void copy(const Box& a_intBox,
            const Interval& a_destInterval,
            const Box& a_toBox,
            const MiniIFFAB<T>& a_src,
            const Interval& a_srcInterval);


  void clear();

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


  bool
  isDefined() const
  {
    return (m_isDefined);
  }



  int nComp() const
  {
    return m_nComp;
  }


  T& operator() (const FaceIndex& a_face,const int& varlocin);
  const T& operator()(const FaceIndex& a_face,const int& varlocin) const;


  T* dataPtr(const int& a_comp);

  T* getIndex(const FaceIndex& a_face, const int& a_comp) const;
  const T* dataPtr(const int& a_comp) const;

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

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

  int numFaces() const
  {
    return m_nFaces;
  }

  const IntVectSet& getIVS() const
  {
    return m_ivs;
  }
  
  static bool s_verbose;
private:
  void   setDefaultValues();
  
  void   getFaceSubset(Vector<FaceIndex>& a_subFaces, const Box& a_subBox) const;
protected:
  
  T* m_data;
  Vector<FaceIndex> m_faces;
  IntVectSet m_ivs;
  int m_nFaces;
  int m_nComp;
  int m_direction;

  EBGraph m_ebgraph;
  bool    m_isDefined;
private:
  //we disallow copy construction and assigment for data holders
  void operator= (const MiniIFFAB<T>&)
  {
    MayDay::Error("invalid operator");
  }
  MiniIFFAB (const MiniIFFAB<T>&)
  {
    MayDay::Error("invalid operator");
  }
};

#include "NamespaceFooter.H"

#ifndef CH_EXPLICIT_TEMPLATES
#include "MiniIFFABI.H"
#endif

#endif

---