PetscGaussLobattoLegendreElementLaplacianCreate#
computes the Laplacian for a single 1d GLL element
Synopsis#
#include "petscdt.h"
PetscErrorCode PetscGaussLobattoLegendreElementLaplacianCreate(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 weights
Output Parameter#
A - the stiffness element
Notes#
Destroy this with PetscGaussLobattoLegendreElementLaplacianDestroy()
You can access entries in this array with AA[i][j] but in memory it is stored in contiguous memory, row oriented (the array is symmetric)
See Also#
PetscDTGaussLobattoLegendreQuadrature(), PetscGaussLobattoLegendreElementLaplacianDestroy()
Level#
beginner
Location#
Examples#
src/ksp/ksp/tutorials/ex68.c
src/ksp/ksp/tutorials/ex69.c
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