26   PetscInt           dim, pdim, spdim, j, s;  in PetscFESetUp_Composite()  local
44   for (s = 0; s < cmp->numSubelements; ++s) {  in PetscFESetUp_Composite()
49     PetscCall(DMPlexGetTransitiveClosure(K, s, PETSC_TRUE, &closureSize, &closure));  in PetscFESetUp_Composite()
56       for (k = 0; k < dof; k++) cmp->embedding[s * spdim + sd++] = off + k;  in PetscFESetUp_Composite()
58     PetscCall(DMPlexRestoreTransitiveClosure(K, s, PETSC_TRUE, &closureSize, &closure));  in PetscFESetUp_Composite()
59 …ement %" PetscInt_FMT " has %" PetscInt_FMT " dual basis vectors != %" PetscInt_FMT, s, sd, spdim);  in PetscFESetUp_Composite()
70   for (s = 0; s < cmp->numSubelements; ++s) {  in PetscFESetUp_Composite()
77       PetscCall(PetscDualSpaceGetFunctional(fem->dualSpace, cmp->embedding[s * spdim + j], &f));  in PetscFESetUp_Composite()
83         invVscalar[(s * spdim + j) * spdim + k] = 0.0;  in PetscFESetUp_Composite()
84 …for (q = 0; q < Nq; ++q) invVscalar[(s * spdim + j) * spdim + k] += Bf[q * spdim + k] * weights[q];  in PetscFESetUp_Composite()
89 …PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&n, &n, &invVscalar[s * spdim * spdim], &n, pivots, &inf…  in PetscFESetUp_Composite()
90 …PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&n, &invVscalar[s * spdim * spdim], &n, pivots, work, &n…  in PetscFESetUp_Composite()
93 …for (s = 0; s < cmp->numSubelements * spdim * spdim; s++) fem->invV[s] = PetscRealPart(invVscalar[…  in PetscFESetUp_Composite()
114   PetscInt           p, s, d, e, j, k;  in PetscFEComputeTabulation_Composite()  local
127     for (s = 0; s < cmp->numSubelements; ++s) {  in PetscFEComputeTabulation_Composite()
133 …= 0; e < dim; ++e) subpoint[d] += cmp->invjac[(s * dim + d) * dim + e] * (points[p * dim + e] - cm…  in PetscFEComputeTabulation_Composite()
137         subpoints[p] = s;  in PetscFEComputeTabulation_Composite()
141 …PetscCheck(s < cmp->numSubelements, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Point %" PetscInt_FMT "…  in PetscFEComputeTabulation_Composite()
154     const PetscInt s = subpoints[p];  in PetscFEComputeTabulation_Composite()  local
159         const PetscInt i = (p * pdim + cmp->embedding[s * spdim + j]) * comp;  in PetscFEComputeTabulation_Composite()
162 …  for (k = 0; k < spdim; ++k) B[i] += fem->invV[(s * spdim + k) * spdim + j] * tmpB[p * spdim + k];  in PetscFEComputeTabulation_Composite()
169           const PetscInt i = ((p * pdim + cmp->embedding[s * spdim + j]) * comp + 0) * dim + d;  in PetscFEComputeTabulation_Composite()
172 …for (k = 0; k < spdim; ++k) D[i] += fem->invV[(s * spdim + k) * spdim + j] * tmpD[(p * spdim + k) …  in PetscFEComputeTabulation_Composite()
180 …        const PetscInt i = ((p * pdim + cmp->embedding[s * spdim + j]) * comp + 0) * dim * dim + d;  in PetscFEComputeTabulation_Composite()
183 …for (k = 0; k < spdim; ++k) H[i] += fem->invV[(s * spdim + k) * spdim + j] * tmpH[(p * spdim + k) …  in PetscFEComputeTabulation_Composite()