31   PetscInt           Ns, Nc, i, Nv, deg;  in PetscSpaceSetFromOptions_Tensor()  local
37   PetscCall(PetscSpaceGetNumComponents(sp, &Nc));  in PetscSpaceSetFromOptions_Tensor()
68     Ncs = (PetscInt)PetscPowReal((PetscReal)Nc, 1. / Ns);  in PetscSpaceSetFromOptions_Tensor()
69 …Nc % PetscPowInt(Ncs, Ns) == 0, PetscObjectComm((PetscObject)sp), PETSC_ERR_ARG_WRONG, "Cannot use…  in PetscSpaceSetFromOptions_Tensor()
133   PetscInt           Nc, Nv, Ns;  in PetscSpaceSetUp_Tensor()  local
141   PetscCall(PetscSpaceGetNumComponents(sp, &Nc));  in PetscSpaceSetUp_Tensor()
168       Ncs = (PetscInt)(PetscPowReal((PetscReal)Nc, 1. / Ns));  in PetscSpaceSetUp_Tensor()
169 …Nc % PetscPowInt(Ncs, Ns) == 0, PetscObjectComm((PetscObject)sp), PETSC_ERR_ARG_WRONG, "Cannot use…  in PetscSpaceSetUp_Tensor()
198 …Nc % Ncprod == 0, PetscObjectComm((PetscObject)sp), PETSC_ERR_ARG_WRONG, "Product of subspace comp…  in PetscSpaceSetUp_Tensor()
200     if (Ncprod != Nc) {  in PetscSpaceSetUp_Tensor()
201       PetscInt    Ncopies = Nc / Ncprod;  in PetscSpaceSetUp_Tensor()
302   PetscInt           Nc   = sp->Nc;  in PetscSpaceEvaluate_Tensor()  local
317   if (B || D || H) PetscCall(DMGetWorkArray(dm, npoints * pdim * Nc, MPIU_REAL, &sB));  in PetscSpaceEvaluate_Tensor()
318   if (D || H) PetscCall(DMGetWorkArray(dm, npoints * pdim * Nc * Nv, MPIU_REAL, &sD));  in PetscSpaceEvaluate_Tensor()
319   if (H) PetscCall(DMGetWorkArray(dm, npoints * pdim * Nc * Nv * Nv, MPIU_REAL, &sH));  in PetscSpaceEvaluate_Tensor()
321     for (PetscInt i = 0; i < npoints * pdim * Nc; i++) B[i] = 1.;  in PetscSpaceEvaluate_Tensor()
324     for (PetscInt i = 0; i < npoints * pdim * Nc * Nv; i++) D[i] = 1.;  in PetscSpaceEvaluate_Tensor()
327     for (PetscInt i = 0; i < npoints * pdim * Nc * Nv * Nv; i++) H[i] = 1.;  in PetscSpaceEvaluate_Tensor()
337 …Nc % cstep) && !(Nc % sNc), PETSC_COMM_SELF, PETSC_ERR_PLIB, "Bad tensor loop: Nv %" PetscInt_FMT …  in PetscSpaceEvaluate_Tensor()
339     cskip = Nc / (cstep * sNc);  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()
376 …for (PetscInt p = 0; p < npoints; p++) D[((pdim * p + i) * Nc + c) * Nv + der] *= sD[((spdim * p +…  in PetscSpaceEvaluate_Tensor()
378 …for (PetscInt p = 0; p < npoints; p++) D[((pdim * p + i) * Nc + c) * Nv + der] *= sB[(spdim * p + …  in PetscSpaceEvaluate_Tensor()
402 …for (PetscInt p = 0; p < npoints; p++) H[(((pdim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sH…  in PetscSpaceEvaluate_Tensor()
404 …for (PetscInt p = 0; p < npoints; p++) H[(((pdim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sD…  in PetscSpaceEvaluate_Tensor()
406 …for (PetscInt p = 0; p < npoints; p++) H[(((pdim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sD…  in PetscSpaceEvaluate_Tensor()
408 …for (PetscInt p = 0; p < npoints; p++) H[(((pdim * p + i) * Nc + c) * Nv + der) * Nv + der2] *= sB…  in PetscSpaceEvaluate_Tensor()
423   if (H) PetscCall(DMRestoreWorkArray(dm, npoints * pdim * Nc * Nv * Nv, MPIU_REAL, &sH));  in PetscSpaceEvaluate_Tensor()
424   if (D || H) PetscCall(DMRestoreWorkArray(dm, npoints * pdim * Nc * Nv, MPIU_REAL, &sD));  in PetscSpaceEvaluate_Tensor()
425   if (B || D || H) PetscCall(DMRestoreWorkArray(dm, npoints * pdim * Nc, MPIU_REAL, &sB));  in PetscSpaceEvaluate_Tensor()
575   PetscInt           Nc, dim, order, i;  in PetscSpaceGetHeightSubspace_Tensor()  local
581   PetscCall(PetscSpaceGetNumComponents(sp, &Nc));  in PetscSpaceGetHeightSubspace_Tensor()
595       PetscCall(PetscSpaceSetNumComponents(sub, Nc));  in PetscSpaceGetHeightSubspace_Tensor()