83   PetscSEMOperators SEMop;  member
143 …PetscCall(PetscMalloc2(appctx.param.N, &appctx.SEMop.gll.nodes, appctx.param.N, &appctx.SEMop.gll.…  in main()
144 …param.N, PETSCGAUSSLOBATTOLEGENDRE_VIA_LINEAR_ALGEBRA, appctx.SEMop.gll.nodes, appctx.SEMop.gll.we…  in main()
145   appctx.SEMop.gll.n = appctx.param.N;  in main()
168   PetscCall(VecDuplicate(u, &appctx.SEMop.grid));  in main()
169   PetscCall(VecDuplicate(u, &appctx.SEMop.mass));  in main()
174   PetscCall(DMDAVecGetArray(appctx.da, appctx.SEMop.grid, &wrk_ptr1));  in main()
175   PetscCall(DMDAVecGetArray(appctx.da, appctx.SEMop.mass, &wrk_ptr2));  in main()
188 … x             = (appctx.param.Le / 2.0) * (appctx.SEMop.gll.nodes[j] + 1.0) + appctx.param.Le * i;  in main()
191       wrk_ptr2[ind] = .5 * appctx.param.Le * appctx.SEMop.gll.weights[j];  in main()
192       if (j == 0) wrk_ptr2[ind] += .5 * appctx.param.Le * appctx.SEMop.gll.weights[j];  in main()
195   PetscCall(DMDAVecRestoreArray(appctx.da, appctx.SEMop.grid, &wrk_ptr1));  in main()
196   PetscCall(DMDAVecRestoreArray(appctx.da, appctx.SEMop.mass, &wrk_ptr2));  in main()
202   PetscCall(DMCreateMatrix(appctx.da, &appctx.SEMop.stiff));  in main()
203   PetscCall(DMCreateMatrix(appctx.da, &appctx.SEMop.advec));  in main()
210   PetscCall(RHSLaplacian(appctx.ts, 0.0, u, appctx.SEMop.stiff, appctx.SEMop.stiff, &appctx));  in main()
211   PetscCall(RHSAdvection(appctx.ts, 0.0, u, appctx.SEMop.advec, appctx.SEMop.advec, &appctx));  in main()
212   PetscCall(MatAXPY(appctx.SEMop.stiff, -1.0, appctx.SEMop.advec, DIFFERENT_NONZERO_PATTERN));  in main()
213   PetscCall(MatDuplicate(appctx.SEMop.stiff, MAT_COPY_VALUES, &appctx.SEMop.keptstiff));  in main()
217   PetscCall(MatSetNullSpace(appctx.SEMop.stiff, nsp));  in main()
218   PetscCall(MatNullSpaceTest(nsp, appctx.SEMop.stiff, NULL));  in main()
237 …PetscCall(TSSetRHSJacobian(appctx.ts, appctx.SEMop.stiff, appctx.SEMop.stiff, TSComputeRHSJacobian…  in main()
265   PetscCall(MatDestroy(&appctx.SEMop.advec));  in main()
266   PetscCall(MatDestroy(&appctx.SEMop.stiff));  in main()
267   PetscCall(MatDestroy(&appctx.SEMop.keptstiff));  in main()
273   PetscCall(VecDestroy(&appctx.SEMop.grid));  in main()
274   PetscCall(VecDestroy(&appctx.SEMop.mass));  in main()
276   PetscCall(PetscFree2(appctx.SEMop.gll.nodes, appctx.SEMop.gll.weights));  in main()
330   PetscCall(DMDAVecGetArrayRead(appctx->da, appctx->SEMop.grid, (void *)&xg));  in InitialConditions()
341   PetscCall(DMDAVecRestoreArrayRead(appctx->da, appctx->SEMop.grid, (void *)&xg));  in InitialConditions()
369   PetscCall(DMDAVecGetArrayRead(appctx->da, appctx->SEMop.grid, (void *)&xg));  in TrueSolution()
376   PetscCall(DMDAVecRestoreArrayRead(appctx->da, appctx->SEMop.grid, (void *)&xg));  in TrueSolution()
400   PetscCall(DMDAVecGetArrayRead(appctx->da, appctx->SEMop.grid, (void *)&xg));  in ComputeReference()
410   PetscCall(DMDAVecRestoreArrayRead(appctx->da, appctx->SEMop.grid, (void *)&xg));  in ComputeReference()
419   PetscCall(MatMult(appctx->SEMop.keptstiff, globalin, globalout));  in RHSFunction()
428   PetscCall(MatCopy(appctx->SEMop.keptstiff, A, DIFFERENT_NONZERO_PATTERN));  in RHSJacobian()
462 …GaussLobattoLegendreElementLaplacianCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->S…  in RHSLaplacian()
488   PetscCall(VecReciprocal(appctx->SEMop.mass));  in RHSLaplacian()
489   PetscCall(MatDiagonalScale(A, appctx->SEMop.mass, 0));  in RHSLaplacian()
490   PetscCall(VecReciprocal(appctx->SEMop.mass));  in RHSLaplacian()
492 …aussLobattoLegendreElementLaplacianDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->S…  in RHSLaplacian()
513 …GaussLobattoLegendreElementAdvectionCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->S…  in RHSAdvection()
539   PetscCall(VecReciprocal(appctx->SEMop.mass));  in RHSAdvection()
540   PetscCall(MatDiagonalScale(A, appctx->SEMop.mass, 0));  in RHSAdvection()
541   PetscCall(VecReciprocal(appctx->SEMop.mass));  in RHSAdvection()
543 …aussLobattoLegendreElementAdvectionDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->S…  in RHSAdvection()
600   PetscCall(VecDot(temp, appctx->SEMop.mass, f));  in FormFunctionGradient()
605   PetscCall(VecPointwiseMult(G, G, appctx->SEMop.mass));  in FormFunctionGradient()
625   PetscCall(VecDot(temp, appctx->SEMop.mass, &nrm));  in MonitorError()
629   PetscCall(VecDot(temp, appctx->SEMop.mass, &gnorm));  in MonitorError()