[petsc-dev] Issues with DMComplexGetTransitiveClosure in Fortran

2012-08-30 Thread Matthew Knepley
On Thu, Aug 30, 2012 at 11:49 AM, Chris Eldred wrote: > That is what I have already been doing- so that is good! As far as the > ordering convention goes, for a 2D mesh I should have vertices, faces > (what I have been calling cells) and then edges? Ex. for the doublet > mesh it would be 0,1,2,3 (

[petsc-dev] Issues with DMComplexGetTransitiveClosure in Fortran

2012-08-30 Thread Chris Eldred
Oh ok that makes sense. The nomenclature was confusing me. On Thu, Aug 30, 2012 at 10:53 AM, Matthew Knepley wrote: > On Thu, Aug 30, 2012 at 11:49 AM, Chris Eldred > wrote: >> >> That is what I have already been doing- so that is good! As far as the >> ordering convention goes, for a 2D mesh I

[petsc-dev] Issues with DMComplexGetTransitiveClosure in Fortran

2012-08-30 Thread Chris Eldred
That is what I have already been doing- so that is good! As far as the ordering convention goes, for a 2D mesh I should have vertices, faces (what I have been calling cells) and then edges? Ex. for the doublet mesh it would be 0,1,2,3 (vertices); 4,5 (faces); 6,7,8,9,10 (edges)? On Wed, Aug 29, 20

[petsc-dev] Issues with DMComplexGetTransitiveClosure in Fortran

2012-08-29 Thread Matthew Knepley
On Wed, Aug 29, 2012 at 1:35 PM, Chris Eldred wrote: > DMComplexGetTransitiveClosure (and Restore) work fine now. What are > the new restrictions on numbering mentioned in > src/dm/impls/complex/examples/tests/ex2f90.F? > In the original Sieve, you could use any id for any point. This is very fle

[petsc-dev] Issues with DMComplexGetTransitiveClosure in Fortran

2012-08-29 Thread Chris Eldred
DMComplexGetTransitiveClosure (and Restore) work fine now. What are the new restrictions on numbering mentioned in src/dm/impls/complex/examples/tests/ex2f90.F? On Tue, Aug 28, 2012 at 8:16 PM, Matthew Knepley wrote: > On Tue, Aug 28, 2012 at 5:01 PM, Chris Eldred > wrote: >> >> On Tue, Aug 28,

[petsc-dev] Issues with DMComplexGetTransitiveClosure in Fortran

2012-08-28 Thread Matthew Knepley
On Tue, Aug 28, 2012 at 5:01 PM, Chris Eldred wrote: > On Tue, Aug 28, 2012 at 3:49 PM, Matthew Knepley > wrote: > > On Tue, Aug 28, 2012 at 4:39 PM, Chris Eldred > > wrote: > >> > >> DMComplexGetTransitiveClosure does not appear to have the calling > >> sequence (DM dm, PetscInt p, PetscBool us

[petsc-dev] Issues with DMComplexGetTransitiveClosure in Fortran

2012-08-28 Thread Matthew Knepley
On Tue, Aug 28, 2012 at 4:39 PM, Chris Eldred wrote: > DMComplexGetTransitiveClosure does not appear to have the calling > sequence (DM dm, PetscInt p, PetscBool useCone, PetscInt *numPoints, > PetscInt *points[]) indicated in the documentation- instead it has DM > dm, PetscInt p, PetscBool useCon

[petsc-dev] Issues with DMComplexGetTransitiveClosure in Fortran

2012-08-28 Thread Chris Eldred
On Tue, Aug 28, 2012 at 3:49 PM, Matthew Knepley wrote: > On Tue, Aug 28, 2012 at 4:39 PM, Chris Eldred > wrote: >> >> DMComplexGetTransitiveClosure does not appear to have the calling >> sequence (DM dm, PetscInt p, PetscBool useCone, PetscInt *numPoints, >> PetscInt *points[]) indicated in the

[petsc-dev] Issues with DMComplexGetTransitiveClosure in Fortran

2012-08-28 Thread Chris Eldred
DMComplexGetTransitiveClosure does not appear to have the calling sequence (DM dm, PetscInt p, PetscBool useCone, PetscInt *numPoints, PetscInt *points[]) indicated in the documentation- instead it has DM dm, PetscInt p, PetscBool useCone, PetscInt *points[]. In addition, it does not appear to be r