PolynomialQuadrature::Polynomial Class Reference

General polynomial class with function for evaluating a polynomial of a given order at a given point. More...

#include <PolynomialQuadrature.h>

Inheritance diagram for PolynomialQuadrature::Polynomial:

Public Member Functions
	Polynomial ()
virtual	~Polynomial ()=default
virtual void	store (std::ostream &stream, void *context) const
virtual void	store (nlohmann::json &json) const
virtual Real	compute (const unsigned int order, const Real x, const bool normalize=true) const
virtual Real	computeDerivative (const unsigned int order, const Real x, const unsigned int m=1) const
	Computes the mth derivative of polynomial: d^mP_n/dx^m.
virtual Real	innerProduct (const unsigned int order) const =0
virtual void	gaussQuadrature (const unsigned int order, std::vector< Real > &points, std::vector< Real > &weights) const =0
virtual void	clenshawQuadrature (const unsigned int order, std::vector< Real > &points, std::vector< Real > &weights) const
Real	productIntegral (const std::vector< unsigned int > order) const

Detailed Description

General polynomial class with function for evaluating a polynomial of a given order at a given point.

Definition at line 33 of file PolynomialQuadrature.h.

Constructor & Destructor Documentation

Polynomial()

PolynomialQuadrature::Polynomial::Polynomial ( )

inline

Definition at line 36 of file PolynomialQuadrature.h.

{}

~Polynomial()

virtual PolynomialQuadrature::Polynomial::~Polynomial ( )

virtualdefault

Member Function Documentation

clenshawQuadrature()

void PolynomialQuadrature::Polynomial::clenshawQuadrature ( const unsigned int  order, std::vector< Real > &  points, std::vector< Real > &  weights  ) const

virtual

Reimplemented in PolynomialQuadrature::Legendre.

Definition at line 81 of file PolynomialQuadrature.C.

{
  ::mooseError("Clenshaw quadrature has not been implemented for this polynomial type.");
}

compute()

Real PolynomialQuadrature::Polynomial::compute ( const unsigned int  order, const Real  x, const bool  normalize = true  ) const

virtual

Reimplemented in PolynomialQuadrature::Legendre, and PolynomialQuadrature::Hermite.

Definition at line 65 of file PolynomialQuadrature.C.

{
  ::mooseError("Polynomial type has not been implemented.");
  return 0;
}

Referenced by productIntegral().

computeDerivative()

Real PolynomialQuadrature::Polynomial::computeDerivative ( const unsigned int  order, const Real  x, const unsigned int  m = 1  ) const

virtual

Computes the mth derivative of polynomial: d^mP_n/dx^m.

Reimplemented in PolynomialQuadrature::Legendre, and PolynomialQuadrature::Hermite.

Definition at line 72 of file PolynomialQuadrature.C.

{
  ::mooseError("Polynomial type has not been implemented.");
  return 0;
}

gaussQuadrature()

virtual void PolynomialQuadrature::Polynomial::gaussQuadrature ( const unsigned int  order, std::vector< Real > &  points, std::vector< Real > &  weights  ) const

pure virtual

Implemented in PolynomialQuadrature::Legendre, and PolynomialQuadrature::Hermite.

Referenced by productIntegral().

innerProduct()

virtual Real PolynomialQuadrature::Polynomial::innerProduct ( const unsigned int  order ) const

pure virtual

Implemented in PolynomialQuadrature::Legendre, and PolynomialQuadrature::Hermite.

Referenced by productIntegral().

productIntegral()

Real PolynomialQuadrature::Polynomial::productIntegral

(

const std::vector< unsigned int >

order

)

const

Definition at line 89 of file PolynomialQuadrature.C.

{
  const unsigned int nprod = order.size();
 
  if (nprod == 1)
    return (order[0] == 0 ? 1.0 : 0.0);
  else if (nprod == 2)
    return (order[0] == order[1] ? innerProduct(order[0]) : 0.0);
 
  unsigned int poly_order = std::accumulate(order.begin(), order.end(), 0);
  unsigned int quad_order = (poly_order % 2 == 0 ? poly_order : poly_order - 1) / 2;
 
  std::vector<Real> xq;
  std::vector<Real> wq;
  gaussQuadrature(quad_order, xq, wq);
 
  Real val = 0.0;
  for (unsigned int q = 0; q < xq.size(); ++q)
  {
    Real prod = wq[q];
    for (unsigned int i = 0; i < nprod; ++i)
      prod *= compute(order[i], xq[q], false);
    val += prod;
  }
 
  return val / std::accumulate(wq.begin(), wq.end(), 0.0);
}

store() [1/2]

void PolynomialQuadrature::Polynomial::store ( nlohmann::json &  json ) const

virtual

Reimplemented in PolynomialQuadrature::Legendre, and PolynomialQuadrature::Hermite.

Definition at line 58 of file PolynomialQuadrature.C.

{
  // Cannot be pure virtual because for dataLoad operations the base class must be constructed
  ::mooseError("Polynomial child class must override this method.");
}

store() [2/2]

void PolynomialQuadrature::Polynomial::store ( std::ostream &  stream, void *  context  ) const

virtual

Reimplemented in PolynomialQuadrature::Legendre, and PolynomialQuadrature::Hermite.

Definition at line 51 of file PolynomialQuadrature.C.

{
  // Cannot be pure virtual because for dataLoad operations the base class must be constructed
  ::mooseError("Polynomial child class must override this method.");
}

---

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

• stochastic_tools/include/utils/PolynomialQuadrature.h
• stochastic_tools/src/utils/PolynomialQuadrature.C