packages/MC/geometry/RTK_Cell.cc
//---------------------------------------------------------------------------//
// PUBLIC INTERFACE
//---------------------------------------------------------------------------//
/*!
 * \brief Initialize a geometric state in the pin.
 *
 * \param r a point somewhere in this pin-cell
 * \param state geometric state that gets set
 */
PUBLIC INTERFACE
\brief Initialize a geometric state in the pin.

\param r a point somewhere in this pin-cell
\param state geometric state that gets set
PUBLIC INTERFACE \ brief Initialize a geometric state in the pin . \ param r a point somewhere in this pin - cell \ param state geometric state that gets set
void RTK_Cell::initialize(const Space_Vector &r,
                          Geo_State_t        &state) const
{
    REQUIRE(r[0] >= d_extent[0][LO]);
    REQUIRE(r[0] <= d_extent[0][HI]);
    REQUIRE(r[1] >= d_extent[1][LO]);
    REQUIRE(r[1] <= d_extent[1][HI]);
    REQUIRE(r[2] >= 0.0);
    REQUIRE(r[2] <= d_z);
    state.region = region(r[0], r[1]);
    state.segment = segment(r[0], r[1]);
    state.face = Geo_State_t::NONE;
    state.exiting_face = Geo_State_t::NONE;
}
void RTK_Cell :: initialize ( const Space_Vector & r , Geo_State_t & state ) const { REQUIRE ( r [ 0 ] >= d_extent [ 0 ] [ LO ] ) ; REQUIRE ( r [ 0 ] <= d_extent [ 0 ] [ HI ] ) ; REQUIRE ( r [ 1 ] >= d_extent [ 1 ] [ LO ] ) ; REQUIRE ( r [ 1 ] <= d_extent [ 1 ] [ HI ] ) ; REQUIRE ( r [ 2 ] >= 0.0 ) ; REQUIRE ( r [ 2 ] <= d_z ) ; state . region = region ( r [ 0 ] , r [ 1 ] ) ; state . segment = segment ( r [ 0 ] , r [ 1 ] ) ; state . face = Geo_State_t :: NONE ; state . exiting_face = Geo_State_t :: NONE ; }
PUBLIC INTERFACE
 \brief Initialize a geometric state in the pin.
PUBLIC INTERFACE \ brief Initialize a geometric state in the pin .
// get the current region // get the current segment // initialize to no-face (this should probably be fixed to allow particle // born on a face) // we are not on an exiting face