3792     const PetscInt cind = c - cStart;  in DMPlexGetCellFields()  local
3797     for (i = 0; i < totDim; ++i) ul[cind * totDim + i] = x[i];  in DMPlexGetCellFields()
3801       for (i = 0; i < totDim; ++i) ul_t[cind * totDim + i] = x_t[i];  in DMPlexGetCellFields()
3808       for (i = 0; i < totDimAux; ++i) al[cind * totDimAux + i] = x[i];  in DMPlexGetCellFields()
3894     const PetscInt  cind = c - cStart;  in DMPlexGetHybridCellFields()  local
3896 …PetscScalar    *ul = &(*u)[cind * totDim], *ul_t = PetscSafePointerPlusOffset(*u_t, cind * totDim);  in DMPlexGetHybridCellFields()
3959       PetscScalar *al = &(*a)[cind * totDimAux];  in DMPlexGetHybridCellFields()
4025     const PetscInt  cind = c - cStart;  in DMPlexGetHybridFields()  local
4033       PetscScalar   *closure = NULL, *xl = &x[s][cind * tdX];  in DMPlexGetHybridFields()
4541         const PetscInt cind = c - cStart;  in DMPlexComputeResidual_Patch_Internal()  local
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()
5352         const PetscInt cind = c - cStart;  in DMPlexComputeResidualByKey()  local
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()
5799       const PetscInt cind = c - cStart;  in DMPlexComputeResidualHybridByKey()  local
5804         PetscInt  Nb, off = cind * totDim, soff = cind * totDimScale[0];  in DMPlexComputeResidualHybridByKey()
5823 …for (i = cind * totDim; i < (cind + 1) * totDim; ++i) elemVecCoh[i] += elemVecNeg[i] + elemVecPos[…  in DMPlexComputeResidualHybridByKey()
5825 …if (mesh->printFEM > 1) PetscCall(DMPrintCellVector(cell, name, totDim, &elemVecCoh[cind * totDim]…  in DMPlexComputeResidualHybridByKey()
5832 …PetscCall(DMPlexVecSetClosure(dm, section, locF, cell, &elemVecCoh[cind * totDim], ADD_ALL_VALUES)…  in DMPlexComputeResidualHybridByKey()
6223     const PetscInt cind = c - cStart;  in DMPlexComputeJacobianByKey()  local
6228     for (i = 0; i < totDim; ++i) u[cind * totDim + i] = x[i];  in DMPlexComputeJacobianByKey()
6232       for (i = 0; i < totDim; ++i) u_t[cind * totDim + i] = x_t[i];  in DMPlexComputeJacobianByKey()
6239       for (i = 0; i < totDimAux; ++i) a[cind * totDimAux + i] = x[i];  in DMPlexComputeJacobianByKey()
6320         const PetscInt cind    = c - cStart;  in DMPlexComputeJacobianByKey()  local
6321         const PetscInt eOffset = cind * totDim * totDim;  in DMPlexComputeJacobianByKey()
6344     const PetscInt cind = c - cStart;  in DMPlexComputeJacobianByKey()  local
6347 …sformPointTensor_Internal(dm, tdm, tv, cell, PETSC_TRUE, totDim, &elemMat[cind * totDim * totDim])…  in DMPlexComputeJacobianByKey()
6350 …FEM > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, totDim, &elemMat[cind * totDim * totDim])…  in DMPlexComputeJacobianByKey()
6351 …ernal(dm, section, globalSection, mesh->useMatClPerm, Jac, cell, &elemMat[cind * totDim * totDim],…  in DMPlexComputeJacobianByKey()
6353 …EM > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, totDim, &elemMatP[cind * totDim * totDim])…  in DMPlexComputeJacobianByKey()
6354 …nal(dm, section, globalSection, mesh->useMatClPerm, JacP, cell, &elemMatP[cind * totDim * totDim],…  in DMPlexComputeJacobianByKey()
6357 …FEM > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, totDim, &elemMat[cind * totDim * totDim])…  in DMPlexComputeJacobianByKey()
6358 …rnal(dm, section, globalSection, mesh->useMatClPerm, JacP, cell, &elemMat[cind * totDim * totDim],…  in DMPlexComputeJacobianByKey()
6446     const PetscInt cind = c - cStart;  in DMPlexComputeJacobianByKeyGeneral()  local
6451     for (i = 0; i < totDim; ++i) u[cind * totDim + i] = x[i];  in DMPlexComputeJacobianByKeyGeneral()
6455       for (i = 0; i < totDim; ++i) u_t[cind * totDim + i] = x_t[i];  in DMPlexComputeJacobianByKeyGeneral()
6462       for (i = 0; i < totDimAux; ++i) a[cind * totDimAux + i] = x[i];  in DMPlexComputeJacobianByKeyGeneral()
6543         const PetscInt cind    = c - cStart;  in DMPlexComputeJacobianByKeyGeneral()  local
6544         const PetscInt eOffset = cind * totDim * ctotDim;  in DMPlexComputeJacobianByKeyGeneral()
6567     const PetscInt cind = c - cStart;  in DMPlexComputeJacobianByKeyGeneral()  local
6570 …formPointTensor_Internal(dmr, tdm, tv, cell, PETSC_TRUE, totDim, &elemMat[cind * totDim * ctotDim]…  in DMPlexComputeJacobianByKeyGeneral()
6573 …if (printFEM > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, ctotDim, &elemMat[cind * totDim …  in DMPlexComputeJacobianByKeyGeneral()
6574 …ction, clPerm, dmc, csection, cglobalSection, clPerm, Jac, cell, &elemMat[cind * totDim * ctotDim]…  in DMPlexComputeJacobianByKeyGeneral()
6576 …if (printFEM > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, ctotDim, &elemMatP[cind * totDim…  in DMPlexComputeJacobianByKeyGeneral()
6577 …ion, clPerm, dmc, csection, cglobalSection, clPerm, JacP, cell, &elemMatP[cind * totDim * ctotDim]…  in DMPlexComputeJacobianByKeyGeneral()
6580 …if (printFEM > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, ctotDim, &elemMat[cind * totDim …  in DMPlexComputeJacobianByKeyGeneral()
6581 …tion, clPerm, dmc, csection, cglobalSection, clPerm, JacP, cell, &elemMat[cind * totDim * ctotDim]…  in DMPlexComputeJacobianByKeyGeneral()
6898       const PetscInt cind = c - cS, coff = cind * totDim * totDim;  in DMPlexComputeJacobianHybridByKey()  local
6903         PetscInt  Nb, soff = cind * totDimScale[0], off = 0;  in DMPlexComputeJacobianHybridByKey()
6936 … > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, totDim, &elemMatCoh[cind * totDim * totDim])…  in DMPlexComputeJacobianHybridByKey()
6937 …(plex, section, globalSection, mesh->useMatClPerm, Jac, cell, &elemMatCoh[cind * totDim * totDim],…  in DMPlexComputeJacobianHybridByKey()
6939 …> 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, totDim, &elemMatCohP[cind * totDim * totDim])…  in DMPlexComputeJacobianHybridByKey()
6940 …PetscCall(DMPlexMatSetClosure(plex, section, globalSection, JacP, cell, &elemMatCohP[cind * totDim…  in DMPlexComputeJacobianHybridByKey()
6942 … > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, totDim, &elemMatCoh[cind * totDim * totDim])…  in DMPlexComputeJacobianHybridByKey()
6943 …plex, section, globalSection, mesh->useMatClPerm, JacP, cell, &elemMatCoh[cind * totDim * totDim],…  in DMPlexComputeJacobianHybridByKey()
7046     const PetscInt cind = c - cStart;  in DMPlexComputeJacobianActionByKey()  local
7051     for (i = 0; i < totDim; ++i) u[cind * totDim + i] = x[i];  in DMPlexComputeJacobianActionByKey()
7055       for (i = 0; i < totDim; ++i) u_t[cind * totDim + i] = x_t[i];  in DMPlexComputeJacobianActionByKey()
7062       for (i = 0; i < totDimAux; ++i) a[cind * totDimAux + i] = x[i];  in DMPlexComputeJacobianActionByKey()
7066     for (i = 0; i < totDim; ++i) y[cind * totDim + i] = x[i];  in DMPlexComputeJacobianActionByKey()
7122     const PetscInt     cind = c - cStart;  in DMPlexComputeJacobianActionByKey()  local
7128 …cCallBLAS("BLASgemv", BLASgemv_("N", &M, &M, &a, &elemMat[cind * totDim * totDim], &M, &y[cind * t…  in DMPlexComputeJacobianActionByKey()
7130       PetscCall(DMPrintCellMatrix(c, name, totDim, totDim, &elemMat[cind * totDim * totDim]));  in DMPlexComputeJacobianActionByKey()
7131       PetscCall(DMPrintCellVector(c, "Y", totDim, &y[cind * totDim]));  in DMPlexComputeJacobianActionByKey()