Home
last modified time | relevance | path

Searched refs:gll (Results 1 – 5 of 5) sorted by relevance

/petsc/src/ksp/ksp/tutorials/
H A Dex69.c22 PetscErrorCode ComputeSolution(DM da, PetscGLL *gll, Vec u) in ComputeSolution() argument
48 PetscErrorCode ComputeRhs(DM da, PetscGLL *gll, Vec b) in ComputeRhs() argument
50 PetscInt i, j, xs, xn, n = gll->n; in ComputeRhs()
69 …bb[j * (n - 1) + i] += -gll->weights[i] * (-20. * PETSC_PI * xd * PetscSinReal(5. * PETSC_PI * xd)… in ComputeRhs()
90 PetscGLL gll; in main() local
133 PetscCall(PetscMalloc2(n, &gll.nodes, n, &gll.weights)); in main()
134 …obattoLegendreQuadrature(n, PETSCGAUSSLOBATTOLEGENDRE_VIA_LINEAR_ALGEBRA, gll.nodes, gll.weights)); in main()
135 gll.n = n; in main()
136 PetscCall(DMDASetGLLCoordinates(da, gll.n, gll.nodes)); in main()
141 PetscCall(PetscGaussLobattoLegendreElementLaplacianCreate(gll.n, gll.nodes, gll.weights, &A)); in main()
[all …]
H A Dex59.c708 GLLData gll; in ComputeMatrix() local
715 PetscCall(GLLStuffs(dd, &gll)); in ComputeMatrix()
717 PetscCall(ComputeSubdomainMatrix(dd, gll, &local_mat)); in ComputeMatrix()
755 PetscCall(PetscFree(gll.zGL)); in ComputeMatrix()
756 PetscCall(PetscFree(gll.rhoGL)); in ComputeMatrix()
757 PetscCall(PetscFree(gll.A[0])); in ComputeMatrix()
758 PetscCall(PetscFree(gll.A)); in ComputeMatrix()
759 PetscCall(MatDestroy(&gll.elem_mat)); in ComputeMatrix()
/petsc/src/ts/tutorials/
H A Dex50.c62 PetscGLL gll; member
118 …PetscCall(PetscMalloc2(appctx.param.N, &appctx.SEMop.gll.nodes, appctx.param.N, &appctx.SEMop.gll.… in main()
119 ….N, PETSCGAUSSLOBATTOLEGENDRE_VIA_LINEAR_ALGEBRA, appctx.SEMop.gll.nodes, appctx.SEMop.gll.weights… in main()
120 appctx.SEMop.gll.n = appctx.param.N; 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()
232 PetscCall(PetscFree2(appctx.SEMop.gll.nodes, appctx.SEMop.gll.weights)); in main()
339 …LobattoLegendreElementLaplacianCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.… in MatMult_Laplacian()
351 …obattoLegendreElementLaplacianDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.… in MatMult_Laplacian()
[all …]
/petsc/src/tao/unconstrained/tutorials/
H A Dburgers_spectral.c70 PetscGLL gll; member
132 …PetscCall(PetscMalloc2(appctx.param.N, &appctx.SEMop.gll.nodes, appctx.param.N, &appctx.SEMop.gll.… in main()
133 ….N, PETSCGAUSSLOBATTOLEGENDRE_VIA_LINEAR_ALGEBRA, appctx.SEMop.gll.nodes, appctx.SEMop.gll.weights… in main()
134 appctx.SEMop.gll.n = appctx.param.N; in main()
176 … x = (appctx.param.Le / 2.0) * (appctx.SEMop.gll.nodes[j] + 1.0) + appctx.param.Le * i; 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()
267 PetscCall(PetscFree2(appctx.SEMop.gll.nodes, appctx.SEMop.gll.weights)); in main()
437 …LobattoLegendreElementLaplacianCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.… in RHSMatrixLaplaciangllDM()
468 …obattoLegendreElementLaplacianDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.… in RHSMatrixLaplaciangllDM()
[all …]
H A Dspectraladjointassimilation.c78 PetscGLL gll; member
143 …PetscCall(PetscMalloc2(appctx.param.N, &appctx.SEMop.gll.nodes, appctx.param.N, &appctx.SEMop.gll.… in main()
144 ….N, PETSCGAUSSLOBATTOLEGENDRE_VIA_LINEAR_ALGEBRA, appctx.SEMop.gll.nodes, appctx.SEMop.gll.weights… in main()
145 appctx.SEMop.gll.n = appctx.param.N; 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()
276 PetscCall(PetscFree2(appctx.SEMop.gll.nodes, appctx.SEMop.gll.weights)); in main()
462 …LobattoLegendreElementLaplacianCreate(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.… in RHSLaplacian()
492 …obattoLegendreElementLaplacianDestroy(appctx->SEMop.gll.n, appctx->SEMop.gll.nodes, appctx->SEMop.… in RHSLaplacian()
[all …]