Lines Matching refs:grid_aniso_proj
12 NodalProjectionData grid_aniso_proj; in GridAnisotropyTensorProjectionSetupApply() local
21 PetscCall(PetscNew(&grid_aniso_proj)); in GridAnisotropyTensorProjectionSetupApply()
26 PetscCall(DMClone(honee->dm, &grid_aniso_proj->dm)); in GridAnisotropyTensorProjectionSetupApply()
27 PetscCall(DMSetMatrixPreallocateSkip(grid_aniso_proj->dm, PETSC_TRUE)); in GridAnisotropyTensorProjectionSetupApply()
28 …PetscCall(PetscObjectSetName((PetscObject)grid_aniso_proj->dm, "Grid Anisotropy Tensor Projection"… in GridAnisotropyTensorProjectionSetupApply()
31 grid_aniso_proj->num_comp = 7; in GridAnisotropyTensorProjectionSetupApply()
32 …RUE, PETSC_TRUE, honee->app_ctx->degree, 1, honee->app_ctx->q_extra, 1, &grid_aniso_proj->num_comp, in GridAnisotropyTensorProjectionSetupApply()
33 grid_aniso_proj->dm)); in GridAnisotropyTensorProjectionSetupApply()
35 PetscCall(DMGetLocalSection(grid_aniso_proj->dm, §ion)); in GridAnisotropyTensorProjectionSetupApply()
47 …PetscCall(DMPlexCeedElemRestrictionCreate(ceed, grid_aniso_proj->dm, DMLABEL_DEFAULT, DMLABEL_DEFA… in GridAnisotropyTensorProjectionSetupApply()
49 …PetscCall(DMPlexCeedBasisCreate(ceed, grid_aniso_proj->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE,… in GridAnisotropyTensorProjectionSetupApply()
50 …PetscCall(QDataGet(ceed, grid_aniso_proj->dm, DMLABEL_DEFAULT, DMLABEL_DEFAULT_VALUE, &elem_restr_… in GridAnisotropyTensorProjectionSetupApply()
58 …PetscCallCeed(ceed, CeedQFunctionAddOutput(qf_rhs_assemble, "v", grid_aniso_proj->num_comp, CEED_E… in GridAnisotropyTensorProjectionSetupApply()
64 …PetscCall(OperatorApplyContextCreate(honee->dm, grid_aniso_proj->dm, ceed, op_rhs_assemble, CEED_V… in GridAnisotropyTensorProjectionSetupApply()
65 &grid_aniso_proj->l2_rhs_ctx)); in GridAnisotropyTensorProjectionSetupApply()
76 PetscCall(HoneeMassQFunctionCreate(ceed, grid_aniso_proj->num_comp, q_data_size, &qf_mass)); in GridAnisotropyTensorProjectionSetupApply()
82 PetscCall(MatCreateCeed(grid_aniso_proj->dm, grid_aniso_proj->dm, op_mass, NULL, &mat_mass)); in GridAnisotropyTensorProjectionSetupApply()
84 PetscCall(KSPCreate(comm, &grid_aniso_proj->ksp)); in GridAnisotropyTensorProjectionSetupApply()
85 PetscCall(KSPSetOptionsPrefix(grid_aniso_proj->ksp, "grid_anisotropy_tensor_projection_")); in GridAnisotropyTensorProjectionSetupApply()
88 PetscCall(KSPGetPC(grid_aniso_proj->ksp, &pc)); in GridAnisotropyTensorProjectionSetupApply()
91 PetscCall(KSPSetType(grid_aniso_proj->ksp, KSPCG)); in GridAnisotropyTensorProjectionSetupApply()
92 PetscCall(KSPSetNormType(grid_aniso_proj->ksp, KSP_NORM_NATURAL)); in GridAnisotropyTensorProjectionSetupApply()
93 …PetscCall(KSPSetTolerances(grid_aniso_proj->ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, PETSC_DEFAUL… in GridAnisotropyTensorProjectionSetupApply()
95 PetscCall(KSPSetFromOptions_WithMatCeed(grid_aniso_proj->ksp, mat_mass)); in GridAnisotropyTensorProjectionSetupApply()
105 PetscCall(DMGetGlobalVector(grid_aniso_proj->dm, &Grid_Anisotropy)); in GridAnisotropyTensorProjectionSetupApply()
106 PetscCall(ApplyCeedOperatorLocalToGlobal(NULL, Grid_Anisotropy, grid_aniso_proj->l2_rhs_ctx)); in GridAnisotropyTensorProjectionSetupApply()
107 PetscCall(KSPSolve(grid_aniso_proj->ksp, Grid_Anisotropy, Grid_Anisotropy)); in GridAnisotropyTensorProjectionSetupApply()
110 PetscCall(DMGetLocalVector(grid_aniso_proj->dm, &grid_anisotropy_loc)); in GridAnisotropyTensorProjectionSetupApply()
112 …PetscCall(DMGlobalToLocal(grid_aniso_proj->dm, Grid_Anisotropy, INSERT_VALUES, grid_anisotropy_loc… in GridAnisotropyTensorProjectionSetupApply()
114 PetscCall(DMRestoreLocalVector(grid_aniso_proj->dm, &grid_anisotropy_loc)); in GridAnisotropyTensorProjectionSetupApply()
115 PetscCall(DMRestoreGlobalVector(grid_aniso_proj->dm, &Grid_Anisotropy)); in GridAnisotropyTensorProjectionSetupApply()
118 PetscCall(NodalProjectionDataDestroy(&grid_aniso_proj)); in GridAnisotropyTensorProjectionSetupApply()