330     PetscInt sNv, sNc, spdim;  in PetscSpaceEvaluate_Tensor()  local
334     PetscCall(PetscSpaceGetNumComponents(tens->tensspaces[s], &sNc));  in PetscSpaceEvaluate_Tensor()
337 …PetscCheck(!(Nc % cstep) && !(Nc % sNc), PETSC_COMM_SELF, PETSC_ERR_PLIB, "Bad tensor loop: Nv %" …  in PetscSpaceEvaluate_Tensor()
339     cskip = Nc / (cstep * sNc);  in PetscSpaceEvaluate_Tensor()
351               for (PetscInt sc = 0; sc < sNc; sc++) {  in PetscSpaceEvaluate_Tensor()
353                   PetscInt c = (l * sNc + sc) * cstep + m;  in PetscSpaceEvaluate_Tensor()
355 …r (PetscInt p = 0; p < npoints; p++) B[(pdim * p + i) * Nc + c] *= sB[(spdim * p + si) * sNc + sc];  in PetscSpaceEvaluate_Tensor()
370               for (PetscInt sc = 0; sc < sNc; sc++) {  in PetscSpaceEvaluate_Tensor()
372                   PetscInt c = (l * sNc + sc) * cstep + m;  in PetscSpaceEvaluate_Tensor()
376 …ts; p++) D[((pdim * p + i) * Nc + c) * Nv + der] *= sD[((spdim * p + si) * sNc + sc) * sNv + der -…  in PetscSpaceEvaluate_Tensor()
378 … = 0; p < npoints; p++) D[((pdim * p + i) * Nc + c) * Nv + der] *= sB[(spdim * p + si) * sNc + sc];  in PetscSpaceEvaluate_Tensor()
395               for (PetscInt sc = 0; sc < sNc; sc++) {  in PetscSpaceEvaluate_Tensor()
397                   PetscInt c = (l * sNc + sc) * cstep + m;  in PetscSpaceEvaluate_Tensor()
402 …dim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sH[(((spdim * p + si) * sNc + sc) * sNv + der -…  in PetscSpaceEvaluate_Tensor()
404 …pdim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sD[((spdim * p + si) * sNc + sc) * sNv + der -…  in PetscSpaceEvaluate_Tensor()
406 …pdim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sD[((spdim * p + si) * sNc + sc) * sNv + der2 …  in PetscSpaceEvaluate_Tensor()
408 …nts; p++) H[(((pdim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sB[(spdim * p + si) * sNc + sc];  in PetscSpaceEvaluate_Tensor()
421     cstep *= sNc;  in PetscSpaceEvaluate_Tensor()