137   DMCeed            sd;  in DMCeedCreateGeometry()  local
147   PetscCall(PetscCalloc1(1, &sd));  in DMCeedCreateGeometry()
156   PetscCall(PetscFEGetCeedBasis(fe, &sd->basis));  in DMCeedCreateGeometry()
157   PetscCall(CeedBasisGetNumQuadraturePoints(sd->basis, &Nq));  in DMCeedCreateGeometry()
158   PetscCall(DMPlexGetCeedRestriction(dm, NULL, 0, 0, 0, &sd->er));  in DMCeedCreateGeometry()
173   PetscCallCEED(CeedQFunctionCreateInterior(ceed, 1, geom, geomName, &sd->qf));  in DMCeedCreateGeometry()
174   PetscCallCEED(CeedQFunctionAddInput(sd->qf, "x", cdim, CEED_EVAL_INTERP));  in DMCeedCreateGeometry()
175   PetscCallCEED(CeedQFunctionAddInput(sd->qf, "dx", cdim * dim, CEED_EVAL_GRAD));  in DMCeedCreateGeometry()
176   PetscCallCEED(CeedQFunctionAddInput(sd->qf, "weight", 1, CEED_EVAL_WEIGHT));  in DMCeedCreateGeometry()
177   PetscCallCEED(CeedQFunctionAddOutput(sd->qf, "qdata", *Nqdata, CEED_EVAL_NONE));  in DMCeedCreateGeometry()
179 …PetscCallCEED(CeedOperatorCreate(ceed, sd->qf, CEED_QFUNCTION_NONE, CEED_QFUNCTION_NONE, &sd->op));  in DMCeedCreateGeometry()
180   PetscCallCEED(CeedOperatorSetField(sd->op, "x", sd->er, sd->basis, CEED_VECTOR_ACTIVE));  in DMCeedCreateGeometry()
181   PetscCallCEED(CeedOperatorSetField(sd->op, "dx", sd->er, sd->basis, CEED_VECTOR_ACTIVE));  in DMCeedCreateGeometry()
182 …PetscCallCEED(CeedOperatorSetField(sd->op, "weight", CEED_ELEMRESTRICTION_NONE, sd->basis, CEED_VE…  in DMCeedCreateGeometry()
183   PetscCallCEED(CeedOperatorSetField(sd->op, "qdata", *erq, CEED_BASIS_NONE, CEED_VECTOR_ACTIVE));  in DMCeedCreateGeometry()
186   *soldata = sd;  in DMCeedCreateGeometry()
201   DMCeed   sd;  in DMCeedCreate_Internal()  local
207   PetscCall(PetscCalloc1(1, &sd));  in DMCeedCreate_Internal()
212   PetscCall(PetscFEGetCeedBasis(fe, &sd->basis));  in DMCeedCreate_Internal()
214   PetscCall(CeedBasisGetNumQuadraturePoints(sd->basis, &Nq));  in DMCeedCreate_Internal()
215   PetscCall(DMPlexGetCeedRestriction(dm, NULL, 0, 0, 0, &sd->er));  in DMCeedCreate_Internal()
221     PetscCall(DMCeedCreateGeometry(cdm, cellIS, &Nqdata, &sd->erq, &sd->qd, &sd->geom));  in DMCeedCreate_Internal()
224   if (sd->geom) {  in DMCeedCreate_Internal()
228     PetscCallCEED(CeedBasisGetNumQuadraturePoints(sd->geom->basis, &Nqx));  in DMCeedCreate_Internal()
235   PetscCallCEED(CeedQFunctionCreateInterior(ceed, 1, func, func_source, &sd->qf));  in DMCeedCreate_Internal()
236   PetscCallCEED(CeedQFunctionAddInput(sd->qf, "u", Nc, CEED_EVAL_INTERP));  in DMCeedCreate_Internal()
237   PetscCallCEED(CeedQFunctionAddInput(sd->qf, "du", Nc * dim, CEED_EVAL_GRAD));  in DMCeedCreate_Internal()
238   PetscCallCEED(CeedQFunctionAddInput(sd->qf, "qdata", Nqdata, CEED_EVAL_NONE));  in DMCeedCreate_Internal()
239   PetscCallCEED(CeedQFunctionAddOutput(sd->qf, "v", Nc, CEED_EVAL_INTERP));  in DMCeedCreate_Internal()
240   PetscCallCEED(CeedQFunctionAddOutput(sd->qf, "dv", Nc * dim, CEED_EVAL_GRAD));  in DMCeedCreate_Internal()
242 …PetscCallCEED(CeedOperatorCreate(ceed, sd->qf, CEED_QFUNCTION_NONE, CEED_QFUNCTION_NONE, &sd->op));  in DMCeedCreate_Internal()
243   PetscCallCEED(CeedOperatorSetField(sd->op, "u", sd->er, sd->basis, CEED_VECTOR_ACTIVE));  in DMCeedCreate_Internal()
244   PetscCallCEED(CeedOperatorSetField(sd->op, "du", sd->er, sd->basis, CEED_VECTOR_ACTIVE));  in DMCeedCreate_Internal()
245   PetscCallCEED(CeedOperatorSetField(sd->op, "qdata", sd->erq, CEED_BASIS_NONE, sd->qd));  in DMCeedCreate_Internal()
246   PetscCallCEED(CeedOperatorSetField(sd->op, "v", sd->er, sd->basis, CEED_VECTOR_ACTIVE));  in DMCeedCreate_Internal()
247   PetscCallCEED(CeedOperatorSetField(sd->op, "dv", sd->er, sd->basis, CEED_VECTOR_ACTIVE));  in DMCeedCreate_Internal()
250   sd->func = func;  in DMCeedCreate_Internal()
251   PetscCall(PetscStrallocpy(func_source, &sd->funcSource));  in DMCeedCreate_Internal()
254   *soldata = sd;  in DMCeedCreate_Internal()
277   DMCeed          sd;  in DMCeedCreateGeometryFVM()  local
283   PetscCall(PetscCalloc1(1, &sd));  in DMCeedCreateGeometryFVM()
296   *soldata = sd;  in DMCeedCreateGeometryFVM()
329   DMCeed   sd;  in DMCeedCreateFVM_Internal()  local
334   PetscCall(PetscCalloc1(1, &sd));  in DMCeedCreateFVM_Internal()
340   PetscCall(DMPlexCreateCeedRestrictionFVM(dm, &sd->erL, &sd->erR));  in DMCeedCreateFVM_Internal()
346     PetscCall(DMCeedCreateGeometryFVM(cdm, faceIS, &Nqdata, &sd->erq, &sd->qd, &sd->geom));  in DMCeedCreateFVM_Internal()
353     PetscCall(DMCeedCreateInfoFVM(cdm, faceIS, &Nqinfo, &sd->eri, &sd->qi, &sd->info));  in DMCeedCreateFVM_Internal()
354     PetscCall(DMCeedComputeInfo(dm, sd));  in DMCeedCreateFVM_Internal()
357   PetscCallCEED(CeedQFunctionCreateInterior(ceed, 1, func, func_source, &sd->qf));  in DMCeedCreateFVM_Internal()
358   PetscCallCEED(CeedQFunctionAddInput(sd->qf, "uL", Nc, CEED_EVAL_NONE));  in DMCeedCreateFVM_Internal()
359   PetscCallCEED(CeedQFunctionAddInput(sd->qf, "uR", Nc, CEED_EVAL_NONE));  in DMCeedCreateFVM_Internal()
360   PetscCallCEED(CeedQFunctionAddInput(sd->qf, "geom", Nqdata, CEED_EVAL_NONE));  in DMCeedCreateFVM_Internal()
361   if (createInfo) PetscCallCEED(CeedQFunctionAddInput(sd->qf, "info", Nqinfo, CEED_EVAL_NONE));  in DMCeedCreateFVM_Internal()
362   PetscCallCEED(CeedQFunctionAddOutput(sd->qf, "cL", Nc, CEED_EVAL_NONE));  in DMCeedCreateFVM_Internal()
363   PetscCallCEED(CeedQFunctionAddOutput(sd->qf, "cR", Nc, CEED_EVAL_NONE));  in DMCeedCreateFVM_Internal()
365   PetscCallCEED(CeedQFunctionSetContext(sd->qf, qfCtx));  in DMCeedCreateFVM_Internal()
367 …PetscCallCEED(CeedOperatorCreate(ceed, sd->qf, CEED_QFUNCTION_NONE, CEED_QFUNCTION_NONE, &sd->op));  in DMCeedCreateFVM_Internal()
368   PetscCallCEED(CeedOperatorSetField(sd->op, "uL", sd->erL, CEED_BASIS_NONE, CEED_VECTOR_ACTIVE));  in DMCeedCreateFVM_Internal()
369   PetscCallCEED(CeedOperatorSetField(sd->op, "uR", sd->erR, CEED_BASIS_NONE, CEED_VECTOR_ACTIVE));  in DMCeedCreateFVM_Internal()
370   PetscCallCEED(CeedOperatorSetField(sd->op, "geom", sd->erq, CEED_BASIS_NONE, sd->qd));  in DMCeedCreateFVM_Internal()
371 …if (createInfo) PetscCallCEED(CeedOperatorSetField(sd->op, "info", sd->eri, CEED_BASIS_NONE, sd->q…  in DMCeedCreateFVM_Internal()
372   PetscCallCEED(CeedOperatorSetField(sd->op, "cL", sd->erL, CEED_BASIS_NONE, CEED_VECTOR_ACTIVE));  in DMCeedCreateFVM_Internal()
373   PetscCallCEED(CeedOperatorSetField(sd->op, "cR", sd->erR, CEED_BASIS_NONE, CEED_VECTOR_ACTIVE));  in DMCeedCreateFVM_Internal()
376   sd->func = func;  in DMCeedCreateFVM_Internal()
377   PetscCall(PetscStrallocpy(func_source, &sd->funcSource));  in DMCeedCreateFVM_Internal()
380   *soldata = sd;  in DMCeedCreateFVM_Internal()
426 PetscErrorCode DMCeedComputeGeometry(DM dm, DMCeed sd)  in DMCeedComputeGeometry()  argument
439 …if (sd->geom->op) PetscCallCEED(CeedOperatorApply(sd->geom->op, ccoords, sd->qd, CEED_REQUEST_IMME…  in DMCeedComputeGeometry()
440   else PetscCall(DMPlexCeedComputeGeometryFVM(dm, sd->qd));  in DMCeedComputeGeometry()
447 PetscErrorCode DMCeedComputeInfo(DM dm, DMCeed sd)  in DMCeedComputeInfo()  argument
455   PetscCallCEED(CeedVectorGetArrayWrite(sd->qi, CEED_MEM_HOST, &a));  in DMCeedComputeInfo()
495   PetscCallCEED(CeedVectorRestoreArray(sd->qi, &a));  in DMCeedComputeInfo()