452   PetscReal **temp;  in RHSLaplacian()  local
462 …ntLaplacianCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.gll.weights, &temp));  in RHSLaplacian()
467     for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv;  in RHSLaplacian()
483 …PetscCall(MatSetValuesLocal(A, appctx->param.N, rowsDM, appctx->param.N, rowsDM, &temp[0][0], ADD_…  in RHSLaplacian()
492 …tLaplacianDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.gll.weights, &temp));  in RHSLaplacian()
503   PetscReal **temp;  in RHSAdvection()  local
513 …ntAdvectionCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.gll.weights, &temp));  in RHSAdvection()
518     for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv;  in RHSAdvection()
534 …PetscCall(MatSetValuesLocal(A, appctx->param.N, rowsDM, appctx->param.N, rowsDM, &temp[0][0], ADD_…  in RHSAdvection()
543 …tAdvectionDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.gll.weights, &temp));  in RHSAdvection()
583   Vec     temp;  in FormFunctionGradient()  local
598   PetscCall(VecDuplicate(G, &temp));  in FormFunctionGradient()
599   PetscCall(VecPointwiseMult(temp, G, G));  in FormFunctionGradient()
600   PetscCall(VecDot(temp, appctx->SEMop.mass, f));  in FormFunctionGradient()
601   PetscCall(VecDestroy(&temp));  in FormFunctionGradient()
616   Vec       temp, grad;  in MonitorError()  local
622   PetscCall(VecDuplicate(appctx->dat.ic, &temp));  in MonitorError()
623   PetscCall(VecWAXPY(temp, -1.0, appctx->dat.ic, appctx->dat.true_solution));  in MonitorError()
624   PetscCall(VecPointwiseMult(temp, temp, temp));  in MonitorError()
625   PetscCall(VecDot(temp, appctx->SEMop.mass, &nrm));  in MonitorError()
628   PetscCall(VecPointwiseMult(temp, temp, temp));  in MonitorError()
629   PetscCall(VecDot(temp, appctx->SEMop.mass, &gnorm));  in MonitorError()
631   PetscCall(VecDestroy(&temp));  in MonitorError()