BoundaryConditions.hpp

Go to the documentation of this file.

// Class to provide description of boundary conditions
//
//
// Author: Jerome Droniou (jerome.droniou@monash.edu)
//
 
/*
*
*   This library was developed around HHO methods, although some parts of it have a more
* general purpose. If you use this code or part of it in a scientific publication, 
* please mention the following book as a reference for the underlying principles
* of HHO schemes:
*
* The Hybrid High-Order Method for Polytopal Meshes: Design, Analysis, and Applications.
* D. A. Di Pietro and J. Droniou. 2019, 516p. 
* url: https://hal.archives-ouvertes.fr/hal-02151813.
*
*/
 
 
#ifndef _BOUNDARY_CONDITIONS_HPP
#define _BOUNDARY_CONDITIONS_HPP
 
#include <iostream>
#include <string>
#include <mesh.hpp>
 
using namespace HArDCore3D;
 
// ----------------------------------------------------------------------------
//                            Class definition
// ----------------------------------------------------------------------------
 
// @addtogroup TestCases
 
class BoundaryConditions {
 
public:
  BoundaryConditions(
    const std::string bc_id,  
    Mesh& mesh          
  );
 
    const std::string type(
        Face& face    
    ) const ;
  inline const size_t n_dir_faces() const {
    return m_n_dir_faces;
  };
 
  inline const std::string name() const {
    if (m_bc_id == "D"){
      return "Dirichlet";
    }else if (m_bc_id == "N"){
      return "Neumann";
    }else if (m_bc_id == "M0"){
      return "Mixed #0";
    }
    std::cout << "Unknown boundary conditions: " << m_bc_id << "\n";
    exit(1);
  };
 
  void reorder_faces(const std::string pos = "end");
 
private:
  // Parameters: id of boundary condition, reference to mesh
  const std::string m_bc_id;
  Mesh& m_mesh;
 
  // Number of Dirichlet faces
  size_t m_n_dir_faces;
};
 
 
 
#endif //_BOUNDARY_CONDITION_HPP