BEBilinearFormLaplace1Layer.h 1.95 KB
Newer Older
Jan Zapletal's avatar
Jan Zapletal committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90
/*!
 * @file    BEBilinearFormLaplace1Layer.h
 * @author  Michal Merta 
 * @date    July 12, 2013
 * @brief   Header file for class BEBilinearFormLaplace1Layer
 * 
 */

#ifndef BEBILINEARFORMLAPLACE1LAYER_H
#define	BEBILINEARFORMLAPLACE1LAYER_H

#include "BEBilinearForm.h"

namespace bem4i {

/*! 
 * Class representing the bilinear form for the Laplace single layer operator
 * 
 * Provides methods for system matrix assembly.
 * 
 */
template<class LO, class SC>
class BEBilinearFormLaplace1Layer : public BEBilinearForm<LO, SC> {
  typedef typename GetType<LO, SC>::SCVT SCVT;



public:
  using BEBilinearForm<LO, SC>::assemble;

  //! default constructor
  BEBilinearFormLaplace1Layer( );

  //! copy constructor
  BEBilinearFormLaplace1Layer(
      const BEBilinearFormLaplace1Layer& orig
      );

  /*! 
   * Constructor taking the boundary element space on which the form acts
   */
  BEBilinearFormLaplace1Layer(
      BESpace<LO, SC>* space,
      int* quadratureOrder = nullptr,
      quadratureType quadrature = SauterSchwab,
      int* quadratureOrderDisjointElems = nullptr,
      bool employSymmetricity = false
      );

  //! destructor
  virtual ~BEBilinearFormLaplace1Layer( );

  /*!
   * The method assembles the Galerkin matrix for the Laplace single layer operator
   */
  virtual void assemble(
      FullMatrix<LO, SC> & matrix
      ) const;

  /*!
   * The method assembles the Galerkin matrix for the Laplace single layer operator
   */
  void assembleP0P0(
      FullMatrix<LO, SC> & matrix
      ) const;

  void assembleP1P1(
      FullMatrix<LO, SC> & matrix
      ) const;
  
  void assembleP1DisP1Dis(
      FullMatrix<LO, SC> & matrix
      ) const;

protected:

private:

  //! default quadrature rule for this class (specific values set in .cpp file)
  static int defaultQuadratureOrder[2];

};

}

// include .cpp file to overcome linking problems due to templates
#include "BEBilinearFormLaplace1Layer.cpp"

#endif	/* BEBILINEARFORMLAPLACE1LAYER_H */