00001 #ifdef CH_LANG_CC
00002
00003
00004
00005
00006
00007
00008
00009 #endif
00010
00011
00012
00013 #ifndef _EBISBOX_H_
00014 #define _EBISBOX_H_
00015
00016 #include "REAL.H"
00017 #include "IntVect.H"
00018 #include "IntVectSet.H"
00019 #include "VolIndex.H"
00020 #include "Interval.H"
00021 #include "FaceIndex.H"
00022 #include "Tuple.H"
00023 #include "IntVectSet.H"
00024 #include "RealVect.H"
00025 #include "BaseFab.H"
00026 #include "RefCountedPtr.H"
00027 #include "SPMD.H"
00028 #include "EBGraph.H"
00029 #include "EBData.H"
00030 #include "NamespaceHeader.H"
00031
00034
00046 class EBISBox
00047 {
00048 public:
00049
00051
00054 EBISBox();
00055
00057 ~EBISBox();
00058
00060
00064 IntVectSet getMultiCells(const Box& a_subbox) const;
00065
00067
00071 IntVectSet getIrregIVS(const Box& a_subbox) const;
00072
00074
00077 IntVectSet boundaryIVS(const Box& a_subbox) const;
00078
00080
00083 int numVoFs(const IntVect& a_iv) const;
00084
00086
00089 Vector<VolIndex> getVoFs(const IntVect& a_iv) const;
00090
00092
00103 Vector<VolIndex> getVoFs(const VolIndex& a_vof,
00104 const int& a_dir,
00105 const Side::LoHiSide& a_sd,
00106 const int& a_steps) const;
00107
00109
00112 bool isAllRegular() const;
00113
00115
00118 bool isAllCovered() const;
00119
00121
00124 bool isRegular(const IntVect& a_iv) const;
00125
00127
00130 bool isIrregular(const IntVect& a_iv) const;
00131
00133
00136 bool isCovered(const IntVect& a_iv) const;
00137
00139
00143 bool isCovered(const Box& a_box) const;
00144
00146
00150 bool isRegular(const Box& a_box) const;
00151
00153
00157 Vector<FaceIndex> getFaces(const VolIndex& a_vof,
00158 const int& a_idir,
00159 const Side::LoHiSide& a_sd) const;
00160
00162 Vector<FaceIndex> getAllFaces(const IntVect& a_iv,
00163 const int& a_idir,
00164 const Side::LoHiSide& a_sd) const;
00166
00176 RealVect centroid(const FaceIndex& facein) const;
00177
00179
00183 int numFaces(const VolIndex& a_vof,
00184 const int& a_idir,
00185 const Side::LoHiSide& a_sd) const;
00186
00188
00191 Real volFrac(const VolIndex& a_vof) const;
00192
00194
00197 bool isConnected(const VolIndex& a_vof1,
00198 const VolIndex& a_vof2) const;
00199
00201
00206 Real areaFrac(const FaceIndex& a_face1) const;
00207
00209
00213 Real sumArea(const VolIndex& a_vof,
00214 const int& a_idir,
00215 const Side::LoHiSide& a_sd) const;
00216
00218
00224 RealVect centroid(const VolIndex& a_vof) const;
00225
00227
00232 Vector<VolIndex> refine(const VolIndex& a_coarVoF) const;
00233
00235
00240 VolIndex coarsen(const VolIndex& a_fineVoF) const;
00241
00243
00245 RealVect bndryCentroid(const VolIndex& a_vof) const;
00246 RealVect bndryCentroid(const VolIndex& a_vof, int face) const;
00247
00249
00251 Real bndryArea(const VolIndex& a_vof) const;
00252 Real bndryArea(const VolIndex& a_vof, int face) const;
00253
00255 int numFacePhase(const VolIndex& a_vof) const ;
00256
00258 int facePhase(const VolIndex& a_vof, int face) const ;
00259
00261 const VolIndex& faceIndex(const VolIndex& a_vof, int face) const ;
00263
00268 RealVect normal(const VolIndex& a_vof) const;
00269 RealVect normal(const VolIndex& a_vof, int face) const;
00270
00272 const Box& getRegion() const;
00273
00275 const ProblemDomain& getDomain() const;
00276
00278 void setDomain(const ProblemDomain& a_domain);
00279
00281
00286 Vector<FaceIndex> refine(const FaceIndex& a_coarFace,const EBISBox& a_fineEBISBox) const;
00287
00289
00295 FaceIndex coarsen(const FaceIndex& a_fineFace) const;
00296
00298
00301 void setToAllRegular();
00302
00304
00307 void setToAllCovered();
00308
00310 void
00311 define(const BaseFab<int>& a_regIrregCovered,
00312 const Vector<IrregNode>& a_irregGraph,
00313 const Box& a_validRegion,
00314 const ProblemDomain& a_domain);
00315
00317
00319 const EBGraph& getEBGraph() const;
00320
00322
00324 const EBData& getEBData() const;
00325
00327 void
00328 define(const EBGraph& a_graph,
00329 const EBData& a_data);
00330
00332 EBISBox& operator=(const EBISBox& ebiin);
00333
00335 EBISBox(const EBISBox& ebiin);
00336
00338 bool operator==(const EBISBox& ebiin);
00339
00340 private:
00341
00343 EBGraph m_graph;
00344
00346 EBData m_data;
00347
00348 };
00349
00350
00351
00352
00353 inline bool
00354 EBISBox::isRegular(const IntVect& a_iv) const
00355 {
00356 return m_graph.isRegular(a_iv);
00357 }
00358
00359 inline bool
00360 EBISBox::isIrregular(const IntVect& a_iv) const
00361 {
00362 return m_graph.isIrregular(a_iv);
00363 }
00364
00365 inline bool
00366 EBISBox::isAllCovered() const
00367 {
00368 return m_graph.isAllCovered();
00369 }
00370
00371 inline bool
00372 EBISBox::isAllRegular() const
00373 {
00374 return m_graph.isAllRegular();
00375 }
00376
00377 inline bool
00378 EBISBox::isCovered(const IntVect& a_iv) const
00379 {
00380 return m_graph.isCovered(a_iv);
00381 }
00382
00383 inline bool
00384 EBISBox::isCovered(const Box& a_box) const
00385 {
00386 return m_graph.isCovered(a_box);
00387 }
00388
00389 inline bool
00390 EBISBox::isRegular(const Box& a_box) const
00391 {
00392 return m_graph.isRegular(a_box);
00393 }
00394
00395 inline Vector<VolIndex>
00396 EBISBox::getVoFs(const IntVect& a_iv) const
00397 {
00398 return m_graph.getVoFs(a_iv);
00399 }
00400
00401
00402 inline Vector<FaceIndex>
00403 EBISBox::getAllFaces(const IntVect& a_iv,
00404 const int& a_idir,
00405 const Side::LoHiSide& a_sd) const
00406 {
00407 return m_graph.getAllFaces(a_iv, a_idir, a_sd);
00408 }
00409
00410 inline Vector<FaceIndex>
00411 EBISBox::getFaces(const VolIndex& a_vof,
00412 const int& a_idir,
00413 const Side::LoHiSide& a_sd) const
00414 {
00415 return m_graph.getFaces(a_vof, a_idir, a_sd);
00416 }
00417
00418 inline int
00419 EBISBox::numVoFs(const IntVect& a_iv) const
00420 {
00421 return m_graph.numVoFs(a_iv);
00422 }
00423 #include "NamespaceFooter.H"
00424 #endif
