ex50.c - OpenGrok cross reference for /petsc/src/ts/tutorials/ex50.c

Lines Matching refs:SEMop
67   PetscSEMOperators SEMop;  member
118 …PetscCall(PetscMalloc2(appctx.param.N, &appctx.SEMop.gll.nodes, appctx.param.N, &appctx.SEMop.gll.…  in main()
119 …param.N, PETSCGAUSSLOBATTOLEGENDRE_VIA_LINEAR_ALGEBRA, appctx.SEMop.gll.nodes, appctx.SEMop.gll.we…  in main()
120   appctx.SEMop.gll.n = appctx.param.N;  in main()
140   PetscCall(VecDuplicate(appctx.dat.curr_sol, &appctx.SEMop.grid));  in main()
141   PetscCall(VecDuplicate(appctx.dat.curr_sol, &appctx.SEMop.mass));  in main()
144   PetscCall(DMDAVecGetArray(appctx.da, appctx.SEMop.grid, &wrk_ptr1));  in main()
145   PetscCall(DMDAVecGetArray(appctx.da, appctx.SEMop.mass, &wrk_ptr2));  in main()
158 … x             = (appctx.param.Le / 2.0) * (appctx.SEMop.gll.nodes[j] + 1.0) + appctx.param.Le * i;  in main()
161       wrk_ptr2[ind] = .5 * appctx.param.Le * appctx.SEMop.gll.weights[j];  in main()
162       if (j == 0) wrk_ptr2[ind] += .5 * appctx.param.Le * appctx.SEMop.gll.weights[j];  in main()
165   PetscCall(DMDAVecRestoreArray(appctx.da, appctx.SEMop.grid, &wrk_ptr1));  in main()
166   PetscCall(DMDAVecRestoreArray(appctx.da, appctx.SEMop.mass, &wrk_ptr2));  in main()
172   PetscCall(DMCreateMatrix(appctx.da, &appctx.SEMop.stiff));  in main()
173   PetscCall(DMCreateMatrix(appctx.da, &appctx.SEMop.grad));  in main()
179 …atrixLaplaciangllDM(appctx.ts, 0.0, appctx.dat.curr_sol, appctx.SEMop.stiff, appctx.SEMop.stiff, &…  in main()
180 …atrixAdvectiongllDM(appctx.ts, 0.0, appctx.dat.curr_sol, appctx.SEMop.grad, appctx.SEMop.grad, &ap…  in main()
187   PetscCall(MatDuplicate(appctx.SEMop.stiff, MAT_COPY_VALUES, &appctx.SEMop.keptstiff));  in main()
191   PetscCall(MatSetNullSpace(appctx.SEMop.stiff, nsp));  in main()
192   PetscCall(MatSetNullSpace(appctx.SEMop.keptstiff, nsp));  in main()
193   PetscCall(MatNullSpaceTest(nsp, appctx.SEMop.stiff, NULL));  in main()
197   PetscCall(MatSetNullSpace(appctx.SEMop.grad, nsp));  in main()
198   PetscCall(MatNullSpaceTest(nsp, appctx.SEMop.grad, NULL));  in main()
215 …PetscCall(TSSetRHSJacobian(appctx.ts, appctx.SEMop.stiff, appctx.SEMop.stiff, RHSJacobian, &appctx…  in main()
226   PetscCall(MatDestroy(&appctx.SEMop.stiff));  in main()
227   PetscCall(MatDestroy(&appctx.SEMop.keptstiff));  in main()
228   PetscCall(MatDestroy(&appctx.SEMop.grad));  in main()
229   PetscCall(VecDestroy(&appctx.SEMop.grid));  in main()
230   PetscCall(VecDestroy(&appctx.SEMop.mass));  in main()
232   PetscCall(PetscFree2(appctx.SEMop.gll.nodes, appctx.SEMop.gll.weights));  in main()
264   PetscCall(DMDAVecGetArrayRead(appctx->da, appctx->SEMop.grid, (void *)&xg));  in TrueSolution()
270   PetscCall(DMDAVecRestoreArrayRead(appctx->da, appctx->SEMop.grid, (void *)&xg));  in TrueSolution()
279   PetscCall(MatMult(appctx->SEMop.grad, globalin, globalout)); /* grad u */  in RHSFunction()
282   PetscCall(MatMultAdd(appctx->SEMop.keptstiff, globalin, globalout, globalout));  in RHSFunction()
302   PetscCall(MatCopy(appctx->SEMop.grad, A, SAME_NONZERO_PATTERN));  in RHSJacobian()
307   PetscCall(MatMult(appctx->SEMop.grad, globalin, Gglobalin));  in RHSJacobian()
313   PetscCall(MatAXPY(A, 0.0, appctx->SEMop.keptstiff, SAME_NONZERO_PATTERN));  in RHSJacobian()
339 …GaussLobattoLegendreElementLaplacianCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->S…  in MatMult_Laplacian()
351 …aussLobattoLegendreElementLaplacianDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->S…  in MatMult_Laplacian()
357   PetscCall(VecPointwiseDivide(y, y, appctx->SEMop.mass));  in MatMult_Laplacian()
379 …GaussLobattoLegendreElementAdvectionCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->S…  in MatMult_Advection()
387 …aussLobattoLegendreElementAdvectionDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->S…  in MatMult_Advection()
393   PetscCall(VecPointwiseDivide(y, y, appctx->SEMop.mass));  in MatMult_Advection()
429 …GaussLobattoLegendreElementLaplacianCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->S…  in RHSMatrixLaplaciangllDM()
456     PetscCall(VecReciprocal(appctx->SEMop.mass));  in RHSMatrixLaplaciangllDM()
457     PetscCall(MatDiagonalScale(A, appctx->SEMop.mass, 0));  in RHSMatrixLaplaciangllDM()
458     PetscCall(VecReciprocal(appctx->SEMop.mass));  in RHSMatrixLaplaciangllDM()
460 …aussLobattoLegendreElementLaplacianDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->S…  in RHSMatrixLaplaciangllDM()
500 …GaussLobattoLegendreElementAdvectionCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->S…  in RHSMatrixAdvectiongllDM()
515     PetscCall(VecReciprocal(appctx->SEMop.mass));  in RHSMatrixAdvectiongllDM()
516     PetscCall(MatDiagonalScale(A, appctx->SEMop.mass, 0));  in RHSMatrixAdvectiongllDM()
517     PetscCall(VecReciprocal(appctx->SEMop.mass));  in RHSMatrixAdvectiongllDM()
519 …aussLobattoLegendreElementAdvectionDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->S…  in RHSMatrixAdvectiongllDM()