FreeFEM Documentation on GitHub

823 stars - 195 forks

Quadrature formulae

The quadrature formula is like the following:

Df(x)=1Lωf(ξ)

int1d

Quadrature formula on an edge.

Notations

|D| is the measure of the edge D.

For a shake of simplicity, we denote:

f(x)=g(t)

with 0t1; x=(1t)x0+tx1.

qf1pE

1int1d(Th, qfe=qf1pE)( ... )

or

1int1d(Th, qforder=2)( ... )

This quadrature formula is exact on P1.

Df(x)|D|g(12)

qf2pE

1int1d(Th, qfe=qf2pE)( ... )

or

1int1d(Th, qforder=3)( ... )

This quadrature formula is exact on P3.

Df(x)|D|2(g(1+1/32)+g(11/32))

qf3pE (default)

1int1d(Th, qfe=qf3pE)( ... )

or

1int1d(Th, qforder=6)( ... )

This quadrature formula is the default one and be exact on P5.

Df(x)|D|18(5g(1+3/52)+8g(12)+5g(13/52))

qf4pE

1int1d(Th, qfe=qf4pE)( ... )

or

1int1d(Th, qforder=8)( ... )

This quadrature formula is exact on P7.

Df(x)|D|72((1830)g(1525+7030352)+(1830)g(1+525+7030352)+(18+30)g(15257030352)+(18+30)g(1+5257030352))

qf5pE

1int1d(Th, qfe=qf5pE)( ... )

or

1int1d(Th, qforder=10)( ... )

This quadrature formula is exact on P9.

Df(x)|D|((3321370)1800g(1245+1470212)+(3321370)1800g(1+245+1470212)+64225g(12)+(332+1370)1800g(12451470212)+(332+1370)1800g(1+2451470212))

qf1pElump

1int1d(Th, qfe=qf1pElump)( ... )

This quadrature formula is exact on P1.

Df(x)|D|2(g(0)+g(1))

int2d

Note

Complete formulas are no longer detailed

qf1pT

1 int2d(Th, qft=qf1pT)( ... )

or

1 int2d(Th, qforder=2)( ... )

This quadrature formula is exact on P1.

qf2pT

1 int2d(Th, qft=qf2pT)( ... )

or

1 int2d(Th, qforder=3)( ... )

This quadrature formula is exact on P2.

qf5pT (default)

1 int2d(Th, qft=qf5pT)( ... )

or

1 int2d(Th, qforder=6)( ... )

This quadrature formula is the default and be exact on P5.

qf1pTlump

1 int2d(Th, qft=qf1pTlump)( ... )

This quadrature formula is exact on P1.

qf2pT4P1

1 int2d(Th, qft=qf2pT4P1)( ... )

This quadrature formula is exact on P1.

qf7pT

1 int2d(Th, qft=qf7pT)( ... )

or

1 int2d(Th, qforder=8)( ... )

This quadrature formula is exact on P7.

qf9pT

1 int2d(Th, qft=qf9pT)( ... )

or

1 int2d(Th, qforder=10)( ... )

This quadrature formula is exact on P9.

int3d

qfV1

1 int3d(Th, qfV=qfV1)( ... )

or

1 int3d(Th, qforder=2)( ... )

This quadrature formula is exact on P1.

qfV2

1 int3d(Th, qfV=qfV2)( ... )

or

1 int3d(Th, qforder=3)( ... )

This quadrature formula is exact on P2.

qfV5 (default)

1 int3d(Th, qfV=qfV5)( ... )

or

1 int3d(Th, qforder=6)( ... )

This quadrature formula is the default one and be exact on P5.

qfV1lump

1 int3d(Th, qfV=qfV1lump)( ... )

This quadrature formula is exact on P1.

Table of content