Panzer_CellData.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//           Panzer: A partial differential equation assembly
//       engine for strongly coupled complex multiphysics systems
//
// Copyright 2011 NTESS and the Panzer contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER

#ifndef PANZER_CELL_DATA_HPP
#define PANZER_CELL_DATA_HPP

#include "PanzerDiscFE_config.hpp"

#include "Teuchos_Assert.hpp"
#include "Teuchos_RCP.hpp"

#include "Shards_CellTopology.hpp"
#include "Shards_BasicTopologies.hpp"

namespace panzer {
  
  class CellData {
    
  public:

    explicit CellData() 
    { }
    
    explicit CellData(std::size_t num_cells,
                      const Teuchos::RCP<const shards::CellTopology> & ct) :
      m_num_cells(num_cells),
      m_is_side(false),
      m_side(-1),
      m_cell_topo(ct)
    { }

    explicit CellData(std::size_t num_cells,
          int local_side_id,const Teuchos::RCP<const shards::CellTopology> & ct) :
      m_num_cells(num_cells),
      m_is_side(local_side_id >= 0),
      m_side(local_side_id),
      m_cell_topo(ct)
    { }
  
    bool isSide() const 
    {return m_is_side;}
    
    int side() const 
    {
      TEUCHOS_TEST_FOR_EXCEPTION(!m_is_side, std::logic_error,
       "Cannot return side index, CellData is not a side!");
      return m_side;
    }
      
    std::size_t numCells() const 
    {return m_num_cells;}
      
    int baseCellDimension() const 
    {return m_cell_topo->getDimension();}
    
    // TODO BWR this terminology is a bit confusing... this is not the base cell topology
    Teuchos::RCP<const shards::CellTopology> getCellTopology() const
    { return m_cell_topo; }

  private:
    std::size_t m_num_cells;
    bool m_is_side;
    int m_side;
      
 
    Teuchos::RCP<const shards::CellTopology> m_cell_topo;
  };

}

#endif