Lines Matching refs:pdim
306 PetscInt pdim; in PetscSpaceEvaluate_Tensor() local
315 PetscCall(PetscSpaceGetDimension(sp, &pdim)); in PetscSpaceEvaluate_Tensor()
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()
336 …pdim % vstep) && !(pdim % spdim), PETSC_COMM_SELF, PETSC_ERR_PLIB, "Bad tensor loop: Nv %" PetscIn… in PetscSpaceEvaluate_Tensor()
338 vskip = pdim / (vstep * spdim); 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()