Lines Matching refs:npoints
298 static PetscErrorCode PetscSpaceEvaluate_Tensor(PetscSpace sp, PetscInt npoints, const PetscReal po… in PetscSpaceEvaluate_Tensor() argument
311 PetscCall(PetscSpaceEvaluate(sp, npoints, points, B, D, H)); in PetscSpaceEvaluate_Tensor()
316 PetscCall(DMGetWorkArray(dm, npoints * Nv, MPIU_REAL, &lpoints)); 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()
340 for (PetscInt p = 0; p < npoints; ++p) { in PetscSpaceEvaluate_Tensor()
343 PetscCall(PetscSpaceEvaluate(tens->tensspaces[s], npoints, lpoints, sB, sD, sH)); 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()
426 PetscCall(DMRestoreWorkArray(dm, npoints * Nv, MPIU_REAL, &lpoints)); in PetscSpaceEvaluate_Tensor()