20   DMField_DS *dsfield = (DMField_DS *)field->data;  in DMFieldDestroy_DS()  local
24   PetscCall(VecDestroy(&dsfield->vec));  in DMFieldDestroy_DS()
25   for (i = 0; i < dsfield->height; i++) PetscCall(PetscObjectDereference(dsfield->disc[i]));  in DMFieldDestroy_DS()
26   PetscCall(PetscFree(dsfield->disc));  in DMFieldDestroy_DS()
27   PetscCall(VecDestroy(&dsfield->vecDG));  in DMFieldDestroy_DS()
28   if (dsfield->discDG)  in DMFieldDestroy_DS()
29     for (i = 0; i < dsfield->height; i++) PetscCall(PetscObjectDereference(dsfield->discDG[i]));  in DMFieldDestroy_DS()
30   PetscCall(PetscFree(dsfield->discDG));  in DMFieldDestroy_DS()
31   PetscCall(PetscFree(dsfield));  in DMFieldDestroy_DS()
37   DMField_DS *dsfield = (DMField_DS *)field->data;  in DMFieldView_DS()  local
43   disc = dsfield->disc[0];  in DMFieldView_DS()
45 … PetscCall(PetscViewerASCIIPrintf(viewer, "PetscDS field %" PetscInt_FMT "\n", dsfield->fieldNum));  in DMFieldView_DS()
51 …PetscCheck(!dsfield->multifieldVec, PetscObjectComm((PetscObject)field), PETSC_ERR_SUP, "View of s…  in DMFieldView_DS()
52   PetscCall(VecView(dsfield->vec, viewer));  in DMFieldView_DS()
88   DMField_DS        *dsfield = (DMField_DS *)field->data;  in DMFieldGetClosure_Internal()  local
89   DM                 fdm     = dsfield->dmDG;  in DMFieldGetClosure_Internal()
109   if (!dsfield->vecDG) goto cg;  in DMFieldGetClosure_Internal()
111   PetscCall(VecGetArrayRead(dsfield->vecDG, array));  in DMFieldGetClosure_Internal()
115   PetscCall(VecRestoreArrayRead(dsfield->vecDG, array));  in DMFieldGetClosure_Internal()
122   PetscCall(PetscSectionGetField(s, dsfield->fieldNum, &s));  in DMFieldGetClosure_Internal()
123   PetscCall(DMPlexVecGetClosure(fdm, s, dsfield->vec, cell, Nc, values));  in DMFieldGetClosure_Internal()
129   DMField_DS  *dsfield = (DMField_DS *)field->data;  in DMFieldRestoreClosure_Internal()  local
135     PetscCall(DMRestoreWorkArray(dsfield->dmDG, *Nc, MPIU_SCALAR, values));  in DMFieldRestoreClosure_Internal()
139     PetscCall(DMPlexVecRestoreClosure(fdm, s, dsfield->vec, cell, Nc, values));  in DMFieldRestoreClosure_Internal()
147   DMField_DS      *dsfield = (DMField_DS *)field->data;  in DMFieldEvaluateFE_DS()  local
164   PetscCall(DMFieldDSGetHeightDisc(field, dsfield->height - 1 - dim, dsfield->disc, &disc));  in DMFieldEvaluateFE_DS()
167   PetscCall(PetscSectionGetField(section, dsfield->fieldNum, &section));  in DMFieldEvaluateFE_DS()
229   DMField_DS        *dsfield = (DMField_DS *)field->data;  in DMFieldEvaluate_DS()  local
252   PetscCall(DMFieldDSGetHeightDisc(field, 0, dsfield->disc, &cellDisc));  in DMFieldEvaluate_DS()
451   DMField_DS      *dsfield = (DMField_DS *)field->data;  in DMFieldEvaluateFV_DS()  local
472   for (h = 0; h < dsfield->height; h++) {  in DMFieldEvaluateFV_DS()
479   PetscCall(DMFieldDSGetHeightDisc(field, h, dsfield->disc, &disc));  in DMFieldEvaluateFV_DS()
658   DMField_DS  *dsfield;  in DMFieldGetDegree_DS()  local
664   dsfield = (DMField_DS *)field->data;  in DMFieldGetDegree_DS()
669     for (h = 0; h < dsfield->height; h++) {  in DMFieldGetDegree_DS()
676   PetscCall(DMFieldDSGetHeightDisc(field, h, dsfield->disc, &disc));  in DMFieldGetDegree_DS()
720   DMField_DS  *dsfield;  in DMFieldCreateDefaultQuadrature_DS()  local
727   dsfield = (DMField_DS *)field->data;  in DMFieldCreateDefaultQuadrature_DS()
739   if (h < dsfield->height) {  in DMFieldCreateDefaultQuadrature_DS()
740     PetscCall(DMFieldDSGetHeightDisc(field, h, dsfield->disc, &disc));  in DMFieldCreateDefaultQuadrature_DS()
754   DMField_DS  *dsfield;  in DMFieldCreateDefaultFaceQuadrature_DS()  local
761   dsfield = (DMField_DS *)field->data;  in DMFieldCreateDefaultFaceQuadrature_DS()
773   if (h < dsfield->height) {  in DMFieldCreateDefaultFaceQuadrature_DS()
776     PetscCall(DMFieldDSGetHeightDisc(field, h, dsfield->disc, &disc));  in DMFieldCreateDefaultFaceQuadrature_DS()
904     DMField_DS    *dsfield = (DMField_DS *)field->data;  in DMFieldComputeFaceData_DS()  local
923     PetscCall(DMFieldDSGetHeightDisc(field, 1, dsfield->disc, &faceDisc));  in DMFieldComputeFaceData_DS()
924     PetscCall(DMFieldDSGetHeightDisc(field, 0, dsfield->disc, &cellDisc));  in DMFieldComputeFaceData_DS()
1148   DMField_DS *dsfield;  in DMFieldCreate_DS()  local
1151   PetscCall(PetscNew(&dsfield));  in DMFieldCreate_DS()
1152   field->data = dsfield;  in DMFieldCreate_DS()
1160   DMField_DS  *dsfield;  in DMFieldCreateDSWithDG()  local
1200   dsfield           = (DMField_DS *)b->data;  in DMFieldCreateDSWithDG()
1201   dsfield->fieldNum = fieldNum;  in DMFieldCreateDSWithDG()
1202   PetscCall(DMGetDimension(dm, &dsfield->height));  in DMFieldCreateDSWithDG()
1203   dsfield->height++;  in DMFieldCreateDSWithDG()
1204   PetscCall(PetscCalloc1(dsfield->height, &dsfield->disc));  in DMFieldCreateDSWithDG()
1205   dsfield->disc[0] = disc;  in DMFieldCreateDSWithDG()
1207   dsfield->vec = vec;  in DMFieldCreateDSWithDG()
1209     dsfield->dmDG = dmDG;  in DMFieldCreateDSWithDG()
1210     PetscCall(PetscCalloc1(dsfield->height, &dsfield->discDG));  in DMFieldCreateDSWithDG()
1211     dsfield->discDG[0] = discDG;  in DMFieldCreateDSWithDG()
1213     dsfield->vecDG = vecDG;  in DMFieldCreateDSWithDG()