Lines Matching refs:pt
5 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data;
9 PetscCall(PetscOptionsInt("-petscspace_ptrimmed_form_degree", "form degree of trimmed space", "PetscSpacePTrimmedSetFormDegree", pt->formDegree, &pt->formDegree, NULL));
16 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data;
20 f = pt->formDegree;
40 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data;
45 if (pt->subspaces) {
48 for (d = 0; d < sp->Nv; ++d) PetscCall(PetscSpaceDestroy(&pt->subspaces[d]));
50 PetscCall(PetscFree(pt->subspaces));
51 PetscCall(PetscFree(pt));
57 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data;
61 if (pt->setupcalled) PetscFunctionReturn(PETSC_SUCCESS);
62 PetscCheck(pt->formDegree >= -sp->Nv && pt->formDegree <= sp->Nv, PetscObjectComm((PetscObject)sp), PETSC_ERR_ARG_OUTOFRANGE, "Form degree %" PetscInt_FMT " not in valid range [%" PetscInt_FMT ",%" PetscInt_FMT "]", pt->formDegree, sp->Nv, sp->Nv);
63 PetscCall(PetscDTBinomialInt(sp->Nv, PetscAbsInt(pt->formDegree), &Nf));
70 PetscInt formDegree = pt->formDegree;
102 sp->maxDegree = (pt->formDegree == 0 || PetscAbsInt(pt->formDegree) == sp->Nv) ? sp->degree : sp->degree + 1;
103 if (pt->formDegree == 0 || PetscAbsInt(pt->formDegree) == sp->Nv) {
109 pt->setupcalled = PETSC_TRUE;
115 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data;
120 f = pt->formDegree;
123 PetscCall(PetscDTPTrimmedSize(sp->Nv, f == 0 ? sp->degree : sp->degree + 1, pt->formDegree, dim));
124 PetscCall(PetscDTBinomialInt(sp->Nv, PetscAbsInt(pt->formDegree), &Nf));
136 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data;
146 if (!pt->setupcalled) {
158 f = pt->formDegree;
282 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data;
285 pt->formDegree = formDegree;
291 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data;
296 *formDegree = pt->formDegree;
302 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data;
308 if (!pt->subspaces) PetscCall(PetscCalloc1(dim, &pt->subspaces));
309 if ((dim - height) <= PetscAbsInt(pt->formDegree)) {
310 if (!pt->subspaces[height - 1]) {
316 PetscCall(PetscDTBinomialInt(dim, PetscAbsInt(pt->formDegree), &Nf));
317 PetscCall(PetscDTBinomialInt(dim - height, PetscAbsInt(pt->formDegree), &Nfsub));
328 PetscCall(PetscSpacePTrimmedSetFormDegree(sub, pt->formDegree));
330 pt->subspaces[height - 1] = sub;
332 *subsp = pt->subspaces[height - 1];
390 PetscSpace_Ptrimmed *pt;
394 PetscCall(PetscNew(&pt));
395 sp->data = pt;
397 pt->subspaces = NULL;