BaseIFFAB.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

#ifndef _BASEIFFAB_H_
#define _BASEIFFAB_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 BaseIFFAB
{
public:

  BaseIFFAB();


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

  ~BaseIFFAB();


  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 BaseIFFAB<T>& a_src,
            const Interval& a_srcInterval);


  void clear();

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


  bool
  isDefined() const;


  int numFaces() const;


  int nComp() const;



  int direction() const;


  const IntVectSet& getIVS() const;


  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;

  const EBGraph& getEBGraph() 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);

  static void setVerbose(bool a_verbose);
  static bool s_verbose;
private:
  void
  setDefaultValues();

protected:

  static Arena* s_Arena;

  //get index vector<int> at the face location
  //using lexographic ordering of vof's indicies
  int getLocalVecIndex(const FaceIndex& a_face) const;

  T* m_data;
  int m_truesize;
  int m_nComp;
  int m_nFaces;
  int m_direction;

  //object to store pointers to facilitate fast indexing
  BaseFab<T*> m_fab;
  IntVectSet m_ivs;
  EBGraph m_ebgraph;
  bool m_isDefined;
private:
  //we disallow copy construction and assigment for data holders
  void operator= (const BaseIFFAB<T>&)
  {
    MayDay::Error("invalid operator");
  }
  BaseIFFAB (const BaseIFFAB<T>&)
  {
    MayDay::Error("invalid operator");
  }
};

#include "NamespaceFooter.H"

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

#endif

---