324   PetscReal        **temp, vv;  in MatMult_Laplacian()  local
339 …ntLaplacianCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.gll.weights, &temp));  in MatMult_Laplacian()
342     for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv;  in MatMult_Laplacian()
348 …x->param.N - 1) PetscCallBLAS("BLASgemv", BLASgemv_("N", &n, &n, &_DOne, &temp[0][0], &n, &xl[j], …  in MatMult_Laplacian()
351 …tLaplacianDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.gll.weights, &temp));  in MatMult_Laplacian()
364   PetscReal        **temp;  in MatMult_Advection()  local
379 …ntAdvectionCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.gll.weights, &temp));  in MatMult_Advection()
384 …x->param.N - 1) PetscCallBLAS("BLASgemv", BLASgemv_("N", &n, &n, &_DOne, &temp[0][0], &n, &xl[j], …  in MatMult_Advection()
387 …tAdvectionDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.gll.weights, &temp));  in MatMult_Advection()
415   PetscReal **temp;  in RHSMatrixLaplaciangllDM()  local
429 …ntLaplacianCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.gll.weights, &temp));  in RHSMatrixLaplaciangllDM()
436       for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv;  in RHSMatrixLaplaciangllDM()
451 …PetscCall(MatSetValuesLocal(A, appctx->param.N, rowsDM, appctx->param.N, rowsDM, &temp[0][0], ADD_…  in RHSMatrixLaplaciangllDM()
460 …tLaplacianDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.gll.weights, &temp));  in RHSMatrixLaplaciangllDM()
487   PetscReal **temp;  in RHSMatrixAdvectiongllDM()  local
500 …ntAdvectionCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.gll.weights, &temp));  in RHSMatrixAdvectiongllDM()
509 …PetscCall(MatSetValuesLocal(A, appctx->param.N, rowsDM, appctx->param.N, rowsDM, &temp[0][0], ADD_…  in RHSMatrixAdvectiongllDM()
519 …tAdvectionDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.gll.weights, &temp));  in RHSMatrixAdvectiongllDM()