39   PetscInt     i, j, f, c;  in LoadData2D()  local
49       for (f = 0; f < user->Nf; ++f) {  in LoadData2D()
50         PetscInt ioff = i * user->k[f], joff = j * user->k[f];  in LoadData2D()
52         for (kj = 0; kj <= user->k[f]; ++kj) {  in LoadData2D()
53           for (ki = 0; ki <= user->k[f]; ++ki) {  in LoadData2D()
54 …for (c = 0; c < user->Nc[f]; ++c) closure[o++] = ((kj + joff) * (Ni * user->k[f] + 1) + ki + ioff)…  in LoadData2D()
67   PetscInt     i, j, k, f, c;  in LoadData3D()  local
78         for (f = 0; f < user->Nf; ++f) {  in LoadData3D()
79           PetscInt ioff = i * user->k[f], joff = j * user->k[f], koff = k * user->k[f];  in LoadData3D()
81           for (kk = 0; kk <= user->k[f]; ++kk) {  in LoadData3D()
82             for (kj = 0; kj <= user->k[f]; ++kj) {  in LoadData3D()
83               for (ki = 0; ki <= user->k[f]; ++ki) {  in LoadData3D()
84 … c < user->Nc[f]; ++c) closure[o++] = (((kk + koff) * (Nj * user->k[f] + 1) + kj + joff) * (Ni * u…  in LoadData3D()
127     PetscInt     closureSize, ki, kj, f, c, foff = 0;  in ReadData2D()  local
131     for (f = 0; f < user->Nf; ++f) {  in ReadData2D()
132       PetscCall(PetscPrintf(PETSC_COMM_SELF, "  Field %" PetscInt_FMT "\n", f));  in ReadData2D()
133       for (kj = user->k[f]; kj >= 0; --kj) {  in ReadData2D()
134         for (ki = 0; ki <= user->k[f]; ++ki) {  in ReadData2D()
136           for (c = 0; c < user->Nc[f]; ++c) {  in ReadData2D()
138 …MM_SELF, "%2.0f", (double)PetscRealPart(closure[(kj * (user->k[f] + 1) + ki) * user->Nc[f] + c + f…  in ReadData2D()
144       foff += PetscSqr(user->k[f] + 1);  in ReadData2D()
160     PetscInt     closureSize, ki, kj, kk, f, c, foff = 0;  in ReadData3D()  local
164     for (f = 0; f < user->Nf; ++f) {  in ReadData3D()
165       PetscCall(PetscPrintf(PETSC_COMM_SELF, "  Field %" PetscInt_FMT "\n", f));  in ReadData3D()
166       for (kk = user->k[f]; kk >= 0; --kk) {  in ReadData3D()
167         for (kj = user->k[f]; kj >= 0; --kj) {  in ReadData3D()
168           for (ki = 0; ki <= user->k[f]; ++ki) {  in ReadData3D()
170             for (c = 0; c < user->Nc[f]; ++c) {  in ReadData3D()
172 …0f", (double)PetscRealPart(closure[((kk * (user->k[f] + 1) + kj) * (user->k[f] + 1) + ki) * user->…  in ReadData3D()
180       foff += PetscSqr(user->k[f] + 1);  in ReadData3D()
190   PetscInt dim, f, o, i, j, k, c, d;  in SetSymmetries()  local
196   for (f = 0; f < user->Nf; f++) {  in SetSymmetries()
199 …if (user->k[f] < 3) continue; /* No symmetries needed for order < 3, because no cell, facet, edge …  in SetSymmetries()
204         PetscInt   numDof  = user->k[f] - 1;  in SetSymmetries()
205         PetscInt   numComp = user->Nc[f];  in SetSymmetries()
226         PetscInt   perEdge = user->k[f] - 1;  in SetSymmetries()
228         PetscInt   numComp = user->Nc[f];  in SetSymmetries()
300     PetscCall(PetscSectionSetFieldSym(s, f, sym));  in SetSymmetries()
313   PetscInt     dim, size = 0, f;  in main()  local
328     for (f = 0; f < user.Nf; ++f) {  in main()
329 … for (d = 0; d <= dim; ++d) numDof[f * (dim + 1) + d] = PetscPowInt(user.k[f] - 1, d) * user.Nc[f];  in main()
330       size += PetscPowInt(user.k[f] + 1, d) * user.Nc[f];  in main()