PetscDTGaussQuadrature#

create Gauss-Legendre quadrature

Synopsis#

#include "petscdt.h" 
PetscErrorCode PetscDTGaussQuadrature(PetscInt npoints, PetscReal a, PetscReal b, PetscReal *x, PetscReal *w)

Not Collective

Input Parameters#

  • npoints - number of points

  • a - left end of interval (often-1)

  • b - right end of interval (often +1)

Output Parameters#

  • x - quadrature points

  • w - quadrature weights

References#

  • **** -*** Golub and Welsch, Calculation of Quadrature Rules, Math. Comp. 23(106), 1969.

See Also#

PetscDTLegendreEval(), PetscDTGaussJacobiQuadrature()

Level#

intermediate

Location#

src/dm/dt/interface/dt.c

Examples#

src/ksp/ksp/tutorials/ex74.c
src/snes/tutorials/ex31.c

Edit on GitLab

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