[petsc-dev] Issues with DMComplex Cone Orientations in Fortran

[Chris Eldred]

Building on the previous example code (for the Doublet mesh), I am
getting an error when trying to set Cone orientations:

[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Argument out of range!
[0]PETSC ERROR: Cone orientation 0 is not in the valid range [-1. 0)!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Development HG revision:
b6fe82991deee4a0d3f9a20654bc7750a6b1fe0f  HG Date: Mon Aug 27 13:17:10
2012 -0500
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: ./test on a arch-linu named Puget-101334 by user Tue
Aug 28 15:54:36 2012
[0]PETSC ERROR: Libraries linked from
/home/user/Desktop/LIBRARIES/petsc-dev/arch-linux2-cxx-debug/lib
[0]PETSC ERROR: Configure run at Mon Aug 27 12:55:14 2012
[0]PETSC ERROR: Configure options --download-boost --download-chaco
--download-ctetgen --download-f-blas-lapack --download-fiat
--download-generator --download-metis --download-ml --download-mpich
--download-parmetis --download-scientificpython --download-triangle
--with-clanguage=cxx --with-dynamic-loading --with-shared-libraries
--with-sieve PETSC_ARCH=arch-linux2-cxx-debug
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: DMComplexSetConeOrientation() line 1207 in
/home/user/Desktop/LIBRARIES/petsc-dev/src/dm/impls/complex/complex.c
application called MPI_Abort(MPI_COMM_WORLD, 63) - process 0
[unset]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 63) - process 0

The code is:

Orientations(1) = 0
Orientations(2) = 1
pOrientations => Orientations
call DMComplexSetConeOrientation(dm,4,pOrientations,ierr)
CHKERRQ(ierr)
call DMComplexGetConeOrientation(dm,4,pOrientations,ierr)
CHKERRQ(ierr)
write(*,*) pOrientations

Mesh point 4 is an edge (so it has a cone of size 2)- It is unclear to
me why the cone orientation range is [-1,0). My understanding of cone
orientations is that they specified the order in which cone's were
traversed for a point- ie for mesh point 4 (an edge) I could give it a
cone Orientation of (1,0) which would mean look at point 1 of the cone
first, then point 0.

-- 
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