#include <ScdInterface.hpp>

Collaboration diagram for moab::ScdBox:

Public Member Functions
	~ScdBox ()
	Destructor. More...

ScdInterface *	sc_impl () const
	Return the ScdInterface responsible for this box. More...

ErrorCode	add_vbox (ScdBox *vbox, HomCoord from1, HomCoord to1, HomCoord from2, HomCoord to2, HomCoord from3, HomCoord to3, bool bb_input=false, const HomCoord &bb_min=HomCoord::getUnitv(0), const HomCoord &bb_max=HomCoord::getUnitv(0))
	Add a vertex box to this box. More...

bool	boundary_complete () const
	Return whether this box has all its vertices defined. More...

int	box_dimension () const
	Return highest topological dimension of box. More...

EntityHandle	start_vertex () const
	Starting vertex handle for this box. More...

EntityHandle	start_element () const
	Starting entity handle for this box. More...

int	num_elements () const
	Return the number of elements in the box. More...

int	num_vertices () const
	Return the number of vertices in the box. More...

const int *	box_dims () const
	Return the parametric coordinates for this box. More...

HomCoord	box_min () const
	Return the lower corner parametric coordinates for this box. More...

HomCoord	box_max () const
	Return the upper corner parametric coordinates for this box. More...

HomCoord	box_size () const
	Return the parameter extents for this box. More...

void	box_size (int *ijk) const
	Return the parametric extents for this box. More...

void	box_size (int &i, int &j, int &k) const
	Return the parametric extents for this box. More...

EntityHandle	get_element (const HomCoord &ijk) const
	Get the element at the specified coordinates. More...

EntityHandle	get_element (int i, int j=0, int k=0) const
	Get the element at the specified coordinates. More...

EntityHandle	get_vertex (const HomCoord &ijk) const
	Get the vertex at the specified coordinates. More...

EntityHandle	get_vertex (int i, int j=0, int k=0) const
	Get the vertex at the specified coordinates. More...

ErrorCode	get_params (EntityHandle ent, int &i, int &j, int &k, int &dir) const
	Get parametric coordinates of the specified entity. More...

ErrorCode	get_params (EntityHandle ent, int &i, int &j, int &k) const
	Get parametric coordinates of the specified entity, intermediate entities not allowed (no dir parameter) More...

ErrorCode	get_params (EntityHandle ent, HomCoord &ijkd) const
	Get parametric coordinates of the specified entity. More...

ErrorCode	get_adj_edge_or_face (int dim, int i, int j, int k, int dir, EntityHandle &ent, bool create_if_missing=true) const
	Get the adjacent edge or face at a parametric location This function gets the left (i=0), front (j=0), or bottom (k=0) edge or face for a parametric element. Left, front, or bottom is indicated by dir = 0, 1, or 2, resp. All edges and faces in a structured mesh block can be accessed using these parameters. More...

bool	contains (int i, int j, int k) const
	Return whether the box contains the parameters passed in. More...

bool	contains (const HomCoord &ijk) const
	Return whether the box contains the parameters passed in. More...

void	box_set (EntityHandle this_set)
	Set/Get the entity set representing the box. More...

EntityHandle	box_set ()

ErrorCode	get_coordinate_arrays (double &xc, double &yc, double *&zc)
	Get coordinate arrays for vertex coordinates for a structured block. More...

ErrorCode	get_coordinate_arrays (const double &xc, const double &yc, const double *&zc) const
	Get read-only coordinate arrays for vertex coordinates for a structured block. More...

bool	locally_periodic_i () const
	Return whether box is locally periodic in i. More...

bool	locally_periodic_j () const
	Return whether box is locally periodic in j. More...

bool	locally_periodic_k () const
	Return whether box is locally periodic in k. More...

void	locally_periodic (bool lperiodic[3])
	Set local periodicity. More...

const int *	locally_periodic () const
	Get local periodicity. More...

ScdParData &	par_data ()
	Return parallel data. More...

const ScdParData &	par_data () const
	Return parallel data. More...

void	par_data (const ScdParData &par_datap)
	set parallel data More...

Private Member Functions
	ScdBox (ScdInterface sc_impl, EntityHandle box_set, EntitySequence seq1, EntitySequence *seq2=NULL)
	Constructor. More...

EntityHandle	get_vertex_from_seq (int i, int j, int k) const
	function to get vertex handle directly from sequence More...

ErrorCode	vert_dat (ScdVertexData *vert_dat)
	set the vertex sequence More...

ScdVertexData *	vert_dat () const
	get the vertex sequence More...

ErrorCode	elem_seq (EntitySequence *elem_seq)
	set the element sequence More...

StructuredElementSeq *	elem_seq () const
	get the element sequence More...

void	start_vertex (EntityHandle startv)
	Set the starting vertex handle for this box. More...

void	start_element (EntityHandle starte)
	Set the starting entity handle for this box. More...

Private Attributes
ScdInterface *	scImpl
	interface instance More...

EntityHandle	boxSet
	entity set representing this box More...

ScdVertexData *	vertDat
	vertex sequence this box represents, if there's only one, otherwise they're retrieved from the element sequence More...

StructuredElementSeq *	elemSeq
	element sequence this box represents More...

EntityHandle	startVertex
	starting vertex handle for this box More...

EntityHandle	startElem
	starting element handle for this box More...

int	boxDims [6]
	lower and upper corners More...

int	locallyPeriodic [3]
	is locally periodic in i or j or k More...

ScdParData	parData
	parallel data associated with this box, if any More...

HomCoord	boxSize
	parameter extents More...

int	boxSizeIJ
	convenience parameters, (boxSize[1]-1)*(boxSize[0]-1) and boxSize[0]-1 More...

int	boxSizeIJM1

int	boxSizeIM1

Friends
class	ScdInterface

Detailed Description

Definition at line 506 of file ScdInterface.hpp.

Constructor & Destructor Documentation

◆ ~ScdBox()

moab::ScdBox::~ScdBox ( )

Destructor.

Definition at line 548 of file ScdInterface.cpp.

 {
     // Reset the tag on the set
     if( boxSet )
     {
         // It is possible that the box set entity has been deleted (e.g. by
         // Core::clean_up_failed_read)
         Core* mbcore = dynamic_cast< Core* >( scImpl->mbImpl );
         assert( mbcore != NULL );
         if( mbcore->is_valid( boxSet ) )
         {
             ScdBox* tmp_ptr = NULL;
             scImpl->mbImpl->tag_set_data( scImpl->box_set_tag(), &boxSet, 1, &tmp_ptr );
         }
         else
             boxSet = 0;
     }
  
     scImpl->remove_box( this );
 }

References moab::ScdInterface::box_set_tag(), boxSet, moab::Core::is_valid(), moab::ScdInterface::mbImpl, moab::ScdInterface::remove_box(), scImpl, and moab::Interface::tag_set_data().

◆ ScdBox()

moab::ScdBox::ScdBox	(	ScdInterface *	sc_impl,
		EntityHandle	box_set,
		EntitySequence *	seq1,
		EntitySequence *	seq2 = `NULL`
	)

private

Constructor.

Create a structured box instance; this constructor is private because it should only be called from ScdInterface, a friend class. This constructor takes two sequences, one of which can be NULL. If both sequences come in non-NULL, the first should be a VertexSequence* corresponding to a structured vertex sequence and the second should be a StructuredElementSeq*. If the 2nd is NULL, the first can be either of those types. The other members of this class are taken from the sequences (e.g. parametric space) or the box set argument. Tags on the box set should be set from the caller.

Parameters

sc_impl	A ScdInterface instance
box_set	Entity set representing this rectangle
seq1	An EntitySequence (see ScdBox description)
seq2	An EntitySequence (see ScdBox description), or NULL

Definition at line 474 of file ScdInterface.cpp.

     : scImpl( impl ), boxSet( bset ), vertDat( NULL ), elemSeq( NULL ), startVertex( 0 ), startElem( 0 )
 {
     for( int i = 0; i < 6; i++ )
         boxDims[i] = 0;
     for( int i = 0; i < 3; i++ )
         locallyPeriodic[i] = false;
     VertexSequence* vseq = dynamic_cast< VertexSequence* >( seq1 );
     if( vseq ) vertDat = dynamic_cast< ScdVertexData* >( vseq->data() );
     if( vertDat )
     {
         // retrieve the parametric space
         for( int i = 0; i < 3; i++ )
         {
             boxDims[i]     = vertDat->min_params()[i];
             boxDims[3 + i] = vertDat->max_params()[i];
         }
         startVertex = vertDat->start_handle();
     }
     else if( impl->boxDimsTag )
     {
         // look for parametric space info on set
         ErrorCode rval = impl->mbImpl->tag_get_data( impl->boxDimsTag, &bset, 1, boxDims );
         if( MB_SUCCESS == rval )
         {
             Range verts;
             impl->mbImpl->get_entities_by_dimension( bset, 0, verts );
             if( !verts.empty() ) startVertex = *verts.begin();
         }
     }
  
     elemSeq = dynamic_cast< StructuredElementSeq* >( seq2 );
     if( !elemSeq ) elemSeq = dynamic_cast< StructuredElementSeq* >( seq1 );
  
     if( elemSeq )
     {
         if( vertDat )
         {
             // check the parametric space to make sure it's consistent
             assert( elemSeq->sdata()->min_params() == HomCoord( boxDims, 3 ) &&
                     ( elemSeq->sdata()->max_params() + HomCoord( 1, 1, 1 ) ) == HomCoord( boxDims, 3 ) );
         }
         else
         {
             // get the parametric space from the element sequence
             for( int i = 0; i < 3; i++ )
             {
                 boxDims[i]     = elemSeq->sdata()->min_params()[i];
                 boxDims[3 + i] = elemSeq->sdata()->max_params()[i];
             }
         }
  
         startElem = elemSeq->start_handle();
     }
     else
     {
         Range elems;
         impl->mbImpl->get_entities_by_dimension(
             bset, ( boxDims[2] == boxDims[5] ? ( boxDims[1] == boxDims[4] ? 1 : 2 ) : 3 ), elems );
         if( !elems.empty() ) startElem = *elems.begin();
         // call the following w/o looking at return value, since it doesn't really need to be there
         if( impl->boxPeriodicTag ) impl->mbImpl->tag_get_data( impl->boxPeriodicTag, &bset, 1, locallyPeriodic );
     }
  
     assert( vertDat || elemSeq || boxDims[0] != boxDims[3] || boxDims[1] != boxDims[4] || boxDims[2] != boxDims[5] );
  
     boxSize     = HomCoord( boxDims + 3, 3 ) - HomCoord( boxDims, 3 ) + HomCoord( 1, 1, 1 );
     boxSizeIJ   = ( boxSize[1] ? boxSize[1] : 1 ) * boxSize[0];
     boxSizeIM1  = boxSize[0] - ( locallyPeriodic[0] ? 0 : 1 );
     boxSizeIJM1 = ( boxSize[1] ? ( boxSize[1] - ( locallyPeriodic[1] ? 0 : 1 ) ) : 1 ) * boxSizeIM1;
  
     scImpl->add_box( this );
 }

References moab::ScdInterface::add_box(), moab::Range::begin(), boxDims, moab::ScdInterface::boxDimsTag, moab::ScdInterface::boxPeriodicTag, boxSize, boxSizeIJ, boxSizeIJM1, boxSizeIM1, moab::EntitySequence::data(), elemSeq, moab::Range::empty(), ErrorCode, moab::Interface::get_entities_by_dimension(), locallyPeriodic, moab::ScdElementData::max_params(), moab::ScdVertexData::max_params(), MB_SUCCESS, moab::ScdInterface::mbImpl, moab::ScdElementData::min_params(), moab::ScdVertexData::min_params(), scImpl, moab::StructuredElementSeq::sdata(), moab::EntitySequence::start_handle(), moab::SequenceData::start_handle(), startElem, startVertex, moab::Interface::tag_get_data(), and vertDat.

Member Function Documentation

◆ add_vbox()

ErrorCode moab::ScdBox::add_vbox	(	ScdBox *	vbox,
		HomCoord	from1,
		HomCoord	to1,
		HomCoord	from2,
		HomCoord	to2,
		HomCoord	from3,
		HomCoord	to3,
		bool	bb_input = `false`,
		const HomCoord &	bb_min = `HomCoord::getUnitv( 0 )`,
		const HomCoord &	bb_max = `HomCoord::getUnitv( 0 )`
	)

Add a vertex box to this box.

Definition at line 580 of file ScdInterface.cpp.

 {
     if( !vbox->vertDat ) return MB_FAILURE;
     ScdVertexData* dum_data = dynamic_cast< ScdVertexData* >( vbox->vertDat );
     ErrorCode rval =
         elemSeq->sdata()->add_vsequence( dum_data, from1, to1, from2, to2, from3, to3, bb_input, bb_min, bb_max );
     return rval;
 }

References moab::ScdElementData::add_vsequence(), elemSeq, ErrorCode, moab::StructuredElementSeq::sdata(), and vertDat.

Referenced by moab::ScdInterface::construct_box().

◆ boundary_complete()

bool moab::ScdBox::boundary_complete ( ) const

Return whether this box has all its vertices defined.

Tests whether vertex boxs added with add_vbox have completely defined the vertex parametric space for this box.

Definition at line 598 of file ScdInterface.cpp.

 {
     return elemSeq->boundary_complete();
 }

References moab::StructuredElementSeq::boundary_complete(), and elemSeq.

Referenced by moab::ScdNCHelper::create_mesh().

◆ box_dimension()

int moab::ScdBox::box_dimension ( ) const

inline

Return highest topological dimension of box.

Definition at line 575 of file ScdInterface.cpp.

 {
     return ( startElem ? scImpl->mbImpl->dimension_from_handle( startElem ) : -1 );
 }

References moab::Interface::dimension_from_handle(), moab::ScdInterface::mbImpl, scImpl, and startElem.

Referenced by get_params(), and moab::ScdInterface::tag_shared_vertices().

◆ box_dims()

const int * moab::ScdBox::box_dims ( ) const

inline

Return the parametric coordinates for this box.

Returns: IJK parameters of lower and upper corners

Definition at line 1520 of file ScdInterface.hpp.

 {
     return boxDims;
 }

References boxDims.

Referenced by moab::ScdInterface::assign_global_ids(), and moab::ScdInterface::get_shared_vertices().

◆ box_max()

HomCoord moab::ScdBox::box_max ( ) const

inline

Return the upper corner parametric coordinates for this box.

Definition at line 1530 of file ScdInterface.hpp.

 {
     return HomCoord( boxDims + 3, 3 );
 }

References boxDims.

Referenced by moab::Skinner::skin_box().

◆ box_min()

HomCoord moab::ScdBox::box_min ( ) const

inline

Return the lower corner parametric coordinates for this box.

Definition at line 1525 of file ScdInterface.hpp.

 {
     return HomCoord( boxDims, 3 );
 }

References boxDims.

Referenced by moab::Skinner::skin_box().

◆ box_set() [1/2]

EntityHandle moab::ScdBox::box_set ( )

inline

Definition at line 1594 of file ScdInterface.hpp.

 {
     return boxSet;
 }

References boxSet.

◆ box_set() [2/2]

void moab::ScdBox::box_set ( EntityHandle this_set )

inline

Set/Get the entity set representing the box.

Definition at line 1589 of file ScdInterface.hpp.

 {
     boxSet = this_set;
 }

References boxSet.

Referenced by moab::ScdInterface::construct_box(), moab::ScdNCHelper::create_mesh(), and moab::ScdInterface::tag_shared_vertices().

◆ box_size() [1/3]

HomCoord moab::ScdBox::box_size ( ) const

inline

Return the parameter extents for this box.

Definition at line 1535 of file ScdInterface.hpp.

 {
     return boxSize;
 }

References boxSize.

Referenced by moab::ScdInterface::construct_box().

◆ box_size() [2/3]

void moab::ScdBox::box_size	(	int &	i,
		int &	j,
		int &	k
	)		const

inline

Return the parametric extents for this box.

Parameters

i	I extent of this box
j	J extent of this box
k	K extent of this box

Definition at line 1547 of file ScdInterface.hpp.

 {
     i = boxSize[0];
     j = boxSize[1];
     k = boxSize[2];
 }

References boxSize.

◆ box_size() [3/3]

void moab::ScdBox::box_size ( int * ijk ) const

inline

Return the parametric extents for this box.

Parameters

ijk	IJK extents of this box

Definition at line 1540 of file ScdInterface.hpp.

 {
     ijk[0] = boxSize[0];
     ijk[1] = boxSize[1];
     ijk[2] = boxSize[2];
 }

References boxSize.

◆ contains() [1/2]

bool moab::ScdBox::contains ( const HomCoord & ijk ) const

inline

Return whether the box contains the parameters passed in.

Parameters

i	Parametric coordinates being evaluated
j	Parametric coordinates being evaluated
k	Parametric coordinates being evaluated

Definition at line 1579 of file ScdInterface.hpp.

 {
     return ( ijk >= HomCoord( boxDims, 3 ) && ijk <= HomCoord( boxDims + 3, 3 ) );
 }

References boxDims.

◆ contains() [2/2]

bool moab::ScdBox::contains	(	int	i,
		int	j,
		int	k
	)		const

inline

Return whether the box contains the parameters passed in.

Parameters

i	Parametric coordinates being evaluated
j	Parametric coordinates being evaluated
k	Parametric coordinates being evaluated

Definition at line 1584 of file ScdInterface.hpp.

 {
     return contains( HomCoord( i, j, k ) );
 }

◆ elem_seq() [1/2]

StructuredElementSeq * moab::ScdBox::elem_seq ( ) const

inlineprivate

get the element sequence

Definition at line 1604 of file ScdInterface.hpp.

 {
     return elemSeq;
 }

References elemSeq.

◆ elem_seq() [2/2]

ErrorCode moab::ScdBox::elem_seq ( EntitySequence * elem_seq )

private

set the element sequence

Definition at line 628 of file ScdInterface.cpp.

 {
     elemSeq = dynamic_cast< StructuredElementSeq* >( elem_sq );
     if( elemSeq ) elemSeq->is_periodic( locallyPeriodic );
  
     if( locallyPeriodic[0] ) boxSizeIM1 = boxSize[0] - ( locallyPeriodic[0] ? 0 : 1 );
     if( locallyPeriodic[0] || locallyPeriodic[1] )
         boxSizeIJM1 = ( boxSize[1] ? ( boxSize[1] - ( locallyPeriodic[1] ? 0 : 1 ) ) : 1 ) * boxSizeIM1;
  
     return ( elemSeq ? MB_SUCCESS : MB_FAILURE );
 }

References boxSize, boxSizeIJM1, boxSizeIM1, elemSeq, moab::StructuredElementSeq::is_periodic(), locallyPeriodic, and MB_SUCCESS.

Referenced by moab::ScdInterface::construct_box().

◆ get_adj_edge_or_face()

ErrorCode moab::ScdBox::get_adj_edge_or_face	(	int	dim,
		int	i,
		int	j,
		int	k,
		int	dir,
		EntityHandle &	ent,
		bool	create_if_missing = `true`
	)		const

Get the adjacent edge or face at a parametric location This function gets the left (i=0), front (j=0), or bottom (k=0) edge or face for a parametric element. Left, front, or bottom is indicated by dir = 0, 1, or 2, resp. All edges and faces in a structured mesh block can be accessed using these parameters.

Get the entity of specified dimension adjacent to parametric element.

Parameters

dim	Dimension of adjacent entity being requested
i	Parametric coordinates of cell being evaluated
j	Parametric coordinates of cell being evaluated
k	Parametric coordinates of cell being evaluated
dir	Direction (0, 1, or 2), for getting adjacent edges (2d, 3d) or faces (3d)
ent	Entity returned from this function
create_if_missing	If true, creates the entity if it doesn't already exist
dim	Dimension of adjacent entity being requested
i	Parametric coordinates of cell being evaluated
j	Parametric coordinates of cell being evaluated
k	Parametric coordinates of cell being evaluated
dir	Direction (0, 1, or 2), for getting adjacent edges (2d, 3d) or faces (3d)
ent	EntityHandle of adjacent entity
create_if_missing	If true, creates the entity if it doesn't already exist

Definition at line 669 of file ScdInterface.cpp.

 {
     // describe connectivity of sub-element in static array
     // subconnect[dim-1][dir][numv][ijk] where dimensions are:
     // [dim-1]: dim=1 or 2, so this is 0 or 1
     // [dir]: one of 0..2, for ijk directions in a hex
     // [numv]: number of vertices describing sub entity = 2*dim <= 4
     // [ijk]: 3 values for i, j, k
     int subconnect[2][3][4][3] = { { { { 0, 0, 0 }, { 1, 0, 0 }, { -1, -1, -1 }, { -1, -1, -1 } },    // i edge
                                      { { 0, 0, 0 }, { 0, 1, 0 }, { -1, -1, -1 }, { -1, -1, -1 } },    // j edge
                                      { { 0, 0, 0 }, { 0, 0, 1 }, { -1, -1, -1 }, { -1, -1, -1 } } },  // k edge
  
                                    { { { 0, 0, 0 }, { 0, 1, 0 }, { 0, 1, 1 }, { 0, 0, 1 } },      // i face
                                      { { 0, 0, 0 }, { 1, 0, 0 }, { 1, 0, 1 }, { 0, 0, 1 } },      // j face
                                      { { 0, 0, 0 }, { 1, 0, 0 }, { 1, 1, 0 }, { 0, 1, 0 } } } };  // k face
  
     // check proper input dimensions and lower bound
     if( dim < 1 || dim > 2 || i < boxDims[0] || j < boxDims[1] || k < boxDims[2] ) return MB_FAILURE;
  
     // now check upper bound; parameters must be <= upper corner, since edges/faces
     // follow element parameterization, not vertex parameterization
     else if( ( boxDims[3] != boxDims[0] && i > ( locallyPeriodic[0] ? boxDims[3] + 1 : boxDims[3] ) ) ||
              ( boxDims[4] != boxDims[1] && j > ( locallyPeriodic[1] ? boxDims[4] + 1 : boxDims[4] ) ) ||
              ( boxDims[5] != boxDims[2] && k > boxDims[5] ) )
         return MB_FAILURE;
  
     // get the vertices making up this entity
     EntityHandle verts[4];
     for( int ind = 0; ind < 2 * dim; ind++ )
     {
         int i1 = i + subconnect[dim - 1][dir][ind][0];
         int j1 = j + subconnect[dim - 1][dir][ind][1];
         // if periodic in i and i1 is boxDims[3]+1, wrap around
         if( locallyPeriodic[0] && i1 == boxDims[3] + 1 ) i1 = boxDims[0];
         // if periodic in j and j1 is boxDims[4]+1, wrap around
         if( locallyPeriodic[1] && j1 == boxDims[4] + 1 ) j1 = boxDims[1];
         verts[ind] = get_vertex( i1, j1, k + subconnect[dim - 1][dir][ind][2] );
         if( !verts[ind] ) return MB_FAILURE;
     }
  
     Range ents;
     ErrorCode rval = scImpl->impl()->get_adjacencies( verts, 2 * dim, dim, false, ents );
     if( MB_SUCCESS != rval ) return rval;
  
     if( ents.size() > 1 )
         return MB_FAILURE;
  
     else if( ents.size() == 1 )
     {
         ent = *ents.begin();
     }
     else if( create_if_missing )
         rval = scImpl->impl()->create_element( ( 1 == dim ? MBEDGE : MBQUAD ), verts, 2 * dim, ent );
  
     return rval;
 }

References moab::Range::begin(), boxDims, moab::Interface::create_element(), dim, ErrorCode, moab::Interface::get_adjacencies(), get_vertex(), moab::ScdInterface::impl(), locallyPeriodic, MB_SUCCESS, MBEDGE, MBQUAD, scImpl, and moab::Range::size().

Referenced by moab::Skinner::skin_box().

◆ get_coordinate_arrays() [1/2]

ErrorCode moab::ScdBox::get_coordinate_arrays	(	const double *&	xc,
		const double *&	yc,
		const double *&	zc
	)		const

Get read-only coordinate arrays for vertex coordinates for a structured block.

Returns error if there isn't a single vertex sequence associated with this structured block

Parameters

xc	X coordinate array pointer returned
yc	Y coordinate array pointer returned
zc	Z coordinate array pointer returned

Definition at line 613 of file ScdInterface.cpp.

 {
     if( !vertDat ) return MB_FAILURE;
     xc = reinterpret_cast< const double* >( vertDat->get_sequence_data( 0 ) );
     yc = reinterpret_cast< const double* >( vertDat->get_sequence_data( 1 ) );
     zc = reinterpret_cast< const double* >( vertDat->get_sequence_data( 2 ) );
     return MB_SUCCESS;
 }

References moab::SequenceData::get_sequence_data(), MB_SUCCESS, and vertDat.

◆ get_coordinate_arrays() [2/2]

ErrorCode moab::ScdBox::get_coordinate_arrays	(	double *&	xc,
		double *&	yc,
		double *&	zc
	)

Get coordinate arrays for vertex coordinates for a structured block.

Returns error if there isn't a single vertex sequence associated with this structured block

Parameters

xc	X coordinate array pointer returned
yc	Y coordinate array pointer returned
zc	Z coordinate array pointer returned

Definition at line 603 of file ScdInterface.cpp.

 {
     if( !vertDat ) return MB_FAILURE;
  
     xc = reinterpret_cast< double* >( vertDat->get_sequence_data( 0 ) );
     yc = reinterpret_cast< double* >( vertDat->get_sequence_data( 1 ) );
     zc = reinterpret_cast< double* >( vertDat->get_sequence_data( 2 ) );
     return MB_SUCCESS;
 }

References moab::SequenceData::get_sequence_data(), MB_SUCCESS, and vertDat.

Referenced by moab::ScdInterface::construct_box(), and moab::ScdNCHelper::create_mesh().

◆ get_element() [1/2]

EntityHandle moab::ScdBox::get_element ( const HomCoord & ijk ) const

inline

Get the element at the specified coordinates.

Parameters

ijk	Parametric coordinates being evaluated

Definition at line 1561 of file ScdInterface.hpp.

 {
     return get_element( ijk[0], ijk[1], ijk[2] );
 }

Referenced by moab::Core::create_scd_sequence(), and main().

◆ get_element() [2/2]

EntityHandle moab::ScdBox::get_element	(	int	i,
		int	j = `0`,
		int	k = `0`
	)		const

inline

Get the element at the specified coordinates.

Parameters

i	Parametric coordinates being evaluated
j	Parametric coordinates being evaluated
k	Parametric coordinates being evaluated

Definition at line 1554 of file ScdInterface.hpp.

 {
     return ( !startElem
                  ? 0
                  : startElem + ( k - boxDims[2] ) * boxSizeIJM1 + ( j - boxDims[1] ) * boxSizeIM1 + i - boxDims[0] );
 }

References boxDims, boxSizeIJM1, boxSizeIM1, and startElem.

◆ get_params() [1/3]

ErrorCode moab::ScdBox::get_params	(	EntityHandle	ent,
		HomCoord &	ijkd
	)		const

Get parametric coordinates of the specified entity.

This function returns MB_ENTITY_NOT_FOUND if the entity is not in this ScdBox.

Parameters

ent	Entity being queried
ijkd	Parametric coordinates returned (including direction, in case of getting adjacent edges (2d, 3d) or faces (3d))

Definition at line 640 of file ScdInterface.cpp.

 {
     // check first whether this is an intermediate entity, so we know what to do
     int dimension = box_dimension();
     int this_dim  = scImpl->impl()->dimension_from_handle( ent );
  
     if( ( 0 == this_dim && !vertDat ) || ( this_dim && this_dim == dimension ) )
     {
         assert( elemSeq );
         return elemSeq->get_params( ent, ijkd[0], ijkd[1], ijkd[2] );
     }
  
     else if( !this_dim && vertDat )
         return vertDat->get_params( ent, ijkd[0], ijkd[1], ijkd[2] );
  
     else
         return MB_NOT_IMPLEMENTED;
 }

References box_dimension(), moab::Interface::dimension_from_handle(), elemSeq, moab::StructuredElementSeq::get_params(), moab::ScdVertexData::get_params(), moab::ScdInterface::impl(), MB_NOT_IMPLEMENTED, scImpl, and vertDat.

◆ get_params() [2/3]

ErrorCode moab::ScdBox::get_params	(	EntityHandle	ent,
		int &	i,
		int &	j,
		int &	k
	)		const

inline

Get parametric coordinates of the specified entity, intermediate entities not allowed (no dir parameter)

This function returns MB_ENTITY_NOT_FOUND if the entity is not in this ScdBox, or MB_FAILURE if the entity is an intermediate-dimension entity.

Parameters

ent	Entity being queried
i	Parametric coordinates returned
j	Parametric coordinates returned
k	Parametric coordinates returned

Definition at line 1624 of file ScdInterface.hpp.

 {
     HomCoord hc;
     ErrorCode rval = get_params( ent, hc );
     if( MB_SUCCESS == rval )
     {
         i = hc[0];
         j = hc[1];
         k = hc[2];
     }
  
     return rval;
 }

References ErrorCode, get_params(), and MB_SUCCESS.

◆ get_params() [3/3]

ErrorCode moab::ScdBox::get_params	(	EntityHandle	ent,
		int &	i,
		int &	j,
		int &	k,
		int &	dir
	)		const

inline

Get parametric coordinates of the specified entity.

This function returns MB_ENTITY_NOT_FOUND if the entity is not in this ScdBox.

Parameters

ent	Entity being queried
i	Parametric coordinates returned
j	Parametric coordinates returned
k	Parametric coordinates returned
dir	Parametric coordinate direction returned (in case of getting adjacent edges (2d, 3d) or faces (3d); not modified otherwise

Definition at line 1609 of file ScdInterface.hpp.

 {
     HomCoord hc;
     ErrorCode rval = get_params( ent, hc );
     if( MB_SUCCESS == rval )
     {
         i   = hc[0];
         j   = hc[1];
         k   = hc[2];
         dir = hc[3];
     }
  
     return rval;
 }

References ErrorCode, and MB_SUCCESS.

Referenced by get_params().

◆ get_vertex() [1/2]

EntityHandle moab::ScdBox::get_vertex ( const HomCoord & ijk ) const

inline

Get the vertex at the specified coordinates.

Parameters

ijk	Parametric coordinates being evaluated

Definition at line 1574 of file ScdInterface.hpp.

 {
     return get_vertex( ijk[0], ijk[1], ijk[2] );
 }

Referenced by moab::Core::create_scd_sequence(), get_adj_edge_or_face(), and main().

◆ get_vertex() [2/2]

EntityHandle moab::ScdBox::get_vertex	(	int	i,
		int	j = `0`,
		int	k = `0`
	)		const

inline

Get the vertex at the specified coordinates.

Parameters

i	Parametric coordinates being evaluated
j	Parametric coordinates being evaluated
k	Parametric coordinates being evaluated

Definition at line 1566 of file ScdInterface.hpp.

 {
     return ( vertDat
                  ? startVertex + ( boxDims[2] == -1 && boxDims[5] == -1 ? 0 : ( k - boxDims[2] ) ) * boxSizeIJ +
                        ( boxDims[1] == -1 && boxDims[4] == -1 ? 0 : ( j - boxDims[1] ) ) * boxSize[0] + i - boxDims[0]
                  : get_vertex_from_seq( i, j, k ) );
 }

References boxDims, boxSize, boxSizeIJ, get_vertex_from_seq(), startVertex, and vertDat.

◆ get_vertex_from_seq()

EntityHandle moab::ScdBox::get_vertex_from_seq	(	int	i,
		int	j,
		int	k
	)		const

private

function to get vertex handle directly from sequence

Parameters

i	Parameter being queried
j	Parameter being queried
k	Parameter being queried

Definition at line 569 of file ScdInterface.cpp.

 {
     assert( elemSeq );
     return elemSeq->get_vertex( i, j, k );
 }

References elemSeq, and moab::StructuredElementSeq::get_vertex().

Referenced by get_vertex().

◆ locally_periodic() [1/2]

const int * moab::ScdBox::locally_periodic ( ) const

inline

Get local periodicity.

Returns: Vector of ijk periodicities for this box

Definition at line 1659 of file ScdInterface.hpp.

 {
     return locallyPeriodic;
 }

References locallyPeriodic.

◆ locally_periodic() [2/2]

void moab::ScdBox::locally_periodic ( bool lperiodic[3] )

inline

Set local periodicity.

Parameters

lperiodic Vector of ijk periodicities to set this box to

Definition at line 1653 of file ScdInterface.hpp.

 {
     for( int i = 0; i < 3; i++ )
         locallyPeriodic[i] = lperiodic[i];
 }

References locallyPeriodic.

Referenced by moab::ScdInterface::assign_global_ids().

◆ locally_periodic_i()

bool moab::ScdBox::locally_periodic_i ( ) const

inline

Return whether box is locally periodic in i.

Return whether box is locally periodic in i

Returns: True if box is locally periodic in i direction

Definition at line 1638 of file ScdInterface.hpp.

 {
     return locallyPeriodic[0];
 }

References locallyPeriodic.

◆ locally_periodic_j()

bool moab::ScdBox::locally_periodic_j ( ) const

inline

Return whether box is locally periodic in j.

Return whether box is locally periodic in j

Returns: True if box is locally periodic in j direction

Definition at line 1643 of file ScdInterface.hpp.

 {
     return locallyPeriodic[1];
 }

References locallyPeriodic.

◆ locally_periodic_k()

bool moab::ScdBox::locally_periodic_k ( ) const

inline

Return whether box is locally periodic in k.

Return whether box is locally periodic in k

Returns: True if box is locally periodic in k direction

Definition at line 1648 of file ScdInterface.hpp.

 {
     return locallyPeriodic[2];
 }

References locallyPeriodic.

◆ num_elements()

int moab::ScdBox::num_elements ( ) const

inline

Return the number of elements in the box.

Definition at line 1493 of file ScdInterface.hpp.

 {
     if( !startElem ) return 0;  // not initialized yet
  
     /* for a structured mesh, total number of elements is obtained by multiplying
         number of elements in each direction
       number of elements in each direction is given by number of vertices in that direction minus 1
       if periodic in that direction, the last vertex is the same as first one, count one more
       element
       */
     int num_e_i = ( -1 == boxSize[0] || 1 == boxSize[0] ) ? 1 : boxSize[0] - 1;
     if( locallyPeriodic[0] ) ++num_e_i;
  
     int num_e_j = ( -1 == boxSize[1] || 1 == boxSize[1] ) ? 1 : boxSize[1] - 1;
     if( locallyPeriodic[1] ) ++num_e_j;
  
     int num_e_k = ( -1 == boxSize[2] || 1 == boxSize[2] ) ? 1 : boxSize[2] - 1;
     if( locallyPeriodic[2] ) ++num_e_k;
  
     return num_e_i * num_e_j * num_e_k;
 }

References boxSize, locallyPeriodic, and startElem.

Referenced by moab::ScdInterface::construct_box(), moab::ScdNCHelper::create_mesh(), moab::ScdInterface::create_scd_sequence(), and moab::ScdInterface::tag_shared_vertices().

◆ num_vertices()

int moab::ScdBox::num_vertices ( ) const

inline

Return the number of vertices in the box.

Definition at line 1515 of file ScdInterface.hpp.

 {
     return boxSize[0] * ( !boxSize[1] ? 1 : boxSize[1] ) * ( !boxSize[2] ? 1 : boxSize[2] );
 }

References boxSize.

Referenced by moab::ScdInterface::assign_global_ids(), moab::ScdNCHelper::create_mesh(), and moab::ScdInterface::create_scd_sequence().

◆ par_data() [1/3]

ScdParData& moab::ScdBox::par_data ( )

inline

Return parallel data.

Return parallel data, if there is any

Returns: par_data Parallel data set on this box

Definition at line 749 of file ScdInterface.hpp.

     {
         return parData;
     }

Referenced by moab::ScdInterface::assign_global_ids(), moab::ScdInterface::construct_box(), moab::ScdInterface::get_shared_vertices(), and moab::ScdInterface::tag_shared_vertices().

◆ par_data() [2/3]

const ScdParData& moab::ScdBox::par_data ( ) const

inline

Return parallel data.

Return parallel data, if there is any

Returns: par_data Parallel data set on this box

Definition at line 758 of file ScdInterface.hpp.

     {
         return parData;
     }

◆ par_data() [3/3]

void moab::ScdBox::par_data ( const ScdParData & par_datap )

inline

set parallel data

Set parallel data for this box

Parameters

par_data Parallel data to be set on this box

Definition at line 767 of file ScdInterface.hpp.

     {
         parData = par_datap;
     }

◆ sc_impl()

ScdInterface * moab::ScdBox::sc_impl ( ) const

inline

Return the ScdInterface responsible for this box.

Definition at line 1468 of file ScdInterface.hpp.

 {
     return scImpl;
 }

References scImpl.

◆ start_element() [1/2]

EntityHandle moab::ScdBox::start_element ( ) const

inline

Starting entity handle for this box.

If this is a vertex box, the start vertex handle is returned.

Definition at line 1483 of file ScdInterface.hpp.

 {
     return startElem;
 }

References startElem.

Referenced by moab::ScdInterface::construct_box(), moab::ScdNCHelper::create_mesh(), and moab::ScdInterface::tag_shared_vertices().

◆ start_element() [2/2]

void moab::ScdBox::start_element ( EntityHandle starte )

inlineprivate

Set the starting entity handle for this box.

Definition at line 1488 of file ScdInterface.hpp.

 {
     startElem = starte;
 }

References startElem.

◆ start_vertex() [1/2]

EntityHandle moab::ScdBox::start_vertex ( ) const

inline

Starting vertex handle for this box.

Definition at line 1473 of file ScdInterface.hpp.

 {
     return startVertex;
 }

References startVertex.

Referenced by moab::ScdInterface::assign_global_ids(), moab::ScdNCHelper::create_mesh(), and moab::ScdInterface::tag_shared_vertices().

◆ start_vertex() [2/2]

void moab::ScdBox::start_vertex ( EntityHandle startv )

inlineprivate

Set the starting vertex handle for this box.

Definition at line 1478 of file ScdInterface.hpp.

 {
     startVertex = startv;
 }

References startVertex.

◆ vert_dat() [1/2]

ScdVertexData * moab::ScdBox::vert_dat ( ) const

inlineprivate

get the vertex sequence

Definition at line 1599 of file ScdInterface.hpp.

 {
     return vertDat;
 }

References vertDat.

◆ vert_dat() [2/2]

ErrorCode moab::ScdBox::vert_dat ( ScdVertexData * vert_dat )

private

set the vertex sequence

Definition at line 622 of file ScdInterface.cpp.

 {
     vertDat = vert_dt;
     return MB_SUCCESS;
 }

References MB_SUCCESS, and vertDat.

Friends And Related Function Documentation

◆ ScdInterface

friend class ScdInterface

friend

Definition at line 508 of file ScdInterface.hpp.

Member Data Documentation

◆ boxDims

int moab::ScdBox::boxDims[6]

private

lower and upper corners

Definition at line 832 of file ScdInterface.hpp.

Referenced by box_dims(), box_max(), box_min(), contains(), get_adj_edge_or_face(), get_element(), get_vertex(), and ScdBox().

◆ boxSet

EntityHandle moab::ScdBox::boxSet

private

entity set representing this box

Definition at line 816 of file ScdInterface.hpp.

Referenced by box_set(), and ~ScdBox().

◆ boxSize

HomCoord moab::ScdBox::boxSize

private

parameter extents

Definition at line 841 of file ScdInterface.hpp.

Referenced by box_size(), elem_seq(), get_vertex(), num_elements(), num_vertices(), and ScdBox().

◆ boxSizeIJ

int moab::ScdBox::boxSizeIJ

private

convenience parameters, (boxSize[1]-1)*(boxSize[0]-1) and boxSize[0]-1

Definition at line 844 of file ScdInterface.hpp.

Referenced by get_vertex(), and ScdBox().

◆ boxSizeIJM1

int moab::ScdBox::boxSizeIJM1

private

Definition at line 845 of file ScdInterface.hpp.

Referenced by elem_seq(), get_element(), and ScdBox().

◆ boxSizeIM1

int moab::ScdBox::boxSizeIM1

private

Definition at line 846 of file ScdInterface.hpp.

Referenced by elem_seq(), get_element(), and ScdBox().

◆ elemSeq

StructuredElementSeq* moab::ScdBox::elemSeq

private

element sequence this box represents

Definition at line 823 of file ScdInterface.hpp.

Referenced by add_vbox(), boundary_complete(), elem_seq(), get_params(), get_vertex_from_seq(), and ScdBox().

◆ locallyPeriodic

int moab::ScdBox::locallyPeriodic[3]

private

is locally periodic in i or j or k

Definition at line 835 of file ScdInterface.hpp.

Referenced by elem_seq(), get_adj_edge_or_face(), locally_periodic(), locally_periodic_i(), locally_periodic_j(), locally_periodic_k(), num_elements(), and ScdBox().

◆ parData

ScdParData moab::ScdBox::parData

private

parallel data associated with this box, if any

Definition at line 838 of file ScdInterface.hpp.

◆ scImpl

ScdInterface* moab::ScdBox::scImpl

private

interface instance

Definition at line 813 of file ScdInterface.hpp.

Referenced by box_dimension(), get_adj_edge_or_face(), get_params(), sc_impl(), ScdBox(), and ~ScdBox().

◆ startElem

EntityHandle moab::ScdBox::startElem

private

starting element handle for this box

Definition at line 829 of file ScdInterface.hpp.

Referenced by box_dimension(), get_element(), num_elements(), ScdBox(), and start_element().

◆ startVertex

EntityHandle moab::ScdBox::startVertex

private

starting vertex handle for this box

Definition at line 826 of file ScdInterface.hpp.

Referenced by get_vertex(), ScdBox(), and start_vertex().

◆ vertDat

ScdVertexData* moab::ScdBox::vertDat

private

vertex sequence this box represents, if there's only one, otherwise they're retrieved from the element sequence

Definition at line 820 of file ScdInterface.hpp.

Referenced by add_vbox(), get_coordinate_arrays(), get_params(), get_vertex(), ScdBox(), and vert_dat().

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

Public Member Functions

Private Member Functions

Private Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ ~ScdBox()

◆ ScdBox()

Member Function Documentation

◆ add_vbox()

◆ boundary_complete()

◆ box_dimension()

◆ box_dims()

◆ box_max()

◆ box_min()

◆ box_set() [1/2]

◆ box_set() [2/2]

◆ box_size() [1/3]

◆ box_size() [2/3]

◆ box_size() [3/3]

◆ contains() [1/2]

◆ contains() [2/2]

◆ elem_seq() [1/2]

◆ elem_seq() [2/2]

◆ get_adj_edge_or_face()

◆ get_coordinate_arrays() [1/2]

◆ get_coordinate_arrays() [2/2]

◆ get_element() [1/2]

◆ get_element() [2/2]

◆ get_params() [1/3]

◆ get_params() [2/3]

◆ get_params() [3/3]

◆ get_vertex() [1/2]

◆ get_vertex() [2/2]

◆ get_vertex_from_seq()

◆ locally_periodic() [1/2]

◆ locally_periodic() [2/2]

◆ locally_periodic_i()

◆ locally_periodic_j()

◆ locally_periodic_k()

◆ num_elements()

◆ num_vertices()

◆ par_data() [1/3]

◆ par_data() [2/3]

◆ par_data() [3/3]

◆ sc_impl()

◆ start_element() [1/2]

◆ start_element() [2/2]

◆ start_vertex() [1/2]

◆ start_vertex() [2/2]

◆ vert_dat() [1/2]

◆ vert_dat() [2/2]

Friends And Related Function Documentation

◆ ScdInterface

Member Data Documentation

◆ boxDims

◆ boxSet

◆ boxSize

◆ boxSizeIJ

◆ boxSizeIJM1

◆ boxSizeIM1

◆ elemSeq

◆ locallyPeriodic

◆ parData

◆ scImpl

◆ startElem

◆ startVertex

◆ vertDat