00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029 #ifndef _BaseEBCellFAB_H_
00030 #define _BaseEBCellFAB_H_
00031
00032 #include <cmath>
00033 #include <cstdlib>
00034 #include "Vector.H"
00035 #include "EBISBox.H"
00036 #include "BaseIVFAB.H"
00037 #include "IntVectSet.H"
00038 #include "BaseFab.H"
00039
00040
00042
00054 template <class T>
00055 class BaseEBCellFAB
00056 {
00057 public:
00059 BaseEBCellFAB();
00060
00062 BaseEBCellFAB(const EBISBox& a_ebisBox, const Box& a_region, int a_nVar);
00063
00065 void define(const EBISBox& a_ebisBox, const Box& a_region, int a_nVar);
00066
00068 virtual ~BaseEBCellFAB();
00069
00070
00072 void clear();
00073
00075 void
00076 setVal(T value);
00077
00079 bool isDefined() const;
00080
00081
00083 int nComp() const ;
00084
00086 const Box& getRegion() const ;
00087
00089 const Box& box() const ;
00090
00092 const IntVectSet& getMultiCells() const ;
00093
00095 const BaseIVFAB<T>& getMultiValuedFAB() const;
00096
00098 BaseIVFAB<T>& getMultiValuedFAB() ;
00099
00101 const BaseFab<T>& getSingleValuedFAB() const;
00102
00104 BaseFab<T>& getSingleValuedFAB();
00105
00107 const T& operator() (const VolIndex& a_ndin,int a_nVarLoc) const;
00108
00110 T& operator() (const VolIndex& a_ndin,int a_nVarLoc);
00111
00112
00114 void copy(const Box& RegionFrom,
00115 const Interval& destInt,
00116 const Box& RegionTo,
00117 const BaseEBCellFAB<T>& source,
00118 const Interval& srcInt);
00119
00121
00124 static int preAllocatable()
00125 {
00126 return 1;
00127 }
00128
00130 int size(const Box& R, const Interval& comps) const ;
00131
00133 int regSize(const Box& R, const Interval& comps) const ;
00134
00136 void linearOut(void* buf, const Box& R, const Interval& comps) const ;
00137
00139 void regLinearOut(void* buf, const Box& R, const Interval& comps) const ;
00140
00142 void linearIn(void* buf, const Box& R, const Interval& comps);
00143
00145 void regLinearIn(void* buf, const Box& R, const Interval& comps);
00146
00147
00149 BaseEBCellFAB(const Box& a_region, int a_nVar)
00150 {
00151 MayDay::Error("invalid constructor called for baseebcellfab");
00152 }
00153
00155 const EBISBox& getEBISBox() const;
00156
00158 void
00159 setCoveredCellVal(const T& a_val,
00160 const int& a_comp);
00161
00162 protected:
00163
00165 BaseIVFAB<T> m_irrFAB;
00166
00168 BaseFab<T> m_regFAB;
00169
00170 EBISBox m_ebisBox;
00171
00172 int m_nComp;
00173
00175 Box m_region;
00176
00178 BaseFab<bool> m_isMultiValued;
00179
00181 bool m_isDefined;
00182
00183 protected:
00184 void
00185 setDefaultValues();
00186
00187 private:
00188
00189 void operator= (const BaseEBCellFAB<T>& ebcin)
00190 {
00191 MayDay::Error("operator= for baseebcellfab undefined");
00192 }
00193 BaseEBCellFAB (const BaseEBCellFAB<T>& ebcin)
00194 {
00195 MayDay::Error("copy constructor for baseebcellfab undefined");
00196 }
00197 };
00198
00199 #include "BaseEBCellFABI.H"
00200
00201 #endif