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 HArDCore2D;
 
// ----------------------------------------------------------------------------
//                            Class definition
// ----------------------------------------------------------------------------
 
class BoundaryConditions {
 
public:
  BoundaryConditions(
                     const std::string bc_id,  
                     Mesh& mesh          
                     );
 
  const std::string type(
                         const Edge& edge    
                         ) const ;
  const std::string type(
                         const Vertex& vertex    
                         ) const ;
 
  inline const size_t n_dir_edges() const {
    return m_n_dir_edges;
  };
 
  inline const size_t n_dir_vertices() const {
    return m_n_dir_vertices;
  };
 
  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_edges(const std::string pos = "end");
 
  void reorder_vertices(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 edges and vertices
  size_t m_n_dir_edges;
  size_t m_n_dir_vertices;
 
};
 
 
 
#endif //_BOUNDARY_CONDITION_HPP