ElementIntegrator.hpp

// This file is part of the imaging2 class library.
//
// University of Innsbruck, Infmath Imaging, 2009.
// http://infmath.uibk.ac.at
//
// All rights reserved.


#ifndef FEM_ELEMENTINTEGRATOR_H
#define FEM_ELEMENTINTEGRATOR_H

#include <core/imaging2.hpp>
#include <core/vector_utilities.hpp>
#include <core/matrix_utilities.hpp>


namespace imaging
{
  template<std::size_t N, std::size_t N_NODES>
  class ElementIntegrator
  {
  protected:
    ublas::fixed_vector< ublas::fixed_vector<float_t, N>, N_NODES> _nodes;
    
    ublas::fixed_vector<float_t, N_NODES> _weights;
  public:
    static const std::size_t n_nodes = N_NODES;

    const ublas::fixed_vector<float_t, N> & node(std::size_t i) const { return _nodes(i); }
    
    float_t weight(std::size_t i) const { return _weights(i); }
  };
  
  


  template<std::size_t N_NODES>
  class SquareIntegrator : public ElementIntegrator<2, N_NODES>
  {
  public:
    SquareIntegrator();
  };

  template<>
  SquareIntegrator<4>::SquareIntegrator();
  

  template<std::size_t N_NODES>
  class CubeIntegrator : public ElementIntegrator<3, N_NODES>
  {
  public:
    CubeIntegrator();
  };

  template<>
  CubeIntegrator<8>::CubeIntegrator();
  
  template<std::size_t N_NODES>
  class TriangleIntegrator : public ElementIntegrator<2, N_NODES>
  {
  public:
    TriangleIntegrator();
  };

  template<>
  TriangleIntegrator<1>::TriangleIntegrator();

  template<>
  TriangleIntegrator<4>::TriangleIntegrator();
  
  template<>
  TriangleIntegrator<7>::TriangleIntegrator();


  template<std::size_t N_NODES>
  class TetrahedraIntegrator : public ElementIntegrator<3, N_NODES>
  {
  public:
    TetrahedraIntegrator();
  };

  template<>
  TetrahedraIntegrator<1>::TetrahedraIntegrator();

  
  template<std::size_t N_NODES>
  class IntervalIntegrator : public ElementIntegrator<1, N_NODES>
  {
  public:
    IntervalIntegrator();
  };

  template<>
  IntervalIntegrator<2>::IntervalIntegrator();

  template<>
  IntervalIntegrator<1>::IntervalIntegrator();
  
  class PointIntegrator : public ElementIntegrator<0, 1>
  {
  public:
    PointIntegrator() : ElementIntegrator<0, 1>() { _weights.assign(1.0); }
  };

}


#endif

---