330     PetscInt sNv, sNc, spdim;  in PetscSpaceEvaluate_Tensor()  local
333     PetscCall(PetscSpaceGetNumVariables(tens->tensspaces[s], &sNv));  in PetscSpaceEvaluate_Tensor()
341       for (PetscInt i = 0; i < sNv; i++) lpoints[p * sNv + i] = points[p * Nv + d + i];  in PetscSpaceEvaluate_Tensor()
375                     if (der >= d && der < d + sNv) {  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()
401                       if (der >= d && der < d + sNv && der2 >= d && der2 < d + sNv) {  in PetscSpaceEvaluate_Tensor()
402 …) * Nv + der) * Nv + der2] *= sH[(((spdim * p + si) * sNc + sc) * sNv + der - d) * sNv + der2 - d];  in PetscSpaceEvaluate_Tensor()
403                       } else if (der >= d && der < d + sNv) {  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()
405                       } else if (der2 >= d && der2 < d + sNv) {  in PetscSpaceEvaluate_Tensor()
406 … * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sD[((spdim * p + si) * sNc + sc) * sNv + der2 - d];  in PetscSpaceEvaluate_Tensor()
419     d += sNv;  in PetscSpaceEvaluate_Tensor()