[petsc-dev] Issues with DMComplexSetCone

[Matthew Knepley]

On Fri, Aug 24, 2012 at 3:48 PM, Blaise Bourdin <bourdin at lsu.edu> wrote:

> I am attaching a debugged version.
>

Checked in as DMComplex ex1f90.F

   Matt


> Blaise
>
>
>
> On Aug 24, 2012, at 3:30 PM, Matthew Knepley <knepley at gmail.com> wrote:
>
> On Fri, Aug 24, 2012 at 2:16 PM, Chris Eldred <chris.eldred at gmail.com>wrote:
>
>> The following code snippet does not appear to be working:
>>
>
> 1) I hate Fortran. A lot. It is not debuggable.
>
> 2) Here is my code that should work, but the do loop somehow does not work:
>
>        program main
>       implicit none
> !
> #include <finclude/petsc.h90>
> #include <finclude/petscdmcomplex.h90>
> !
>       DM dm
>       PetscInt, target, dimension(4) :: EC
>       PetscInt, pointer :: pEC(:)
>       PetscInt c, firstCell, numCells
>       PetscErrorCode ierr
>
>       call PetscInitialize(PETSC_NULL_CHARACTER,ierr)
>       call DMComplexCreate(PETSC_COMM_WORLD, dm, ierr)
>       firstCell = 0
>       numCells = 5
>       call DMComplexSetChart(dm, 0, numCells, ierr)
>       do c=firstCell,numCells
>          write(*,*) 'c',c
>          call DMComplexSetConeSize(dm, c, 4)
>       end do
>       call DMSetUp(dm, ierr)
>
>       EC(1) = 1
>       EC(2) = 2
>       EC(3) = 3
>       EC(4) = 4
>       pEC => EC
>       write(*,*) 'cell',c,pEC
>       c = 0
>       call DMComplexSetCone(dm, c , pEC, ierr)
>       CHKERRQ(ierr)
>       call DMComplexGetCone(dm, c , pEC, ierr)
>       CHKERRQ(ierr)
>       write(*,*) 'cell',c,pEC
>
>       call DMDestroy(dm,ierr)
>       call PetscFinalize(ierr)
>       end
>
>
>    Matt
>
>
>> 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
>>
>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
>
> --
> Department of Mathematics and Center for Computation & Technology
> Louisiana State University, Baton Rouge, LA 70803, USA
> Tel. +1 (225) 578 1612, Fax  +1 (225) 578 4276
> http://www.math.lsu.edu/~bourdin
>
>
>
>
>
>
>
>
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120824/59464a9b/attachment.html>