46       PetscInt ki, kj, o = 0;  in LoadData2D()  local
52         for (kj = 0; kj <= user->k[f]; ++kj) {  in LoadData2D()
54 …for (c = 0; c < user->Nc[f]; ++c) closure[o++] = ((kj + joff) * (Ni * user->k[f] + 1) + ki + ioff)…  in LoadData2D()
75         PetscInt ki, kj, kk, o = 0;  in LoadData3D()  local
82             for (kj = 0; kj <= user->k[f]; ++kj) {  in LoadData3D()
84 …< user->Nc[f]; ++c) closure[o++] = (((kk + koff) * (Nj * user->k[f] + 1) + kj + joff) * (Ni * user…  in LoadData3D()
127     PetscInt     closureSize, ki, kj, f, c, foff = 0;  in ReadData2D()  local
133       for (kj = user->k[f]; kj >= 0; --kj) {  in ReadData2D()
138 …PetscCall(PetscPrintf(PETSC_COMM_SELF, "%2.0f", (double)PetscRealPart(closure[(kj * (user->k[f] + …  in ReadData2D()
160     PetscInt     closureSize, ki, kj, kk, f, c, foff = 0;  in ReadData3D()  local
167         for (kj = user->k[f]; kj >= 0; --kj) {  in ReadData3D()
172 …OMM_SELF, "%2.0f", (double)PetscRealPart(closure[((kk * (user->k[f] + 1) + kj) * (user->k[f] + 1) …  in ReadData3D()