[petsc-dev] Issues with DMComplexSetCone

[Chris Eldred]

The following code snippet does not appear to be working:

PetscInt, dimension(4) :: EC

EC(1) = 1
EC(2) = 2
EC(3) = 3
EC(4) = 4
write(*,*) 'cell',EC
call DMComplexSetCone(model_mesh, lcell , EC, ierr)
CHKERRQ(ierr)
call DMComplexGetCone(model_mesh, lcell , EC, ierr)
CHKERRQ(ierr)
write(*,*) 'cell', EC

I get the following as output:

cell           1           2           3        4
cell    18673792           0          -1          -1

DMComplexGetConeSize returns 4 as expected. I am using the latest
version of petsc-dev.

Later on there is a segmentation fault violation, probably related to
this issue.

Any ideas?

-- 
Chris Eldred
DOE Computational Science Graduate Fellow
Graduate Student, Atmospheric Science, Colorado State University
B.S. Applied Computational Physics, Carnegie Mellon University, 2009
chris.eldred at gmail.com