PetscGaussLobattoLegendreElementAdvectionCreate#
computes the advection operator for a single 1d GLL element
Synopsis#
#include "petscdt.h"
PetscErrorCode PetscGaussLobattoLegendreElementAdvectionCreate(PetscInt n, PetscReal *nodes, PetscReal *weights, PetscReal ***AA)
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#
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