packages/MC/geometry/Mesh_Geometry.cc
//---------------------------------------------------------------------------//
/*!
 * \brief Get bounding box for cell
 */
\brief Get bounding box for cell
\ brief Get bounding box for cell
Bounding_Box Mesh_Geometry::get_cell_extents(geometry::cell_type cell) const
{
    REQUIRE( cell < d_mesh.num_cells() );
    using def::I;  using def::J;  using def::K;
    auto ijk = d_mesh.cardinal(cell);
    const auto &edges = d_mesh.edges();
    return Bounding_Box( edges[I][ijk[I]],
                         edges[I][ijk[I]+1],
                         edges[J][ijk[J]],
                         edges[J][ijk[J]+1],
                         edges[K][ijk[K]],
                         edges[K][ijk[K]+1]);
}
Bounding_Box Mesh_Geometry :: get_cell_extents ( geometry :: cell_type cell ) const { REQUIRE ( cell < d_mesh . num_cells ( ) ) ; using def :: I ; using def :: J ; using def :: K ; auto ijk = d_mesh . cardinal ( cell ) ; const auto & edges = d_mesh . edges ( ) ; return Bounding_Box ( edges [ I ] [ ijk [ I ] ] , edges [ I ] [ ijk [ I ] + 1 ] , edges [ J ] [ ijk [ J ] ] , edges [ J ] [ ijk [ J ] + 1 ] , edges [ K ] [ ijk [ K ] ] , edges [ K ] [ ijk [ K ] + 1 ] ) ; }
\brief Get bounding box for cell
\ brief Get bounding box for cell