351               for (PetscInt sc = 0; sc < sNc; sc++) {  in PetscSpaceEvaluate_Tensor()  local
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()  local
372                   PetscInt c = (l * sNc + sc) * cstep + m;  in PetscSpaceEvaluate_Tensor()
376 … p++) D[((pdim * p + i) * Nc + c) * Nv + der] *= sD[((spdim * p + si) * sNc + sc) * sNv + der - d];  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()  local
397                   PetscInt c = (l * sNc + sc) * cstep + m;  in PetscSpaceEvaluate_Tensor()
402 …p + i) * Nc + c) * Nv + der) * Nv + der2] *= sH[(((spdim * p + si) * sNc + sc) * sNv + der - d) * …  in PetscSpaceEvaluate_Tensor()
404 …m * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sD[((spdim * p + si) * sNc + sc) * sNv + der - d];  in PetscSpaceEvaluate_Tensor()
406 … * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sD[((spdim * p + si) * sNc + sc) * sNv + der2 - d];  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()