GeomCellJunction.cpp

Go to the documentation of this file.

/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   This file is part of CEPS.
 
   CEPS is free software: you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.
 
   CEPS is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.
 
   You should have received a copy of the GNU General Public License
   along with CEPS (see file LICENSE at root of project).
   If not, see <https://www.gnu.org/licenses/>.
 
 
   Copyright 2019-2024 Inria, Universite de Bordeaux
 
   Authors, in alphabetical order:
 
     Pierre-Elliott BECUE, Florian CARO, Yves COUDIERE(*), Andjela DAVIDOVIC, 
     Charlie DOUANLA-LONTSI, Marc FUENTES, Mehdi JUHOOR, Michael LEGUEBE(*), 
     Pauline MIGERDITICHAN, Valentin PANNETIER(*), Nejib ZEMZEMI.
     * : currently active authors
 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
#include "geometry/base/GeomCellJunction.hpp"
 
GeomCellJunction::
GeomCellJunction () :
    m_type (JunctionType::TYPE_UNKNOWN)
{
}
 
GeomCellJunction::
~GeomCellJunction () 
{
}
 
const CepsVector<GeomCell *> & GeomCellJunction::
get3dcells () const
{
  return m_3d;
}
 
const CepsVector<GeomCell *> & GeomCellJunction::
get2dcells () const
{
  return m_2d;
}
 
const CepsVector<GeomCell *> & GeomCellJunction::
get1dcells () const
{
  return m_1d;
}
 
JunctionType GeomCellJunction::
getType () const
{
  return m_type;
}
 
void GeomCellJunction::
addCell(GeomCell* cell)
{
  CEPS_ABORT_IF(ceps::isNullPtr(cell),"Cannot add nullptr cell to cell junction");
 
  CepsUInt dim = cell->getDimension();
  switch (dim)
  {
    case 3:
      m_3d.push_back(cell);
      break;
    case 2:
      m_2d.push_back(cell);
      break;
    case 1:
      m_1d.push_back(cell);
      break;
    default:
      CEPS_ABORT("Cannot add cell to dimension for dimension " << dim);
  }
  checkIntegrity();
  return;
}
 
void GeomCellJunction::
addBoundaryCell(GeomCell* cell)
{
  CEPS_ABORT_IF(ceps::isNullPtr(cell),"Cannot add nullptr boundary cell to cell junction");
 
  CepsUInt dim = cell->getDimension();
  switch (dim)
  {
    case 2:
      m_3dBoundary.push_back(cell);
      break;
    case 1:
      m_2dBoundary.push_back(cell);
      break;
    default:
      CEPS_ABORT("Cannot add boundary cell to dimension for dimension " << dim);
  }
  checkIntegrity();
  return;
}
 
void GeomCellJunction::
checkIntegrity ()
{
  // Determine junction cell type
  CepsBool has1d = (m_1d.size () != 0);
  CepsBool has2d = (m_2d.size () != 0) or (m_2dBoundary.size () != 0);
  CepsBool has3d = (m_3d.size () != 0) or (m_3dBoundary.size () != 0);
 
  if (has1d)
  {
    if (has3d)
      m_type = has2d ? JunctionType::TYPE_1D2D3D : JunctionType::TYPE_1D3D;
    else
    {
      if (has2d)
        m_type = JunctionType::TYPE_1D2D;
      else
      {
        m_type = JunctionType::TYPE_UNKNOWN;
        CEPS_WARNS ("unknown junction type for this junction cell.");
      }
    }
  }
  else
  {
    if (has3d and has2d)
      m_type = JunctionType::TYPE_2D3D;
    else
    {
      m_type = JunctionType::TYPE_UNKNOWN;
      CEPS_WARNS ("unknown junction type for this junction cell.");
    }
  }
  return;
}