11 #ifndef _LEVELDATAI_H_ 12 #define _LEVELDATAI_H_ 21 #include "NamespaceHeader.H" 47 : m_disjointBoxLayout(dp), m_ghost(ghost)
51 this->numReceives = 0;
64 MayDay::Error(
"non-disjoint DisjointBoxLayout: LevelData<T>::LevelData(const DisjointBoxLayout& dp, int comps)");
72 this->
setVector(*
this, interval, interval);
86 CH_TIME(
"LevelData<T>::define(dbl,comps,ghost,factory)");
96 MayDay::Error(
"non-disjoint DisjointBoxLayout: LevelData<T>::define(const DisjointBoxLayout& dp,....)");
100 MayDay::Error(
"LevelData::LevelData(const BoxLayout& dp, int comps) comps<=0");
129 CH_TIME(
"LevelData<T>::define(LevelData<T>,factory)");
136 if (
this == &da)
return;
160 CH_TIME(
"LevelData<T>::define(LevelData<T>,comps,factory)");
171 MayDay::Error(
" LevelData<T>::define(const LevelData<T>& da, const Interval& comps) called with 'this'");
213 int nbox = it.
size();
214 #pragma omp parallel for if(this->m_threadSafe) 215 for (
int ibox = 0; ibox < nbox; ibox++)
217 dest[it[ibox]].copy(this->
box(it[ibox]),
220 this->
operator[](it[ibox]),
227 copyTo(srcComps, dest, destComps, copier);
240 int nbox = it.
size();
241 #pragma omp parallel for if(this->m_threadSafe) 242 for (
int ibox = 0; ibox < nbox; ibox++)
244 Box srcBox = this->
box(it[ibox]);
245 Box dstBox = this->
box(it[ibox]);
250 dest[it[ibox]].copy(minBox,
253 this->
operator[](it[ibox]),
291 MayDay::Error(
"src == dest in copyTo function. Perhaps you want exchange ?");
298 int nbox = it.
size();
299 #pragma omp parallel for if(this->m_threadSafe) 300 for (
int ibox = 0; ibox < nbox; ibox++)
302 dest[it[ibox]].copy(this->
box(it[ibox]),
305 this->
operator[](it[ibox]),
312 copyTo(srcComps, dest, destComps, copier);
341 this->
makeItSo(srcComps, *
this, dest, destComps, copier, a_op);
371 this->
makeItSo(srcComps, *
this, dest, destComps, copier, a_op);
436 this->
makeItSo(comps, *
this, *
this, comps, copier, a_op);
482 MayDay::Error(
"LevelData<T>::define called with BoxLayout input");
490 MayDay::Error(
"LevelData<T>::define called with BoxLayout input");
498 MayDay::Error(
"LevelData<T>::define called with BoxLayout input");
507 MayDay::Error(
"LevelData<T>::define called with BoxLayout input");
516 int nbox = it.
size();
517 #pragma omp parallel for 518 for (
int ibox = 0; ibox < nbox; ibox++)
531 int nbox = it.
size();
532 #pragma omp parallel for 533 for (
int ibox = 0; ibox < nbox; ibox++)
542 template<
class T>
void 565 template<
class T>
void 587 for (toDit.
begin(); toDit.
ok(); ++toDit)
599 Box intersectBox = toFAB.
box();
600 intersectBox &= fromFAB.
box();
602 toFAB.
copy(fromFAB, intersectBox);
611 if (!fromDit.
ok()) done =
true;
620 #include "NamespaceFooter.H" int m_comps
Definition: BoxLayoutData.H:387
const ProblemDomain & physDomain() const
bool isDefined() const
Definition: Copier.H:306
Interval interval() const
Definition: BoxLayoutData.H:312
#define CH_SPACEDIM
Definition: SPACE.H:51
#define CH_assert(cond)
Definition: CHArray.H:37
Vector< T * > m_vector
Definition: LayoutData.H:124
void degenerate(LevelData< T > &a_to, const SliceSpec &a_ss) const
Definition: LevelDataI.H:543
LevelData()
Definition: LevelDataI.H:30
DisjointBoxLayout m_disjointBoxLayout
Definition: LevelData.H:297
void makeItSoLocalCopy(const Interval &a_srcComps, const BoxLayoutData< T > &a_src, BoxLayoutData< T > &a_dest, const Interval &a_destComps, const Copier &a_copier, const LDOperator< T > &a_op=LDOperator< T >()) const
Definition: BoxLayoutDataI.H:290
A not-necessarily-disjoint collective of boxes.
Definition: BoxLayout.H:145
int size() const
Definition: Interval.H:75
A strange but true thing to make copying from one boxlayoutdata to another fast.
Definition: Copier.H:145
int begin() const
Definition: Interval.H:97
Definition: SliceSpec.H:41
IntVect m_ghost
Definition: LevelData.H:299
const BoxLayout & boxLayout() const
Definition: LayoutData.H:107
void begin()
initialize this iterator to the first Box in its Layout
Definition: LayoutIterator.H:122
DataIterator dataIterator() const
Definition: LayoutDataI.H:78
void makeItSoBegin(const Interval &a_srcComps, const BoxLayoutData< T > &a_src, BoxLayoutData< T > &a_dest, const Interval &a_destComps, const Copier &a_copier, const LDOperator< T > &a_op=LDOperator< T >()) const
Definition: BoxLayoutDataI.H:249
void setVector(const BoxLayoutData< T > &da, const Interval &srcComps, const Interval &destComps)
Definition: BoxLayoutDataI.H:46
Box box(const DataIndex &a_index) const
Definition: LayoutDataI.H:66
virtual bool ok() const
return true if this iterator is still in its Layout
Definition: LayoutIterator.H:117
Definition: DataIterator.H:190
Copier m_exchangeCopier
Definition: LevelData.H:305
const Box & box() const
Definition: BaseFabImplem.H:226
BoxLayout m_boxLayout
Definition: LayoutData.H:118
virtual void exchange(void)
Simplest case – do all components.
Definition: LevelDataI.H:417
void degenerateLocalOnly(LevelData< T > &a_to, const SliceSpec &a_ss) const
version of degenerate which does strictly local copying
Definition: LevelDataI.H:566
Definition: EBInterface.H:45
virtual void exchangeBegin(const Copier &copier)
asynchronous exchange start. load and fire off messages.
Definition: LevelDataI.H:461
Box minBox(const Box &b1, const Box &b2)
virtual void apply(void(*a_Function)(const Box &, int, T &))
int size() const
Definition: DataIterator.H:218
void degenerate(DisjointBoxLayout &a_to, const SliceSpec &a_ss, bool a_maintainProcAssign=false) const
#define CH_TIME(name)
Definition: CH_Timer.H:82
Structure for passing component ranges in code.
Definition: Interval.H:23
virtual void localCopyTo(const Interval &srcComps, LevelData< T > &dest, const Interval &destComps) const
only works if source and dest have the same disjointboxlayout
Definition: LevelDataI.H:231
void allocateGhostVector(const DataFactory< T > &factory, const IntVect &ghost=IntVect::Zero)
Definition: BoxLayoutDataI.H:171
Data on a BoxLayout.
Definition: BoxLayoutData.H:97
virtual void define(const DisjointBoxLayout &dp, int comps, const IntVect &ghost=IntVect::Zero, const DataFactory< T > &a_factory=DefaultDataFactory< T >())
Definition: LevelDataI.H:83
BaseFab< T > & copy(const BaseFab< T > &a_src, const Box &a_srcbox, int a_srccomp, const Box &a_destbox, int a_destcomp, int a_numcomp)
Definition: BaseFabImplem.H:381
void makeItSo(const Interval &a_srcComps, const BoxLayoutData< T > &a_src, BoxLayoutData< T > &a_dest, const Interval &a_destComps, const Copier &a_copier, const LDOperator< T > &a_op=LDOperator< T >()) const
Definition: BoxLayoutDataI.H:236
virtual void copyTo(const Interval &srcComps, BoxLayoutData< T > &dest, const Interval &destComps) const
Definition: LevelDataI.H:201
A BoxLayout that has a concept of disjointedness.
Definition: DisjointBoxLayout.H:30
Box get(const LayoutIndex &it) const
Definition: BoxLayout.H:716
bool isPeriodic(int a_dir) const
Returns true if BC is periodic in direction a_dir.
Definition: ProblemDomain.H:894
const DisjointBoxLayout & getBoxes() const
Definition: LevelData.H:203
const IntVect & ghostVect() const
Definition: LevelData.H:170
static void Error(const char *const a_msg=m_nullString, int m_exitCode=CH_DEFAULT_ERROR_CODE)
Print out message to cerr and exit with the specified exit code.
virtual bool threadSafe() const
Definition: BoxLayoutData.H:49
static const IntVect Zero
Definition: IntVect.H:654
A Rectangular Domain on an Integer Lattice.
Definition: Box.H:465
void makeItSoEnd(BoxLayoutData< T > &a_dest, const Interval &a_destComps, const LDOperator< T > &a_op=LDOperator< T >()) const
Definition: BoxLayoutDataI.H:317
int nComp() const
Definition: BoxLayoutData.H:306
const DisjointBoxLayout & disjointBoxLayout() const
Definition: LevelData.H:209
const T & operator[](const DataIndex &a_index) const
const accessor function
Definition: LayoutDataI.H:25
bool m_isdefined
Definition: BoxLayoutData.H:389
virtual ~LevelData()
Definition: LevelDataI.H:37
An integer Vector in SpaceDim-dimensional space.
Definition: CHArray.H:42
Definition: FArrayBox.H:45
Factory object to data members of a BoxLayoutData container.
Definition: BoxLayoutData.H:30
Box & grow(int i)
grow functions
Definition: Box.H:2247
bool isClosed() const
Definition: BoxLayout.H:729
Definition: BoxLayoutData.H:173
bool m_threadSafe
Definition: BoxLayoutData.H:388
virtual void define(const DisjointBoxLayout &a_level, const BoxLayout &a_dest, bool a_exchange=false, IntVect a_shift=IntVect::Zero)
virtual void exchangeEnd()
finish asynchronous exchange
Definition: LevelDataI.H:471
int direction
Definition: SliceSpec.H:48
virtual void exchangeNoOverlap(const Copier &copier)
Definition: LevelDataI.H:450