Subproblem3D.cpp

#include "Subproblem3D.h"

using gmshfem::equation::dof;
using gmshfem::equation::tf;
using gmshfem::function::operator-;
using gmshfem::function::operator+;
using gmshfem::function::operator*;

namespace D3 {


  // **********************************
  // Boundary
  // **********************************

  Boundary::Boundary(const unsigned int boundary) : _boundary(boundary)
  {
  }

  Boundary::~Boundary()
  {
  }

  std::string Boundary::orientation() const
  {
    switch (_boundary) {
      case 0:
        return "E";
        break;
      case 1:
        return "N";
        break;
      case 2:
        return "W";
        break;
      case 3:
        return "S";
        break;
      case 4:
        return "D";
        break;
      case 5:
        return "U";
        break;
      default:
        break;
    }
    return "null";
  }
  
  // **********************************
  // Edge
  // **********************************

  Edge::Edge(const unsigned int edge, Boundary *first, Boundary *second) : _edge(edge), _bnd {first, second}
  {
  }

  Edge::~Edge()
  {
  }

  std::string Edge::orientation() const
  {
    switch (_edge) {
      case 0:
        return "ED";
        break;
      case 1:
        return "ND";