[petsc-users] DMPlexComputeCellGeometryFVM: "Cannot handle faces with 1 vertices"

[finnkochinski at keemail.me]

Hello,
with the code below, I create a tetrahedron using DMPlexCreateFromDAG, then I try to run DMPlexComputeCellGeometryFVM on this cell. The latter call fails with output:
...
[0]PETSC ERROR: --------------------- Error Message -----------------
[0]PETSC ERROR: Argument out of range
[0]PETSC ERROR: Cannot handle faces with 1 vertices...
(full output is attached).

What am I doing wrong?
regards
Chris

Here is the code:
static char help[] = "No help \n";

#include <petscdmplex.h>

#undef __FUNCT__
#define __FUNCT__ "main"
int main(int argc,char **args)
{
  PetscErrorCode ierr;
  PetscInitialize(&argc,&args,(char*)0,help);

  DM                dm;
  int cStart,cEnd;
  int fStart,fEnd;

  int depth = 3;
  int dim = 3;

  PetscInt    numPoints[4]        = {1,4,6,4};
  PetscInt    coneSize[15]         = {4,3,3,3,3,2,2,2,2,2,2,0,0,0,0};
  PetscInt    cones[28]            = {1,2,3,4, 5,9,8, 9,6,10, 10,8,7, 5,6,7, 11,12, 12,13, 13,11, 11,14, 12,14, 13,14};
  PetscInt    coneOrientations[28] = {0 };
  PetscScalar vertexCoords[12]     = {0,0,0, 1,0,0, 0,1,0, 0,0,1};

  DMCreate(PETSC_COMM_WORLD, &dm);
  DMSetType(dm, DMPLEX);
  DMSetDimension(dm,dim);
  DMPlexCreateFromDAG(dm, depth, numPoints, coneSize, cones, coneOrientations, vertexCoords);

  DMPlexGetHeightStratum(dm, 0, &cStart, &cEnd);
  DMPlexGetHeightStratum(dm, 1, &fStart, &fEnd);

  for (int k =cStart;k<cEnd;k++){
    double vol;
    double centroid[3];
    double normal[3];
    ierr = DMPlexComputeCellGeometryFVM(dm, k, &vol, centroid,NULL);CHKERRQ(ierr);
    printf("FVM: V=%f c=(%f %f %f) n=(%f %f %f)\n",vol,centroid[0],centroid[1],centroid[2],
      normal[0],normal[1],normal[2]);
  }
  ierr = PetscFinalize();
  return 0;
}


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190312/0cbd41c5/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: log
Type: application/octet-stream
Size: 4241 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190312/0cbd41c5/attachment-0001.obj>