oomph::GaussLobattoLegendre< 2, NPTS_1D > Class Template Reference

2D Gauss Lobatto Legendre integration class More...

#include <integral.h>

Inheritance diagram for oomph::GaussLobattoLegendre< 2, NPTS_1D >:

Public Member Functions
	GaussLobattoLegendre ()
	Deafault constructor. Calculates and stores GLL nodes. More...
unsigned	nweight () const
	Number of integration points of the scheme. More...
double	knot (const unsigned &i, const unsigned &j) const
	Return coordinate s[j] (j=0) of integration point i. More...
double	weight (const unsigned &i) const
	Return weight of integration point i. More...
Public Member Functions inherited from oomph::Integral
	Integral ()
	Default constructor (empty) More...
	Integral (const Integral &dummy)=delete
	Broken copy constructor. More...
void	operator= (const Integral &)=delete
	Broken assignment operator. More...
virtual	~Integral ()
	Virtual destructor (empty) More...
virtual Vector< double >	knot (const unsigned &i) const
	Return local coordinates of i-th intergration point. Broken virtual. More...

Private Attributes
double	Knot [NPTS_1D *NPTS_1D][2]
	Array to hold weight and knot points. More...
double	Weight [NPTS_1D *NPTS_1D]

Static Private Attributes
static const unsigned long int	Npts = NPTS_1D * NPTS_1D
	Number of integration points in scheme. More...

Detailed Description

template<unsigned NPTS_1D>
class oomph::GaussLobattoLegendre< 2, NPTS_1D >

2D Gauss Lobatto Legendre integration class

Constructor & Destructor Documentation

GaussLobattoLegendre()

template<unsigned NPTS_1D>

oomph::GaussLobattoLegendre< 2, NPTS_1D >::GaussLobattoLegendre

Deafault constructor. Calculates and stores GLL nodes.

Calculate positions and weights for the 2D Gauss Lobatto Legendre integration class

  {
    // Tempoarary storage for the 1D knots and weights
    Vector<double> s(NPTS_1D), w(NPTS_1D);
    // Call the function to populate the array
    Orthpoly::gll_nodes(NPTS_1D, s, w);
    int i_fast = 0, i_slow = 0;
    for (unsigned i = 0; i < NPTS_1D * NPTS_1D; i++)
    {
      if (i_fast == NPTS_1D)
      {
        i_fast = 0;
        i_slow++;
      }
      Knot[i][0] = s[i_fast];
      Knot[i][1] = s[i_slow];
      Weight[i] = w[i_fast] * w[i_slow];
      i_fast++;
    }
  }

References oomph::Orthpoly::gll_nodes(), i, s, w, and ProblemParameters::Weight.

Member Function Documentation

knot()

template<unsigned NPTS_1D>

double oomph::GaussLobattoLegendre< 2, NPTS_1D >::knot ( const unsigned &  i, const unsigned &  j  ) const

inlinevirtual

Return coordinate s[j] (j=0) of integration point i.

Implements oomph::Integral.

    {
      return Knot[i][j];
    }

References i, and j.

nweight()

template<unsigned NPTS_1D>

unsigned oomph::GaussLobattoLegendre< 2, NPTS_1D >::nweight ( ) const

inlinevirtual

Number of integration points of the scheme.

Implements oomph::Integral.

    {
      return Npts;
    }

weight()

template<unsigned NPTS_1D>

double oomph::GaussLobattoLegendre< 2, NPTS_1D >::weight ( const unsigned &  i ) const

inlinevirtual

Return weight of integration point i.

Implements oomph::Integral.

    {
      return Weight[i];
    }

References i, and ProblemParameters::Weight.

Member Data Documentation

Knot

template<unsigned NPTS_1D>

double oomph::GaussLobattoLegendre< 2, NPTS_1D >::Knot[NPTS_1D *NPTS_1D][2]

private

Array to hold weight and knot points.

Npts

template<unsigned NPTS_1D>

const unsigned long int oomph::GaussLobattoLegendre< 2, NPTS_1D >::Npts = NPTS_1D * NPTS_1D

staticprivate

Number of integration points in scheme.

Weight

template<unsigned NPTS_1D>

double oomph::GaussLobattoLegendre< 2, NPTS_1D >::Weight[NPTS_1D *NPTS_1D]

private

---

The documentation for this class was generated from the following file:

• integral.h