BEBilinearFormLaplace2Layer.h 2.04 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 91 92 93
/*!
 * @file    BEBilinearFormLaplace2Layer.h
 * @author  Michal Merta
 * @date    July 18, 2013
 * @brief   Header file for class BEBilinearFormLaplace2Layer
 *
 */

#ifndef BEBILINEARFORMLAPLACE2LAYER_H
#define BEBILINEARFORMLAPLACE2LAYER_H

#include "BEBilinearForm.h"
#include "Quadratures.h"

namespace bem4i {

/*!
 * Class representing the bilinear form for the Laplace double layer operator
 *
 * Provides methods for system matrix assembly.
 *
 */
template<class LO, class SC>
class BEBilinearFormLaplace2Layer : public BEBilinearForm<LO, SC> {
  // to get inner type of complex numbers (for Helmholtz)
  typedef typename GetType<LO, SC>::SCVT SCVT;

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

  //! default constructor
  BEBilinearFormLaplace2Layer( );

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

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

  //! destructor
  virtual ~BEBilinearFormLaplace2Layer( );

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

  /*!
   * The method assembles the Galerkin matrix for the Laplace double layer
   * operator with p0p1 elements
   */
  void assembleP0P1(
    FullMatrix<LO, SC> & matrix
    ) const;

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

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

  void assembleP0P0(
    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 "BEBilinearFormLaplace2Layer.cpp"

#endif /* BEBILINEARFORMLAPLACE2LAYER_H */