51   PetscInt      pdim, j;  in PetscFESetUp_Basic()  local
54   PetscCall(PetscDualSpaceGetDimension(fem->dualSpace, &pdim));  in PetscFESetUp_Basic()
55   PetscCall(PetscMalloc1(pdim * pdim, &fem->invV));  in PetscFESetUp_Basic()
56   for (j = 0; j < pdim; ++j) {  in PetscFESetUp_Basic()
64     PetscCall(PetscMalloc1(Nc * Nq * pdim, &Bf));  in PetscFESetUp_Basic()
66     for (k = 0; k < pdim; ++k) {  in PetscFESetUp_Basic()
68       fem->invV[j * pdim + k] = 0.0;  in PetscFESetUp_Basic()
71 …for (c = 0; c < Nc; ++c) fem->invV[j * pdim + k] += Bf[(q * pdim + k) * Nc + c] * weights[q * Nc +…  in PetscFESetUp_Basic()
77   PetscCall(PetscMalloc2(pdim, &pivots, pdim, &work));  in PetscFESetUp_Basic()
78   PetscCall(PetscBLASIntCast(pdim, &n));  in PetscFESetUp_Basic()
125   PetscInt   pdim; /* Dimension of FE space P */  in PetscFEComputeTabulation_Basic()  local
136   PetscCall(PetscDualSpaceGetDimension(fem->dualSpace, &pdim));  in PetscFEComputeTabulation_Basic()
139   if (K >= 0) PetscCall(DMGetWorkArray(dm, npoints * pdim * Nc, MPIU_REAL, &tmpB));  in PetscFEComputeTabulation_Basic()
140   if (K >= 1) PetscCall(DMGetWorkArray(dm, npoints * pdim * Nc * dim, MPIU_REAL, &tmpD));  in PetscFEComputeTabulation_Basic()
141   if (K >= 2) PetscCall(DMGetWorkArray(dm, npoints * pdim * Nc * dim * dim, MPIU_REAL, &tmpH));  in PetscFEComputeTabulation_Basic()
146     PetscCall(TensorContract_Private(npoints, pdim, Nc, pdim, tmpB, fem->invV, B));  in PetscFEComputeTabulation_Basic()
150     PetscCall(TensorContract_Private(npoints, pdim, Nc * dim, pdim, tmpD, fem->invV, D));  in PetscFEComputeTabulation_Basic()
154     PetscCall(TensorContract_Private(npoints, pdim, Nc * dim * dim, pdim, tmpH, fem->invV, H));  in PetscFEComputeTabulation_Basic()
156   if (K >= 0) PetscCall(DMRestoreWorkArray(dm, npoints * pdim * Nc, MPIU_REAL, &tmpB));  in PetscFEComputeTabulation_Basic()
157   if (K >= 1) PetscCall(DMRestoreWorkArray(dm, npoints * pdim * Nc * dim, MPIU_REAL, &tmpD));  in PetscFEComputeTabulation_Basic()
158   if (K >= 2) PetscCall(DMRestoreWorkArray(dm, npoints * pdim * Nc * dim * dim, MPIU_REAL, &tmpH));  in PetscFEComputeTabulation_Basic()