Lines Matching refs:cell
146 static PetscErrorCode CheckFEMGeometry(DM dm, PetscInt cell, PetscInt spaceDim, PetscReal v0Ex[], P… in CheckFEMGeometry() argument
152 PetscCall(DMPlexComputeCellGeometryFEM(dm, cell, NULL, v0, J, invJ, &detJ)); in CheckFEMGeometry()
175 static PetscErrorCode CheckFVMGeometry(DM dm, PetscInt cell, PetscInt spaceDim, PetscReal centroidE… in CheckFVMGeometry() argument
182 …PetscCall(DMPlexComputeCellGeometryFVM(dm, cell, volEx ? &vol : NULL, centroidEx ? centroid : NULL… in CheckFVMGeometry()
185 …" PetscInt_FMT ", Invalid centroid[%" PetscInt_FMT "]: %g != %g diff %g", cell, d, (double)centroi… in CheckFVMGeometry()
186 …B, "Cell %" PetscInt_FMT ", Invalid normal[%" PetscInt_FMT "]: %g != %g", cell, d, (double)normal[… in CheckFVMGeometry()
188 …SC_ERR_PLIB, "Cell %" PetscInt_FMT ", Invalid volume = %g != %g diff %g", cell, (double)vol, (doub… in CheckFVMGeometry()
192 static PetscErrorCode CheckGaussLaw(DM dm, PetscInt cell) in CheckGaussLaw() argument
204 PetscCall(DMPlexGetCellType(dm, cell, &ct)); in CheckGaussLaw()
206 PetscCall(DMPlexGetConeSize(dm, cell, &coneSize)); in CheckGaussLaw()
207 PetscCall(DMPlexGetCone(dm, cell, &cone)); in CheckGaussLaw()
208 PetscCall(DMPlexGetConeOrientation(dm, cell, &ornt)); in CheckGaussLaw()
216 …ponent %" PetscInt_FMT ": %g != 0. as it should be for a constant field", cell, d, (double)integra… in CheckGaussLaw()
220 static PetscErrorCode CheckCell(DM dm, PetscInt cell, PetscBool transform, PetscReal v0Ex[], PetscR… in CheckCell() argument
230 if (v0Ex) PetscCall(CheckFEMGeometry(dm, cell, cdim, v0Ex, JEx, invJEx, detJEx)); in CheckCell()
232 PetscCall(CheckFVMGeometry(dm, cell, cdim, centroidEx, normalEx, volEx)); in CheckCell()
233 PetscCall(CheckGaussLaw(dm, cell)); in CheckCell()
235 PetscCall(DMPlexGetConeSize(dm, cell, &coneSize)); in CheckCell()
236 PetscCall(DMPlexGetCone(dm, cell, &cone)); in CheckCell()
251 PetscCall(DMPlexVecGetClosure(dm, coordSection, coordinates, cell, &coordSize, &coords)); in CheckCell()
256 PetscCall(DMPlexVecRestoreClosure(dm, coordSection, coordinates, cell, &coordSize, &coords)); in CheckCell()
388 PetscCall(CheckFVMGeometry(dm, cell, cdim, centroidExT, normalExT, volExT)); in CheckCell()
389 PetscCall(CheckGaussLaw(dm, cell)); in CheckCell()
391 PetscCall(DMPlexGetConeSize(dm, cell, &coneSize)); in CheckCell()
392 PetscCall(DMPlexGetCone(dm, cell, &cone)); in CheckCell()