NodeCoupling.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/NodeCoupling.hpp"
 
NodeCoupling::
NodeCoupling()
{
}
 
NodeCoupling::
~NodeCoupling()
{
}
 
const std::array<CepsString,2>& NodeCoupling::
getMeshes() const
{
  return m_meshes;
}
 
const CepsString& NodeCoupling::
getMesh(CepsUInt id) const
{
  CEPS_ABORT_IF(id>1,"A node coupling has only two meshes, but getMesh(i) was called with i="<< id);
  return m_meshes[id];
}
 
const CepsSet<std::pair<CepsNodeGlobalIndex,CepsNodeGlobalIndex>>& NodeCoupling::
getNodeCouples() const
{
  return m_couples;
}
 
CepsBool NodeCoupling::
isValidIn(const CepsVector<CepsString>& meshes)
{
 
  if (m_meshes[1].empty())
  {
    CEPS_WARNS("Node coupling must have exactly two meshes registered to be valid");
    return false;
  }
 
  return ceps::contains(meshes,m_meshes[0]) and ceps::contains(meshes,m_meshes[1]);
}
 
void NodeCoupling::
addMesh(const CepsString& mesh)
{
  if (m_meshes[0].empty())
    m_meshes[0] = mesh;
  else if (m_meshes[1].empty())
    m_meshes[1] = mesh;
  else
    CEPS_ABORT("Cannot register more than two meshes in node coupling");
  return;
}
 
void NodeCoupling::
addCouple(const CepsNodeGlobalIndex& id1, const CepsNodeGlobalIndex& id2)
{
  std::pair<CepsNodeGlobalIndex,CepsNodeGlobalIndex> p(id1,id2);
  m_couples.insert(p);
}