CoupledNodesReader.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 "CoupledNodesReader.hpp"
 
CoupledNodesReader::CoupledNodesReader (const CepsString &fileName) :
    FileReader (fileName)
{
  CEPS_ABORT_IF(not open(),"Unable to open file " << m_fileName);
}
 
CoupledNodesReader::~CoupledNodesReader ()
{}
 
CepsVector<NodeCoupling>
CoupledNodesReader::getCoupledNodesBetween (const CepsVector<CepsString> &meshes)
{
  CepsVector<NodeCoupling> res;
  CEPS_ABORT_IF (meshes.size()<2,"Coupled nodes can be registered only between two meshes");
 
  this->reset ();
 
  CepsUInt            n;
  CepsString          mesh;
  CepsNodeGlobalIndex i1, i2;
 
  while (findNext ("COUPLING"))
  {
    NodeCoupling c;
    for (CepsUInt i = 0u; i < 2; i++)
    {
      CEPS_ABORT_IF(not (m_file >> mesh),
        "Impossible to read mesh names after CONNECTION keyword in file " << m_fileName
      );
      CEPS_ABORT_IF (not ceps::fileExists (mesh),
        "Cannot find mesh given after CONNECTION keyword" << mesh
      );
      c.addMesh (mesh);
    }
 
    // Is the connection valid ?
    if (c.isValidIn (meshes))
    {
      n = ceps::readInt(m_file,"Impossible to read number of couples after COUPLING keyword in file " + m_fileName);
      for (CepsUInt i = 0u; i < n; i++)
      {
        i1 = ceps::readInt(m_file,"Impossible to read node IDs after COUPLING keyword in file " + m_fileName);
        i2 = ceps::readInt(m_file,"Impossible to read node IDs after COUPLING keyword in file " + m_fileName);
        c.addCouple (i1, i2);
      }
      // Couples read, add them
      res.push_back (c);
    }
  }
 
  return res;
}