PetscGaussLobattoLegendreElementAdvectionCreate#

computes the advection operator for a single 1d GLL element

Synopsis#

Not Collective

Input Parameters#

  • n - the number of GLL nodes

  • nodes - the GLL nodes

  • weights - the GLL weightss

Output Parameter#

  • AA - the stiffness element

Notes#

Destroy this with PetscGaussLobattoLegendreElementAdvectionDestroy()

This is the same as the Gradient operator multiplied by the diagonal mass matrix

You can access entries in this array with AA[i][j] but in memory it is stored in contiguous memory, row oriented

See Also#

PetscDTGaussLobattoLegendreQuadrature(), PetscGaussLobattoLegendreElementLaplacianCreate(), PetscGaussLobattoLegendreElementAdvectionDestroy()

Level#

beginner

Location#

src/dm/dt/interface/dt.c

Examples#

src/tao/unconstrained/tutorials/burgers_spectral.c
src/tao/unconstrained/tutorials/spectraladjointassimilation.c
src/ts/tutorials/ex50.c


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