Main Page Namespace List Class Hierarchy Alphabetical List Compound List File List Compound Members File Members

ProblemDomain Class Reference

A class to facilitate interaction with physical boundary conditions. More...

#include <ProblemDomain.H>

Collaboration diagram for ProblemDomain:

Public Methods
	ProblemDomain ()
	The default constructor. The constructed domain box is empty.
	ProblemDomain (const Box &a_domBox)
	Construct ProblemDomain with a_domBox* as computational domain.*
	ProblemDomain (const Box &a_domBox, const bool *a_isPeriodic)
	Construct ProblemDomain with a_domBox as computational domain.
	ProblemDomain (const IntVect &small, const IntVect &big)
	Construct a ProblemDomain.
	ProblemDomain (const IntVect &small, const IntVect &big, const bool *a_isPeriodic)
	Construct a ProblemDomain.
	ProblemDomain (const IntVect &small, const int *vec_len)
	Construct ProblemDomain with specified lengths.
	ProblemDomain (const IntVect &small, const int vec_len, const bool a_isPeriodic)
	Construct ProblemDomain with specified lengths.
	ProblemDomain (const ProblemDomain &a_src)
	The copy constructor.
const Box &	domainBox () const
	Returns the logical computational domain.
bool	isPeriodic (int a_dir) const
	Returns true if BC is periodic in direction a_dir.
bool	isPeriodic () const
	Returns true is BC is periodic in _any_ direction.
ShiftIterator	shiftIterator () const
	Returns the shiftIterator for this ProblemDomain.
bool	isEmpty () const
	Returns true if this ProblemDomain is empty or undefined.
bool	contains (const IntVect &p) const
	Returns true if argument is contained within this ProblemDomain.
bool	image (IntVect &p) const
	Returns the periodic image of this IntVect inside of the ProblemDomain.
bool	contains (const Box &b) const
	Returns true if argument is contained within this ProblemDomain.
bool	contains_box (const Box &b) const
	Equivalent to contains()* function.*
bool	intersects (const Box &a_box) const
	Returns true if this ProblemDomain and the argument intersect.
bool	intersectsNotEmpty (const Box &a_box) const
	Returns true if this ProblemDomain and the argument intersect.
bool	intersects (const Box &box1, const Box &box2) const
	Returns true if box1 and box2 or any of their periodic images intersect.
ProblemDomain &	operator= (const ProblemDomain &b)
	The assignment operator.
void	setPeriodic (int a_dir, bool a_isPeriodic)
	Sets whether BC is periodic in direction a_dir (true == periodic).
ProblemDomain &	grow (int i)
	Grows (or shrinks) the domain Box by i in all directions.
ProblemDomain &	grow (const IntVect &v)
	Grow this ProblemDomain.
ProblemDomain &	grow (int idir, int n_cell)
	Grow this ProblemDomain.
ProblemDomain &	growLo (int idir, int n_cell=1)
	Grow this ProblemDomain on the low side.
ProblemDomain &	growHi (int idir, int n_cell=1)
	Grow this ProblemDomain on the high side.
Box	operator & (const Box &a_b) const
	Returns the Box intersection of this ProblemDomain and a_b.
ProblemDomain &	refine (int a_refinement_ratio)
	Refine this problem domain.
ProblemDomain &	refine (const IntVect &a_refinement_ratio)
	Refine this ProblemDomain.
ProblemDomain &	coarsen (int a_refinement_ratio)
	Coarsen this ProblemDomain.
ProblemDomain &	coarsen (const IntVect &refinement_ratio)
	Coarsen this ProblemDomain.
void	dumpOn (std::ostream &strm) const
	Gives more detail than printOn.
Protected Attributes
bool	m_isPeriodic [SpaceDim]
Box	m_domainBox
ShiftIterator	m_shiftIt
Friends
class	HDF5Handle
void	operator &= (Box &a_box, const ProblemDomain &a_probdomain)
	Modifies a_box to be the intersection of a_box and a_probdomain.
Box	operator & (const Box &a_box, const ProblemDomain &a_probdomain)
	Returns a Box which is the interesection of a_box and a_probdomain.
ProblemDomain	grow (const ProblemDomain &pd, int i)
	Returns a ProblemDomain with a domainBox grown by the given amount.
ProblemDomain	grow (const ProblemDomain &pd, const IntVect &v)
	Returns a grown version of pd.
Box	bdryLo (const ProblemDomain &a_pd, int a_dir, int a_len=1)
	Returns a face-centered Box at the low side of a_pd.
Box	bdryHi (const ProblemDomain &a_pd, int a_dir, int a_len=1)
	Returns a face-centered Box at the high side of a_pd.
Box	adjCellLo (const ProblemDomain &a_pd, int a_dir, int a_len=1)
	Returns the cell-centered Box adjacent to the low side of a_pd.
Box	adjCellHi (const ProblemDomain &a_pd, int a_dir, int a_len=1)
	Returns the cell-centered Box adjacent to the low side of a_pd.
ProblemDomain	refine (const ProblemDomain &a_probdomain, int a_refinement_ratio)
	Return a ProblemDomain which is a refinement of a_probdomain.
ProblemDomain	refine (const ProblemDomain &a_probdomain, const IntVect &a_refinement_ratio)
	Refinement function.
ProblemDomain	coarsen (const ProblemDomain &a_probdomain, int a_refinement_ratio)
	Coarsening function.
ProblemDomain	coarsen (const ProblemDomain &a_probdomain, const IntVect &a_refinement_ratio)
	Coarsening function.
std::ostream &	operator<< (std::ostream &os, const ProblemDomain &bx)
	Write an ASCII representation to the ostream.
std::istream &	operator>> (std::istream &is, ProblemDomain &bx)
	Read from istream.

Detailed Description

A class to facilitate interaction with physical boundary conditions.

ProblemDomain is a class which facilitates the application of physical boundary conditions, both periodic and non-periodic. This class contains much of the functionality of the Box class, since logically the computational domain is generally a Box.

Intersection with a ProblemDomain object will result in only removing regions which are outside the physical domain in non-periodic directions. Regions outside the logical computational domain in periodic directions will be treated as ghost cells which can be filled with an exchange() function or through suitable interpolation from a coarser domain.

Since ProblemDomain will contain a Box, it is a dimension dependent class, so SpaceDim must be defined as either 1, 2, or 3 when compiling.

Note that this implementation of ProblemDomain is inherently cell-centered.

Constructor & Destructor Documentation

ProblemDomain::ProblemDomain ( ) [inline]

The default constructor. The constructed domain box is empty.

ProblemDomain::ProblemDomain ( const Box & a_domBox )

Construct ProblemDomain with a_domBox as computational domain.
This constructor defaults to non-periodic domain

ProblemDomain::ProblemDomain ( const Box & a_domBox,

const bool * a_isPeriodic

)

Construct ProblemDomain with a_domBox as computational domain.
a_isPeriodic is a SpaceDim array of bools; if true, the physical boundary condition is periodic in the coordinate direction. False means a non-periodic BC.

ProblemDomain::ProblemDomain ( const IntVect & small,

const IntVect & big

)

Construct a ProblemDomain.
It is an error if small is greater than big. Defaults to non-periodic domain.

ProblemDomain::ProblemDomain ( const IntVect & small,

const IntVect & big,

const bool * a_isPeriodic

)

Construct a ProblemDomain.
It is an error if small is greater than big. a_isPeriodic is a SpaceDim array of bools; if true, the physical boundary condition is periodic in the coordinate direction. False means a non-periodic BC.

ProblemDomain::ProblemDomain ( const IntVect & small,

const int * vec_len

)

Construct ProblemDomain with specified lengths.
It is an error if the lengths are negative. Defaults to non-periodic domain.

ProblemDomain::ProblemDomain ( const IntVect & small,

const int * vec_len,

const bool * a_isPeriodic

)

Construct ProblemDomain with specified lengths.
It is an error if the lengths are negative. a_isPeriodic is a SpaceDim array of bools; if true, the physical boundary condition is periodic in the coordinate direction. False means a non-periodic BC.

ProblemDomain::ProblemDomain ( const ProblemDomain & a_src ) [inline]

The copy constructor.

Member Function Documentation

ProblemDomain& ProblemDomain::coarsen ( const IntVect & refinement_ratio )

Coarsen this ProblemDomain.
Modifies this ProblemDomain by coarsening by given (positive) refinement ratio. The Empty ProblemDomain is not modified by this function.

ProblemDomain& ProblemDomain::coarsen ( int a_refinement_ratio )

Coarsen this ProblemDomain.
Modifies this ProblemDomain by coarsening it by given (positive) refinement ratio. The Empty ProblemDomain is not modified by this function.

bool ProblemDomain::contains ( const Box & b ) const [inline]

Returns true if argument is contained within this ProblemDomain.
It is an error if the Box is not cell-centered. An empty ProblemDomain does not contain any Box. An entirely periodic domain contains any Box.

bool ProblemDomain::contains ( const IntVect & p ) const [inline]

Returns true if argument is contained within this ProblemDomain.
An empty ProblemDomain does not contain and is not contained by any ProblemDomain, including itself. Note also that in a periodic Direction, any index is valid, since the domain is infinite in that direction.

bool ProblemDomain::contains_box ( const Box & b ) const [inline]

Equivalent to contains() function.

const Box & ProblemDomain::domainBox ( ) const [inline]

Returns the logical computational domain.

void ProblemDomain::dumpOn ( std::ostream & strm ) const

Gives more detail than printOn.
Useful for exiting due to an error.

ProblemDomain & ProblemDomain::grow ( int idir,

int n_cell

) [inline]

Grow this ProblemDomain.
Modifies this ProblemDomain by growing it on the low and high end by n_cell cells in direction idir.

ProblemDomain & ProblemDomain::grow ( const IntVect & v ) [inline]

Grow this ProblemDomain.
Modifies this ProblemDomain by growing the domainBox in each direction by the specified amount

ProblemDomain & ProblemDomain::grow ( int i ) [inline]

Grows (or shrinks) the domain Box by i in all directions.

ProblemDomain & ProblemDomain::growHi ( int idir,

int n_cell = 1

) [inline]

Grow this ProblemDomain on the high side.
Modifies this ProblemDomain by growing it on the high end by n_Cell cells in direction idir

ProblemDomain & ProblemDomain::growLo ( int idir,

int n_cell = 1

) [inline]

Grow this ProblemDomain on the low side.
Modifies this ProblemDomain by growing it on the low end by n_cell cells in direction idir.

bool ProblemDomain::image ( IntVect & p ) const [inline]

Returns the periodic image of this IntVect inside of the ProblemDomain.
Return true if the domain contains this IntVect, returning the image in 'p', otherwise returns false

bool ProblemDomain::intersects ( const Box & box1,

const Box & box2

) const

Returns true if box1 and box2 or any of their periodic images intersect.
(useful for checking disjointness)

bool ProblemDomain::intersects ( const Box & a_box ) const

Returns true if this ProblemDomain and the argument intersect.
It is an error if a_box is not cell-centered. An empty ProblemDomain does not intersect any Box. This will do nothing in periodic directions (since a periodic domain is infinite in the periodic direction. If periodic in all dimensions, this will always return true.

bool ProblemDomain::intersectsNotEmpty ( const Box & a_box ) const

Returns true if this ProblemDomain and the argument intersect.
It is an error if a_box is not cell-centered. This routine does not perform the check to see if *this or b are empty Boxes. It is the callers responsibility to ensure that this never happens. If you are unsure, the use the .intersects(..) routine. In periodic directions, will always return true.

bool ProblemDomain::isEmpty ( ) const [inline]

Returns true if this ProblemDomain is empty or undefined.

bool ProblemDomain::isPeriodic ( ) const [inline]

Returns true is BC is periodic in _any_ direction.

bool ProblemDomain::isPeriodic ( int a_dir ) const [inline]

Returns true if BC is periodic in direction a_dir.

Box ProblemDomain::operator & ( const Box & a_b ) const

Returns the Box intersection of this ProblemDomain and a_b.
The Box MUST be cell-centered. The intersection of the Empty ProblemDomain and any Box is the Empty Box. This operator does nothing in periodic directions (since a periodic domain is an infinite domain).

ProblemDomain & ProblemDomain::operator= ( const ProblemDomain & b ) [inline]

The assignment operator.

ProblemDomain& ProblemDomain::refine ( const IntVect & a_refinement_ratio )

Refine this ProblemDomain.
Modifies this ProblemDomain by refining it by given (positive) refinement ratio. The Empty ProblemDomain is not modified by this function.

ProblemDomain& ProblemDomain::refine ( int a_refinement_ratio )

Refine this problem domain.
Modifies this ProblemDomain by refining it by given (positive) refinement ratio. The Empty ProblemDomain is not modified by this function.

void ProblemDomain::setPeriodic ( int a_dir,

bool a_isPeriodic

)

Sets whether BC is periodic in direction a_dir (true == periodic).

ShiftIterator ProblemDomain::shiftIterator ( ) const [inline]

Returns the shiftIterator for this ProblemDomain.
The shiftIterator is defined based on the periodicity of this ProblemDomain.

Friends And Related Function Documentation

Box adjCellHi ( const ProblemDomain & a_pd,

int a_dir,

int a_len = 1

) [friend]

Returns the cell-centered Box adjacent to the low side of a_pd.
Returns the cell centered Box of the given length adjacent to the argument ProblemDomain on the high end along the given coordinate direction. The return Box is identical to the argument ProblemDomain in the other directions. The return Box and the argument ProblemDomain have an empty intersection.
NOTES:

len >= 1.

Box retval = adjCellHi(b,dir,len) is equivalent to the following set of operations:
Box retval(b); retval.convert(dir,ProblemDomain::CELL); retval.setrange(dir,retval.bigEnd(dir)+1,len);
Directions are zero-based. It is an error if not 0 <= dir < SpaceDim. The neighbor of an Empty ProblemDomain is an Empty Box of the appropriate type. If a_dir is a periodic direction, will return an empty Box.

Box adjCellLo ( const ProblemDomain & a_pd,

int a_dir,

int a_len = 1

) [friend]

Returns the cell-centered Box adjacent to the low side of a_pd.
Returns the cell centered Box of the given length adjacent to the argument ProblemDomain on the low end along the given coordinate direction. The return Box is identical to the argument ProblemDomain in the other directions. The return ProblemDomain and the argument ProblemDomain have an empty intersection.
NOTES:

len >= 1.
Box retval = adjCellLo(b,dir,len) is equivalent to the following set of operations:
Box retval(b); retval.convert(dir,ProblemDomain::CELL); retval.setrange(dir,retval.smallEnd(dir)-len,len);
Directions are zero-based. It is an error if not 0 <= dir < SpaceDim. The neighbor of an Empty ProblemDomain is an Empty Box of the appropriate type. If a_dir is a periodic direction, will return an empty Box as well.

Box bdryHi ( const ProblemDomain & a_pd,

int a_dir,

int a_len = 1

) [friend]

Returns a face-centered Box at the high side of a_pd.
Returns the edge-centered Box (in direction a_dir) defining the high side of the argument ProblemDomain. The return Box will have the given length in the given direction. Directions are zero-based. It is an error if not 0 <= dir < SpaceDim. The neighbor of an Empty ProblemDomain is an Empty Box of the appropriate type. If a_dir is a periodic direction, will return an empty Box.

Box bdryLo ( const ProblemDomain & a_pd,

int a_dir,

int a_len = 1

) [friend]

Returns a face-centered Box at the low side of a_pd.
Returns the edge-centered Box (in direction a_dir) defining the low side of the argument ProblemDomain. The output Box will have the given length in the given direction. Directions are zero-based. It is an error if not 0 <= dir < SpaceDim. The neighbor of an Empty ProblemDomain is an Empty Box of the appropriate type. If a_dir is a periodic direction, will return an empty Box.

ProblemDomain coarsen ( const ProblemDomain & a_probdomain,

const IntVect & a_refinement_ratio

) [friend]

Coarsening function.
Returns a ProblemDomain that is the argument ProblemDomain coarsened by given (positive) refinement ratio. The Empty ProblemDomain is not modified by this function.

ProblemDomain coarsen ( const ProblemDomain & a_probdomain,

int a_refinement_ratio

) [friend]

Coarsening function.
Returns a ProblemDomain that is the argument ProblemDomain coarsened by given (positive) refinement ratio. The Empty ProblemDomain is not modified by this function.

ProblemDomain grow ( const ProblemDomain & pd,

const IntVect & v

) [friend]

Returns a grown version of pd.
Returns a ProblemDomain that is the argument ProblemDomain with a DomainBox grown by the given amount.

ProblemDomain grow ( const ProblemDomain & pd,

int i

) [friend]

Returns a ProblemDomain with a domainBox grown by the given amount.

friend class HDF5Handle [friend]

Box operator & ( const Box & a_box,

const ProblemDomain & a_probdomain

) [friend]

Returns a Box which is the interesection of a_box and a_probdomain.

void operator &= ( Box & a_box,

const ProblemDomain & a_probdomain

) [friend]

Modifies a_box to be the intersection of a_box and a_probdomain.

std::ostream& operator<< ( std::ostream & os,

const ProblemDomain & bx

) [friend]

Write an ASCII representation to the ostream.

std::istream& operator>> ( std::istream & is,

ProblemDomain & bx

) [friend]

Read from istream.

ProblemDomain refine ( const ProblemDomain & a_probdomain,

const IntVect & a_refinement_ratio

) [friend]

Refinement function.
Returns a ProblemDomain that is the argument ProblemDomain refined by given (positive) refinement ratio. The Empty ProblemDomain is not modified by this function.

ProblemDomain refine ( const ProblemDomain & a_probdomain,

int a_refinement_ratio

) [friend]

Return a ProblemDomain which is a refinement of a_probdomain.
Returns a ProblemDomain that is the argument ProblemDomain refined by given (positive) refinement ratio. The Empty ProblemDomain is not modified by this function.

Member Data Documentation

Box ProblemDomain::m_domainBox [protected]

Domain index extents

bool ProblemDomain::m_isPeriodic[SpaceDim] [protected]

Periodicity info

ShiftIterator ProblemDomain::m_shiftIt [protected]

Shift iterator for this ProblemDomain

The documentation for this class was generated from the following file:

ProblemDomain.H

Generated on Wed Jan 19 17:56:34 2005 for Chombo&INSwithParticles by

1.2.16

ProblemDomain Class Reference

Public Methods

Protected Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation