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
