Searched refs:valsum (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/dm/impls/plex/ |
| H A D | plexfem.c | 2151 PetscScalar *interpolant, *valsum; in DMPlexComputeClementInterpolant() local
2188 …PetscCall(PetscMalloc6(Nc * 2, &valsum, Nc, &interpolant, cdim * Nq, &coords, Nq, &fegeom.detJ, cd… in DMPlexComputeClementInterpolant()
2196 PetscCall(PetscArrayzero(valsum, Nc)); in DMPlexComputeClementInterpolant()
2200 PetscScalar *val = &valsum[Nc]; in DMPlexComputeClementInterpolant()
2240 for (fc = 0; fc < Nc; ++fc) valsum[fc] += val[fc]; in DMPlexComputeClementInterpolant()
2251 for (fc = 0; fc < Nc; ++fc) valsum[fc] /= volsum; in DMPlexComputeClementInterpolant()
2253 PetscCall(DMPlexVecSetClosure(dmc, NULL, locC, v, valsum, INSERT_VALUES)); in DMPlexComputeClementInterpolant()
2255 PetscCall(PetscFree6(valsum, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMPlexComputeClementInterpolant()
|