PetscErrorCode PetscGLLElementAdvectionCreate(PetscGLL *gll,PetscReal ***AA)Not Collective
Notes: Destroy this with PetscGLLElementAdvectionDestroy()
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