4477     PetscScalar     *elemVec, *fluxL = NULL, *fluxR = NULL;  in DMPlexComputeResidual_Patch_Internal()  local
4487       PetscCall(DMGetWorkArray(dm, numCells * totDim, MPIU_SCALAR, &elemVec));  in DMPlexComputeResidual_Patch_Internal()
4488       PetscCall(PetscArrayzero(elemVec, numCells * totDim));  in DMPlexComputeResidual_Patch_Internal()
4523 …     PetscCall(PetscFEIntegrateResidual(prob, key, Ne, chunkGeom, u, u_t, probAux, a, t, elemVec));  in DMPlexComputeResidual_Patch_Internal()
4525 …erPlusOffset(u_t, offset * totDim), probAux, &a[offset * totDimAux], t, &elemVec[offset * totDim])…  in DMPlexComputeResidual_Patch_Internal()
4543 … if (mesh->printFEM > 1) PetscCall(DMPrintCellVector(cell, name, totDim, &elemVec[cind * totDim]));  in DMPlexComputeResidual_Patch_Internal()
4550 …  PetscCall(DMPlexVecSetClosure(dm, section, locF, cell, &elemVec[cind * totDim], ADD_ALL_VALUES));  in DMPlexComputeResidual_Patch_Internal()
4590       PetscCall(DMRestoreWorkArray(dm, numCells * totDim, MPIU_SCALAR, &elemVec));  in DMPlexComputeResidual_Patch_Internal()
4927   PetscScalar    *u = NULL, *u_t = NULL, *a = NULL, *elemVec = NULL;  in DMPlexComputeBdResidualSingleByKey()  local
4966 … &u, (locX_t ? (size_t)numFaces * totDim : 0), &u_t, numFaces * totDim, &elemVec, (locA ? (size_t)…  in DMPlexComputeBdResidualSingleByKey()
5001     PetscCall(PetscArrayzero(elemVec, numFaces * totDim));  in DMPlexComputeBdResidualSingleByKey()
5023 …etscCall(PetscFEIntegrateBdResidual(prob, wf, key, Ne, chunkGeom, u, u_t, probAux, a, t, elemVec));  in DMPlexComputeBdResidualSingleByKey()
5026 …totDim), probAux, PetscSafePointerPlusOffset(a, offset * totDimAux), t, &elemVec[offset * totDim])…  in DMPlexComputeBdResidualSingleByKey()
5032 …if (mesh->printFEM > 1) PetscCall(DMPrintCellVector(point, name, totDim, &elemVec[face * totDim]));  in DMPlexComputeBdResidualSingleByKey()
5034 …PetscCall(DMPlexVecSetClosure(plex, NULL, locF, support[0], &elemVec[face * totDim], ADD_ALL_VALUE…  in DMPlexComputeBdResidualSingleByKey()
5040     PetscCall(PetscFree4(u, u_t, elemVec, a));  in DMPlexComputeBdResidualSingleByKey()
5280     PetscScalar     *elemVec, *fluxL, *fluxR;  in DMPlexComputeResidualByKey()  local
5290       PetscCall(DMGetWorkArray(dm, numCells * totDim, MPIU_SCALAR, &elemVec));  in DMPlexComputeResidualByKey()
5291       PetscCall(PetscArrayzero(elemVec, numCells * totDim));  in DMPlexComputeResidualByKey()
5334         PetscCall(PetscFEIntegrateResidual(ds, key, Ne, chunkGeom, u, u_t, dsAux, a, t, elemVec));  in DMPlexComputeResidualByKey()
5336 …* totDim), dsAux, PetscSafePointerPlusOffset(a, offset * totDimAux), t, &elemVec[offset * totDim])…  in DMPlexComputeResidualByKey()
5354 … if (mesh->printFEM > 1) PetscCall(DMPrintCellVector(cell, name, totDim, &elemVec[cind * totDim]));  in DMPlexComputeResidualByKey()
5361 …  PetscCall(DMPlexVecSetClosure(dm, section, locF, cell, &elemVec[cind * totDim], ADD_ALL_VALUES));  in DMPlexComputeResidualByKey()
5450       PetscCall(DMRestoreWorkArray(dm, numCells * totDim, MPIU_SCALAR, &elemVec));  in DMPlexComputeResidualByKey()