Searched refs:cintegral (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/snes/utils/ |
| H A D | dmplexsnes.c | 248 PetscScalar *cintegral; in DMPlexSNESComputeObjectiveFEM() local
254 PetscCall(PetscCalloc1((cEnd - cStart) * Nf, &cintegral)); in DMPlexSNESComputeObjectiveFEM()
256 PetscCall(DMPlexComputeIntegral_Internal(dm, X, cStart, cEnd, cintegral, ctx)); in DMPlexSNESComputeObjectiveFEM()
260 for (PetscInt f = 0; f < Nf; ++f) *obj += PetscRealPart(cintegral[(c - cStart) * Nf + f]); in DMPlexSNESComputeObjectiveFEM()
262 PetscCall(PetscFree(cintegral)); in DMPlexSNESComputeObjectiveFEM()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexfem.c | 2401 …ral_Internal(DM dm, Vec locX, PetscInt cStart, PetscInt cEnd, PetscScalar *cintegral, PetscCtx ctx) in DMPlexComputeIntegral_Internal() argument
2543 PetscCall(PetscFEIntegrate(prob, f, Ne, chunkGeom, u, probAux, a, cintegral)); in DMPlexComputeIntegral_Internal()
2545 …t * totDim], probAux, PetscSafePointerPlusOffset(a, offset * totDimAux), &cintegral[offset * Nf])); in DMPlexComputeIntegral_Internal()
2561 cintegral[c * Nf + f] += PetscRealPart(lint) * cgeomFVM[c].volume; in DMPlexComputeIntegral_Internal()
2603 PetscScalar *cintegral, *lintegral; in DMPlexComputeIntegralFEM() local
2617 PetscCall(PetscCalloc2(Nf, &lintegral, (cEnd - cStart) * Nf, &cintegral)); in DMPlexComputeIntegralFEM()
2623 PetscCall(DMPlexComputeIntegral_Internal(dm, locX, cStart, cEnd, cintegral, ctx)); in DMPlexComputeIntegralFEM()
2630 if (printFEM > 1) PetscCall(DMPrintCellVector(cell, "Cell Integral", Nf, &cintegral[c * Nf])); in DMPlexComputeIntegralFEM()
2631 for (f = 0; f < Nf; ++f) lintegral[f] += cintegral[c * Nf + f]; in DMPlexComputeIntegralFEM()
2639 PetscCall(PetscFree2(lintegral, cintegral)); in DMPlexComputeIntegralFEM()
[all …]
|