PetscDTGaussLobattoJacobiQuadrature#

quadrature for the interval [a, b] with the weight function \((x-a)^\alpha (x-b)^\beta\), with endpoints a and b included as quadrature points.

Synopsis#

#include "petscdt.h" 
PetscErrorCode PetscDTGaussLobattoJacobiQuadrature(PetscInt npoints, PetscReal a, PetscReal b, PetscReal alpha, PetscReal beta, PetscReal x[], PetscReal w[])

Not Collective

Input Parameters#

  • npoints - the number of points in the quadrature rule

  • a - the left endpoint of the interval

  • b - the right endpoint of the interval

  • alpha - the left exponent

  • beta - the right exponent

Output Parameters#

  • x - array of length npoints, the locations of the quadrature points

  • w - array of length npoints, the weights of the quadrature points

Note#

This quadrature rule is exact for polynomials up to degree 2*npoints - 3.

See Also#

PetscDTGaussJacobiQuadrature()

Level#

intermediate

Location#

src/dm/dt/interface/dt.c

Edit on GitLab

Index of all DT routines
Table of Contents for all manual pages
Index of all manual pages