76   PetscParam        param;  member
112   appctx.param.N     = 10;   /* order of the spectral element */  in main()
113   appctx.param.E     = 10;   /* number of elements */  in main()
114   appctx.param.L     = 4.0;  /* length of the domain */  in main()
115   appctx.param.mu    = 0.01; /* diffusion coefficient */  in main()
117   appctx.param.steps = PETSC_INT_MAX;  in main()
118   appctx.param.Tend  = 4;  in main()
120   PetscCall(PetscOptionsGetInt(NULL, NULL, "-N", &appctx.param.N, NULL));  in main()
121   PetscCall(PetscOptionsGetInt(NULL, NULL, "-E", &appctx.param.E, NULL));  in main()
122   PetscCall(PetscOptionsGetReal(NULL, NULL, "-Tend", &appctx.param.Tend, NULL));  in main()
123   PetscCall(PetscOptionsGetReal(NULL, NULL, "-mu", &appctx.param.mu, NULL));  in main()
124   appctx.param.Le = appctx.param.L / appctx.param.E;  in main()
127 …PetscCheck((appctx.param.E % size) == 0, PETSC_COMM_WORLD, PETSC_ERR_ARG_WRONG, "Number of element…  in main()
132 …PetscCall(PetscMalloc2(appctx.param.N, &appctx.SEMop.gll.nodes, appctx.param.N, &appctx.SEMop.gll.…  in main()
133 …PetscCall(PetscDTGaussLobattoLegendreQuadrature(appctx.param.N, PETSCGAUSSLOBATTOLEGENDRE_VIA_LINE…  in main()
134   appctx.SEMop.gll.n = appctx.param.N;  in main()
135   lenglob            = appctx.param.E * (appctx.param.N - 1);  in main()
167   xs = xs / (appctx.param.N - 1);  in main()
168   xm = xm / (appctx.param.N - 1);  in main()
175     for (j = 0; j < appctx.param.N - 1; j++) {  in main()
176 … x             = (appctx.param.Le / 2.0) * (appctx.SEMop.gll.nodes[j] + 1.0) + appctx.param.Le * i;  in main()
177       ind           = i * (appctx.param.N - 1) + j;  in main()
179       wrk_ptr2[ind] = .5 * appctx.param.Le * appctx.SEMop.gll.weights[j];  in main()
180       if (j == 0) wrk_ptr2[ind] += .5 * appctx.param.Le * appctx.SEMop.gll.weights[j];  in main()
226   PetscCall(TSSetMaxSteps(appctx.ts, appctx.param.steps));  in main()
227   PetscCall(TSSetMaxTime(appctx.ts, appctx.param.Tend));  in main()
239   PetscCall(ComputeObjective(appctx.param.Tend, appctx.dat.obj, &appctx));  in main()
304 …for (i = xs; i < xs + xn; i++) s[i] = 2.0 * appctx->param.mu * PETSC_PI * PetscSinScalar(PETSC_PI …  in InitialConditions()
332 …for (i = xs; i < xs + xn; i++) s[i] = 2.0 * appctx->param.mu * PETSC_PI * PetscSinScalar(PETSC_PI …  in TrueSolution()
358 …param.mu * PETSC_PI * PetscSinScalar(PETSC_PI * xg[i]) * PetscExpScalar(-PETSC_PI * PETSC_PI * t *…  in ComputeObjective()
439 …PetscCheck(appctx->param.N > 1, PETSC_COMM_WORLD, PETSC_ERR_ARG_WRONG, "Spectral element order sho…  in RHSMatrixLaplaciangllDM()
442   for (i = 0; i < appctx->param.N; i++) {  in RHSMatrixLaplaciangllDM()
443     vv = -appctx->param.mu * 2.0 / appctx->param.Le;  in RHSMatrixLaplaciangllDM()
444     for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv;  in RHSMatrixLaplaciangllDM()
450   xs = xs / (appctx->param.N - 1);  in RHSMatrixLaplaciangllDM()
451   xn = xn / (appctx->param.N - 1);  in RHSMatrixLaplaciangllDM()
453   PetscCall(PetscMalloc1(appctx->param.N, &rowsDM));  in RHSMatrixLaplaciangllDM()
458     for (l = 0; l < appctx->param.N; l++) rowsDM[l] = 1 + (j - xs) * (appctx->param.N - 1) + l;  in RHSMatrixLaplaciangllDM()
459 …PetscCall(MatSetValuesLocal(A, appctx->param.N, rowsDM, appctx->param.N, rowsDM, &temp[0][0], ADD_…  in RHSMatrixLaplaciangllDM()
503   xs = xs / (appctx->param.N - 1);  in RHSMatrixAdvectiongllDM()
504   xn = xn / (appctx->param.N - 1);  in RHSMatrixAdvectiongllDM()
506   PetscCall(PetscMalloc1(appctx->param.N, &rowsDM));  in RHSMatrixAdvectiongllDM()
508     for (l = 0; l < appctx->param.N; l++) rowsDM[l] = 1 + (j - xs) * (appctx->param.N - 1) + l;  in RHSMatrixAdvectiongllDM()
509 …PetscCall(MatSetValuesLocal(A, appctx->param.N, rowsDM, appctx->param.N, rowsDM, &temp[0][0], ADD_…  in RHSMatrixAdvectiongllDM()