oomph::GaussLegendre< 3, NPTS_1D > Class Template Reference

3D Gauss Legendre integration class More...

#include <integral.h>

Inheritance diagram for oomph::GaussLegendre< 3, NPTS_1D >:

Public Member Functions
	GaussLegendre ()
	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 *NPTS_1D][3]
	Array to hold weight and knot points. More...
double	Weight [NPTS_1D *NPTS_1D *NPTS_1D]

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

Detailed Description

template<unsigned NPTS_1D>
class oomph::GaussLegendre< 3, NPTS_1D >

3D Gauss Legendre integration class

Constructor & Destructor Documentation

GaussLegendre()

template<unsigned NPTS_1D>

oomph::GaussLegendre< 3, NPTS_1D >::GaussLegendre

Deafault constructor. Calculates and stores GLL nodes.

Calculate positions and weights for the 3D Gauss 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::gl_nodes(NPTS_1D, s, w);
    for (unsigned k = 0; k < NPTS_1D; k++)
    {
      for (unsigned j = 0; j < NPTS_1D; j++)
      {
        for (unsigned i = 0; i < NPTS_1D; i++)
        {
          unsigned index = NPTS_1D * NPTS_1D * k + NPTS_1D * j + i;
          Knot[index][0] = s[i];
          Knot[index][1] = s[j];
          Knot[index][2] = s[k];
          Weight[index] = w[i] * w[j] * w[k];
        }
      }
    }
  }

References oomph::Orthpoly::gl_nodes(), i, j, k, s, w, and ProblemParameters::Weight.

Member Function Documentation

knot()

template<unsigned NPTS_1D>

double oomph::GaussLegendre< 3, 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::GaussLegendre< 3, NPTS_1D >::nweight ( ) const

inlinevirtual

Number of integration points of the scheme.

Implements oomph::Integral.

    {
      return Npts;
    }

weight()

template<unsigned NPTS_1D>

double oomph::GaussLegendre< 3, 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::GaussLegendre< 3, NPTS_1D >::Knot[NPTS_1D *NPTS_1D *NPTS_1D][3]

private

Array to hold weight and knot points.

Npts

template<unsigned NPTS_1D>

const unsigned long int oomph::GaussLegendre< 3, NPTS_1D >::Npts = NPTS_1D * NPTS_1D * NPTS_1D

staticprivate

Number of integration points in scheme.

Weight

template<unsigned NPTS_1D>

double oomph::GaussLegendre< 3, NPTS_1D >::Weight[NPTS_1D *NPTS_1D *NPTS_1D]

private

---

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

• integral.h