You are here: Concepts>Concepts Web>Class documentation

Class documentation of Concepts

Loading...
Searching...
No Matches
Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
concepts::Quadrature< type > Class Template Reference

#include <quadrature.hh>

Inheritance diagram for concepts::Quadrature< type >:
concepts::OutputOperator

Public Member Functions

 Quadrature (uint n)
 
const Real * abscissas () const
 Returns a pointer into the array of the abscissas.
 
const Real * weights () const
 Returns a pointer into the array of the weights.
 
uint n () const
 Returns the number of quadrature points.
 

Static Public Member Functions

static std::string storedPoints ()
 

Protected Member Functions

virtual std::ostream & info (std::ostream &os) const
 Returns information in an output stream.
 

Protected Attributes

Real * abscissas_
 Abscissas.
 
Real * weights_
 Weights.
 

Detailed Description

template<int type>
class concepts::Quadrature< type >

Basic class for numerical integration. This class provides the datastructures for the weights and abscissas and the methods to get them.

The values are computed in the template instantiations.

Parameters
typeTemplate parameter: type of the quadrature rule. Can take the values
  • 0: Gauss Lobatto quadrature. This rule includes both endpoints.

    \[ \int_{-1}^1 f(x) \, dx \approx \sum_{i=0}^p w_i f(x_i) \]

    is exact for $f \in P_{2p-1}$ and n = p+1 points. n must be greater or equal to 2.
    The abscissas $x_i$ are the zeros of $(1-x^2) P_{p-1}^{(1,1)}(x)$ and the weights are $w_i = 2/(p(p+1) (P_p^{(0,0)}(x_i))^2)$.
  • 1: Gauss Jacobi Lobatto quadrature. This rule includes both endpoints.

    \[ \int_{-1}^1 f(x) (1-x) \, dx \approx \sum_{i=0}^p w_i f(x_i) \]

    is exact for $f \in P_{2p-1}$ and n = p+1 points. n must be greater or equal to 2.
    The abscissas $x_i$ are the zeros of $(1-x^2) P_{p-1}^{(2,1)}(x)$ and the weights are $w_i = 4/(p(p+2) (P_p^{(1,0)}(x_i))^2)$ and $w_p = 8/(p(p+2) (P_p^{(1,0)}(x_i))^2)$.
  • 2: Gauss Jacobi Lobatto quadrature. This rule includes both endpoints.

    \[ \int_{-1}^1 f(x) (1-x)^2 \, dx \approx \sum_{i=0}^p w_i f(x_i)\]

    is exact for $f \in P_{2p-1}$ and n = p+1 points. n must be greater or equal to 2.
    The abscissas $x_i$ are the zeros of $(1-x^2) P_{p-1}^{(3,1)}(x)$ and the weights are $w_i = 8/(p(p+3) (P_p^{(2,0)}(x_i))^2)$ and $w_p = 24/(p(p+3) (P_p^{(2,0)}(x_i))^2)$.
  • 3: Gauss Radau Jacobi quadrature. This rule includes only one of the endpoints, ie. -1.

    \[ \int_{-1}^1 f(x) (1-x) \, dx \approx \sum_{i=0}^p w_i f(x_i) \]

    is exact for $f \in P_{2p}$ and n = p+1 points. n must be greater or equal to 2.
    The abscissas $x_i$ are the zeros of $(1+x) P_p^{(1,1)}(x)$ and the weights are $w_i = 2(1-x_i)/(p(p+1) (P_p^{(1,0)}(x_1))^2)$.
  • 4: Gauss Jacobi quadrature. This rule does not include the endpoints -1 and 1.

    \[ \int_{-1}^1 f(x) \, dx \approx \sum_{i=0}^p w_i f(x_i) \]

    is exact for $f \in P_{2p+1}$ and n = p+1 points. n must be greater or equal to 1.
    The abscissas $x_i$ are the zeros of $P_{p+1}^{(0,0)}(x)$ and the weights are

    \[w_i = \frac{2}{1-x_i^2} \left( \frac{d}{dx} \left. P^{(0,0)}_{p+1}(x) \right|_{x=x_i} \right)^{-2}.\]

  • 5: Trapeze quadrature. This rule include the endpoints -1 and 1.

    \[ \int_{-1}^1 f(x) \, dx \approx \sum_{i=0}^{n-1} w_i f(x_i) \]

    is only exact for $f \in P_{1}$. n must be greater or equal to 2.
    The abscissas $x_i = -1 + 2i/(n-1)$ are equidistant and the weights are

    \[w_i = \frac{1}{n-1}\left\{\begin{array}{ll}1 & i \in \{0,n\}\\2 & \mbox{otherwise}\end{array}\right.\]

Requesting the same quadrature rule with the same amount of integration points several times does not harm: the values are only stored once (internally).

Test:

test::QuadratureTest

test::KarniadakisTest

Author
Philipp Frauenfelder, 2000
Examples
linearFEM1d-simple.cc, and linearFEM1d.cc.

Definition at line 97 of file quadrature.hh.

Constructor & Destructor Documentation

◆ Quadrature()

template<int type>
concepts::Quadrature< type >::Quadrature ( uint  n)

Constructor. Computes the quadrature points.

Parameters
nNumber of entries to be computed

Member Function Documentation

◆ abscissas()

template<int type>
const Real * concepts::Quadrature< type >::abscissas ( ) const
inline

Returns a pointer into the array of the abscissas.

Examples
linearFEM1d-simple.cc, and linearFEM1d.cc.

Definition at line 105 of file quadrature.hh.

◆ info()

template<int type>
virtual std::ostream & concepts::Quadrature< type >::info ( std::ostream &  os) const
protectedvirtual

Returns information in an output stream.

Reimplemented from concepts::OutputOperator.

◆ n()

template<int type>
uint concepts::Quadrature< type >::n ( ) const
inline

Returns the number of quadrature points.

Examples
linearFEM1d-simple.cc, and linearFEM1d.cc.

Definition at line 109 of file quadrature.hh.

◆ storedPoints()

template<int type>
static std::string concepts::Quadrature< type >::storedPoints ( )
static

Lists the integration orders (points) which have been computed so far.

Returns
a string containing the computed orders.

◆ weights()

template<int type>
const Real * concepts::Quadrature< type >::weights ( ) const
inline

Returns a pointer into the array of the weights.

Examples
linearFEM1d-simple.cc, and linearFEM1d.cc.

Definition at line 107 of file quadrature.hh.

Member Data Documentation

◆ abscissas_

template<int type>
Real* concepts::Quadrature< type >::abscissas_
protected

Abscissas.

Definition at line 119 of file quadrature.hh.

◆ weights_

template<int type>
Real* concepts::Quadrature< type >::weights_
protected

Weights.

Definition at line 121 of file quadrature.hh.

The documentation for this class was generated from the following file:
Generated on Wed Sep 13 2023 21:06:43 for Concepts by doxygen 1.9.8