330     PetscInt sNv, sNc, spdim;  in PetscSpaceEvaluate_Tensor()  local
335     PetscCall(PetscSpaceGetDimension(tens->tensspaces[s], &spdim));  in PetscSpaceEvaluate_Tensor()
336 …spdim), PETSC_COMM_SELF, PETSC_ERR_PLIB, "Bad tensor loop: Nv %" PetscInt_FMT ", Ns %" PetscInt_FM…  in PetscSpaceEvaluate_Tensor()
337 …", s %" PetscInt_FMT ", cstep %" PetscInt_FMT ", sNc %" PetscInt_FMT, Nv, Ns, Nc, s, cstep, spdim);  in PetscSpaceEvaluate_Tensor()
338     vskip = pdim / (vstep * spdim);  in PetscSpaceEvaluate_Tensor()
346         for (PetscInt si = 0; si < spdim; si++) {  in PetscSpaceEvaluate_Tensor()
348             PetscInt i = (k * spdim + si) * vstep + j;  in PetscSpaceEvaluate_Tensor()
355 …for (PetscInt p = 0; p < npoints; p++) B[(pdim * p + i) * Nc + c] *= sB[(spdim * p + si) * sNc + s…  in PetscSpaceEvaluate_Tensor()
365         for (PetscInt si = 0; si < spdim; si++) {  in PetscSpaceEvaluate_Tensor()
367             PetscInt i = (k * spdim + si) * vstep + j;  in PetscSpaceEvaluate_Tensor()
376 … p = 0; p < npoints; p++) D[((pdim * p + i) * Nc + c) * Nv + der] *= sD[((spdim * p + si) * sNc + …  in PetscSpaceEvaluate_Tensor()
378 …t p = 0; p < npoints; p++) D[((pdim * p + i) * Nc + c) * Nv + der] *= sB[(spdim * p + si) * sNc + …  in PetscSpaceEvaluate_Tensor()
390         for (PetscInt si = 0; si < spdim; si++) {  in PetscSpaceEvaluate_Tensor()
392             PetscInt i = (k * spdim + si) * vstep + j;  in PetscSpaceEvaluate_Tensor()
402 …ints; p++) H[(((pdim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sH[(((spdim * p + si) * sNc + …  in PetscSpaceEvaluate_Tensor()
404 …oints; p++) H[(((pdim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sD[((spdim * p + si) * sNc + …  in PetscSpaceEvaluate_Tensor()
406 …oints; p++) H[(((pdim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sD[((spdim * p + si) * sNc + …  in PetscSpaceEvaluate_Tensor()
408 …points; p++) H[(((pdim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sB[(spdim * p + si) * sNc + …  in PetscSpaceEvaluate_Tensor()
420     vstep *= spdim;  in PetscSpaceEvaluate_Tensor()