FEType.hpp

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
 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
#pragma once
 
#include "geometry/base/GeomCell.hpp"
 
using CepsFeLagrangeType = std::pair<CepsCellType,CepsUInt>;
 
inline CepsString
toString(CepsCellType cat, CepsUInt order)
{
  if (cat == CepsCellType::Quad)
    return "Q" + ceps::toString(order);
  return "P" + ceps::toString(order);
}
 
inline CepsFeLagrangeType
fromString(CepsString type)
{
  CepsChar k = ceps::toKey (type).at (0);
  CepsUInt i = ceps::toUInt(type.substr (1));
 
  CepsBool ok = (k=='P' or k=='Q');
  CEPS_ABORT_IF(not ok,
    "Unknown finite element type. Valid types are P[1,2,...] and Q[1,2,...]."
  );
 
  if (k == 'P')
    return {CepsCellType::Simplex,i};
  return {CepsCellType::Quad,i};
 
}
 
inline CepsCellType
getCategory (const CepsFeLagrangeType &type)
{
  return type.first;
}
 
inline CepsUInt
getOrder (const CepsFeLagrangeType &type)
{
  return type.second;
}