4630 …PetscScalar    *work, *u = NULL, *u_t = NULL, *a = NULL, *elemMat = NULL, *elemMatP = NULL, *elemM…  in DMPlexComputeJacobian_Patch_Internal()  local
4699   elemMatP = hasPrec ? (sz = chunkSize * totDim * totDim, off += sz, work + off - sz) : NULL;  in DMPlexComputeJacobian_Patch_Internal()
4752 …rob, prob, PETSCFE_JACOBIAN_PRE, key, Ncell, cgeomFEM, u, u_t, probAux, a, t, X_tShift, elemMatP));  in DMPlexComputeJacobian_Patch_Internal()
4767             if (hasPrec) elemMatP[eOffset + i * totDim + i] = 1.0;  in DMPlexComputeJacobian_Patch_Internal()
4781 …if (hasPrec) PetscCall(DMPrintCellMatrix(cell, nameP, totDim, totDim, &elemMatP[(c - cStart) * tot…  in DMPlexComputeJacobian_Patch_Internal()
5915   PetscScalar    *u = NULL, *u_t = NULL, *a = NULL, *elemMat = NULL, *elemMatP = NULL;  in DMPlexComputeBdJacobianSingleByLabel()  local
5969 …im : 0), &elemMat, (hasPrec ? (size_t)numFaces * totDim * totDim : 0), &elemMatP, (locA ? (size_t)…  in DMPlexComputeBdJacobianSingleByLabel()
6004     if (elemMatP) PetscCall(PetscArrayzero(elemMatP, numFaces * totDim * totDim));  in DMPlexComputeBdJacobianSingleByLabel()
6028 …cobian(ds, wf, PETSCFE_JACOBIAN_PRE, key, Ne, chunkGeom, u, u_t, dsAux, a, t, X_tShift, elemMatP));  in DMPlexComputeBdJacobianSingleByLabel()
6036 …dsAux, PetscSafePointerPlusOffset(a, offset * totDimAux), t, X_tShift, &elemMatP[offset * totDim *…  in DMPlexComputeBdJacobianSingleByLabel()
6046 …sformPointTensor_Internal(dm, tdm, tv, support[0], PETSC_TRUE, totDim, &elemMatP[face * totDim * t…  in DMPlexComputeBdJacobianSingleByLabel()
6052 … > 1) PetscCall(DMPrintCellMatrix(point, "BdJacobian", totDim, totDim, &elemMatP[face * totDim * t…  in DMPlexComputeBdJacobianSingleByLabel()
6053 …al(plex, section, globalSection, mesh->useMatClPerm, JacP, support[0], &elemMatP[face * totDim * t…  in DMPlexComputeBdJacobianSingleByLabel()
6065     PetscCall(PetscFree5(u, u_t, elemMat, elemMatP, a));  in DMPlexComputeBdJacobianSingleByLabel()
6182   PetscScalar    *elemMat, *elemMatP, *elemMatD, *u, *u_t, *a = NULL;  in DMPlexComputeJacobianByKey()  local
6219 …im : 0), &elemMat, (hasPrec ? (size_t)numCells * totDim * totDim : 0), &elemMatP, (hasDyn ? (size_…  in DMPlexComputeJacobianByKey()
6244   if (hasPrec) PetscCall(PetscArrayzero(elemMatP, numCells * totDim * totDim));  in DMPlexComputeJacobianByKey()
6290 …(prob, prob, PETSCFE_JACOBIAN_PRE, key, Ne, chunkGeom, u, u_t, probAux, a, t, X_tShift, elemMatP));  in DMPlexComputeJacobianByKey()
6291 …obAux, PetscSafePointerPlusOffset(a, offset * totDimAux), t, X_tShift, &elemMatP[offset * totDim *…  in DMPlexComputeJacobianByKey()
6330               elemMatP[eOffset + i * totDim + i] = vol;  in DMPlexComputeJacobianByKey()
6353 …if (mesh->printFEM > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, totDim, &elemMatP[cind * t…  in DMPlexComputeJacobianByKey()
6354 …e_Internal(dm, section, globalSection, mesh->useMatClPerm, JacP, cell, &elemMatP[cind * totDim * t…  in DMPlexComputeJacobianByKey()
6364   PetscCall(PetscFree5(u, u_t, elemMat, elemMatP, elemMatD));  in DMPlexComputeJacobianByKey()
6397   PetscScalar    *elemMat, *elemMatP, *elemMatD, *u, *u_t, *a = NULL;  in DMPlexComputeJacobianByKeyGeneral()  local
6442 …m : 0), &elemMat, (hasPrec ? (size_t)numCells * totDim * ctotDim : 0), &elemMatP, (hasDyn ? (size_…  in DMPlexComputeJacobianByKeyGeneral()
6467   if (hasPrec) PetscCall(PetscArrayzero(elemMatP, numCells * totDim * ctotDim));  in DMPlexComputeJacobianByKeyGeneral()
6513 …bian(rds, cds, PETSCFE_JACOBIAN_PRE, key, Ne, chunkGeom, u, u_t, dsAux, a, t, X_tShift, elemMatP));  in DMPlexComputeJacobianByKeyGeneral()
6514 …dsAux, PetscSafePointerPlusOffset(a, offset * totDimAux), t, X_tShift, &elemMatP[offset * totDim *…  in DMPlexComputeJacobianByKeyGeneral()
6553               elemMatP[eOffset + i * ctotDim + i] = vol;  in DMPlexComputeJacobianByKeyGeneral()
6576 …if (printFEM > 1) PetscCall(DMPrintCellMatrix(cell, name, totDim, ctotDim, &elemMatP[cind * totDim…  in DMPlexComputeJacobianByKeyGeneral()
6577 …balSection, clPerm, dmc, csection, cglobalSection, clPerm, JacP, cell, &elemMatP[cind * totDim * c…  in DMPlexComputeJacobianByKeyGeneral()
6587   PetscCall(PetscFree5(u, u_t, elemMat, elemMatP, elemMatD));  in DMPlexComputeJacobianByKeyGeneral()